Suite

Réflectance Landsat8 DN à TOA (valeurs négatives)


J'apprends à effectuer la correction atmosphérique manuellement, donc sans utiliser les méthodes automatisées disponibles dans ENVI

les métadonnées (bande 1) pour ma scène Landsat 8 sont les suivantes ;

RADIANCE_MULT_BANDE_1 = 1.2344E-02

RADIANCE_ADD_BAND_1 = -61.72016

REFLECTANCE_MULT_BAND_1 = 2.0000E-05

REFLECTANCE_ADD_BAND_1 = -0,100000

SOLEIL_ELEVATION = 48.12237303

Donc, quand j'utilise les mathématiques de bande pour calculer TOA Radiance, je tape ceci:

(0.012344*B1)-61.72016

Et cela donne des statistiques spectrales qui ont l'air plutôt correctes

Lorsque j'utilise l'une des formules mathématiques de bande suivantes ;

Réflectance TOA (hors réglage de l'élévation du soleil)

(B1*0.00002)-0.1

Réflectance TOA (y compris le réglage de l'élévation du soleil)

((B1*0.00002)-0,1)/0,74457226676389733207607359928648

Ou calculez directement à partir du DN brut vers la réflectance TOA (y compris l'ajustement de l'élévation du soleil)

((((0,012344*B1)-61,72016)*0,00002)-0,1)/0,74457226676389733207607359928648

Je finis avec quelque chose comme ça :

ce qui semble incorrect.

J'ai cherché sur Google et soit mes termes de recherche sont nuls, soit personne n'utilise les mathématiques manuelles pour le faire.

[ÉDITER]

La réponse de Mikkel ci-dessous avec le calcul mathématique de bande révisé (en omettant les étapes intermédiaires) a donné les statistiques suivantes :

qui semblent un peu plus réalistes.

Vous vous demandez si je dois ajouter 13.4305 à chaque bande ou ignorer le minimum négatif ?

[FIN DE MODIFIER]

[EDIT #2]

Mikkel était encore une fois sur place, en suggérant que les valeurs -13,4305 étaient isolées des zones "pas de données" autour de l'image;

comme vous pouvez le voir, isolé seulement au bord nord de la scène.

[FIN EDIT #2]


Un défaut dans votre approche. Vous n'avez pas besoin de passer par DN pour rayonner.
Vous pouvez aller directement au DN pour la réflectance.
Il suffit de s'en tenir à((B1*0.00002)-0,1)/0,74457226676389733207607359928648.


Voici une bonne lecture sur la physique de la correction atmosphérique, et pourrait aider avec le problème… Je suis un gars qui pousse des boutons, je ressens votre douleur, au moins vous pouvez simplement cliquer sur ce magnifique bouton QUAC ou FLAASH une fois que vous êtes en dehors du milieu universitaire .

http://modis.gsfc.nasa.gov/data/atbd/atbd_mod08.pdf


Voir la vidéo: How to convert DN value to Reflectance using ArcGIS (Octobre 2021).