<font face="courier new,monospace">Gavin,<br>Is GDAL_DATA environment variable set on the Fedora system?<br clear="all"></font><br># ============================<br>Kyle Shannon<br>Physical Science Technician<br>RMRS Fire Sciences Lab<br>
Fire, Fuels & Smoke - RWU 4405<br>5775 Highway 10 W.<br>Missoula, MT 59808<br>(406)829-6954<br><a href="mailto:kshannon@fs.fed.us" target="_blank">kshannon@fs.fed.us</a><br># ============================<br>
<br><br><div class="gmail_quote">On Tue, Jan 18, 2011 at 04:15, Gavin <span dir="ltr"><<a href="mailto:gavin@afrispatial.co.za">gavin@afrispatial.co.za</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div>
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><font color="#888888">
<br>
Gavin
</font></div>
<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" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br></blockquote></div><br>