Landsat8 DN zu TOA Reflexion (Minuswerte)

8

Ich lerne, wie man die atmosphärische Korrektur manuell durchführt, daher verwende ich nicht die in ENVI verfügbaren automatisierten Methoden

Die Metadaten (Band 1) für meine Landsat 8-Szene lauten wie folgt:

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

Wenn ich also Bandmathematik verwende, um die TOA-Strahlung zu berechnen, gebe ich Folgendes ein:

(0,012344 * B1) -61,72016

Und es liefert eine Spektralstatistik, die irgendwie richtig aussieht SpectralStats

Wenn ich eine der folgenden Bandmathematikformeln verwende;

TOA-Reflexion (ohne Anpassung der Sonnenhöhe)

(B1 * 0,00002) -0,1

TOA-Reflexion (einschließlich Anpassung der Sonnenhöhe)

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

Oder berechnen Sie direkt vom Roh-DN zum TOA-Reflexionsvermögen (einschließlich Anpassung der Sonnenhöhe).

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

Am Ende habe ich so etwas:

TOA Reflexion

das sieht falsch aus.

Ich habe gegoogelt und entweder sind meine Suchbegriffe Müll oder niemand verwendet manuelle Bandmathematik, um dies zu tun.

[BEARBEITEN]

Mikkels Antwort unten mit der überarbeiteten Berechnung der Bandmathematik (ohne die Zwischenschritte) ergab die folgenden Statistiken:

Geben Sie hier die Bildbeschreibung ein

die etwas realistischer aussehen.

Sie fragen sich, ob ich jedem Band 13.4305 hinzufügen oder das negative Minimum ignorieren soll?

[END EDIT]

[EDIT # 2]

Mikkel war wieder genau richtig und schlug vor, dass die -13,4305-Werte auf die Bereiche "keine Daten" um das Bild herum isoliert wurden.

Geben Sie hier die Bildbeschreibung ein

Wie Sie sehen können, nur bis zum nördlichen Rand der Szene isoliert.

[END EDIT # 2]

Kyuss
quelle

Antworten:

7

Ein Fehler in Ihrem Ansatz. Sie müssen nicht durch DN gehen, um zu strahlen.
Sie können direkt zum DN zum Reflexionsgrad gehen.
Bleib einfach bei ((B1*0.00002)-0.1)/0.74457226676389733207607359928648.

Mikkel Lydholm Rasmussen
quelle
das hat ganz gut funktioniert! Ich nehme an, es wurde auch nie erklärt, dass TOA-Strahlung / TOA-Reflexionsvermögen excl. Sonnenaufgang waren keine notwendigen Schritte in der Korrektursequenz.
Kyuss
Ist es gemäß dem oben bearbeiteten Beitrag erforderlich, jedem Band 13.4305 hinzuzufügen, um sie auf min = 0 zu bringen, oder ist das unnötig?
Kyuss
@kyuss, versuche zu überprüfen, wo sich diese -0.134305 befinden. Ich vermute, dass es sich um die Bereiche ohne Daten am Bildrand handelt.
Mikkel Lydholm Rasmussen
Genial, ich habe nicht einmal etwas anderes überprüft als die südlichen Grenzen, an denen es keine negativen Werte gab. Vielen Dank für all Ihre Hilfe.
Kyuss
@kyuss jetzt ist der nächste Schritt in Ihrem Prozess "einfach" die atmosphärische Korrektur durchzuführen.
Mikkel Lydholm Rasmussen
2

Hier ist eine gute Lektüre über die Physik der atmosphärischen Korrektur und könnte bei dem Problem helfen ... Ich bin ein Typ, der einen Knopf drückt, ich fühle deinen Schmerz, zumindest kannst du einfach auf diesen schönen QUAC- oder FLAASH-Knopf klicken, wenn du draußen bist der Wissenschaft.

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

geodranisch
quelle
danke Steven. Ich habe - wie Sie gefolgert haben - während meines Studiums viel gelernt, aber dieses modis-Dokument kann durchaus dazu beitragen, meine Verständnislücken in Bezug auf den Strahlungstransfer zu schließen.
Kyuss