Extrahera information från rasterlager
Nu ska vi extrahera information från rasterlagret utifrån en vektorfil. Den här funktionaliteten tillhandahålls av paketet rasterstats. I den här övningen ska vi bestämma vegetationstypen vid samtliga gruvplatser genom att hämta det närmaste rasterpixelvärdet för varje punkt i gruvplatsdatamängden.
En delmängd av gruvplatsdatamängden (mining_sites) är redan inläst, och GeoPandas och matplotlib är redan importerade.
Den här övningen är en del av kursen
Arbeta med geospatial data i Python
Övningsinstruktioner
- Importera paketet
rasterstats. - Extrahera det närmaste rastervärdet för gruvplatsernas positioner (
mining_sites) med funktionenrasterstats.point_query(), och tilldela resultatet till en ny kolumn'vegetation'. - Skapa ett diagram över gruvplatsdatan där du använder kolumnen
'vegetation'för att färgsätta punkterna. Se till att lägga till en legend.
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# 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()