Identificar el sistema de coordenadas de Shapefile cuando se desconoce

Tengo un Shapefile pero su sistema de coordenadas es Desconocido y no hay un archivo * .prj.

¿Cómo puedo identificarlo ahora?

Siempre existe el método de "fuerza bruta":

1. Tome una capa con un sistema de coordenadas conocido que se supone que se superpone con su capa desconocida.

2. Ahora haga algunas conjeturas sobre qué proyección podría ser la capa desconocida. (UTM, Placa Carree, etc.). Proyecte su capa de sistema de coordenadas conocida en cada proyección hasta que encuentre una que coincida con la capa desconocida tanto como sea posible.

Estoy promoviendo los comentarios de Mapperz y Brandon Copeland agregando una respuesta que utiliza su técnica.

1. Si tiene ArcGIS Desktop, agregue algunos datos de referencia que cubran la misma área a ArcMap. Estos datos de referencia deben tener una definición de sistema de coordenadas (proyección) válida.
2. Agregue los datos en un sistema de coordenadas desconocido.
3. Investigue qué posibles sistemas de coordenadas se utilizan en esta área. Un lugar para verificar es el Registro EPSG .
4. Establezca el sistema de coordenadas del marco de datos a las posibilidades. Si los datos con un coordsys conocido se alinean con los datos sin un coordsys conocido, probablemente haya identificado el sistema de coordenadas.

Esta técnica aprovecha el hecho de que ArcMap puede convertir datos en un nuevo sistema de coordenadas en memoria. Los datos sin un coordsys conocido no se pueden convertir, por lo que solo se muestran. Al proyectar los datos conocidos en la memoria, podemos ver rápidamente qué coordenadas tiene en los diversos sistemas de coordenadas posibles.

Como esta pregunta nunca pasa de moda, construí un sitio que utiliza el método de la fuerza bruta . Si arrastra un shp + shx comprimido al mapa, lo asignará en todos los sistemas de coordenadas disponibles en PostGIS. Suponiendo que sabe cómo se ve "correcto", puede hacer zoom en esa área y hacer clic en el polígono para obtener el archivo .prj de epsg.io.

Hay dos excelentes enlaces de Esri que entran en detalles sobre esto:

• Identifique un sistema de coordenadas proyectadas desconocido utilizando ArcMap
• Identificar la referencia espacial, proyección o sistema de coordenadas de datos.

En mi último trabajo, tenía una capa de datos de Geología (polígono) llamada "FSU_Geol.shp". Mi jefe me lo dio y me pidió que descubriera varias cosas. En primer lugar, el cliente le entregó este archivo de forma, y ​​no había un archivo .prj, por lo que quería que lo resolviera. También quería que descubriera cuáles eran las categorías de Geología. Podría seguir y seguir, pero vamos al grano ... Lo busqué en Google y terminé AQUÍ .

El acrónimo en la convención de nombres significaba "Geología de superficie de la antigua Unión Soviética", y la búsqueda en Google me llevó "directamente" a la fuente (USGS). Todo lo que podría necesitar saber sobre este archivo de forma estaba en el enlace superior que golpeé. No estoy diciendo que Google pueda encontrar nada y todo , pero estaba recién salido de la universidad y solo estaba tomando una "foto en la oscuridad", ¡y mira los comentarios que recibí!

En otro caso, alguien en la oficina había descargado una serie de archivos de forma usando una función por lotes. No recuerdo el nombre del paquete de software en este momento, pero faltaban los archivos .prj. Simplemente ingresé a los metadatos en ArcCatalog, y realmente encontré mi respuesta allí. Nuevamente, no me cite, pero creo que había una URL a la fuente allí.

¡Desde ese día "comienzo" con una búsqueda en Google si no hay pistas en los metadatos!

Pregúntale al productor.

Si conoce CRS comunes para su región geográfica, puede probar algunos de ellos. Pero preguntar es mejor.

Realmente amo esta herramienta web: http://projfinder.com/ . Revise su archivo y busque algunas coordenadas. Zoom a aproximadamente el mismo lugar en el mundo y dejarlo adivinar.

No es una herramienta (no sé de una existente que te permita hacer eso), pero mira la respuesta de @ mkennedy a ¿Cómo convertiría este punto a WKID 4326? . Ella explica cómo llegó a la referencia espacial correcta. SpatialReference.org y paciencia serán tus amigos.

Además, Esri proporciona una guía sobre cómo adivinar un sistema de coordenadas (aunque prefiero el método de mkennedy si sabes un poco más sobre los datos).

La mejor respuesta que he encontrado para esta pregunta no es técnica: averigüe de dónde provienen sus datos. Las agencias y organizaciones tienden a ser consistentes con su uso de proyecciones. ¿Sabes que vino de tu estado DOT? Mire el resto de sus datos y vea lo que le dice. ¿No sabes de dónde vino? Es muy probable que una suposición educada lo envíe por el camino correcto.

¡Al menos hace que abordar el problema con la fuerza bruta sea un poco más factible!

Seriamente desactualizado en este punto, pero Werner Flacke y Birgit Klaus publicaron Find Projection en ArcScripts en 2007. No creo que el código fuente esté allí, desafortunadamente. Está basado en VBA, por lo que solo se puede usar en ArcGIS Desktop v9.2 y posiblemente 9.3. Incluye dos archivos de forma con las áreas de interés del conjunto de datos de parámetros geodésicos EPSG que podrían utilizarse para reducir las posibilidades.

Blue Marble Geographic Calculator y Geographic Transformer tienen herramientas de recuperación del sistema de coordenadas.

Pruebe el programa ogrinfo suministrado como parte de GDAL.

Consulte ¿Cómo acceder a los metadatos de Shapefile usando OGR?

Por lo tanto, es posible que ogrinfo no le brinde la información de proyección sin un archivo .prj, pero sigue siendo una herramienta útil para ayudarlo a investigar la lista de posibles proyecciones.

Por ejemplo:

Geometry: Polygon
Feature Count: 269 
Extent: (320000.000000, 505000.000000) - (323000.000000, 511000.000000) 

Esta respuesta me sugiere que el shapefile está usando un sistema de referencia basado en metros y no en grados.

Suponiendo que sabe de dónde provienen los datos, ahora tiene una lista más corta de posibles proyecciones.

Otras mejores prácticas pueden ser:

• Pregunte a las personas que le dieron los datos en primer lugar, incluso si no saben exactamente, deberían tener una lista de posibles proyecciones que podrían analizar en primera instancia.
• Mire para ver si hay algún metadato asociado.

Probablemente sea mejor comenzar desde alguna posición de conocimiento, en lugar de intentar y la fuerza bruta desde cero.

Para agregar a las otras respuestas en esta publicación, agregaría lo siguiente:

Cómo: Identificar un sistema de coordenadas proyectadas desconocido utilizando ArcMap

Proporcionar una verificación específica en NAD1927

Si las coordenadas están en grados decimales, como entre longitud -180 y +180, y latitud -90 y +90, identifique el Sistema de coordenadas geográficas (dato) utilizado para los datos. Antes de la versión 9.2, ArcMap asigna GCS_Assumed_Geographic_1 a los datos, de manera predeterminada. Esto coloca los datos en el dato NAD_1927

Adicionalmente

Si los datos se encuentran en los Estados Unidos y muestran una extensión en la que las coordenadas a la izquierda del decimal son 6, 7 u 8 dígitos, los datos probablemente se proyecten en los sistemas de coordenadas del plano estatal o UTM.

Preguntas frecuentes: Conceptos básicos de proyección: lo que el profesional de SIG necesita saber

Esta página proporciona varias verificaciones para determinar las proyecciones, lo cual es importante para ayudar a tratar de identificar un sistema de coordenadas desconocido.

Los siguientes conceptos son fundamentales para comprender el uso de las proyecciones de mapas en ArcGIS. Sin embargo, tenga en cuenta que el tema de las proyecciones es extremadamente amplio, y este artículo no puede hacer más que tocar algunos temas importantes.

1. Los sistemas de coordenadas, también conocidos como proyecciones de mapas, son designaciones arbitrarias para datos espaciales. Su propósito es proporcionar una base común para la comunicación sobre un lugar o área particular en la superficie de la tierra. El problema más crítico al tratar con las proyecciones de mapas es saber cuál es la proyección y tener la información correcta del sistema de coordenadas asociado con un conjunto de datos.

2. Cuando se idearon las primeras proyecciones del mapa, se asumió, incorrectamente, que la tierra era plana. Más tarde, la suposición fue revisada, y se supuso que la tierra era una esfera perfecta. En el siglo XVIII, la gente comenzó a darse cuenta de que la tierra no era perfectamente redonda. Este fue el comienzo del concepto del esferoide cartográfico.

3. Para representar con mayor precisión las ubicaciones en la superficie de la tierra, los creadores de mapas estudiaron la forma de la tierra (geodesia) y crearon el concepto del esferoide. Luego se diseñaron los sistemas de coordenadas geográficas (GCS), que incluyen un dato, unidades de medida y un meridiano principal. Un dato vincula un esferoide a una porción particular de la superficie de la tierra. Los datos recientes están diseñados para adaptarse bien a toda la superficie terrestre.

4. Los datos más utilizados en América del Norte son:

• NAD 1927 (North American Datum 1927) con el esferoide Clarke 1866
• NAD 1983 (North American Datum 1983) con el esfero GRS 1980
• WGS 1984 (World Geodetic Survey 1984) con el esferoide WGS 1984

Los esferoides más nuevos se desarrollan a partir de mediciones satelitales y son más precisos que los desarrollados por Clarke en 1866. Los términos 'sistema de coordenadas geográficas' y 'dato' se usan indistintamente, pero como se señaló anteriormente, un GCS incluye un dato, esferoide, unidades de medida y un primer meridiano.

5. Las coordenadas para los datos cambian según el datum y el esferoide en el que se basan esas coordenadas, incluso si están utilizando la misma proyección de mapa y parámetros.

Por ejemplo, las coordenadas geográficas a continuación son para un solo punto ubicado dentro de la ciudad de Bellingham, Washington, usando 3 datos diferentes:

Código: DATUM X-Coordinate Y-Coordinate NAD_1927 -122.466903686523 48.7440490722656 NAD_1983 -122.46818353793 48.7438798543649 WGS_1984 -122.46818353793 48.7438798534299

6. Un principio de una buena gestión de datos es obtener los parámetros de proyección de la fuente de datos que proporciona los datos. No haga una conjetura sobre la proyección de datos, porque el resultado será una base de datos SIG inexacta. Los parámetros necesarios son los siguientes:

• Proyección
• Unidades de medida
• ZONA (para UTM)
• Zona FIPS (para plano de estado)
• Datum

Se pueden requerir otros parámetros, dependiendo de la proyección. Por ejemplo, las proyecciones de Albers y Lambert requieren los siguientes parámetros:

• 1er paralelo estándar, en grados, minutos y segundos (DMS)
• 2do paralelo estándar (DMS)
• Meridiano central (DMS)
• Latitud del origen de las proyecciones (DMS)
• Falso este y unidades de medida
• Falso norte y unidades de medida
• Desplazamiento X y unidades de medida
• Desplazamiento Y y unidades de medida

Tuve el mismo problema mientras trabajaba en mi proyecto (DEM). Una de las DEM era coordenadas desconocidas y no se puede proyectar, así que hice la 'georreferenciación' utilizando la herramienta Georreferenciación en ArcMap . Después de Georeferenciar, obtendrá su archivo prj.

Mira este sitio web:

http://www.egger-gis.at/shapefile-projectionfinder/

Soy el desarrollador de esta herramienta gratuita. Tal vez esta herramienta pueda ayudarlo a encontrar y definir la proyección correcta de su shapefile. Se basa en la idea de http://projfinder.com/ por Aaron Racicot.

También puede probar estas soluciones especiales para:

Austria (alemán): https://www.data.gv.at/anwendungen/checkaustrianprojection/

Australia: https://maegger.github.io/map_australia.html

UTM - Zonas: https://maegger.github.io/map_utm.html

Puede discriminar entre coordenadas geográficas y de cartografía, el primer lat / lon está en grado y el segundo norte / este está en metros. Si el archivo tiene forma, puede volcar las coordenadas de la geometría con shapelib-1.3.0 y ver los valores. ..

Puedes hacerlo con Fiona.

    import fiona as f
    a = fiona.open("C:\QGIS_ShapeFile\qgis\shafile_XXX.shp")
    print(a.crs)

Obtendrá respuesta como

    {'init': 'epsg:4326'}