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


8

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 SpectralStats

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:

Réflexion TOA

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:

entrez la description de l'image ici

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;

entrez la description de l'image ici

comme vous pouvez le voir, isolé uniquement à la limite nord de la scène.

[FIN EDIT # 2]

Réponses:


7

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.
Tenez-vous-en ((B1*0.00002)-0.1)/0.74457226676389733207607359928648.


cela a plutôt bien fonctionné! je suppose que soit il n'a jamais été expliqué que la radiance TOA / réflectance TOA excl. l'élévation du soleil n'était pas une étape nécessaire dans la séquence de correction.
kyuss

d'après l'article édité ci-dessus, est-il nécessaire d'ajouter 13.4305 à chaque bande pour les amener à min = 0, ou est-ce inutile?
kyuss

@kyuss, essayez de vérifier où se trouvent ces -0.134305. Je soupçonne que ce sont les zones sans données autour du bord de l'image.
Mikkel Lydholm Rasmussen du

génial, je n'ai même pas vérifié autre chose que les frontières sud où il n'y avait pas de valeurs négatives. Merci pour toute votre aide.
kyuss

@kyuss maintenant la prochaine étape de votre processus est "simplement" de faire la correction atmosphérique.
Mikkel Lydholm Rasmussen du

2

Voici une bonne lecture sur la physique de la correction atmosphérique, et pourrait aider avec le problème ... Je suis un gars poussant le bouton, je ressens votre douleur, au moins vous pouvez simplement cliquer sur ce beau bouton QUAC ou FLAASH une fois que vous êtes à l'extérieur du monde universitaire.

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


merci Steven. comme vous l'avez déduit, j'ai beaucoup appris au cours de mon étude, mais ce document modis pourrait bien aider à combler mes lacunes en matière de transfert radiatif.
kyuss
En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.