點的緩衝區
考慮剛果(金)北基伍省的省會戈馬(Goma),鄰近盧安達邊境。它的座標為 1.66°S、29.22°E(Point 已以 UTM 座標給定,存在變數 goma 中)。
有多少採礦據點位於戈馬 50 km 範圍內?國家公園涵蓋多少面積?讓我們用緩衝區(buffer)運算來判定。記得距離要以 CRS 的單位表示(本例為公尺)。
注意:如果你有一個布林 Series(例如空間關係方法的結果),要計算 True 的個數(也就是通過檢查的幾何數),可以對該布林值取總和,因為在這種情況下 True 和 False 會分別視為 1 和 0。
本練習屬於課程
在 Python 中處理地理空間資料
練習說明
- 以
goma建立 50 km 的緩衝區,並指定為變數goma_buffer。 - 檢查
mining_sites中哪些點位於此緩衝區內(將其稱為mask)。對布林mask取總和以計算這些據點數量,並印出結果。 - 計算國家公園與該緩衝區的交集。將這些交集的面積加總,並以 km² 印出結果。
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
# goma is a Point
print(type(goma))
# Create a buffer of 50km around Goma
goma_buffer = ____
# The buffer is a polygon
print(type(goma_buffer))
# Check how many sites are located within the buffer
mask = ____
print(____)
# Calculate the area of national park within the buffer
print(national_parks.____ / (1000**2))