Extrakce informací z rastrové vrstvy
Teď z rastrové vrstvy extrahujeme informace na základě vektorového souboru. Tuto funkcionalitu poskytuje balíček rasterstats. V tomto cvičení konkrétně zjistíme typ vegetace ve všech těžebních lokalitách – a to tak, že pro každý bod v datasetu těžebních lokalit získáme hodnotu nejbližšího rastrového pixelu.
Podmnožina datasetu těžebních lokalit (mining_sites) je již načtena a GeoPandas spolu s matplotlib jsou také k dispozici.
Toto cvičení je součástí kurzu
Práce s geoprostorovými daty v Pythonu
Pokyny k cvičení
- Importuj balíček
rasterstats. - Pomocí funkce
rasterstats.point_query()extrahuj hodnotu nejbližšího rastrového pixelu pro lokality těžebních míst (mining_sites) a výsledek ulož do nového sloupce'vegetation'. - Vytvoř graf dat těžebních lokalit, kde barva bodů bude odpovídat sloupci
'vegetation'. Nezapomeň přidat legendu.
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
# 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()