2017年01月29日

幹材積計算プログラム

自分で胸高直径や樹高を計測した,かわいい杉の木達の材積を求めるために,他人が作成した計算式を自分のExcel表に入れ,セル番地に書き込まれていた定数を直接式の中に書き込んで計算させています。

理屈も何もわからず,ただD(直径)とH(樹高)から計算するだけです。

同僚が作成した計算式が正しいのかどうかもわかりません。

と思ったら,この分野のオーソリティーである,国立研究開発法人の森林総合研究所さんがExcelのマクロを公開していました↓。
http://www.ffpri.affrc.go.jp/database/stemvolume/

早速ダウンロードして使ってみました。

マクロを有効にした後の画面です↓。
材積計算例0.tif.gif

富山県のスギは,”名古屋スギ”に含まれます↓。
材積計算例1.tif.gif

名古屋スギの計算結果の例です↓。
材積計算例2.tif.gif

困ったことになりました。
D=31.0cm,H=23.3mのスギの幹材積は,
(1)同僚の計算式では,V=0.8397m3
(2)私の計算式では,V=0.8381m3
(3)森総研マクロでは,V=0.8305m3
と,それぞれ違った答えになりました。

元の計算式の根拠はこれ↓に載っているはずなのですが,
http://www.ffpri.affrc.go.jp/labs/shukakushiken/03gyomu/03gyomu_11.pdf
(1),(2)が(3)より過大になっています。

どうしましょうか?

追伸)
(1)と(2)は一致しました。
私の計算式の定数項の桁数が小数点以下3桁になっていたのですが,同僚のものはよく見ると点以下6桁になっていました。
私の式も6桁にしたら(1)と(2)の答えは完全に一致しました↓。
材積計算例3.gif

一歩前進です。

追伸2)
3者の答えが一致しました↓。
材積_小林.gif

材積_総研.gif

3者の元の計算式は,以下の式に基づいています。
材積式.gif

この元式の直径範囲の境目(10-12,30-32,40-42)で誤差が出るようで,森林総研のプログラムは,この境目で何らかの補正をしているらしい?とのことでした。

既に退職された大先輩には,違いはゴミのようなものなので気にしないようにと言われました。

我が社方式で前進します。

ただ,常識の定数(例えば5.770734,5の上にマイナス記号)の意味は大先輩も忘れたそうで,とにかく絶対値を足したら約10.0になるということ(-4.2293)だそうです。

う〜〜ん。

追伸3)
足しても10.0にならないのが4行目の式です。
誤植?があるそうで,正しくはこんな感じだそうです↓。
s34幹材積式.gif

エクセルマクロを開発された,森林総研の細田さんの論文でも,4行目の定数が,
-3.923549
となっていることを確認しました。
細田和男・光田 靖・家原敏郎「現行立木幹材積表と材積式による計算値との相違およびその修正方法」 森林計画学会誌44巻2号:23〜39ページ、2010年12月発行

材積を計算するのも大変ですね。


posted by Dr.koba at 10:13| Comment(0) | TrackBack(0) | 林業 | このブログの読者になる | 更新情報をチェックする

2016年04月26日

青焼きコピーのカシミール3Dでの表示について

コメント欄に投稿できない(ちょくちょくエラーが出て本当に困ります)ので新たな記事として質問にお答えします。

画像を直接カシミール3Dで読むと,誤差が出てまずいです。
前にも書いたように,画像データの測地系,座標系を変換するには,全ての画素(ピクセル)を並べ替える(リサンプリングとかワープとか呼びます)必要があります。

拙ブログの以下の記事
・森林計画図のgeoref(1)
http://gpsrsgis.seesaa.net/article/386679119.html
・森林計画図のgeoref(2)
・森林計画図のwarp
・森林計画図のbmp出力
・森林計画図をカシミール3Dで開く
を参考に,やってみてください。

QGISを使うのは厳しいでしょうか?


posted by Dr.koba at 08:09| Comment(1) | TrackBack(0) | 林業 | このブログの読者になる | 更新情報をチェックする

2016年03月24日

林業研究グループ総会

3/22(火)に富山林研グループの総会がありました。

林研総会1.gif

林研総会2.gif

総会後約50分間にわたり、「森林境界明確化支援システム」に関する講演をさせてもらいました。

林研総会3.gif

森林所有者さんたちに直接お目にかかれる貴重な機会をいただき、ありがとうございました。


posted by Dr.koba at 08:35| Comment(0) | TrackBack(0) | 林業 | このブログの読者になる | 更新情報をチェックする