Eu tenho tentado calcular o NDVI para uma cena usando as bandas Sentinel 2 4 (vermelho) e 8 (NIR). Eu o fiz primeiro no ArcMap e depois no QGIS. Nos dois softwares, usei a calculadora raster e incluí a fórmula básica NDVI = (NIR-Red) / (NIR + Red).
O ArcMap me forneceu uma faixa de 0-9221 (isso varia com o alongamento aplicado), enquanto o QGIS me forneceu valores mais plausíveis de 0,0258 a 0,5255.
Li brevemente que o Sentinel 2 é uma refletância TOA e possui um fator de escala de 10000, que aparentemente pode ser encontrado no arquivo xml que é baixado com as imagens. Nunca encontrei esse fator de escala mencionado no arquivo xml que baixei.
Por que existe essa diferença? Como o cálculo está sendo feito no QGIS?
Respostas:
O Arcmap é notório por manter o tipo de varredura e não tratar os valores como decimais / flutuação quando calculado a partir de rasters inteiros de 16 bits.
Para forçar o Arcmap a tratar o resultado do cálculo de varredura como uma flutuação, você deve multiplicar um dos argumentos por 1,0 ou simplesmente anexar um '.0' a qualquer número constante ou real da fórmula.
No seu caso
(NIR-Red)*1.0/(NIR+Red)
será suficiente.fonte
Você checou os histogramas iniciais das respectivas bandas que você está usando? É possível que cada software esteja interpretando diferentes valores sem dados que podem distorcer os resultados.
fonte