[postgis-tickets] [PostGIS] #3020: [raster] ST_AddBand out-db bug with height

[PostGIS]

#3020: [raster] ST_AddBand out-db bug with height
-----------------------+----------------------------------------------------
 Reporter:  dustymugs  |       Owner:  dustymugs    
     Type:  defect     |      Status:  new          
 Priority:  critical   |   Milestone:  PostGIS 2.1.6
Component:  raster     |     Version:  2.1.x        
 Keywords:             |  
-----------------------+----------------------------------------------------
 From an email:

 http://lists.osgeo.org/pipermail/postgis-devel/2015-January/024713.html

 Hello
 I’m writing this message, because of bug i found in posgis raster.
 (I was unable to create user id on osgeo trac, and as result to create
 ticket)

 The bug occurs in ST_AddBand function for out-db raster.
 I used this function in fallowing way:

 select ST_AddBand(null,'D:\\Mapy\\hrpt_test_12_01_2015.tif',ARRAY[1, 2, 3,
 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14]);

 And I received such error:
 rt_raster_add_band: Can't add a 1200x1200 band to a 1200x1600 raster

 Where hrpt_test_12_01_2015.tif is 1200px x 1600px valid GeoTiff with 14
 bands, and it gdalinfo is as fallows:
 Driver: GTiff/GeoTIFF
 Files: hrpt_test_12_01_2015.tif
 Size is 1200, 1600
 Coordinate System is:
 PROJCS["Mercator, WGS 84, Bands: Solar zenith°,Solar azimuth°,Nadir to
 pixel°,Satellite inclination°,AVHRR channel 1 - raw,AVHRR channel 1 -
 albedo%,AVHRR channel 2 - raw,AVHRR channel 2 - albedo%,AVHRR channel 3 -
 raw,AVHRR channel 3 - temperature°C,AVHRR channel 4 - raw,AVHRR channel 4
 - temperature°C,AVHRR channel 5 - raw,AVHRR channel 5 - temperature°C",
     GEOGCS["Unknown datum based upon the WGS 84 ellipsoid",
         DATUM["Not_specified_based_on_WGS_84_ellipsoid",
             SPHEROID["WGS 84",6378137,298.257223563,
                 AUTHORITY["EPSG","7030"]],
             AUTHORITY["EPSG","6030"]],
         PRIMEM["Greenwich",0],
         UNIT["degree",0.0174532925199433],
         AUTHORITY["EPSG","4030"]],
     PROJECTION["Mercator_1SP"],
     PARAMETER["central_meridian",0],
     PARAMETER["scale_factor",1],
     PARAMETER["false_easting",0],
     PARAMETER["false_northing",0],
     UNIT["unknown",1]]
 Origin = (1500000.000000000500000,7799999.999998852600000)
 Pixel Size = (1000.000000000000000,-1000.000000000000000)
 Metadata:
   AREA_OR_POINT=Area
   TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
   TIFFTAG_XRESOLUTION=96
   TIFFTAG_YRESOLUTION=96
 Image Structure Metadata:
   INTERLEAVE=BAND
 Corner Coordinates:
 Upper Left  ( 1500000.000, 7800000.000) ( 13d28'29.03"E, 57d22'10.14"N)
 Lower Left  ( 1500000.000, 6200000.000) ( 13d28'29.03"E, 48d44'51.33"N)
 Upper Right ( 2700000.000, 7800000.000) ( 24d15'16.25"E, 57d22'10.14"N)
 Lower Right ( 2700000.000, 6200000.000) ( 24d15'16.25"E, 48d44'51.33"N)
 Center      ( 2100000.000, 7000000.000) ( 18d51'52.64"E, 53d16'34.72"N)
 Band 1 Block=1200x1 Type=Float32, ColorInterp=Gray
 Band 2 Block=1200x1 Type=Float32, ColorInterp=Undefined
 Band 3 Block=1200x1 Type=Float32, ColorInterp=Undefined
 Band 4 Block=1200x1 Type=Float32, ColorInterp=Undefined
 Band 5 Block=1200x1 Type=Float32, ColorInterp=Undefined
 Band 6 Block=1200x1 Type=Float32, ColorInterp=Undefined
 Band 7 Block=1200x1 Type=Float32, ColorInterp=Undefined
 Band 8 Block=1200x1 Type=Float32, ColorInterp=Undefined
 Band 9 Block=1200x1 Type=Float32, ColorInterp=Undefined
 Band 10 Block=1200x1 Type=Float32, ColorInterp=Undefined
 Band 11 Block=1200x1 Type=Float32, ColorInterp=Undefined
 Band 12 Block=1200x1 Type=Float32, ColorInterp=Undefined
 Band 13 Block=1200x1 Type=Float32, ColorInterp=Undefined
 Band 14 Block=1200x1 Type=Float32, ColorInterp=Undefined


 As you can see bug occurs for non square rasters.

 I tired to track this bug in source code, and I found somethig. It is
 located in rtpg_create.c file in RASTER_addBandOutDB() function, line 776.
 http://trac.osgeo.org/postgis/browser/trunk/raster/rt_pg/rtpg_create.c#L776
 Width of raster is written to height variable, which causes problem.

 I hope this will be helpful.

 I’m using postgis:
 “POSTGIS="2.1.5 r13152" GEOS="3.4.2-CAPI-1.8.2 r3924" PROJ="Rel. 4.8.0, 6
 March 2012" GDAL="GDAL 1.11.1, released 2014/09/24" LIBXML="2.7.8"
 LIBJSON="UNKNOWN" RASTER”,
 and postgres:
 "PostgreSQL 9.4.0, compiled by Visual C++ build 1800, 64-bit"

 Best Regards
 Tomasz Bieliński
 Gdańsk University of Technology

-- 
Ticket URL: <http://trac.osgeo.org/postgis/ticket/3020>
PostGIS <http://trac.osgeo.org/postgis/>
The PostGIS Trac is used for bug, enhancement & task tracking, a user and developer wiki, and a view into the subversion code repository of PostGIS project.