He estado tratando de calcular el NDVI para una escena usando Sentinel 2 bandas 4 (rojo) y 8 (NIR). Primero lo hice en ArcMap y luego en QGIS. En ambos softwares, utilicé la calculadora ráster e incluí la fórmula básica de NDVI = (NIR-Red) / (NIR + Red).
ArcMap me dio un rango de 0-9221 (esto varía con el estiramiento aplicado) mientras que QGIS me dio valores más plausibles de 0.0258-0.5255.
He leído brevemente que Sentinel 2 es una reflectancia TOA y tiene un factor de escala de 10000, que aparentemente se puede encontrar en el archivo xml que se descarga con las imágenes. Nunca encontré este factor de escala mencionado en el archivo xml que descargué.
¿Por qué hay esta diferencia? ¿Cómo se realiza el cálculo en QGIS?
Respuestas:
Arcmap es conocido por retener el tipo de ráster y no tratar los valores como decimales / flotantes cuando se calcula a partir de rásteres enteros de 16 bits.
Para forzar a Arcmap a tratar el resultado del cálculo de ráster como flotante, debe multiplicar uno de los argumentos por 1.0, o simplemente agregar un '.0' a cualquier número constante o real en la fórmula.
En su caso
(NIR-Red)*1.0/(NIR+Red)
será suficiente.fuente
¿Verificó los histogramas iniciales de las bandas respectivas que está utilizando? Es posible que cada software esté interpretando diferentes valores sin datos que pueden sesgar los resultados ...
fuente