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@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@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_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@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
