[gdal-dev] gdalwarp differences between v1.5.2 and v1.4.4.0

[Stephen Woodbridge]

Hi all,

I have a serious problem that you might be able to shed some light on.

I have a landsat geotif (actually lots of them) and when I projection 
them on a systems with v1.5.2 and v1.4.4.0 I get very different results!

The image sizes are different, the projection info is different, the 
pixel size is different.

And I think the imagery had shifted significantly in the vertical direction.

Thoughts on how to fix this using 1.4.4.0, I'm not sure I can upgrade 
the system because of package dependencies.

-Steve


Original image same on both systems:

$ gdalinfo n-41-25_0-0.tif
Driver: GTiff/GeoTIFF
Files: n-41-25_0-0.tif
Size is 27192, 18524
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 = (59.228401534356898,27.479955940770701)
Pixel Size = (0.000138699063483,-0.000138699063483)
Metadata:
   AREA_OR_POINT=Area
Image Structure Metadata:
   INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  59.2284015,  27.4799559) ( 59d13'42.25"E, 27d28'47.84"N)
Lower Left  (  59.2284015,  24.9106945) ( 59d13'42.25"E, 24d54'38.50"N)
Upper Right (  62.9999065,  27.4799559) ( 62d59'59.66"E, 27d28'47.84"N)
Lower Right (  62.9999065,  24.9106945) ( 62d59'59.66"E, 24d54'38.50"N)
Center      (  61.1141540,  26.1953252) ( 61d 6'50.95"E, 26d11'43.17"N)
Band 1 Block=256x256 Type=Byte, ColorInterp=Red
   NoData Value=0
   Overviews: 13596x9262, 6798x4631, 3399x2316, 1700x1158, 850x579, 
425x290, 213x145, 107x73, 54x37
Band 2 Block=256x256 Type=Byte, ColorInterp=Green
   NoData Value=0
   Overviews: 13596x9262, 6798x4631, 3399x2316, 1700x1158, 850x579, 
425x290, 213x145, 107x73, 54x37
Band 3 Block=256x256 Type=Byte, ColorInterp=Blue
   NoData Value=0
   Overviews: 13596x9262, 6798x4631, 3399x2316, 1700x1158, 850x579, 
425x290, 213x145, 107x73, 54x37


I run this command on the two systems:

$ gdalwarp -s_srs EPSG:4326 -t_srs epsg:900913 -co TILED=YES 
landsat-af2/n-41-25_0-0.tif junk/n-41-25_0-0.tif

I get results on GDAL 1.5.2, released 2008/05/29:

$ gdalinfo junk/n-41-25_0-0.tif
Driver: GTiff/GeoTIFF
Files: junk/n-41-25_0-0.tif
Size is 26204, 19897
Coordinate System is:
PROJCS["Google Maps Global Mercator",
     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"]],
     PROJECTION["Mercator_1SP"],
     PARAMETER["central_meridian",0],
     PARAMETER["scale_factor",1],
     PARAMETER["false_easting",0],
     PARAMETER["false_northing",0],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]]]
Origin = (6593275.499304148368537,3183564.941239751409739)
Pixel Size = (16.021935325381026,-16.021935325381026)
Metadata:
   AREA_OR_POINT=Area
Image Structure Metadata:
   INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  ( 6593275.499, 3183564.941) ( 59d13'42.25"E, 27d38'15.97"N)
Lower Left  ( 6593275.499, 2864776.494) ( 59d13'42.25"E, 25d 3'28.72"N)
Upper Right ( 7013114.293, 3183564.941) ( 62d59'59.56"E, 27d38'15.97"N)
Lower Right ( 7013114.293, 2864776.494) ( 62d59'59.56"E, 25d 3'28.72"N)
Center      ( 6803194.896, 3024170.718) ( 61d 6'50.90"E, 26d21'18.52"N)
Band 1 Block=256x256 Type=Byte, ColorInterp=Red
Band 2 Block=256x256 Type=Byte, ColorInterp=Green
Band 3 Block=256x256 Type=Byte, ColorInterp=Blue


I get results on GDAL 1.4.4.0, released 2007/11/23:

$ gdalinfo junk/n-41-25_0-0.tif
Driver: GTiff/GeoTIFF
Size is 26256, 19828
Coordinate System is:
PROJCS["unnamed",
     GEOGCS["unnamed ellipse",
         DATUM["unknown",
             SPHEROID["unnamed",6378137,0]],
         PRIMEM["Greenwich",0],
         UNIT["degree",0.0174532925199433]],
     PROJECTION["Mercator_1SP"],
     PARAMETER["latitude_of_origin",0],
     PARAMETER["central_meridian",0],
     PARAMETER["scale_factor",1],
     PARAMETER["false_easting",0],
     PARAMETER["false_northing",0],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]]]
Origin = (6593275.499304148368537,3163848.570181989576668)
Pixel Size = (15.990354027000551,-15.990354027000551)
Metadata:
   AREA_OR_POINT=Area
Corner Coordinates:
Upper Left  ( 6593275.499, 3163848.570) ( 59d13'42.25"E, 27d19'21.76"N)
Lower Left  ( 6593275.499, 2846791.831) ( 59d13'42.25"E, 24d45'50.60"N)
Upper Right ( 7013118.235, 3163848.570) ( 62d59'59.69"E, 27d19'21.76"N)
Lower Right ( 7013118.235, 2846791.831) ( 62d59'59.69"E, 24d45'50.60"N)
Center      ( 6803196.867, 3005320.200) ( 61d 6'50.97"E, 26d 3'1.31"N)
Band 1 Block=256x256 Type=Byte, ColorInterp=Red
Band 2 Block=256x256 Type=Byte, ColorInterp=Green
Band 3 Block=256x256 Type=Byte, ColorInterp=Blue