Intersectie van twee polygonen

In deze oefening gebruiken we 2 afzonderlijke polygonen: het district Muette uit de gegevensset districts, en het groene stedelijke gebied Boulogne, een groot openbaar park in het westen van Parijs, uit de gegevensset land_use. De twee polygonen zijn al toegewezen aan de variabelen muette en park_boulogne.

We visualiseren eerst de twee polygonen. Je zult zien dat ze overlappen, maar het park ligt niet volledig in het district Muette. Laten we het overlappende deel bepalen.

GeoPandas en matplotlib zijn al geïmporteerd.

Deze oefening maakt deel uit van de cursus

Werken met georuimtelijke data in Python

Cursus bekijken

Oefeninstructies

Plot de twee polygonen in één kaart om visueel te bekijken in hoeverre ze overlappen.
Bereken de intersectie van de polygonen park_boulogne en muette.
Print het aandeel van de oppervlakte van het district dat door het park wordt ingenomen.

Praktische interactieve oefening

Probeer deze oefening eens door deze voorbeeldcode in te vullen.

# Plot the two polygons
geopandas.GeoSeries([park_boulogne, muette]).plot(alpha=0.5, color=['green', 'blue'])
plt.show()

# Calculate the intersection of both polygons
intersection = ____

# Plot the intersection
geopandas.GeoSeries([intersection]).plot()
plt.show()

# Print proportion of district area that occupied park
print(____ / ____)

Code bewerken en uitvoeren