Cuando se trabaja con datos espaciales en software SIG, es habitual encontrarse con el problema de archivos que no tienen un sistema de coordenadas claramente definido.
Cuando un archivo carece de su sistema de coordenadas (SCR), puede dar lugar a capas desalineadas, mediciones incorrectas y otros muchos problemas que dificultan la visualización y el análisis precisos de los datos.
En este artículo intentaremos darte los consejos y herramientas disponibles en 2023 para identificar sistemas de coordenadas de archivos con proyecciones desconocidas.
Antes de abordar el problema, pongamos las cartas sobre la mesa: usted puede ser parcial o totalmente responsable del problema. Habrá dos capítulos en este artículo: el primero para los que no estén acostumbrados a los sistemas de coordenadas, el segundo para los que ya hayan hecho por sí mismos las comprobaciones del primer capítulo y que, por tanto, pueden pasar directamente al segundo.
Recuerda que las soluciones sugeridas aquí nunca resolverán el 100% de los problemas…
Eres principiante en sistemas de coordenadas
El problema puede surgir en dos ocasiones
Has creado una nueva capa de datos con el procesamiento espacial de QGis y el resultado no se solapa con otros datos.
Intenta cargar nuevos datos y éstos no se solapan en absoluto con los demás datos del proyecto.
La capa es el resultado de un procesamiento
En este caso, su preocupación es que el programa de procesamiento no asignó, o no pudo asignar, el SCR de su proyecto a su salida.
Para averiguar el SCR del proyecto, vaya al menú Proyecto -> Propiedades, y a la pestaña SCR
El panel inferior le muestra el nombre y la configuración del SCR de su proyecto, pero fíjese especialmente en el código que aparece al final de la línea en la lista de sistemas utilizados recientemente, en este ejemplo EPSG:2154.
Para indicar a QGis que su capa utiliza este mismo sistema, haga clic con el botón derecho del ratón sobre la capa -> SCR de capa
La primera línea, Ningún SCR, corresponde al estado de descripción del SCR conocido por QGis. Para asignarle el SCR del proyecto, seleccione la línea correspondiente al SCR en la ventana anterior.
Si, por el contrario, la primera línea no es Ningún SCR, como en este caso:
Esto significa que el programa de tratamiento ha asignado a su capa un SCR incorrecto. Esto puede ocurrir cuando su proyecto y las capas que contiene no están en el mismo sistema de coordenadas. QGis sabe cómo manejar esto, pero el procesamiento tiene varios orígenes y puede tener errores.
Mire el SCR de la capa de entrada. Si no es el mismo que el del proyecto, la solución suele ser cambiar el SCR de la capa problemática por el de la otra.
La capa procede de datos externos que no conoce.
Compruebe los metadatos y los archivos que los acompañan. En la mayoría de los casos, los archivos están en formato shapefile. Un shapefile consta de varios archivos obligatorios y opcionales, a saber
- .shp: el archivo principal que contiene los datos geométricos. (obligatorio)
- .shx: el archivo índice para acceder rápidamente a los datos geométricos. (obligatorio)
- .dbf: La tabla dBASE que contiene los datos de atributos (obligatorio)
- .prj: El fichero de proyección que contiene la información del sistema de coordenadas (si está disponible)
En general, el problema de la definición del SCR viene de la ausencia del fichero .prj. Si este fichero falta o está vacío, tendrá que identificar usted mismo el sistema de coordenadas. Si el fichero está presente y cumplimentado, puede abrirlo con un editor de texto para visualizar la descripción del sistema de coordenadas.
He aquí un ejemplo de un fichero .prj completo:
Si falta o está vacío, tendrá que fijarlo en las propiedades de la capa. Pero para ello tendrá que averiguar de qué sistema se trata: busque archivos opcionales, como . txt o . xml, que puedan contener metadatos o información adicional sobre los datos; a veces los shapefiles van acompañados de archivos de documentación, como un README.txt o un informe en PDF, que proporcionan información.
Si no encuentra la información que necesita en los metadatos o en los archivos relacionados, póngase en contacto con el proveedor o la fuente de los datos. Es posible que ellos puedan proporcionarle la información necesaria sobre el sistema de coordenadas.
Si encuentra la respuesta, deberá indicar a QGis qué SIR debe utilizar. Haga clic en la capa, menú Capa SCR ->Capa SCR :
Se abre la ventana de definición del SCR.
Utilizando el campo de filtro, introduzca partes del nombre del SCR buscado. Obtendrá una lista de SCR que contienen esta parte en la lista Sistemas de coordenadas de referencia predefinidos.
Si a estas alturas no ha resuelto su problema, tendrá que armarse de valor y pasar al capítulo siguiente.
Sabes que el SCR de los datos es desconocido y quieres averiguarlo.
En 2015 había publicado un documento «Determinación de la proyección desconocida de una capa para QGis 2.X», en francés, que sigue disponible pinchando aquí. Solo es válido para las proyecciones actuales de Francia Metropolitana y se basa en el análisis de los valores XY de los datos desconocidos. Afortunadamente, QGis y sus desarrolladores han avanzado mucho desde entonces, ya que QGis 3.x ofrece una herramienta para encontrar un sistema desconocido, en los tratamientos -> Herramientas vectoriales generales -> Encontrar proyección.
La herramienta de búsqueda de proyecciones en QGIS está diseñada para crear una lista de sistemas de referencia de coordenadas (SRC) candidatos para una capa con una proyección desconocida, basada en un área objetivo especificada utilizando otra capa de referencia, con un SRC bien definido.
Por ejemplo, si sabe o cree que sus datos se encuentran en algún lugar alrededor de Arcachon, cargue una capa de datos que incluya las Landes :
- Límites administrativos (por ejemplo, fronteras nacionales, regionales o departamentales)
- Redes de carreteras
- Imágenes por satélite o fotografía aérea (por ejemplo, Google Maps, OpenStreetMap)
- etc.
En este ejemplo tenemos una capa nbm con un SCR desconocido, pero sabemos que está relacionada con la cuenca de Arcachon.
La herramienta Buscar Proyección tomará los datos nbm e iterará a través de todos los SCRs conocidos: Para cada SCR, lo asignará a la capa desconocida, y comprobará si los datos así rellenados están dentro o cerca del área objetivo definida por la capa de referencia. En caso afirmativo, el tratamiento añade el SCR a la lista de posibles candidatos. El proceso continúa recorriendo todos los SCR hasta crear la lista completa de candidatos.
El algoritmo genera así una lista corta de SCR candidatos que, si se utilizan para la capa desconocida, la situarían cerca de la zona objetivo. El usuario deberá probar manualmente cada SCR candidato utilizando QGis, definiendo el SCR para la capa y comparándolo con la capa de referencia con un SCR conocido, basándose en la alineación y la compatibilidad con la capa de referencia.
Debe tenerse en cuenta que el algoritmo puede producir varios candidatos o ningún candidato, dependiendo de la precisión de la entrada del área objetivo y de la extensión de la capa. Cuanto mayor sea la extensión de la capa de referencia, mayor será el número de SCR candidatos. También es posible que el SCR correcto no esté entre los candidatos preseleccionados, por lo que es crucial comprobar los resultados y combinar este enfoque con otros métodos para garantizar una evaluación completa.
Para continuar con nuestro ejemplo, abrimos la herramienta Buscar proyección:
Pour poursuivre notre exemple, nous ouvrons l’outil Trouver une projection:
Puede definir el Área de destino de la capa haciendo clic en el icono situado a la derecha del campo. En este caso el área será la que se muestre en la ventana del mapa QGis. También puede abrir la lista desplegable que le ofrece otros métodos de definición. Haga clic en Ejecutar y espere a que QGis proyecte su capa desconocida en cada uno de los SCR que conoce.
Cuando termine, mostrará la preselección en la pestaña de registro, y habrá añadido una tabla en el panel de capas con un registro por SCR candidato.
Para ponértelo más difícil, y QGis no tiene nada que ver con esto, el mismo sistema de coordenadas tendrá diferentes códigos dependiendo de la organización de referencia. Son las 4 letras que aparecen al principio de cada código. EPSG significa European Petroleum Survey Group, la primera organización que estandarizó los sistemas de coordenadas en todo el mundo, IGNF significa IGN Francia, y ESRI por ESRI, por supuesto. Cada uno tiene un código diferente para el mismo sistema.
Si nos fijamos en las definiciones de todos los SIR en este ejemplo:
EPSG:27573
EPSG:27583
IGNF:LAMB3C
ESRI:102583
EPSG:7403
EPSG:7413
EPSG:7422
ESRI:102587
IGNF:NTFLAMB3C
IGNF:NTFLAMB3C.IGN69
Todos ellos corresponden al mismo sistema de proyección zonal Lambert 3, utilizado hasta la década de 2000 para el sur de Francia. Sea cual sea el SCR utilizado, el resultado será el mismo:
Pero no siempre es tan simple. Si tienes datos no proyectados, en coordenadas geográficas Latitud-Longitud tendrás el problema de elegir DATUM, en términos menos complicados, el centro de la tierra de referencia. Las mismas coordenadas Lat/Long, según el DATUM utilizado pueden encontrarse en posiciones distantes de 350 metros…
En resumen, si una vez seguidos todos los pasos anteriores no has podido resolver tu problema, te queda la opción de contactarme a través del formulario de contacto. Trata de describir tu problema lo mejor posible, nos ahorrará tiempo a los dos 😉.
