En entradas anteriores hablamos del plugin SEN-ET de SNAP y de su complejo funcionamiento. Hoy, os presentamos el resultado de aplicarlo a una zona agrícola para mostrar el tipo de información que se produce. Hemos utilizado para ello la famosa zona de Cal/Val situada en las cercanías de Barrax (Albacete, España).

En la imagen podéis ver la salida final (ET en mm/día) junto a una composición RGB en infrarrojo color de la zona, que facilita su interpretación.

Un examen visual muestra el resultado esperado: valores muy altos (tonalidades claras) sobre zonas de cultivo y bastante más bajos (tonalidades oscuras) sobre suelos desnudos. Hay que advertir que la herramienta está pensada sobre todo para gestión del riego y estimación de balances hídricos en vegetación y, por tanto, los resultados en zonas sin vegetación son menos relevantes.

En el análisis visual llama la atención un artefacto rectangular (señalado con la flecha amarilla), que se repite en otras partes de la escena no mostradas en la figura. Este se corresponde a una clasificación incorrecta en el mapa de usos del suelo ESA CCI land-cover, que es el recomendado en el manual de SEN-ET. Para ese mapa, la zona de estudio alterna entre "cultivos herbáceos" y "cultivos agrícolas", pero en algunas celdas aparece la clase "Tree / shrub". La estimación de los parámetros estructurales de la vegetación para estas celdas resulta muy diferente a la del resto de la escena (e incorrecta en este caso), por lo que el valor final de ET estimado difiere significativamente de su entorno, causando el artefacto. Afortunadamente, SEN-ET permite introducir un mapa de usos del suelo diferente al de ESA CCI, y esta sería la solución para este problema.

También se aprecia un cierto patrón de ruido, aunque de pequeña amplitud. Este patrón aparece en el sharpening del producto LST de Sentinel-3, paso que los propios autores del plugin consideran la mayor causa de incertidumbre del resultado.

Para mostrar los resultados cuantitativos obtenidos, hemos preparado un gráfico bivariante con los valores de ET frente al índice de área foliar (LAI) estimado por el plugin. Se aprecia el esperado aumento de ET con LAI, pero organizado según diferente categorías, que se corresponden probablemente a clases de parámetros estructurales de vegetación.

En la publicación "Modelling High-Resolution Actual Evapotranspiration through Sentinel-2 and Sentinel-3 Data Fusion"  se discuten detalladamente la exactitud y las limitaciones del plugin. Os recomendamos estudiarla a fondo si, después de la breve presentación que os hemos ofrecido en las últimas entradas, tenéis interés en esta herramienta tan compleja y completa.