Estrai informazioni dal layer raster
Ora estraiamo informazioni dal layer raster a partire da un file vettoriale. Questa funzionalità è fornita dal pacchetto rasterstats. In particolare, in questo esercizio determineremo il tipo di vegetazione in tutti i siti minerari, recuperando il valore del pixel raster più vicino per ciascun punto dell'insieme di dati sui siti minerari.
Un sottoinsieme dell'insieme di dati dei siti minerari (mining_sites) è già caricato, e GeoPandas e matplotlib sono già importati.
Questo esercizio fa parte del corso
Lavorare con i dati geospaziali in Python
Istruzioni dell'esercizio
- Importa il pacchetto
rasterstats. - Estrai il valore raster più vicino per le posizioni dei siti minerari (
mining_sites) usando la funzionerasterstats.point_query()e assegna il risultato a una nuova colonna'vegetation'. - Crea un grafico dei dati dei siti minerari usando la colonna
'vegetation'per colorare i punti. Assicurati di aggiungere una legenda.
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
# 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()