[gdal-dev] Fwd: Bounding box change in recent versions of gdal?
Ole Nielsen
ole.moller.nielsen at gmail.com
Sat Sep 3 06:23:55 EDT 2011
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 at 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 at shiva:~/dev/risiko_test_data$ gdalinfo --version
GDAL 1.6.3, released 2009/11/19
nielso at 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 at 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 at shakti:~/dev/risiko_test_data$ gdalinfo --version
GDAL 1.8.0, released 2011/01/12
nielso at 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20110903/38eb3456/attachment.html
More information about the gdal-dev
mailing list