Trích xuất thông tin từ lớp raster
Giờ hãy trích xuất thông tin từ lớp raster dựa trên một tệp vector. Chức năng này được cung cấp bởi gói rasterstats. Cụ thể trong bài này, bạn sẽ xác định loại thảm thực vật tại tất cả các điểm khai thác mỏ bằng cách lấy giá trị pixel raster gần nhất tại mỗi điểm trong bộ dữ liệu điểm khai thác mỏ.
Một tập con của bộ dữ liệu điểm khai thác mỏ (mining_sites) đã được nạp sẵn, và GeoPandas cùng matplotlib đã được nhập sẵn.
Bài tập này là một phần của khóa học
Làm việc với Dữ liệu Không gian địa lý trong Python
Hướng dẫn bài tập
- Import gói
rasterstats. - Trích xuất giá trị raster gần nhất cho các vị trí điểm khai thác mỏ (
mining_sites) bằng hàmrasterstats.point_query()và gán vào cột mới'vegetation'. - Vẽ biểu đồ dữ liệu điểm khai thác mỏ, dùng cột
'vegetation'để tô màu các điểm. Nhớ thêm chú giải (legend).
Bài tập tương tác thực hành trực tiếp
Hãy thử làm bài tập này bằng cách hoàn thành đoạn mã mẫu này.
# 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()