ogr2ogr de NTF_Lambert_II_etendu es de aproximadamente 30 metros

Question

Me estoy convirtiendo RPG archivos en NTF_Lambert_II_etendu a wgs_84 el uso de:

ogr2ogr -t_srs EPSG:4326 ilot_2008_027_wgs84_2.shp ilot_2008_027.shp

pero el resultado es de aproximadamente 30 metros. Sorprendentemente, la misma conversión de uso de Qgis funciona bien. Podría alguien decirme cómo corregir eso o cómo investigar la exacta transformación de qgis está haciendo?

El proj archivo del archivo de salida es :

PROJCS["NTF_Lambert_II_etendu",GEOGCS["GCS_NTF",DATUM["D_NTF",SPHEROID["Clarke_1880_IGN",6378249.2,293.46602]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["False_Easting",600000.0],PARAMETER["False_Northing",2200000.0],PARAMETER["Central_Meridian",2.3372291667],PARAMETER["Standard_Parallel_1",45.8989188889],PARAMETER["Standard_Parallel_2",47.6960144444],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",46.8],UNIT["Meter",1.0]]

Answer 1

NTF Lambert II tiene un towgs84 el dato de cambio, pero que no está incluido en el .prj archivo.

El .prj archivo utiliza una nomenclatura diferente para la proyección, por lo que el código EPSG finder puede fallar.

Supongo QGIS asigna el correcto código EPSG (tal vez en el .qpj archivo), y la convierte en un estándar de la transformación de un código EPSG a otro con plena towgs84 parámetros.

Usted debe obtener el mismo resultado con ogr2ogr el uso de los códigos EPSG 27572, 7411, 7412 o 7421 como -s_srs.

Answer 2

Como se ha mencionado por @AndreJ especificando el código EPSG facilita el trabajo. En concreto, los siguientes trabajos:

ogr2ogr -t_srs EPSG:4326 -s_srs EPSG:7421 ilot_2008_027_wgs84_2.shp ilot_2008_027.shp

Alternativamente, yo estaba tratando de usar el qgis API de python para hacerlo todo en python usando qgis, pero no podía llegar a trabajar.