gdalwarp with -t_srs 'EPSG:3857' produces expected output on one system
but not on another. It seems gdal can't find the projection definition
on System 2 (below). System 2 is the one where I need it to work and
where I ultimately need to publish as 900913 through geoserver.
What am I missing?
This command was run on both systems (both have gdal 1.7 and proj
4.9.1):
gdalwarp -s_srs "+proj=tmerc +lat_0=0 +lon_0=27 +k=1 +x_0=0 +y_0=0
+ellps=WGS84 +datum=WGS84 +units=m +no_defs" -t_srs 'EPSG:3857' -r near
-multi -dstnodata "0 0 0" -dstnodata "0 0 0" -dstalpha -co
COMPRESS="DEFLATE" -co "TILED=YES" -of GTiff
2627BB_05_1117_RGB_RECT.TIF /tmp/test900913.tif
Geotiff header output:
System1 (Ubuntu 10.10):
Driver: GTiff/GeoTIFF
Files: /tmp/test.tif
Size is 10749, 11713
Coordinate System is:
PROJCS["WGS 84 / Pseudo-Mercator",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
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],
AUTHORITY["EPSG","3857"],
EXTENSION["PROJ4","+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 +wktext +no_defs"],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (3111167.042467535007745,-2998912.428198689594865)
Pixel Size = (0.557611099232479,-0.557611099232479)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
COMPRESSION=DEFLATE
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( 3111167.042,-2998912.428) ( 27d56'53.12"E,
25d59'55.10"S)
Lower Left ( 3111167.042,-3005443.727) ( 27d56'53.12"E, 26d
3'4.90"S)
Upper Right ( 3117160.804,-2998912.428) ( 28d 0'6.95"E,
25d59'55.10"S)
Lower Right ( 3117160.804,-3005443.727) ( 28d 0'6.95"E, 26d
3'4.90"S)
Center ( 3114163.923,-3002178.078) ( 27d58'30.04"E, 26d
1'30.01"S)
System 2 (Fedora 10):
Driver: GTiff/GeoTIFF
Files: test3857.tif
Size is 10749, 11713
Coordinate System is:
PROJCS["unnamed",
GEOGCS[,
DATUM["unknown",
SPHEROID["unretrievable - using
WGS84",6378137,298.257223563]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","3857"]]
Origin = (3111167.042467534542084,-2998912.428198689594865)
Pixel Size = (0.557611099232459,-0.557611099232459)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
COMPRESSION=DEFLATE
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( 3111167.042,-2998912.428)
Lower Left ( 3111167.042,-3005443.727)
Upper Right ( 3117160.804,-2998912.428)
Lower Right ( 3117160.804,-3005443.727)
Center ( 3114163.923,-3002178.078)
thanks
Gavin
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
