[gdal-dev] gdaladdo bus errors for large, LZW-compressed GTiff

[David Shean]

While the topic of gdaladdo resampling methods is active...

I get bus errors when I use gdaladdo's gauss resampling method for large, float32, LZW-compressed GTiff.  This is repeatable.
Using a different gdaladdo resampling method (e.g., average) on the same LZW-compressed GTiff works.
When I decompress first, gdaladdo with gauss works.
For smaller, LZW-compressed GTiff input, gdaladdo with gauss works.

I'm using gdal 1.9.1 with internal libtiff, libgeotiff, libz on OS X 10.7.

Hopefully this will involve a simple fix.  Thanks.
-David



gdalinfo WV02_10JUL091534204-P1BS-1030010005B3D100__WV02_10JUL091533199-P1BS-103001000612DC00-DEM.tif
Driver: GTiff/GeoTIFF
Files: WV02_10JUL091534204-P1BS-1030010005B3D100__WV02_10JUL091533199-P1BS-103001000612DC00-DEM.tif
Size is 17719, 14955
Coordinate System is:
PROJCS["unnamed",
   GEOGCS["WGS 84",
       DATUM["WGS_1984",
           SPHEROID["WGS 84",6378137,298.257223563,
               AUTHORITY["EPSG","7030"]],
           AUTHORITY["EPSG","6326"]],
       PRIMEM["Greenwich",0],
       UNIT["degree",0.0174532925199433],
       AUTHORITY["EPSG","4326"]],
   PROJECTION["Polar_Stereographic"],
   PARAMETER["latitude_of_origin",70],
   PARAMETER["central_meridian",-45],
   PARAMETER["scale_factor",1],
   PARAMETER["false_easting",0],
   PARAMETER["false_northing",0],
   UNIT["metre",1,
       AUTHORITY["EPSG","9001"]]]
Origin = (-201111.404963798559038,-2284012.343544036149979)
Pixel Size = (1.000000000000000,-1.000000000000000)
Metadata:
 AREA_OR_POINT=Area
Image Structure Metadata:
 COMPRESSION=LZW
 INTERLEAVE=BAND
Corner Coordinates:
Upper Left  ( -201111.405,-2284012.344) ( 50d 1'55.27"W, 69d 3'39.06"N)
Lower Left  ( -201111.405,-2298967.344) ( 49d59'58.03"W, 68d55'39.79"N)
Upper Right ( -183392.405,-2284012.344) ( 49d35'26.36"W, 69d 4'26.88"N)
Lower Right ( -183392.405,-2298967.344) ( 49d33'39.31"W, 68d56'27.28"N)
Center      ( -192251.905,-2291489.844) ( 49d47'44.82"W, 69d 0' 3.78"N)
Band 1 Block=256x256 Type=Float32, ColorInterp=Gray
 NoData Value=0

LZW-compressed, -r gauss:

gdaladdo -r gauss WV02_10JUL091534204-P1BS-1030010005B3D100__WV02_10JUL091533199-P1BS-103001000612DC00-DEM.tif 2 4 8 16 32
0Bus error

LZW-compressed, -r average:

gdaladdo -r average WV02_10JUL091534204-P1BS-1030010005B3D100__WV02_10JUL091533199-P1BS-103001000612DC00-DEM.tif 2 4 8 16 32
0...10...20...30...40...50...60...70...80...90...100 - done.

Decompressed, -r gauss:

gdaladdo -r gauss WV02_10JUL091534204-P1BS-1030010005B3D100__WV02_10JUL091533199-P1BS-103001000612DC00-DEM_nolzw.tif 2 4 8 16 32
0...10...20...30...40...50...60...70...80...90...100 - done.