Dear all! This may be related to my recent posting about SetField behaving differently in different gdal versions. In any case I find that the bounding box of this test geotif file<http://www.keepandshare.com/doc/3116689/population-2010-clip-tif-september-3-2011-5-20-pm-162k?da=y> is different on the two installations (details below).
In brief, the bounding box using gdal 1.6 is Upper Left ( 99.3641696, -0.0041806) ( 99d21'51.01"E, 0d 0'15.05"S) Lower Left ( 99.3641696, -2.2031806) ( 99d21'51.01"E, 2d12'11.45"S) Upper Right ( 102.2411696, -0.0041806) (102d14'28.21"E, 0d 0'15.05"S) Lower Right ( 102.2411696, -2.2031806) (102d14'28.21"E, 2d12'11.45"S) Center ( 100.8026696, -1.1036806) (100d48'9.61"E, 1d 6'13.25"S) whereas using gdal 1.8 it is Upper Left ( 99.3600000, 0.0000000) ( 99d21'36.00"E, 0d 0' 0.01"N) Lower Left ( 99.3600000, -2.1990000) ( 99d21'36.00"E, 2d11'56.40"S) Upper Right ( 102.2370000, 0.0000000) (102d14'13.20"E, 0d 0' 0.01"N) Lower Right ( 102.2370000, -2.1990000) (102d14'13.20"E, 2d11'56.40"S) Center ( 100.7985000, -1.0995000) (100d47'54.60"E, 1d 5'58.20"S) It'd be great if someone can 1. Reproduce this using the test file (download using link above) 2. Let us know which bounding box is correct 3. Raise this as a ticket if it is a genuine bug Cheers and thanks Ole Nielsen -------------------------------------------- DETAILED OUTPUTS -------------------------------------------- --------------------------------------------- On Maverick 10.10, GDAL 1.6 --------------------------------------------- nielso@shiva:~/dev/risiko_test_data$ dpkg -l | grep gdal ii gdal-bin 1.6.3-4build1 Geospatial Data Abstraction Library - Utility programs ii libgdal1-1.6.0 1.6.3-4build1 Geospatial Data Abstraction Library ii libgdal1-1.6.0-grass 1.6.3-3 GRASS extension for the GDAL library ii libgdal1-dev 1.6.3-4build1 Geospatial Data Abstraction Library - Development files ii python-gdal 1.6.3-4build1 Python bindings to the Geospatial Data Abstraction Library nielso@shiva:~/dev/risiko_test_data$ gdalinfo --version GDAL 1.6.3, released 2009/11/19 nielso@shiva:~/dev/risiko_test_data$ gdalinfo Earthquake_Ground_Shaking_clip.tif Driver: GTiff/GeoTIFF Files: Earthquake_Ground_Shaking_clip.tif Size is 345, 263 Coordinate System is: GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.2572235630016, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4326"]] Origin = (99.364169565217395,-0.004180608365019) Pixel Size = (0.008339130434783,-0.008361216730038) Metadata: AREA_OR_POINT=Point TIFFTAG_XRESOLUTION=1 TIFFTAG_YRESOLUTION=1 TIFFTAG_RESOLUTIONUNIT=1 (unitless) Image Structure Metadata: COMPRESSION=LZW INTERLEAVE=BAND Corner Coordinates: Upper Left ( 99.3641696, -0.0041806) ( 99d21'51.01"E, 0d 0'15.05"S) Lower Left ( 99.3641696, -2.2031806) ( 99d21'51.01"E, 2d12'11.45"S) Upper Right ( 102.2411696, -0.0041806) (102d14'28.21"E, 0d 0'15.05"S) Lower Right ( 102.2411696, -2.2031806) (102d14'28.21"E, 2d12'11.45"S) Center ( 100.8026696, -1.1036806) (100d48'9.61"E, 1d 6'13.25"S) Band 1 Block=256x256 Type=Float64, ColorInterp=Gray ----------------------------------------- On Natty 11.04, GDAL 1.8.0 ----------------------------------------- nielso@shakti:~/dev/risiko_test_data$ dpkg -l | grep gdal ii gdal-bin 1.8.0-2~natty2 Geospatial Data Abstraction Library - Utility programs pi libgdal1-1.6.0 1.6.3-4build4 Geospatial Data Abstraction Library pi libgdal1-1.8.0 1.8.0-2~natty2 Geospatial Data Abstraction Library ii python-gdal 1.8.0-2~natty2 Python bindings to the Geospatial Data Abstraction Library nielso@shakti:~/dev/risiko_test_data$ gdalinfo --version GDAL 1.8.0, released 2011/01/12 nielso@shakti:~/dev/risiko_test_data$ gdalinfo Earthquake_Ground_Shaking_clip.tif Driver: GTiff/GeoTIFF Files: Earthquake_Ground_Shaking_clip.tif Size is 345, 263 Coordinate System is: 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"]] Origin = (99.359999999999999,0.000000000000000) Pixel Size = (0.008339130434783,-0.008361216730038) Metadata: TIFFTAG_XRESOLUTION=1 TIFFTAG_YRESOLUTION=1 TIFFTAG_RESOLUTIONUNIT=1 (unitless) AREA_OR_POINT=Point Image Structure Metadata: COMPRESSION=LZW INTERLEAVE=BAND Corner Coordinates: Upper Left ( 99.3600000, 0.0000000) ( 99d21'36.00"E, 0d 0' 0.01"N) Lower Left ( 99.3600000, -2.1990000) ( 99d21'36.00"E, 2d11'56.40"S) Upper Right ( 102.2370000, 0.0000000) (102d14'13.20"E, 0d 0' 0.01"N) Lower Right ( 102.2370000, -2.1990000) (102d14'13.20"E, 2d11'56.40"S) Center ( 100.7985000, -1.0995000) (100d47'54.60"E, 1d 5'58.20"S) Band 1 Block=256x256 Type=Float64, ColorInterp=Gray
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev