Primera pregunta:
¿Que productos mantener en Nor?
Se han buscado imágenes Landsat 8 y 9 posteriores a la fecha de la primeras imágenes de Landsat 9 (2021-10-31) que estuvieran completamente libres de nubes.
Una vez seleccionadas las escenas se procedió a ver el histograma de las bandas 1 y 2 (Coastal blue y blue).
Las siguientes escenas cumplían los requisitos:
LC08_L2SP_202034_20220802_20220806_02_T1
LC08_L2SP_202034_20220818_20220824_02_T1
LC08_L2SP_202034_20230805_20230812_02_T1
LC09_L2SP_202034_20230712_20230714_02_T1
LC09_L2SP_202034_20230728_20230802_02_T1
Finalmente la decisión ha sido la escena Landsat 8 del 6 de agosto de 2022.
Se ha preferido una Landsat 8 porque parecen radiométricamente más coherentes que Landsat 9 y porque están disponibles desde 2013, habiendo por tanto muchas más escenas en ese periodo.
20220806l8oli_202_34
Imagen de referencia:
Nomenclatura:
Formato:
El código está hecho y funcionando, a falta de normalizar las escenas para probarlo con las imágenes normalizadas
El mayor esfuerzo de tiempo se ha empleado en la mejora de la máscara de máxima inundación. Se trata de un enfoque dinámico basado tanto en donde puede haber agua como en donde no puede haber agua:
P10 >= 0.4
Mean >= 0.4
(rs < 0.13 & Fmask >= 5 | dtm <= 5) | (rs < 0.13 & (ndvip10 < 0.4 & slope < 10 & cobveg < 90)