miércoles, 30 de septiembre de 2020

Asignar coordenadas y alturas a puntos en QGIS.

Hoy vamos a describir un par de operaciones básicas que podemos necesitar realizar sobre una capa de puntos vectorial. Esta capa es posible que se encuentre desnuda de atributos y podamos enriquecer sus datos con unos campos tan importantes como las coordenadas geográficas y la altitud de esos puntos.
Partimos de un ejemplo práctico en el que tenemos cargado en QGIS un archivo de relieve sobre cuya extensión generaremos un determinado número de puntos de forma aleatoria. Usamos para ello la herramienta Random points in extent (podemos buscarla en la cajita inferior izquierda del programa):
Rellenamos el Input con la extensión de nuestra capa MDT.
Establecemos en 30 el número de puntos a crear a lo largo de la extensión de la capa y ejecutamos el algoritmo.
Tenemos nuestra nueva capa temporal de puntos llamada Random points con el único atributo de su valor numérico correlativo ID.

El siguiente paso será darle los valores de altitud a cada punto, proveniente de nuestra capa de relieve MDT, utilizando el complemento Point Sampling Tool.
Con nuestra capa de puntos seleccionada en el complemento marcamos el MDT como fuente de datos para los valores y le damos al archivo una ubicación y un nombre (en el ejemplo puntos Z.shp). En la pestaña Fields podemos personalizar el nombre del campo atributo para el dato de altitud (en el ejemplo Z_MDT25).

Ejecutamos el complemento y la capa puntos Z.shp se añade a nuestro proyecto. El siguiente paso será darle unas coordenadas a cada punto para lo que utilizaremos el algoritmo Saga Add coordinates to points

Una nueva capa temporal Points with coordinates se añade a nuestro proyecto. En este caso ya conteniendo todos los atributos que buscábamos: la altitud y las coordenadas de cada punto.

Solamente nos restaría guardar nuestra capa en el formato que deseemos antes de cerrar el proyecto.

PD: en este caso hemos usado un random points para generar una capa de puntos de ejemplo sobre la que trabajar, pero obviamente podemos partir de una capa de puntos ya dada.
PD2: los algoritmos disponibles en QGIS son muchos y a veces dependen de los proveedores instalados, por lo que hay otras formas de insertar datos en los puntos. Ésta es solamente una.

miércoles, 9 de septiembre de 2020

Mosaicos nacionales históricos Landsat.

Disponible nuevo producto en el Centro de Descargas del Centro Nacional de Información Geográfica (CNIG). Bajo la categoría Ortofotos e Imágenes de Satélite se ofrece la serie de mosaicos nacionales históricos de ortoimágenes por satélite Landsat 5 y 8, a 30 y 15 metros de resolución de pixel respectivamente. La unidad de descarga es un mosaico para la Península, Islas Baleares, Ceuta y Melilla, y otro mosaico para las Islas Canarias. Desde 1986 se generan cada década dos modalidades: color natural y falso color infrarrojo.
El Sistema Geográfico de Referencia (SGR) es ETRS89 en la Península, Islas Baleares, Ceuta y Melilla, y WGS84 en las Islas Canarias. Proyección UTM en su huso correspondiente.

Las imágenes empleadas en los mosaicos son las procesadas por el United States Geological Survey (USGS) a nivel L1TP. Por lo tanto, tienen corrección geométrica pero no corrección atmosférica. Cada mosaico Landsat 5 se presenta en falso color infrarrojo (bandas 4-3-2), y una composición de color natural y pseudocolor (bandas 4, 3, 2 y 1). Cada mosaico Landsat 8 se presenta en falso color infrarrojo (bandas 8, 6, 5 y 4), y una composición de color natural y pseudocolor (bandas 8, 4, 3 y 2). Los datos están corregidos de nubosidad y la periodicidad de la serie es de un producto por década.
Listado completo de ficheros disponibles a fecha de este artículo. Quizá se echa de menos especificar si se trata de imagen Landsat 5 u 8, cosa que no sabemos hasta comenzar la descarga del archivo y ver su nomenclatura, aunque por el tamaño de los archivo tampoco es difícil deducir que las imágenes Landsat 8 son las de un peso mucho mayor al tener más resolución (2014).
Los archivos se ofrecen según el IGN en formato COG. Y si os ha pasado lo mismo que a mi que al leerlo me quedé igual de extrañado por no reconocerlo, se trata de una especie de acrónimo de Cloud Optimized GeoTIFF, un formato de imágenes para el procesamiento geoespacial nativo de la nube. Se trata de un archivo GeoTIFF normal, destinado a ser alojado en un servidor de archivos HTTP, con una organización interna que permite flujos de trabajo más eficientes en la nube. Presentan una nomenclatura como la de este ejemplo:
pnt_landsat8_2014_mosaico_islas_canarias_b432_hu28_8bits_COG_90.tif

La información auxiliar del producto nos proporciona la distribución de las teselas Landsat, tanto en formato PDF como SHP.
En el caso de la versión de las teselas en SHP las mismas ofrecen una buena cantidad de atributos que podrían sernos útiles a la hora de trabajar con la teselas territoriales de las huellas de imágenes Landsat.
Un vistazo por último a una de las imágenes de las Islas Canarias, en este caso Landsat 8 de 2014 a 15 metros de resolución y color natural. Debajo detalle a zoom 100% o pixel real sobre la zona de Santa Cruz de Tenerife.
Esto es todo por hoy. Saludos y hasta la próxima entrada.