Prozessierungsablauf

Die SAR Geophysical Retrieval Toolbox (SGRT)¹ erlaubt die Benützung von Sentinel-1A and Sentinel-1B Rohdaten (Level 1, SAFE Format) als Eingangsdaten, um höherwertige Produkte zu erzeugen, wie zum Beispiel vorprozessierte (Level 2), Parameter und Produkt-Datensätze (Level 3-4) (die Namenskonvention folgt der SGRT Terminologie). SGRT wird innerhalb der Mikrowellen-Fernerkundungs Forschungsgruppe am Geo Departement der TU Wien entwickelt.

Vorprozessierung

Die vorprozessierten Daten werden von den Sentinel-1 Rohdaten mithilfe von Sentinel Application Platform (SNAP) abgeleitet. Der Interferometric Wide (IW) Streifenmodus (cf. Aufnahme Modi) wird hauptsächlich für die Überwachung von Landgebieten benützt und ist daher der meistverwendetste Aufnahmemodus, welcher Daten für die Vorprozessierung der Sentinel-1 Daten an der TU Wien liefert.  Eine Roh-/.SAFE Datei (d.h. ein Verzeichnis) beinhaltet zwei Bilder, eines mit  VV und eines mit VH Polarisation, und Metadaten, welche für die weiteren Prozessierungschritte essentiell sind. Diese Bilder sind in der Streifen-/Sensorgeometrie verfügbar, das heißt das Koordinatensystem ist an den Sensor gebunden (in Flugrichtung (Azimut) und normal zur Flugrichtung (Distanz)), und nicht an die Erde (z.b. durch geographische Länge und Breite). Die folgenden Schritte sind nötig, um diese Rohdaten in eine zuverlässige Darstellung der Rückstreuung und anderer Produkte umzuwandeln:

• Aktualisierung des Positionsvektors in den Metadaten anhand von präzisen Bahndaten.
• Entfernung von thermischem Rauschen (tritt in Regionen mit niederer Rückstreuung auf und wird aufgrund von zufälligen Bewegungen geladener Teilchen verursacht).
• Eliminierung von Randartefakten, welche bei vorhergehenden Prozessierungsschritten erzeugt werden.
• Radiometrische Kalibrierung zur Konvertierung  von einem digitalen Pixelwert in eine kalibrierte Rückstreuung.
• Verwendung von Multi-looking, um Rauschen in SAR Bildern zu reduzieren (e.g. Speckle).
• Auf Anfrage kann eine radiometrische Geländeeinebnung erfolgen. Diese Prozedur korrigiert Layover-Effekte, welche wegen der seitlichen Aufnahmegeomtrie entstehen (ein digitales Geländemodell (DEM) ist ein verpflichtender Basisdatensatz für diese Korrektur).
• Eine Geländekorrektur adjustiert verzerrte Gebiete (aufgrund der seitlichen Aufnahmegeometrie) and reprojiziert die Daten von der Sensorgeometrie in eine georeferenzierte Darstellung/Projektion.
• Im letzten Schritt werden die vorprozessierten Daten in das EQUI7 Gittersystem reprojiziert und gekachelt (cf. Datenspezifikationen und Formate).

Nach dieser Routine und in Abhängigkeit von den gewählten Einstellung, sind folgende Level-2 Produkte verfügbar:

SIG0:

SIG0 entspricht der Rückstreuung von einer Einheitsfläche über Grund. Es ist hauptsächlich in VV und VH Polarisation verfügbar.

GAM0:

GAM0 repräsentiert die radiometrisch geländekorrigierte Rückstreuung von einer Einheitsfläche orthogonal zur Aufnahmerichtung. Die radiometrische Geländekorrektur berücksichtigt und trennt überlappende Pulse (in "Layover" und "Foreshortening" Regionen), welche aufgrung der seitlichen Aufnahmegeometrie entstehen. Dieses Produkt ist von einem minimalen Einfluss des Geländes auf die Rückstreuung gekennzeichnet mit dem Nachteil eines höheren Rechenaufwandes (cf ² für weitere Details). Es ist hauptsächlich in VV und VH Polarisation verfügbar.

PLIA:

PLIA ist die Abkürzung für Projected Local Incidence Angle und ist der Winkel zwischen der Oberflächennormalen und dem einfallenden Strahl  (auch bekannt unter lokaler Einfallswinkel (LIA)),  welcher in die Distanzebene projiziert wird. Die Oberflächennormale wird entweder über das Ellipsoid (vereinfachtes Modell der  Erde) oder über die lokale Nachbarschaft, welche mithilfe eines DEMs ermittelst werden kann, festgelegt.

Folglich stehen für eine Sentinel-1 Aufnahme (genauer, für einen Zeitpunkt) zwei SIG0 Bilder (VV und VH), zwei GAM0 Bilder (VV und VH) und ein PLIA Bild zur Verfügung.

Produkte

SGRT Produkte werden entweder von einem Bild (genauer, ein Zeitstempel) oder von aggregierten vorprozessierten Daten (über ein bestimmte Zeitspanne) abgeleitet.  Es muss bemerkt werden, dass das radiometrisch geländekorrigierte GAM0 eine relativ neue Darstellung der Rückstreuung ist und daher im Moment nur auf der Vorprozessierungsebene verwendet wird. Eine Liste von höherwertigen Produkten sind im Folgenden aufgelistet:

• Relative Bodenfeuchte (Surface Soil Moisture (SSM)): Der TU Wien Bodenfeuchte Algorithmus mit Sentinel-1 Daten kann unter anderem in ³ gefunden werden. SSM lieferte die relative Bodenfeuchte von den obersten Schichten des Bodens und stell damit einen unverzichtbaren Parameter für die Landwirtschaft, Dürre- und Vegetationsüberwachung und Klimaforschung dar.
• Unsicherheit der Bodenfeuchteschätzung (Surface Soil Moisture Noise (SSM-NOISE)): SSM-NOISE zeigt das "Rauschen", oder in anderen Worten, die Verlässlichkeit der Bodenfeuchteschätzung.
• Wasser (WATER): WATER besteht aus einer Kartierung von Wasserflächen.
• Nässe (Wetness (WWS)): Die Nässe ist ein kombiniertes Produkt, welches Informationen von SSM und WATER beinhaltet und eine Klassifikation von Wasserflächen, feuchten Gebieten und nicht feuchten Gebieten darstellt.
• Saisonale RGB Komposite (S-COMP): Die Komposite vermitteln bzw. visualisieren Änderungen der Landbedeckung basierend auf monatlichen Mitteln der Rückstreuung, die für die drei unterschiedlichen Kanäle zu unterschiedlichen Zeitpunkten eingesetzt werden. Dies beinhaltet auch eine Variation der Polarisation für jeden Kanal. Aufgrund ihrer schönen visuellen Impression zeigen diese Komposite die Möglichkeit der Detektion unterschiedlicher Landbedeckungstypen mithilfe des Sensors.

Parameter

Parameter-Produkte  beruhen entweder auf einer Zeitserie von vorprozessierten (SIG0, GAM0, PLIA) oder höherwertigen Produkten (WWS, SSM), und zielen auf die Aggregation der jeweilgen Produkte über die zeit ab. Die Zeitintervalle der Aggregation können folgendermaßen festgelegt werden:

• eine benutzerdefinierte Zeitspanne (zeitliche Kennzeichnung "T")
• saisonal (zeitliche Kennzeichnung "S")
• monatlich (zeitliche Kennzeichnung "M")
• täglich, z.b. über eine bestimmte Anzahl von Tagen (zeitliche Kennzeichnung "D")

Für eine detailierte Beschreibung der Namenskonvention dieser Produkte sei der Leser auf Dateinamenskonvention verwiesen.  Die anschließende Auflistung zeigt die meisten verfügbaren, statistischen Maße:

• Minimum ("MIN")
• Maximum ("MAX")
• Mittelwert ("MEN")
• Standardabweichung ("STD")
• Perzentile (5^th, 10^th, 25^th, 50^th (Median), 75^th, 90^th, 95^th) ("P05",  "P10", ...)
• Häufigkeit (Wahrscheinlichkeit des Auftretens einer Klasse in einer gewissen Zeitspanne) ("FRQ")

Beispiele

Bodenfeuchte (SSM)

Sentinel-1 SSM product at 500m spatial resolution over Austria. The depicted acquisition ranges from 02.08.2016 to 10.08.2016. One can identify dry regions in the eastern part and wet regions in the center part. Water bodies, mountains and cities are depicted as white, since a soil moisture retrieval can not be performed for this types of land cover (at the moment). Additionally, forest regions also do not yield a reliable estimate of soil moisture, but have not been masked out due to the lack of a forest mask.

Saisonale Komposite (S-COMP)

Left: Landsat-8 RGB image of London area, UK. Middle: False colour composite (RGB) of Sentinel-1 backscatter data over London area, UK. Red band shows mean of backscatter (VH) during summer (Jun-Jul-Aug), Blue band shows mean of backscatter (VH) during winter (Dec-Jan-Feb), and Green band shows the ratio of the Red and Blue band (mean_summer/mean_winter). Right: False colour composite (RGB) of Sentinel-1 backscatter data over London area, UK. Red band shows mean of backscatter (VV) during winter (Dec-Jan-Feb), Blue band shows mean of backscatter (VH) during winter (Dec-Jan-Feb), and Green band shows the ratio of the Red and Blue band (mean_summer/mean_winter). This composite highlights the variation of backscatter due to different polarizations. (data composites)

Referenzen

¹ Naeimi, Vahid et al. (2016). “Geophysical parameters retrieval from sentinel-1 SAR data: a case study for high performance computing at EODC”. In: Proceedings of the 24th High Performance Computing Symposium. Society for Computer Simulation International, p. 10.

² Small, David (2011). "Flattening gamma: Radiometric terrain correction for SAR imagery". In: IEEE Transactions on Geoscience and Remote Sensing 49.8, pp. 3081-3093.

³Gruber, Alexander et al. (2013). “Potential of Sentinel-1 for high-resolution soil moisture monitoring”. In: Geoscience and Remote Sensing Symposium (IGARSS), 2013 IEEE International. IEEE, pp. 4030–4033.