Ekstrak informasi dari layer raster
Sekarang mari kita ekstrak informasi dari layer raster berdasarkan berkas vektor. Fungsionalitas ini disediakan oleh paket rasterstats. Secara khusus untuk latihan ini, kita akan menentukan jenis vegetasi di semua lokasi pertambangan dengan mengambil nilai piksel raster terdekat pada setiap titik dalam himpunan data lokasi pertambangan.
Subset dari himpunan data lokasi pertambangan (mining_sites) sudah dimuat, dan GeoPandas serta matplotlib sudah diimpor.
Latihan ini adalah bagian dari kursus
Bekerja dengan Data Geospasial di Python
Petunjuk latihan
- Impor paket
rasterstats. - Ekstrak nilai raster terdekat untuk lokasi pertambangan (
mining_sites) menggunakan fungsirasterstats.point_query(), lalu tetapkan hasilnya ke kolom baru'vegetation'. - Buat plot data lokasi pertambangan dengan menggunakan kolom
'vegetation'untuk mewarnai titik. Pastikan untuk menambahkan legenda.
Latihan interaktif praktis
Cobalah latihan ini dengan menyelesaikan kode contoh berikut.
# 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()