Raster katmanından bilgi çıkar
Şimdi bir vektör dosyasına dayanarak raster katmanından bilgi çıkaralım. Bu işlevsellik rasterstats paketi tarafından sağlanır. Bu egzersizde özel olarak, her madencilik sahası noktasında en yakın raster piksel değerini alarak tüm sahalardaki bitki örtüsü türünü belirleyeceğiz.
Madencilik sahaları veri kümesinin (mining_sites) bir alt kümesi zaten yüklendi ve GeoPandas ile matplotlib zaten içe aktarılmış durumda.
Bu egzersiz, kursun bir parçasıdır
Python ile Coğrafi Verilerle Çalışma
Egzersiz talimatları
rasterstatspaketini içe aktar.rasterstats.point_query()fonksiyonunu kullanarak madencilik sahası konumları (mining_sites) için en yakın raster değerini çıkar ve bunu yeni bir'vegetation'sütununa ata.- Noktaları renklendirmek için
'vegetation'sütununu kullanarak madencilik sahası verisinin bir grafiğini oluştur. Açıklama (legend) eklediğinden emin ol.
Uygulamalı etkileşimli egzersiz
Bu egzersizi bu örnek kodu tamamlayarak deneyin.
# 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()