Hello, I do not understand why there seems to be a drift (about 100 meters) of the polygons contained in a shapefile after reprojecting it to 4326 using the following command :
ogr2ogr -f 'ESRI Shapefile' -wrapdateline -t_srs EPSG:4326 destinationfile.shp sourcefile.shp The original prj file is the following. It is supposed to describe the projection EPSG:28992 ( http://spatialreference.org/ref/epsg/28992/ ) PROJCS["RD_New",GEOGCS["GCS_Amersfoort",DATUM["D_Amersfoort", SPHEROID["Bessel_1841",6377397.155,299.1528128]], PRIMEM["Greenwich",0.0], UNIT["Degree",0.0174532925199433]], PROJECTION["Double_Stereographic"], PARAMETER["False_Easting",155000.0], PARAMETER["False_Northing",463000.0], PARAMETER["Central_Meridian",5.38763888888889], PARAMETER["Scale_Factor",0.9999079], PARAMETER["Latitude_Of_Origin",52.15616055555555], UNIT["Meter",1.0]] Apparently QGIS could understand the projection file correctly and reproject it to WGS84 with no problem. But it seems that for some specific projections, GDAL does not reproject as expected. Is it possible to force it by adding the projection specifications in one of the files included in the following folder (linux debian install) : /usr/share/gdal and how ? For example, I could see a file named gdal_datum.csv containing the datum "D_Amersfoort" specified in the prj file. However I could not find in any other files the projection GCS_Amersfoort or RD_New as specified in the prj file. Thanks for your help
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
