Een GeoDataFrame kruisen met een polygoon
Door de gegevensset met grondgebruik te combineren met de gegevensset met wijken, kunnen we nu onderzoeken wat het grondgebruik is in een bepaalde wijk.
Daarvoor moeten we eerst de doorsnede bepalen van de gegevensset met grondgebruik met een gegeven wijk. Laten we opnieuw de wijk Muette als voorbeeld nemen.
De gegevenssets voor grondgebruik en wijken zijn al geïmporteerd als land_use en districts, en de wijk Muette is geëxtraheerd in de shapely-polygoon muette. Verder zijn GeoPandas en matplotlib geïmporteerd.
Deze oefening maakt deel uit van de cursus
Werken met georuimtelijke data in Python
Oefeninstructies
- Bereken de doorsnede van de
land_use-polygonen met de enkelemuette-polygoon. Noem het resultaatland_use_muette. - Maak snel een plot van deze doorsnede en geef
edgecolor='black'mee om de grenzen van de verschillende polygonen duidelijker te zien. - Print de eerste vijf rijen van
land_use_muette.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# 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.____)