Croiser un GeoDataFrame avec un polygone
En combinant le jeu de données d’occupation du sol et celui des arrondissements, nous pouvons maintenant analyser l’occupation du sol dans un arrondissement donné.
Pour cela, nous devons d’abord déterminer l’intersection du jeu de données d’occupation du sol avec un arrondissement donné. Prenons de nouveau l’arrondissement de la Muette comme exemple.
Les jeux de données d’occupation du sol et des arrondissements ont déjà été importés sous les noms land_use et districts, et l’arrondissement de la Muette a été extrait dans le polygone shapely muette. Par ailleurs, GeoPandas et matplotlib sont importés.
Cet exercice fait partie du cours
Travailler avec des données géospatiales en Python
Instructions
- Calculez l’intersection des polygones
land_useavec le polygone uniquemuette. Nommez le résultatland_use_muette. - Réalisez un tracé rapide de cette intersection et passez
edgecolor='black'pour voir plus clairement les limites des différents polygones. - Affichez les cinq premières lignes de
land_use_muette.
Exercice interactif pratique
Essayez cet exercice en complétant cet exemple de code.
# Print the land use datset and Notre-Dame district polygon
print(land_use.head())
print(type(muette))
# Calculate the intersection of the land use polygons with Notre Dame
land_use_muette = ____
# Plot the intersection
land_use_muette.____
plt.show()
# Print the first five rows of the intersection
print(land_use_muette.____)