探す範囲
文献：Mitsuda, Y., Ito, S.（2015）J. For. Plann. 20: 45–51。J-STAGE で無料公開：J-STAGE
1. 定義
Introduction（p.46）より：ラスター DEM で地形因子を計算するアルゴリズムは、対象のセルの現象に影響する地形をどこまで探すかを決める任意のパラメータ（半径またはウィンドウの大きさ）を必要とする。そのパラメータで決まる範囲（円または矩形）を search range と呼ぶ。
つまり、あるセルの値を計算するときに、周りの地形を何 m まで見るかという距離。セルの細かさ（解像度）と並ぶ、地形解析の「スケール」のもう一つの軸と位置づけられている（Wilson et al. 2000 を引用）。
2. 3つの要素ごとの使われ方（p.47）
| 要素 | 探す範囲で決まること |
|---|---|
| SRI（＝SRAD） | 太陽がそのセルから見て周りの稜線より上にあるか（日陰の判定）を、範囲内の標高で判断する |
| UCA（集水面積） | 範囲内で、そのセルに水が集まってくる面積を求める（ZoningKit には含まれない） |
| VTEX | 8方向それぞれで範囲内の最も低い点を探し、そこへの俯角を合計する |
試した範囲：100・250・500・1000 m。DEM：国土地理院50m（GSI50）、国土地理院10m（GSI10）、等高線のデジタル化（CNT）。解像度：10・12.5・25・50 m。対象地：宮崎大学田野演習林、18地点（伐倒・樹幹解析による地位指数）。
3. 結果
- VTEX：範囲が広いほど、地位指数との相関がどの DEM でも悪くなった。適した範囲は「明らかに100m」。VTEX は「尾根らしさ」を表し、範囲が広いと尾根の向こうの谷底まで拾って尾根らしさを過大に見積もるため
- SRI・UCA：範囲が広がると相関はすぐ頭打ちになり、範囲の影響ははっきりしなかった
最良モデルで採用された組み合わせ（Table 2）
| モデル | SRI | UCA | VTEX |
|---|---|---|---|
| 全データ | 1000 m（CNT-10） | 100 m（CNT-12.5） | 100 m（GSI50-50） |
| 国土地理院 DEM（GSI） | 1000 m（GSI50-50） | 100 m（GSI10-10） | 100 m（GSI50-50） |
| 国土地理院50m のみ（GSI50） | 1000 m（50 m） | — | 100 m（50 m） |
Discussion の本文には「UCA・SRI・VTEX の範囲はそれぞれ 1000 m、100 か 250 m、100 m」とあり、表と SRI・UCA の順番が逆になっている。表のほうが正しいと考える。
DEM の細かさとの相性
- UCA は、情報の多い DEM（CNT、GSI10）の細かい解像度（10〜12.5 m）が向く。地下水の流れという細かい現象のため
- SRI は、情報の少ない DEM（GSI50）の粗い解像度（50 m）が向く。日射は広い範囲の現象で、細かい DEM は「細かすぎて」かえって表しにくい
- 粗い元データを細かくしたもの（GSI50-10 など）や、細かい元データを粗くしたもの（CNT-50 など）は役に立たなかった。元データの質に合った解像度を選ぶ必要がある
4. 計画への意味
- ZoningKit の設定の手がかり：ZoningKit は国土地理院50m DEM から作った50mのデータなので、この論文の「50m のみ」の最良モデルと同じく SRAD 1000 m・VTEX 100 m で作られた可能性がある（推測）。手元の DEM で VTEX を計算して ZoningKit の値と比べれば確かめられる
- 外国での設定値：範囲の効き方は対象地の地形（試験地は起伏の激しい土地）に左右されると論文は書いている。外国では探す範囲をプラグインで変えられる設定にし、既定値を VTEX 100 m・SRAD 1000 m にするのが自然
- DEM の細かさ：外国で LiDAR 1 m の DEM を使う場合も、SRAD は粗くしてから計算する必要があるかもしれない。設計に効く点
5. 残ること
| # | 項目 | 状態 |
|---|---|---|
| 1 | ZoningKit を作ったときの探す範囲（SRAD 1000 m・VTEX 100 m という推測の検証） | 手がかりあり |
| 2 | SRAD で、どの日時の太陽の位置を合計したか | 未 |
| 3 | NPP をどの樹種のモデルで計算したか | 未 |