Einen GeoDataFrame mit einem Polygon schneiden
Wenn wir den Datensatz zur Flächennutzung mit dem Datensatz der Bezirke kombinieren, können wir untersuchen, wie die Flächennutzung in einem bestimmten Bezirk aussieht.
Dazu müssen wir zuerst die Schnittmenge des Flächennutzungs-Datensatzes mit einem gegebenen Bezirk bestimmen. Nehmen wir erneut den Bezirk Muette als Beispielsfall.
Die Datensätze zur Flächennutzung und zu den Bezirken wurden bereits als land_use bzw. districts importiert, und der Bezirk Muette wurde als shapely-Polygon in muette extrahiert. Außerdem sind GeoPandas und matplotlib importiert.
Diese Übung ist Teil des Kurses
<Kurs>Arbeiten mit Geodaten in Python</Kurs>Übungsanweisungen
- Berechne die Schnittmenge der
land_use-Polygone mit dem einzelnenmuette-Polygon. Nenne das Ergebnisland_use_muette. - Erstelle eine schnelle Visualisierung dieser Schnittmenge und übergib
edgecolor='black', um die Grenzen der verschiedenen Polygone klarer zu sehen. - Gib die ersten fünf Zeilen von
land_use_muetteaus.
Interaktive praktische Übung
Versuche dich an dieser Übung, indem du diesen Beispielcode vervollständigst.
# Print the land use datset and Notre-Dame district polygon
print(land_use.head())
print(type(muette))
# Calculate the intersection of the land use polygons with Notre Dame
land_use_muette = ____
# Plot the intersection
land_use_muette.____
plt.show()
# Print the first five rows of the intersection
print(land_use_muette.____)