Intersection de deux polygones
Dans cet exercice, nous allons utiliser 2 polygones individuels : le quartier de Muette, extrait du jeu de données districts, et la zone urbaine verte de Boulogne, un grand parc public à l’ouest de Paris, extrait du jeu de données land_use. Les deux polygones ont déjà été affecté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 le quartier de Muette. Déterminons la partie qui se superpose.
GeoPandas et matplotlib sont déjà importés.
Cet exercice fait partie du cours
Travailler avec des données géospatiales en Python
Instructions
- Tracez les deux polygones sur une seule carte pour évaluer visuellement le chevauchement.
- Calculez l’intersection des polygones
park_boulogneetmuette. - Affichez la proportion de la surface du quartier occupée par le parc.
Exercice interactif pratique
Essayez cet exercice en complétant cet exemple de code.
# 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(____ / ____)