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.

De QGIS al mundo Google.

Hoy traemos una rápida entrada con dos complementos muy sencillos que conectan nuestro lienzo de QGIS con el  mundo cartográfico Google, lo que siempre puede resultar útil e interesante.

El primero de ellos es Send2GE (Send to Google Earth), simple herramienta que graba la posición del mapa en la que pinchemos y la envía a Google Earth, creando una marca de posición. Si tenemos Google Earth cerrado el complemento lo arranca automáticamente. Y si marcamos varias posiciones éstas se van acumulando en los lugares temporales de Earth.

El segundo complemento es Street View, que abre la legendaria vista de pie de calle de Google en nuestro navegador, justo en el punto donde hayamos pinchado en QGIS. Si al pinchar estiramos con el ratón, una línea roja nos permite escoger la dirección de la vista.

Ambos complementos generan al instalarse tanto una entrada en el menú como un icono en la barra de botones.

Por cierto que probando en lugares donde no había calle física (grupo de edificios en este ejemplo) me he llevado alguna sorpresa y Street View me ha abierto vistas como esta....

Pero obviamente si pinchamos en lugares alejados de toda cobertura de Street View el navegador nos devolverá una imagen en negro. Esto se echa un poco en falta, algún indicador que nos asegure o no si Street View cubre la zona en la que vamos a pinchar.

En resumen dos complementos muy sencillos y sin ninguna pretensión ni complicación que unen nuestro escritorio de QGIS con el mundillo Google a través de Earth o Street View, y que muchas veces pueden abrirnos las puertas a trabajar con diferentes puntos de vista desde una misma posición. Saludos.

Mapa LIDAR

El Centro de Descargas del CNIG contiene en su categoría Modelos Digitales de Elevaciones los llamados mapas LIDAR. Se trata de unos modelos digitales de superficies sombreados a los que se les han superpuesto dos capas procedentes de la rasterización de las clases edificación y vegetación de las nubes de puntos en formato .las y una tercera de hidrografía. 

Tenemos los derivados de la primera cobertura LIDAR que cubren todo el territorio nacional y ya llevan un tiempo disponibles para descarga. Se ofrecen en el Sistema Geográfico de Referencia WGS84, proyección Web Mercator (EPSG:3857). Y el motivo de este artículo es que el IGN ha comenzado la publicación del mapa LIDAR derivado de la segunda cobertura, para el que en esta ocasión ha optado por un Sistema Geográfico de Referencia ETRS89 en la Península, Islas Baleares, Ceuta y Melilla, y REGCAN95 en las Islas Canarias (ambos sistemas compatibles con WGS84) y proyección UTM en el huso correspondiente.
Datos comunes para ambas coberturas son el formato de archivo de descarga ECW y la unidad de dicha descarga (hojas del MTN50). También el proceso de elaboración que el mismo IGN nos especifica:

Este producto se genera por composición de cuatro capas, debidamente superpuestas. A partir de las nubes de puntos LiDAR correspondientes a la segunda cobertura del proyecto PNOA-LiDAR, debidamente clasificadas, se genera un modelo digital de superficies al que se le aplica un sombreado. Sobre este modelo se superponen tres capas: dos modelos digitales de superficies normalizados, uno de la clase vegetación y otro de la clase edificación, y la rasterización de la clase "agua" de las nubes de puntos LiDAR. Se aplica una paleta de color a cada capa, estableciendo un degradado en la información correspondiente a la vegetación y a la edificación en función de su altura sobre el terreno.

A fecha de este artículo las hojas disponibles de la segunda cobertura del mapa LIDAR son 669, con esta extensión geográfica (similar en cierta lógica a la del MDT02 en misma fecha):

La información auxiliar proporciona una cuadrícula en formato shape con las hojas disponibles, su nomenclatura y las fechas de vuelo.

En lo que respecta al mapa propiamente dicho se trata de una composición muy atractiva. La combinación del relieve sombreado con el coloreado de los elementos más importantes del paisaje (edificaciones, vegetación y aguas) resultan en un lienzo realista que no deja indiferente y con el que es más sencillo realizar escenas 3D visualmente buenas.

El mapa LIDAR cuenta así mismo con sendos servicios web disponibles que podemos usar en nuestros softwares favoritos. Además se puede consultar la altura de los edificios y vegetación a partir de cierta escala (1:20.000):

WMS........... http://wms-mapa-lidar.idee.es/lidar?request=GetCapabilities&service=WMS

WMTS......... http://wmts-mapa-lidar.idee.es/lidar?request=GetCapabilities&service=WMTS

Gracias a estos servicios web y SASPlanet pudimos elaborar por ejemplo los fondos JNX para Garmin que tenemos subidos en nuestra sección Garmin JNX del blog y que son una buena alternativa para complementar con nuestros mapas vectoriales en BaseCamp.

Ejemplo con TopoValencia para Garmin y su fondo LIDAR en BaseCamp.

Por último decir sobre los archivos físicos a descarga, los ECW's, que están elaborados a una resolución de 2,5 metros y con una compresión 1:6. La nomenclatura de estos archivos es un libro abierto, por ejemplo: MapaLidar-COB2-H30-0509.ecw Es decir nombre de producto, cobertura 2, Huso 30, hoja del MTN50 y formato. Una imagen de su aspecto a zoom 100%:

En definitiva un mapa LIDAR fantástico como fondo o mapa base para aplicaciones web y para el trabajo en entornos SIG. Un mapa además que cuenta con todas las versiones deseables en un producto TOP: servicios web y archivos físicos, para que cada usuario elija la forma de trabajo que mejor le convenga. Enhorabuena y gracias una vez más a nuestro IGN.