[gdal-dev] gdalwarp EPSG not being applied

[Gavin]

On Tue, 2011-01-18 at 10:06 -0700, Kyle Shannon wrote:

> Gavin,
> Is GDAL_DATA environment variable set on the Fedora system?


yip: GDAL_DATA=/opt/gdal_data


> 
> 
> 
> On Tue, Jan 18, 2011 at 04:15, Gavin <gavin at afrispatial.co.za> wrote:
> 
>         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 at lists.osgeo.org
>         http://lists.osgeo.org/mailman/listinfo/gdal-dev
> 
> 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20110118/03eda5fb/attachment-0001.html