Absolute Korrekturmodelle
Algorithmus 1:
Input:
- N Pixelwerte mit folgenden Informationen pro Pixel:
- Wellenlänge
- Korrektionsparameter der Atmosphäre
- lokale Information
- Tag der Aufnahme
Output: Reflexionsgrad der Oberfläche pro Pixel
Beschreibung:
Schritt: Lesen der Input-Parameter und Erstellen der Parameter-Tabelle
(LUT = Look-Up Table)
Schritt: Ausarbeiten einer Initialisierung und der benötigten Normalisierung
Schritt: Berechnung von mithilfe Interpolation der Tabellenwerte
Schritt: Berechnung der Oberflächenreflexion mithilfe der oben genannten Strahlungsbilanz
zu 2.:
In diesem Schritt wird die solare Stralungsstärke in Abhängigkeit der Wellenlänge mittels Interpolation berechnet (aus Strahlungsintensität und Entfernung Erde – Sonne). Ebenso kann man die Absorption von Kohlendioxid und Ozon ermitteln und somit die Sonnenstrahlung korrigieren.
Die verwendeten Interpolationen sind Spline-Interpolationen vom Grad 1. Vereinfacht dargestellt sind das stückweise lineare Funktionen.
Zu 3.:
Hierbei werden sowohl lineare als auch nicht-lineare Interpolationen angewandt. Zur Berechnung von L0 sind beispielweise folgende Schritte notwendig:
Interpolation von L0 für die gemessenen Höhe (linear)
Interpolation von L0 für die entsprechende Wellenlänge sowie für die Absorptionen von Gasen und Wasserdampf (stückweise exponentiell)
Interpolation von L0 für den Zenitwinkel der Sonne (Spline-Interpolation)
Interpolation von L0 für den Azimutwinkel des Beobachtungsortes (Spline-Interpolation)
Interpolation von L0 für den Zenitwinkel des Beobachtungsortes (Spline-Interpolation)
Interpolation von L0 für die optische Dicke (nicht-linear)
Der Algorithmus für die direkte Berechnung der atmosphärischen Korrektion benötigt also insgesamt 8 Look-Up Tables. Zwei der AVHRR Bänder (Advanced Very High Resolution Radiometer) und sechs der TM Bänder [siehe: Sensoren] müssen korrigiert werden.