Croiser un GeoDataFrame avec un polygone
En combinant l'ensemble de données sur l'utilisation du sol et celui des arrondissements, nous pouvons maintenant examiner l'utilisation du sol dans un arrondissement donné.
Pour cela, nous devons d'abord déterminer l'intersection de l'ensemble sur l'utilisation du sol avec un arrondissement précis. Reprenons l'arrondissement de la Muette comme exemple.
Les ensembles de données sur l'utilisation du sol et les arrondissements ont déjà été importés sous land_use et districts, et l'arrondissement de la Muette a été extrait dans le polygone shapely muette. De plus, GeoPandas et matplotlib sont importés.
Cette activité fait partie du cours
Travailler avec des données géospatiales en Python
Instructions de l’exercice
- Calculez l'intersection des polygones
land_useavec le polygone uniquemuette. Nommez le résultatland_use_muette. - Faites un tracé rapide de cette intersection et passez
edgecolor='black'pour mieux voir 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 ce code d’exemple.
# 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.____)