Informationen aus einer Rasterebene extrahieren
Lass uns nun Informationen aus der Rasterebene auf Basis einer Vektordatei extrahieren. Diese Funktionalität wird vom Paket rasterstats bereitgestellt. In dieser Übung bestimmen wir den Vegetationstyp an allen Abbaustandorten, indem wir für jeden Punkt im Datensatz der Abbaustandorte den nächstgelegenen Raster-Pixelwert abrufen.
Ein Teil des Datensatzes der Abbaustandorte (mining_sites) ist bereits geladen, und GeoPandas sowie matplotlib sind bereits importiert.
Diese Übung ist Teil des Kurses
Arbeiten mit Geodaten in Python
Anleitung zur Übung
- Importiere das Paket
rasterstats. - Extrahiere den nächstgelegenen Rasterwert für die Standorte der Abbaustandorte (
mining_sites) mit der Funktionrasterstats.point_query()und weise das Ergebnis einer neuen Spalte'vegetation'zu. - Erstelle einen Plot der Daten der Abbaustandorte und verwende die Spalte
'vegetation', um die Punkte einzufärben. Achte darauf, eine Legende hinzuzufügen.
Interaktive Übung
Vervollständige den Beispielcode, um diese Übung erfolgreich abzuschließen.
# Import the rasterstats package
____
# Extract the nearest value in the raster for all mining sites
vegetation_raster = "central_africa_vegetation_map_foraf.tif"
mining_sites['vegetation'] = ____(____, ____, interpolate='nearest')
print(mining_sites.head())
# Replace numeric vegation types codes with description
mining_sites['vegetation'] = mining_sites['vegetation'].replace(vegetation_types)
# Make a plot indicating the vegetation type
____
plt.show()