Extraire des informations à partir d’une couche raster
Extrayons maintenant des informations de la couche raster à partir d’un fichier vecteur. Cette fonctionnalité est fournie par le paquet rasterstats. Dans 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 proche 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 ainsi que matplotlib sont déjà importés.
Cet exercice fait partie du cours
Travailler avec des données géospatiales en Python
Instructions
- Importez le paquet
rasterstats. - Extrayez la valeur raster la plus proche pour les emplacements des sites miniers (
mining_sites) à l’aide de la fonctionrasterstats.point_query(), et affectez-la à une nouvelle colonne'vegetation'. - Réalisez une visualisation des données des sites miniers en utilisant la colonne
'vegetation'pour colorer les points. Veillez à ajouter une légende.
Exercice interactif pratique
Essayez cet exercice en complétant cet exemple de code.
# 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()