[gdal-dev] GDAL/OGR projection problem only with epsg:3785/900913

Smith, Michael ERDC-CRREL-NH michael.smith at usace.army.mil
Fri Jul 1 10:16:15 EDT 2011


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 at 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 at 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 at 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 at 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 at 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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20110701/98ac1071/attachment.html


More information about the gdal-dev mailing list