Průnik dvou polygonů
V tomto cvičení budeme pracovat se dvěma samostatnými polygony: čtvrtí Muette vyextrahovanou z datasetu districts a zelenou plochou Boulogne – rozsáhlým veřejným parkem na západě Paříže – vyextrahovanou z datasetu land_use. Oba polygony jsou už přiřazeny do proměnných muette a park_boulogne.
Nejprve si oba polygony vizualizujeme. Uvidíš, že se překrývají, ale park neleží celý ve čtvrti Muette. Pojďme zjistit, která část se překrývá.
GeoPandas a matplotlib jsou už naimportované.
Toto cvičení je součástí kurzu
Práce s geoprostorovými daty v Pythonu
Pokyny k cvičení
- Vykresli oba polygony do jedné mapy a vizuálně posuď míru jejich překryvu.
- Vypočítej průnik polygonů
park_boulogneamuette. - Vypiš, jakou část plochy čtvrti park zaujímá.
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
# 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(____ / ____)