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

Fri Apr 17 18:46:08 EDT 2009

Stephen,

Which version gives the right result ?
I've checked that GDAL SVN trunk (1.7.0dev) gives the same result as the one 
you get with GDAL 1.5.2.
It looks like EPSG:900913 definition was added in data/cubewerx_extra.wkt for 
GDAL 1.5.0. I've tested again with GDAL 1.4.4 and it doesn't even manage to 
start the warping :
"ERROR 6: EPSG PCS/GCS code 900913 not found in EPSG support files.  Is this a 
valid
EPSG coordinate system?
ERROR 1: Translating source or target SRS failed:
EPSG:900913
"

So I don't understand how you can get it to work with GDAL 1.4.4.

Even

Le Saturday 18 April 2009 00:14:09 Stephen Woodbridge, vous avez écrit :
> 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
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev