ぎょーむ日誌 2013-10-09

2013 年 10 月 09 日 (水)

• 0430 とら泣きおき …… 30 分ほどかけて再ねかしつけ．
• 0550 起床 …… とら坊はまだぐーぐー寝ている． コーヒー． 朝飯・弁当の準備． とらめし． 朝飯． 0820 自宅発． 曇天． 0830 保育所着． 0900 研究室着．

[いきなり破断]

自転車が走行中にいきなりこわれた． head tube と seat tube をナナメにつないでいたパイプがぱきっと …… まあ， たしかにチカラがかかりそうなところに， こういう弱そうな部品を使っているというのが …… ともあれ走行不可能状態． とりあえず， 小野サイクル に交換用の部品が手に入るかどうか尋ねるメイルをだしてみた．

• また R の library(rgl) 調査．
• 昼飯．
• ヒノキデータをひねくってみる．
• だめだ， 調査地の三次元座標がまったくわからん …… 隅田さんの論文では傾斜 37 度とか書いてあるけれど， それとはどうやっても矛盾してしまう． また後日に問い合わせてみますか ……
• ということで， かなりらんぼうながら， 傾斜はゼロということにする． とりあえず， の措置 ……
• library(rgl) を使った三次元作図プログラミングは初めてなのだが …… いろいろ工夫してみる．
• とりあえず， ヒノキの幹については， こんな図示かな．

• 根本が太く， 先端ほど細くなるのは観測データを反映しているため …… つまりこのデータは 1977 年時点における， 142 本のヒノキについて高さ 1m おきに太さを測定したもの． このデータが 21 年ぶんある．
• 幹については， library(rgl) の shade3d と cylinder3d() をこんなふうに組み合わせて作図している．

  shade3d(
  	cylinder3d(center = xzy, radius = radius),
  	color = "#a04040" # alpha = 0.5 などと透過 (semi-transparent) 指定可能
  )
  

  円筒といっても下の図 (上下転倒した視点から見ている図) でわかるように， ホントの円筒ではなくポリゴンで近似していることがわかる． 円周を分割するポリゴン数は cylinder3d() のオプション sides で指定できる． 下の図の場合は， デフォルトである sides = 8．

• おっと， 幹がむやみに太い理由は， まちがって幹の半径ではなく周囲長を使って作図してしまったため ……
• 1700 研究室発． 買いもの． 1730 帰宅． 晩飯の準備． とらめし． 晩飯．
• また買いものに出たり …… しかしばてていたので， いろいろとしくじる．
• [今日の運動]
  • とらちゃんえくささいず．
• [今日の食卓]
  • 朝 (0730): 米 0.7 合． ネギ・豆腐の味噌汁． キュウリ・モズク酢．
  • 昼 (1200): 研究室． 米 0.7 合． イカあげ． ソーセイジ． カボチャ． ミニトマト．
  • 晩 (1830): 米 0.7 合． ネギ・豆腐の味噌汁． ピーマン・タマネギ・ラムの炒めもの． ホウレンソウごまあえ．