COMO EMPLEAR IMÁGENES SENTINEL
Parte 2
A continuación haremos índices NDVI para las cuatro imágenes descargados,
empezamos con la zona RT_T29SQC de verano:
La fórmula de cálculo es la siguiente:
("RT_T29SQC_A032910_20230625T110839_B08" -
"RT_T29SQC_A032910_20230625T110839_B04") /
("RT_T29SQC_A032910_20230625T110839_B08" +
"RT_T29SQC_A032910_20230625T110839_B04")
Para la imagen de invierno es:
("RT_T29SQC_A030908_20230205T111142_B08" -
"RT_T29SQC_A030908_20230205T111142_B04")/
("RT_T29SQC_A030908_20230205T111142_B08" +
"RT_T29SQC_A030908_20230205T111142_B04")
Para la zona zona RT_T29SQC los cálculos son los siguientes:
Verano
("RT_T29SQD_A032910_20230625T110839_B08" -
"RT_T29SQD_A032910_20230625T110839_B04")/
("RT_T29SQD_A032910_20230625T110839_B08" +
"RT_T29SQD_A032910_20230625T110839_B04")
Invierno
("RT_T29SQD_A030908_20230205T111142_B08" -
"RT_T29SQD_A030908_20230205T111142_B04") /
("RT_T29SQD_A030908_20230205T111142_B08" +
"RT_T29SQD_A030908_20230205T111142_B04")
"NDVI_RTS9SQC_VERANO" - "NDVI_RTS9SQC_INVIERNO"
Para dar simbología:
https://www.youtube.com/watch?v=HgvyBtuGj3k
Repetimos la misma simbología para las cuatro imágenes NDVI
También se pueden calcular los índices NDVI con la calculadora ráster, para el
caso de la imagen RT_T29SQC de verano sería:
("RT_T29SQC_A032910_20230625T110839_B08@1" -
"RT_T29SQC_A032910_20230625T110839_B04@1" /)
(S"RT_T29SQC_A032910_20230625T110839_B08@1" +
"RT_T29SQC_A032910_20230625T110839_B04@1")
Y para la de invierno:
("RT_T29SQC_A030908_20230205T111142_B08@1" -
"RT_T29SQC_A030908_20230205T111142_B04@1") /
("RT_T29SQC_A030908_20230205T111142_B08@1" +
"RT_T29SQC_A030908_20230205T111142_B04@1")
Para calcular el Índice de Vegetación Estandarizado (SVI), con solo dos fechas,
se puede hacer de dos maneras:
Empleando la calculadora ráster, siendo las fechas de verano e invierno
correspondientes a la zona RTS9SQC (1):
("NDVI_RTS9SQC_VERANO@1"-("NDVI_RTS9SQC_VERANO@1"
+"NDVI_RTS9SQC_INVIERNO@1")/2)/sqrt((("NDVI_RTS9SQC_VERANO@1"-
("NDVI_RTS9SQC_VERANO@1" +"NDVI_RTS9SQC_INVIERNO@1")/2)^2)+
(("NDVI_RTS9SQC_INVIERNO@1"-("NDVI_RTS9SQC_VERANO@1"
+"NDVI_RTS9SQC_INVIERNO@1")/2)^2))
En primer lugar, se calcula como un ráster virtual, para que una vez hecho se
exporta el mismo como geotiff:
Me sale una imagen con los valores no coincidentes con la capa virtual .
Probaré la segunda opción, usando el plugin SCP:
El algoritmo será:
("NDVI_RTS9SQC_VERANO"-("NDVI_RTS9SQC_VERANO"
+"NDVI_RTS9SQC_INVIERNO")/2)/sqrt((("NDVI_RTS9SQC_VERANO"-
("NDVI_RTS9SQC_VERANO" +"NDVI_RTS9SQC_INVIERNO")/2)^2)+
(("NDVI_RTS9SQC_INVIERNO"-("NDVI_RTS9SQC_VERANO"
+"NDVI_RTS9SQC_INVIERNO")/2)^2))
Pero me sale la expresión en rojo, por lo que vamos a realizar dos imágenes
numerador y denominador de la expresión y en un segundo paso hacemos el
cociente de las mismas.
Para el numerador, será
Al hacer clic en RUN, sale el desplegable donde indico el nombre y ruta del
fichero, en este caso es numerador. La expresión es:
("NDVI_RTS9SQC_VERANO"-("NDVI_RTS9SQC_VERANO"
+"NDVI_RTS9SQC_INVIERNO")/2)
Con el denominador repito los pasos.
La expresión es:
sqrt((("NDVI_RTS9SQC_VERANO"-("NDVI_RTS9SQC_VERANO"
+"NDVI_RTS9SQC_INVIERNO")/2)^2)+ (("NDVI_RTS9SQC_INVIERNO"-
("NDVI_RTS9SQC_VERANO" +"NDVI_RTS9SQC_INVIERNO")/2)^2))
La operación final será:
"numerador"/"denominador"
El fichero resultado es:
RT_T29SQC DIFERENCIA ESTANDARD SAC.tif
Este cálculo del SVI solo nos da pixeles con dos valores que son -0,07606575 o
su simétrico 0,07606575.
Es un cálculo de prueba, una metodología más correcta sería la indicada en la
publicación de la sociedad española de teledetección: Anomalías de vegetación
asociadas con el fenómeno del ENOS en el valle geográfico del río Cauca,
Colombia (7715-33030-1-PB.pdf):
“El SVI es un índice basado en la probabilidad de verdor de un píxel según su
media histórica y que puede estar altamente influenciado por los efectos del
clima. Para la obtención del SVI se utilizaron las imágenes de EVI del producto
para el periodo comprendido entre los años 2000 y 2015.
Con los datos depurados y filtrados, se procedió a realizar un promedio y
desviación estándar global
y un promedio trimestral para obtener el SVI, (Ecuación 1) (Sousa et al., 2010).”
Otro tipo de cálculo que podría ser de interés es la anomalía del NDVI, que es la
diferencia entre el NDVI medio para un mes en particular de cierto año y el NDVI
promedio para el mismo mes a lo largo de un período específico de años.
(1) C:\Users\cesar.diaz\OneDrive - Junta de
Extremadura\Documentos\ESTUDIOS\SENTINEL\DOCUMENTACION\77
15-33030-1-PB.pdf
(2) 2011_Casagranda-Ma-Elvira.pdf
(3) Operational_drought_monitoring_in_Kenya_using_MODI.pdf
(4) https://www.sciencedirect.com/topics/earth-and-planetary-sciences/
normalized-difference-vegetation-index
(5) https://gis.stackexchange.com/questions/457911/calculate-a-
guassian-standardized-variable-in-google-earth-engine
(6) https://un-spider.org/advisory-support/recommended-practices/
recommended-practice-agricultural-drought-monitoring-svi/in-detail
(7) https://www.walshmedicalmedia.com/open-access/investigate-the-
sensitivity-of-the-satellitebased-agricultural-drought-indices-to-monitor-
the-drought-condition-of-paddy-and-intro-53154.html
Para hacer un mapa de porcentaje de cambio del indice NDVI entre dos fechas
se emplea an la calcukadora de bandas del plugin SCP el siguiente algoritmo:
where ("NDVI_RTS9SQC_INVIERNO" < 0, ((("NDVI_RTS9SQC_VERANO" -
"NDVI_RTS9SQC_INVIERNO") / "NDVI_RTS9SQC_INVIERNO") * 100),
((("NDVI_RTS9SQC_VERANO" - "NDVI_RTS9SQC_INVIERNO") /
"NDVI_RTS9SQC_INVIERNO") * 100))
Para estudir comparaciones entre dos NDVI:
Dónde se aprecia que la primera condicion del where es para cambiar de signo
a los valores negativos
de la imagen de invierno.
Lo represento según el criterio
siguiente:
Siendo este un extracto del resultado:
