Extragerea informațiilor dintr-un strat raster
Acum vom extrage informații din stratul raster pe baza unui fișier vectorial. Această funcționalitate este oferită de pachetul rasterstats. În cadrul acestui exercițiu, vom determina tipul de vegetație din dreptul fiecărui sit minier, obținând valoarea celui mai apropiat pixel raster pentru fiecare punct din setul de date.
Un subset al setului de date cu situri miniere (mining_sites) este deja încărcat, iar GeoPandas și matplotlib sunt deja importate.
Acest exercițiu face parte din cursul
Lucrul cu date geospațiale în Python
Instrucțiuni pentru exercițiu
- Importă pachetul
rasterstats. - Extrage valoarea raster cea mai apropiată pentru locațiile siturilor miniere (
mining_sites) folosind funcțiarasterstats.point_query()și atribuie rezultatul unei coloane noi numite'vegetation'. - Realizează un grafic al datelor despre siturile miniere, folosind coloana
'vegetation'pentru a colora punctele. Asigură-te că adaugi și o legendă.
Exercițiu interactiv practic
Încearcă acest exercițiu completând acest cod de exemplu.
# 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()