All,
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.
A couple of example errors
[u4rt9mds@lidarora1 /] ogr2ogr -f "ESRI Shapefile" -t_srs epsg:3785 test.shp
cincinnati_forests.shp
ERROR 1: no system list, errno: 2
Failed to reproject feature 0 (geometry probably out of source or
destination SRS).
ERROR 1: Terminating translation prematurely after failed
translation of layer cincinnati_forests (use -skipfailures to skip errors)
[u4rt9mds@lidarora1 /] gdalwarp -of vrt -t_srs epsg:900913 -r near
cincinnati_bare_earth_01.img test2.vrt
ERROR 1: no system list, errno: 2
ERROR 1: GDALSuggestedWarpOutput() failed because the passed
transformer failed.
Reprojection to any other projection seems to work fine
[u4rt9mds@lidarora1 /] gdalwarp -of vrt -t_srs epsg:26919
/u02/LASdata/Office_of_Americas/Cincinnati_OH/Delivery/rasters/cincinnati_ba
re_earth_01.img test.vrt
Creating output file that is 20023P x 24321L.
Processing input file
/u02/LASdata/Office_of_Americas/Cincinnati_OH/Delivery/rasters/cincinnati_ba
re_earth_01.img.
Using internal nodata values (eg. -3.40282e+38) for image
/u02/LASdata/Office_of_Americas/Cincinnati_OH/Delivery/rasters/cincinnati_ba
re_earth_01.img.
[u4rt9mds@lidarora1 /] ogr2ogr -progress -f "ESRI Shapefile" -t_srs
epsg:26917 test.shp cincinnati_forests.shp
0...10...20...30...40...50...60...70...80...90...100 - done.
Testepsg reports the correct output
[u4rt9mds@lidarora1 shapes]$ testepsg epsg:3785
Validate Succeeds.
WKT[epsg:3785] =
PROJCS["Popular Visualisation CRS / Mercator (deprecated)",
GEOGCS["Popular Visualisation CRS",
DATUM["Popular_Visualisation_Datum",
SPHEROID["Popular Visualisation Sphere",6378137,0,
AUTHORITY["EPSG","7059"]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY["EPSG","6055"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4055"]],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
PROJECTION["Mercator_1SP"],
PARAMETER["central_meridian",0],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
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"],
AUTHORITY["EPSG","3785"],
AXIS["X",EAST],
AXIS["Y",NORTH]]
Simplified WKT[epsg:3785] =
PROJCS["Popular Visualisation CRS / Mercator (deprecated)",
GEOGCS["Popular Visualisation CRS",
DATUM["Popular_Visualisation_Datum",
SPHEROID["Popular Visualisation Sphere",6378137,0],
TOWGS84[0,0,0,0,0,0,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]]
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_no
rthing",0]]
ESRI'ified WKT[epsg:3785] =
PROJCS["Popular_Visualisation_CRS_Mercator_deprecated",
GEOGCS["GCS_Popular Visualisation CRS",
DATUM["D_",
SPHEROID["Popular_Visualisation_Sphere",6378137,0]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.017453292519943295]],
PROJECTION["Mercator"],
PARAMETER["central_meridian",0],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["Meter",1],
PARAMETER["standard_parallel_1",0.0]]
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
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
