報告要求如上
首先我們先取得大安森林公園的圖資,然後進行分層的作業
這邊插入一段超級抱怨文,你可以直接跳到下一段
這是我第一手拿到的圖資,想說人家已經把大安森林公園分好區了
我可以粗估覆蓋率,計算公園總面積等等
但程式怎麼也不讓我算面積,我眉頭一皺發現案情並不單純
你看看阿你看看,這張GIS的圖資竟然是線圖!!!
大姊阿不是把線框起來就是多邊形(polygon)阿!!
你有沒有修過GIS的課程啊???
虧你今年還考上林業技術的高普考=口=!?!?!?
你什麼都沒看到,我們繼續來討論如何分層
主要分成四層0-無法測量的地區,如水池,建物等等
1-可測量,林木覆蓋率佔10%左右之地區
2-林木覆蓋率在30%~70%的地區
3-林木覆蓋率90%以上的地區
分層的取決以目測為主
輸出圖資如下
隨便輸出的,只是要讓大家看一下大致上如何分配
這個圖可以看出森林的分布情況和計算面積
基本上如果包含建築物接近25公頃,符合網路上查的資料
再來是利用GIS軟體依照面積比例隨機佈上點位
我想佈20個樣區,因此照面積比例分成8,6,6下去
輸出如下
我有標上TWD97的座標,雖然在現場時我還是用經緯度的GPS找點位
然後就是今天早上實地去測量
發現這張航照圖太舊了,森林北邊的幾點可能受到捷運施工影響無法進入
只好另外挑點位
另外航照圖有"偏移"的現象,會和GPS經緯度搭不上,尤其是小範圍的公園內
這點我們事後討論
將20個點的樹高先繪成樹高曲線,我是用R做的
> summary(treeH.lm)
Call:
lm(formula = Height ~ D.B.H, data = treeH)
Residuals:
Min 1Q Median 3Q Max
-5.0637 -1.7138 -0.1226 1.9490 4.7767
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.4153 1.0064 2.400 0.0226 *
D.B.H 0.2765 0.0391 7.071 6.08e-08 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 2.441 on 31 degrees of freedom
Multiple R-squared: 0.6173, Adjusted R-squared: 0.605
F-statistic: 50.01 on 1 and 31 DF, p-value: 6.078e-08
不能說很漂亮,畢竟 R-squared只有0.6173
但我先做殘差分析
總之還算合格,細節我會再討論
我接著測試如果樹高曲線是二次函數的話...
Call:
lm(formula = Height ~ D.B.H + I(D.B.H^2), data = treeH)
Residuals:
Min 1Q Median 3Q Max
-4.9497 -1.6394 -0.2133 1.9068 4.6348
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.643454 2.488276 0.660 0.5140
D.B.H 0.342113 0.196964 1.737 0.0927 .
I(D.B.H^2) -0.001146 0.003370 -0.340 0.7361
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 2.477 on 30 degrees of freedom
Multiple R-squared: 0.6188, Adjusted R-squared: 0.5934
F-statistic: 24.35 on 2 and 30 DF, p-value: 5.219e-07
更糟了,或許是多次函數(笑),我們暫且使用一次吧
所以樹高曲線的公式是
y=0.2765x+ 2.4153
(我要是看到哪組敢用這個公式小心我扁你,自己組的自己算)
(何況這誤差頗大)
接著算各樣區的材積統計
這東西好醜,總之最後出來誤差率高達55%
我想請組員們重算看看
也許我哪個步驟有出錯
至於要降低到10%的話,取樣點要高達700才行= =
這有不少要討論的
樹高曲線誤差
航照圖誤差
劃分區域不當
公園內樹木生長情況特殊
樹木分支增加誤差
明天見= =
沒有留言:
張貼留言