Extraire de l'information d'une couche raster
Extrayons maintenant de l'information à partir de la couche raster, en nous basant sur un fichier vecteur. Cette fonctionnalité est offerte par le paquet rasterstats. Plus précisément pour cet exercice, nous allons déterminer le type de végétation à tous les sites miniers, en récupérant la valeur du pixel raster le plus près pour chaque point du jeu de données des sites miniers.
Un sous-ensemble du jeu de données des sites miniers (mining_sites) est déjà chargé, et 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
- Importez le paquet
rasterstats. - Extrayez la valeur raster la plus proche pour les emplacements des sites miniers (
mining_sites) en utilisant la fonctionrasterstats.point_query(), et assignez-la à une nouvelle colonne'vegetation'. - Tracez les données des sites miniers en utilisant la colonne
'vegetation'pour colorer les points. Assurez-vous d'ajouter une légende.
Exercice interactif pratique
Essayez cet exercice en complétant ce code d’exemple.
# Import the rasterstats package
____
# Extract the nearest value in the raster for all mining sites
vegetation_raster = "central_africa_vegetation_map_foraf.tif"
mining_sites['vegetation'] = ____(____, ____, interpolate='nearest')
print(mining_sites.head())
# Replace numeric vegation types codes with description
mining_sites['vegetation'] = mining_sites['vegetation'].replace(vegetation_types)
# Make a plot indicating the vegetation type
____
plt.show()