ぎょーむ日誌 2001-06-06

2001 年 06 月 06 日 (水)

• 0730 起床． 札幌は今朝も快晴．
• 午前中費して Pasoh シミュレイターの改善につとめる． Voxels の初期化を簡単化して， パラメーターファイルの数をひとつ減じることができたのは 良いのだが …… すべての計算が終った仮想森林を細部・末端から 逐次的に「あの世」から delete していくときに， Voxel の持つ近傍 Voxel のリストを delete すると core を吐くことがある． うーん， わざわざ安全を確認してから「消す」 ように配慮してるんだが．
• 結局これは， free の STL である libstdc++ に問題あったようだ． 最新版に交換するととりあえず core は出なくなった． 時刻はすでに 1300． うう． ちょっと昼飯にしますか． 甲山家の庭でとれたイタリアンパセリを恵んでいただいたので， それを炒めてパスタソースにする．
• パセリパスタを食べてると M2 の長谷川さんが 石鹸をめぐんでくださる． 献血のときにもらったそうだ． うう， ありがとうございます． そして今日こそ銭湯に行きたい．
• 昼からも Pasoh 解析の続き． 新ライブラリの元で開空度ごとき指標の計算， というのをまだやっていなかったことにようやく気づく． 前のライブラリとだいたい同じ結果を出そうとすると， あれこれ細かいバグが見つかる． ライブラリまで遡って直したりする．
• 1515 ようやく以前と比較可能な計算の試験開始． 1ha 区画というかその中の 50x50m 区画について計算するのに 15 秒ぐらいかかってるから， 171 区画を計算するのに 40 分ぐらいかかるのか． やれやれ．
• えーと， これは樹冠を構成するタマの直径が 1m の場合で， 以前との比較のためにやっているのである． 直径を 2m ぐらいにすると 処理速度 4 倍ぐらいにならないだろうか? そのための新ライブラリなんだが ……
• 計算待ちのあいだ， ちょっと Numerical Recipes を見て統計的手法の おさらい． うーん， 二次元のτ検定みたいなことはやりたくないなぁ． どうしよう．
• これから計算せんといかんのは ……
  • 樹種ごとの 成長パラメーター推定． これは従来どーり．
  • 樹種ごとの 新規加入速度の推定量を出す． 甲山さん方式． 区画など取ったりするんで， いまいちロバストネスがわからん統計量だ． さーてと， 莫大な本数の樹木データ流れる 処理のどこに分岐を作って そういう計算に回せばいいのやら ……
  • 他にもいくつかあったな …… やはりいままで単線的な処理ばかりしてたので， 途中で分岐作るのちと注意せんといかんな． 今のような Makefile 使った方法では そのへんの管理が難しいか?
• 1600 ごろ開空度の計算終ってるのに気づく． やはり 45 分ぐらいかかるな． じゃあ， これをもとに樹種ごとにパラメーター推定． 時間かかるんだよね， これも …… そうか， 3 パラメーターモデルだけじゃなくて， 2 パラメーターモデルの場合も 推定するようにしていたんだった． ますます時間かかるな．
• 38 万本の樹木個体とか 470 樹種とか …… 何をやってるのか だんだんわからなくなってくるところがある． 後から簡単に全処理を再現できるようにしなければならない， というきつい制約もある (としておいたおかげで， ライブラリ変えようがモデルを変更しようが 何度でもやりなおしがきいてるわけだが)．
• 1630 より Trendy セミナー， 永光さん (森林総研北海道) の シウリザクラ というクローナルな樹木の話である． 苫小牧演習林の毎木と遺伝子の調査から， Genet の境目がはっきりしてるとか， 養子の多い親樹は苦労してるみたいだ， といった楽しい内容であった．
• セミナー後， 演者の永光さんほか 内田さん・浦口さん・かとー大先生と 北15西4 の 「中華美食楼」 で会食． おお， 軟禁されていた建物の外に出るのは 48 時間ぶりだ．
• 帰りがけに銭湯 「江札湯」 に入る． じつは札幌に来てから初めての入浴である． 客少なく浴槽独占状態． ああ， さっぱり．
• もどってから Pasoh 開空度計算のつづき． 樹冠モジュールなるつぶつぶの直径を 1m→2m にすると 計算時間が 45 分から 28 分に短縮された． ファイル読み書きの時間は以前と変わらぬから 半減とまではいかないのである．

[樹冠模型] 4 月に描いたやつを再掲． 50ha 調査区の中の 1ha を抽出して 三次元モデル化．

• つぶつぶの大きさだけ変えてもパラメーター推定結果は そんなに変わらんな (樹冠の形状だのつぶつぶ表面積合計は同等)． よしよし．

[パラメーター推定の結果の例] 赤は「つぶつぶ 1m」 青は「つぶつぶ 2m」 の場合． 横軸は成長率の定数項の係数の最尤推定値， 縦軸は明る依存項の係数の最尤推定値． 一点一樹種である．

• 今日の食卓
  • 朝 (0750): 米 0.6 合． タマネギと豆腐の炒めもの． ジンギスカンソース味． 書き忘れていたけど「お茶部屋料理」な 日々である．
  • 昼 (1300): スパゲッティー． ソースは甲山さんにいただいたパセリを 炒めてトマトケチャップ・キムチの素で味つけしたもの．
  • 晩 (1830): 北15西4 の 中華美食楼で餃子定食 (630 円)．