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
