Intersection de deux polygones
Pour cet exercice, nous allons utiliser 2 polygones distincts : l'arrondissement de Muette extrait de l'ensemble de données districts, et la zone verte urbaine de Boulogne, un grand parc public à l'ouest de Paris, extraite de l'ensemble de données land_use. Les deux polygones ont déjà été assignés aux variables muette et park_boulogne.
Nous allons d'abord visualiser les deux polygones. Vous verrez qu'ils se chevauchent, mais le parc n'est pas entièrement situé dans l'arrondissement de Muette. Déterminons la partie de chevauchement.
GeoPandas et matplotlib sont déjà importés.
Cette activité fait partie du cours
Travailler avec des données géospatiales en Python
Instructions de l’exercice
- Tracez les deux polygones sur une même carte pour examiner visuellement le degré de chevauchement.
- Calculez l'intersection des polygones
park_boulogneetmuette. - Affichez la proportion de la superficie de l'arrondissement occupée par le parc.
Exercice interactif pratique
Essayez cet exercice en complétant ce code d’exemple.
# 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(____ / ____)