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

[Even Rouault]

Le vendredi 01 juillet 2011 16:16:15, Smith, Michael ERDC-CRREL-NH a écrit :
> 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

--> this error comes from proj.4, so I'd suspect something is wrong with the 
setup of the vertical datum data files. Perhaps activating PROJ_DEBUG 
environmenet variable might give more info.

> 
> 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_b
> a 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_b
> a re_earth_01.img.
> Using internal nodata values (eg. -3.40282e+38) for image
> /u02/LASdata/Office_of_Americas/Cincinnati_OH/Delivery/rasters/cincinnati_b
> a 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