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

Sat, 18 Apr 2009 06:07:52 -0700 

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
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev





_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev






















					Reply via email to