ぎょーむ日誌 2001-04-17

2001 年 04 月 17 日 (火)

• 0640 起床． 真夜中から全然寝られなかった． まぁ， 晩飯後に 3 時間ぐらい寝たけど．
• 朝飯・弁当の準備． シャワー． 朝飯． コーヒー．
• 0740 自宅発． 晴れ． 0759 東京モノレイル流通センター発． 浜松町の気温 16 ℃． 0820 研究所着．
• 早くも眠くなってきた． うう． 脳内ホワイトアウト．
• とりあえず Pasoh シミュレイター起動する → 50ha のうち 1ha だけ読ませる → 計算はまだいいや → 読んだものをそのまま出す， という作業やらせようとする．
• 自分で準備しておいたエラーメッセイジを何度も眺めたり， pasoh_simulator の吐いた core ファイルも解析しつつ， バグ取りというか OFR ⇔ Pasoh のすりあわせ作業を行う．
• 1015 必要なファイルを全部読ませることはできた． しかし出力がヘンだ．
• すでにこの独房群内は引越し屋さん多数が徘徊している． ばたばたして落ち着かない．
• Pasoh シミュレイターの出力． 書式がヘンだとわかった． 無意味な空白行が数千も出ていた．
• 計算はともかく， 読む → 出す， というのはできた． 時刻は 1130． じゃ， あとは途中の計算をちゃんとする， ということで． 具体的には新造した Canopy クラスの関数の 中身を書いてしまえばよい．
• 1155 アヤしげな開空度ごときものも計算されるようになった …… しかし， なんかやけにすかすかだな． そして， 1ha 中心の個体たちのうち約 5% の開空度 (らしき指標) が完全にゼロなのはなぜだろう …… あ， 「無眼」指定樹種かな? 「無眼」なら「無眼」と出力するようにしてみよう …… やはりそうか． Pasoh 第 155 (暫定) 計算区画 960 個体中 61 個体は 調査区内稀少樹種指定なので計算対象外， 開空度観測点を設置せず， というワケね．
• 「すかすか」なのは林冠木も含んでいるからだろうか …… 899 個体の中央値は約 0.55 か． なんか明るい． まぁ， このあたりは「樹冠の形状」と 葉群モジュール密度の調整で変わってくるだろう． 幸か不幸か， 小川より Pasoh のほうがこのあたりの パラメーターに関して すでに決められてしまってる部分が多い．
• まぁ， 午後は Canopy クラスをもうちょっと修正しよう． その前に昼飯にしますか．
• 昼飯食うとやはり眠くなったので， でれでれできる場所を求めて浜松町界隈を散歩したんだが …… そういう場所がない． ちょっと休めそうな場所は 喫煙者が群をなしている． 杉田 (横浜市金沢区) に移ったら 少しはマシになるのだろうか．
• チャンドラ部品 Yahoo! オークション落札者 4 名のうち， 連絡の取れなかったヒトがようやくメイル送ってくる． すでに送金ずみとのこと． ネット経由で口座の状況を確認． あ， 入ってるね． とはいえ， 発送は帰宅してからになるんだが．
• いきなり昏睡してしまいそうな眠気こらえつつ， 熱帯林の樹木の数値樹冠をでっちあげる作業． 概形は (これまた信用ならない) 円筒形として与えられているんだよね． ありがたいことに葉面積まで むりやり計算されている． どう使えばいいんだろうね．
• …… う． 座ったまま 10 分ほどサスペンドしてた． いやはや． しかし眠気はおさまった．
• さーてと， 樹冠問題の続き． さっきも書いたように樹冠の概形は すでに甲山さんがモデリングしており， dbh (胸高直径) がわかれば
  • 樹高
  • 樹冠直径
  • 枝下高
  • 幹重
  • 葉重
  • 葉面積
  までわかってしまうというような， なんというか魔法のようなハナシだ． ちなみに変換式は全樹種共通．
• 他に代替案ないので， これを信用したフリをして使うわけだが …… さーて， どうしてくれよう． 何も考えてないソリッドモデルは使いもんにならんし． とりあえず， (こちらも何の進歩もなく) つぶつぶ (粒状樹冠モジュール) モデルを使うとして， 上の怪しげな憶測値とどう対応つけるか， だな．
• 出張してたんで入金遅れると連絡してくれてた 落札者からも入金通知． これでお金のやりとり関係は全部終わったはず． あとは宅急便で発送するだけだ． 今晩， 自宅に帰ってからだけど．
• 樹冠モデリング問題 …… まずは樹冠モジュールの個数が 葉面積に比例， とするか． ちょっと気色悪いが． えーと， どうスケイリングすればいいんだ． 個数は (葉面積) / (Voxel の体積)^(2/3) …… とでもするかな．
• 次に， このつぶつぶの三次元空間分布だ． 「外枠」たる甲山円筒モデルとどう対応させるか …… これも面白くはないが， 円筒に内接する回転半楕円体上に一様分布， とでもしてみるかな．

一様分布なのは角度ね (座標変換前の)．
• …… とゆーふうに 樹冠モデルを変えて， 計算やりなおし． 時刻 1505． ふーん， 開空度の中央値が 0.1 以下になったね． かなり暗いというべきか． 結局このあたりは樹冠モジュールの個数だな． というか， 「透過率」みたいに翻訳できる 勝手なパラメーターがどうしても一個入る．
• 空間配置モデリングは小川シミュレイターの方式を 引きずっているので， 樹冠モジュールの配置に乱数を使っている． 今回は推定したパラメーター眺めてどうこう言おう， という心算らしいから， 三次元モデリングに乱数使わないほうがいいかなぁ． 独立反復試行がたいへんだ． 小川とちがって， この 30 数万本とか尋常じゃない個体数だしねー．
• 樹冠形状アイデアも尽きたんで， どこかにお伺いたててみようか． しかしソリッドモデルとか言われるだろうな． ひとつ真に受けたふりをしてソリッドでやるか? しかし使ってるライブラリはそういうことに 向いてない．
• 計算量が少なくて 確率論的じゃなくて決定論的で 見てくれもヘンではない樹冠モデルって ないもんかね． そもそも上のような値のセット渡されたときに， 何も考えずに概形というか境界が円筒 と信じてしまったあたりが莫迦なんだろうな．
• 何の問題解決にもなってないけど， もうちょっと確率論的樹冠形成で押してみよう． このモデルは計算量少なくて形がそんなにヘンではない (かも)， という利点だけはある． Pasoh についても視覚化してみるか， とりあえず． ヘンかどうかは見てから考えよう．
• 1635， 作画例できた．

[Pasoh 50ha 調査区の一部] 第 155 計算区画 (第 0 から第 170 まである)． クリックで拡大図． でかいのが一本あると思ってたら， エマージェント (林冠突出個体) だったのかな?

…… つらつらと眺めてるうちに， なんか計算間違いしてるような気がしてきた． 再点検 …… エマージェント (みたいな個体) の「スカート丈」が 短いような気がしたんだけど， あれでいいみたい． 甲山アロメトリーモデルによると， すべての個体において「スカート」の長さは 樹高の 36% ということになってるのである．
• よーし， このへんの妥当性はまた後から考えることにして 先に進もう．
• あとはこういう 1ha 計算を 171 回やっていけば， 50ha (500m × 1000m) 調査区の端から端まで 計算できるな． 数があってないのは， 1ha プロット真ん中の 50m 区画にいる 樹木個体にしか「眼」を与えていないからである． つまり 0.25 ha × 171 個．
• ちなみにひとつの計算区画 (「眼」をもつ個体数 900-1000 ぐらい) をクリアーするのに要する時間は …… BookPC (PentiumIII 733MHz) で 6.4-6.5 秒ぐらいか． 開空度 (ごとき指標) 計算するより， 2500-3000 樹木個体ぶんの データ読み込みのほうに時間かかってるかんぢだ． えーと， ならば 50ha 全体だと …… 約 18 分， というところか．
• この計算時間って， 下手な新ライブラリ作ったりすると 数倍にハネ上がりそう …… やる気失せるな．
• ともかく， いまやることは終わりの見えない新ライブラリじゃなくて， 171 個の計算区画を順番に呼び出していって， pasoh_simulator に投げ込んでいき， 結果をまとめてくれる「進行とりまとめ役」の作成だな．
• おっと考えてみれば， 25m 区画を 16 個つなぎあわせて 1ha 計算区画を 生成するところも時間かかりそうだ …… 測定してみると， だいたい 2.2-2.3 秒ぐらいか． とすると 1 計算区画 8.7 秒ぐらいかかるだろうから， 全体で 24 分ぐらいね． さすがに 50 ha 調査区はデカい．
• 発注者に報告メイル書いて撤退． 1803 研究所発． 1810 都営浅草線大門発． 京急平和島まで 360 円． 1840 帰宅． 途中のローソンでクロネコヤマト宅急便の 荷札を二枚もらう．
• ここで晩飯食うと昏睡するだろうから， 晩飯の時刻を遅らせる実験をやってみる． 紅茶を淹れてビスケット少々かじり空腹をしのぐ．
• これでもやはり眠くなったんだけど， 1 時間ぐらい寝たり起きたりを繰り返してたら， 眠気がさめた．
• Yahoo! オークション落札者たちに「商品」発送に出る． これで作業は全部片付いたか． あとは先方に届いてからだが …… 買物などして帰宅．
• 遅くに晩飯食っても眠くならない． まことに人体とは不思議である． うう， 何時になったら眠くなるんだろうか．
• 今日の食卓
  • 朝 (0650): 米 0.7 合． 乾物味噌汁の残り．
  • 昼 (1215): 弁当． 米 0.7 合． 豚肉麦酒煮． 冷えてもうまいが， 弁当の菜にするには工夫がいるな． タッパーウェアの中で凝結してる ラードのかたまりってのは ……
  • 晩 (2350): 米 0.3 合． 豚肉麦酒煮を分派して雑炊．