Protocolo v2

Collection 2 images

Landsat bands

Landsat bands

Primera pregunta:

¿Que productos mantener en Nor?

Tareas

  • Descarga Automática
  • Selección Escena de referencia
  • Actualización de MongoDB e integración con PostgreSQL
  • Instanciación
  • Geo (Recorte con el shape de la escena 202_34)
  • Rad (Aplicación de la correción para pasar a SR de 0 a 1)
  • Normalización
  • Productos

 

  • Selección Escena de referencia

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
  • Selección Escena de referencia

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

  • Selección Escena de referencia

Imagen de referencia:

  1. ¿Se aplica la nueva imagen a toda la serie histórica?
  2. ¿Se normaliza la nueva con la antigua escena de referencia (o viceversa)?
  3. ¿Se usa una escena de referencia por sensor y/o periodo?
  • Selección Escena de referencia

Nomenclatura:

  1. ¿Ori por SR?
  2. ¿Geo y Rad existen pero no se guardan?
  3. ¿Nor pasa a Nor2?
  4. ¿Se mantiene el antiguo formato?
  • Selección Escena de referencia

Formato:

  1. Valores: Float/Integer (ahora Float 0-1)
  2. COG/GTIFF

Productos

Productos

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:

  • Zonas por debajo de 5 m s.n.m.
  • Zonas con Σ Fmask >= 5
  • Zonas con NDVI alto todo el año
  • Zonas con altura de la vegetación elevada (no podemos ver debajo de los árboles)
  • Zonas con mucha pendiente
  • Zonas con sombra orográfica
  • Zonas de nubes y sombra de nubes (*)

Productos

Productos

Productos

P10 >= 0.4

Mean >= 0.4

Productos Cob Veg

Productos Slope

Productos (zona No Inundable)

Productos (DTM <= 5 m s.n.m.)

Productos (Fmask)

Productos (Fmask)

Productos (zona inundable)

Productos (Water Mask)

(rs < 0.13 & Fmask >= 5 | dtm <= 5) | (rs < 0.13 & (ndvip10 < 0.4 & slope < 10 & cobveg < 90)

  • ¿ Que sumamos a Fmask la máscara de Fmask o la máscara de agua resultante?
  • ¿Hacemos otro sumatorio con nuestras máscaras de agua?
  • ¿Merece la pena calcular esos parámetros por ciclos?

Otras cuestiones

  • ¿Actualización de máquina virtual?
  • Landsat 7 gapfill
  • Nuevos productos (Fenología, GPP, etc...)
  • Metadatos

No Data:               -9999

Tierra limpia:                0

Agua:                              1

Nubes:                           2

Sombra Nubes:            3

Sombra Orográfica:    4

Valores Propuestos

Protocolo v2

By Diego García Díaz

Protocolo v2

  • 64