<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 TRANSITIONAL//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; CHARSET=UTF-8">
<META NAME="GENERATOR" CONTENT="GtkHTML/3.30.3">
</HEAD>
<BODY>
<BR>
<BR>
On Tue, 2011-01-18 at 10:06 -0700, Kyle Shannon wrote:<BR>
<BLOCKQUOTE TYPE=CITE>
Gavin,<BR>
Is GDAL_DATA environment variable set on the Fedora system?<BR>
</BLOCKQUOTE>
<BR>
yip: GDAL_DATA=/opt/gdal_data<BR>
<BR>
<BLOCKQUOTE TYPE=CITE>
<BR>
<BR>
</BLOCKQUOTE>
<BLOCKQUOTE TYPE=CITE>
On Tue, Jan 18, 2011 at 04:15, Gavin <<A HREF="mailto:gavin@afrispatial.co.za">gavin@afrispatial.co.za</A>> wrote:
</BLOCKQUOTE>
<BLOCKQUOTE TYPE=CITE>
<BLOCKQUOTE>
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. <BR>
<BR>
What am I missing?<BR>
<BR>
This command was run on both systems (both have gdal 1.7 and proj 4.9.1):<BR>
<BR>
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<BR>
<BR>
Geotiff header output:<BR>
<BR>
System1 (Ubuntu 10.10):<BR>
<BR>
<BLOCKQUOTE>
Driver: GTiff/GeoTIFF<BR>
Files: /tmp/test.tif<BR>
Size is 10749, 11713<BR>
Coordinate System is:<BR>
PROJCS["WGS 84 / Pseudo-Mercator",<BR>
GEOGCS["WGS 84",<BR>
DATUM["WGS_1984",<BR>
SPHEROID["WGS 84",6378137,298.257223563,<BR>
AUTHORITY["EPSG","7030"]],<BR>
AUTHORITY["EPSG","6326"]],<BR>
PRIMEM["Greenwich",0],<BR>
UNIT["degree",0.0174532925199433],<BR>
AUTHORITY["EPSG","4326"]],<BR>
PROJECTION["Mercator_1SP"],<BR>
PARAMETER["central_meridian",0],<BR>
PARAMETER["scale_factor",1],<BR>
PARAMETER["false_easting",0],<BR>
PARAMETER["false_northing",0],<BR>
AUTHORITY["EPSG","3857"],<BR>
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"],<BR>
UNIT["metre",1,<BR>
AUTHORITY["EPSG","9001"]]]<BR>
Origin = (3111167.042467535007745,-2998912.428198689594865)<BR>
Pixel Size = (0.557611099232479,-0.557611099232479)<BR>
Metadata:<BR>
AREA_OR_POINT=Area<BR>
Image Structure Metadata:<BR>
COMPRESSION=DEFLATE<BR>
INTERLEAVE=PIXEL<BR>
Corner Coordinates:<BR>
Upper Left ( 3111167.042,-2998912.428) ( 27d56'53.12"E, 25d59'55.10"S)<BR>
Lower Left ( 3111167.042,-3005443.727) ( 27d56'53.12"E, 26d 3'4.90"S)<BR>
Upper Right ( 3117160.804,-2998912.428) ( 28d 0'6.95"E, 25d59'55.10"S)<BR>
Lower Right ( 3117160.804,-3005443.727) ( 28d 0'6.95"E, 26d 3'4.90"S)<BR>
Center ( 3114163.923,-3002178.078) ( 27d58'30.04"E, 26d 1'30.01"S)<BR>
</BLOCKQUOTE>
<BR>
System 2 (Fedora 10):<BR>
<BR>
<BLOCKQUOTE>
Driver: GTiff/GeoTIFF<BR>
Files: test3857.tif<BR>
Size is 10749, 11713<BR>
Coordinate System is:<BR>
PROJCS["unnamed",<BR>
GEOGCS[,<BR>
DATUM["unknown",<BR>
SPHEROID["unretrievable - using WGS84",6378137,298.257223563]],<BR>
PRIMEM["Greenwich",0],<BR>
UNIT["degree",0.0174532925199433]],<BR>
UNIT["metre",1,<BR>
AUTHORITY["EPSG","9001"]],<BR>
AUTHORITY["EPSG","3857"]]<BR>
Origin = (3111167.042467534542084,-2998912.428198689594865)<BR>
Pixel Size = (0.557611099232459,-0.557611099232459)<BR>
Metadata:<BR>
AREA_OR_POINT=Area<BR>
Image Structure Metadata:<BR>
COMPRESSION=DEFLATE<BR>
INTERLEAVE=PIXEL<BR>
Corner Coordinates:<BR>
Upper Left ( 3111167.042,-2998912.428) <BR>
Lower Left ( 3111167.042,-3005443.727) <BR>
Upper Right ( 3117160.804,-2998912.428) <BR>
Lower Right ( 3117160.804,-3005443.727) <BR>
Center ( 3114163.923,-3002178.078)<BR>
</BLOCKQUOTE>
<BR>
<BR>
thanks<BR>
<BR>
<FONT COLOR="#888888">Gavin</FONT>
</BLOCKQUOTE>
</BLOCKQUOTE>
<BLOCKQUOTE TYPE=CITE>
<BLOCKQUOTE>
<BR>
_______________________________________________<BR>
gdal-dev mailing list<BR>
<A HREF="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</A><BR>
<A HREF="http://lists.osgeo.org/mailman/listinfo/gdal-dev">http://lists.osgeo.org/mailman/listinfo/gdal-dev</A>
</BLOCKQUOTE>
</BLOCKQUOTE>
<BLOCKQUOTE TYPE=CITE>
<BR>
</BLOCKQUOTE>
</BODY>
</HTML>