Patricio Soriano Castro

En numerosas ocasiones el resultado de la toma de datos en campo con un GPS, una estación total o capturando datos sobre bases cartográficas es un listado de coordenadas que identifican una serie de localizaciones, como podrían ser elementos arqueológicos, un elemento poligonal como un entorno de protección o un corte, o una estructura lineal tipo infraestructura o camino.

Visor Zona Arqueológica de Ategua.BIC y entorno de protección. Ver mapa más grande

Este tipo de información geográfica, me refiero al listado de coordenadas, suele ser solicitado por administraciones públicas como parte de los informes finales de proyectos y pueden encontrarse en boletines y publicaciones de la administración. Saber extraer esta información de nuestras entidades geográficas digitales o convertir estos datos en información que pueda ser tratada con una herramienta informática, como un SIG, puede ser un aprendizaje interesante y de gran utilidad desde el punto de vista profesional.

El Sistema de Información Geográfica QGIS, posee herramientas directas para convertir listados de coordenadas en formato tabla (CSV o TXT) en entidades geográficas puntuales. Si en las tablas existe otra columna que nos sirva para agrupar los datos,por ejemplo un código de yacimiento, estos datos pueden pasar a convertir de forma automática en líneas y estas a su vez en polígonos.

ATEGUA. Banco imágenes IAPH. Autor: J.C. Cazalla

Vista panorámica de Ategua. Fuente: Banco de imágenes IAPH. Autor: Cazalla, Juan Carlos

Supongamos que tenemos un listado de coordenadas que identifican un bien arqueológico y también su entorno de protección. En este caso vamos a obtener los datos del yacimiento arqueológico de Ategua en la provincia de Córdoba. Los datos de las coordenadas del BIC y su entorno están disponibles en el Boletín Oficial de la Junta de Andalucía (BOJA)  nº 244 de 16/12/2005 . Lamentablemente  estos datos no están en el PDF en formato imagen por lo que tendremos que ir copiando manualmente las coordenadas en una hoja de cálculo en la que se incluya  el número de orden la coordenada, el valor de X y el valor de Y. Como tenemos la información tanto del BIC como del entorno , vamos a añadir una nueva columna al principio que denominaremos entidad, donde quedará reflejado si la coordenadas es de la Zona Arqueológica o del Entorno.

Una vez obtenidos los datos, pasaremos la tabla a formato CSV, utilizando las opción de Guardar como de nuestro programa de hoja de cálculo, en mi caso Calc de LibreOffice. El formato csv no es más que un texto plano donde cada fila es un registro de la tabla, y en que los valores por columnas se encuentra separado por comas. Es un formato estándar que suele poder cargarse en la mayoría de los programas que trabajan con datos tabulares como las bases de datos o las hojas de cálculo, o en nuestro caso los SIG.

Listado de coordenadas en hoja de cálculo y archivo CSV.

Una vez que tengamos nuestro csv, abrimos QGIS y seleccionamos el Sistema de Coordenadas de Referencia en el que se encuentren los vértices. Aunque en el BOJA, no se hace referencia al SRC, aunque debería, por la fecha estos datos están en ED50 UTM30N.

Para añadir la tabla de coordenadas, usamos la opción “Añadir capa de texto delimitado” del menú Añadir Capa. Indicaremos la ubicación del archivo CSV, indicaremos el tipo y comprobaremos que QGIS a seleccionado las columnas correctas para las coordenadas X e Y. Una vez revisado las opciones, pulsamos en aceptar y tendremos representados las coordenadas de  nuestra tabla en puntos. La capa se genera de forma temporal por lo que si queremos trabajar con ella podemos guardarla en Shape o almacenarla en nuestra base de datos PostGIS o SpatiaLite. Para una mejor compresión de los datos podemos cargar cartografía base o alguna ortofoto como las del PNOA.

Para pasar los datos puntuales a lineales o poligonales y mejorar la compresión de la información, vamos a utilizar dos procesos de la caja de Herramientas.

  • Point to path nos generará una línea según el orden de los puntos (Campo de Orden) y utilizará los datos de la columna entidad para agrupar los trazados (Campo de Grupo).
  • Con Lines to Polygon generaremos una nueva capa poligonal a partir de las líneas. Tendremos que tener en cuenta que nos generará un polígono para la delimitación de la Zona Arqueológica y otro para el Entorno de Protección  por lo que una vez ejecutado será interesado obtener dos capas poligonales con cada una de las entidades.

Este manual forma parte del temario del curso  on line “Sistemas de Información Geográfica y Arqueología” organizado por Almagre. Para más información puede consultarse la página del curso en la web del curso.

Este artículo continua la entrada anterior en la se expusieron algunos de los usos de los comandos gdalinfo y gdaltranslate de la librería GDAL. En el siguiente texto nos vamos a centrar en el comando gdalwarp. Este comando se usa para hacer reproyecciones del Sistema de Referencia de Coordenadas de una imagen georeferenciada.

Reproyección de una imagen georeferenada con gdalwarp

El uso básico es sencillo: partimos de la base de que la imagen ya posee un SRC definido, tras el comando indicaremos el SRC de salida usando el código EPSG (ej. -t_srs "EPSG:4326"), el fichero a reproyectar y el nombre la imagen proyectada.

$ gdalwarp  -t_srs "EPSG:4326"  -of GTiff img.tif img4326.tif

Las opciones de gdalwarp son muchos más variadas establecer un valor sin datos (- dstnodata) o usar una archivo vectorial para recortar la imagen de salida (-cutline). Todas las opciones pueden consultarse en la siguiente dirección http://www.gdal.org/gdalwarp.html

De ED50 a ETRS89

En España, a partir del 1 de enero de 2015  y según  el Real Decreto 1071/2007 , de 27 de julio, por el que se regula el sistema geodésico de referencia oficial en España, "...toda la cartografía y bases de datos de información geográfica y cartográfica producida o actualizada por las Administraciones Públicas deberá compilarse y publicarse conforme a lo que se dispone en este real decreto” o lo que es lo mismos debe estar encontrarse en el ETRS89. Más información en esta página del IGN.

Para facilitar la transformación del datum ED50 a ETRS89, el Instituto Geográfico Nacional ha generado dos rejillas (Península y Baleares) en formato NTV2. Estas rejillas pueden ser descargadas desde la web del IGN.  Las rejillas pueden añadidas  usadas en la mayoría de los SIG pero también podemos utilizarla directamente usando GDAL con el comando gdalwarp. En la web de Mappingis podréis encontrar un artículo de su uso en QGIS.

Si queremos trasformar con gdalwarp una imagen en ED50 UTM30N (EPGS:23030) a ETRS89 (EPGS:25830), lo primero será descargar la correspondiente rejilla a nuestro ordenador. Tras el comando debemos añadir los SRC de entrada (-s_srs) y salida (-t_srs), pero en esta ocasión añadiendo los parámetros de la librería PROJ.4 e indicando la ruta de la rejilla del IGN tras la opción nadgrids.

Suponiendo que hemos almacenado la rejilla de la península en nuestro disco duro C: en la carpeta "rejillas", el comando quedaría así:

$ gdalwarp -s_srs "+init=epsg:23030 +nadgrids=C/:rejilla/PENR2009.gsb +wktext" -t_srs "+init=epsg:25830 +nadgrids=null +wktext" -of GTiff img23030.tif img25830.tif

Un poco de programación: procesando de múltiples imágenes con gdalwarp

Como he comentado al principio, todo este textazo tiene en primer lugar una función didáctica para mi. El segundo aspecto por el que merece utilizar los comando de estas librerías es el hecho de poder crear nuestros propios fragmentos de código o scripts que nos ayuden a mejorar nuestra vida..al menos desde el punto de vista informático. Creo que en esta búsqueda de la felicidad, podría ser interesante hacer un un pequeño script que por ejemplo convirtiera a formato GeoTiff y mejorara el tamaño de, ¿porqué no?, 200 ficheros.

Creamos un script utilizando el shell de Linux con la extensión sh (ej. ed50etrs89.sh). El script está pensado para ser utilizado desde el terminal de Linux. Si alguno se lo trabaja para Windows o Mac, y quiere que lo ponga en la entrada avisad por correo.

El código hace lo siguiente:

  • Localiza las imágenes con la extensión JPG que se encuentran en el mismo directorio del archivo  ed50etrs89.sh
  • Reproyecta  y comprime las imágenes a ETRS89 UTM 30N usando la rejilla del IGN localizada en en mismo directorio del script.
  • Salva las imágenes con el sufijo “_25830”  y en formato GeoTiff
  • Usa gdalinfo para generar una archivo de texto con los metadatos de las imágenes reproyectadas.

 

Si es necesario le asignamos, los correspondientes permisos de ejecución con chmod

$ chmod +x ed50etrs89.sh

Y a continuación, lo ejecutamos

$ ./ed50etrs89.sh

Sobre mí

SIGdeletras es Patricio Soriano y Patricio Soriano es SIGdeletras. Trabajo el campo las Tecnologías de la Información Geográfica y especialmente su aplicación en el ámbito del Administración Pública y el Patrimonio Cultural...  ¿Quieres saber más sobre mí?

 

Buscar