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

Stephen Woodbridge woodbri at swoodbridge.com
Fri Apr 17 19:40:58 EDT 2009


Even,

gdal 1.5.0 gives the correct result.

I added the 900913 definition to /usr/share/proj/epsg files as:

# Spherical Mercator - Google
<900913> +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 
+x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs

-Steve

Even Rouault wrote:
> 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
> 
> 



More information about the gdal-dev mailing list