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

Stephen Woodbridge Fri, 17 Apr 2009 16:41:51 -0700

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 +nadgri...@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 oneyou get with GDAL 1.5.2.It looks like EPSG:900913 definition was added in data/cubewerx_extra.wkt forGDAL 1.5.0. I've tested again with GDAL 1.4.4 and it doesn't even manage tostart the warping :"ERROR 6: EPSG PCS/GCS code 900913 not found in EPSG support files. Is this avalid

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

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

Reply via email to