Intersezione di due poligoni
Per questo esercizio useremo 2 poligoni singoli: il distretto di Muette estratto dall'insieme di dati districts, e l'area verde urbana di Boulogne, un grande parco pubblico a ovest di Parigi, estratta dall'insieme di dati land_use. I due poligoni sono già stati assegnati alle variabili muette e park_boulogne.
Per prima cosa visualizziamo i due poligoni. Vedrai che si sovrappongono, ma il parco non ricade interamente nel distretto di Muette. Determiniamo la parte in sovrapposizione.
GeoPandas e matplotlib sono già importati.
Questo esercizio fa parte del corso
Lavorare con i dati geospaziali in Python
Istruzioni dell'esercizio
- Rappresenta i due poligoni in un'unica mappa per esaminare visivamente il grado di sovrapposizione
- Calcola l'intersezione dei poligoni
park_boulogneemuette. - Stampa la proporzione dell'area del distretto occupata dal parco.
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
# 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(____ / ____)