Snittet av två polygoner
I den här övningen använder vi 2 enskilda polygoner: distriktet Muette, hämtat från datamängden districts, och det gröna stadsområdet Boulogne – en stor offentlig park i västra Paris – hämtat från datamängden land_use. De två polygonerna har redan tilldelats variablerna muette och park_boulogne.
Vi visualiserar först de två polygonerna. Du ser att de överlappar varandra, men att parken inte är helt belägen inom distriktet Muette. Låt oss ta reda på vilken del som överlappar.
GeoPandas och matplotlib är redan importerade.
Den här övningen är en del av kursen
Arbeta med geospatial data i Python
Övningsinstruktioner
- Rita de två polygonerna i en och samma karta för att visuellt undersöka hur mycket de överlappar.
- Beräkna snittet av polygonerna
park_boulogneochmuette. - Skriv ut hur stor andel av distriktets yta som upptas av parken.
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# 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(____ / ____)