[gdal-dev] gdalwarp EPSG not being applied

Kyle Shannon ksshannon at gmail.com
Tue Jan 18 12:06:34 EST 2011


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

# ============================
Kyle Shannon
Physical Science Technician
RMRS Fire Sciences Lab
Fire, Fuels & Smoke - RWU 4405
5775 Highway 10 W.
Missoula, MT 59808
(406)829-6954
kshannon at fs.fed.us
# ============================


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/2188d40b/attachment.html


More information about the gdal-dev mailing list