Featured image of post 一般森林管理與 GIS — 定位最佳砍伐區

一般森林管理與 GIS — 定位最佳砍伐區

在指定的條件下,找出可以砍伐樹木的區域。

目標

從 forest 與 shoreline 這兩張圖之中找出可以砍樹的地方。

條件

  1. Only White Pine
  2. 距離水源超過 1000 公尺
  3. 面積至少要 1000 公頃
  4. 與道路(roads)的距離最多 500 公尺

Marco Modeler 操作

Only White Pine

從 forest 中分類出我們要的 White Pine,使用 Edit 建立一個 avl. 檔,取名為 white pine。內容如下:

1
2
3
1 1
2 0
3 0

再將 white pine 與 forest 使用 assign output 出只有 White Pine 的像素檔,取名為 wpine,結果如下:

1 的部分為我們需要的 White Pine 區域。

滿足第一個條件的模型

距離水源超過 1000 公尺

要從 shoreline 中分類出水源區,再建立一個 avl. 檔,取名為 water。

1
2
1 0
2 1

也一樣再進行 assign,output 的檔案取名為 water。

1 的部分為水源區。

但這個圖是水源區,要用 buffer 取出距離水源區 1000 公尺的區域。buffer 內 Additional parameters 的設定為:

  • Value for target:在指定的區域內 → 我們沒有要在水源區內,所以設定為 0
  • Value for buffer zone:在指定的區域帶上 → 在水源區上我們也不要,所以也是設定為 0
  • Value for area outside buffer:在指定的區域外 → 我們要在水源區外,所以設定為 1
  • Buffer distance:設定距離,為 1000

再連接剛剛的 water,output 的檔案取名為 waterbuffer,結果跑出來如下:

黑色部分為 buffer 內設定的前兩個部分;紅色為離水源區 1000 公尺外的區域。

滿足第二個條件的模型

面積至少要 1000 公頃

wpine 與 water buffer 使用 overlay 相乘,output 的檔案取名為 pinebuff。

再使用 group 進行分組編號,output 的檔案取名為 classification。

再使用 area 計算面積,Unit type 要選擇 Hectares(公頃),計算後的檔案取名為 pinebuffarea。

再用 reclass 把 pinebuffarea 面積小於 1000 公頃的區域刪去,先使用 Edit 建立一個 rcl. 檔,取名為 pinebuffclass。

1
2
3
1 1000 6100
0 0 999
0 6100 999999

放入 reclass 裡,output 的檔案取名為 area0n1。

最後再把 area0n1 與 pinebuffarea 使用 overlay 相乘,最後 output 的檔案取名為 lease。

滿足第三個條件的模型

與道路(roads)的距離最多 500 公尺

使用 intital 連接任何一張像素檔,output 一個空白的像素圖,取名為 blank。再與 roads 連接 lineras,把 roads 這個向量檔轉為像素檔,output 的像素檔也取名為 roads。

再算圖上的點到道路的距離,使用 distance 連接,output 的檔案取名為 dis2road。

使用 extract 功能萃取出離道路較遠的區域。在 extract 的設定中,Output type 選擇 Image file,Summary type 選擇 Min。Feature definition image 連接 lease,Image to be processed 連接 dist2road,output 出的檔案取名為 far。extract 內的資料如下:

output 出的 far 如下圖:

這兩塊地離 roads 太遠,所以要從原本的圖去掉這兩塊地。

使用 reclass 把圖分為 0 空白與 1 有顏色的區域。reclass 需要的 rcl. 檔可以如下編輯,因為兩塊地的數值分別為 2241 跟 2128,所以只要 1 的範圍有包含這兩個數值就好。

1
2
0 0 2000
1 2001 2300

進行 reclass 後 output 出的檔案取名為 farclass,從原本有數值的圖分成 0 與 1。

這兩塊是我們不要的,所以要原本的圖(area0n1)減去只有這兩塊的圖(farclass),使用 overlay 中的減法,要先連接 area0n1 再連 farclass,因為 overlay 的減法是 first - second。減完後 output 的檔案取名為 -over,就是我們需要的。

再與面積(pinebuffarea)使用 overlay 相乘,最後出來的圖取名為 final,就是滿足四個條件的可以砍樹的地方了。

滿足第四個條件的模型

結果

結果圖

final 去掉黑色背景、加入 roads 與圖例,使可以使用的區域看起來更明確。

Marco Modeler 的所有流程

comments powered by Disqus
使用 Hugo 建立
主題 StackJimmy 設計