[gdal-dev] gdalwarp EPSG not being applied

[Gavin]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20110118/c5920405/attachment.html