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