Identificar el sistema de coordenadas de Shapefile cuando se desconoce

Respuestas:

49

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.

brenth
fuente
1
Este es mi método también. Utilizo dos programas diferentes para "probarlo". También puede leer esta pregunta para obtener una orientación básica sobre las proyecciones. gis.stackexchange.com/questions/2769/…
Brad Nesom el
66
Puede simplificar un poco el paso 2 simplemente cambiando la proyección del marco de datos (ArcMap) y dejando que la capa conocida se vuelva a proyectar sobre la marcha.
Brandon Copeland
13
Tenga en cuenta que no debe proyectar los datos. Deberías estar definiendo la proyección. Proyectar los datos cambia las coordenadas. Definir la proyección solo asigna una proyección a esos datos y le dice al software en qué proyección están los datos. Estas son dos herramientas diferentes en ArcMap, "Definir proyección" es la que desea usar.
Sean
Prefiero leer los valores de coordenadas, generalmente puedo decir si son grados decimales, metros o pies (más difícil de distinguir) y luego pasar al mejor sistema. Comience con WGS 84 si aparece decimal, UTM si metros, Estado plano si pies. Una vez tuve uno que nadie podía arreglar (por anomalías gravitacionales). Era un avión del estado de EE. UU. Con una unidad lineal de ...... kilómetros. Me di cuenta por el hecho de que estaba cerca, tenía el sistema correcto pero el estado era pequeño, ahh kilómetros era la unidad lineal.
Si no lo sabes, solo GIS
33

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.

mkennedy
fuente
44
Nota al margen: en QGIS, la reproyección en memoria se denomina "sobre la marcha".
oscuro
3
¡También se llama 'sobre la marcha' en ArcGIS! Sin embargo, sigo encontrando muchas personas que no entienden esa frase, así que he estado tratando de usar 'en memoria'.
mkennedy
¿Sigue siendo válido este enlace? Parece un poco más complicado de lo que recuerdo, sin embargo, existe este epsg-registry.org
GISKid
1
@GISKid Tienes razón. Creo que el sitio web epsg-registry.org es posterior a la fecha en que escribí esta respuesta. Voy a actualizar
mkennedy
1
Otro gran recurso: epsg.io
GISKid
29

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.

jmk
fuente
15
idea excelente, pero creo que la página web no está realmente trabajando
Fábio Dias
22

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!

Dano
fuente
18

Pregúntale al productor.

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

bajo oscuro
fuente
El único problema es que no puedo preguntarle ... Este fue también el primer pensamiento que tuve.
com
Entonces, tendrás que usar la fuerza bruta como dijo Brenth.
oscuro
18

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.

johanvdw
fuente
Esta es en realidad la mejor herramienta que he encontrado. Las conjeturas en ArcMap no son realmente útiles a largo plazo o si no sabes cuál debería ser. No puedo ayudar con los sistemas de coordenadas homebrew y la falta de uno por completo, ¡pero quién puede!
Martin
Lástima que ya no funciona, ¿alguien sabe a quién contactar?
sventechie
13

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!

Jason Baker
fuente
11

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.

mkennedy
fuente
Gracias por el dato Melita. El código está ahí (* .frm). La licencia es de dominio público. PrjFinder es un compañero del libro "Trabajar con proyecciones y transformaciones de datos en ArcGIS" (2006).
Matt Wilkie
4

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.

nmtoken
fuente
1
Esto no responde la pregunta que se hizo.
Dan C
1
Si no hay un archivo prj, ogrinfo no devolverá resultados para la información del sistema de coordenadas.
Devdatta Tengshe
4

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.

  1. 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

  1. 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

whyzar
fuente
¿Qué debo hacer si tengo una capa de cuadrícula que solo muestra un sistema de coordenadas desconocido y no puedo asignarle nada? Gracias.
AndrewLebron
2

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.

ndr
fuente
0

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

Manfred Egger
fuente
0

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. ..

Natan
fuente
0

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'}
Sandip Nepal
fuente
¿Estás seguro de que funcionará sin un archivo .prj? ¿Cómo lo sabrá Fiona?
JR
No sin el archivo .prj, no funcionará
Sandip Nepal