Suite

Valeurs Landsat8 NDVI non comprises entre -1 et 1

Valeurs Landsat8 NDVI non comprises entre -1 et 1


J'essayais de créer un NDVI à partir d'une image Landsat8 OLI de juillet 2013. J'ai exécuté DOS1 et TOA. Les valeurs vont jusqu'à 3. Les fichiers originaux ont été téléchargés depuis EarthExplorer (EE), exécutés sur QGIS et Multispec. Les données d'origine ont été téléchargées à nouveau à partir d'EE, puis exécutées uniquement dans QGIS (calculateur de bande).

Qu'est-ce que je fais mal?


Je suis passé à une image Landsat 5 et j'ai un problème. Après le TOA et le DOS1, voici les résultats pour le rouge visible : 'RT_LT50010862008320COA00_B3.sta' Nombre de lignes incorrectes = 0 Nombre total de pixels = 2 969 703 Tableau récapitulatif de l'histogramme Première dernière ligne Première colonne Intervalle de ligne Intervalle de colonne Intervalle de colonne 1 7251 20 1 8181 1 Plage de données de canal Moyenne Médiane Norme Description Écart 1 0,0075 à 0,6560 0,0409 0,0432 0,0294


Les valeurs NDVI en dehors de la plage valide -1 : 1 sont généralement causées par des valeurs de pixels négatives dans la bande rouge. Une fois votre pré-traitement terminé, vous devez toujours corriger les valeurs de pixels négatives avant de calculer le NDVI.

Un exemple:

Le NDVI est calculé comme

Supposons maintenantNIR = 100etVIS = -50. Votre calcul devient maintenant

NDVI = (100 - (-50)) / (100 + (-50)) = 3

Certains pixels sont saturés. C'est un problème L8 connu. Vous pouvez examiner les bandes originales - les valeurs de ces pixels doivent être inhabituelles.

Pour résoudre ce problème, vous pouvez définir des règles pour gérer ces pixels.

Par exemple:

si valeur > 1 : valeur = 1 si valeur < -1 : valeur = -1 si valeur < -1 OU valeur > 1 : valeur = 0