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

Stephen Woodbridge woodbri at swoodbridge.com
Sat Apr 18 09:06:48 EDT 2009 

Even,

I copied the newer data/cubewerx_extra.wkt over to the gdal 1.4.4 system 
and that seems to have fixed the problem.

I'm curious why just adding 900913 to epsg did not work for gdal, but 
did work for the mapserver which is projection the vectors on the fly.

-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