Część wspólna dwóch wielokątów
W tym ćwiczeniu skorzystamy z 2 pojedynczych wielokątów: dzielnicy Muette wyodrębnionej ze zbioru danych districts oraz zielonego terenu miejskiego Lasku Bulońskiego – dużego parku publicznego na zachodzie Paryża – wyodrębnionego ze zbioru danych land_use. Oba wielokąty zostały już przypisane do zmiennych muette i park_boulogne.
Na początek zwizualizujemy oba wielokąty. Zobaczysz, że zachodzą na siebie, ale park nie leży w całości w granicach dzielnicy Muette. Sprawdźmy, jaki obszar mają wspólny.
GeoPandas i matplotlib są już zaimportowane.
To ćwiczenie jest częścią kursu
Praca z danymi geoprzestrzennymi w Pythonie
Instrukcje do ćwiczenia
- Narysuj oba wielokąty na jednej mapie, aby wizualnie ocenić stopień ich nakładania się.
- Oblicz część wspólną wielokątów
park_boulogneimuette. - Wyświetl, jaką część obszaru dzielnicy zajmuje park.
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# 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(____ / ____)