Le vendredi 05 octobre 2012 22:21:27, David Shean a écrit : > 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.
Repeatable on various large float32 LZW-compressed GTiff, or that particular one ? > 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. This is something very likely related to update-in-place TIFF support, so generally not in the "simple fix" category. You might want to re-check with GDAL trunk whose internal libtiff has received per http://trac.osgeo.org/gdal/ticket/4771 a fix in that area (no promise that it fixes your particular issue). Otherwise, if you can provide a link to the dataset, I can try to see if I reproduce. That kind of bug really depends on how the imagery compress. I have created a fake dataset with similar characteristics to yours with : $ python from osgeo import gdal ds = gdal.GetDriverByName('GTiff').Create('buglzw.tif', 17719, 14955, \ 1, gdal.GDT_Float32, options = ['TILED=YES', 'COMPRESS=LZW']) ds.GetRasterBand(1).SetNoDataValue(0) ds = None And : $ gdaladdo -r gauss buglzw.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