Intersecția a două poligoane
În acest exercițiu vom folosi 2 poligoane individuale: districtul Muette, extras din setul de date districts, și zona verde urbană Boulogne, un mare parc public din vestul Parisului, extras din setul de date land_use. Cele două poligoane au fost deja atribuite variabilelor muette și park_boulogne.
Începe prin a vizualiza cele două poligoane. Vei observa că se suprapun, însă parcul nu se află în întregime în districtul Muette. Să determinăm zona de suprapunere.
GeoPandas și matplotlib sunt deja importate.
Acest exercițiu face parte din cursul
Lucrul cu date geospațiale în Python
Instrucțiuni pentru exercițiu
- Reprezentați grafic cele două poligoane pe o singură hartă pentru a examina vizual gradul de suprapunere.
- Calculează intersecția poligoanelor
park_boulogneșimuette. - Afișează proporția din suprafața districtului ocupată de parc.
Exercițiu interactiv practic
Încearcă acest exercițiu completând acest cod de exemplu.
# 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(____ / ____)