Wyodrębnianie informacji z warstwy rastrowej
Teraz wyodrębnimy informacje z warstwy rastrowej na podstawie pliku wektorowego. Tę funkcjonalność zapewnia pakiet rasterstats. W tym ćwiczeniu określimy typ roślinności dla wszystkich miejsc wydobycia, pobierając wartość najbliższego piksela rastrowego dla każdego punktu z naszego zbioru danych.
Podzbiór zbioru danych o miejscach wydobycia (mining_sites) jest już wczytany, a biblioteki GeoPandas i matplotlib są już zaimportowane.
To ćwiczenie jest częścią kursu
Praca z danymi geoprzestrzennymi w Pythonie
Instrukcje do ćwiczenia
- Zaimportuj pakiet
rasterstats. - Wyodrębnij wartość najbliższego piksela rastrowego dla lokalizacji miejsc wydobycia (
mining_sites), używając funkcjirasterstats.point_query(), i przypisz wynik do nowej kolumny'vegetation'. - Utwórz wykres danych o miejscach wydobycia, używając kolumny
'vegetation'do kolorowania punktów. Pamiętaj, aby dodać legendę.
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# 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()