CORINE Land Cover 2018

Como anunciamos hace unos días en nuestra página de Facebook, ya se encuentra disponible la nueva versión de la base europea de ocupación del suelo Corine Land Cover 2018 (CLC2018). Podemos obtenerla a través del Centro de Descargas del CNIG, y se trata de la actualización pertinente de Corine 2012 de la que ya hablamos en el blog. Analizamos su contenido.

Listado de productos por año Corine disponibles en el Centro de Descargas, incluida la versión 2018 objeto de este artículo. Incluye las versiones de 1990, 2000, 2006, 2012 y 2018, así como los mapas de cambios de ocupación del suelo entre los años 1990-2000, 2000-2006, 2006-2012 y 2012-2018.

Lo primero que llama la atención de esta nueva versión de Corine es los formatos en los que se ofrece. Al igual que sucedió con el reciente SIOSE 2014, el IGN ha abandonado el formato shape SHP en favor de los formatos vectoriales .GDB (ArcGIS) y GeoPackage. Es este último formato abierto .GPKG el que nos interesa en Cartografía Digital, no sin antes reiterar nuestra extrañeza por el abandono de un formato actualmente universalizado como .shp en favor de uno comercial.

El IGN ofrece también acceso a la descarga de información auxiliar del producto (con exhaustiva documentación técnica sobre su creación) y enlace a los metadatos de cada versión, en los que encontramos información relevante como esta definición general: Mapa de ocupación del suelo en España correspondiente al proyecto europeo CORINE Land Cover (CLC) con una nomenclatura de 44 clases, versión del año de referencia 2018. En España el Sistema geodésico de referencia utilizado es ETRS89 (en Canarias, WGS84/REGCAN95) y proyección UTM huso 30 (en Canarias, huso 28). El proyecto CORINE Land Cover (CLC) es responsabilidad desde 1995 de la Agencia Europea del Medio Ambiente con el objetivo fundamental de obtener una base de datos europea de ocupación del suelo a escala 1:100.000, útil para el análisis territorial y la gestión de políticas europeas. En la actualidad este proyecto está incluido en el Programa «Copernicus», cuyo principal objetivo es establecer un Sistema de Observación de la Tierra bajo el mandato de la Comisión Europea.

O este otro interesante párrafo: El CORINE Land Cover se obtuvo mediante fotointerpretación sobre imágenes de satélite Sentinel 2, Landsat en falso color (combinación de bandas 453) y Spot como referencia básica del proyecto, utilizándose como apoyo otras bases de datos geográficas, como mapas topográficos o temáticos relativos a la ocupación del suelo y ortofotografías aéreas, así como otra información exógena de apoyo a la fotointerpretación. La producción se produjo de manera coordinada entre diferentes organismos públicos. Sobre las imágenes se delimitaron áreas o zonas susceptibles de representar una unidad de ocupación del suelo según CORINE Land Cover, asociando a cada zona una clase de nomenclatura CLC y extendiendo esta delimitación a un conjunto de zonas de la imagen que presentaran características similares, consiguiendo una base de datos continua para todo el territorio.

La descarga del paquete comprimido 2018 - CLC2018-GPKG.ZIP correspondiente a todo el territorio nacional tiene un tamaño de 1,23 Gigas, y en su interior encontramos el siguiente contenido:

- CLC2018_ES.gpkg: es la base vectorial Corine propiamente dicha en formato GeoPackage.

- LEER_DatosCLC2018_gpkg.pdf: archivo PDF explicativo de la descarga y que contiene los imprescindibles códigos de la nomenclatura Corine.

- LEYENDA_CLC.lyr y LEYENDA_CLC.sld: archivos con la simbología de capas CLC en esos dos formatos (.lyr para ArcGIS y .sld para software abierto).

Una vez descomprimido nuestro archivo GeoPackage vemos que tiene el nada desdeñable tamaño de 2,83 Gigas, por lo que un ordenador de cierta potencia se agradecerá en su manejo. Con toda seguridad las diferentes regiones irán ofreciendo próximamente recortes autonómicos de un tamaño más comedido del producto. Cargamos Corine en QGIS para ver lo que nos ofrece:

De nuevo observamos las bondades de este formato .gpkg abierto en QGIS: la simbología de capas se muestra por defecto, lo que hace que la primera impresión del producto ya sea inmejorable. En cuanto al contenido, nuestro archivo CLC2018_ES.gpkg nos ofrece cuatro capas perfectamente diferenciadas (QGIS nos da la opción previa de cargarlas todas o de forma individual).

- CHA18_ES: capa de cambios 2012-2018 para Península e Islas Baleares.
- CHA18_ES_Canarias: idem para las Islas Canarias.
- CLC18_ES: Corine Land Cover 2018 para Península e Islas Baleares.
- CLC18_ES_Canarias: idem para las Islas Canarias.



¿Qué campos nos ofrece la tabla de atributos de cada capa?

Las capas correspondientes a la versión 2018: CODE_18: código CLC. AREA_HA: superficie del polígono en hectáreas. Las capas correspondientes a los cambios 2012-2018: CODE_12: código CLC en 2012. CODE_18: código CLC en 2018. CHANGE: código 2012-código 2018. 

La documentación adjunta a CLC2018 afirma que las capas de cambios 2012-2018 llevan también el atributo AREA_HA, pero parece que no. A ellos hay que sumarle el campo ID que otorga un identificador único en la tabla a cada polígono.
Obviamente para trabajar con los datos CLC se hace imprescindible conocer los códigos que identifican el uso de cada parcela, y que como hemos comentado se encuentran en el archivo PDF de la descarga.

Tabla de códigos CLC por categoría correspondientes al atributo CODE del archivo vectorial.

De esta forma sabremos rápidamente que un CODE 312 (3.1.2.) corresponde a una parcela con bosque de coníferas, o que el CODE 222 (2.2.2.) identifica una zona con frutales. Igualmente si el atributo CHANGE es 333-211, sabremos de inmediato que la parcela en cuestión era un espacio con vegetación escasa en 2012 y ha pasado a ser en 2018 tierra de labor en secano. La herramienta Identificar objetos espaciales de QGIS es para ello de gran comodidad.

Incluso la comodidad sería aún mayor si agregáramos un campo a las tablas con el texto de cada CODE, lo que nos evitaría tener que consultar la documentación para conocer el uso descriptivo del polígono.

De la misma forma que trabajábamos con los datos SIOSE, aunque la complejidad de códigos es aquí mucho menor al haber un único uso de suelo por parcela, podemos efectuar consultas espaciales sobre los datos CLC. Un ejemplo sencillo sería buscar las parcelas catalogadas como viñedos y con una extensión mayor de 3 hectáreas y menor de 50:

Resultados de la selección por expresión reflejados en la tabla de atributos: 17 parcelas.

Respecto a la capa de cambios producidos entre CLC2012 y CLC2018 (casi 20000 en la capa peninsular) no hay mucho que añadir. Esta capa no contiene simbología integrada y nos discrimina aquellos polígonos que han sufrido algún cambio en los usos del suelo (o simplemente en su extensión manteniendo el uso) durante dicho periodo.

Capa de cambios (negro) mostrando aquellas parcelas CLC2012 que en 2018 (derecha) han cambiado de uso o extensión.

No queremos terminar sin hablar del archivo con la simbología de capas presente en la descarga, el mencionado LEYENDA_CLC.sld de formato abierto. Os preguntaréis para qué sirve si el archivo GeoPackage ya se abre en QGIS con la simbología correcta cargada... Aparte de que el archivo puede servirnos para otros softwares que no actúen como lo hace QGIS, imaginaos que yo he hecho un recorte de CLC2018 de la zona de Valencia y lo he guardado en formato SHP. Cuando cargue dicho archivo en QGIS lo hará con una simbología única:

Peeeero gracias a nuestro archivo SLD de simbología nos será muy fácil cargarle el estilo CLC2018. Nos vamos a las Propiedades de la capa, pestaña Simbología, botón Estilo en la parte inferior. 

Seleccionamos Cargar estilo... y buscamos la ruta a nuestro archivo LEYENDA_CLC.sld

Pulsamos en Cargar estilo y.....

La simbología basada en reglas de CLC2018 se carga en nuestra capa.

Simplemente aplicamos y aceptamos y nuestro recorte SHP queda configurado con la simbología CORINE del GeoPackage original.

En realidad trabajando con QGIS podríamos haber salvado nosotros mismos el estilo del archivo GeoPackage en un archivo SLD, pero desconozco si otros softwares cargan el GPKG con la simbología integrada, por lo que este archivo de estilo sería muy importante.

El producto CLC2018 ya se encuentra así mismo disponible a través de los servicios web de ocupación del suelo del IGN (WMS, WMTS y WFS) con la misma simbología vista en este artículo.

Y ya que teníamos todo cargado no nos hemos resistido a jugar un poco con el diseñador de impresión para finalizar.... Saludos!!

Nuevos servicios de Ocupación del Suelo.

Nuestro IGN, en su sección de Noticias, informa de la puesta en marcha de nuevos servicios de mapas correspondientes a la información sobre Cubierta Terrestre y Usos del Suelo, y de que estos servicios sustituyen a los que hasta la fecha venía publicando el propio IGN sobre esta temática. Se trata concretamente de servicios de visualización WMS/WMTS y del servicio de descarga WFS.

Os dejamos sus direcciones URL enlazadas a su correspondiente archivo de capacidades:

                                             WMS:    http://servicios.idee.es/wms-inspire/ocupacion-suelo?

                                             WMTS:  http://servicios.idee.es/wmts/ocupacion-suelo?

                                             WFS:     http://servicios.idee.es/wfs-inspire/ocupacion-suelo?

En dichos archivos de capacidades encontramos mucha información sobre estos nuevos servicios conformes con las especificaciones de la Directiva Inspire. Así el resumen del servicio WMS nos dice: "Información de Ocupación de Suelo de España (SIOSE) 2011 y CORINE Land Cover 2012. La denominación de las capas es conforme con las especificaciones de la Directiva Inspire 2007/2/EC (nombre, título) así como estilo Inspire por defecto. . Servicio de visualización WMS 1.3.0 conforme al perfil Inspire de ISO 19128:2005". Consta de dos capas cuya estructura y árbol vemos en la siguiente imagen de QGIS:

La primera capa lleva por título "Superficies de Cubierta terrestre" y por nombre de layer LC.LandCoverSurfaces. El propio servicio nos resume su contenido: "Representación de datos espaciales correspondientes al tema Cubierta terrestre del Anexo II de la Directiva Inspire y Anexo I de LISIGE. Datos procedentes CORINE Land Cover 2012 y SIOSE 2011. Para escalas menores de 1:50.000 se distinguen las clases correspondientes al nivel 3 de CORINE Land Cover; para mayores, la clasificación de cubierta terrestre correspondiente a la nomenclatura establecida por el CODIIGE para SIOSE." Y también encontramos enlace a la imagen de la Leyenda.

Distintos estilos de representación con las diferentes clases atendiendo a las diferentes escalas.

La segunda capa del servicio WMS lleva por título "Usos de suelo existentes" y por nombre técnico de layer LU.ExistingLandUse. Su correspondiente resumen (Representación de datos espaciales correspondientes al tema "Usos del Suelo" del Anexo III de la Directiva Inspire y Anexo II de LISIGE. Datos procedentes de CORINE 2012 y Sistema de Información de Ocupación del Suelo en España (SIOSE) 2011 . Se utiliza la nomenclatura de clases de Usos del Suelo HILUCS, distinguiendo las siguientes clases: 1_1 Uso agrícola, 1_2 Uso forestal, 1_3 Minas y canteras, 1_4 Pesca y acuicultura, 2 Sector secundario, 3 Sector terciario, 3_1 Servicios comerciales, 3_3 Servicios comunitarios, 3_4 Servicios recreativos y culturales, 4 Redes de transporte y logística, 4_1 Redes de transporte, 4_3 Utilidades, 5 Uso residencial, 6_1 Áreas de transición, 6_2 Áreas abandonadas, 6_3 Áreas naturales sin uso económico, 6_3_1 Zonas terrestres sin uso económico, 6_3_2 Zonas de agua sin uso económico, 6_6 Uso desconocido.) y su enlace a la Leyenda.

Nos ofrece también detalles alternativos y más específicos sobre los códigos de colores HTML empleados en cada superficie: El color del relleno depende del valor de la clase de Uso del Suelo de HILUCS (1_1 Uso agrícola, #E6E66E; 1_2 Uso forestal, #6EE66E; 1_3 Minas y canteras, #A64D00; 1_4 Pesca y acuicultura, #8CE6FF; 2 Sector secundario, #646464; 3 Sector terciario, #969696; 3_1 Servicios comerciales, #9696C8; 3_3 Servicios comunitarios, #96B496; 3_4 Servicios recreativos y culturales, #B49696; 4 Redes de transporte y logística, #B478F0; 4_1 Redes de transporte, #E678F0; 4_3 Utilidades, #B450F0; 5 Uso residencial, #F07864; 6_1 áreas de transición, #F0D2D2; 6_2 áreas abandonadas, #D2FAD2; 6_3 áreas naturales sin uso económico, #D2D2F0; 6_3_1 Zonas terrestres sin uso económico, #D2F0FA; 6_3_2 Zonas de agua sin uso económico #F0F0D2; 6_3 Uso desconocido, #E6E6E6). Se visualiza a escalas mayores de 1:350.000.
El servicio cuenta también con capacidades GetFeatureInfo, con lo que pinchando con la herramienta información dispondremos de una ventana con la información del polígono sobre el que hayamos pedido los datos.

Segunda capa del servicio WMS con previsualización de ventana de datos GetFeatureInfo.

Respecto al servicio WMTS poco más que añadir. Simplemente se trata de la misma información servida a la manera teselada propia de este servicio, lo que agiliza la respuesta a la hora de pedir las imágenes. Su resumen: Información de Ocupación de Suelo de España (SIOSE) 2011 y CORINE Land Cover 2012. La denominación de las capas es conforme con las especificaciones de la Directiva Inspire 2007/2/EC. Servicio de visualización Teselado conforme al perfil INSPIRE de Web Map Tile Service (WMTS) 1.0.0. Y la estructura del servicio con sus posibles conjuntos de teselas y proyecciones:

Por último disponemos también del servicio de descarga WFS. Sus capacidades nos dicen que es compatible con las versiones 1.0, 1.1 y 2.0, y cargado en el administrador de QGIS nos ofrece esta estructura:

El cliente por defecto de QGIS para servicios WFS nos arroja error en todos los conjuntos y no nos muestra la info, por lo que hemos recurrido al complemento WFS 2.0 que ya nos ha sacado del apuro en otras ocasiones. No hemos conseguido encontrar en el archivo de capacidades la limitación máxima de elementos para cada petición.

Tabla de atributos de info servida desde el WFS.

En definitiva este es el repaso que queríamos darle a estos nuevos servicios mostrando sus características y funcionamiento. La Ocupación del Suelo es uno de los grandes productos cartográficos, fuente y origen de innumerables posibilidades en cuanto a estudios del terreno, desde la faceta medioambiental hasta la agrícola. Un saludo.

CORINE Land Cover 2012.

Ya en el mes de mayo de este año hicimos un primer artículo general sobre el sistema de ocupación del suelo europeo denominado CORINE Land Cover. El proyecto CORINE Land Cover (CLC) es responsabilidad desde 1995 de la Agencia Europea del Medio Ambiente con el objetivo fundamental de obtener una base de datos europea de ocupación del suelo a escala 1:100.000, útil para el análisis territorial y la gestión de políticas europeas. En la actualidad este proyecto está incluido en el Programa «Copernicus», cuyo principal objetivo es establecer un Sistema de Observación de la Tierra bajo el mandato de la Comisión Europea.

Si en aquella ocasión los datos disponibles eran de 2006, ahora el CNIG acaba de publicar a través de su Centro de Descargas la última actualización de este excelente producto, correspondiente al año 2012.

Bajo el enlace Documentación técnica accedemos a un par de archivos PDF (en inglés) con documentación exhaustiva y pormenorizada del proyecto: método de elaboración, guías técnicas, validación del producto, explicación exhaustiva de las categorías, etc, etc... Una lectura de nivel para aquellos interesados en los pormenores más técnicos del asunto.

Y bajo en enlace Descargar accedemos al catálogo de productos CORINE disponibles.

Además de las diferentes versiones del producto CORINE de cada año y demás archivos de metadatos y cambios
 (por ejemplo los cambios de usos del suelo entre la última versión 2006 y la actual 2012 -en amarillo-), encontramos
 la descarga del conjunto de datos CORINE 2012 para España objeto de este artículo (en verde).

Lo primero que nos llama la atención es el aumento de tamaño de la descarga en esta nueva versión 2012 (951 megas frente a los 341 de 2006), por lo que hay que tener claro que estamos hablando de archivos de manejo pesado y gran cantidad de datos. De hecho, el contenido de CORINE 2012 ya descomprimido en carpeta, supera los 2 Gigas.

Dicho contenido es el siguiente:

- Archivo .lyr correspondiente a la simbología de los vectoriales para ArcGIS.

- Vectorial clc12_ca.shp con los datos CORINE 2012 para Canarias.

- Vectorial clc12_es_1.shp con los datos CORINE 2012 correspondientes a Galicia, Asturias, Cantabria, País Vasco, Navarra, La Rioja, Aragón, Cataluña, Castilla y León, Comunitat Valenciana e Illes Balears.

- Vectorial clc12_es_2.shp con los datos CORINE 2012 correspondientes a Extremadura, Madrid, Castilla-La Mancha, Andalucía, Murcia y Ciudades Autónomas de Ceuta y Melilla.

- Archivo PDF resumen del producto y conteniendo las 44 clases CORINE.

Tres consideraciones queremos hacer sobre este contenido:

1.- Constatar, como podéis ver en la imagen, el enorme aumento de tamaño que han sufrido los datos CORINE respecto a la versión 2006 y que han obligado a dividir la Península en dos archivos al alcanzar las limitaciones de tamaño del propio formato shape. Esto puede hacer que el manejo de la información sea engorroso, o directamente imposible, sin un PC de cierta potencia.

2.- Lamentar la posición de privilegio de ArcGIS, única plataforma para la que se ofrece archivo contenedor de simbología de los vectoriales (formato lyr), ignorando a softwares libres de gran calado en el mercado y rompiendo la lógica filosofía open data en beneficio de programas comerciales. Un archivo SLD compatible con ellos, y que además es el estándar OGC para servicios WMS, no hubiera estado de más.

3.- Agradecer la inclusión del PDF resumen de las clases de suelo CORINE, simplificado y traducido al castellano, pues los archivos técnicos en inglés no están al alcance de todos.

Estos 44 códigos de suelo, o clases CORINE, son el contenido fundamental para poder clasificar, categorizar los vectoriales y poder de esta forma trabajar con ellos. Abramos el vectorial clc12_ca.shp correspondiente a Canarias en Global Mapper para ver su aspecto:

CORINE 2012 de Canarias desagrupado por los códigos de las diferentes clases de suelo.

De esta forma, y ayudados por los correspondientes códigos del archivo PDF, podemos saber que los polígonos con código 133 pertenecen a zonas en construcción, o los del código 243 a "Terrenos principalmente agrícolas, pero con importantes espacios de vegetación natural". Otro ejemplo: detalle de la isla de Gran Canaria y todas las clases de suelo categorizadas con colores aleatorios, esta vez en QGIS.

Comentar también que los datos CORINE 2012 ya están incluidos en el servicio WMS del IGN sobre ocupación del suelo en España, cuya URL es http://www.ign.es/wms-inspire/ocupacion-suelo? Dichos datos se visualizan por defecto en la capa llamada LC.LandCoverSurfaces y con la simbología correspondiente a la leyenda del siguiente enlace para cada escala. Os mostramos un ejemplo del servicio abierto en Global Mapper y usando la herramienta Get Feature Info que nos permite obtener datos directos sobre el polígono.

En definitiva, una esperada actualización de uno de los grandes productos cartográficos que, si bien adolece del inevitable desfase con la realidad de estos productos complejos (pienso por ejemplo en las clases zonas en construcción o zonas quemadas cinco años después de la recogida de datos), supone una grandísima noticia de la que seguramente volveremos a hablar muy pronto. Saludos.

CORINE Land Cover.

Si el mes pasado hablábamos del SIOSE, el producto estrella en cuanto a ocupación del suelo del IGN, no podemos dejar de rematar la faena hablando de CORINE Land Cover (CLC), el otro producto destacable dentro de esta temática en nuestro imprescindible Centro de Descargas.

Lo primero que se nos ocurre pensar es ¿qué diferencias tiene con el SIOSE? Podéis echar un vistazo a la web del IGN sobre el asunto, o también a la web de la EEA para una completa información. Nosotros, en nuestro estilo de andar por casa, resumiríamos las principales diferencias así:

• La escala de referencia de CORINE es 1:100000 frente al 1:25000 del SIOSE.
• La unidad mínima cartogafiable por defecto son 25 hectáreas (0,5 a 2 hectáreas en SIOSE).
• La clasificación del terreno y sus niveles de jerarquía son más simplificados que en SIOSE.

Digamos que, básicamente, CORINE Land Cover es un proyecto a nivel europeo que trabaja la información a mayores escalas, por lo que es un producto más simplificado tanto en polígonos como en categorías de datos. Como solemos hacer, una vez ofrecidos los enlaces para interesados, vamos a lo práctico. Un primer vistazo a las posibles descargas del producto nos arroja estos resultados:

Junto a los archivos de metadatos de cada serie (naranja), lo cierto es que aparecen varias descargas de varios productos y actualizaciones (amarillo) y el producto de 2006 (verde) que da la impresión de ser el único unificado e ideal para descargar.

Las descargas en este producto, a diferencia del SIOSE que son por Comunidades Autónomas, son a nivel nacional. Descargado y descomprimido el archivo ZIP correspondiente a 2006 este es su contenido:

- Carpeta con simbología (códigos de colores RGB empleados en la clasificación en formato mdb Access).

- Vectoriales con los polígonos para las Islas Canarias (amarillo).

- Vectoriales con los polígonos para toda España. Canarias incluido (verde).

- Diversos archivos TXT y PDF con información del producto y su metodología, que también contienen los imprescindibles códigos numéricos de clasificación del terreno (Anexo del pdf).

La verdad que este tipo de productos tienen un abordaje complicado, casi siempre debido a la codificación numérica de los diferentes tipos de coberturas. Abrimos el vectorial de toda España -CLC06_es.shp- en QGIS para hacernos la siempre recomendable visión global (542 megas, no es poca cosa):

CORINE Land Cover abierto en QGIS con detalle de zona.

Un vistazo siempre igual de recomendable a la tabla de atributos del archivo vectorial nos proporciona información muy importante; en este caso los campos que contiene dicha tabla son los siguientes:

155770 polígonos en la tabla. Area_ha: extensión en hectáreas del polígono. Shape_leng: longitud del vectorial del polígono (equivalente a su perímetro en metros). Shape_Area: extensión del polígono, esta vez en metros cuadrados. code_06: código identificativo de la ocupación del suelo. ID: identificador único de cada polígono.

Aunque todos los campos tienen su importancia, obviamente el campo code_06 es vital. En él se contiene el código que clasifica la cobertura del suelo para dicho polígono y que podemos encontramos en el anexo del PDF de la descarga (44 clases en total). Clasificamos y categorizamos una zona al azar para verlo mejor:

La utilidad de CORINE Land Cover es evidente, sobre todo a escalas intermedias, y su facilidad de uso en comparación al SIOSE toda una ventaja. Así por ejemplo, en pocos minutos, sobre el sombreado de relieve del Mapa Base del IGN, activamos las coberturas que nos interesen y conseguimos un interesante mapa de ejemplo con abundante y utílisima información:

Coberturas con su código CORINE sobre relieve de una zona de Huesca.

No olvidemos para terminar que este servicio está disponible vía WMS (también WFS), y en su correspondiente capa con su correspondiente leyenda dependiente de la escala de visualización, en la misma URL que el SIOSE: http://www.ign.es/wms-inspire/ocupacion-suelo? y que a través de una petición GetFeatureInfo podemos acceder a los datos básicos de cada polígono. También contiene capas con las variaciones históricas de la cobertura de los polígonos en cada actualización CORINE.

Consulta en Global Mapper de información sobre polígono desde el servicio WMS de CORINE.

Si os interesan los aspectos más técnicos sobre el futuro de los productos relacionados con la ocupación del suelo, se han publicado guías al respecto (final de página bajo el epígrafe Implementación española) de cómo los datos deben ajustarse a la directiva INSPIRE para fomentar la máxima interoperabilidad entre organismos. Saludos y a disfrutar de este extraordinario producto, para mi básico en la elaboración de cartografía topográfica, pues nos permite aportar un útil e interesante contenido extra.