[gdal-dev] Reprojection shapefile using ogr2ogr from EPSG:28992 to EPSG:4326

[Pierric de Laborie]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20150122/4d74eec8/attachment.html>