¿Cómo puedo enseñar a ogr2ogr sobre una proyección?

Question

He JSON generado a partir de ArcGIS que quiero convertir a GeoJSON en EPSG 4326. Esta falla:

$ ogr2ogr -f GeoJSON -t_srs EPSG:4326 out.geo.json in.arcgis.json OGRGeoJSON
ERROR 6: EPSG PCS/GCS code 102100 not found in EPSG support files.  Is this a valid
EPSG coordinate system?

He encontrado un post del foro que dice 102100 puede ser cambiado a 3857, que es matemáticamente idéntico. Entonces, esto funciona, con el mismo mensaje de error:

$ ogr2ogr -f GeoJSON -s_srs EPSG:3857 -t_srs EPSG:4326 out.geo.json in.arcgis.json OGRGeoJSON 
ERROR 6: EPSG PCS/GCS code 102100 not found in EPSG support files.  Is this a valid
EPSG coordinate system?

Sin embargo, tengo varios archivos en el 102100 sistema de proyección. Me gustaría enseñar ogr2ogr acerca de la 102100 de proyección, por lo que puedo secuencia de comandos de la conversión. Hay una manera de actualizar la configuración para que 102100 es reconocida y tratada de la misma manera 3857?

La actualización de estos datos es de la API de REST de ArgGIS Servidor 9.3.1. Si no especifica una referencia espacial de salida, luego me sale 102100, que ellos definen como el sistema de coordenadas proyectadas "WGS_1984_Web_Mercator_Auxiliary_sphere". Una manera mucho más fácil para lograr mi objetivo es pedir la proyección deseada, con la outSR=4326 parámetro de consulta. La pregunta sigue siendo válida si usted no tiene acceso a los datos originales.

Answer 1

Como 102100 y 3857 son matemáticamente idénticas, creo que se pueden ignorar los errores y continuar con la invalidación de la fuente srs usando el -s_srs bandera.

Usted puede tratar de añadir la definición de 102100 en su$GDAL_DATA$, pero no es realmente recomendable. Si usted realmente quiere hacer de ella, aunque, aquí está el método y el acompañamiento de la advertencia de Frank Warmerdam:

El .los archivos csv son principalmente para EPSG sistemas de coordenadas, por lo que no es realmente apropiado para extender manualmente.

Pero usted puede agregar sus propios diccionarios personalizados de una manera similar a gdal/datos/cubewerx_extra.wkt. Estos archivos de diccionario consisten en líneas con un nombre corto, una coma y luego el WKT definición de la sistema de coordenadas.

Puede hacer referencia a una definición de SetFromUserInput() y los programas que uso esa función (como s_srs, a_srs, t_srs para gdal_translate, ogr2ogr y gdalwarp) utilizando el formulario:

gdal_translate -a_srs DICT:cubewerx_extra.wkt,41001

Espero que esto ayude. Buena suerte!