<HTML>
<HEAD>
<TITLE>GDAL/OGR projection problem only with epsg:3785/900913</TITLE>
</HEAD>
<BODY>
<FONT SIZE="2"><FONT FACE="Consolas, Courier New, Courier"><SPAN STYLE='font-size:10pt'>All,<BR>
<BR>
I’m using GDAL (1.9dev r22508) and getting errors only when reprojecting (raster or vector) to any of the Google Mercator projections. On a number of linux boxes I have, this works fine, but, on a couple its a problem. The only difference on the problem boxes that I can see is that they have the configuration files for vertical datum reprojection and are using external libgeotiff (svn r2002) for use with libLAS and PDAL. Reprojection using either the epsg codes or the full proj string fails.<BR>
<BR>
A couple of example errors<BR>
<BR>
[u4rt9mds@lidarora1 /] ogr2ogr -f "ESRI Shapefile" -t_srs epsg:3785 test.shp cincinnati_forests.shp<BR>
ERROR 1: no system list, errno: 2<BR>
<BR>
Failed to reproject feature 0 (geometry probably out of source or destination SRS).<BR>
ERROR 1: Terminating translation prematurely after failed<BR>
translation of layer cincinnati_forests (use -skipfailures to skip errors)<BR>
<BR>
[u4rt9mds@lidarora1 /] gdalwarp -of vrt -t_srs epsg:900913 -r near cincinnati_bare_earth_01.img test2.vrt<BR>
ERROR 1: no system list, errno: 2<BR>
<BR>
ERROR 1: GDALSuggestedWarpOutput() failed because the passed<BR>
transformer failed.<BR>
<BR>
Reprojection to any other projection seems to work fine<BR>
<BR>
[u4rt9mds@lidarora1 /] gdalwarp -of vrt -t_srs epsg:26919 /u02/LASdata/Office_of_Americas/Cincinnati_OH/Delivery/rasters/cincinnati_bare_earth_01.img test.vrt<BR>
Creating output file that is 20023P x 24321L.<BR>
Processing input file /u02/LASdata/Office_of_Americas/Cincinnati_OH/Delivery/rasters/cincinnati_bare_earth_01.img.<BR>
Using internal nodata values (eg. -3.40282e+38) for image /u02/LASdata/Office_of_Americas/Cincinnati_OH/Delivery/rasters/cincinnati_bare_earth_01.img.<BR>
<BR>
<BR>
[u4rt9mds@lidarora1 /] ogr2ogr -progress -f "ESRI Shapefile" -t_srs epsg:26917 test.shp cincinnati_forests.shp<BR>
0...10...20...30...40...50...60...70...80...90...100 - done.<BR>
<BR>
Testepsg reports the correct output<BR>
<BR>
[u4rt9mds@lidarora1 shapes]$ testepsg epsg:3785<BR>
Validate Succeeds.<BR>
WKT[epsg:3785] =<BR>
PROJCS["Popular Visualisation CRS / Mercator (deprecated)",<BR>
GEOGCS["Popular Visualisation CRS",<BR>
DATUM["Popular_Visualisation_Datum",<BR>
SPHEROID["Popular Visualisation Sphere",6378137,0,<BR>
AUTHORITY["EPSG","7059"]],<BR>
TOWGS84[0,0,0,0,0,0,0],<BR>
AUTHORITY["EPSG","6055"]],<BR>
PRIMEM["Greenwich",0,<BR>
AUTHORITY["EPSG","8901"]],<BR>
UNIT["degree",0.0174532925199433,<BR>
AUTHORITY["EPSG","9122"]],<BR>
AUTHORITY["EPSG","4055"]],<BR>
UNIT["metre",1,<BR>
AUTHORITY["EPSG","9001"]],<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>
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>
AUTHORITY["EPSG","3785"],<BR>
AXIS["X",EAST],<BR>
AXIS["Y",NORTH]]<BR>
<BR>
Simplified WKT[epsg:3785] =<BR>
PROJCS["Popular Visualisation CRS / Mercator (deprecated)",<BR>
GEOGCS["Popular Visualisation CRS",<BR>
DATUM["Popular_Visualisation_Datum",<BR>
SPHEROID["Popular Visualisation Sphere",6378137,0],<BR>
TOWGS84[0,0,0,0,0,0,0]],<BR>
PRIMEM["Greenwich",0],<BR>
UNIT["degree",0.0174532925199433]],<BR>
UNIT["metre",1],<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>
<BR>
Old Style WKT[epsg:3785] = PROJCS["Popular Visualisation CRS / Mercator (deprecated)",GEOGCS["Popular Visualisation CRS",DATUM["Popular_Visualisation_Datum",SPHEROID["Popular Visualisation Sphere",6378137,0]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],UNIT["metre",1],PROJECTION["Mercator_1SP"],PARAMETER["central_meridian",0],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0]]<BR>
ESRI'ified WKT[epsg:3785] = <BR>
PROJCS["Popular_Visualisation_CRS_Mercator_deprecated",<BR>
GEOGCS["GCS_Popular Visualisation CRS",<BR>
DATUM["D_",<BR>
SPHEROID["Popular_Visualisation_Sphere",6378137,0]],<BR>
PRIMEM["Greenwich",0],<BR>
UNIT["Degree",0.017453292519943295]],<BR>
PROJECTION["Mercator"],<BR>
PARAMETER["central_meridian",0],<BR>
PARAMETER["false_easting",0],<BR>
PARAMETER["false_northing",0],<BR>
UNIT["Meter",1],<BR>
PARAMETER["standard_parallel_1",0.0]]<BR>
PROJ.4 rendering of [epsg:3785] = +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>
</SPAN></FONT></FONT>
</BODY>
</HTML>