Nuevo producto IGN: cartografía antigua.

Disponible a descarga el nuevo producto Cartografía antigua con los mapas que se conservan en la Cartoteca del IGN desde el siglo XV hasta el siglo XIX de España y distintas partes del Mundo. Puede realizarse la descarga de los datos a través de la agrupación Documentación geográfica y cartografía antiguas del Centro de Descargas del CNIG.

Se trata de ficheros digitales, imágenes JPG y documentos en PDF, resultado del escaneado de cartografía y atlas antiguos (entre los siglos XV y XIX) pertenecientes a los fondos cartográficos del Instituto Geográfico Nacional. Son mapas, cartas náuticas, vistas topográficas y atlas de múltiples autores, tanto españoles como extranjeros, la mayoría originales y en algunos casos reproducciones facsímiles. Hay mapamundis, mapas de España y de distintos lugares del mundo.

Son ficheros sin georreferenciar con una resolución de 250 ppp para los JPG y cuya unidad de descarga es el documento individual. A fecha de este artículo la colección cuenta con un total de nada más y nada menos que 1752 documentos (1675 en JPG y 77 PDF's). Pensamos que se trata de trasladar los fondos de la Cartoteca al Centro de Descargas, quizá para unificar aún más el material disponible dependiente del mismo organismo.

Ante semejante aluvión de documentos podemos al menos establecer algún pequeño filtro por siglo en la búsqueda. También por tipo de fichero (JPG o PDF) e incluso el habitual buscador genérico que puede delimitarnos la búsqueda inicial (por ejemplo: "África"). El localizador espacial con el área que ocupa el mapa también puede sernos especialmente útil para este producto.

PDF con Plano del Real Sitio de Aranjuez. 1810

En el Centro de Descargas simplemente tendremos la imagen o documento, sin ninguna información accesoria más, por lo que para todos los detalles de cada documento recomendamos seguir accediendo a su referencia en la Cartoteca (ejemplo anterior).

JPG con Cartas Naúticas de Cádiz. 1584

Desde luego para los amantes de la cartografía antigua se configura un acceso sencillo, potente y rápido a una buena colección de documentos que harán la delicia de cualquier aficionado. Somos de los que compartimos esa opinión de que los mapas antiguos tienen un alma diferente, algo que atrae por enigmático, remoto y trazos particularmente bellos. Agradecemos al IGN esta nueva iniciativa en la que seguramente "perderemos" unas cuantas horas en el futuro.

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.

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.

Modelo Digital de Pendientes MDP05

Disponible en el Centro de Descargas CNIG el Modelo Digital de Pendientes con paso de malla de 5 m, generador a partir del MDT-LIDAR 1ª cobertura. Ficheros cortados por hojas del Mapa Topográfico Nacional a escala 1: 50.000. En las zonas de solape entre husos se ha calculado cada hoja en los dos husos. Cada celda de la malla almacena el valor de la pendiente en grados y fracción de grados sexagesimales. Realice la descarga a través de la agrupación Modelos Digitales de Elevaciones.

Según la información auxiliar del producto se han seguido los siguientes pasos para obtener el modelo digital de pendientes: 

1. Generación de un catálogo con las hojas del modelo digital de elevaciones MDE05, con 5 metros de paso de malla, cortado por hojas 50 mil y obtenido a partir de la clasificación de los datos lídar 1ª cobertura. 
2. Cálculo de la pendiente para cada celda. Para ello se calcula la pendiente entre cada celda y cada una de las celdas colindantes, tomándose de todas ellas la mayor en valor absoluto.
3. El mosaico formado por todas las hojas del mapa de pendientes se recorta por hojas 50 mil para cada huso, manteniendo un solape entre husos de 1 hoja.

Representación del MDP05 con un degradado de color

Formato de archivo ASCII matriz ESRI (asc). Sistema geodésico de referencia ETRS89 (en Canarias REGCAN95, compatible con ETRS89) y proyección UTM en el huso correspondiente a cada hoja del MTN50. En Canarias el huso UTM es el 28. 

El MDP es un producto derivado del MDT el cual fue obtenido a partir de la clasificación de la nube de puntos LIDAR de la 1ª cobertura del proyecto PNOA-LIDAR. En el proceso de cálculo del MDP a cada celda del fichero asc se le asigna el valor de la mayor de las pendientes calculadas entre la celda homóloga del MDT, con iguales coordenadas, y cada una de las celdas colindantes. No se dispone de MDP de Ceuta, Melilla ni Gibraltar ya que no existen datos LIDAR en estas zonas.

Obtención del dato de pendiente en una coordenada concreta

Uso de la calculadora ráster para extraer las zonas con pendientes entre 15 y 30 grados

WMS de los Modelos Digitales de Andalucía 2016.

El artículo de hoy quiere desglosar el servicio web interoperable de la Junta de Andalucía que, a través del Instituto de Estadística y Cartografía de Andalucía, ofrece el Modelo Digital del Terreno (MDT) y los Modelos Derivados generados a partir de datos de dos vuelos fotogramétricos realizados en el año 2016. Uno de ellos de la mitad norte de Andalucía con resolución de 0.5 metros/pixel y otro para la zona Sur de Andalucía con resolución 0.25 metros/pixel.

Decir antes de nada que este MDT 2016 de Andalucía solamente parece estar disponible a través del servicio web, ya que en el portal cartográfico andaluz Line@ solamente se ofrece para descarga en formato TIFF el MDT 2008-2009 (10 metros/pixel). Imaginamos que está en fase de producción y que su resolución estimada dada la de sus vuelos origen será de 5 metros/pixel.

Se trata pues de un servicio WMS cuya URL es https://www.ideandalucia.es/wms/mdt_2016

Cargamos el WMS en QGIS y según observamos en la imagen superior dicho servicio consta de cuatro capas diferenciadas. La primera de ellas, denominada modelo_digital_terreno_2016_color, muestra el Modelo Digital del Terreno del año 2016 representado mediante un gradiente de tintas hipsométricas. Esta clásica representación de colores por alturas ofrece una paleta muy atractiva.

La segunda de las capas se denomina modelo_digital_terreno_2016_escala_grises y muestra el Modelo Digital del Terreno del año 2016 representado mediante un gradiente de escala de grises.

La tercera capa es la llamada modelo_digital_pendientes_2016 y muestra el Modelo Digital de Pendientes obtenido a partir del MDT.

Y la última de las capas, denominada sombreado_orografico_2016, nos muestra el sombreado orográfico de Andalucía obtenido a partir del modelo digital del terreno del año 2016.

Es ésta la única de las capas totalmente opaca, ya que las tres primeras vienen predefinidas con cierto grado de transparencia, ideal para ser superpuestas sobre la del sombreado y conseguir un efecto 3D mucho más realista en cada caso. En ninguna de las capas hemos detectado función GetFeatureInfo por la que podamos obtener la combinación RGB, el dato de altitud o alguna otra información.

En el visor geográfico de la IDE Andalucía está presente, como capa de fondo, una combinación de las capas 1 y 4 de este servicio (MDT color+sombreado), sin duda la combinación más clásica y aparente para representaciones cartográficas, por lo que es también la que nosotros hemos escogido para agregar a nuestro lote de mapas para SASPlanet. Ya es cuestión de los gustos de cada usuario utilizarla como capa, como mapa o superponer ésta o aquella información para conseguir productos más atractivos o útiles.