Nuevos servicios WCS de descarga del MDP y MDS de España.

El Centro Nacional de Información Geográfica (CNIG) ha puesto en marcha dos nuevos servicios de descarga de coberturas (WCS): el Modelo Digital de Pendientes (MDP) y el Modelo Digital de Superficies (MDS). Éste último contiene tres coberturas: el MDS propiamente dicho y dos Modelos Digitales de Superficies Normalizados de las clases Vegetación y Edificación de todos los cuales ya hemos hablado en el blog.

La resoluciones de los modelos son de 5 metros para superficie y pendientes, y 2,5 metros para las coberturas de las clases vegetación y edificación. Gracias a estos nuevos servicios se puede contar con información actualizada del valor de la pendiente en grados sexagesimales y la altura en metros de cada clase. Ambos modelos han sido realizados a partir de las nubes de puntos LIDAR de la primera Cobertura del proyecto PNOA-LiDAR.

Las URL de conexión para dichos servicios WCS son:

• Modelo Digital de Pendientes (MDP): https://wcs-pendientes.idee.es/pendientes?
• Modelo Digital de Superficies (MDS): https://wcs-mds.idee.es/mds?

Estructura de las capas de cada servicio en el asistente de QGIS.

Los servicios WCS nos permiten tanto conseguir una imagen como extraer datos de ella. Por tanto, cargada la capa de cobertura tendremos la información numérica correspondiente de cada punto.

Cargado el MDS05 en QGIS (y con simbología de Mapa de sombras para una mejor visualización) solamente debemos clickar
en un punto con la herramienta Info para obtener la altura en metros.

Y lo mismo con la cobertura de Vegetación.

Usando estos servicios web cualquier operación que requiera extraer el dato que proporciona el WCS es muy sencilla de realizar, como por ejemplo la asignación de alturas a una capa de puntos. Y ello sin necesidad de cargar pesados archivos físicos o de importar una buena colección de ellos para cubrir la misma zona. Es más, si necesitamos extraer a archivo una zona podemos hacerlo sin problema.

Y lo mismo para cualquier tipo de análisis ráster.

La descarga de la imagen desde el servicio se puede realizar en varias proyecciones y formatos, por lo que deberemos elegir el formato más apropiado a nuestra tarea, ya que parece que hay formatos que cargan como simple imagen RGB -la herramienta info nos dará la combinación RGB del pixel pinchado pero no el dato de altitud o pendiente-.

Unos servicios que se antojan imprescindibles por rapidez de uso, accesibilidad y la multitud de utilidades que todo lo relacionado con los archivos de relieve generan en el mundo cartográfico. En un ordenador potente y buena conexión de internet su manejo es toda una gozada, aunque no conviene abusar con la carga de zonas demasiado amplias. Saludos.

Esri 10-Meter Land Cover

La multinacional Esri publica un mapa de coberturas de la superficie terrestre con una resolución de 10 metros y fecha de referencia 2020. El conjunto de datos se desarrolló a partir de imágenes satelitales Sentinel-2 de la Agencia Espacial Europea (ESA) mediante la clasificación de imágenes en diez clases de cobertura terrestre diferentes utilizando un modelo de aprendizaje automático. En esta página de inicio tenemos toda la información de presentación del producto.

Desde allí se nos ofrece acceso a una estupenda Storie del proceso y datos más detallados, así como al mapa del producto, desde el que podremos acceder a una App de descarga de datos.

Podremos muy fácilmente descargar el global del planeta o pinchar en la cuadrícula de nuestro interés y descargar el fichero de datos en formato TIF de la zona.

La precisión de estos datos reclasificados automáticamente se estima por la propia empresa en un 86%, por lo que dependiendo de la zona es posible esperar algunos errores...

Ejemplo de error de clasificación como "áreas construidas" zonas que son simplemente piedra caliza.

Tanto Europa como España en concreto disponen de grandes productos de coberturas de suelo (CORINE Land Cover ó SIOSE son dos buenos ejemplos) que se nos antojan más efectivos y profesionales, sobre todo para estudios de detalle, y no han sido pocas las críticas a este proceso basado únicamente en la IA... Pero en todo caso este mapa de coberturas ofrece las innegables ventajas de ser un producto global disponible en cualquier lugar del planeta y de una potencial actualización más frecuente.

Aspecto de varias cuadrículas del producto ESRI cargadas en QGIS con simbología.

Tampoco conocemos ninguna herramienta abierta para visualizar los datos a través de algún servicio web WMS/WMTS o similar, por lo que suponemos que esa posibilidad queda limitada a los usuarios de los programas de la empresa. A pesar de ello ya sabemos de magníficas iniciativas para poder cargar automáticamente los datos del producto en QGIS a través de las acciones. Saludos.

Modelo Digital de Superficies normalizado de Vegetación

Si hace pocos meses hablamos del Modelo Digital de Superficies normalizado de Edificaciones, ahora le llega el turno al Modelo Digital de Superficies normalizado de Vegetación (MDSnV2,5) correspondiente a los vuelos de la primera cobertura del proyecto PNOA-LiDAR, y con un paso de malla de 2,5 metros. A diferencia de los tradicionales MDT y MDS en los que se almacenan altitudes ortométricas, en este producto se rasterizan las alturas relativas al suelo (considerando el suelo aquellos puntos de la nube clasificados como clase 2) de las clases que contienen información relativa a la vegetación (clases 3, 4 y 5) de los archivos LAS correspondientes a la primera cobertura del proyecto PNOA-LiDAR. El producto se ofrece en formato ASCII matriz ESRI (.asc).

Publicados en el SGR ETRS89 en la Península y en las Islas Baleares, y REGCAN95 en las Islas Canarias (ambos sistemas compatibles con WGS84) y en proyección UTM (en el huso correspondiente), la unidad de descarga son las hojas del MTN50 (a excepción de las hojas de Ceuta, Melilla e Isla de Alborán (1110, 1111, 1078B)).

Son necesarios para conocer las alturas de la vegetación y las masas forestales respecto del suelo y para tener un modelo digital de elevaciones normalizadas homogéneo y preciso para todo el territorio español. De utilidad, por ejemplo, para estudios en los que se requieran las altitudes de la vegetación, para analizar la evolución de las masas forestales, estudio de impacto de incendios forestales o el análisis de zonas verdes en entornos urbanos. 

La descarga puede hacerse a través de la agrupación Modelos Digitales de Elevaciones y se acompaña de información auxiliar con conversor entre formatos ASC y XYZ, descripción del formato ASC y archivo Excel con fechas y resoluciones de cada hoja MTN50 del producto. Abrimos en QGIS una hoja para mostrar su contenido, concretamente la hoja 105 del MTN50 cuyo archivo descargado lleva por nombre NDSM-Vegetacion-ETRS89-H30-0105-COB1.asc

El tamaño del ASC de esta hoja 105 son 376 megas.

Como hemos contado el ráster muestra información solamente en aquellas zonas que los puntos LIDAR clasifican como vegetación, quedando el resto en blanco. A partir de aquí podemos interactuar con el archivo como con cualquier otro modelo de relieve, teniendo siempre en cuenta que las alturas que nos ofrece cada pixel son las relativas al suelo, nunca las absolutas.

Obviamente suponemos que estos datos no presentan una exactitud impecable ni exenta de errores en todos los casos. Por ejemplo, en nuestra hoja de prueba para este artículo, estamos casi seguros que no existe ningún punto de vegetación con una altura de 100 metros respecto al suelo, tal como nos "canta" el archivo. Igualmente ciertos puntos del ráster arrojan una altitud negativa.

Comparativa sobre ortofoto de las zonas captadas como vegetación por la nube de puntos.

Pero en su conjunto se trata de un gran producto normalizado, de extensión nacional y de enorme validez y utilidad. Y es que las opciones de generar productos derivados de la información obtenida en un vuelo LIDAR multiplican las posibilidades de generar estudios geográficos en muy diversos ámbitos. Y lo mejor: al alcance de todos. Gracias de nuevo al IGN por hacerlo posible.

Servicio ATOM del FEGA de descarga del SIGPAC Nacional, con datos de cultivos.

Es posible que los términos sobre los que trata este artículo necesiten una breve explicación introductoria, así que, tirando de las propias web oficiales, vamos con ello.

¿Qué es el SIGPAC?

El Sistema de Información Geográfica de Parcelas Agrícolas (SIGPAC) permite identificar geográficamente las parcelas declaradas por los agricultores y ganaderos, en cualquier régimen de ayudas relacionado con la superficie cultivada o aprovechada por el ganado.

Concebido inicialmente con el propósito de facilitar a los agricultores la presentación de solicitudes, con soporte gráfico, así como para facilitar los controles administrativos y sobre el terreno, el SIGPAC se ha convertido en una herramienta de enorme utilidad en campos diferentes del agrario (geología, infraestructuras, urbanismo...), lo que obedece a su concepción y desarrollo, en el que se hace uso continuo y permanente de las tecnologías más avanzadas en información geográfica automatizada.

¿Qué es el FEGA?

El Fondo Español de Garantía Agraria es el organismo autónomo adscrito al Ministerio de Agricultura, Pesca y Alimentación que tiene como misión principal hacer que los fondos europeos asignados a España se apliquen estrictamente a lograr dichos objetivos.

Con fecha 20 de diciembre de 2014 se ha publicado el Real Decreto 1077/2014, de 19 de diciembre, por el que se regula el sistema de información geográfica de parcelas agrícolas, que atribuye al FEGA el carácter de autoridad competente, a nivel nacional, al estar facultado para la coordinación de los controles a realizar en el ámbito del sistema de gestión y control de las ayudas, así como para llevar a cabo el seguimiento de las actuaciones que, en esta materia, realicen las Comunidades Autónomas, con objeto de garantizar la aplicación armonizada en el territorio nacional de la reglamentación comunitaria

¿Qué es un servicio ATOM?

Es un formato de redifusión web de contenido y metadatos basado en XML. La redifusión web permite publicar información desde una fuente web que luego se podrá consultar, para determinar y obtener la información publicada, así como suscribir, para conocer sus actualizaciones posteriores. Existen otros formatos de redifusión web como RSS. Sin embargo, ATOM es el indicado por la directiva INSPIRE de la UE para los servicios de descargas.

Siendo ATOM un formato de uso en la web, dicho cliente solía estar integrado en los navegadores web, de los que MS Internet Explorer 11 sigue siendo un caso aún vigente. Sin embargo, los navegadores actuales, como Chrome, Edge, Firefox o Safari, ya no lo incorporan, por lo que se requiere la instalación de una extensión que añada esa funcionalidad al navegador.

Por lo tanto, dicho lo anterior y resumiendo, hablamos en este artículo de conseguir los archivos vectoriales del SIGPAC a nivel nacional desde el servicio de descargas ATOM proporcionado por el FEGA.

La puerta de entrada al servicio ATOM es https://www.fega.es/orig/

Web principal en la que tenemos un enlace "¿Cómo funciona este servicio?" muy recomendable que nos explica cómo instalar extensiones para ATOM en otros navegadores, una guía completa con las especificaciones del formato o el texto del Acuerdo de distribución del FEGA.

Nosotros hemos usado Internet Explorer por los motivos que se han citado más arriba, y una vez aceptadas las condiciones, este es el aspecto de la página con los enlaces ATOM.

Junto a los enlaces de descarga para todas las provincias se proporcionan las imprescindibles Tablas de códigos usadas por los datos SIGPAC.

Las tablas con los códigos contienen un lote de archivos CSV de texto que nos serán de gran utilidad para establecer las diferentes relaciones especiales entre las geometrías y los datos mediante uniones.

Ejemplo del contenido de la tabla con los códigos de uso de suelo SIGPAC.

Tras ubicarnos en la provincia de nuestro interés, el enlace nos adentra en una nueva página con los datos desglosados por sus municipios. Nos situamos por ejemplo sobre el de Boca de Huérgano (24021), perteneciente a la provincia de León (24).

La descarga principal consiste en un archivo único en formato GeoPackage cuya nomenclatura -24021_20210104.zip- es el código INE del municipio seguido de la fecha de actualización de los archivos. Abramos el archivo en QGIS para ver su contenido.

El archivo GPKG nos da la opción de cargar alguno o todos los elementos de SIGPAC, tablas de códigos ya incluidas, que deberíamos cargar de forma individual si utilizamos los archivos en formato Shapefile.

A partir de aquí es sencillo unir las geometrías con sus respectivas tablas de códigos para enriquecer las tablas de atributos con descripciones o nuevos campos que nos permitan por ejemplo una rápida categorización de las parcelas por su uso...

... o consultas espaciales más complejas a través de la calculadora de campos:

Selección (en amarillo) mediante expresión de las parcelas compuestas por pasto con arbolado con una pendiente media mayor de 40% y una superficie mayor a 500 metros cuadrados.

En definitiva, las posibilidades de combinación de la información en estos grandes productos geográficos son innumerables, y el partido que se puede sacar a un buen conocimiento de los datos es inmenso. En este sentido una concienzuda lectura del documento PDF describiendo la estructura de las geometrías y de las tablas de atributos se antoja imprescindible. Y los conocimientos sobre expresiones para extraer datos cruzados en un SIG muy recomendables. Saludos.

Ortofotos históricas PNOA 10

Si en nuestro anterior artículo sobre las ortofotos expeditas acabábamos anunciando una nueva capa en el servicio WMS de ortofotos históricas con el nombre PNOA10, pocos días después se nos anuncia la disponibilidad de este producto en su formato físico para descarga.

Se trata de mosaicos de ortofotos de diferentes años de vuelos realizados por diversos organismos de Administraciones Públicas a 0,10 m de resolución de píxel, en el formato COG (Cloud Optimized GeoTIFF) del que ya hablamos a cuenta de los mosaicos Landsat. Las imágenes se sirven en cortes de hoja 1:2.000 (subdivisión de cada MTN50 en 20 x 20 hojas) o en cortes de superficies de 1x1 km de extensión. Por el momento se dispone de imágenes de Galicia, Castilla-La Mancha, Islas Baleares y Madrid, pero se irán incorporando nuevas zonas.

La información auxiliar del producto nos proporciona los vectores con las zonas, fechas y resoluciones de las imágenes disponibles a fecha de este artículo.

Este producto tiene como objetivo la recopilación de ortofotografías aéreas digitales producidas por diversas Administraciones Públicas, realizadas en zonas de interés (básicamente núcleos de población y zonas de costa), para ponerlas a disposición de los usuarios. Se realiza un vuelo fotogramétrico único y un tratamiento riguroso de los datos cumpliendo con unas especificaciones técnicas consensuadas entre todas las Administraciones Públicas productoras. Este planteamiento es acorde con el espíritu de la Directiva Inspire para el establecimiento de una Infraestructura de datos geográficos en Europa, que persigue que el nivel de detalle máximo de la información geográfica se capture una sola vez y que se comparta abiertamente entre los distintos agentes sociales que tengan necesidad de la misma. Los recubrimientos de las ortoimágenes forman parte de los "datos básicos de referencia" recogidos en dicha Directiva.

La colección PNOA 10 histórico se estrena con un total de 3840 imágenes.
El buscador nos permite discriminarlas por año de vuelo.

Así por ejemplo, para el entorno de la ortofoto de Toledo del vuelo de octubre 2007, se nos sitúa en el marco de su hoja MTN50 (gris) la zona cubierta por las ortofotos PNOA 10 (naranja), compuesta en este caso por los 39 archivos COG correspondientes.

Recordamos que los archivos COG son a efectos de visualización un archivo formato TIF, por lo que programas habituales para nosotros como Global Mapper o QGIS los abren sin ningún problema. Sus tamaños son muy heterogéneos dependiendo del corte de imagen aplicado y podemos encontrar imágenes en rangos de 10 hasta 200 megas. La calidad de imagen es también variopinta.

Calidad y resolución de la ortofoto PNOA10 2018 visionada a pixel real (izquierda) frente a la ortofoto máxima actualidad 25 cms 2017 (derecha) en la ciudad de Lugo.

Comentar por último como dijimos al comienzo la posibilidad de acceder a estas imágenes PNOA10 a través de sus respectivas capas del servicio WMS de ortofotos históricas y PNOA anual (que no parece ir muy fino últimamente, dicho sea de paso...).

En definitiva un nuevo producto disponible en esta escalada de novedades en la que está inmerso nuestro IGN este comienzo de año 2021. Enhorabuena una vez más por el trabajo y el acercamiento de todos estos interesantes productos al público general.