J'apprends à effectuer manuellement la correction atmosphérique, donc je n'utilise pas les méthodes automatisées disponibles dans ENVI
les métadonnées (bande 1) de ma scène Landsat 8 sont les suivantes;
RADIANCE_MULT_BAND_1 = 1.2344E-02
RADIANCE_ADD_BAND_1 = -61.72016
REFLECTANCE_MULT_BAND_1 = 2.0000E-05
REFLECTANCE_ADD_BAND_1 = -0.100000
SUN_ELEVATION = 48.12237303
Donc, lorsque 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 semblent plutôt correctes
Lorsque j'utilise l'une des formules mathématiques de bande suivantes;
Réflexion TOA (hors ajustement de l'élévation du soleil)
(B1 * 0,00002) -0,1
Réflexion TOA (y compris l'ajustement de l'élévation du soleil)
((B1 * 0,00002) -0,1) /0,74457226676389733207607359928648
Ou calculez directement du DN brut à 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 me retrouve avec quelque chose comme ça:
qui semble incorrect.
J'ai googlé et soit mes termes de recherche sont des ordures, soit personne n'utilise le calcul manuel des bandes pour ce 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 EDIT]
[EDIT # 2]
Mikkel était de nouveau sur place, suggérant que les valeurs de -13,4305 étaient isolées dans les zones "sans données" autour de l'image;
comme vous pouvez le voir, isolé uniquement à la limite nord de la scène.
[FIN EDIT # 2]