Copernicus FPCUP IGNCNIG
Copernicus FPCUP IGNCNIG
Copernicus FPCUP IGNCNIG
la información geoespacial
Cada uno de los ODS informará sobre la manera en la que cada país y la comunidad internacional mide,
gestiona y monitoriza el progreso y comunicación sobre el crecimiento económico, inclusión social y
sostenibilidad medioambiental como pilares del desarrollo sostenible.
En apoyo de las mediciones y la monitorización del progreso hacia los 17 ODS, las Naciones Unidas
estableció un marco global de 232 indicadores que significan los medios que tienen los gobiernos de las
naciones para monitorizar de manera práctica y reportar hacia las 169 metas de la Agenda 2030.
Figura 3 Reporte del progreso hacia los ODS a partir de los indicadores. Importancia de la
información geográfica y de los datos de observación de la Tierra. Fuente: CEOS
La Agenda 2030 (resolución 70/1 aprobada por la Asamblea General de Naciones Unidas el 25 de
septiembre de 2015), establece en su artículo 76 que:
Ayudaremos a los países en desarrollo, particularmente los países africanos, los países menos adelantados,
los pequeños Estados insulares en desarrollo y los países en desarrollo sin litoral, a fortalecer la capacidad
de sus oficinas de estadística y sus sistemas de datos estadísticos para asegurar el acceso a datos de gran
calidad, oportunos, fiables y desglosados. Promoveremos el aumento, con transparencia y rendición de
cuentas, de una cooperación adecuada entre los sectores público y privado para aprovechar una amplia
gama de datos, incluidos los de observación de la Tierra e información geoespacial, garantizando al mismo
tiempo que sean los propios países quienes se encarguen de apoyar y vigilar los progresos conseguidos.
Por lo tanto, al más alto nivel de las Naciones Unidas se reconoce la importancia que para la monitorización
de los ODS tiene este tipo de datos y servicios.
Para realizar el seguimiento hacia estos objetivos y metas asociadas, el marco de indicadores debe capturar
las ambiciosas y variadas aspiraciones que reflejen el desarrollo continuo de las naciones y las sociedades.
Una adecuada política de informes sobre estos indicadores requerirá el uso de diferentes tipos de datos, tanto
los que ya están disponibles como nuevas fuentes de datos ajenas a los sistemas estadísticos tradicionales, en
concreto datos de observación de la Tierra e información geográfica, empleando técnicas de procesamiento
de datos apropiados para grandes volúmenes de datos, el llamado big data.
En conjunto con lo anterior y en la misma línea, el 27 de julio de 2011 el Consejo Económico y Social de las
Naciones Unidas (ECOSOC) reconoció la necesidad de promover la cooperación internacional en el campo
de la información geoespacial global y decidió, a este respecto, establecer el Comité de Expertos sobre la
Gestión Mundial de la Información Geoespacial de las Naciones Unidas (Global Geospatial Information
Management – UN-GGIM) mediante Resolución ECOSOC 2011/24. En los Términos de Referencia, el
Comité estaría compuesto por expertos de todos los Estados Miembros, así como expertos de organizaciones
internacionales como observadores.
UN-GGIM tiene como objetivo desempeñar un papel de liderazgo en el desarrollo de la información
geoespacial mundial y promover su uso para abordar desafíos globales clave. Proporciona un foro para la
comunicación y la coordinación entre los estados miembros y las organizaciones internacionales. UN-GGIM,
junto con el Banco Mundial, consciente de que la información geoespacial muestra dónde se están
produciendo las circunstancias sociales, medioambientales y económicas, se esfuerza en establecer un marco
mediante el cual una información geográfica fiable y acreditada sirva para medir y gestionar los ODS para
mejorar la vida de las personas y proteger el planeta. Esta información, según UN-GGIM, integra todos los
datos digitales que tienen un componente de localización, desde características topográficas básicas que
ofrecen los mapas hasta conjuntos de datos complejos en múltiples capas.
Este apartado se tratará con más detalle en otro punto del curso, para poder tratar más en concreto la
contribución de la teledetección a los ODS.
Un resumen visual para ver en qué ODS la observación de la Tierra tiene más peso puede verse en la Figura
6.
Tabla 1 Desglose de algunos ODS, metas e indicadores asociados donde la observación de la Tierra
tiene un gran potencial. Fuente: CEOS
Otro de los aspectos fundamentales es el hecho de que las agencias cartográficas y estadísticas viven en los
últimos años unas transformaciones importantes debidas a las nuevas tecnologías y las nuevas metodologías
de adquisición y depuración de datos. Existe una tentación a la hora de adoptar nuevos métodos de descartar
lo anterior por poco preciso. Es fundamental que estas instituciones desarrollen metodologías para garantizar
la continuidad a lo largo de los años y décadas para garantizar la continuidad en la serie de datos. Esto
también aplica a las agencias espaciales a la hora de diseñar las nuevas misiones, las cuales deben garantizar
continuidad con las anteriores y evitar, en la medida de lo posible, periodos temporales con ausencia de
datos.
Tabla 2 Contribuciones significativas a los ODS por parte del programa espacial europeo Copernicus.
Fuente: UNOOSA
Por ejemplo: incendios forestales, inundaciones, vertidos de petróleo, terremotos…etc. Son aplicaciones en
las que requieren imágenes con alta frecuencia temporal en el periodo de tiempo que se están produciendo el
fenómeno. En estos casos es donde más se utilizan imágenes radar, ya que pueden tomar imágenes en
cualquier momento del día o de la noche y en condiciones meteorológicas. Para estos casos se buscan
también plataformas satelitales que posean “constelaciones” de satélites, iguales o complementarios, para
recoger la información demandada con la frecuencia requerida.
Figura 2 Imagen del sensor MODIS. Incendios declarados en la Comunidad Valenciana, España
La teledetección permite clasificar los tipos de cultivos que aparecen en las imágenes y ayuda a conocer el
estado fenológico de cada cultivo mediante las firmas espectrales, por lo que puede ser utilizada para:
Obtener estadísticas de superficies de cultivos.
Evolución de superficies ocupadas por cultivos.
Seguimiento del estado vegetativo de los cultivos y detección de anomalías.
Gestión del agua de riego y control de su aplicación.
Agricultura de precisión: predicción del rendimiento de los cultivos y del momento óptimo para su
cosecha.
La teledetección permite clasificar los tipos de suelo o coberturas vegetales que aparecen en las imágenes,
por lo que puede ser utilizada para:
Seguimiento de especies y poblaciones animales.
Realizar cartografía e inventarios de la cobertura vegetal del suelo.
Evaluar las condiciones de estrés en la vegetación por los efectos de la sequía o la deforestación.
La siguiente imagen es el llamado “índice de vegetación” (NDVI, Normalized Difference Vegetation Index)
obtenido para toda la Tierra a partir de imágenes del sensor MODIS, con 1 Km de tamaño de píxel. A golpe
de vista permite hacerse una buena idea de dónde se encuentran localizadas las zonas de mayor vegetación
en el mundo y donde los grandes desiertos.
Figura 7 Ventanas atmosféricas y su relación con el máximo de energía emitida por el Sol
La cantidad de energía que es reflejada, transmitida o absorbida por la superficie terrestre depende de la
longitud de onda de la radiación incidente y de las características de cada cubierta (propiedades químicas,
rugosidad, etc.).
La siguiente ilustración representa las distintas técnicas existentes en teledetección según el tipo de radiación
capturada por los sensores:
Hay que aclarar que los sensores a bordo de satélites no miden directamente la radiación reflejada por las
cubiertas de la superficie terrestre, puesto que cuando llega al sensor se ha visto afectada por la atmósfera
que modifica los valores originales. La energía reflejada sufre procesos de dispersión y absorción al
interactuar con los componentes gaseosos y los aerosoles de la atmósfera.
b) Órbitas heliosíncronas
Son órbitas polares o cuasi polares (casi pasan por los polos). Muchos satélites de teledetección están
equipados con sensores pasivos que dependen de la iluminación solar. Al ir midiendo la reflexión de la luz
solar procedente de la Tierra se deben ajustar sus órbitas al ritmo del día y de la noche. Esto es importante
para poder comparar imágenes registradas a lo largo de un periodo de tiempo, para lo que las condiciones de
luz deben ser lo más parecidas posible. Los registros deben tener lugar a la misma hora local del día para que
el plano de la órbita del satélite debe permanecer a un ángulo constante a la luz del Sol. Estos requisitos
pueden cumplirse situando el satélite en una órbita heliosíncrona.
Mientras el satélite avanza en su órbita, la Tierra gira sobre su eje. Cada vez que el satélite completa una
vuelta se escanea una nueva franja de la superficie de la Tierra y, pasado un cierto número de vueltas, se
habrá obtenido toda la superficie de la Tierra. Algunos satélites escanean una franja ancha cada vez y pueden
de este modo cubrir la totalidad de la superficie de la Tierra en unas pocas vueltas. Por contra, los satélites de
alta resolución que escanean sólo tiras finas tardan varios días en completar la cobertura de la Tierra.
Los sensores pasivos se clasifican según del número de bandas o canales espectrales de las imágenes
adquiridas:
Sensores pancromáticos: capturan información en una sola banda, normalmente en el rango entre la
radiación visible y el infrarrojo próximo. Las imágenes obtenidas se visualizan en escala de grises.
Sensores multiespectrales: capturan información en varias bandas, asignando una banda para cada
rango del espectro electromagnético del que pueden capturar información.
Sensores hiperespectrales: cuando el número de bandas en el que capturan información es muy
elevado, llegando a varias centenas, se habla de sensores hiperespectrales. Hoy en día son los menos
Los sensores activos se clasifican dependiendo de la longitud de onda en que trabajan. Para trabajos de
teledetección se distinguen dos tipos de sensores activos:
Sensores RADAR (RAdio Detection And Ranging) en el rango de las microondas.
Sensores LiDAR (LIght Detection and Ranging) emite pulsos de luz láser entre el ultravioleta y el
infrarrojo cercano.
Los sensores LIDAR permiten obtener mayor detalle en la observación que los sensores RADAR porque
trabajan con longitudes de onda más pequeñas pero, por el contrario, el láser no puede atravesar las nubes.
Los sensores LiDAR no se suelen utilizar en teledetección espacial, y normalmente se colocan a bordo de
aviones.
Los sistemas radar se verán con detalle más adelante en el tema de teledetección RADAR.
Los Landsat son una serie de satélites construidos y puestos en órbita por Estados Unidos para la
monitorización de la superficie terrestre. Todos con órbita circular heliosincrónica, a 705 km de altura,y con
una inclinación de 98.2º respecto del Ecuador. La órbita de los satélites está diseñada de tal modo que cada
vez que éstos cruzan el Ecuador lo hacen de Norte a Sur sobre las diez de la mañana hora local. Los Landsat
están equipados con sensores específicos para la teledetección multiespectral en el rango óptico del espectro
electromagnético.
El Landsat 1 se lanzó al espacio en 1972; el siguiente el Landsat 2 en 1975, el Landsat 3 en 1978, el Landsat
4 en 1982, el Landsat 5 en 1985, el Landsat 6 tuvo el lanzamiento fallido en 1993, el Landsat 7 en 1999, y el
último de la serie por ahora, el Landsat 8, fue lanzado en febrero de 2013. Los dos últimos siguen en
funcionamiento, adquiriendo imágenes con un tamaño de pixel de 30 metros (más un pancromático de 15 m)
y un periodo de revisita de 16 días cada satélite y desfasados 8 días entre ellos.
La NASA mantiene desde hace varios años una política de distribución gratuita de imágenes. Se pueden
descargar a través del catálogo del USGS: http://earthexplorer.usgs.gov/.
Los Spot son una serie de satélites desarrollados por el CNES (Centro Nacional de Estudios Espaciales de
Francia) para la monitorización de la superficie terrestre. Su órbita alrededor de la Tierra es heliosíncrona,
casi polar y circular, a 822 km de altura.
El Spot 1 se lanzó en 1986; Spot 5 en 2002, Spot 6 en 2012 y el último lanzado fue el Spot 7 en junio de
2014. Estos dos últimos siguen en funcionamiento, y conjuntamente permiten la observación diaria de
prácticamente cualquier punto del globo terráqueo con una resolución de 1.5 metros.
El satélite SPOT 5, que se mantuvo en funcionamiento desde 2002 a 2015, capturaba 2 imágenes
pancromáticas simultáneas con una resolución de 5 metros que convertían a 2,5 metros en postproceso e
imágenes de 10 metros de resolución espacial con 3 bandas espectrales: Verde, Rojo, NIR, además de un
SWIR de 20 metros.
Los satélites SPOT6 y 7 garantizan la continuidad de los productos de alta resolución hasta el año 2024,
capturando imágenes pancromáticas con 1,5 metros de resolución espacial y multiespectrales con 6 metros y
4 bandas espectrales: Azul, Verde, Rojo, NIR.
Los formatos de fichero más utilizados para almacenar imágenes de satélite (imágenes raster multibanda que
requieren tamaños grandes) son: los que se especifican a continuación:
Formato TIFF: permite almacenar imágenes raster de hasta 4 GB de espacio de almacenamiento; y
hasta 32 bits por píxel. Las imágenes se pueden almacenar sin compresión, comprimidas sin
pérdidas (algoritmo LZW) o comprimidas con pérdidas (algoritmos: PackBits / CCITT G3y4 / RLE
/ JPEG / UIT-T). En teledetección se suelen almacenar las imágenes sin comprimir. La información
sobre la georreferenciación de la imagen se almacena en un fichero de texto aparte con extensión
TFW y otro con extensión .PRJ.
Formato GeoTIFF: el mismo que el anterior, pero almacena la información de la
georrefereanciación de la imagen en la cabecera del fichero.
Formato BigTIFF: es una versión mejorada del formato TIFF. Permite almacenar imágenes
mayores de 4 GB (hasta 143 Gb). Existe la variante GeoBigTIFF.
Formato ECW: comprime las imágenes con poca pérdida de calidad. Reduce considerablemente el
tamaño de los archivos, manteniendo una aceptable calidad visual y permitiendo una rápida
compresión y descompresión. Utiliza un algoritmo “quad-tree”, que subdivide sucesivamente la
imagen en cuadrantes. La información de la georreferenciación de la imagen se puede almacenar en
la cabecera del archivo o en un fichero de texto aparte con extensión ERS.
Formato IMG: Formato utilizado por el programa ERDAS para almacenar imágenes raster, pero
bastante extendido y que pueden leer muchos otros programas.
Formato JPEG: Técnica estándar de compresión para el almacenamiento de imágenes a color y en
escala de grises. Es uno de los formatos más conocidos. Es uno de los pocos formatos que se
soporta en Internet (Web). Almacena las imágenes comprimidas con pérdidas. Se puede escoger
diferentes niveles de compresión: a más alta compresión menor calidad. Se puede optar por una
compresión alta para crear archivos que ocupen poco espacio, que se usarán en la Web o para ser
enviadas por correo electrónico. JPEG puede llegar a comprimir una imagen hasta sólo un 10 % de
su tamaño original, sin que el ojo humano pueda percibir grandes diferencias. El teledetección se
utiliza mayoritariamente para almacenar “quicklooks” o para publicar una imagen de satélite en
Las imágenes multiespectrales disponen de varias bandas donde se almacena información de diversas
regiones del espectro electromagnético.
Los monitores y tarjetas de video disponen de 3 canales RGB: rojo, verde y azul. Por tanto puede utilizarse
cada canal del monitor para representar los ND de una banda y obtener así una composición de color.
Por tanto, los datos multiespectrales se pueden representar por medio de composiciones de color o
combinaciones de tres bandas.
En general, se trata de aprovechar que podemos visualizar tres canales a la vez para introducir las tres bandas
que más nos van a ayudar a discriminar visualmente los elementos que nos interesan. Estas composiciones
servirán entonces para resaltar los elementos que mayor variación de reflectividad presentan en las bandas
utilizadas, además de obtener visualizaciones más o menos estéticas.
La denominación de las combinaciones de color puede variar de un autor a otro, pero generalmente se llaman
“color natural” y “falso color infrarrojo” a dos de las combinaciones más utilizadas:
Color natural consiste en combinar las bandas rojo, verde y azul del sensor, sacándolas por los
canales RGB del monitor respectivamente.
Falso color infrarrojo: consiste en combinar las bandas NIR, rojo y verde del sensor, enviándolas a
los cañones R, G y B del monitor respectivamente. En esta composición de color la vegetación se
visualiza en tonalidades de rojo y se utiliza para discriminar con mayor facilidad distintos tipos de
vegetación, ya que en la banda 5 la vegetación posee alta reflectividad.
A continuación, pueden explicarse las distintas combinaciones de bandas utilizando una imagen del satélite
Landsat 8.
Figura 20 Combinaciones color natural (4-3-2) y falso color (6-5-4) de una imagen Landsat
En el siguiente cuadro se listan las combinaciones de bandas para los sensores más habituales en
teledetección de recursos terrestres:
Espectro electromagnético
Satélite Coastal Azul Verde Rojo NIR SWIR 1 SWIR 2
Landsat 8 1 2 3 4 5 6 7
Sentinel 2 1 2 3 4 8 11 12
Teniendo en cuenta la información de la Tabla 1, la Tabla 2 nos proporciona la información necesaria para
poder visualizar una misma imagen Landsat 8 o Sentinel 2 en diferentes combinaciones de bandas,
Tabla 2 Configuración de los cañones del monitor para visualizar diferentes combinaciones de color
con Sentinel 2 y Landsat 8
Cuando se trabaja con varias imágenes que han sido obtenidas en fechas distintas, poseerán una radiometría
distinta debido a la variación en las condiciones atmosféricas, y esto hace que el aspecto y los colores no
sean iguales visualmente para todas las imágenes. Los equilibrados radiométricos de realizan para igualar la
radiometría de varias imágenes con una de referencia.
Se pueden consideran como un paso previo a la obtención de mosaicos de imágenes, aunque también se
realizan equilibrados para calibrar series temporales de imágenes o para otras aplicaciones. El equilibrado
radiométrico consiste en igualar las estadísticas (media y desviación típica) de las imágenes para que todas
tengan un aspecto visual similar. Una de las imágenes se tomará como referencia para igualar la radiometría
de las demás a la de esta.
El método más sencillo para igualar las estadísticas de una imagen a otra, y que funciona bastante bien, es
mediante una transformación lineal. Para ello las imágenes a equilibrar deben poseer una zona de solape o
recubrimiento de donde se obtienen las estadísticas, media y desviación típica, de ambas imágenes para cada
banda y se calculan los coeficientes de la recta de transformación. Posteriormente se aplican las ecuaciones
de esa recta a todos los ND de los píxeles de la imagen a equilibrar para igualar su radiometría con la otra.
4.4.2 Mosaicos
Es frecuente que las escenas de satélite no cubran la totalidad del área de estudio. En este caso se pueden
adquirir escenas adyacentes y realizar un mosaico con todas ellas de forma que quede una imagen continua.
Para ello habrá que ajustar previamente la radiomatería para que todas las imágenes tengan el mismo aspecto
y a la hora de mosaicar no aparezcan saltos de color, este proceso conocido como equilibrado radiométrico
es el que se ha explicado en un apartado anterior.
El realce se aplica para producir un mejor aspecto visual, mediante una transformación del histograma de la
imagen, que habitualmente implica mejorar el contraste de la imagen. Los realces se realizan píxel a píxel sin
tener en cuenta los de alrededor. La modificación píxel a píxel, significará que para un valor ND de entrada
le corresponderá un NDmodificado de salida.
El histograma de una imagen es la representación gráfica (o en tabla) del número de píxeles asignado a cada
nivel de gris o valor digital (ND) distinto en la imagen. La distribución de Niveles Digitales representada en
el histograma aporta una información a partir de la cual pueden definirse estas operaciones de realce a nivel
de píxel.
En la práctica, los realces se aplican modificando el histograma de las imágenes. La modificación o
transformación más empleada en la expansión lineal del histograma para mejorar el contrate de la imagen.
Figura 26 Ejemplo de expansión del histograma en una imagen con 3 bandas RGB
Los realces son de fácil implementación por los programas de tratamiento de imágenes, confeccionándose
tablas de codificación o Look Up Tables (LUT), donde se establecen las transformaciones adecuadas.
Las siguientes imágenes de la isla de Mallorca, muestran el resultado de aplicar el proceso de realce anterior:
4.4.4 Filtros
Los filtros se utilizan básicamente para destacar algunos elementos de la imagen. Consiste en aplicar a la
matriz de pixeles de la imagen una matriz de filtrado “deslizante” (por ejemplo de 3x3 o de 5x5) que genera
un nuevo valor mediante una media ponderada del valor original y el valor de los pixeles circundantes.
Mediante diferentes combinaciones de los factores de ponderación se pueden conseguir diferentes efectos.
Los filtros más utilizados son:
Los filtros de paso bajo que asignan a cada pixel el valor medio de los pixeles circundantes: efecto
suavizado.
Los filtros de paso alto para resaltar zonas de gran variabilidad: efecto realce de bordes
Existen muchos satélites de observación de la Tierra que llevan a bordo un sensor pancromático
multiespectral, que capturan la misma imagen de la Tierra en el mismo instante, de forma que la banda
pancromática es perfectamente superponible o coincidente con las bandas multiespectrales de la imagen. Por
condiciones de diseño, los sensores pancromáticos poseen siempre mayor resolución espacial que los
multiespectrales, es decir que proporcionan imágenes con menor tamaño de píxel, pero por el contrario
capturan información en una sola banda que abarca un rango espectral más amplio que las multiespectrales.
La fusión o pan-sharpening es un producto derivado de las imágenes de satélite. Consiste en combinar las
bandas multiespectrales de una imagen con su banda pancromática (capturada en el mismo instante que las
multiespectrales), para mejorar la resolución espacial de las bandas multiespectrales extrayendo la
información espacial que les falta de la banda pancromática.
En el siguiente ejemplo se puede ver cómo mejora la calidad espacial de las bandas multiespectrales de una
imagen SPOT al fusionarlas con la banda pancromática:
El componente espacio del programa está liderado técnicamente por la ESA, y a su vez se subdivide en 2
elementos, dependiendo de la propiedad de los satélites empleados en la captura de datos:
1. Las misiones Sentinel, operadas de manera conjunta por la propia ESA y la Organización Europea
para la Explotación de Satélites Meteorológicos (EUMETSAT).
2. Las denominadas como misiones contributivas, las cuales incluyen misiones propias de la ESA,
EUMETSAT y otros proveedores privados de datos que ponen a disposición de Copernicus una
amplia variedad de información satelital.
En la actualidad existen alrededor de 30 misiones contributivas ya existentes o previstas, las cuales
se incluyen en las siguientes categorías:
Radares de apertura sintética (SAR) para observar los océanos tanto de día como de noche.
Sensores ópticos para monitorizar actividades sobre el territorio y dinámicas oceánicas.
Altímetros para medir el nivel del mar.
Radiómetros para monitorizar la temperatura de las superficies emergidas y de los océanos.
Sentinel-2A y Sentinel-2B ocupan la misma órbita, con un desfase de 180º. La altitud media de la órbita
sobre la superficie terrestre es de 786 km. La inclinación de la órbita es de 98,62° y la Hora Solar Local
Media (MLST, Mean Local Solar Time) en el nodo descendente se produce a las 10:30 (am). Este valor de
MLST se eligió como un compromiso entre un nivel adecuado de iluminación solar y la minimización de la
cobertura potencial de nubes. El valor de MLST es cercano al tiempo de paso local de Landsat y casi
idéntico al de SPOT-5, lo que facilita la integración de los datos de Sentinel-2 con estas misiones, existentes
e históricas, y contribuye a la recopilación de datos de series temporales.
El archivo de datos KML que informa sobre la posición de las órbitas relativas Sentinel-2A y Sentinel-2B
para un ciclo completo (143 órbitas) con un paso de tiempo de 10 segundos está disponible en formato kml
en el enlace Sentinel-2 Relative Orbits.
2.4 Resoluciones
En este apartado se especifican todas las resoluciones de la misión Sentinel-2: temporal, radiométrica,
espectral y espacial.
La resolución temporal de cada satélite de la misión Sentinel-2 por separado es de 10 días, duplicándose la
frecuencia al estar los dos operativos, alcanzándose por tanto una resolución temporal de 5 días bajo las
mismas condiciones de iluminación, salvo para las regiones únicamente observadas por un satélite, tal y
como se indica en la Figura 3.
Figura 3 Cobertura y tiempo de revisita previsto para las adquisiciones MSI de Sentinel 2
En la Figura 17 muestra las regiones donde se incrementa la frecuencia gracias a la superposición entre las
franjas de las órbitas adyacentes, si bien cambian las condiciones de la geometría de observación.
La resolución espacial de cada una de las bandas y que se resume en las siguientes líneas y en la Figura 5:
4 bandas a 10 m: banda 2, banda 3, banda 4 y banda 8.
6 bandas a 20 m: banda 5, banda 6, banda 7, banda 8a, banda 11 y banda 12.
3 bandas a 60 m: banda 1, banda 9 y banda 10.
La resolución radiométrica del instrumento MSI es de 12 bits, lo que permite adquirir la imagen en un rango
de 0 a 4095 valores de intensidad de luz potencial.
Los datos de Sentinel-2 se adquieren en 13 bandas espectrales tal y como se muestra en la Tabla 1.
2.4.5 Aplicaciones
En este apartado se describe la aportación de la misión Sentinel-2 a varios servicios del programa Copernicus
y se describen algunas aplicaciones realizadas en España.
A través de la información suministrada por sus satélites, Sentinel-2 aporta información útil para prácticas
agrícolas y forestales, la gestión de la seguridad alimentaria, la supervisión del estado de la superficie
terrestre y control de cambios, alertas producidas por la contaminación en lagos y aguas costeras, soporte en
tareas de ayuda humanitaria y gestión de consecuencias de desastres naturales, tales como inundaciones,
erupciones volcánicas y deslizamientos.
a) Gestión del territorio: Los datos adquiridos, la amplia cobertura de la misión y la alta frecuencia de
revisita de Sentinel-2 permiten el suministro de información geográfica a escala local, regional,
nacional e internacional. Estos datos son complementarios a los aportados por otras misiones y están
diseñados para ser modificados y adaptados por usuarios interesados en áreas temáticas tales como:
Supervisión del recurso agua: La utilización de datos de observación de la Tierra de las
misiones Sentinel-2 y su procesamiento rápido de productos derivados permite obtener
resultados rápidos de una forma muy sencilla. La alta resolución temporal de los dos
satélites Sentinel-2A y Sentinel-2B permite inspeccionar posibles incidencias cuando el
cultivo está implantado, incluso con evidencias de haber sido regado recientemente.
Figura 7 Evolución de NDVI a lo largo de una campaña a partir de Sentinel 2. Fuente: Instituto de
Desarrollo Regional, Castilla la Mancha, España
En este caso, a partir del cálculo de determinados índices de vegetación derivados de las
imágenes de satélite se pretende generar mapas que se correlacionen con el rendimiento real.
Figura 9 Mapas de valoración de daños del incendio que se inició el 24 de junio 2017 en Moguer
(Huelva). Rapid Mapping. Fuente: Copernicus EMS
Figura 11 Mapa de la nube tóxica producida por el accidente industrial del 12 de febrero de 2015 en
Igualada (Barcelona). Risk and Recovery Mapping. Fuente: EMS
En este apartado se describen algunas convenciones utilizadas en Sentinel. En este documento se ha decidido
mantener la denominación original utilizada en los productos Sentinel. De este modo, en este marco, se
denomina datatake a la franja de adquisición continua de imágenes de Sentinel-2 en un modo MSI
determinado y que forma la base de los productos que se estudiarán posteriormente. La longitud máxima de
un datatake de imágenes es de 15.000 km (por ejemplo, observación continua desde el norte de Rusia hasta
el sur de África). Si un datatake es adquirido en dos estaciones terrestres receptoras separadas, el datatake
puede subdividirse en datastrips.
A continuación, se establecen las unidades en las que se dividen los productos Seninel-2 y que guardan
relación directa con los distintos niveles de procesamiento de cada producto.
Los productos son una compilación de unidades elementales de tamaño fijo, dentro una única órbita. Existen
unidades mínimas indivisibles en cada producto que contiene todas las bandas espectrales posibles. El
tamaño de esta unidad depende del nivel del producto.
Proyecto cofinanciado por la Comisión Europea mediante acuerdo 2018/SI2.810140/04 15
Para productos de Nivel-0, Nivel-1A y Nivel-1B: las unidades de división se denominan gránulos,
que hacen referencia a una subimagen de un número determinado de líneas a lo largo del recorrido
del satélite. Tienen un tamaño de 25 km en disposición transversal al recorrido y 23 km en
disposición longitudinal.
Para productos ortorectificados de Nivel-1C y Nivel-2A: las unidades de división se denominan tiles. Una
tile tiene una extensión geográfica de 100 km x 100 km en el sistema de referencia geodésico WGS84 y en la
proyección cartográfica UTM, en el huso correspondiente a la ubicación geográfica. Se puede descargar la
malla de tiles a nivel mundial en formato kml (Figura 14) en el siguiente enlace de la ESA:
https://sentinel.esa.int/web/sentinel/missions/sentinel-2/data-products.
2.5.2 Nivel-1C
Dato Información
Metadatos A nivel de producto y tiles
Datos de imagen Una imagen en reflectancias TOA por cada banda espectral (JPEG2000)
Derivados del Nivel-1B y Nivel-1C: calidad radiométrica, calidad
Datos indicadores de calidad geométrica, calidad del contenido de las imágenes, información de la
verificación de los controles de calidad.
Referencias GIPP, DEM y GRI utilizadas.
Parámetros meteorológicos ECMWF: Total Column Ozone (TCO3),
Datos auxiliares
Total Column Water Vapour (TCWV) y Mean Sea Level Atmospheric
Pressure (MSLP)
2.5.3 Nivel-2A
Desde mediados de marzo de 2018, el Nivel-2A se convirtió en un producto operacional, comenzando con la
cobertura de la región euromediterránea, con un aumento gradual a la cobertura mundial sistemática en
diciembre de 2018. Además, el Nivel-2A también puede ser generado por el usuario a partir del producto
Nivel-1C utilizando la caja de herramientas de Sentinel-2 o la versión independiente del procesador
Sen2Cor.
El Nivel-2A de Sentinel-2, L2A Product Definition Document (Louis and Devignot, 2018), es una
ortoimagen expresada en valores de reflectancias a nivel de la superficie (BOA), en tanto por unidad, con
valores entre 0 y 1, con una resolución radiométrica de 12 bits. Los píxeles defectuosos, sin datos y saturados
se excluyen de los pasos de procesamiento del Nivel-2A. La reflectancia de la superficie se realiza a la
resolución de los productos nativos en función de la resolución de la banda (10, 20 y 60 m). Los algoritmos
empleados en el proceso de generación del Nivel-2A requieren información derivada de bandas con
diferentes resoluciones. Algunos canales son remuestreados a 60 m según las necesidades de los algoritmos.
Estas imágenes se encuentran corregidas atendiendo a los parámetros relativos a la atmósfera (visibilidad,
vapor de agua, presencia de nubes o cirros y contenido de aerosoles y ozono), a la geometría de la
iluminación-observación (ángulo cenital solar, ángulo acimutal solar y topografía del terreno) y, previa
activación en las opciones del archivo de cabecera del algoritmo, del comportamiento reflectivo que tienen
las distintas cubiertas (Bidirectional Reflectance Distribution Function, BRDF). El modelo atmosférico que
utiliza Sen2Cor se basa en el modelo ATCOR (Atmospheric and Topographic CORrection), el cual se
fundamenta en interpretar los parámetros citados con anterioridad mediante tablas de transformación (Look
Up Tables, LUT) que, a su vez, dependen del cálculo de funciones de transferencia radiativa para: diferentes
sensores, geometrías solares, topografía del terreno y parámetros atmosféricos (Mueller-Wilm, 2016)
Figura 15 Ráster Aerosol Optical Thickness (izquierda); ráster Water Vapour (derecha)
A partir del 6 de diciembre de 2016 se distribuye en un nuevo formato para todos los productos de Nivel-1C.
Este formato se introdujo para superar la limitación de 256 caracteres en las rutas de acceso impuestas por
las plataformas del sistema operativo Windows. Esto se ha logrado mediante la compactación de los nombres
de archivo, incluida la denominación de carpetas y archivos internos. Los productos de Nivel-2A ya se
empezaron a distribuir en este formato.
El formato actualizado del producto se describe en detalle en el Documento de especificación de productos
(PSD, ProductsSpecificationDocument de Sentinel-2, versión 14.5).
La convención de nomenclatura compacta es la siguiente:
MMM_MSILNN_YYYYMMDDHHMMSS_Nxxyy_ROOO_Txxxxx_ <Product Discriminator> .SAFE
Los productos contienen dos fechas:
La primera fecha (YYYYMMDDHHMMSS) corresponde a la adquisición.
La segunda fecha es el campo "<Product Discriminator>", que tiene una longitud de 15 caracteres
y se utiliza para distinguir entre diferentes productos para usuarios finales.
Los otros componentes del nombre de archivo son:
MMM: es el ID de misión y satélite (S2A / S2B).
MSILNN: denota el nivel del producto nivel (1C/2A).
YYYYMMDDHHMMSS: la hora de inicio de la detección del Datatake.
Nxxyy: el número de línea de base de procesamiento.
ROOO: Número de órbita relativa (R001 - R143).
Txxxxx: Código de identificación de la tile.
SAFE: Formato de producto (Formato de archivo estándar para Europa)
A modo de ejemplo:
S2A_MSIL1C_20180206T111231_N0206_R137_T30TUL_20180206T132826.SAFE
Se identifica:
Satélite: Sentinel 2A.
Producto Nivel 1C.
Fecha de adquisición: 2018/02/06 hora: 11h 12min 31seg.
Procesado con Baseline N0206.
Órbita: R137.
Tile: T30TUL.
Fecha de generación del producto: 2018/02/06 13h 28m 26seg.
El producto Nivel-2A tiene la siguiente estructura de carpetas. Por ejemplo, para el producto:
S2B_MSIL2A_20180906T105639_N0208_R094_T30TVL_20180906T155143.SAFE
Los productos Sentinel se distribuyen en formato SAFE standard (SAFE Archive Format for Europe) e
incluye la siguiente información:
Archivos de metadatos en formato XML. Los metadatos del producto se encuentran en la carpeta
raíz. El fichero se denomina MTD_MSILNP.xml donde NN es el nivel de procesamiento. También
se incluye el metadato del producto de acuerdo con las especificaciones INSPIRE y se denomina
INSPIRE.XML. Se muestra un ejemplo de producto 1C (Figura 19) y 2A (Figura 20).
S2B_MSIL1C_20180820T110609_N0206_R137_T30TUL_20180820T154632.SAFE
En esta carpeta se incluyen archivos que expresan la calidad en forma máscaras en formato GML:
o Máscaras que hacen referencia a la información radiométrica de las imágenes, como por
ejemplo: píxeles defectuosos (MASK_DEFECT), saturados (MASK_SATURA), sin
información (MASK_NODATA), etc.
o Máscaras de nubes (MSK_CLOUDS), clasificadas en densas y cirros, como se muestra en la
Figura 23.
En el producto 2A se incluyen: máscaras de probabilidad de nubes (MSK_CLDPRB), máscaras de
probabilidad de nieves (MSK_SNWPRB) y, de forma opcional, máscaras de vegetación densa
(MSK_DDVPXL).
Huellas de los detectores para cada banda (MSK_DETFOO).
En la carpeta AUX_DATA se encuentran los parámetros meteorológicos y físicos
(AUX_ECMWFT) y, de forma opcional con los procesamientos de la Toolbox, el DEM a la
resolución de la imagen. Además, puede contener ficheros GIPP y otros datos auxiliares.
En este apartado se describen los cuatro modos de adquisición, mostrados en la Figura 24, exclusivos en los
que opera el instrumento SAR de Sentinel-1:
Stripmap (SM).
Interpherometric Wide swath (IW).
Extra-Wide swath (EW).
Wave (WV).
El haz de antena de radar de banda C Sentinel-1 ilumina el terreno hacia el lado derecho del satélite. Debido
al movimiento del satélite y al ancho del haz de la antena a lo largo de la trayectoria (azimut), cada objeto
observado sólo permanece dentro del haz de iluminación durante un corto período de tiempo. Como parte del
procesamiento en tierra, las complejas señales de eco recibidas durante este tiempo se agregan de manera
coherente. De esta forma, se sintetiza una antena larga, con la longitud de apertura sintética igual a la
distancia recorrida por el satélite durante el tiempo de integración. En principio, la resolución geométrica en
azimut que se puede obtener con SAR es aproximadamente la mitad de la longitud real de la antena.
Para mejorar la resolución radiométrica se puede emplear el procesamiento multi-look de azimut, lo que
implica una reducción en la resolución geométrica en azimut en un factor igual al número de looks. La
resolución geométrica en distancia (range resolution) es una función del ancho de banda de la señal radar
transmitida.
El hecho de que el sistema funcione de forma coherente de extremo a extremo significa que tanto la relación
de amplitud como la de fase entre las señales complejas transmitidas y recibidas se mantienen a lo largo de
los instrumentos y la cadena de procesamiento. Esto facilita la síntesis de apertura, así como la
interferometría de radar de múltiples pasadas utilizando pares de imágenes tomadas sobre la misma área en
diferentes momentos.
Las características de los diferentes modos de adquisición se explican a continuación:
a) Modo de franja o StripMap (SM): Este modo está previsto para dar continuidad a las misiones ERS
y Envisat. En este modo, el instrumento proporciona cobertura interrumpida con una resolución
ligeramente mejor de 5 m x 5 m en un ancho de swath de 80 km. Seis pulsos generan seis swaths
superpuestos que cubren un rango de 375 km. Para cada swath la antena está configurada para
generar un haz con un azimut fijo y un ángulo de incidencia diferente.
b) Modo de franja de cobertura interferométrica o Interpherometric Wide Swath (IW): Es el modo
de adquisición predeterminado y principal de Sentinel-1 sobre tierra con resoluciones geométrica y
radiométrica adecuadas para la mayoría de las aplicaciones. Los datos se adquieren en tres sub-
swaths usando la Técnica de Observación de la Tierra con Escaneo Progresivo (TOPSAR, Terrain
Característica SM IW EW WV
Ancho swath 80 km 250 km 410 km -
Cobertura en terreno de una - -
- 20 x 20 km
vignette
Dual HH+HV, Dual HH+HV, Dual HH+HV,
VV+VH VV+VH VV+VH Single HH,
Opciones de polarización
VV
Single HH, VV Single HH, VV Single HH, VV
Resolución espacial 5x5m 5 x 20 m 20 x 40 m 5x5m
En relación con las aplicaciones de los productos derivados, los modos principales sin conflicto de Sentinel-
1 son:
Supervisión de la tierra: modo IW, con polarización VV+VH.
Supervisión mar abierto: modo WV, con polarización VV.
El modo EW se usa principalmente para la supervisión de zonas amplías de costa en mares europeos
seleccionados, áreas del Ártico y del océano Austral, en aplicaciones relacionadas con el tráfico marítimo,
posibles derrames de petróleo y para la supervisión del hielo marino.
El modo SM solo se usa para islas pequeñas y bajo petición, para eventos extraordinarios como la gestión de
emergencias.
Disponer del modo IW como el modo operativo principal sobre tierra, satisface la mayoría de los requisitos
de servicio actuales, evita conflictos y preserva el rendimiento de visitas, simplifica la planificación de
misiones, reduce los costos operativos y crea un archivo coherente a largo plazo.
En la Figura 25 se muestra el escenario de observación Sentinel-1 en términos de modo de adquisición,
polarización y geometría de observación.
3.3.1 Planificación
Este apartado describe la frecuencia de revisita y la cobertura proporcionada por los modos principales de
Sentinel-1. La constelación Sentinel-1 ofrece una mejora en el tiempo de revisita respecto a las misiones
ERS-1/2 y ENVISAT ASAR y una continuidad en la cobertura de ENVISAT ASAR, aportando una mayor
resolución espacial y una cobertura de doble polarización potencialmente global sobre masas terrestres.
Sentinel-1A y Sentinel-1B comparten el mismo plano orbital, casi polar, sincrónico con el Sol, con una
diferencia de fase orbital de 180°, con un ciclo de repetición de 12 días y 175 órbitas por ciclo.
Por tanto, un solo satélite Sentinel-1 es potencialmente capaz de cartografiar masas terrestres globales en el
modo IW una vez cada 12 días, en una sola pasada (ascendente o descendente). La constelación combinada
de los dos satélites ofrece por tanto un ciclo de repetición exacto de 6 días en el ecuador. Dado que el
espaciamiento de la trayectoria de la órbita varía con la latitud, la tasa de revisita es significativamente
mayor en las latitudes más altas que en el ecuador (Figura 121).
Los productos Nivel-1 SLC son datos SAR georeferenciados utilizando información de la órbita y de la
actitud del satélite.
Cada píxel de la imagen está representado por un valor de magnitud compleja (I y Q) y, por lo tanto,
contiene información de amplitud y fase. Cada valor I y Q es de 16 bits por píxel. El procesamiento de todos
los productos SLC produce como resultado una única vista (single look) en cada dimensión utilizando el
ancho de banda completo de la señal disponible.
Los productos SLC se generan para todos los modos de adquisición:
Los productos Nivel-1 GRD consisten en datos SAR proyectados utilizando un modelo elipsoidal WGS84.
Los píxeles incluyen el valor de amplitud, perdiéndose la información de la fase. El producto original tiene
píxeles aproximadamente cuadrados. Tras la aplicación del procesamiento multi-look para la reducción del
efecto de moteado o efecto “sal y pimienta” (speckle) se obtendrán píxeles cuadrados con pérdida de
resolución.
Los productos de GRD están disponibles en tres resoluciones, que se caracterizan por el modo de
adquisición: Stripmap (SM), Interferometric Wide swath (IW) o Extra-Wide swath (EW), y por el nivel de
procesamiento multi-looking aplicado:
Resolución completa (FR).
Alta Resolución (HR).
Resolución media (MR).
En la siguiente Tabla se proporciona un resumen de las resoluciones espaciales que dependen del modo de
adquisición y del tipo de producto. Para el caso de productos Nivel-1 GRD se proporciona además una
subclasificación por clase de resolución (completa, alta y media). Nótese que, en este último caso, la
resolución completa (FR) solo está disponible en el modo de adquisición SM.
Todos los productos se procesan directamente desde el producto Nivel-0 SAR. Cada modo de adquisición
puede generarse en los niveles de procesamiento Nivel-1 SLC, Nivel-1 GRD y Nivel-2 OCN.
Las operaciones del segmento Tierra Sentinel-1 implementan un escenario de producción predefinido de
referencia. Este escenario prevé el procesamiento sistemático y la difusión en línea de todos los datos
Sentinel-1 adquiridos en modos IW, EW y SM en productos GRD de Nivel-0 y Nivel-1.
Además, se complementa con el procesamiento también sistemático y la difusión en línea de:
Para el modo WV, los productos de Nivel-0 SAR no se distribuyen.
Productos Nivel-1 SLC para todos los datos adquiridos del modo WV (activo desde octubre de 2016
para Sentinel-1B y desde mayo de 2017 para Sentinel-1A).
Productos Nivel-1 SLC para todos los datos adquiridos en modo IW en áreas regionales específicas.
Las áreas geográficas donde los productos Nivel-1 SLC se generan sistemáticamente han ido
creciendo desde la apertura del acceso a datos en línea el 03/10/2014.
Productos Nivel-2 OCN para todos los datos adquiridos del modo WV (activo desde julio de 2015
para Sentinel-1A y desde octubre de 2016 para Sentinel-1B).
Productos Nivel-2 OCN para todos los datos adquiridos en modo IW y EW en regiones específicas.
En este caso, para mayor claridad, se acompañan las descripciones con un ejemplo de un producto concreto:
S1A_IW_GRDH_1SDV_20180413T175451_20180413T175516_021452_024F10_6834.SAFE
El nombre de la carpeta del producto Sentinel-1 de nivel superior se compone de caracteres alfanuméricos en
MAYÚSCULAS separados por un guion bajo (_).
3.5 Aplicaciones
Sentinel-1 proporciona soporte a tres servicios prioritarios del programa Copernicus:
Copernicus Marine Environment Monitoring Service.
Copernicus Land Monitoring Service.
Copernicus Emergency Management Service.
Se puede encontrar más información sobre las áreas temáticas de Copernicus en la guía de áreas temáticas.
Cada aplicación puede atenderse mejor con modos de adquisición y tipos de productos de Sentinel-1
específicos, como se muestra en la Tabla 6. El modo de IW es el modo de operación principal para la
mayoría de las aplicaciones sobre Tierra.
Modo de adquisición
Aplicación
SM IW EW WV
Zonas árticas y hielo marino
Vigilancia de embarcaciones en mar abierto
Vigilancia de la contaminación por hidrocarburos
Vientos marinos
Ámbito forestal
Agricultura
Cartografiado de deformaciones urbanas
Cartografiado de inundaciones
Análisis de terremotos
Vigilancia de deslizamientos y volcanes
Tabla 6 Aplicaciones asignadas a distintos modos de adquisición
Nota: para varias islas volcánicas pequeñas en todo el mundo, desde mayo de 2015, se utiliza el modo SM para la cobertura regular.
Esto se implementa en el caso de que una swath en modo SM cubra por completo la isla en una pasada, y no se lleven a cabo
actividades de vigilancia marítima alrededor de la isla. El modo SM presenta swath más pequeño que el modo IW (80 km frente a
250 km), pero garantiza una resolución espacial más alta, de particular interés para la monitorización de volcanes.
Figura 35 Patrón clásico de la estela de una embarcación en imágenes SAR (S-1A Modo SM,
Copernicus Data [2015])
Figura 36 Derrame de petróleo observado en la costa belga cerca de Zeebrugge (Bélgica) el 10/08/2015
después de una colisión entre dos buques
Con objeto de determinar las fuentes y encausar a los infractores, las detecciones pueden además
correlacionarse con las transmisiones de información del Sistema de Identificación Automática (AIS,
Automatic Identification System) o de Identificación y Seguimiento a Largo Plazo (LRIT, Long-Range
Identification and Tracking) de los buques. Los datos relevantes para la supervisión de derrames de petróleo
son transmitidos por el satélite en tiempo real para su recepción por estaciones terrestres colaboradoras
locales que prestan su apoyo a servicios europeos y nacionales.
Otra aplicación está relacionada con la detección de filtraciones de petróleo de origen natural desde el fondo
del océano, que puede proporcionar pistas para futuras prospecciones de petróleo.
c) Vientos marinos
El SAR es sensible a los patrones de rugosidad de la superficie que varían espacialmente y que son causados
por los vientos en la superficie del océano (Figura 37). Los productos oceánicos Sentinel-1, junto con los
modelos de ola oceánica global, ayudan a determinar la dirección, longitud de onda y alturas de las olas en
océanos abiertos, proporcionando soporte en la predicción del clima, el transporte de embarcaciones y las
aplicaciones relacionadas con la energía de las olas.
Las estimaciones de los campos de viento también juegan un papel importante en la supervisión de derrames
de hidrocarburos, permitiendo discriminar causas accidentales de los vertidos voluntarios reales.
Además, Sentinel-1 proporciona datos de las interacciones entre las olas del océano y las corrientes,
permitiendo la visualización de corrientes oceánicas a gran escala, masas de agua fría / caliente, corrientes de
agua costeras y ondas internas.
a) Aplicaciones forestales
Sentinel-1 puede desempeñar un papel importante en la gestión forestal sostenible con la clear-cut and
partial-cut detection, clasificación del tipo de bosque, estimación de biomasa y detección de perturbaciones.
Para aplicaciones relacionadas con el cambio climático, el estudio y la delimitación de incendios forestales
tiene un elevado interés medioambiental, dado que estos pueden cambiar ciertos ecosistemas, y son una
importante fuente de aerosoles y gases de efecto invernadero. Su detección automática empleando series
MultiTemporales de Coherencia (MTC) calculadas a partir de imágenes radar Sentinel-1, permiten disponer
de un historial de la presencia de carbono de un bosque, jugando un papel fundamental en la estimación de
las emisiones de carbono.
En este ámbito de aplicación, se han obtenido buenos resultados en la delimitación de zonas quemadas
mediante el análisis de imagen radar MTC, lo cual abre la posibilidad de su uso en evaluación de daños
(Donezar Hoyos et al., 2017).
Finalmente, la cartografía de la cobertura terrestre se puede utilizar para apoyar la gestión forestal y la
supervisión de extracción ilegal de madera en todo el mundo.
b) Agricultura
La supervisión de las condiciones de los cultivos, las propiedades del suelo y el cartografiado de actividades
de labranza tillage, ayudan a evaluar el uso de la tierra, predicen las cosechas, supervisan los cambios
estacionales y apoyan la implementación de políticas para el desarrollo sostenible. Sentinel-1 también se
puede utilizar para cartografiar los cambios en la producción agrícola y la productividad de los pastizales
La clasificación de cultivos es una de las aplicaciones agrícolas más importantes de teledetección. Aunque la
teledetección radar se comenzó a desarrollar hace ya cuatro décadas, sus aplicaciones operacionales son aún
mucho menores que las basadas en teledetección óptica. Con la puesta en marcha de misiones como
Sentinel-1, se están desarrollando nuevas aplicaciones basadas en este tipo de observaciones, aprovechando
las ventajas de que su señal no se ve afectada por la nubosidad, además de proporcionar información sobre
características geométricas y dieléctricas de las cubiertas que puede ser útil para discriminar algunos cultivos
que no se podrían distinguir sólo con imágenes ópticas.
En cuanto a las características técnicas de los productos a utilizar, se han realizado diversos estudios con
diferentes bandas respecto a las longitudes de onda apropiadas a la hora de clasificar cultivos, y se ha
demostrado que la banda C es muy útil para aplicaciones agrícolas. Respecto a la polarización, se concluye
que la opción de polarización dual (VV-VH) puede ser una alternativa interesante para estas aplicaciones,
destacando además el hecho de que el breve periodo de revisita favorece la disponibilidad de series
multitemporales.
Figura 39 Subsidencias en Venecia a partir de interferometría del satélite ERS-2. Fuente: ESA
a) Cartografiado de inundaciones
Más del 75% de los desastres naturales que ocurren en todo el mundo están relacionados con inundaciones.
Las condiciones climáticas adversas durante los eventos de inundación hacen que el radar de apertura
sintética sea un enfoque adecuado debido a su independencia con respecto al clima y la iluminación solar.
Sentinel-1 proporciona una gran cantidad de datos con tiempos de revisita muy cortos, lo que mejora las
capacidades de las actividades de respuesta a inundaciones.
Entre las aplicaciones en este ámbito se encuentra la evaluación del alcance de áreas inundadas mediante el
cartografiado posterior al desastre (Figura 128). Además, los DEM de alta resolución generados a través de
los modos interferométricos de Sentinel-1 se pueden utilizar para realizar análisis de escorrentía y de
inundación en áreas que anteriormente carecían de ellos, alimentando los sistemas de alerta temprana.
Figura 41 Interferograma coseísmico SAR en el terremoto del año 2016 en Italia central procesado por
TRE ALTAMIRA. Cada ciclo de color (rojo-azul-amarillo) representa aproximadamente 3 cm de
movimiento en la línea de visión del radar
InSAR proporciona la capacidad única de producir mapas de deformaciones de terremotos de media y alta
resolución. A través de la supervisión constante de las áreas propensas a terremotos, se pueden descubrir las
fallas activas y estudiar los riesgos potenciales. El modo de adquisición EW facilita el seguimiento de
terremotos a gran escala.
Figura 42 Mediciones DInSAR sobre el condado de Maoxian para el período de octubre de 2014 a
junio de 2017 (arriba) y mediciones sobre el área de deslizamiento de tierra (abajo). La serie temporal
de desplazamiento superpuesta revela un claro movimiento precursor meses antes del evento
Obsérvese que la herramienta es común para los productos Sentinel-1, Sentinel-2, Sentinel-3, en este punto
nos centraremos en las opciones referentes a la misión Sentinel-2:
Sensing period: Rango de fechas para la búsqueda basada en la fecha de captura de los datos (paso
del satélite).
Ingestion period: Rango de fechas para la búsqueda basada en la fecha de introducción del producto
procesado en el sistema.
Satellite Platform: Pudiendo elegir entre las dos plataformas actualmente en órbita (S2A y S2B).
Product Type: Tipo de producto o nivel (S2MSI1C, S2MSI2A, S2MSI2Ap).
Relative Orbit Number: Identificador de orbita relativa.
Cloud Cover: Porcentaje de cobertura nubosa. Para expresar un rango de búsqueda lo haremos con
el siguiente formato [min To max], por ejemplo, utilizando la expresión [0 TO 20] se mostrarán
resultados con un máximo del 20 % de cobertura nubosa.
Situando el puntero sobre uno de los elementos del listado podemos acceder a distintas opciones:
Select product: Selección el producto para poder añadir más de uno a la vez a la “cesta” ( ).
Zoom to product: Centra la vista del mapa en el producto y hace zoom al mismo.
Add Product to chart: Añade el producto a la cesta para poder descargarlo más tarde, incluso en una
sesión distinta. Podemos acceder a la lista de productos que tenemos en la cesta mediante el botón de
usuario de la parte superior derecha, pulsando en el icono Cart:
También es posible descargar ficheros individuales, como los de metadatos o las imágenes de las bandas
separadas:
Además de utilizando el listado de resultados, también podemos seleccionar productos en el mapa, utilizando
el botón izquierdo o derecho del ratón según estemos en Navigation Mode o en Area Mode respectivamente.
Existen herramientas, disponibles bajo licencias de código abierto, que, haciendo uso de estas APIs, nos
facilitan la automatización de procesos de descarga si necesidad de programar herramientas propias o
reduciendo esa programación a la elaboración de sencillos scripts escritos en lenguajes como Python.
1.2.1 Sentinelsat
Sentinelsat es herramienta de código abierto (licencia GPL v3) que facilita la búsqueda, descarga y acceso a
los metadatos de imágenes de Sentinel-1, Sentinel-2 y Sentinel-3. Esta desarrollada en Python y ofrece dos
interfaces de uso: interfaz de línea de comandos (CLI) e interfaz de programación de aplicaciones (API).
Código fuente: https://github.com/sentinelsat/sentinelsat.
Documentación del proyecto: https://sentinelsat.readthedocs.io/.
Sentinelsat permite realizar búsquedas mediante diversos criterios, tales como el rango de fechas, el área
geográfica, el límite de cobertura nubosa, nombre del producto. Incluyendo el parámetro -d se inicia
automáticamente la descarga de los productos encontrados, en caso contrario sólo se muestran los resultados
de la búsqueda. Es necesario especificar el usuario y contraseña del sistema Copernicus Open Access Hub en
todas las llamadas al comando utilizando los parámetros -u y -p.
La geometría de las consultas se indica mediante un fichero que contenga un polígono en formato geojson.
Existe gran variedad de herramientas GIS, entre ellas QGIS, capaces de exportar una geometría en formato
geojson.
A continuación se muestra un fragmento de código con el cual se obtiene el número de imágenes Sentinel 2
en nivel 1C sobre el territorio de Guatemala, siendo necesario proporcionar:
Las claves de acceso al Sentinel data hub (usuario y contraseña).
Un archivo geojson con el área de interés.
Criterios de búsqueda (fecha, constelación, producto, etc.)
# conectar a la API
import geojson
from sentinelsat import SentinelAPI, read_geojson, geojson_to_wkt
from collections import OrderedDict
products = api.query(footprint,
date = ('20140101','NOW'),
platformname = 'Sentinel-2',
output_filename = 'imagenes_S2L1C_guatemala.geojson'
with open(output_filename,'w') as output_file:
geojson.dump(datasets,output_file)
total_size_images = 0
n_images = 0
En tipos de ficheros vemos All Files, se podría filtrar por Sentinel-2 MSI productor tile (*.xml).
En el panel de explorador de productos aparecerá una nueva entrada, desplegaremos el árbol para poder ver
como se estructura la información del producto y aquella que, aunque en principio puede no existir, se puede
ir generando en los distintos procesos que se lleven a cabo en el propio entorno de la aplicación. Está
agrupada en cuatro contenedores: Metadata, Vector Data, Bands y Masks, donde se accede a la información
relativa a los metadatos del producto, información vectorial, información de la imagen y las máscaras,
respectivamente.
Si se despliega Bands se accede a las bandas de la imagen. Haciendo doble clic en la banda que se desee
visualizar o sobre la banda en el explorador presionar botón derecho abriremos la imagen en el visor central.
En la parte derecha, existe una pestaña donde se mostrará el panel de propiedades de la banda seleccionada:
nombre, descripción, resolución espacial, espectral, etc. y en la parte inferior la información resumida.
También se puede acceder se puede abrir este panel en botón derecho -> Properties.
La información asociada a cada píxel, localizado en el cursor en el visor en el que se está navegando, se
visualiza en la pestaña Pixel Info junto al explorador de productos. Muestra la posición del cursor en
coordenadas imagen, coordenadas geográficas y coordenadas proyectadas; el nivel digital en este producto es
el valor de reflectancia BOA de cada píxel e información temporal.
Para visualizar la imagen en RGB, en el explorador de productos se debe pulsar el botón derecho sobre el
producto y elegir la opción correspondiente. A continuación, se podrá seleccionar composición en color
natural, composición en falso color infrarrojo y penetración atmosférica que como combinaciones definidas
por defecto, pero se pueden modificar por el usuario y salvar como nuevos perfiles. Además, se pueden
almacenar como imágenes virtuales en el propio producto.
Figura 20 Selección de bandas para una visualización RGB (izquierda) y falso color (derecha)
De este modo, se pueden generar tantas composiciones como requieran los análisis que se estén realizando
Para manipular la forma de visualización podemos acceder al panel de Colour Manipulation situado en la
parte inferior izquierda de la aplicación. Primero procederemos a calcular las estadísticas precisas de la
imagen.
Se crea una banda nueva done se puede visualizar el DEM descargado en este ejemplo SRTM 3Sec.
Aunque es opcional, se recomienda la actualización de plugins, para ello se accede a Tools -> Plugins. Por
defecto tras la instalación completa accedemos a la pestaña Installed.
Si hubiera pasado tiempo desde la instalación, la aplicación detecta si existen actualizaciones disponibles,
aunque siempre podemos consultarlas en la pestaña Updates, y proceder a su actualización en caso de ser
necesario. Si el programa solicita reiniciar el sistema para completar el proceso de actualización, debemos
hacerlo.
A continuación en este apartado se verá la instalación de Sen2Cor integrado en SNAP. En primer lugar, se
instalará el plugin que controlará Sen2Cor. En la pestaña Available Plugins podemos consultar los plugins
que no tenemos instalados. Se selecciona la última versión disponible del procesador Sen2Cor. En el
momento de escribir esta guía el plpugin era el: Sentinel-2 SEN2COR280 Processor, tal y como se muestra
en la Figura 25.
Tras realizar esta tarea, SNAP nos pedirá reiniciar. Tras hacerlo, si accedemos a la pestaña de plugins
instalados, el plugin se encontrará activo en la pestaña Installed. Si no hemos reiniciado la aplicación, se
debe selccionar y pulsar en activar.
El siguiente paso será instalar Sen2Cor, porque en la fase previa se ha instalado el plugin pero no la
aplicación en sí misma. Para manejar las utilidades externas se accede al Menú Tools->ManageExternal
Tools.
Si se accede de nuevo a las variables del sistema, se puede observar que SEN2COR_BIN apunta a donde se
ha instalado la aplicación, la cual podemos comprobar.
Se puede consultar el fichero L2A_Process y comprobar las variables de sistema y la llamada al programa en
Python, que realizará el proceso y la ubicación del intérprete de Python que integra la aplicación Sen2Cor.
La selección de LUT se puede configurar a través del archivo de la configuración del usuario
(L2A_GIPP.xml) ubicado dentro de la carpeta cfg del directorio que se encuentra en la carpeta donde apunta
la variable de entorno $SEN2COR_HOME. En la selección Look_Up_Tables del archivo de configuración.
Sen2Cor puede configurarse por el usuario, teniendo en cuenta los siguientes criterios:
Tipo de aerosoles: por defecto se establece RURAL. Si el usuario establece el tipo de aerosol en
AUTO, se determinará automáticamente por la aplicación, aplicando un algoritmo descrito
anteriormente.
Tipo de atmósfera: por defecto la aplicación establece SUMMER para latitudes medias, si el
usuario indica AUTO se determinará en función de la fecha de adquisición de la imagen.
El valor de la concentración de ozono por defecto se establece en 331 unidades Dobson (DU). Si el
usuario establece Ozone_Content ‘0’, el valor se determinará automáticamente por la aplicación,
accediendo a la información almacenada en el fichero denominado AUX_ECMWFT, ubicado en la
carpeta AUX_DATA de la tile del producto de entrada Nivel-1C, y se utiliza la LUT que mejor se
ajuste a la concentración de ozono medida. Puede variar a escala global, valores extremos pueden
presentarse (200-500 DU), no obstante, el rango típico de valores se encuentra en 250 - 400 DU.
El espesor óptico de aerosoles en la atmósfera, el contenido de vapor de agua atmosférico y el
tratamiento de los cirros es calculado internamente a partir de la propia imagen.
El valor por defecto empleado para la visibilidad por Sen2Cor son 40km, lo cual es una buena
estimación para el interior de la Península Ibérica. Debería ajustarse según cada caso.
La ejecución de la aplicación Sen2Cor en el interfaz de SNAP puede ser configurada editando la
herramienta. Accediendo a la pestaña Operator Parameters dentro del menú External tolos se pueden
establecer los valores por defecto que pueden ser configurados por el usuario, tal y como se ha comentado
previamente. También se puede modificar la descripción del parámetro.
Desde el menú superior de SNAP puede accederse a la opción Optical>ThematicLand Processing >
Sen2Cor para realizar la corrección atmosférica. También puede ejecutarse Sen2Cor en el Menú Tools-
>Manage External Tools.
Del mismo modo, haciendo clic sobre el lienzo con el botón derecho, se podrá conectar la entrada, el
operador y la salida de forma automática, aunque también se pueden conectar de forma manual moviendo las
flechas que surgen cuando se flota sobre la caja.
Si se compara esta llamada con la que se hacía previamente, la única diferencia es que en este caso se
generará un fichero plantilla y se establecen por defecto 2 ficheros de parámetros (GIPP).
Se selecciona la imagen de entrada y podemos guardar el modelo creado. Se lanza el proceso pulsando en el
botón Run y cuando finaliza lanza un mensaje informando del tiempo de procesamiento.
En la carpeta donde se encuentra el producto de entrada Nivel-1C se habrá generado el producto Nivel-2A.
El modelo sen2cor. xml guarda la configuración del proceso para ser reutilizada. La guardaremos en nuestro
sistema para emplearlo en el siguiente apartado del tutorial.
A continuación, se muestra paso a paso, ventana a ventana, la generación de un .bat para procesar con
Sen2Cor.
Se accede en el menú Tools -> Batch Processing.
A continuación cargamos el flujo de proceso que guardamos en el paso anterior en File > Load Graph, y
ponemos en la tabla las imágenes que vamos a querer procesar mediante el flujo definido. Con esto
conseguiremos que se nos realice la corrección atmosférica a 3 imágenes Sentinel 2 nivel 1C sobre la zona
que queremos.
Una vez abiertas anteriormente las imágenes, en Raster -> Geometric Operations -> Resampling podemos
configurar los parámetros que queremos obtener. En nuestro caso, en la pestaña de Resampling Parameters
elegimos igualar la resolución por píxel a 10 metros (podría hacerse por banda de referencia también) y el
método bilinear. Dejamos el nombre de salida por defecto y clicamos Run.
Una vez tenemos las 2 imágenes resampleadas se procede a realizar el subset con la finalidad de reducir el
número de bandas de 12 a 4 (bandas 2, 3, 4 y 8).
En Raster -> Subset nos aparece un menú con 3 pestañas. En la pestaña Band Subset nos aparecen la
totalidad de bandas. Solo dejaremos activas las bandas 2, 3, 4 y 8 y desactivamos el resto. Realizaremos este
paso con las dos imágenes.
3.2.3 Mosaicado
El mosaico de cada par de fotos es el proceso final del pre-procesamiento. Teniendo ya las imágenes a 10
metros de resolución (resampling) y con sólo 4 bandas (Subset), se procede a obtener el mosaico de cada par
de fotos contiguas de la misma fecha.
En la pestaña Variables & conditions, simplemente se añaden las 4 bandas que van a intervenir en el proceso
con el primer botón que aparece en el menú.
El procesado del mosaico, con una máquina con procesador Intel-core i7-4790 y 8 Gb de RAM demoró un
total de 18 minutos.
La imagen mosaicada final tiene 4 bandas, con una resolución espacial de 10 metros, que nos permite hacer
visualizaciones en RGB y en falso color natural con la banda del infrarrojo cercano.
Figura 42 Visualización en falso color natural y en color natural del mosaico realizado
Figura 43 Zonas afectadas por un incendio en el departamento del Petén. Imagen pre-incendio
(arriba) y post-incendio (abajo)
Podemos incluso realizar una mejor visualización mediante otra combinación de bandas para discriminar
mejor las áreas quemadas de las que no. Vamos a cambiar la visualización de la clásica de color real anterior
(4-3-2) a otra diferente (12-11-8A). Para ello, cambiamos las bandas en la función Open RGB Image
Aplicamos esta combinación de bandas a las dos imágenes y obtenemos una visualización como la siguiente:
Figura 45 Comparación entre visualización con combinación 4-3-2 y 12-11-8ª. Pre-incendio abajo y
post-incendio arriba
Con esta visualización ya se discriminan de manera visual las zonas quemadas de las que no, y qué zonas
han ardido con mayor severidad que otras. Para emplear más herramientas de SNAP y automatizar las tareas,
vamos a implementar un flujo de trabajo que consistirá en lo siguiente:
La manera de implementar este flujo de trabajo es mediante el graph builder tal y como hemos hecho en la
práctica anterior en Sen2Cor.
Realizaremos esto en graph builder y lo guardamos como un archivo xml en el directorio que estimemos
conveniente. A continuación emplearemos el módulo batch processing para poder realizar esta serie de
operaciones secuencialmente a las dos imágenes descargadas. El primer paso es por lo tanto cargar las dós
imágenes en la pestaña I/O parameters.
A continuación tenemos que escoger las bandas que queremos emplear y especificar si queremos realizar el
procesado en toda la imagen o únicamente en un subconjunto de ella. Haremos lo siguiente en la pestaña
subset:
Elegiremos de todas las bandas, la 3, la 8 y la 12, que emplearemos posteriormente para calcular el
índice NBR y otro producto adicional.
En dimensiones, únicamente vamos a procesar una pequeña parte de la imagen para facilitar el
cálculo a las máquinas. Estudiaremos el incendio presente en la zona definida por las coordenadas
pixel (4.500,6.000) y (7.000,8.500).
Esta expresión la vamos a implementar dentro de la pestaña BandMath como se indica en la Figura 51.
Lo que se va a realizar a continuación es el fusionado de los dos productos pre-procesados para calcular
rápidamente los valores de NBR pre y post incendio.
Para ello vamos a acceder a Raster -> Geometric operations -> Collocation
Los cuerpos de agua pueden mostrar comportamientos análogos a zonas quemadas, por lo que es necesario
enmascararlas. En este apartado vamos a aprender a enmascarar zonas acuáticas. Para ello emplearemos el
Índice de Diferencia de Agua Normalizado o NDWI, el cual se calcula con la siguiente expresión:
El índice NDWI está diseñado para aprovecharse de la máxima reflectancia de los cuerpos de agua en la
banda del verde y la mínima reflectancia en el infrarrojo cercano (NIR), tal y como se muestra en la Figura
55.
Para realizar el cálculo, hacemos clic derecho sobre la imagen generada y seleccionar BandMaths.
Estableceremos los siguientews valores,
Nombre: mascara_agua
En edit expression introduciremos lo siguiente:
if ($5.B3_pre- $5.B8_pre)/($5.B3_pre+ $5.B8_pre) >= 0.0 then 1 else 0
De esta manera obtenemos una nueva imagen con como se muestra en la que 1 corresponde a pixeles de agua
y 0 a píxeles no de agua.
Para identificar las zonas recientemente quemadas y diferenciarlas del suelo desnudo u otras zonas no
cubiertas por la vegetación, se empleará la diferencia entre el valor del NBR calculado antes y después del
incendio:
Este valor es una diferencia absoluta que puede presentar problemas en áreas con bajas cubiertas vegetales
antes del incendio. En estas zonas, el cambio absoluto entre el NBR antes y después del incendio será
bastante pequeño. En estos casos, conviene emplear una versión relativizada de este nuevo índice:
De esta manera se obtiene una imagen nueva en la que los píxeles quemados aparecen mucho más claros ya
que la diferencia entre los valores de los píxeles anteriores y posteriores al incendio son mucho mayores.
Ahora vamos a crear una nueva banda en la que únicamente se incluyan áreas quemadas. Para ello vamos a
establecer un umbral de clasificación para área quemada de 0.5. Cualquier pixel con valor mayor de 0.5 de
NBR será clasificado como área afectada por el incendio.
Volvemos de nuevo a emplear BandMaths y aplicamos la siguiente expresión:
if $5.RBR > 0.5 then $5.RBR else NaN
Obteniéndose el siguiente resultado:
Ahora vamos a exportar esta imagen para visualizarla en otros softwares como puede ser QGIS o Google
Earth.
Para ello vamos a File -> Export -> GeoTIFF
Escogemos el directorio donde queremos guardar la imagen de salida y en la derecha del menú de diálogo
que aparece vemos un botón que indica subset. Haciendo clic en el botón llegamos a otro menú donde
podemos escoger las dimensiones espaciales de la imagen a exportar y un subconjunto de bandas. En nuestro
caso, elegiremos las bandas RBR y RBR_ para exportar, dejando el resto deseleccionadas.
Aceptamos y generamos la imagen en formato GeoTIFF.
Abrimos QGIS y cargamos la imagen. Para una mejor visualización de la severidad de las zonas afectadas,
podemos emplear la clasificación que hace el United States Geological Survey (USGS) para interpretar este
NBR.
Añadiendo un mapa base como el de OpenStreetMaps a través del plugin OpenLayersPlugin podemos
obtener visualizaciones para la generación de cartografía.
Figura 2 Ventanas atmosféricas y su relación con el máximo de energía emitida por el Sol
Debido a estas propiedades, la radiación microondas puede emplearse para caracterizar cubiertas bajo
prácticamente todas las condiciones medioambientales imaginables.
Las diferentes frecuencias en las que se trabaja en radar tienen una capacidad mayor o menor de penetrar en
las cubiertas, en función de su longitud de onda. Las principales aplicaciones según las bandas se resumen en
la siguiente tabla.
Ka Ku X C S L P
Uso general, más usado
Penetración de cubierta vegetal
Captura subsuelo
Estimación de biomasa
Monitoreo de agricultura
Monitoreo de océanos y hielos
Estudios de subsidencias
Monitoreo de nieve
Muy alta resolución
Figura 3 Fundamentos de emisión y receptor de un pulso radar. Fuente: John Wiley & sons, Inc.
(1994)
En este caso es importante señalar la diferencia entre el fundamento del sensor LiDAR y el del radar. En el
caso del LiDAR, todos los retornos de un pulso vienen de la misma dirección, mientras que en el caso del
radar, todos provienen del mismo acimut, pero en distintos alcances.
Las medidas radar pueden describirse mediante una ecuación básica que relaciona la potencia, la distancia, la
reflectividad y las características de la antena:
Donde:
es la potencia del pulso recibida.
es la potencia del pulso emitida.
es la distancia entre la antena emisora y el objeto en el terreno.
es la distancia entre el objeto en el terreno y la antena emisora.
es la ganancia de la antena.
es el área de la antena.
es la retro-dispersión del objeto en el terreno.
En un sistema radar de iluminación lateral se produce una imagen de alta resolución de la superficie bajo
observación. Al moverse a lo largo de su trayectoria, el radar ilumina hacia un lado de la dirección de vuelo
franjas continuas y paralelas entre sí de la superficie en estudio y acumula la información de las microondas
que se reflejan. La señal que se graba a bordo, se procesa apropiadamente para formar una imagen digital. La
distancia entre el radar y el objetivo en la superficie en la dirección perpendicular al vuelo se llama alcance,
mientras que se conoce como acimut la distancia a lo largo de la trayectoria. Por lo tanto, en un sistema
radar la resolución tiene dos dimensiones, una en la dirección del alcance y otra en la dirección del acimut.
Mediante procesamiento digital de la señal, la imagen puede enfocarse y obtener así una resolución mejor
que la de un radar convencional.
Donde L es la longitud de la antena y la longitud de onda en la que opera el radar. Esta expresión puede
transformarse en resolución en acimut resultando:
Siendo R la distancia de la antena al objeto a observar. Tomando ejemplos numéricos: un radar de 10 metros
de longitud de antena, operando en banda C (longitud de onda 5 centímetros) y observando objetos a 1000
kms de distancia, proporciona una resolución en acimut de 5 kilómetros. Estas características serían las
comunes en sensores radar embarcados en satélites, y puesto que antenas de centenares de metros no son
tecnológicamente viables, es necesario emplear otra técnica.
Los sistemas radar de apertura sintética o SAR son más sofisticados que los de apertura real. Se puede
aumentar considerablemente la resolución en acimut simulando una antena de dimensiones mayores a la que
realmente se emplea. Un array de antenas es equivalente a una antena que se mueve a lo largo de una línea si
las señales se miden de manera coherente (amplitud y fase como una función del tiempo) y se añaden las
unas a las otras, asumiendo que el objeto está estático durante el momento de la captura.
Mediante procesamiento digital de la señal, la imagen puede enfocarse y obtener una resolución mejor que la
de un radar convencional.
Los efectos geométricos que se producen en una imagen radar y el relieve son los siguientes:
Figura 8 Ejemplos de escorzo (izquierda) y de inversión por relieve (derecha) en una imagen SAR. La
toma de la imagen está realizada desde arriba hacia abajo. Fuente: ESA
Figura 9 Ejemplo de diferente huella del satélite Sentinel 1A a su paso por Madrid
(ascendente/descendente)
La dispersión a una determinada longitud de onda variará por una serie de condiciones, como el tamaño
físico de los dispersores, las propiedades eléctricas de estos y su contenido en humedad. La propia longitud
de onda de los pulsos SAR y la polarización de estos (de este asunto trataremos más adelante), también
afectarán a la retro-dispersión de la señal.
Se producen fenómenos de dispersión por superficie cuando la radiación incide sobre una superficie
determinada. Este fenómeno depende de la rugosidad o textura de dicha superficie y de su constante
dieléctrica, resultando:
A mayor constante dieléctrica de una superficie o material, mayor será la retro-dispersión y por lo
tanto más brillante se observará un pixel en la imagen radar. Los objetos metálicos y el agua tienen
una constante dieléctrica alta (80) mientras que el resto de materiales tienen una más baja. En
condiciones secas la constante también es baja (3-8), lo que implica que la humedad en suelos o en
cubiertas vegetales puede producir un aumento significativo de reflectividad de la señal radar (25-
30). Por lo tanto, este tipo de teledetección puede usarse para obtener el contenido en humedad de
los suelos.
Además se producirán fenómenos de dispersión por volumen principalmente en cubiertas vegetales.
Estas propiedades definirán el valor del coeficiente de retro-dispersión ( ) ya visto en apartados anteriores.
Valores típicos de este coeficiente en diferentes superficies se muestran en la siguiente tabla.
Valores de retro-dispersión Superficies típicas
Objetos antropogénicos (ambientes urbanos).
Muy alto (más de -5 dB) Pendientes orientadas hacia el sensor.
Superficies muy rugosas.
Alto (de -10 dB a 0 dB) Superficies rugosas.
Vegetación densa (bosques).
Niveles medios de vegetación.
Moderado (de -20 a -10 dB) Cultivos.
Superficies moderadamente rugosas.
Superficies lisas.
Bajo (menos de -20 dB) Aguas en calma.
Carreteras.
Suelos muy secos (arenas).
Tabla 3 Valores típicos del coeficiente de retro-dispersión
El resultado de este fenómeno es una imagen con el denominado efecto de “sal y pimienta”, con una textura
de imagen que se asemeja precisamente a este fenómeno. Es un fenómeno que ocurre únicamente con
sensores coherentes como el SAR, en los que se producen interferencias de la señal en dimensiones
inferiores a la longitud de onda.
Es por esta razón que en segmentación y clasificación de imágenes SAR se incluyan criterios de textura
además de la intensidad de las señales.
Para reducir los efectos de speckle en una imagen se usan diferentes criterios:
Multi-looking: divide la señal en frecuencias mínimamente solapadas, se procesan estas frecuencias
a una resolución imagen reducida y se añaden finalmente para obtener una imagen final.
Filtros espaciales: detecta la resolución total de la imagen, efectúa medias locales sobre los píxeles
y finalmente realiza un remuestreo para crear resoluciones más pequeñas pero con el efecto speckle
reducido.
1.7 Polarización
Independientemente de la longitud de onda, las señales radar pueden transmitirse como campos
electromagnéticos vectoriales horizontales (H) o verticales (V), y recibirlos de vuelta como señales
horizontales (H), verticales (V) o ambas. Por lo tanto, existen las siguientes posibilidades en la emisión-
recepción de la radiación polarizada:
HH: transmitida horizontalmente y recibida horizontalmente.
HV: transmitida horizontalmente y recibida verticalmente.
VH: transmitida verticalmente y recibida horizontalmente.
VV: transmitida verticalmente y recibida verticalmente.
La diferencia de fase entre dispersores de la onda incidente que viaja desde el radar hasta el dispersor y
vuelve cambia en función de:
Donde es la diferencia en la distancia de viaje del pulso. Dado que el pixel de una imagen SAR contiene
un número muy grande de dispersores, la fase de cada uno de los píxeles se distribuye aparentemente de
manera aleatoria.
Si ya estamos dados de alta en el sistema, podremos realizar búsquedas de los distintos productos. En este
punto del curso nos centraremos en las opciones referentes a la misión Sentinel-1:
Sensing period: Rango de fechas para la búsqueda basada en la fecha de captura de los datos (paso
del satélite).
Ingestion period: Rango de fechas para la búsqueda basada en la fecha de introducción del producto
procesado en el sistema.
Satellite Platform: Pudiendo elegir entre las dos plataformas actualmente en órbita (S1A y S1B).
Product Type: Tipo de producto (SLC, GRD y OCN).
Polarisation: Polarización de emisión y recepción de la señal (HH, VV, HV, VH, HH+HV,
VV+VH).
Sensor mode: Se puede escoger entre los diferentes métodos de captura de imágenes de Sentinel 1
(SM, IW, EW, WV).
Relative Orbit Number: Identificador de orbita relativa.
El área geográfica de búsqueda se especificará dibujando un polígono en el mapa. Para ello podemos hacer
clic y arrastras para dibujar un rectángulo o hacer clic sucesivamente en distintos puntos para definir otro
tipo de polígono.
Situando el puntero sobre alguno de los elementos del listado podemos acceder a distintas opciones:
Figura 3 Información relativa a uno de los productos listados. Se muestra por ejemplo un producto
SLC, un GRD y un OCN
1.2.1 Sentinelsat
A continuación se muestra un fragmento de código con el cual se obtiene el número de imágenes Sentinel 1
con producto sobre el territorio de Guatemala, siendo necesario proporcionar:
Las claves de acceso al Sentinel data hub (usuario y contraseña).
products = api.query(footprint,
date = ('20140101','NOW'),
platformname = 'Sentinel-1',
producttype = 'GRD')
total_size_images = 0
n_images = 0
En el caso de productos SLC, bastaría con sustituir ‘GRD’ por ‘SLC’ en producttype, obteniendo los
siguientes resultados:
Número total de imágenes Sentinel 1 SLC sobre Guatemala hasta el día en el que se ejecuta el script
(30 de octubre): 2667 imágenes.
Tamaño total que ocupan dichas imágenes en disco (comprimidas): 3,57 terabytes.
Fichero geojson con las huellas y metadatos de dichas imágenes, el cual se muestra en la Figura 5.
Figura 5 Visualización del fichero geojson y tabla de metadatos asociada en QGIS para S1 SLC
Figura 8 Tipos de visualizaciones posibles por defecto en SNAP para imágenes SAR
Una visualización de las 3 imágenes nos arroja los siguientes resultados:
La imagen calibrada radiométricamente se habrá creado con el nombre indicado y en el directorio de salida
elegido. Se abre en el explorador de productos y se puede ver que en la carpeta de metadata se ha generado
Los filtros speckle se encuentran en la aplicación SNAP en el menú Radar-> Speckle Filtering. En esta
primera fase al tratarse una única imagen se procesará con Single Product Speckle Filter.
Existen distintos filtros, el más común es Filtro Lee Sigma. En este caso se opta por utilizar el filtro Refined
Lee, que aporta algunas mejoras respecto al tradicional.
Figura 15 Diferencia entre imágenes sin filtro (arriba) y con filtro speckle (abajo)
3.4.2 Multilooking
Se trata de otra manera de terminar con el efecto de sal y pimienta. En general o empleamos un filtro de
speckle o empleamos un criterio de multilloking. En este apartado explicaremos cómo realizar el
multilloking.
A continuación, se procederá a crear como bandas físicas, para ello nuevamente botón derecho y hacer click
en Convert Band.
Lo que vamos a hacer es guardar el gráfico de procesamiento sin particularizar entrada y salida para
posteriormente realizar un Batch Processing. Para ello, siempre que queramos hacer uso de este grafo de
procesado, iremos al menú Tools -> Batch Processing y se procede a abrir el gráfico xml guardado en el paso
anterior, en File -> Load Graph o en el botón correspondiente, de la ventana del Batch Processing.
Una vez abierta la imagen en cuestión, se accede al histograma de los valores de retrodispersión en ambas
polarizaciones. En el panel de Color Manipulation se modifica la visualización. Las masas de agua y las
áreas inundadas se pueden apreciar que son áreas de poca retrodispersión, zonas oscuras, corresponden a
valores aproximadamente inferiores a -20.3 dB. Como la proporción entre zonas con agua es significativa en
la imagen aparecen claramente dos picos en el histograma. Las zonas con tierra son más claras.
Figura 29 Expresión matemática para el álgebra de bandas con el objetivo de discriminar agua
El resultado es una imagen donde el valor 255 representaría presencia de agua.
Figura 30 Imagen clasificada en pixeles con agua (blanco, 255) y no agua (negro, 0)
A continuación, se procede a generar un tipo de visualización muy común en radar que permitirá
fotointerpretar las consecuencias de la inundación. Para ello, se empleará una imagen resultado de calcular
las diferencias entre las dos polarizaciones (VV-VH). Esta operación se puede hacer virtualmente para que se
pueda visualizar o físicamente creando la banda. En cualquier caso, posteriormente se podría guardar y/o
exportar. Se procede a generar una visualización en RGB empleando las dos polarizaciones y la diferencia en
los canales RGB. Para ello, se pulsa con el botón derecho Open RGB Image Window, se elige un perfil de
visualización de partida, tal y como muestra la siguiente Figura:
En la siguiente pestaña se seleccionará Product Geolocation en Initial Offset Method. Adviértase que, tal y
como ilustra la siguiente Figura, la primera imagen actuará como imagen maestra.
El resultado será una imagen compuesta de 8 bandas con las polarizaciones de ambas fechas:
Aunque es opcional, se puede guardar el perfil de selección de canales para su uso posterior. También se
pueden guardar como bandas virtuales el modo de visualización RGB.
En la imagen se visualiza en color negro las zonas de agua permanente, como lagunas y el cauce del río; en
color azul oscuro la zona inundada; en color fucsia se asocia a vegetación con agua; y en color blanco se
observan las zonas urbanas y de construcciones artificiales.
Haciendo clic con el botón derecho sobre el visor de la imagen se presentan distintas opciones, una de ellas
es la opción Export View as Image. Para exportar toda la imagen a la máxima resolución se deben elegir las
opciones Full scene en el grupo de opciones Image Region y Full resolution en el grupo de opciones Image
Resolution.
Es posible exportar la imagen en RGB a otros formatos ráster. En este caso de uso se va a exportar a formato
GeoTIFF. Si se desea guardar todo el producto, no sólo la visualización RGB, también se podría exportar en
File -> Export -> GeoTIFF/BigTIFF.
Este ejemplo se realizará con la misma pila de imágenes empleada en el proceso anterior.
Para realizar una clasificación no supervisada se accede al menú Raster -> Classification -> Unsupervised
Classification -> EM Cluster Analysis.
Símbolo Clase B G B
Agua no permanente 0 255 255
Agua permanente 0 0 255
Cultivos herbáceos 0 255 0
Cultivos ribereños 124 242 142
Asfalto 128 128 128
Zona urbana 255 0 0
Roquedo 80 41 22
Suelo desnudo 219 191 48
Forestal 27 132 34
Matorral 83 169 93
Se irán creando en el Vector Data del producto activo, los nuevos ítems con la nomenclatura establecida por
el usuario.
La aplicación interrogará (aunque a veces lo obvia) acerca del almacén vectorial donde se desee incorporar el
recinto. De esta forma, se van añadiendo las geometrías y pueden ser consultadas haciendo doble clic en el
ítem correspondiente incluido en el Vector Data, desplegándose una tabla donde se visualiza la información
geométrica en WKT y la simbología.
Finalmente, se guardará como otro producto con la opción Save Product o Save Product As.
Entramos ahora en la fase de entrenamiento y asignación de clases. SNAP incorpora diversos algoritmos
para llevar a cabo el entrenamiento del clasificador y para la asignación posterior de clases. En este caso se
utilizará el algoritmo Random Forest.
Antes de realizar la clasificación es necesario renombrar las bandas de la imagen, ya que SNAP puede
arrojar un error por no discriminar diferentes bandas dentro de la imagen. Por lo tanto, es necesario cambiar
los nombres de las bandas haciendo clic derecho en cada una de ellas y luego pulsando en properties. Se
abrirá una ventana y cambiaremos lo que está escrito en el apartado Name. Por ejemplo, para la primera
banda haremos lo siguiente:
Cambiaremos el nombre de Sigma_0_VH_mst_25Sep2018 a S0_VH_25Sep2018.
Usaremos imágenes de agosto del año 2019 para realizarlo. A la hora de seleccionar las imágenes se
pretenderá hacerlo siguiendo el criterio de mantener aquellas que se capturasen en la misma pasada del
satélite. De esta manera se seleccionaron las siguientes imágenes:
S1A_IW_GRDH_1SDV_20190823T120201_20190823T120226_028696_033F8B_27F1
S1A_IW_GRDH_1SDV_20190823T120226_20190823T120252_028696_033F8B_2070
S1B_IW_GRDH_1SDV_20190824T115254_20190824T115319_017727_0215A7_1D47
S1B_IW_GRDH_1SDV_20190824T115319_20190824T115344_017727_0215A7_7E44
S1B_IW_GRDH_1SDV_20190824T115344_20190824T115409_017727_0215A7_4CCD
S1A_IW_GRDH_1SDV_20190825T114540_20190825T114605_028725_03409B_9858
S1A_IW_GRDH_1SDV_20190825T114605_20190825T114630_028725_03409B_CBA9
En la siguiente pestaña podrían seleccionarse las bandas a incorporar, en este caso incluimos todas. En el
tamaño de pixel de salida estableceremos 70 metros, puesto que es la resolución de salida que nos da la
operación de speckle con filtrado de tipo Refined Lee.
Finalmente, se seleccionará la carpeta de salida, el nombre del mosaico y ejecutamos el proceso con run. El
proceso no lleva bastante tiempo. En el caso de las 7 imágenes sobre Guatemala demoró cerca de 1 minuto
con una máquina con procesador inte-core i7-4790 y 8 Gb de RAM.
El mosaico ocupa una extensión como la que se muestra en la Figura 61, donde puede apreciar el espacio
vacío sin imágenes que la imagen tiene que rellenar para poder crearse. Esto causa el incremento importante
en cuanto a tamaño del archivo final. En la Figura 62 puede verse el mosaico en niveles de gris con el
perímetro del país de Guatemala superpuesto.
Figura 62 Visualización del mosaico de imágenes Sentinel 1 sobre Guatemala con la silueta del país
superpuesta
Tras aplicar dicha función se obtiene la siguiente Figura, donde ya tenemos la silueta de Guatemala
recortada. No obstante repárese en que existen muchos píxeles vacíos susceptibles de recortar.
Sus posibles aplicaciones son muy diversas como la planificación espacial, gestión forestal, gestión del agua,
agricultura y seguridad alimentaria, gestión de emergencias.
La implementación conjunta del servicio, en funcionamiento desde 2012, corre por cuenta de la Agencia
Europea de Medio Ambiente y del Centro Común de Investigación de la Comisión Europea.
El servicio del territorio de Copernicus está estructurado en cuatro componentes: Global, Pan-Europeo, Local
y datos de referencia los cuales se describen en los siguientes apartados.
Global: El seguimiento sistemático de los parámetros biofísicos genera principalmente una serie de productos
biogeofísicos cualificados sobre el estado y la evolución de la superficie terrestre. Estos datos a escala mundial
se facilitan cada diez días con una resolución espacial intermedia y se complementan con una serie temporal a
largo plazo. Los productos se utilizan para monitorizar la vegetación, los cultivos, el ciclo del agua, el
presupuesto energético y las variables de la criosfera terrestre.
Paneuropeo: La cartografía de la ocupación y el uso del suelo produce clasificaciones de ocupación del suelo
con distintos niveles de detalle, en un contexto paneuropeo o mundial. A nivel paneuropeo, esta información
se complementa con capas detalladas sobre las características de la cobertura del suelo (como impermeabilidad,
bosques, prados, agua y humedad o pequeños elementos leñosos). A nivel mundial, la cartografía de la
ocupación del suelo se ajusta al sistema modular y jerárquico de clasificación de la cubierta terrestre de la
FAO.
Local: La cartografía temática de "puntos calientes" tiene por objeto proporcionar información específica más
detallada sobre áreas de interés concretas, conocidas como puntos calientes. En el contexto del Servicio de
Vigilancia Terrestre de Copernicus, los puntos calientes suelen corresponder a desafíos medioambientales
específicos.
Las imágenes y datos de referencia: proporcionan mosaicos de imágenes satelitales de alta y muy alta
resolución, así como conjuntos de datos de referencia. Incluyen, por una parte, los mosaicos de imágenes
satelitales procedentes de las misiones participantes que cubren el territorio de Europa, así como la producción
de mosaicos de imágenes Sentinel-2 a nivel mundial. A esto se añaden conjuntos de datos de referencia que
proporcionan una cobertura paneuropea homogénea de algunos aspectos geoespaciales clave, como la
hidrografía y la elevación.
Además de los componentes ya mencionados, se está poniendo en marcha una nueva actividad europea de
movimientos del suelo, que medirá los desplazamientos del terreno (corrimientos y hundimientos de tierras) y
la deformación de las infraestructuras.
De estos cuatro componentes el IGN/CNIG, en lo que se refiere a ocupación del suelo, analiza, verifica y
produce información de los componentes Pan-Europeo y Local.
Los ecosistemas de la Tierra cambian constantemente debido a la naturaleza y las condiciones atmosféricas
bajo la presión de las actividades humanas. Para monitorear los cambios en los biomas continentales, el
servicio Copernicus Global Land proporciona, de manera oportuna, un conjunto de variables biofísicas que
describen el estado, el dinamismo y las perturbaciones de la vegetación terrestre.
El índice de vegetación de diferencia normalizada (NVDI) ampliamente utilizado, y sus condiciones derivadas
e índices de productividad, dan una indicación sobre el verdor actual de los biomas, así como sobre su situación
en comparación con el promedio a largo plazo. Otras variables físicas son el índice de área de la hoja, la
fracción de cobertura vegetal y la fracción de radiación absorbida para la fotosíntesis que cuantifican
respectivamente la densidad, la extensión y la salud de la vegetación. Además, la productividad de la materia
seca presenta el crecimiento de la biomasa en pie para aplicaciones agronómicas específicas. Finalmente, los
mapas de áreas quemadas delinean las zonas del globo afectadas por incendios.
Estos productos biofísicos son útiles para una amplia gama de áreas temáticas como el monitoreo global de
cultivos y las aplicaciones de seguridad alimentaria; manejo de bosques, agua y recursos naturales; modelado
de carbono en la tierra y pronóstico del tiempo y el clima.
Los mapas de cobertura del suelo representan información espacial sobre diferentes tipos (clases) de cobertura
física de la superficie de la Tierra como por ejemplo bosques, praderas, tierras de cultivo, lagos, humedales.
Los mapas dinámicos de la cobertura del suelo incluyen transiciones de clases de cobertura del suelo a lo largo
del tiempo y, por lo tanto, capturan los cambios de la cobertura del suelo. Los mapas de uso de la tierra
contienen información espacial sobre los arreglos, actividades e insumos que las personas realizan en un
determinado tipo de cobertura de la tierra para producirla, cambiarla o mantenerla.
Los mapas anuales de cobertura de la tierra de resolución moderada apuntan principalmente a la detección
de la cubierta de la tierra y sus cambios, aunque no es tan sencillo establecer límites entre las definiciones de
cobertura de la tierra y las clases de uso de la tierra. El mapa se proporciona junto con capas de campo continuo
de vegetación que proporcionan estimaciones proporcionales de la cubierta vegetal para varios tipos de
cobertura del suelo.
Formato de archivo:
Los productos globales son de 100m con teselado de 20x20 de resolución y se proporcionan como archivos
ZIP, que contienen un conjunto de 20 archivos GeoTIFF de banda única, uno por capa.
Sistema de proyección: UTM
Su precisión global general es de 80 %.
Para acceder a la información del Mapa Global de ocupación del suelo es necesario acceder al siguiente link:
https://lcviewer.vito.be/
Esta aplicación te permite realizar estadísticas sobre las diferentes coberturas terrestres de diferentes regiones
y países, además de mostrar en el mapa la información que el usuario solicite.
Además se pueden descargar los datos y verlos, comparar y exportar estadísticas de cobertura de la tierra por
país o región y hacer comparativas entre distintas ciudades.
En la figura 4 se observa un ejemplo del visor Global Land Cover Map del país de Guatemala.
Base de datos vectorial a escala Base de datos raster de las principales Base de datos raster de los
1:100.000 sobre la ocupación del suelo características de cobertura del suelo asentamientos humanos en Europa.
de Europa producida por el IGN a nivel de Europa verificada por el IGN a No es verificada ni producida por el
de España. nivel de España. IGN a nivel de España.
Con fecha del 27 de junio de 1985 y en virtud de una decisión del Consejo de la Unión Europea (CE/338/85),
se inicia “un proyecto experimental para la recopilación de datos, la coordinación y homogeneización de la
información sobre el estado del Medio Ambiente y los recursos naturales en la Comunidad”: este sería el
Programa CORINE (Coordination of Information of the Environment). En 1990, el Consejo de la Unión
Europea decidió crear la Agencia Europea de Medioambiente y el establecimiento de una Red de Información
y Observación del Medioambiente (EIONET). El objetivo de AEMA es proveer a la Unión Europea y a los
estados miembros de información fidedigna y comparable en el ámbito europeo que permita desarrollar
políticas y criterios de protección del Medio Ambiente, así como que dicha información esté disponible para
todos los posibles usuarios. CORINE pasa en 1995 a ser responsabilidad de la AEMA, la cual constituye, en
octubre de ese año, el Centro Temático Europeo “Land Cover” (ETC/LC), cuyo objetivo es el apoyo técnico
a la AEMA para la finalización, evolución, mantenimiento y actualizaciónde las bases de datos en materia de
“Land Cover”. Con la constitución del Programa Copernicus en 2014 y la coordinación de la AEMA del
Servicio de Territorio, CORINE se enmarca dentro del Servicio Territorio y su Componente Pan-Europeo.
El objetivo fundamental del proyecto CORINE Land Cover (CLC) es la creación de una base de datos a escala
1:100.000 sobre la ocupación del suelo de Europa con unas pautas establecidas de homogeneidad para toda
Europa (mismas especificaciones, misma escala y clases temáticas), fecha de referencia comparable para toda
Europa, producido por los estados miembros y coordinado por la AEMA.
Esta producción fue llevada a cabo por fotointerpretación tradicional desde las imágenes de referencia en las
versiones del CLC 1990, 2000 y 2006. Sin embargo, desde la existencia de datos de ocupación del suelo más
detallados a nivel nacional procedentes del proyecto SIOSE, la producción del CLC 2012 y CLC 2018 se
vienen realizando por un enfoque productivo de abajo-arriba (bottom-up) mediante generalización geométrica
y semántica del SIOSE a escala 1:100.000. Esto asegura una mayor coincidencia en los datos y estadísticas de
ambos proyectos, transmisión del conocimiento nacional a los datos continentales y una reducción de esfuerzos
y costes.
Figura 5: SIOSE transformado a nomenclatura CLC sin generalizar (izquierda) y CLC obtenido por generalización
(derecha).
La viabilidad futura del proyecto CORINE Land Cover está asegurada pues proporciona datos fundamentales
para el entendimiento del territorio, además posee una cadena temporal de versiones que se remonta hasta el
año 1990, muy útil para el estudio de la dinámica del territorio. Sin embargo ya desde hace años se detectó,
principalmente por usuarios nacionales, la demanda de un detalle geométrico y temático mayor. Esto ha
desembocado en un nuevo paradigma para el proyecto CORINE Land Cover siendo la AEMA y el grupo de
expertos internacionales EAGLE de su desarrollo y evolución.
Ciertas de estas limitaciones del CORINE Land Cover ya fueron detectadas por los usuarios españoles en la
versión del CLC 2000, por lo que se decidió extender el número de clases de la nomenclatura de 44 a un total
de 85, alcanzando un además de 5 niveles jerárquicos. Sin embargo los resultados de esta experiencia no fueron
tan satisfactorios como se esperaban no volviendo a repetirse. Lo que implicó la creación de un nuevo proyecto
de ocupación del suelo más detallado en geometría y temática, el SIOSE.
Son un juego de imágenes, actualmente, sobre cinco variables de ocupación del suelo: suelo sellado
(Imperviousness), fracción de cabida cubierta arbórea y tipo de hoja (Tree Cover Density and Dominant Leaft
Type), vegetación herbácea natural (Grassland) y masas de agua y zonas húmedas (Water and Wetness). Son
producidas centralizadamente por la AEMA y verificadas y/o mejoradas por los estados miembros.
Dependiendo de la capa, existen versiones de los años 2006, 2009, 2012, y 2015.
El Instituto Geográfico Nacional (Ministerio de Fomento) tiene encomendada la verificación de las HRL a
partir de fuentes de referencia nacionales entre ellas el SIOSE, PNOA, LIDAR e información geográfica de
referencia (hidrografía, poblaciones).
Otro producto dentro de la componente paneuropea que no es producido ni verificado por parte del IGN es el
Mapa Europeo de Poblaciones. El Mapa europeo de poblaciones es un conjunto de datos de tipo raster que
mapea los asentamientos humanos en Europa sobre la base de imágenes de satélite SPOT5 y SPOT6. Su
versión del 2016 representa el porcentaje de cobertura de área construida por unidad espacial y describe el
asentamiento humano.
Base de datos vectorial de áreas urbanas Base de datos vectorial de ocupación del Base de datos vectorial sobre la
de toda Europa con población mayor de suelo en áreas ribereñas de Europa ocupación del suelo de áreas de
50.000 habitantes. Verificado y mejorado verificada por el IGN a nivel de España. protección medioambiental de Europa
por el IGN/CNIG a nivel de España. verificada por el IGN a nivel de España.
Ofrece bases de datos vectoriales de ocupación del suelo y usos del suelo de alto detalle de áreas urbanas de
todas Europa designadas por Copernicus con una población mayor a 50.000 habitantes (Funtional Urban Areas
– FUA). Su unidad mínima de resolución es 0.25 ha para áreas urbanas y 1 ha para áreas rurales, haciendo uso
de una nomenclatura de 27 clases (17 urbanas y 10 rurales). Existen versiones de los años 2006 y 2012. En la
versión 2012 hay un total de 785 FUAs a nivel de Europa. Son producidas centralizadamente por la AEMA y
verificadas y/o mejoradas por los estados miembros, en el caso de España por el IGN/CNIG.
Ofrece bases de datos vectoriales de ocupación del suelo de alto detalle de zonas ribereñas para los cauces más
importantes. Representan áreas de transición que ocurren entre la tierra y los ecosistemas de agua dulce.
Proporcionan una amplia gama de funciones ribereñas (por ejemplo, filtración química, control de
inundaciones, estabilización de bancos, vida acuática y apoyo de la fauna ribereña, etc.) y servicios de
ecosistemas.
Los productos de Zonas ribereñas respaldarán los objetivos de varios actos legales europeos e iniciativas
políticas, como la Estrategia de Biodiversidad de la UE para 2020, las Directivas de Hábitats y Aves y la
Directiva Marco del Agua.
Su unidad mínima de resolución es 0.25 ha y hace uso de una nomenclatura jerárquica de 4 niveles y 80 clases.
Existe la versión de 2012. Son producidas centralizadamente por la AEMA y verificadas y/o mejoradas por
los estados miembros.
Natura 2000 es una red europea de áreas de protección medioambiental establecida por la Directiva de Hábitat
en 1992 y representa una pieza clave para las políticas europeas de naturaleza y biodiversidad. Copernicus
ofrece bases de datos vectoriales de ocupación del suelo de alto detalle de una selección de las áreas de Natura
2000. Su unidad mínima de resolución es 0.25 ha y hace uso de una nomenclatura de MAES y CLC con un
total de 62 clases. Existen la versiones de 2006 y 2012. Son producidas centralizadamente por la AEMA y
verificadas y/o mejoradas por los estados miembros.
Las imágenes de satélite forman la entrada para la creación de muchos productos y servicios de información,
como mapas de cobertura terrestre o capas de alta resolución sobre características de cobertura terrestre. Con
todas las imágenes de satélites disponibles para cubrir 39 países del EEE (EEA39), las escenas de imágenes
individuales se han procesado en mosaicos orto-rectificados paneuropeos sin interrupciones. El Servicio de
Monitoreo de la Tierra de Copernicus también brinda acceso al servicio Mosaico Global Sentinel-2.
EU-DEM es un modelo de superficie digital paneuropeo y EU-Hydro proporciona una red fluvial
fotointerpretada y una red de drenaje modelada para los países del EEA39.
Muchos datos in situ se gestionan y se hacen accesibles a nivel nacional. Sin embargo, debido a problemas
como el acceso a los datos y las restricciones de uso, la calidad de los datos y la disponibilidad en los países
del EEE39, los servicios de Copernicus y, en particular, el Servicio de Monitoreo de la Tierra de Copernicus
también se basa en conjuntos de datos in situ paneuropeos creados y / o coordinados a nivel europeo.
La base de datos LUCAS es una base de datos europea coordinada por la EUROSTAT, que se utiliza para la
verificación y validación de varios servicios de información del Servicio de Monitoreo de la Tierra de
Copernicus. Se centra en el estado y la dinámica de los cambios en el uso y la cobertura de la tierra en la Unión
Europea. Se realizan y registran un gran número de observaciones en toda la UE. En 2012, los 27 países de la
UE han analizado más de 270,000 puntos en diferentes tipos de cobertura terrestre (tierras de cultivo, praderas,
bosques, áreas urbanizadas, red de transporte, etc.). En estos puntos, los agrimensores han examinado la
cobertura del suelo y el uso del suelo, la gestión del riego y los elementos estructurales del paisaje. Además,
se toma una muestra de 500 gr de suelo superior en uno de cada 10 puntos. Estas muestras se analizan en un
laboratorio y se utilizan con fines relacionados con la evaluación de factores ambientales, como la
actualización de mapas europeos del suelo, la validación de modelos de suelo y la medición de la cantidad de
carbono orgánico en el suelo, que es un factor importante que influye en el cambio climático.
Urban Atlas Servicio Territorio, AEMA Estados Miembros Vectorial 2006, 2012 0,25 ha
Componente Local polígonos
Natura 2000 Servicio Territorio, AEMA Estados Miembros Vectorial 2006, 2012 0,25 ha
Componente Local polígonos
Tabla 1: Resumen de los principales productos del Servicio Land Monitoring de Copernicus
https://lcviewer.vito.be/download
Componente Pan-Europeo
Figura 9: Imagen del acceso a cada una de las tipologías de HRL y CLC en la página web del Servicio de
Territorio de Copernicus
HRL:
https://land.copernicus.eu/pan-european/high-resolution-layers
Figura 10: Descarga CNIG de los productos del componente local de España
Natura 2K 2012
http://centrodedescargas.cnig.es/CentroDescargas/buscadorCatalogo.do?codFamilia=NAT2K
Figura 11: Descarga productos del componente local de la página del Servicio del Territorio de Copernicus
Natura 2K 2012
https://land.copernicus.eu/local/natura/natura-2000-2012
El resto de productos de Copernicus europeos los puedes encontrar en la web de cada uno de los
servicios.
Se materializa como un escenario innovador en la producción del siguiente CORINE Land Cover mediante la
armonización de los datos nacionales y una posterior integración a nivel europeo (enfoque de abajo-arriba).
Una herramienta de comparación semántica entre las definiciones de clase de los diferentes sistemas
de ocupación del suelo por la descomposición de las clases en elementos de cubiertas, usos del suelo
y características del paisaje.
Un modelo de datos conceptual, en notación UML. Siguiendo ISO19109 (Información Geográfica -
Reglas para el esquema de aplicación), restantes estándares ISO y como extensión de las
especificaciones INSPIRE dando servicio concreto a la ocupación del suelo en Europa.
Estructura de base de datos de acuerdo al modelo, que permite almacenar y mantener los datos, ser
consultados fácilmente por una amplia gama de usuarios con propósitos variados.
Conjunto de reglas de generalización y agregación con el fin de construir un escenario europeo común
en la generación de bases de datos continentales desde conjuntos de datos nacionales más detallados.
Materialización del enfoque de abajo-arriba.
Los mapas anuales de cobertura de la tierra de resolución moderada apuntan principalmente a la detección de
la cubierta de la tierra y sus cambios, aunque no es tan sencillo establecer límites entre las definiciones de
cobertura de la tierra y las clases de uso de la tierra. El mapa se proporciona junto con capas de campo continuo
de vegetación que proporcionan estimaciones proporcionales de la cubierta vegetal para varios tipos de
cobertura del suelo.
Formato de archivo:
Los productos globales son de 100m con teselado de 20x20 de resolución y se proporcionan como archivos
ZIP, que contienen un conjunto de 20 archivos GeoTIFF de banda única, uno por capa.
Sistema de proyección: UTM
Su precisión global general es de 80 %.
Una vez accedamos a la plataforma investigaremos los widgets de la plataforma para ver qué información nos
aporta la aplicación sobre el Global Land Cover Map.
Navegamos por las herramientas.
1-)¿Puedo visualizar Land Cover de diferentes formas? Si es que si, ¿cómo se muestra la clasificación de Land
Cover en el visor?
Hay 3 formas diferentes:
3-) Calcula las estadísticas de Guatemala región Sacatepéquez para ‘Single Class Forest’
4-) Compara las estadísticas anteriores con las de la Comunidad de Madrid en España y descarga el informe.
5-) Descarga imágenes de la ciudad que elijas. ¿en qué formato se bajan y con qué resolución?
Imágenes tiff de cada una de las clases de Global Land Cover Map con una resolución de 100 m.
En 1990, el Consejo de la Unión Europea decidió crear la Agencia Europea de Medioambiente (EEA) y el
establecimiento de una Red de Información y Observación del Medioambiente (EIONET). El programa
CORINE Land Cover (en adelante CLC) pasa en 1995 a ser responsabilidad de la EEA y de los Centros
Nacionales de Referencia en Ocupación del Suelo en cada uno de los estados miembros, que como hemos
visto, en España es el Instituto Geográfico Nacional.
En sus orígenes el programa CLC englobaba la adquisición de imágenes y la realización del primer sistema
de información sobre ocupación del suelo en Europa, con fecha de referencia 1990. Se han realizado
actualizaciones posteriores, siendo ya las últimas versiones coordinadas por el Servicio de Territorio de
Copernicus (Land Monitoring Service) en su Componente Pan Europeo.
El CORINE Land Cover, proporciona información sobre la cubierta terrestre, el uso del suelo y sus cambios
temporales. Entre sus aplicaciones se encuentra la gestión del territorio, la gestión forestal, la gestión del
agua, la gestión de emergencias, etc.
Principales características:
Es un producto homogéneo para toda Europa, es decir, con las mismas especificaciones y clases y
misma fecha de referencia.
Tiene la escala y clases temáticas necesarias para la toma de decisiones medioambientales, obtención
de variables de expansión urbana, cuantificación de la pérdida de entornos naturales,
deforestaciones e incendios, identificación de zonas de protección, etc a nivel europeo.
Producido por los países miembros y coordinado por la Comisión Europea.
Versiones: 1990, 2000, 2006, 2012 y 2018
Sostenible económicamente y actualizado periódicamente.
Frecuencia actualización de seis años.
El futuro del proyecto CORINE Land Cover está asegurado y actualmente desde Copernicus se están
desarrollando nuevas metodologías y nuevos productos CLC que satisfagan las necesidades actuales de un
mayor detalle geométrico y temático.
2.1 Nomenclatura
La nomenclatura jerárquica CLC contiene 44 clases divididas en tres niveles de detalle temático. Esta
definición de clases es fruto del trabajo conjunto de expertos de distintos países europeos con el propósito de
describir todos los paisajes existentes en Europa.
En el primer nivel se definieron cinco grandes grupos de clases: zonas artificiales, cultivos, vegetación y
zonas naturales, zonas húmedas y zonas de agua. Estas clases se subdividen en un segundo y tercer nivel de
detalle, debiendo clasificar cada polígono CLC al máximo nivel de detalle.
Cada clase se identifica con un código de tres dígitos, el primer dígito es la clase al 1º nivel, el segundo al 2º
nivel y el tercero al 3º nivel. Por ejemplo, el código 312 corresponde a la clase 3 en 1º nivel (vegetación y
zonas naturales), a la clase 1 en 2º nivel (bosques) y a la clase 2 en 3º nivel (bosques de coníferas).
En la siguiente Tabla 1 Nomenclatura CLC se listan todas las clases CLC al tercer nivel.
La nomenclatura CLC se puede seguir extendiendo en nuevas clases si un país así lo necesita, como fue el
caso de España en el CLC del año 2000, donde la nomenclatura se extendió a cinco niveles, hasta que la
aparición del proyecto SIOSE cubrió las necesidades de un mayor detalle temático y no hizo falta volver a
trabajar con más niveles de jerarquía.
La producción del CLC en España también se realizó por fotointerpretación y creación de bases de datos de
cambio sobre imágenes de referencia en las versiones del CLC 1990, 2000 y 2006, pero en CLC 2012, al
existir datos de ocupación del suelo más detallados a nivel nacional del proyecto SIOSE, la producción se
realizó mediante generalización geométrica y semántica del SIOSE a escala 1:100.000. Posteriormente el
CLC 2018 se actualizó sobre esta base de datos generalizada.
Esta metodología de trabajo asegura una mejor interoperabilidad entre ambas bases de datos, una mayor
coherencia en los resultados estadísticos y una reducción de esfuerzos y costes.
https://land.copernicus.eu/pan-european/corine-land-cover
Se selecciona el año para el producto final CLC como “CLC 1990” o “CLC 2000” o se puede seleccionar
también de manera independiente la capa de cambios entre dos fechas de referencia, como “CHA 2000-
2006”, la capa de cambios entre 2000 y el año 2006.
Una vez seleccionado el producto y la versión, se abre la ventana de gestión, donde es posible visualizar la
capa, ver sus metadatos, ver la leyenda, activar capas adicionales, descubrir las URLs de los servicios web
relacionados y descargar los datos.
Accediendo a la opción de descarga, es posible listar cada una de las partes o conjuntos de datos
descargables ligados a los productos. Para la descargar es necesaria la creación de un usuario y
contraseña.
*****
Por su parte, La Agencia Europea de Medio Ambiente pone también a disposición de los usuarios servicios
mediante el portal “DiscoMap” y se describen dentro de la plataforma “CORDA”.
Una vez en el Centro de Descargas, hay que seleccionar el tipo de información a descarga, en este caso CLC
se encuentra dentro del grupo “Información Geográfica Temática”, marcada en rojo en la Figura 6 Imagen
del Centro de Descargas
Se selecciona el producto CORINE Land Cover, y elegimos la opción de visualización de datos para
descarga “Por listado”, marcada en rojo en la Figura 7 Imagen del producto CORINE Land Cover a descarga
Aparece una nueva ventana de “Búsqueda por listado” donde hay que elegir dentro de la pestaña “División
Administrativa” elegimos para la primera casilla la división administrativa de “Comunidad Autónoma” y en
En una nueva ventana aparecen todos los productos CLC disponibles. Podemos filtrar por “Temática”
eligiendo el año y por “Tipo de fichero” eligiendo el tipo de formato.
Podemos dar al icono de “Descargar” para iniciar la descarga de los datos o “Añadir” para añadir a la cesta el
producto y buscar nuevas descargas.
Figura 9 Imagen del filtro de resultados por año 2012 y tipo de fichero GeoPackage.
Adicionalmente el IGN / CNIG ha generado un servicio web de visualización para contener todos aquellos
productos validados en su cometido como Centro Nacionales de Referencia de Cubiertas el Suelo. Su
dirección de acceso es la siguiente.
http://servicios.idee.es/wms/copernicus-landservice-spain
El Componente Continental proporciona productos de ámbito continental sobre la ocupación del suelo y su
evolución temporal, entre los que se encuentra las High Resolution Layers o Capas de Alta Resolución sobre
la cubierta del suelo.
Estas capas ofrecen información específica de características de la cubierta terrestre y funcionan como
complemento a otros productos de Copernicus como el CORINE Land Cover. Las HRL son imágenes de 20m
de resolución, derivadas por procesamientos automáticos e interactivos de imágenes de satélite. Por este
motivo no ofrecen una visión completa de la Ocupación del Suelo como puede ser el CORINE Land Cover, si
no que se cada una de ellas se centra en un aspecto concreto observable.
Son un juego de imágenes sobre la siguientes variables: suelo sellado (impreviourness), fracción de cabida
cubierta arbórea y tipo de bosque (tree cover density and forest type), vegetación herbácea natural (grassland),
zonas húmedas (wetlands), masas de agua (water bodies) o pequeños elementos arbóreos (small woody
features). Son producidas centralizadamente por la AEMA, y verificadas y/o mejoradas por los estados
miembros.
https://land.copernicus.eu/pan-european/high-resolution-layers
Figura 1: Ejemplo de HRL del suelo sellado sobre mosaico SPOT en Valladolid (izquierda), y HRL de
densidad arbórea e imagen SPOT de Portillo, Valladolid (derecha.).
A lo largo de las diferentes versiones que han tenido las HRL, éstas han ido evolucionando en contenido y
metodología de obtención. Por este motivo, es posible encontrar capas nuevas según la versión o algunas otras
que no se han continuado. Un motivo por el cual aparece una nueva capa es la necesidad de tener un inventario
continental de ese tipo de cubierta del suelo. Por el contrario un motivo por el cual las capas dejan de producirse
es el inherente a la fiabilidad de sus resultados. A continuación se muestra un resumen de las capas actualmente
disponibles a descarga de las versiones 2006, 2009, 2012 y 2015.
Figura 2: Ejemplo del aspecto de capa de suelo sellado sobre ortofotografía y elementos topográficos básicos.
Se puede ver como la intensidad del color rojo indica el grado del suelo sellado (1-100).
El producto de fracción o densidad de cubierta arbórea consiste en una capa que muestra la de densidad de esta
variable de la ocupación del suelo un rango de 0-100%. Se encuentra disponible para los años de referencia
2012 y 2015, y un producto de cambio que muestra el aumento o la disminución de los cambios reales de entre
2012-2015.
Esta capa es un producto que proporciona información sobre el tipo de hoja dominante: frondosa o coníferas,
disponible para los años de referencia 2012 y 2015.
Este producto consta de: Una capa de tipo de hoja dominante que tiene una unidad mínima cartografiable de
0.5 ha, así como un umbral de densidad de cobertura de árbol del 10%. Y una segunda capa adicional que se
asigna, según el tipo de hoja dominante, árboles de uso agrícola y en contexto urbano (derivados de los datos
del CORINE Land Cover y la HRL de sellado del suelo. Ambas capas están disponibles para los años de
referencia 2012 y 2015.
Esta capa ha cambiado su definición con respecto a la versión anterior de 2012. El producto de vegetación
herbácea del año de referencia 2015 es un nuevo producto, que reemplaza por completo a todos los productos
del año de referencia 2012 anterior. La capa de 2015 no es directamente comparable con la de 2012 por su
metodología distinta de producción. Como reacción a las lecciones aprendidas de la producción de 2012, el
enfoque se modificó drásticamente y el producto de 2015 se basa en una serie de imágenes más largas de varios
sensores diferentes. La serie temporal más larga permite la detección de eventos de arado en los años anteriores
al año de referencia, lo que ayuda a distinguir los cultivos de los pastizales. Además por primera vez se combina
datos de radar (SAR) y datos ópticos para mejorar la precisión de la clasificación.
También para esta versión 2015 se ha conseguido generar capas adicionales sobre índices de la probabilidad
la vegetación herbácea y un indicador de arado.
Esta capa de agua y humedad del año de referencia 2015 es un nuevo, que reemplaza por completo a los
productos de "agua permanente" y "zonas húmedas" por separado que se generaron en 2012.
El producto combinado Agua y Zonas Húmedas es una capa temática que muestra la ocurrencia de agua y
superficies húmedas durante el período de 2009 a 2015. Esta capa se basa en imágenes de satélites ópticas
multitemporales y multiestacionales. Además, esta capa también se basa en información de radar (imágenes
Sentinel-1). Se utiliza una multitud de imágenes ópticas y radar, que cubren una serie de tiempo prolongada a
lo largo de 7 años, cuyo objetivo es capturar la dinámica interanual dentro de un área determinada y alcanzar
una imagen compuesta por temporada (de 3 meses) y año.
Producto: Índice de probabilidad de agua y zonas húmedas (water and wetness probability
indicator - WWPI)
Resolución: 20m píxel
Sistema de referencia y proyección: ETRS89 LAEA
Formato: GEOTIFF 8bit unsigned, comprimido LZW
Valores de los píxeles
0: áreas sin agua o humedad
1-100: índice de probabilidad de agua y zonas húmedas
254: áreas no clasificadas (falta de dato en la imagen, p.e. nubes, sombras, etc.)
255: áreas fuera de la imagen
Desde la producción del año de referencia 2015, la producción se basa cada vez más en series temporales de
imágenes de satélite de varios sensores, incluida la combinación de datos ópticos y de radar. La fuente principal
son los satélites Sentinel de Copernicus (en particular Sentinel-2 y Sentinel-1), además de los mosaicos de
imágenes de alta resolución (HR 2006-09-12-15) y desde 2015, también se utilizan imágenes de muy alta
resolución (VHR 2015) para algunos de los productos.
Las HRL no ofrecen una visión completa de la Ocupación del Suelo como puede ser el CORINE Land Cover,
si no que se cada una de ellas se centra en un aspecto concreto observable. También es importante conocer que
los tratamientos de imágenes de cada tipo de HRL es diferente, por lo que los resultados no nos integrables en
una única HRL, o lo que es lo mismo, los píxeles de las HRL no son excluyente entre ellas (p.e. un mismo
píxel de forest puede estar también recogido en grassland)
A pesar de la documentación disponible sobre la metodología de obtención de las HRL, sigue siendo escasa
para su conocimiento completo o reproducción por parte de los estados miembros.
El aseguramiento de la calidad en unos productos como las capas HRL es una pieza muy importante. Debido
a que la mayoría de los píxeles vienen derivados de la integración o la clasificación de imágenes, su resultado
es automático en una componente muy alta. Por este motivo son necesarias un conjunto de tareas de evaluación
de la calidad para ofrecer una respuesta completa al usuario.
Por definición y como premisa en la producción de las HRL, existe la necesidad de alcanzar un valor fiabilidad
de 90%. O lo que es lo mismo, las muestras utilizadas para la evaluación deben ser correctas al menos en un
90% de los casos. Este requerimiento y las siguientes tareas de verificación descritas no se realizan sobre todas
las HRL, sólo en aquellas consideradas más importantes: IMD, TCD, DLT, GRA y WAW.
Tareas de verificación:
a) Estudio interno realizado durante la producción de las capas y adjunto al reporte de entrega. Lo realiza
la empresa contratada por Copernicus.
b) Estudio independiente de análisis de la fiabilidad realizado a nivel completo del mosaico europeo.
Realizado por equipos técnicos independientes a la producción.
c) Estudio independiente de análisis de la fiabilidad realizado a nivel nacional y con estadísticas
nacionales. Realizado por los Centro Nacionales de Referencia de Cubiertas el Suelo de la Red
EIONET (EIONET NRC-Land Cover)
Los resultados de la verificación interna “a” de la versión 2012 están accesibles a descarga junto con los datos
de las HRL en la web de Copernicus Land. A día de redacción de este informe no se encuentran accesibles los
informes de la versión 2015.
https://land.copernicus.eu/user-corner/technical-library
Los resultados de la validación interna “a” no son de fácil entendimiento debido a su extensión metodológica
que hace enmascarar los aspectos fundamentales en la validación. Los valores de la principal componente
validad “Thematic accuracy” se distribuye irregularmente entre la tipología de capas y las áreas geográficas.
Por lo que es posible adivinar unos valores de fiabilidad fluctuantes entre un 50% a un 90%. En el caso de que
se desee una información más detallada en este aspecto, se recomienda acceder a los informes de validación
interna.
Por el contrario los resultados de las verificaciones independientes, a sea nivel continental “b” o nacional “c”
no se encuentra a libre disposición en la web, pero pueden ser solicitados a Copernicus o a los NRC-LC
nacionales. Resulta de especial interés explicar brevemente las tareas de verificación nacional para comprender
estas tareas de verificación.
La verificación nacional “c” realizada por los NRC-LC para cada una de las HRL se realizan en 3 pasos:
1) General overview (cualitativo, obligatorio): intersección con la mejor base de datos disponible
temática a nivel nacional y análisis de las áreas con grandes discrepancias. El resultado de este paso
son estadísticas porcentuales globales para todo el país de la correspondencia de los píxeles de las
HRL con los datos nacionales.
2) Look and feel verification (cualitativo, obligatorio): revisión de zonas seleccionadas como críticas para
los errores de comisión y omisión. En este caso se seleccionan áreas de control condicionadas por un
error esperado y se evalúa cualitativamente la fiabilidad de las HRL. El resultado en este paso es una
estimación cualitativa de la calidad.
Figura 8: Ejemplo de área de control condicionada por posibles ausencias de vegetación en dehesas. Esta
área puede ser valorado cualitativamente según el comportamiento de la HRL: coincidencia, omisión o
comisión. Partiendo de unas áreas de control es posible obtener estimaciones cualitativas.
La principal faceta de los datos, imágenes y la mayoría de los servicios de Copernicus es la facilidad de acceso
a los mismos. A continuación se muestran los nodos de acceso a los datos de la HRL disponibles.
Las HRL forman parte del Componente Pan-Europeo y desde la propia página web es posible descargarse los
datos de las capa de alta resolución.
https://land.copernicus.eu/pan-european/high-resolution-layers
Figura 10: Imagen del acceso a cada una de las tipologías de HRL en la página web del Servicio de Territorio
de Copernicus
Seleccionando cada uno de los tipos de HRL se puede acceder a la descarga. Como se ha comentado
previamente, cada tipología puede tener más de un producto. Es precisamente en esta navegación donde pueden
encontrarse.
Figura 11: Ejemplo de navegación por la capa de suelo sellado, donde se identifican 2 tipos de productos,
“Status Maps”, o “Change Maps”. Tras la selección de uno de ellos, es posible seguir accediendo a cada uno
de las versiones del producto.
Una vez seleccionado el producto y la versión a acceder, se abre la ventana de gestión, donde es posible
visualizar la capa, ver sus metadatos, ver la leyenda, activar capas adicionales, descubrir las URLs de los
servicios web relacionados y descargar los datos.
En función de la versión del producto es posible que estén desactivadas algunas de las acciones previamente
comentadas.
Accediendo a la opción de descarga, es posible listar cada una de las partes o conjuntos de datos descargables
ligados a los productos. En el caso concreto de las HRL es posible descargarse los datos fraccionados por una
rejilla de 100km x1000km pan-europea, por lo tanto será previamente necesario elegir la celdilla deseada. Las
HRL pueden encontrarse a resoluciones originales de 20 m y remuestreadas a 100 m. Es muy importante
conocer que los datos de las HRL descargables desde la web oficial del Servicio de Territorio de Copernicus
se encuentran en el sistema de referencia y proyección europea ETRS89 LAEA (EPSG:3035) y que es
necesario tener un usuario en la web para su descarga.
Figura 13: Iquierda: ejemplo de listado de datos HRL a descarga. Derecha: rejilla de 1000 km x 1000 km en
los que se ofrecen los datos HRL.
http://discomap.eea.europa.eu/
https://corda.eea.europa.eu/
El Instituto Geográfico Nacional, como Centro Nacional de Referencia de Cubiertas el Suelo de la Red
EINOET, es el encargado en la validación nacional de las HRL. Adicionalmente el IGN ha publicado mediante
el Centro de Descargas del CNIG los datos de las HRL 2015 comentadas en este documento.
http://centrodedescargas.cnig.es/
Desde la pestaña “Información Geográfica Temática” es posible encontrar todos los productos de ocupación
del suelo en el que el IGN participa (coproduce o valida), entre ellos las HRL. Aspecto a destacar es que desde
este nodo de descarga es posible conseguir los datos de las HRL a 20 m de resolución y en el sistema de
referencia y proyección nacional ETRS89 UMT30 para la Península y Baleares, y WGS84 UTM28 para
Canarias. La documentación anexa a los datos se encuentra directamente contenidas en los ficheros de
descarga.
Figura 14: Ventana de acceso a descarga de los datos y documentación sobre las HRL en el Centro de
Descargas del CNIG
Adicionalmente el IGN / CNIG ha generado un servicio web de visualización para contener todos aquellos
productos validados en su cometido como Centro Nacionales de Referencia de Cubiertas el Suelo. Su dirección
de acceso es la siguiente.
http://servicios.idee.es/wms/copernicus-landservice-spain
Urban Atlas
El Urban Atlas es una iniciativa conjunta de la Dirección General de Política Regional y Urbana de la Comisión
Europea y la Dirección General de Empresa e Industria en el marco del programa Copernicus de la Unión
Europea con el apoyo de la Agencia Espacial Europea y la Agencia Europea de Medio Ambiente.
Este producto se localiza dentro del Servicio del Territorio de Copernicus en el componente local junto con
Riparian Zones y Natura 2000 y futuro Coastal Zones. Es coordinado y producido por la Agencia Europea
de Medio Ambiente (AEMA) y verificada y/o mejorada por los estados miembros, en el caso de España por
el IGN/CNIG.
El objetivo de Urban Atlas es proporcionar específica información detallada que se complementa con la
información obtenida en el componente Pan Europeo del servicio de Territorio de Copernicus.
Existen dos versiones de Urban Atlas: Urban Atlas 2006 y Urban Atlas 2012 y próximamente Urban Atlas
2018.
A partir de Urban Atlas 2012 se han creado los siguientes productos que se explican a continuación: Street
Tree Layer (STL), Cambios de Urban Atlas entre 2006 y 2012 y Alturas de edificios 2012.
El UA 2006 es la primera versión creada de Urban Atlas. Está definido por áreas urbanas funcionales (FUAs)
con más de 100.000 habitantes según lo definido por el Urban Audit. Está formado por un total de 319 FUAs
europeas y por una nomenclatura de 17 clases urbanas y 3 clases rurales. Los datos vectoriales se pueden
descargar en formato shape en la siguiente URL como se explicará en apartados posteriores:
https://land.copernicus.eu/local/urban-atlas/urban-atlas-2006
Su nomenclatura definida es 17 clases urbanas con una unidad mínima de resolución de 0.25 ha y 10 clases
rurales con una unidad mínima de resolución de 1 ha. Están basadas en la nomenclatura del CORINE Land
Cover y en la leyenda del servicio urbano de Copernicus. Su leyenda es oficial definida en la documentación
técnica oficial de Urban Atlas.
El sistema de referencia del Urban Atlas es el ETRS 1989 en proyección LAEA (Lambert Acimutal área).
UA2012
1 Superficies artificiales
1.1 Tejido urbano
1.1.1 Tejido urbano continuo (IMD>=80%)
1.1.2 Tejido ubano discontinuo (IMD 10%-80%)
1.1.2.1 Tejido ubano discontinuo denso (IMD 50%-80%)
1.1.2.2 Tejido ubano discontinuo densidad media (IMD 30%-
50%)
1.1.2.3 Tejido ubano discontinuo densidad baja (IMD 10%-30%)
1.1.2.4 Tejido ubano discontinuo densidad muy baja (IMD
<10%)
1.1.3 Estructuras aisladas
1.2 Unid industriales, comerciales, públicas, militares, priv
y de transporte
1.2.1 Unidades industriales, comerciales, públicas, militares
y privadas
1.2.2 Carreteras, red de ferrocarriles y terrenos asociados
1.2.2.1 Carreteras de tránsito rápido y terrenos asociados
1.2.2.2 Otras carreteras y terrenos asociados
1.2.2.3 Ferrocarriles y terrenos asociados
1.2.3 Puertos
1.2.4 Aeropuertos
1.3 Minas, vertederos y zonas en construcción
1.3.1 Zonas de extracción de minerales y vertederos
1.3.2 Zonas en construcción
1.3.3 Terrenos que no están en uso
1.4 Areas artificiales con vegetación no relacionadas con
agricultura
1.4.1 Zonas verdes urbanas
1.4.1 Instalaciones deportivas y de ocio
2 Agricultura
2.1 Tierras de cultivo (cultivos anuales)
2.2 Cultivos permanentes
Los datos vectoriales de las 785 FUAs se pueden descargar en formato shape en la siguiente URL:
https://land.copernicus.eu/local/urban-atlas/urban-atlas-2012
También se pueden descargar las FUAs españolas en el centro de descargas del IGN/CNIG:
http://centrodedescargas.cnig.es/CentroDescargas/buscadorCatalogo.do?codFamilia=URBAT
La capa Street Tree (STL) es una capa separada de la capa UA2012 producida dentro de la máscara urbana de
nivel 1 para cada FUA. Incluye filas contiguas o parches de árboles que cubren 500 m² o más y con un ancho
mínimo de 10 m sobre “Superficies artificiales” (nomenclatura clase 1) dentro de FUA, es decir, hileras de
árboles a lo largo de la red de carreteras fuera de las áreas urbanas o el bosque adyacente a las áreas urbanas
no debe ser incluido.
La unidad mínima de resolución es de 0.05 ha y de anchura es de 10 metros con algunas excepciones descritas
en las especificaciones del producto. Su nomenclatura es una única clase que indica la presencia de árboles o
de no de árboles.
Esta capa utiliza las imágenes satelitales orto-rectificadas VHR 2006 y 2012 para detección de cambios de
imagen a imagen, y desde el cual compara los cambios con el Urban Atlas 2006 que tiene que ser mapeado.
Un producto intermediario, el “UA 2006-2012 Combined LULC", combina la información de LULC por dos
períodos (por ejemplo, 2006 y 2012). Basado en la primera Capa LULC (por ejemplo, 2006), esta capa se
produce al delinear polígonos donde la clase LULC tiene cambiado en el segundo período (por ejemplo, 2012).
Finalmente, la "capa de cambio LULC" se deriva de esta capa combinada y solo incluye polígonos con cambios
de Uso del suelo / Cobertura del suelo entre los dos periodos. La detección de cambios se aplicará con una
MMU diferente de la capa LULC.
Con el fin de garantizar que los cambios relevantes de LULC sean apropiados, la MMU para la capa de cambio
LULC (por ejemplo, 2006-2012) se define a continuación:
Las excepciones a este límite de MMU se pueden encontrar en los polígonos de cambio donde las redes de
carreteras y ferrocarriles (las clases 12210, 12220, 12230) están presentes en una fecha (por ejemplo, 2006 o
2012) o polígonos directamente conectado a estas redes para mantener la consistencia de estas redes.
Los datos vectoriales de las FUAs con cambios se pueden descargar en formato shape en la siguiente URL:
https://land.copernicus.eu/local/urban-atlas/change-2006-2009
Urban Servicio Territorio, Agencia Estados Miembros Vectorial 2006, 2012 0,25 ha 1/10000
Atlas Componente Local Europea de (España en el caso polígonos
Medio del UA en
Ambiente España)
Street Tree Servicio Territorio, Agencia Estados Miembros Vectorial 2012 0,05 ha 1/10000
Layer Componente Local Europea de (España en el caso
Medio del UA en
Ambiente España)
1. La segmentación y clasificación automatizadas para lograr una diferenciación inicial entre las clases de
cobertura terrestre (urbana, bosque, agua u otras coberturas terrestres) es posible siguiendo una decisión de los
proveedores de servicios.
2. La impermeabilidad de la capa de alta resolución (HRL) se utiliza para clasificar las densidades de sellado
del tejido urbano clase 1.1 para los niveles 3 y 4. La asignación de los niveles de impermeabilidad / sellado (es
decir, las clases 1.1.2.1 - 1.1.2.4) será llevado a cabo utilizando la capa de impermeabilidad HRL.
3. La delineación se debe hacer en los datos de datos de observación de la tierra. Los datos de EO deben ser
considerados como el principal Fuente de datos (guía).
o Los datos de EO, mapas topográficos, datos de navegación (COTS u OSM) y otros datos
relevantes. datos auxiliares;
o Información auxiliar incluida la experiencia local.
5. El área interpretada debe interpretarse con una extensión mínima de 100 m (100 m de búfer) para asegurar
la precisión y continuidad de los polígonos.
6. Si dos o más objetos se superponen en niveles diferentes, se asigna el nivel superior continuamente.
7. En caso de que dos o más objetos se superpongan en el mismo nivel de altura, el visualmente dominante y
el objeto completo (en uso y forma) se mapea continuamente.
Por definición y como premisa en la producción del Urban Atlas, existe la necesidad de alcanzar un valor
fiabilidad de 85-90%. Las muestras utilizadas para la evaluación deben ser correctas al menos en un 85-90%
de los casos. Este requerimiento y las siguientes tareas de verificación descritas se realizan además de para el
Urban Atlas para el Riparian Zones y Natura 2000.
Las pautas establecen las reglas para obtener resultados en la verificación y establecen los criterios para
considerar la calidad de un conjunto de datos ya verificados. Sin embargo, en el caso de que estas pautas no
fueran suficientemente precisas, puede provocar una fiabilidad de datos no óptima.
- Calidad de los datos. El objetivo de la evaluación. Si hay errores geométricos o semánticos, son medidos o
identificados por los procesos de verificación.
- Metodología de verificación. Si hay errores fuera del alcance de la metodología, se omiten en la revisión y
no se miden.
Compartimos el objetivo de la verificación del país para el componente local de Copernicus con el fin de
proporcionar una información complementaria basada en la experiencia regional / local.
El proceso de verificación se basa en una inspección visual de polígonos de muestra que analizan
características geométricas y semánticas.
- Corrección del código LC / LU alrededor del punto de muestra dentro del polígono de muestra, considerando
el valor real de MMU:
o Correcto / No correcto
o Si el código LC / LU se encuentra como incorrecto, se le pide al fotointerprete que proporcione su
opinión sobre el código correcto.
- Corrección de la delimitación del polígono de muestra:
o Área delineada: correctas / partes innecesarias incluidas / partes incompletas / partes incompletas y
partes innecesarias incluidas.
o Detalle de la delineación: Correcto / Demasiado grueso / Demasiado detallado
o Precisión posicional: correcta / desplazada
Sin embargo, su riqueza en contenido de información hace que el servicio sea adecuado para otras aplicaciones
como:
El Urban Atlas forma parte del Componente Local y desde la propia página web es posible descargarse los
datos de las capas en formato vectorial junto con Riparian Zones y Natura 2000.
https://land.copernicus.eu/local
Figura 7: Imagen del acceso a cada producto del componente local en la página web del Servicio de Territorio
de Copernicus.
Seleccionando Urban Atlas se puede acceder a la descarga de cada uno de los tipos de productos de Urban
Atlas disponibles. Como se ha comentado previamente, cada tipología puede tener más de un producto. Es
precisamente en esta navegación donde pueden encontrarse.
Figura 8: Ejemplo de navegación de Urban Atlas, donde se identifican 5 tipos de productos. Tras la selección
de uno de ellos, se accede a la descarga.
En función de la versión del producto es posible que estén desactivadas algunas de las acciones previamente
comentadas.
Accediendo a la opción de descarga, es posible listar cada una de las partes o conjuntos de datos descargables
ligados a los productos. En el caso concreto de Urban Atlas 2012 es posible descargarse los datos en un zip
por ciudad y por país (que contiene el Urban Atlas en formato shape y una guía de usuario).
Es muy importante conocer que los datos de Urban Atlas descargables desde la web oficial del Servicio de
Territorio de Copernicus se encuentran en el sistema de referencia y proyección europea ETRS89 LAEA
(EPSG:3035) y que es necesario tener un usuario en la web para su descarga.
Por otro lado, Copernicus tiene desarrollados servicios de red centralizados mediante el portal DiscoMap de la
Agencia Europea de Medio Ambiente. Se debe recordar que el acceso y conocimiento de servicios web de
productos Copernicus también se encuentran descritos en la plataforma CORDA de la Agencia Europea de
Medio Ambiente.
http://discomap.eea.europa.eu/
https://corda.eea.europa.eu/
El Instituto Geográfico Nacional, como Centro Nacional de Referencia de Cubiertas el Suelo de la Red
EINOET, es el encargado en la validación nacional del Urban Atlas. Adicionalmente el IGN ha publicado
mediante el Centro de Descargas del CNIG los datos del Urban Atlas 2012 comentados en este documento
junto con los otros dos productos del componente local del Servicio del Territorio de Copernicus que son
Riparian Zones 2012 y Natura 2K 2012.
Desde la pestaña “Información Geográfica Temática” es posible encontrar todos los productos de ocupación
del suelo en el que el IGN participa (coproduce o valida), entre ellos el Urban Atlas. En este nodo de descarga
es posible conseguir los datos de las Urban Atlas a 0,25 ha de resolución y en el sistema de referencia y
proyección nacional ETRS89 LAEA para toda España. La documentación anexa a los datos se encuentra
directamente contenidas en los ficheros de descarga.
Figura 11: Ventana de acceso a descarga de los datos y documentación sobre Urban Atlas en el Centro de
Descargas del CNIG
Adicionalmente el IGN / CNIG ha generado un servicio web de visualización para contener todos aquellos
productos validados en su cometido como Centro Nacionales de Referencia de Cubiertas el Suelo. Su dirección
de acceso es la siguiente.
http://servicios.idee.es/wms/copernicus-landservice-spain
Componente In Situ
2.1 Observaciones
El Programa Copernicus utiliza una amplia gama de infraestructuras de observación terrestre, cada una
creada para servir a una serie de propósitos diferentes y diversos, con el fin de producir, enriquecer y validar
sus servicios, mejorando su precisión y confiabilidad.
Los datos de observaciones, son información geoespacial ligada a la observación directa de características
físicas de la atmósfera y del mar. Suelen ser variables captadas por sensores terrestres (p.e. estaciones
meteorológicas, boyas, aeroplanos, globos, etc.). Esto incluye sensores colocados en las orillas de los ríos,
transportados en globos meteorológicos o aviones, arrastrados por el mar en barcos y flotando en el océano
en flotadores o boyas. También abarca datos recopilados por drones, por ciudadanos o científicos.
Algunos de los datos recogidos sobre el entorno físico se han recopilado durante siglos y, en consecuencia,
existe una gran variedad de registros y responsables interesados que participan en la recopilación y gestión
de dichos datos. Este trabajo se lleva a cabo en dominios temáticos, los cuales son transversales con respecto
a las necesidades de los Servicios de Copernicus.
- Meteorología: Datos históricos y diarios de estaciones terrestres y sensores atmosféricos en relación
con la temperatura, precipitación, humedad relativa, viento, índice solar, etc.
- Clima: Datos históricos climatológicos y series temporales armonizadas.
- Oceanografía. Datos físicos sobre la temperatura, salinidad, corrientes marinas, altura del nivel
medio del mar o concentración de hielos. Así como datos biológicos como la clorofila, oxígeno
disuelto en el agua, nutrientes, zooplancton, flujo de radiante o metano.
- Atmosfera: Datos sobre la concentración de polución en el aire (p.e. NOx, PM10), composición
química y concentración vertical de los aerosoles o gases de efecto invernadero (p.e. CO2).
- Terrestres: Datos físicos medidos mediante variables sobre el estado de la superficie terrestre, como
por ejemplo: FAPAR, LAI, biomasa sobre nivel del suelo, huella de carbono, humedad del suelo,
etc.
Muchas fuentes de datos in situ son propiedad y están operadas por instituciones y agencias a nivel nacional.
Un buen ejemplo son los datos nacionales de control de la calidad del aire recopilados por algunos Servicios
Meteorológicos nacionales, que los estados miembros de la UE están obligados a informar a la Agencia
Europea de Medio Ambiente.
Una creciente cantidad de datos in situ provienen de infraestructuras de investigación y otras colaboraciones
internacionales, como la Comisión Técnica Conjunta para Oceanografía y Meteorología Marina (JCOMM),
un organismo intergubernamental que reúne la experiencia y las capacidades tecnológicas de la Organización
Meteorológica Internacional (WMO) y la Comisión Oceanográfica Intergubernamental de la UNESCO
(COI).
Las redes de investigación científica a menudo también sirven como fuentes de datos de observaciones in
situ valiosos para el programa Copernicus. Gran número de estas redes existen en Europa, algunos de las
cuales están financiadas como infraestructuras de investigación europeas. En el área de Medio Ambiente, se
están coordinando esfuerzos y desarrollo en proyectos globales como ENVRIplus. Los ejemplos más
característicos de infraestructuras de investigación europeas en el ámbito de datos de observaciones in situ
son:
- Sistema integrado de observación de carbono (ICOS), una infraestructura de investigación
paneuropea que proporciona datos armonizados y de alta precisión sobre el ciclo del carbono, el
presupuesto de gases de efecto invernadero y sus perturbaciones.
Ilustración 1 Ejemplo de estaciones o sensores que recogen observaciones de datos in situ del mar y la
atmósfera
Es muy importante recalcar que el presente documento está centrado en el estudio de los datos espaciales In
Situ demandados por Copernicus.
Este servicio proporciona datos geográficos y productos basados en imágenes de satélite de muy alta
resolución (VHR) e información in situ, y respalda temas como la agricultura, la seguridad alimentaria, la
gestión forestal y la gestión del agua. Proporciona información general y específica, coherente y detallada
sobre la cobertura terrestre y el uso del suelo (LC / LU) en intervalos de tiempo regulares:
- CORINE Land Cover (CLC), conjunto de datos clave para Europa sobre cubiertas y usos del suelo
que proporciona una serie temporal y sus cambios desde 1990 hasta 2018. Comprende 44 clases con
una unidad geométrica estándar de 25 ha y 5 ha para cambios.
- Capas de Alta Resolución (HRL), actualizadas actualmente cada tres años, que proporcionan
características de cobertura terrestre paneuropeas consistentes derivadas de datos de observación de
la Tierra a una resolución espacial de 20 m. Las HRL actualmente disponibles de los años de
referencia 2012 y 2015 abordan 5 tipos principales de cobertura terrestre: impermeabilidad
(superficies selladas), bosques, pastizales naturales, humedales y cuerpos de agua permanentes.
- El conjunto de datos de Urban Atlas proporciona datos LC / LU coherentes y paneuropeos para áreas
urbanas funcionales (es decir, ciudades con más de 50,000 habitantes y sus alrededores).
- Las zonas ribereñas representan áreas de transición entre los ecosistemas terrestres y de agua dulce,
suministrando una amplia gama de funciones ribereñas y servicios ecosistémicos.
- Natura 2000 (N2K) representa una red comunitaria de áreas protegidas establecida en la Directiva de
Hábitats de 1992. El conjunto de datos de LC / LU de los sitios Natura 2000 seleccionados se
adaptan a las necesidades del seguimiento de la biodiversidad.
- El Componente Global ofrece una serie de productos biogeofísicos para el seguimiento del estado y
la evolución de las superficies terrestres a escala global. La producción y la entrega de los
parámetros se realizan de manera oportuna y se complementan con la generación de series de tiempo
a largo plazo. Los productos se utilizan para controlar la vegetación, el ciclo del agua y el
presupuesto energético.
Una amplia variedad de datos in situ con una cobertura paneuropea preferiblemente (al menos a nivel
nacional) son necesarios para generar productos y servicios dentro del componente CLMS. Los datos in situ
más importantes son, en orden de prioridad: ortoimágenes, mapas topográficos, datos del sistema de
información de parcelas (LPIS), datos de elevaciones, hidrografía, datos de cubierta y usos del suelo o
límites administrativos y datos adicionales de fuentes abiertas como OpenStreetMap, o privativas como
Bing, Google, etc. En detalle:
- Ortoimágenes: imágenes georreferenciadas de la superficie de la Tierra con una resolución espacial
muy alta, ya sea por satélite o por sensores aerotransportados, idealmente como imágenes de color
infrarrojo.
- Elevaciones: modelos digitales del terreno para superficie de tierra y hielo. Incluye elevación
terrestre y línea de costa.
- Cubierta y uso del suelo: cubierta real observada de la superficie de la tierra, tanto en términos de
características naturales como artificiales (p.e. tierras de cultivo, pastizales, bosques, vegetación
natural, suelo desnudo, humedales, área urbana). Los datos de los sistemas de información de
parcelas son también necesarios para una mejor identificación.
- Hidrografía: características naturales y artificiales del agua (p.e, arroyos, ríos, canales, lagos,
embalses) que representan los cuerpos de agua y la red de drenaje. La descarga del río, así como los
datos sobre el transporte de nutrientes y sedimentos, etc., son información in situ importante para la
cartografía de la zona ribereña y el futuro producto de la zona costera.
- Suelo: información sobre el tipo de suelo y subsuelo, que se caracteriza por la profundidad, textura,
estructura y contenido de partículas y material orgánico, pedregosidad, sensibilidad a la erosión y
capacidad de almacenamiento de agua.
Este servicio requiere los datos in situ actualizados sobre: asentamientos (tanto formales como informales),
infraestructuras de transporte, infraestructura portuaria, industria, servicios públicos, puntos de interés,
hidrografía, cubierta del suelo, topónimos, límites administrativos, fisiografía, densidad de población / datos
del censo, datos meteorológicos / oceanográficos, batimetría, datos de contaminantes,
Adicionalmente se requieren datos In Situ de sensores activos y pasivos (radar, electro-ópticos) montados en
activos marítimos (patrulleros) y aéreos (MPA, RPA), así como en estaciones terrestres de vigilancia.
Sistema de identificación automática por satélite (AIS), AIS en tierra, Servicio de tráfico de embarcaciones
(VTS), Sistema de seguimiento de embarcaciones (VMS), Identificación y seguimiento de largo alcance
(LRIT) e Inteligencia de señal (SIGINT).
La cobertura de datos in situ para este servicio, a diferencia de otros, debe ser mundial.
El Instituto Geográfico Nacional de España desempeña el rol de coordinador de este grupo de expertos
EIONET en datos In Situ asegurando satisfacer los objetivos planteados desde la EEA y motivando al
conocimiento de la importancia del acceso a los datos nacionales de observaciones y geoespaciales por parte
de los servicios Copernicus.
Como se puede observar, Copernicus demanda aquella información fundamental para conocer un territorio
que posteriormente sobre él se realizarán tomas de imágenes, o generación de productos derivados por parte
de sus servicios. Esta demanda coincide enormemente con la prioridades de IGR definidas en UN-GGIM y
en INSPIRE. Las autoridades nacionales, cuyo cometido es producir y mantener IGR oficial y confiable, y a
su vez publicar estos datos de manera interoperable como indica INSPIRE, han descubierto en los últimos
años que existe una relación directa entre las necesidades de información por la componente In Situ de
Copernicus y los datos IGR nacionales.
La situación deseada es que los países pongan en servicio su IGR según especificaciones interoperables
INSPIRE, y estos datos sean utilizables por Copernicus In Situ. Desgraciadamente esta situación no siempre
ocurre y se ve motivada por:
- No disponibilidad de datos IGR nacionales válidos. No siempre los países tienen datos espaciales
sobre las temáticas demandadas.
- Limitaciones de licencias, políticas de datos y condiciones de uso de la IGR de los países. Cada país
puede tener sus políticas propias referentes a la IGR.
- Limitaciones técnicas en las IGR de los países (p.e. escaso niveles de detalle geométrico o temático,
escala insuficiente, versiones antiguas de los datos, modelo de datos, deficiencias en la calidad de los
datos, etc.)
- Consistencia y comparabilidad de los datos IGR entre diversos países (p.e. los elementos geográficos
debe ser representados de manera homogénea entre todos los países).
- Y como es lógico voluntad de las agencias cartográficas (u organismos competentes
correspondientes) en poner en servicio su IGR para el componente In Situ de Copernicus.
En concreto en España, el IGN como agencia cartográfica nacional ofrece datos de referencia en la mayoría
de las temáticas de datos geoespaciales consideradas en la Componente In Situ. El IGN pone a disposición
de Copernicus datos de referencia nacionales a alto nivel de detalle: Información Geoespacial de Referencia
en Transportes, Poblaciones e Hidrografía, el Sistema de Información sobre Ocupación del Suelo en España
(SIOSE) que ofrece información de referencia para cubiertas y usos del suelo, así como el Plan Nacional de
Ortofotografía Aérea (PNOA) que ofrece información de referencia sobre elevaciones y ortoimágenes. IGN
trabaja de manera continua para que estos conjuntos de datos garanticen los requerimientos de Copernicus en
términos técnicos, de licencias y de acceso a datos, y puedan ser considerados como datos in situ.
A continuación se describen brevemente los diferentes proyectos del IGN en materia de IGR. Estos
proyectos se encuentran mucho más en detalle en otros apartados de la documentación.
- PNOA Imagen. Proporciona coberturas de imágenes aéreas ortorrectificadas cada 3 años para todo el
territorio nacional a altas resoluciones de 50 y 25 cm, a partir de especificaciones de captura y
procesado consensuadas por los usuarios en organismos públicos.
- PNOA Lidar: Proporciona coberturas de datos altimétricos cada 6 años para todo el territorio
nacional, a partir de especificaciones de captura y procesado consensuadas por los usuarios en
organismos públicos. Los productos generados del proyecto son nubes de puntos Lidar de al menos
0.5p/m2 de densidad, y productos derivados como modelos digitales del terreno (MDT, MDS) de 25,
5 y 2 m de resolución.
- IGR Hidrografía: Red fluvial precisa y actualizada. Proporciona información espacial de alta
resolución, homogénea a nivel nacional para satisfacer a los responsables en información
higrológica. Extraída de manera automáticamente a partir de una cobertura completa de Lidar para
todo el territorio español con una densidad de 0,5 puntos por metro cuadrado. Posteriormente se
generaron modelos de terreno hidrológicos con 2 m de resolución, de los cuales se extrajo la red
combinando criterios hidrográficos (red topográfica) y criterios hidrológicos (red de ríos de
acumulación de flujo). Los resultados de esta producción es una red de hidrografía automática
precisa extraída para todo el país con una mejora significativa para el componente altimétrico del
vector lineal 3D y coherente con los conjuntos de datos de elevación.
- IGR Transportes: Redes de transporta precisa y actualizada. Estas redes se generan por integración y
depuración de los actuales conjuntos de datos IGN, resto de agencias cartográficas españolas,
catastro, correos y agencias nacionales y regionales de transporte. Abarca transportes por carretera
(intra e interurbanos), ferrocarriles, marítimos y aéreos. Posee un modelo de datos y
especificaciones, conformes con INSPIRE, identificando características y atributos consensuados
entre IGN, agencias cartográficas regionales y proveedores de datos oficiales.
Los datos presentes en CORDA son recopilados por Coperncius, pero también se basa en una gran cantidad
de información recopilada en tierra, aire y mar, desde sensores in situ, así como datos de referencia y datos
auxiliares de organismos responsables a lo largo de toda Europa (p.e. Información Geoespacial de Referencia
de los estados miembros)
La EEA trabaja en estrecha colaboración con importantes redes de proveedores de datos europeos para
mejorar el acceso a datos in situ esenciales. Los conjuntos de datos adicionales se agregarán gradualmente a
medida que se establezcan nuevos acuerdos con los proveedores de datos. A nivel paneuropeo, la EEA se
centra actualmente en crear o revisar acuerdos con EuroGeoGraphics para datos de referencia espacial,