Umwelt-, Agrar-Reporting
Learnings
Verbindung mit der API herstellen
Download von Umweltdaten
Vergleich der Landbedeckung von 2 gleichgroßen Gebiete
Visualisierung der Ergebnisse
VORBEREITUNGEN
In diesem Beispiel wollen wir zwei gleichgroße Gebiete hinsichtlich ihrer Anteile an Grünfläche, Laubwald, Nadelwald und Fläche, die mit Gewässern bedeckt ist, analysieren und vergleichen.
Die deeeper.API nutzt den OGC-API Standard, um möglichst kompatibel mit gängigen GIS-Systemen zu sein.
Wir definieren zunächst alle Imports sowie die Basis-URLs von deeeper.Kassandra. Mit OWSLib steht ein Python-Paket zur Verfügung, durch das sich die dem OGC Standard entsprechenden APIs leicht ansprechen lassen. Der API Key ist nur zu Demozwecken und auf Abfragen des Gebiets der Städteregion Aachen beschränkt.
Zusätzlich werden die Bounding-Boxen der Gebiete, die miteinander verglichen werden sollen, festgelegt.
Im folgenden Snippet definieren wir außerdem eine Funktion, durch welche die Daten der verschiedenen Klassen einfach für eine gegebe Bounding-Box abgefragt werden können. Die Funktion erleichtert später den Download der Daten und lässt sich für jede beliebige Klasse von deeeper.Kassandra anwenden.
BERECHNUNG DER FLÄCHE JE KLASSE
Nach Abschluss der Vorbereitungen sind wir nun bereit für die Analyse. Der nachfolgende Code-Ausschnitt demonstriert dafür beispielhaft wie man für 4 verschiedene Landbedeckungsklassen aus deeeper.Kassandra die Fläche innerhalb der angegebenen Bounding-Boxen berechnet. Da in den untersuchten Gebieten nur bestimmte Klassen vorkommen, konzentriert sich dieses Beispiel auf die Klassen Grünland, Laubwald, Nadelwald und Gewässer.
In jeder Zeile wird dabei der gleiche Prozess angewendet. Die Daten werden in einem ersten Schritt durch die vorher definierte Funktion heruntergeladen und anschließend wird mittels "set_crs" ein Koordinatenreferenzsystem für den geladenen GeoDataFrame festgelegt. Danach wird der GeoDataFrame auf das Koordinatenreferenzsystem EPSG:25832 projiziert. Dies ist notwendig, damit die Funktion "area", welche die Fläche jedes im GeoDataFrame vorhandenen Objekts berechnet, auch tatsächlich in m² ausgegeben wird und nicht durch eine falsche Projektion verzerrt wird.
Die berechneten Einzelflächen werden danach für jede Landbedeckungsklasse aufsummiert und in km² umgerechnet. Als Ergebnis kennen wir nun für alle vorkommenden Landbedeckungsklassen die Fläche innerhalb der definierten Gebiete.
VISUALISIERUNG
Die berechneten Ergebnisse können anschließend nach Wunsch visualisiert werden. Das hier dargestellte Diagramm stellt beispielsweise die Fläche der untersuchten Landbedeckungsklassen von Gebiet A (links) den Ergebnissen von Gebiet B (rechts) gegenüber. Obwohl der Anteil an Nadelwald in den gleichgroßen Gebieten annähernd gleichgroß zu sein scheint, gibt es in Gebiet A deutlich weniger Laubwald - dafür jedoch zusätzlich Anteile an Grünland und Gewässern.
