MORIZON NEXT 操作マニュアル
収益性と災害リスクから森林をゾーニングする QGIS プラグイン
はじめに
MORIZON NEXT は、林野庁の森林ゾーニング支援ツール「もりぞん(MORIZON)」を改良したプラグインです。ゾーニングの考え方と計算方法は原版をそのまま受け継いでいます。改良したのは、解析に必要なオープンデータをプラグインが自動で取得するようにした点です。
原版では、解析に入る前に、DEM の変換・結合、建物や道路のダウンロードと座標変換、作業システム表の Excel 編集などをすべて手作業で行う必要がありました。MORIZON NEXT では、QGIS の地図で範囲を表示してボタンを押せば、これらのデータがそろいます。
このマニュアルは、原版の「収益性と災害リスクを考慮した森林ゾーニングの手引き」(令和8年3月、以下「手引」)を置き換えるものではありません。操作に必要なことと、判断に必要な考え方に絞って書いています。詳しい背景を知りたいところには、手引の該当ページを p.00 の形で示しました。
最初の章はデータの準備です。データごとに「なぜこのデータが必要なのか」を説明しながら進めるので、読み終えるころには、この解析が何を計算しているのかもつかめるはずです。
1. データの準備
1.1 何のためのデータか
最終的に作るのは、森林を次の4つに色分けしたゾーニング図です。
| 災害リスクが低い | 災害リスクが高い | |
|---|---|---|
| 収益性が高い | 林業に向いている | 収益は見込めるが、施業に配慮が必要 |
| 収益性が低い | 災害の心配は少ないが、収益は見込みにくい | 林業生産には向かない |
「収益性」と「災害リスク」は、それぞれ3つの要素の点数(1〜3点)を合計して求めます。要素は全部で6つです。
| 軸 | 要素 | 何を表すか | 元になるデータ |
|---|---|---|---|
| 収益性 | 地位 | その土地で木がどれだけよく育つか | 地位データ(NPP・SRAD・VTEX) |
| 収益性 | 集材作業効率 | その地形で使える機械と作業のしやすさ | DEM + 作業システム |
| 収益性 | 地利 | 道からの近さ(運び出しやすさ) | 道路 + DEM |
| 災害リスク | 傾斜 | 斜面の急さ | DEM |
| 災害リスク | 地形の複雑さ | 尾根や谷の入り組み具合(崩れやすさの目安) | DEM |
| 災害リスク | 保全対象を含む流域 | 崩れたとき、下流に人家があるか | DEM + 建物 |
ここから2つのことがわかります。
- 用意するデータは5種類だけです。DEM、地位データ、建物、道路、作業システムです。
- DEM は6つの要素すべてに関わります。DEM の質と範囲が、結果全体の質を決めます。
1.2 準備の前に:範囲と保存場所を決める
QGIS プロジェクトを先に保存する
自動取得したデータは、すべて QGIS プロジェクトファイルと同じフォルダの中の morizon_next フォルダの、プロジェクトのファイル名のフォルダ(プロジェクト内の作業フォルダ)に保存されます。同じフォルダに複数の QGIS プロジェクトを置いても、プロジェクトごとに分かれます。そのため、データを取得する前にプロジェクトを一度保存してください。保存していないと、取得を始めたときに案内が出ます。
ここから「1つの QGIS プロジェクトで扱う解析は1つです」の手前までは、自動で作られるフォルダの参考です。読み飛ばして、必要なときに参照してください。
取得が終わると、フォルダの中は次のようになります。
(プロジェクトフォルダ)/
├─ ○○.qgz … QGIS プロジェクト
└─ morizon_next/
├─ ○○/ … プロジェクト内の作業フォルダ(プロジェクトのファイル名)
│ ├─ DATA/ … 入力データ
│ │ ├─ DEM/ dem_fetched.tif
│ │ ├─ SiteIndex/
│ │ │ ├─ NPP/ npp_clipped.tif
│ │ │ ├─ SRAD/ srad_clipped.tif
│ │ │ └─ VTEX/ vtex_clipped.tif
│ │ ├─ TATEMONO/ building.shp
│ │ ├─ ROAD/ road_edge.shp
│ │ └─ SAGYO-SYSTEM_CSV/ costcsv.csv
│ ├─ YOUSO/ … 要素計算の出力
│ ├─ ZONING/ … 収益性・災害リスク・ゾーニング図の出力
│ └─ AGGREGATE/ … ゾーン統計量の出力
├─ shared/ … 再取得を省くための控え(同じフォルダのプロジェクトで共有。触らなくてよい)
└─ archive/ … 保存ファイル(ZIP)
作業フォルダの構成は原版の「ZoningKit」と同じです(AGGREGATE だけは MORIZON NEXT で追加したフォルダです)。画面右上の「保存ファイル出力」で、入力と出力をこの構成のまま ZIP に保存できます(保存先は morizon_next/archive)。保存した ZIP は要素計算タブの「保存ファイルを読み込む」で読み込むと、各タブの入力欄と出力レイヤーが保存したときの状態に戻ります。QGIS プロジェクトを保存して開き直したときは、「再開しますか?」と確認して、前回の作業フォルダの続きから作業できます。
範囲は「QGIS の地図に表示している範囲」
「DEMブラウザから開始する」を押すと、その時点で QGIS の地図に表示されている範囲がデータの取得範囲になります。計画図や市町村界などのレイヤを表示し、対象地が画面に収まるように拡大・縮小してから始めてください。
- 起伏量と地形の複雑さは、周囲約 490m 四方の標高を使います。
- 保全対象を含む流域は、水が流れ下る先まで追いかけて計算します。
- 地利は、範囲の外にある道のほうが近いこともあります。
広さと解像度の目安
手引が想定しているのは、市町村全域を 10m 解像度で解析する使い方です。
- 5m の DEM は、解析のときに自動で 10m に変換されます。
- 1m の DEM は、約 2,500 万画素(5km 四方程度)を超えると 10m に変換されます。1m のまま解析できるのは、施業団地程度の狭い範囲に限られます。
1.3 DEM(標高)
- 傾斜:急な斜面ほど崩れやすく、作業もしにくくなります。
- 地形の複雑さ:尾根と谷が細かく入り組んだ場所は、侵食が進んでいて崩れやすい傾向があります。
- 起伏量:周囲の最高点と最低点の差です。傾斜と組み合わせて、使える作業システムを決めます。
- 流域:雨水や崩れた土砂がどこへ流れ下るかを求めます。
- QGIS プロジェクトを保存し、対象地を余白付きで地図に表示します。
- 「要素計算」タブの 「DEMブラウザから開始する」を押します。
- 取得元を選んで OK を押します。通常は「国土地理院 DEM1A/5A/10B」を選んでください。
- 取得が終わると、解像度と座標系が表示され、DEM 欄が「設定済」になります。そのまま地位データの取得に進みます(1.4)。
解像度はプラグインが選びます
国土地理院の DEM には 1m(航空レーザ測量、一部地域のみ)、5m、10m(全国)があります。プラグインは取得前に、表示範囲の全体を同じ解像度でカバーできるかを調べます。
- 全体をカバーできる解像度が1つだけなら、その解像度で取得します。
- 複数あれば、どれを使うかを選べます。市町村規模なら 5m か 10m で十分です。
座標系も自動で揃えます
取得した DEM は、その地域の平面直角座標系(例:静岡県なら第8系)に変換してから保存されます。
1.4 地位データ(NPP・SRAD・VTEX)
ただし、全域の樹高を実際に測ることはできません。そこで、全国の森林調査データから作られた推定式を使い、次の3つのデータからスギ・ヒノキ・カラマツの地位指数を計算します。
- NPP(純一次生産力):気候から見た、植物が育つ力
- SRAD(日射係数):斜面の向きや地形による日当たり
- VTEX(凹凸度):谷か尾根か。谷は水が集まり、木がよく育つ傾向があります
- 取得元は、林野庁がもりぞん用に公開している「ZoningKit」(G空間情報センター)です。登録は不要です。
- 表示範囲の中心がどの座標系の地域にあるかを判定し、必要な部分だけを取り出します。
- 一度取得した地域のデータは morizon_next/shared に控えが残り、同じフォルダのプロジェクトで取り直すときにはダウンロードを省きます。
1.5 建物(保全対象)
そこでこの解析では、DEM から流域(雨水が同じところへ集まる範囲)を求め、その中に建物があるかどうかで点数を付けます。これが「保全対象を含む流域」です。
- 地位データの取得が終わると、「続けて建物ポリゴン・道路縁データも自動取得しますか?」と聞かれます。「はい」を選びます。
- 基盤地図情報ダウンロードサービスのログイン画面が出るので、ID とパスワードを入力します。
- 国土地理院の基盤地図情報「建築物の外周線」を取得し、DEM と同じ座標系に変換して保存します。
地域の保全対象を反映する場合は、独自の建物データの作成をおすすめします。作業フォルダの建物データに描き足して計算し直す手順は「2.2 高度な使い方」で説明します。
公式の手引:収益性と災害リスクを考慮した森林ゾーニングの手引き(令和8年3月) 参照ページ:p.23〜24/p.75
- 対象にしているのは建物だけです。道路、鉄道、農地、取水施設など、ほかの保全対象は含まれません。
- 「流域に建物があるか」を見ているだけで、土石流が実際に届くかどうかは計算していません。
- 流域は約 5ha 以上のまとまりで区切っています。
1.6 道路
実際の林道・作業道を反映する場合は、独自の道路データの作成をおすすめします。作業フォルダの道路データに林道・作業道を描き足して計算し直す手順は「2.2 高度な使い方」で説明します。
公式の手引:収益性と災害リスクを考慮した森林ゾーニングの手引き(令和8年3月) 参照ページ:p.42〜43/p.71〜72
- 林道や森林作業道は十分に含まれていません。山の中では、実際より道から遠く計算されがちです。手引でも、林道データは自分で用意するデータとされています。
- 逆に、高速道路や国道など、搬出には使えない道も含まれています。
- 距離は地図上の直線距離です。傾斜や、川を渡れないことは考慮していません。
1.7 作業システム
この対応表は、5種類のデータの中で唯一、地域の人が判断して決めるものです。地域で実際に使われている機械と合っていなければ、集材作業効率は現実とずれます。
- 色分布パネル(右):横軸が傾斜、縦軸が起伏量です。色が、その地形に割り当てる機械を表します。デフォルトは、原版の Excel ひな形の「入力例」シートと同じ内容です。「デフォルトに戻す」でいつでもこの内容に戻せます。
- 機材名の一覧(左):色と機械の対応です。数字が大きい機械ほど作業効率が高い扱いになります。色の四角をクリックすると「保有していない機械」になり、その機械が割り当てられていた地形は「該当なし」(最も効率が低い)になります。
- マージン:プラス側にすると、斜面を実際より最大3度緩いものとみなし、同じ地形でもより効率のよい機械を割り当てます(強気)。マイナス側はその逆(慎重)です。緩斜面と、機械の限界を超える急斜面の割り当ては動きません。
- 「CSVとして反映する」を押すと、作業フォルダの DATA/SAGYO-SYSTEM_CSV/costcsv.csv に保存され、作業システムCSV 欄に自動で設定されます。
地域の作業システムをマスごとに細かく決めたい場合は、原版と同じ Excel のひな形も使えます。
- 「Excelシートを開く」を押すと、ひな形(集材作業効率の設定220113.xlsx)が、作業フォルダの DATA/SAGYO-SYSTEM_CSV に置かれて開きます。すでに置いてある場合は、そのファイルを開きます。
- ひな形の説明どおりに「入力」シートを編集し、上書き保存します。CSV で保存し直す必要はありません。
- 保存すると、その内容がパネルと一覧に表示されます。
- 表示を確かめてから「CSVとして反映する」を押します。
素材生産事業体の意見を聞きながら決めるのが理想です。まずは初期値のまま解析を一度通し、結果を見てから調整しても構いません。
1.8 自分でデータを用意する場合
航空レーザの DEM、独自の地位指数、林道データなど、手元により良いデータがあれば、各入力欄の右端「…」から個別に指定できます。
守ってほしいこと
- DEM は平面直角座標系(単位 m)の GeoTIFF にしてください。ほかのデータはすべて DEM の座標系に合わせて変換されるため、DEM の座標系が解析全体の基準になります。緯度経度や Web メルカトルの DEM では、傾斜・距離・面積が正しく計算されません。
- DEM に穴を開けないでください。民有林だけ、5条森林だけといった切り抜き方をした DEM は縁が多くなり、流域などが正しく計算できません。対象地より広い範囲を、隙間なく覆ってください。
- 海や湖の標高値は取り除いてください。水面の平らな値が、傾斜などの統計を歪めます。
- 建物と道路は Shapefile(.shp)で指定します。座標系が DEM と違っても自動で変換されます。
林道データを道路に加える
既設路網ライン欄に指定できるのは1ファイルだけです。基盤地図情報の道路縁と林道データを両方使うときは、先に QGIS で1つにまとめます。
- QGIS メニュー「ベクタ」→「データ管理ツール」→「ベクタレイヤのマージ」を開きます。
- road_edge.shp と林道データを選び、出力先の座標系には DEM と同じものを指定して、Shapefile として保存します。
- 保存したファイルを既設路網ライン欄に指定します。
フォルダ選択から開始する
原版の ZoningKit のフォルダ(またはその中の DATA フォルダ)や、MORIZON NEXT の作業フォルダ(morizon_next/○○)を選ぶと、DATA の中のサブフォルダ(DEM、SiteIndex/NPP など)から入力ファイルをまとめて設定します。ファイルはコピーせず、選んだ場所のものをそのまま使い、出力もそのフォルダの YOUSO・ZONING・AGGREGATE に書き込みます(作業フォルダ:外部)。開始する前に、使うファイルの一覧を確認できます。選んだフォルダに出力(YOUSO・ZONING・AGGREGATE のファイル)がある場合は、途中まで進んだフォルダとして読み込み、出力をレイヤーに反映します。入力欄の「…」でファイルを個別に選んだ場合は、作業フォルダの DATA にコピーして使います。作業フォルダは QGIS プロジェクトに記録されるので、プロジェクトを保存して開き直すと「再開しますか?」と確認し、同じフォルダで続きから作業できます(「設定をクリアする」では初期状態、作業フォルダ:なし に戻ります)。
1.9 実行前の確認
各入力欄の右端に状態が表示されます。
- 設定済:ファイルが指定されています。
- 未設定:計算にチェックが入った要素に必要なのに、まだ指定されていません。
- 任意:チェックが外れた要素だけに使うので、なくても実行できます。
「計算する要素を選択」のチェックは、必要なデータが揃った要素だけ選べるようになり、揃った時点で自動的にチェックが入ります。
- DEM は対象地より広く、穴がないか
- 地位データ・建物・道路が「設定済」になっているか
- 作業システムは地域の機械に合わせたか(初期値のままでもよいが、そのことを覚えておく)
- 林道データを加えるかどうか決めたか
- 作業フォルダが「なし」になっていないか(出力先は作業フォルダの中に決まります)
すべて揃うと「実行」ボタンが押せるようになります。操作の続きは「2.1 要素計算」で説明します。
1.10 出典と利用条件
自動取得するデータの出典は次のとおりです。成果を資料にするときは出典を明記してください。
| データ | 出典 |
|---|---|
| DEM | 国土地理院(地理院タイル 標高タイル DEM1A / DEM5A / DEM10B) |
| DEM(静岡県内で選択時) | 静岡県「VIRTUAL SHIZUOKA」航空レーザ測量データ LP/Grid(G空間情報センター、CC BY 4.0) |
| DEM(長野県内で選択時) | 長野県建設部砂防課 航空レーザ測量成果 0.5mメッシュDEM(令和3〜7年度)、長野県林務部 航空レーザ測量成果 0.5mメッシュDEM(2013〜2014年度)(いずれも G空間情報センター) |
| 建物・道路縁 | 国土地理院 基盤地図情報(建築物の外周線、道路縁) |
| NPP・SRAD・VTEX | 林野庁「森林ゾーニング支援ツール もりぞん」ZoningKit(G空間情報センター)。地位指数の推定式は光田靖・北原文章(2017)による |
| DEM(予備) | AWS Terrain Tiles(Mapzen Terrarium) |
ゾーニング図を印刷して配布する場合など、使い方によっては国土地理院への利用手続きが必要になることがあります。事前に国土地理院の案内で確認することをお勧めします。
2. 操作
解析は、タブの並びのとおり 要素計算 → スコアリング → ゾーニング → ゾーン統計量 → 印刷 の順に進めます。前の工程の結果を次の工程が使うので、途中の工程をやり直したときは、その先の工程もやり直します。
2.1 要素計算
- 収益性地位(スギ・ヒノキ・カラマツ)、集材作業効率、地利
- 災害リスク地形の複雑さ、傾斜、保全対象を含む流域
- 「計算する要素を選択」で、計算する要素にチェックを入れます。必要なデータが揃った要素だけ選べます(1.9)。最初は全部を計算します。
- 「乗算出力を適用する」は、通常はオンのままにします。出力のレイヤーを乗算で重ねるので、上のレイヤーが下を隠さず、要素を重ねて見比べられます。表示が重いときはオフにできます。
- 「実行」を押します。処理の窓に進み具合が表示されます。「中断」を押すと、それまでに計算できた要素だけをレイヤーに出して止まります。
出力
計算結果は、作業フォルダの YOUSO に保存され、レイヤーパネルの「Morizon Next」グループにまとめて表示されます。
- 「収益性」と「災害リスク」の下に、要素ごとのグループができます。各グループには、計算した値そのもの(生データ)と、点数に分けたもの(名前の末尾が「[スコアリング]」)の2つのレイヤーが入ります。
- 最初は各グループの表示がオフです。見たいグループのチェックを入れると表示されます。地利だけは、最初に点数の方が表示されます。
- 「収益性」の要素どうしは、1つを表示すると他が消えます(同時に1つだけ表示)。「災害リスク」の要素は、傾斜の上に地形の複雑さや保全対象を含む流域を重ねて見られるよう、同時に表示できます。
作業フォルダの行のボタン
「出力先フォルダ」の上の行に、今の作業フォルダ(なし・プロジェクト内・外部)が表示されます。その右のボタンは次のとおりです。
- 道路と建物を出力:取得・読み込みした道路縁と建築物を、「Morizon Next」グループの「災害リスク」の下の「道路・建物」グループに表示します。要素計算の入力に使っている道路と建物を、地図で確かめるときに使います。
- 作業フォルダを開く:作業フォルダを、パソコンのファイル管理の画面で開きます。
2.2 高度な使い方:道路・建物を編集して計算し直す
- 要素計算を一度実行し、結果を確かめます。
- 「作業フォルダを開く」で作業フォルダを開き、DATA/ROAD の道路のファイル(自動取得なら road_edge.shp)を QGIS に読み込みます。建物なら DATA/TATEMONO のファイルです。
- 読み込んだレイヤーを編集します。計画路網を線で描き足す、使えない道を削る、などです。編集が終わったら保存して編集を終えます。
- 「道路と建物を出力」を押します。編集が見つかると、表示が新しい内容に更新され、「要素計算をもう一度実行してください」と案内されます。
- 要素計算をもう一度実行します。道路を編集したなら「地利」、建物を編集したなら「保全対象を含む流域」だけを選べば足ります。
- 要素が変わったので、スコアリング・ゾーニング・ゾーン統計量も、この順にやり直します。
- 要素計算が使うのは、作業フォルダの DATA/ROAD・DATA/TATEMONO にあるファイルです。「道路・建物」グループに表示されるレイヤーは確認用の表示なので、こちらを編集しても計算には反映されません。
- ファイルの名前と場所は変えずに、上書きで保存してください。別の名前で保存した場合は、要素計算タブの入力欄でそのファイルを指定し直します。
- 編集中のまま保存していないと、変更はファイルに書き込まれていないので反映されません。
- 前回と同じ範囲なら、「編集しているデータを使いますか?」と聞かれます。「はい」なら編集したデータをそのまま使い、取得しません。「いいえ」なら、編集したデータを同じフォルダの backup フォルダに日時付きの名前で残してから、取得し直します。
- 前回と違う範囲なら、編集したデータを backup フォルダに残してから、新しい範囲で取得し直します。残した場所は画面で知らせます。
- 編集していなければ、前回と同じ範囲で取得し直したときは、取得済みのデータをそのまま使います(まとめ直さないので早く終わります)。
2.3 スコアリング以降
準備中:スコアリング・ゾーニング・ゾーン統計量・印刷の操作
3. パラメータとしきい値
準備中
4. 凡例の見方
準備中
5. 結果の使い方と注意点
準備中
付録 手引との対応表
準備中
不具合報告
GitHub Issues で受け付けています。個人で開発しているため、すべての報告に対応できるとは限りません。
不具合は、使う環境や設定によって起きることが多く、同じ環境を用意しないと確かめられません。対応を検討するには、次の情報が必要です。
- QGIS のバージョンと OS
- MORIZON NEXT のバージョン
- 操作の手順(どの画面で何をしたか)
- エラーメッセージ(QGIS の「ログメッセージ」パネルの内容)
- 使ったデータの種類(自動取得か、手持ちのデータか)
- 追加している QGIS プラグインの一覧
- 追加した Python ライブラリ。自分で入れたものと、プラグインの機能で自動的に入ったものの両方を書いてください
- 上記の追加プラグインとライブラリを外し、動作要件に書いたものだけの環境で試した結果
頂いた情報で再現しない、判断ができない報告には対応ができないことがあります。