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.
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
