Un buen ejemplo de aplicación de la teledetección es el seguimiento de la superficie inundada en zonas húmedas.

Dentro del programa 4ESO+Empresa de la Consejería de Educación de la Comunidad de Madrid, Irene (IES María Rodrigo) y Sergio (IES Antonio López García) han realizado un ejercicio para comprobar esta capacidad. La zona de estudio elegida ha sido Laguna de Gallocanta, en Aragón. Esta es una laguna salina endorreica a casi 1000 m sobre el nivel del mar declarada como reserva natural. Es un sistema muy interesante desde el punto de vista paisajístico, hidrogeológico y biológico, destacando sus poblaciones de aves acuáticas.

Esta laguna tiene un nivel de agua muy variable, en función de la cantidad y distribución de la precipitación a lo largo de los años y de otros factores más complejos. Para comprobar cómo puede estudiarse la superficie ocupada por el agua mediante teledetección hemos usado imágenes de los satélites Sentinel-2, cuya resolución espacial y espectral son adecuadas al objetivo del estudio.

La identificación de la superficie del agua la hemos hecho en base a la señal recibida en el infrarrojo medio (banda 11 de Sentinel-2, que detecta muy bien las aguas someras). Una identificación de píxeles de agua se puede encontrar en los productos L2a de Sentinel-2, dentro de la información de calidad asociada a la corrección atmosférica (banda  quality_scene_classification), pero hemos comprobado que esta información falla bastante en este caso, tanto por la alta reflectancia del fondo (limos salinos) como por la frecuente presencia de crestas de oleaje (al ser una laguna muy extensa y en zona de bastante viento). En la imagen presentamos un ejemplo de cómo se veía la laguna el pasado 17 de febrero en la banda 11 (1.6 micras), en RGB color natural, y en la clasificación realizada por Sen2Cor.

Hemos utilizado 2 imágenes por año desde 2016, eligiendo escenas sin nubes de agosto y febrero, que representan momentos de probable máximo y mínimo anual de agua.

Para procesarlas hemos usado SNAP. Y lo que hemos hecho es:

1. i) visualizamos la zona de la laguna y la recortamos mediante Raster > subset >, eligiendo en la configuración spatial subset y dejando la selección por defecto (la imagen visible),
2. ii) usamos Raster > Band Maths para crear una nueva banda con la condición Band 11 < X, donde el umbral X lo hemos ajustado visualmente para cada imagen (no había un valor único adecuado para todas las fechas).

Los píxeles por debajo del umbral se contabilizan con Analysis > Statistics. En estas estadísticas, como la imagen creada por Band Math es una imagen lógica, solo existen los puntos con valor 1, y por tanto el número de píxeles de agua es igual al número total de píxeles.

El resultado es el que se muestra en la gráfica. Lo interesante sería por supuesto discutir la evolución frente a los patrones de precipitación; pero eso queda para otra edición del 4ESO+Empresa.