forestfire.html

状　態	記　号	画　像
空き地	0
木が生えている	1
木が燃えている	2

[ 1 ]　初期状態　：　木の初期密度（セル当たりの初期存在確率） r と，初期に木が燃えている確率 k を入力で与え， 0，1，2の個数がその割合になるようにして，nxn 個のセルにランダムに配置する．

[ 2 ]　状態遷移： von Neumann近傍を使い，周期的境界条件を採用．４個の隣接セルの状態に依存して，着目セルの次の時間ステップでの状態が決まる．

状態遷移（延焼）のルール　：

空き地のセルからは，確率 p で木が生える。（生成：0 から 1 へ）
着目する木の風下の隣接セルに燃えている木があれば，着目する木は確率 g1 で燃え始める．（風上への延焼：1 から 2 へ）
着目する木の風上の隣接セルに燃えている木があれば，着目する木は確率 g2 で燃え始める．（風下への延焼：1 から 2 へ）
着目する木の風とは垂直な方向の隣接セルに燃えている木があれば，着目する木は確率 (g1+g2)/2 で燃え始める．（風と垂直な方向への延焼：1 から 2 へ）
４個の隣接セルに燃えている木が全く存在しなければ，その木は確率 h×g1 で燃え始める．(自然発火：1 から 2 へ）
燃えている木は確率 q で燃え尽き，空き地になる．（燃え尽き：2 から 0 へ）

[ 3 ]　森林延焼の時間的空間的変化　：　このルールを nxn 個のすべてのセルに同時に適用して，１時間ステップ後の森林の状態が決まる．新しい森林の状態に同じルールを適用して，次の１時間ステップ後の森林の状態が決まる．この過程を繰り返して時間と共に森林が延焼していく様子が実現される．

n	森林領域の大きさ（n×n個のセルから成る正方格子で表す）
r	木の初期密度（セル当たりの初期存在確率）
k	初期に木が燃えている確率
p	空き地から木が１時間ステップで生える確率
q	燃えている木が１時間ステップで燃え尽きる確率
h	落雷等の確率（自然発火の確率は h×g1 とする）
g1	風上へ延焼する確率
g2	風下へ延焼する確率
w	風向（w = 0：無風，w = 1：北風，w = 2：東風，w = 3：南風，w = 4：西風）
t	シミュレーションの時間ステップ数

森林領域（200 x 200 個のセル）に木が 30% 生え，そのうち0.1% が燃えているとして，ランダムに 0, 1, 2 を配置すると，下の左に示す画像を得る（初期状態　t = 0）． ♣ この画像上にマウスポインタを置くとアニメーションが始まる．はずすと初期画像にもどる．自然発火の確率は 0.0 とする．

森林領域（200 x 200 個のセル）に木が 30% 生え，初めに燃えている木はないとして，ランダムに 0, 1 を配置すると，下の左に示す画像を得る（初期状態　t = 0）． ♣ この画像上にマウスポインタを置くとアニメーションが始まる．はずすと初期画像にもどる．自然発火の確率を 0.0001 とする．
[ 2 ] の場合と初期状態は同じであるが，空き地から木が生え難く，また木が燃え尽き難いとした場合である．

森林領域（200 x 200 個のセル）に木が 30% 生え，そのうち5% が燃えているとして，ランダムに 0, 1, 2 を配置すると，下の左に示す画像を得る（初期状態　t = 0）． ♣ この画像上にマウスポインタを置くとアニメーションが始まる．はずすと初期画像にもどる．自然発火の確率は 0.0 とし，北風が吹いているとする．