Informatie uit een rasterlaag halen

Laten we nu informatie uit de rasterlaag halen op basis van een vectorbestand. Deze functionaliteit wordt aangeboden door het pakket rasterstats. Specifiek voor deze oefening bepalen we het vegetatietype op alle mijnlocaties door voor elk punt in de gegevensset met mijnlocaties de dichtstbijzijnde rasterpixelwaarde op te halen.

Een subset van de gegevensset met mijnlocaties (mining_sites) is al ingeladen, en GeoPandas en matplotlib zijn al geïmporteerd.

Deze oefening maakt deel uit van de cursus

Werken met georuimtelijke data in Python

Cursus bekijken

Oefeninstructies

Importeer het pakket rasterstats.
Haal de dichtstbijzijnde rasterwaarde op voor de mijnlocaties (mining_sites) met de functie rasterstats.point_query() en zet dit in een nieuwe kolom 'vegetation'.
Maak een plot van de mijnlocaties en gebruik de kolom 'vegetation' om de punten te kleuren. Voeg zeker een legenda toe.

Praktische interactieve oefening

Probeer deze oefening eens door deze voorbeeldcode in te vullen.

# 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()

Code bewerken en uitvoeren