7 votos

¿Calcular el EVI a partir de las imágenes de Sentinel 2 da resultados diferentes entre ArcMap y QGIS?

Aunque los coeficientes EVI no se han desarrollado para las imágenes de Sentinel 2, estoy utilizando los mismos coeficientes utilizados para MODIS para calcular el EVI de las imágenes de Sentinel.

Sé que las imágenes de Sentinel 2 son de reflectancia TOA y tienen un factor de escala de 10000. Por lo tanto, puedo utilizar la siguiente fórmula para calcular el EVI

(2.5* (NIR/10000-red/10000))/(NIR/10000+(6*red/10000)- (7.5*blue/10000)+1)

Hice el cálculo anterior en Arcmap y QGIS, pero obtuve resultados completamente diferentes. QGIS me dio un rango plausible de -0,048 a 0,3915, mientras que ArcMap incluyó un rango de -31,68 a 45,17. La media y la norma en ArcMap son 0,1718 y 0,0887, lo que me lleva a pensar que estos valores extremos son sólo unos pocos píxeles. ¿Cómo puedo recortar estos píxeles adicionales en ArcMap?

0 votos

¿Ha comparado los píxeles reales de las imágenes, o sólo ha mirado las estadísticas de las imágenes proporcionadas?

0 votos

@MikkelLydholmRasmussen- He comparado las estadísticas de la imagen, pero también he comparado los valores de los píxeles (utilizando la herramienta de identificación en ambos programas) y están ligeramente desviados. Como mencioné en el post, estoy seguro de que los valores extremos obtenidos usando ArcMap son unos pocos píxeles que necesito recortar muy probablemente. ¿Puede sugerir un método para hacerlo?

1 votos

@MikkelLydholmRasmussen- Además, ¿tienes idea de si puedo usar el factor de escala como el anterior? ¿Dividir cada valor de reflectancia de las bandas por 10000 o se supone que el factor de escala se utiliza hacia el final?

6voto

Por favor, limítese a una pregunta por mensaje en el futuro. Así será mucho más fácil responder.

Pregunta 1 - Cómo calcular el EVI:
La fórmula que está utilizando es "correcta", en el sentido de que está utilizando las constantes calculadas para MODIS.

Pregunta 2 - ¿Por qué QGis y ArcGIS proporcionan resultados diferentes?
Las minúsculas diferencias en los resultados se deben probablemente a los diferentes métodos de remuestreo utilizados al combinar las bandas espectrales. Es difícil descubrir la razón de las diferencias, sin disponer de los datos.

Pregunta 3 - Cómo eliminar los valores extremos de un cálculo raster:
La mejor manera de eliminar esos valores es ajustando el cálculo de la trama. Un enfoque sería simplemente añadir un paréntesis alrededor de toda la ecuación y multiplicarla por un par de cálculos "mayor que" y "menor que".
Ejemplo:

(EVI) * (EVI > -1) * (EVI < 1) + (EVI > 1)*1 + (EVI < -1)*(-1)

El cálculo descrito anteriormente establecería todos los valores siguientes -1 y superiores 1 a -1 y 1 respectivamente.

Pregunta 4 - Cuándo aplicar los factores de escala:
Los factores de escala deben aplicarse cuando se utilizan las bandas. Si los aplicas después, en el EVI, el efecto de las constantes en la ecuación cambiaría, y tus cálculos serían incorrectos. Por lo tanto, el uso que haces es el correcto.

0 votos

Gracias por su respuesta. Y disculpas por las múltiples preguntas en un solo post. Lo publicaré por separado en el futuro.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X