(57)【要約】【課題】立木密度の違いや樹高及び樹冠の大小が混交する森林においても、森林資源情報を高精度に算定する。【解決手段】３次元点群データと地理的情報とに基づいて、調査対象森林域画像データを作成する調査対象森林域画像データ作成処理（ステップＳ３０）と、調査対象森林域に存在する各樹木の樹冠の高さを樹冠高として表す樹冠高画像データの作成を行う樹冠高画像データ作成処理ステップ（ステップＳ４０）と、樹冠高画像データから得られる樹冠高のうち、最大樹冠高が存在する階層を第１階層として、最大樹冠高を起点に所定長さの階層幅を有する階層を前記第１階層から地表面に向かって複数設定したときに、第１階層を含む各階層に樹頂点が存在する樹冠の樹頂点を、階層ごとにそれぞれ抽出する階層別樹頂点抽出処理ステップ（ステップＳ５０）と、を有する。【選択図】図１
2019185449.tif 000002
【特許請求の範囲】
【請求項１】
調査対象森林域に関する森林資源情報を作成する森林資源情報算定方法であって、
前記調査対象森林域を含む地域に対応する３次元点群データと地理的情報とに基づいて調査対象森林域画像データを作成する調査対象森林域画像データ作成処理ステップと、
前記調査対象森林域画像データに基づいて、前記調査対象森林域に存在する各樹木の樹冠の高さを樹冠高として表す樹冠高画像データの作成を行う樹冠高画像データ作成処理ステップと、
前記樹冠高画像データから得られる樹冠高のうち、最大樹冠高が存在する階層を第１階層として、前記最大樹冠高を起点に所定長さの階層幅を有する階層を前記第１階層から地表面に向かって複数設定したときに、前記第１階層を含む各階層に樹頂点が存在する樹冠の樹頂点を、前記階層ごとにそれぞれ抽出する階層別樹頂点抽出処理ステップと、
を有することを特徴とする森林資源情報算定方法。
【請求項２】
請求項１に記載の森林資源情報算定方法において、
前記樹冠高画像データ作成処理ステップは、前記調査対象森林域画像データに基づいて、メッシュ化されたデジタル表層モデルデータとデジタル標高モデルデータとを作成し、メッシュごとに前記デジタル表層モデルデータと前記デジタル標高モデルデータとの差分を取ることによって、メッシュごとに樹冠高が求められた樹冠高画像データを作成することを特徴とする森林資源情報算定方法。
【請求項３】
請求項２に記載の森林資源情報算定方法において、
前記階層別樹頂点抽出処理ステップには、
前記メッシュごとに樹冠高が求められた樹冠高画像データを用いて、前記第１階層を含む各階層に存在する個々の樹冠に対応する階層別樹冠高点群データを、前記階層ごとに作成する階層別樹冠高点群データ作成処理ステップと、
前記階層別樹冠高点群データ作成処理ステップによって作成された前記個々の樹冠に対応する階層別樹冠高点群データから、前記個々の樹冠に対応する階層別樹冠高点群データの最大値を前記階層ごとに算定して、算定された最大値を当該階層別樹冠高点群データに対応する樹冠の樹頂点とする樹頂点算定処理ステップと、
が含まれていることを特徴とする森林資源情報算定方法。
【請求項４】
請求項３に記載の森林資源情報算定方法において、
前記樹頂点算定処理ステップは、前記最大値を算定する際においては、前記個々の樹冠に対応する前記階層別樹冠高点群データのうち、前記階層別樹冠高点群データが既に作成されている樹冠と同じ樹冠に対応する階層別樹冠高点群データは除外して、前記最大値を算定することを特徴とする森林資源情報算定方法。
【請求項５】
請求項３又は４に記載の森林資源情報算定方法において、
前記階層別樹頂点抽出処理ステップには、
前記樹冠高画像データの各メッシュをそれぞれ処理対象メッシュとし、当該処理対象メッシュに平均化フィルターを重ね合わせ、当該平均化フィルターに含まれる各メッシュの樹冠高の平均値を当該処理対象メッシュの樹冠高とする処理を、各処理対象メッシュにおいて順次行うことによって、ノイズ除去済みの樹冠高画像データを作成する平均化処理ステップがさらに含まれており、
前記階層別樹冠高点群データ作成処理ステップは、前記平均化処理ステップによって作成されたノイズ除去済みの樹冠高画像データを、前記メッシュごとに樹冠高が求められた樹冠高画像データとして用いて前記階層別樹冠高点群データを作成することを特徴とする森林資源情報算定方法。
【請求項６】
請求項１〜５のいずれかに記載の森林資源情報算定方法において、
前記階層幅は、０．２ｍ〜２．０ｍの範囲の長さに設定されていることを特徴とする森林資源情報算定方法。
【請求項７】
請求項１〜６のいずれかに記載の森林資源情報算定方法において、
前記３次元点群データは、上空から撮影箇所を重複させ、かつ、ずらしながら撮影して得られた撮影画像データから、３次元形状復元技術であるＳＦＭ(Structure from Motion)によって作成されていることを特徴とする森林資源情報算定方法。
【請求項８】
請求項１〜６のいずれかに記載の森林資源情報算定方法において、
前記３次元点群データは、上空からレーザー光を照射して得られたレーザー計測データに基づいて作成されていることを特徴とする森林資源情報算定方法。
【請求項９】
請求項１〜８のいずれかに記載の森林資源情報算定方法において、
前記階層ごとにそれぞれ抽出された樹冠の樹頂点を基にして、当該樹冠の樹冠高から、隣接して存在する樹冠を領域分割する領域分割アルゴリズムを用いて、個々の樹冠が区分された精密樹冠画像データを作成する精密樹冠画像データ作成処理ステップをさらに有することを特徴とする森林資源情報算定方法。
【請求項１０】
請求項９に記載の森林資源情報算定方法において、
前記領域分割アルゴリズムは、Watershedアルゴリズムであることを特徴とする森林資源情報算定方法。
【請求項１１】
請求項９又は１０に記載の森林資源情報算定方法において、
前記精密樹冠画像データ作成処理ステップによって作成された前記精密樹冠画像データから、前記各樹木に対応する精密樹冠情報を作成する精密樹冠情報作成処理ステップと、
前記精密樹冠情報作成処理ステップによって作成された前記各樹木に対応する精密樹冠情報に基づいて、前記調査対象森林域に関する森林資源情報を作成する森林資源情報算定処理ステップと、
をさらに有し、
前記各樹木に対応する精密樹冠情報には、少なくとも、前記各樹木に対応する樹冠のラベル番号、前記各樹木に対応する樹冠の樹冠位置、前記各樹木に対応する樹冠の樹冠直径及び前記各樹木に対応する樹冠の樹冠面積が含まれており、
前記森林資源情報算定処理ステップには、
前記各樹木において人間のほぼ胸の高さの直径を各樹木の胸高直径として算定する胸高直径算定処理ステップと、
前記各樹木の材積を算定する材積算定処理ステップと、
前記各樹木の胸高直径、前記各樹木の材積、前記精密樹冠情報作成処理ステップで作成された前記精密樹冠情報及び前記階層別樹頂点抽出ステップで抽出された樹頂点に基づいて得られる各樹木の樹高を集計する処理を行う森林資源情報集計処理ステップと、が含まれていることを特徴とする森林資源情報算定方法。
【請求項１２】
調査対象森林域に関する森林資源情報を作成する森林資源情報算定装置であって、
前記調査対象森林域を含む地域に対応する３次元点群データと地理的情報とに基づいて調査対象森林域画像データを作成する調査対象森林域画像データ作成部と、
前記調査対象森林域画像データに基づいて、前記調査対象森林域に存在する各樹木の樹冠の高さを樹冠高として表す樹冠高画像データの作成を行う樹冠高画像データ作成部と、
前記樹冠高画像データから得られる樹冠高のうち、最大樹冠高が存在する階層を第１階層として、前記最大樹冠高を起点に所定長さの階層幅を有する階層を前記第１階層から地表面に向かって複数設定したときに、前記第１階層を含む各階層に樹頂点が存在する樹冠の樹頂点を、前記階層ごとにそれぞれ抽出する階層別樹頂点抽出部と、を有することを特徴とする森林資源情報算定装置。
【請求項１３】
請求項１２に記載の森林資源情報算定装置において、
前記樹冠高画像データ作成部は、前記調査対象森林域画像データに基づいて、メッシュ化されたデジタル表層モデルデータとデジタル標高モデルデータとを作成し、メッシュごとに前記デジタル表層モデルデータと前記デジタル標高モデルデータとの差分を取ることによって、メッシュごとに樹冠高が求められた樹冠高画像データを作成することを特徴とする森林資源情報算定装置。
【請求項１４】
請求項１３に記載の森林資源情報算定装置において、
前記階層別樹頂点抽出部には、
前記メッシュごとに樹冠高が求められた樹冠高画像データを用いて、前記第１階層を含む各階層に存在する個々の樹冠に対応する階層別樹冠高点群データを、前記階層ごとに作成する階層別樹冠高点群データ作成部と、
前記階層別樹冠高点群データ作成部によって作成された前記個々の樹冠に対応する階層別樹冠高点群データから、前記個々の樹冠に対応する階層別樹冠高点群データの最大値を前記階層ごとに算定して、算定された最大値を当該階層別樹冠高点群データに対応する樹冠の樹頂点とする樹頂点算定部と、
が含まれていることを特徴とする森林資源情報算定装置。
【請求項１５】
請求項１４に記載の森林資源情報算定装置において、
前記樹頂点算定部は、前記最大値を算定する際においては、前記個々の樹冠に対応する前記階層別樹冠高点群データのうち、前記階層別樹冠高点群データが既に作成されている樹冠と同じ樹冠に対応する階層別樹冠高点群データは除外して、前記最大値を算定することを特徴とする森林資源情報算定装置。
【請求項１６】
請求項１４又は１５に記載の森林資源情報算定装置において、
前記階層別樹頂点抽出部には、
前記樹冠高画像データの各メッシュをそれぞれ処理対象メッシュとし、当該処理対象メッシュに平均化フィルターを重ね合わせ、当該平均化フィルターに含まれる各メッシュの樹冠高の平均値を当該処理対象メッシュの樹冠高とする処理を、各処理対象メッシュにおいて順次行うことによって、ノイズ除去済みの樹冠高画像データを作成する平均化処理部がさらに含まれており、
前記階層別樹冠高点群データ作成部は、前記平均化処理部によって作成されたノイズ除去済みの樹冠高画像データを、前記メッシュごとに樹冠高が求められた樹冠高画像データとして用いて前記階層別樹冠高点群データを作成することを特徴とする森林資源情報算定装置。
【請求項１７】
請求項１２〜１６のいずれかに記載の森林資源情報算定装置において、
前記階層幅は、０．２ｍ〜２．０ｍの範囲の長さに設定されていることを特徴とする森林資源情報算定装置。
【請求項１８】
請求項１２〜１７のいずれかに記載の森林資源情報算定装置において、
前記３次元点群データは、上空から撮影箇所を重複させ、かつ、ずらしながら撮影して得られた撮影画像データから、３次元形状復元技術であるＳＦＭ(Structure from Motion)によって作成されていることを特徴とする森林資源情報算定装置。
【請求項１９】
請求項１２〜１７のいずれかに記載の森林資源情報算定装置において、
前記３次元点群データは、上空からレーザー光を照射して得られたレーザー計測データに基づいて作成されていることを特徴とする森林資源情報算定装置。
【請求項２０】
請求項１２〜１９のいずれかに記載の森林資源情報算定装置において、
前記階層ごとにそれぞれ抽出された樹冠の樹頂点を基にして、当該樹冠の樹冠高から、隣接して存在する樹冠を領域分割する領域分割アルゴリズムを用いて、個々の樹冠が区分された精密樹冠画像データを作成する精密樹冠画像データ作成部をさらに有することを特徴とする森林資源情報算定装置。
【請求項２１】
請求項２０に記載の森林資源情報算定装置において、
前記領域分割アルゴリズムは、Watershedアルゴリズムであることを特徴とする森林資源情報算定装置。
【請求項２２】
請求項２０又は２１に記載の森林資源情報算定装置において、
前記精密樹冠画像データ作成部によって作成された前記精密樹冠画像データから、前記各樹木に対応する精密樹冠情報を作成する精密樹冠情報作成部と、
前記精密樹冠情報作成部によって作成された前記各樹木に対応する精密樹冠情報に基づいて、前記調査対象森林域に関する森林資源情報を作成する森林資源情報算定部と、
をさらに有し、
前記各樹木に対応する精密樹冠情報には、少なくとも、前記各樹木に対応する樹冠のラベル番号、前記各樹木に対応する樹冠の樹冠位置、前記各樹木に対応する樹冠の樹冠直径及び前記各樹木に対応する樹冠の樹冠面積が含まれており、
前記森林資源情報算定部には、
前記各樹木において人間のほぼ胸の高さの直径を各樹木の胸高直径として算定する胸高直径算定部と、
前記各樹木の材積を算定する材積算定部と、
前記各樹木の胸高直径、前記各樹木の材積、前記精密樹冠情報作成部で作成された前記精密樹冠情報及び前記階層別樹頂点抽出部で抽出された樹頂点に基いて得られる各樹木の樹高を集計する処理を行う森林資源情報集計部と、が含まれていることを特徴とする森林資源情報算定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、調査対象森林域を含む地域に対応する画像データに基づいて、調査対象森林域に関する森林資源情報を作成する森林資源情報算定方法及び森林資源情報算定装置に関する。
【背景技術】
【０００２】
樹高を含む森林資源に関する情報（森林資源情報という。）は、森林管理において重要な基礎情報である。森林資源情報は、従来においては、地上調査によって、調査員が標準的な森林と思われる箇所に小面積の調査地を設定して、当該調査地における樹木の位置、胸高直径、樹高を測定し、その結果に管理面積を乗ずることで、全体の森林資源情報を推定する方法が広く行われている。
【０００３】
このような森林資源情報の推定方法では、森林資源の推定精度を上げるためには、現地調査箇所の数を増やす必要があるが、そうとすると、多大な労力と費用を要する。地上調査で最も誤差が生じやすい樹高算定は、森林内は障害物や枝葉が茂っていることから、地上からは樹木の頂点が見えない場合がほとんどであり、林縁や孤立木などの標本木を使用して樹高を測定し、実測した胸高直径と標本木の樹高との関係から回帰式を作成し、樹高を推定し、内挿することを行っている。さらに、測定者の技量や測定斜面などの地形の影響を受けると測定誤差が大きくなる。
【０００４】
そこで、最近では、航空レーザー計測で得られた３次元点群データから、樹木頂点や樹木抽出を評価することも行われている。また、ドローンや航空機等で上空から森林を撮影した画像を３次元点群データに変換して、樹木頂点や樹木抽出を評価することも行われている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特許第４９４６０７２号公報
【特許文献２】特許第４２７９８９４号公報
【特許文献３】特許第５５０７４１８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
特許文献１及び特許文献２に記載されている技術は、３次元点群データから樹木頂点や樹木位置を公知の局所最大値フィルターといわれる空間フィルタリング法を用いて、一定サイズ（３×３、５×５等）のフィルターを樹冠高点群データ上で移動させながら領域内の最大値となる樹木頂点を一括処理抽出する方法である。
【０００７】
特許文献１及び特許文献２に記載されている技術は、碁盤のように区画された人工林で、かつ、単純で疎な針葉樹では有効な方法であるが、殆どの森林は、立木密度の違いや樹高及び樹冠の大小が混交しているため、「未抽出や誤抽出が多くなる」、「フィルターのサイズに依存し、小さい樹冠と大きい樹冠を同時に抽出することはできない」、「立木密度の高い森林や集団木では１つの樹高しか抽出できず未抽出木が多い」、「森林の樹冠の隙間（ギャップ）では、下層にある植生を樹木頂点として誤抽出してしまう」など、実用面での課題は多い。
【０００８】
一方、特許文献３に記載されている技術は、航空レーザー計測により取得された森林の３次元点群データを地上に平面投影して樹木位置を検出する方法であるが、立木密度の違いや樹高及び樹冠の大小が混交する森林を単木ごとに区分することは難しいといった課題がある。
【０００９】
このように、上記特許文献１〜３においては、立木密度の違いや樹高及び樹冠の大小が混交する森林においては、各樹木の樹高、胸高直径及びこれら樹高と胸高直径とから得られる材積などの森林資源情報を高精度に算定できないといった問題がある。なお、この明細書において、「樹冠」とは森林を構成する各樹木を上部から見たときに各樹木の枝と葉で構成された部分をいう。
【００１０】
そこで本発明は、上記課題を解決するためになされたものであり、立木密度の違いや樹高及び樹冠の大小が混交する森林においても、森林資源情報を高精度に算定できる森林資源情報算定方法及び森林資源算情報定装置を提供することを目的とする。
【課題を解決するための手段】
【００１１】
［１］本発明の森林資源情報算定方法は、調査対象森林域に関する森林資源情報を作成する森林資源情報算定方法であって、前記調査対象森林域を含む地域に対応する３次元点群データと地理的情報とに基づいて調査対象森林域画像データを作成する調査対象森林域画像データ作成処理ステップと、前記調査対象森林域画像データに基づいて、前記調査対象森林域に存在する各樹木の樹冠の高さを樹冠高として表す樹冠高画像データの作成を行う樹冠高画像データ作成処理ステップと、前記樹冠高画像データから得られる樹冠高のうち、最大樹冠高が存在する階層を第１階層として、前記最大樹冠高を起点に所定長さの階層幅を有する階層を前記第１階層から地表面に向かって複数設定したときに、前記第１階層を含む各階層に樹頂点が存在する樹冠の樹頂点を、前記階層ごとにそれぞれ抽出する階層別樹頂点抽出処理ステップと、を有することを特徴とする。
【００１２】
このように、本発明の森林資源情報算定方法においては、調査対象森林域を含む地域に対応する３次元点群データと地理的情報とに基づいて作成された調査対象森林域画像データに基づいて、調査対象森林域に存在する樹木の樹冠高を表す樹冠高画像データを作成し、当該樹冠高画像データから得られる樹冠高のうち、最大樹冠高が存在する階層を第１階層として、当該第１階層を含む各階層に樹頂点が存在する樹冠の樹頂点を、階層ごとにそれぞれ抽出するようにしている。
【００１３】
これにより、本発明の森林資源情報算定方法によれば、調査対象森林域に存在する調査対象樹木の樹高にかかわらず各樹木に対応する樹冠を高精度に抽出することができ、それによって、立木密度の違いや樹高及び樹冠の大小が混交する森林においても、森林資源情報を高精度に算定できる。
【００１４】
［２］本発明の森林資源情報算定方法においては、前記樹冠高画像データ作成処理ステップは、前記調査対象森林域画像データに基づいて、メッシュ化されたデジタル表層モデルデータとデジタル標高モデルデータとを作成し、メッシュごとに前記デジタル表層モデルデータと前記デジタル標高モデルデータとの差分を取ることによって、メッシュごとに樹冠高が求められた樹冠高画像データを作成することが好ましい。
【００１５】
これにより、樹冠高（地表からの樹冠の高さ）をメッシュごとに高精度に算定することができる。なお、この場合、デジタル表層モデルデータというのは、樹木などの表層の標高を表すデータを指し、デジタル標高モデルデータというのは、地表の標高を表すデータを指している。従って、デジタル表層モデルデータと前記デジタル標高モデルデータとの差分を取ることによって樹冠高（地表からの樹冠の高さ）をメッシュごとに求めることができる。
【００１６】
［３］本発明の森林資源情報算定方法においては、前記階層別樹頂点抽出処理ステップには、前記メッシュごとに樹冠高が求められた樹冠高画像データを用いて、前記第１階層を含む各階層に存在する個々の樹冠に対応する階層別樹冠高点群データを、前記階層ごとに作成する階層別樹冠高点群データ作成処理ステップと、前記階層別樹冠高点群データ作成処理ステップによって作成された前記個々の樹冠に対応する階層別樹冠高点群データから、前記個々の樹冠に対応する階層別樹冠高点群データの最大値を前記階層ごとに算定して、算定された最大値を当該階層別樹冠高点群データに対応する樹冠の樹頂点とする樹頂点算定処理ステップと、が含まれていることが好ましい。
【００１７】
このように、本発明の森林資源情報算定方法においては、第１階層を含む各階層に存在する個々の樹冠に対応する階層別樹冠高点群データを階層ごとに作成し、当該階層ごとに作成された個々の樹冠に対応する階層別樹冠高点群データから、個々の樹冠に対応する階層別樹冠高点群データの最大値を算定するようにしている。ここで、階層別樹冠高点群データを階層ごとに作成し、階層ごとに作成された個々の樹冠に対応する階層別樹冠高点群データから、個々の樹冠に対応する階層別樹冠高点群データの最大値を算定する処理は、例えて言えば、個々の樹冠をＣＴ（Computed Tomography）スキャンのように、階層ごとに輪切りにして階層別樹冠高点群データを作成して、階層ごとに求められた階層別樹冠高点群データから最大値を算定するといった処理である。このような処理を行うことによって、各階層に存在する個々の樹冠の最大値を、階層ごとに算定できる。なお、階層ごとに算定された「樹冠の樹頂点」は、当該樹冠を有する樹木の樹頂点であり、当該算定された樹頂点から当該樹木の樹高を得ることができる。
【００１８】
［４］本発明の森林資源情報算定方法においては、前記樹頂点算定処理ステップは、前記最大値を算定する際においては、前記個々の樹冠に対応する前記階層別樹冠高点群データのうち、前記階層別樹冠高点群データが既に作成されている樹冠と同じ樹冠に対応する階層別樹冠高点群データは除外して、前記最大値を算定することが好ましい。
【００１９】
これにより、各階層のうち現時点の処理対象階層に樹頂点が存在する樹冠のみの最大値を算定できる。例えば、現時点の処理対象階層を第２階層とした場合、当該第２階層よりも上の階層（第１階層）に樹頂点が存在する樹冠の最大値が当該階層の樹頂点として誤算出されることを防止できる。それによって、現時点の処理対象階層（例えば、第２階層）に樹頂点が存在する樹冠のみの最大値を算定できる。
【００２０】
［５］本発明の森林資源情報算定方法においては、前記階層別樹頂点抽出処理ステップには、前記樹冠高画像データの各メッシュをそれぞれ処理対象メッシュとし、当該処理対象メッシュに平均化フィルターを重ね合わせ、当該平均化フィルターに含まれる各メッシュの樹冠高の平均値を当該処理対象メッシュの樹冠高とする処理を、各処理対象メッシュにおいて順次行うことによって、ノイズ除去済みの樹冠高画像データを作成する平均化処理ステップがさらに含まれており、前記階層別樹冠高点群データ作成処理ステップは、前記平均化処理ステップによって作成されたノイズ除去済みの樹冠高画像データを、前記メッシュごとに樹冠高が求められた樹冠高画像データとして用いて前記階層別樹冠高点群データを作成することが好ましい。
【００２１】
このような平均化処理を行うことにより、例えば、樹冠に極端な凹凸（枝葉が周囲に比べて疎になっている部分や徒長枝葉など）が存在している場合に、これをノイズとして除去することができ、ノイズ除去済みの樹冠高画像データを作成することができる。そして、階層別樹冠高点群データ作成処理ステップは、ノイズ除去済みの樹冠高画像データを樹冠高画像データとして用いて階層別樹冠高点群データを作成することによって、高精度な階層別樹冠高点群データを作成することができる。なお、徒長枝葉というのは、樹冠において枝葉が周囲の枝葉に比べて突出している部分を指す。
【００２２】
［６］本発明の森林資源情報算定方法においては、前記階層幅は、０．２ｍ〜２．０ｍの範囲の長さに設定されていることが好ましい。
【００２３】
階層幅は、特に限定されるものではないが、あまり細かく（長さを短く）設定すると、階層数が増大して処理回数が増えてデータ量が多くなり、演算に時間がかかるといった課題が生じ、逆に、あまり大きく（長さを長く）設定すると、階層別樹冠高点群データの精度が低くなるといった課題が生じる。このため、階層幅は、０．２ｍ〜２．０ｍの範囲に設定することが好まく、例えば、０．４ｍ、０．８ｍ、１．０ｍ程度の階層幅を例示できるが、調査対象森林域内に存在する樹木の樹高なども考慮して設定することが好ましい。
【００２４】
［７］本発明の森林資源情報算定方法においては、前記３次元点群データは、上空から撮影箇所を重複させ、かつ、ずらしながら撮影して得られた撮影画像データから、３次元形状復元技術であるＳＦＭ(Structure from Motion)によって作成されていることが好ましい。
【００２５】
３次元形状復元技術であるＳＦＭを用いることによって、２次元の撮影画像データから、歪みのない地図補正されたオルソ画像を作成することができる。そして、このようにして作成されたオルソ画像に基づいて、高さが表わされている３次元点群データを作成することができる。
【００２６】
［８］本発明の森林資源情報算定方法においては、前記３次元点群データは、上空からレーザー光を照射して得られたレーザー計測データに基づいて作成されていてもよい。
【００２７】
このように、上空からレーザー光を照射して得られたレーザー計測データによっても、３次元点群データを作成することができる。
【００２８】
［９］本発明の森林資源情報算定方法においては、前記階層ごとにそれぞれ抽出された樹冠の樹頂点を基にして、当該樹冠の樹冠高から、隣接して存在する樹冠を領域分割する領域分割アルゴリズムを用いて、個々の樹冠が区分された精密樹冠画像データを作成する精密樹冠画像データ作成処理ステップをさらに有することが好ましい。
【００２９】
このような領域分割アルゴリズムを用いることによって、隣接する樹冠が個々の樹冠に区分された精密な樹冠画像データ（精密樹冠画像データ）を作成することができる。このようにして作成された精密樹冠画像データは、調査対象樹木の各樹木に対応した個々の樹冠を高精度に抽出可能な精密樹冠画像データとなる。このため、当該精密樹冠画像データに基づいて作成した精密樹冠情報は、調査対象樹木の各樹木に対応した高精度な樹冠情報となり、当該精密樹冠情報に基づいて、森林資源に関する情報（森林資源情報）を算定することにより、高精度な森林資源情報を得ることができる。
【００３０】
［１０］本発明の森林資源情報算定方法においては、前記領域分割アルゴリズムは、Watershedアルゴリズムであることが好ましい。
【００３１】
Watershedアルゴリズムは、隣接して密集する多数の物体同士が繋がった状態で存在する場合において、個々の物体の境界を高精度に抽出して個々の物体を区分させることができるアルゴリズムである。従って、Watershedアルゴリズムを用いることによって、密集して存在する樹冠を、樹冠の大きさ（樹冠直径や高さ）に係わらず個々の樹冠として高精度に区分することができる。
【００３２】
［１１］本発明の森林資源情報算定方法においては、前記精密樹冠画像データ作成処理ステップによって作成された前記精密樹冠画像データから、前記各樹木に対応する精密樹冠情報を作成する精密樹冠情報作成処理ステップと、前記精密樹冠情報作成処理ステップによって作成された前記各樹木に対応する精密樹冠情報に基づいて、前記調査対象森林域に関する森林資源情報を作成する森林資源情報算定処理ステップと、をさらに有し、前記各樹木に対応する精密樹冠情報には、少なくとも、前記各樹木に対応する樹冠のラベル番号、前記各樹木に対応する樹冠の樹冠位置、前記各樹木に対応する樹冠の樹冠直径及び前記各樹木に対応する樹冠の樹冠面積が含まれており、前記森林資源情報算定処理ステップには、前記各樹木において人間のほぼ胸の高さの直径を各樹木の胸高直径として算定する胸高直径算定処理ステップと、前記各樹木の材積を算定する材積算定処理ステップと、前記各樹木の胸高直径、前記各樹木の材積、前記精密樹冠情報作成処理ステップで作成された前記精密樹冠情報及び前記階層別樹頂点抽出ステップで抽出された樹頂点に基づいて得られる各樹木の樹高を集計する処理を行う森林資源情報集計処理ステップと、が含まれていることが好ましい。
【００３３】
精密樹冠画像データに基づいて作成した精密樹冠情報は、調査対象樹木の各樹木に対応した高精度な樹冠情報となる。このため、精密樹冠画像データに基づいて作成した精密樹冠情報に基づいて、森林資源情報算定処理ステップが、森林資源に関する情報（森林資源情報）を算定することにより、高精度な森林資源情報を得ることができる。これにより、計測誤差が少なく客観性の高い高精度な森林資源情報を、調査対象森林域の全域、小班ごと、任意の範囲において提供できることから、間伐などの森林施業立案などの森林管理に活用できる。具体的には、例えば、森林資源情報算定処理ステップが算定した算定結果を森林資源の概要を示す森林資源概要情報として属性データベースに登録する処理を行い、当該属性データベースに登録されている森林資源概要情報を、一覧表としてディスプレイ上に表示することによって、調査対処森林域の森林資源の概要を把握することができる。
【００３４】
［１２］本発明の森林資源情報算定装置は、調査対象森林域に関する森林資源情報を作成する森林資源情報算定装置であって、前記調査対象森林域を含む地域に対応する３次元点群データと地理的情報とに基づいて、前記調査対象森林域画像データを作成する調査対象森林域画像データ作成部と、前記調査対象森林域画像データに基づいて、前記調査対象森林域に存在する各樹木の樹冠の高さを樹冠高として表す樹冠高画像データの作成を行う樹冠高画像データ作成部と、前記樹冠高画像データから得られる樹冠高のうち、最大樹冠高が存在する階層を第１階層として、前記最大樹冠高を起点に所定長さの階層幅を有する階層を前記第１階層から地表面に向かって複数設定したときに、前記第１階層を含む各階層に樹頂点が存在する樹冠の樹頂点を、前記階層ごとにそれぞれ抽出する階層別樹頂点抽出部と、を有することを特徴とする。
【００３５】
［１３］本発明の森林資源情報算定装置においては、前記樹冠高画像データ作成部は、前記調査対象森林域画像データに基づいて、メッシュ化されたデジタル表層モデルデータとデジタル標高モデルデータとを作成し、メッシュごとに前記デジタル表層モデルデータと前記デジタル標高モデルデータとの差分を取ることによって、メッシュごとに樹冠高が求められた樹冠高画像データを作成することが好ましい。
【００３６】
［１４］本発明の森林資源情報算定装置においては、前記階層別樹頂点抽出部には、
前記メッシュごとに樹冠高が求められた樹冠高画像データを用いて、前記第１階層を含む各階層に存在する個々の樹冠に対応する階層別樹冠高点群データを、前記階層ごとに作成する階層別樹冠高点群データ作成部と、前記階層別樹冠高点群データ作成部によって作成された前記個々の樹冠に対応する階層別樹冠高点群データから、前記個々の樹冠に対応する階層別樹冠高点群データの最大値を前記階層ごとに算定して、算定された最大値を当該階層別樹冠高点群データに対応する樹冠の樹頂点とする樹頂点算定部と、が含まれていることが好ましい。
【００３７】
［１５］本発明の森林資源情報算定装置においては、前記樹頂点算定部は、前記最大値を算定する際においては、前記個々の樹冠に対応する前記階層別樹冠高点群データのうち、前記階層別樹冠高点群データが既に作成されている樹冠と同じ樹冠に対応する階層別樹冠高点群データは除外して、前記最大値を算定することが好ましい。
【００３８】
［１６］本発明の森林資源情報算定装置においては、前記階層別樹頂点抽出部には、
前記樹冠高画像データの各メッシュをそれぞれ処理対象メッシュとし、当該処理対象メッシュに平均化フィルターを重ね合わせ、当該平均化フィルターに含まれる各メッシュの樹冠高の平均値を当該処理対象メッシュの樹冠高とする処理を、各処理対象メッシュにおいて順次行うことによって、ノイズ除去済みの樹冠高画像データを作成する平均化処理部がさらに含まれており、前記階層別樹冠高点群データ作成部は、前記平均化処理部によって作成されたノイズ除去済みの樹冠高画像データを、前記メッシュごとに樹冠高が求められた樹冠高画像データとして用いて前記階層別樹冠高点群データを作成することが好ましい。
【００３９】
［１７］本発明の森林資源情報算定装置においては、前記階層幅は、０．２ｍ〜２．０ｍの範囲の長さに設定されていることが好ましい。
【００４０】
［１８］本発明の森林資源情報算定装置においては、前記３次元点群データは、上空から撮影箇所を重複させ、かつ、ずらしながら撮影して得られた撮影画像データから、３次元形状復元技術であるＳＦＭ(Structure from Motion)によって作成されていることが好ましい。
【００４１】
［１９］本発明の森林資源情報算定装置においては、前記３次元点群データは、上空からレーザー光を照射して得られたレーザー計測データに基づいて作成されていてもよい。
【００４２】
［２０］本発明の森林資源情報算定装置においては、前記階層ごとにそれぞれ抽出された樹冠の樹頂点を基にして、当該樹冠の樹冠高から、隣接して存在する樹冠を領域分割する領域分割アルゴリズムを用いて、個々の樹冠が区分された精密樹冠画像データを作成する精密樹冠画像データ作成部をさらに有することが好ましい。
【００４３】
［２１］本発明の森林資源情報算定装置においては、前記領域分割アルゴリズムは、Watershedアルゴリズムであることが好ましい。
【００４４】
［２２］本発明の森林資源情報算定装置においては、前記精密樹冠画像データ作成部によって作成された前記精密樹冠画像データから、前記各樹木に対応する精密樹冠情報を作成する精密樹冠情報作成部と、前記精密樹冠情報作成部によって作成された前記各樹木に対応する精密樹冠情報に基づいて、前記調査対象森林域に関する森林資源情報を作成する森林資源情報算定部と、をさらに有し、前記各樹木に対応する精密樹冠情報には、少なくとも、前記各樹木に対応する樹冠のラベル番号、前記各樹木に対応する樹冠の樹冠位置、前記各樹木に対応する樹冠の樹冠直径及び前記各樹木に対応する樹冠の樹冠面積が含まれており、前記森林資源情報算定部には、前記各樹木において人間のほぼ胸の高さの直径を各樹木の胸高直径として算定する胸高直径算定部と、前記各樹木の材積を算定する材積算定部と、前記各樹木の胸高直径、前記各樹木の材積、前記精密樹冠情報作成部で作成された前記精密樹冠情報及び前記階層別樹頂点抽出部で抽出された樹頂点に基いて得られる各樹木の樹高を集計する処理を行う森林資源情報集計部と、が含まれていることが好ましい。
【００４５】
なお、［１２］〜［２２］に記載した本発明の森林資源情報算定装置においても、[１]〜［１１］に記載したそれぞれ対応する本発明の森林資源情報算定方法で得られる効果と同様の効果が得られる。
【図面の簡単な説明】
【００４６】
【図１】実施形態に係る森林資源情報算定方法を説明するために示すフローチャートである。
【図２】ディスプレイ上に表示された調査対象森林域画像を示す図である。
【図３】ディスプレイ上に表示された樹冠高画像である。
【図４】図１に示されている階層別樹頂点抽出処理（ステップＳ５０）を詳細に説明するフローチャートである。
【図５】ディスプレイ上に表示された第１階層樹冠高点群領域Ａ１，Ａ２，Ａ３，・・・を示す図である。
【図６】図５に示されている第１階層樹冠高点群領域Ａ１，Ａ２，Ａ３，・・・において抽出された樹頂点ａ１，ａ２，ａ３，・・・を示す図である。
【図７】ディスプレイ上に表示された第２階層樹冠高点群領域Ｂ１，Ｂ２，Ｂ３，・・・を示す図である。
【図８】ディスプレイ上に表示された第３階層樹冠高点群領域Ｃ１，Ｃ２，Ｃ３，・・・を示す図である。
【図９】第１階層から第５層までの各階層の階層別樹冠高点群領域と当該階層別樹冠高点群領域における樹頂点とを示す図である。
【図１０】階層別樹冠高点群データ作成処理（ステップＳ５３）及び樹頂点算定処理（ステップＳ５４）を第１階層（最上層）から最下層まで順次行うことによって得られた全階層の階層別樹冠高点群領域と樹頂点を示す図である。
【図１１】ディスプレイ上に表示された精密樹冠画像を示す図である。
【図１２】ディスプレイ上に表示された単木資源情報の一覧表を示す図である。
【図１３】ディスプレイ上に表示された森林資源概要情報の一覧表を示す図である。
【図１４】実施形態に係る森林資源情報算定装置１を説明するために示す図である。
【図１５】図１４に示されている階層別樹頂点抽出部５０の構成を示す図である。
【発明を実施するための形態】
【００４７】
図１は、実施形態に係る森林資源情報算定方法を説明するために示すフローチャートである。実施形態に係る森林資源情報算定方法は、調査対象森林域を含む地域に対応する画像データに基づいて、調査対象森林域に関する森林資源情報を作成する森林資源情報算定方法である。ここで、調査対象森林域を含む地域に対応する画像データは、上空から調査対象森林域を含む地域を空撮することによって取得する。なお、空撮するには、航空機、ドローンなどを用いることができるが、実施形態に係る森林資源情報算定方法においては、普及型のドローンを用いることとする。
【００４８】
また、実施形態に係る森林資源情報算定方法において、調査対象樹木というのは、例えば、カラマツ、アカマツ、スギ、ヒノキなどの高木となる植栽木を指しており、調査対象森林域というのは、これらの樹木によって形成される林が所定の広さを有して存在する森林域であるとする。また、「調査対象樹木」を一本ごとの樹木として説明する場合には、単に「各樹木」と略記する場合もある。このことは、後述する実施形態に係る森林資源情報算定装置においても同様である。
【００４９】
実施形態に係る森林資源情報算定方法において行われる処理は、図１に示すように、ドローンによって撮影した空撮画像データをパーソナルコンピュータなどの情報処理装置に入力するデータ入力処理（ステップＳ１０）と、入力された空撮画像データに基づいて３次元点群データを作成する３次元点群データ作成処理（ステップＳ２０）と、３次元点群データ作成処理（ステップＳ２０）によって作成された３次元点群データと地理情報システム（ＧＩＳ）１００から得られる地理的情報地理的情報とに基づいて、調査対象森林域画像データを作成する調査対象森林域画像データ作成処理（ステップＳ３０）と、調査対象森林域画像データ作成処理（ステップＳ３０）によって作成された調査対象森林域画像データに基づいて、調査対象森林域に存在する各樹木の樹冠の高さを樹冠高として表す樹冠高画像データの作成を行う樹冠高画像データ作成処理（ステップＳ４０）と、樹冠高画像データ作成処理（ステップＳ４０）によって作成された樹冠高画像データから得られる樹冠高のうち、最大樹冠高が存在する階層を第１階層（最上層）として、所定長さの階層幅を有する階層を第１階層から地表面に向かって複数設定したときに、第１階層を含む各階層に樹頂点が存在する樹冠の樹頂点を、階層ごとにそれぞれ抽出する階層別樹頂点抽出処理（ステップＳ５０）と、を有している。
【００５０】
また、これら各ステップＳ１０〜Ｓ５０の処理に加えて、個々の樹冠が区分された精密樹冠画像データを作成する精密樹冠画像データ作成処理（ステップＳ６０）と、精密樹冠画像データ作成処理（ステップＳ６０）によって作成された精密樹冠画像データから、各樹木に対応する精密樹冠情報を作成する精密樹冠情報作成処理（ステップＳ７０）と、当該精密樹冠情報作成処理（ステップＳ７０）によって作成された各樹木に対応する精密樹冠情報に基づいて、調査対象森林域に関する森林資源情報を作成する森林資源情報算定処理（ステップＳ８０）と、を有している。
なお、精密樹冠画像データ作成処理は、階層別樹頂点抽出処理（ステップＳ５０）によって階層ごとにそれぞれ抽出された樹冠の樹頂点を基にして、当該樹冠の樹冠高から、隣接して存在する樹冠を領域分割する領域分割アルゴリズムを用いて、個々の樹冠が区分された精密樹冠画像データを作成する。
【００５１】
また、森林資源情報算定処理（ステップＳ８０）には、各樹木において人間のほぼ胸の高さの直径を各樹木の胸高直径として算定する胸高直径算定処理（ステップＳ８１）と、各樹木の材積を算定する材積算定処理（ステップＳ８２）と、各樹木の胸高直径、各樹木の材積、精密樹冠情報作成処理（ステップＳ７０）で作成された精密樹冠情報、階層別樹頂点抽出処理（ステップＳ５０）で抽出された樹頂点に基いて得られる各樹木の樹高を集計する処理を行う森林資源情報集計処理（ステップＳ８３）と、が含まれている。
【００５２】
また、階層別樹頂点抽出処理（ステップＳ５０）には、複数の処理が含まれており、これらの各処理については、図４により詳細に説明する。
【００５３】
続いて、上述した各ステップＳ１０〜Ｓ８０のうちの３次元点群データ作成処理（ステップＳ２０）以降の各処理について詳細に説明する。３次元点群データ作成処理（ステップＳ２０）は、調査対象地域を含む所定地域をドローンに搭載されているカメラで、撮影箇所を例えば８０％程度を重複させ、かつ、ずらしながら撮影し、その画像データ（空撮画像データ）から、公知の３次元形状復元技術であるＳＦＭ(Structure from Motion)によって、歪みのない地図補正されたオルソ画像を作成して、作成したオルソ画像から３次元点群データを作成する。
【００５４】
なお、調査対象森林域を含む所定地域全体の広さとしては、２〜３ヘクタール程度である。このような広さを有する調査対象森林域を含む所定地域全体をドローンに搭載されているカメラで、８０％程度を重複させ、かつ、ずらしながら数百枚撮影して、その撮影画像から、ＳＦＭによって、歪みのない地図補正されたオルソ画像を作成する。このようにして作成されたオルソ画像に基づいて、高さが表わされている３次元点群データを作成することができる。なお、ＳＦＭについては、下記の公知文献１に記載されている。
【００５５】
公知文献1：Linda G. Shapiro; George C. Stockman (2001). Computer Vision. Prentice Hall.
【００５６】
調査対象森林域画像データ作成処理（ステップＳ３０）は、地理情報システム（ＧＩＳ）１００から得られる地理的情報（この場合、都道府県などで整備されている森林の管理台帳などに登録されている森林域の境界データ）を用いて、３次元点群データから調査対象となる森林域（調査対象森林域）の画像データを抽出して調査対象森林域画像データとする。そして、調査対象森林域画像データに対応する調査対象森林域画像をディスプレイ上に表示する。
【００５７】
図２は、ディスプレイ上に表示された調査対象森林域画像を示す図である。図２は調査対象森林域のうちの一部が抽出されて示されている。図２に示す調査対象森林域画像によれば、各樹木の樹冠を識別することができる。なお、図２はモノクロ画像であるため、各樹木の樹冠を識別しにくいが、図２の元となるカラー画像上では、各樹木の樹冠を容易に識別できる。
【００５８】
その後、調査対象森林域画像データに基づいて、調査対象森林域に存在する樹木の樹冠高を表す樹冠高画像データを作成する樹冠高画像データ作成処理（ステップＳ４０）を行う。当該樹冠高画像データ作成処理（ステップＳ４０）は、調査対象森林域画像データ作成処理（ステップＳ２０）によって作成された調査対象森林域画像データに基づいて、調査対象森林域に存在する樹木の樹冠高（地表からの樹冠の高さ）を表す樹冠高画像データを作成する。
【００５９】
具体的には、調査対象森林域画像データに基づいて、メッシュ化されたデジタル表層モデルデータ（樹木などの表層の標高を表すデータ）とデジタル標高モデルデータ（地表の標高を表すデータ）を作成し、各メッシュにおいてデジタル表層モデルデータとデジタル標高モデルデータとの差分を取ることによって、各メッシュにおいて樹冠高が表わされた樹冠高画像データを作成する。このようにして作成された当該樹冠高画像データに対応する樹冠高画像をディスプレイ上に表示する。
【００６０】
なお、メッシュサイズは、調査対象森林域画像の地上分解能に基づいて設定することができるが、メッシュサイズを大きくし過ぎると、算定される森林資源情報の精度に課題が生じる場合もあり、また、逆にメッシュサイズを小さくし過ぎると、処理すべきデータ量が多くなる。このため、メッシュサイズは、調査対象森林域の面積、当該調査対象森林域における森林の状況及び算定される森林資源情報の精度などを考慮して最適なサイズに設定することが好ましい。実施形態においては、メッシュサイズは、１０ｃｍ×１０ｃｍとする。
【００６１】
図３は、ディスプレイ上に表示された樹冠高画像である。なお、図３に示す樹冠高画像における黒から白のグラデーションは、地表を基準（０．０メートル）としたときの当該基準に対する高さの変化（例えば、０．０メートルを黒とする。）を表している。なお、図３に示す樹冠高画像に対応する樹冠高画像データは、この段階では、各樹木に対応する樹冠が明確なものとはなってはいない。
【００６２】
続いて、樹冠高画像データ作成処理（ステップＳ４０）によって作成された樹冠高画像データに基づいて、階層別樹頂点抽出処理（ステップＳ５０）を行う。
【００６３】
図４は、図１に示されている階層別樹頂点抽出処理（ステップＳ５０）を詳細に説明するフローチャートである。図４に示すように、階層別樹頂点抽出処理（ステップＳ５０）には、平均化処理（ステップＳ５１）、初期値設定処理（ステップＳ５２）、階層別樹冠高点群データ作成処理（ステップＳ５３）、樹頂点算定処理（ステップＳ５４）が含まれている。また、これらの各処理の他に、Ｈｍａｘ−Ｄ×Ｎｉ＝Ｈｉ（ステップＳ５５）、Ｈｉ＞０（ステップＳ５６）、Ｎｉ＝Ｎｉ＋１（ステップＳ５７）といった演算処理が含まれている。なお、Ｈｍａｘは調査対象森林域内の最大樹冠高、Ｄは階層幅、Ｎｉは現時点の処理対象階層、Ｈｉは次の処理対象階層内の最大樹冠項であり、これらについては後述する。
【００６４】
平均化処理（ステップＳ５１）は、図３に示す樹冠高画像において、例えば、樹冠に極端な凹凸（枝葉が周囲に比べて疎になっている部分や徒長枝葉など）が存在している場合に、これらをノイズとして除去して、ノイズが除去された樹冠高画像データを作成する処理である。当該、平均化処理(ステップＳ５１)は、樹冠高画像データ作成処理（ステップＳ４０）によって作成された樹冠高画像データの各メッシュをそれぞれ処理対象メッシュとし、当該処理対象メッシュに平均化フィルターを重ね合わせ、当該平均化フィルターに含まれる各メッシュの樹冠高の平均値を当該処理対象メッシュの樹冠高とする。このような処理を、各処理対象メッシュにおいて順次行うことによって、ノイズ除去済みの樹冠高画像データを作成する。
【００６５】
平均化処理を簡単化して具体的に説明する。１つのメッシュのサイズを一辺が仮に１０ｃｍとした場合においては、平均化処理の対象となる１つのメッシュ（処理対象メッシュという。）が中心となるように平均化フィルター（フィルターサイズを例えば５メッシュ×５メッシュとする。）を重ね合わせ、当該平均化フィルターに含まれる各メッシュの樹冠高の平均値を当該処理対象メッシュの樹冠高とする処理を、各処理対象メッシュにおいて順次行う。
【００６６】
例えば、平均化フィルターに重ね合わせられている範囲（平均化処理範囲という。）に含まれる各メッシュ（例えば５メッシュ×５メッシュの合計２５メッシュとする。）の各樹冠高を積算して、その合計値を２５で割った値を当該平均化処理範囲の平均値として、その平均値を当該処理対象メッシュの樹冠高とする。このような処理を樹冠高画像データに含まれるすべてのメッシュにおいて行う。なお、平均化フィルターのサイズは、調査対象森林域に存在する調査対象樹木の平均的な樹冠直径などに基づいて設定することができる。
【００６７】
なお、このような平均化処理を行うと、平均化処理によって得られたメッシュごとの樹冠高が、本来の樹冠高（メッシュごとの樹冠高）と異なってしまう場合もあるが、メッシュのサイズが一辺１０ｃｍというように小さいサイズであれば、平均化処理によって得られたメッシュごとの樹冠高は、本来の樹冠高に対して大差のない樹冠高を得ることができる。また、このような平均化処理を行うことによって、個々の樹冠に存在する極端な凹凸（枝葉が周囲に比べて疎になっている部分や徒長枝葉など）が除去されて、個々の樹冠が平滑化されるため、後述するWatershedアルゴリズによる領域分割を行う際においては、高精度な領域分割が行える。
【００６８】
初期値設定処理（ステップＳ５２）は、階層別樹冠高点群データ作成処理（ステップＳ５３）を行う際の初期値を設定する処理であり、ここでの初期値としては、平均化処理（ステップＳ５１）によって得られたメッシュごとの樹冠高のうちの最大樹冠高（調査対象森林域内の最大樹冠高）Ｈｍａｘと、当該最大樹冠高Ｈｍａｘを起点に下方向（地面）に向かって階層を設定したときの階層幅Ｄと、現時点の処理対象階層Ｎｉとを例示できる。
【００６９】
なお、各メッシュの樹冠高のうちの最大樹冠高Ｈｍａｘというのは、各メッシュにおいて得られた樹冠高のうち、あるメッシュにおいて得られた樹冠高が例えば２５ｍであって、当該２５ｍの樹冠高がすべてのメッシュにおいて得られた樹冠高のうちの最大樹冠高であったとすると、Ｈｍａｘ＝２５ｍとする。また、階層幅Ｄは、最大樹冠高Ｈｍａｘを起点に下方向（地面）に向かって、例えば、１ｍごとに階層幅を設定する場合には、Ｄ＝１とする。また、処理対象階層Ｎｉは、初期において処理対象となる階層は、第１階層（最上層）であるため、Ｎｉ＝１とする。
【００７０】
このように、初期値としては、Ｈｍａｘ＝２５、Ｄ＝１、Ｎｉ＝１を設定して、階層別樹冠高点群データ作成処理（ステップＳ５３）を行う。この場合、Ｎｉ＝１であるため、第１階層（最上層）における階層別樹冠高点群データ作成処理を行う。なお、Ｈｍａｘ＝２５ｍであり、また、Ｄ＝１であるため、当該第１階層は、２５ｍ〜２４ｍである。
【００７１】
すなわち、階層別樹冠高点群データ作成処理（ステップＳ５３）は、平均化処理によって算出された各メッシュの樹冠高のうちの最大樹冠高（Ｈｍａｘ＝２５ｍ）が存在する層を第１階層（最上層）として、当該最大樹冠高（Ｈｍａｘ）を起点に所定長さの階層幅（Ｄ＝１ｍ）を設定し、最大樹冠高Ｈｍａｘ＝２５ｍから階層幅ごとに、第１階層を含む各階層に存在する個々の樹冠に対応する階層別樹冠高点群データを、前記階層ごとに順次作成する。
【００７２】
階層別樹冠高点群データ作成処理（ステップＳ５３）について具体的に説明する。この場合、各メッシュの樹冠高のうちの最大樹冠高Ｈｍａｘは２５ｍであり、また、階層幅が１ｍであるため、初期値としては、Ｈｍａｘ＝２５、Ｎｉ＝１（第１階層）を設定する。この場合、まず、各メッシュの樹冠高のうち当該第１階層（２４ｍ〜２５ｍ）に属する樹冠に対応する樹冠高点群データ（第１階層樹冠高点群データとする。）を作成する。なお、第１階層は、２４ｍ〜２５ｍであるが、厳密には、２４．０１ｍ〜２５ｍというように、２４ｍをわずかに超える値から２５ｍまでとする。これは、第１階層よりも下位の階層においても同様であり、例えば、第２階層を２３ｍ〜２４ｍとしたときは、２３．０１ｍ〜２４ｍというように、２３ｍをわずかに超える値から２４ｍまでとする。
そして、階層別樹冠高点群データ作成処理（ステップＳ５３）によって作成された第１階層樹冠高点群データに対応する第１階層樹冠高点群領域をディスプレイ上に表示する。
【００７３】
なお、ここでは、第１階層は２４．０１ｍ〜２５ｍ、第２階層は２３．０１ｍ〜２４ｍというように設定したが、このような設定の仕方に限らず、例えば、第１階層は２４ｍ〜２４．９９ｍ、第２階層は２３ｍ〜２３．９９ｍというように設定してもよい。このことは、他の階層においても同様である。但し、第１階層については、当該第１階層は最上層であり、それより上の階層が存在しないため、２４ｍ〜２５ｍとしてもよい。
【００７４】
図５は、ディスプレイ上に表示された第１階層樹冠高点群領域を示す図である。図５において、黒色で塗りつぶされている領域が第１階層樹冠高点群領域Ａ１，Ａ２，Ａ３，・・・であり、当該第１階層樹冠高点群領域Ａ１，Ａ２，Ａ３，・・・は、樹頂点付近が第１階層（２４ｍ〜２５ｍ）に達している樹木の存在を示している。例えば、図５において、黒枠内に含まれる樹木群の様子（横断面図）を見ると、これらの樹木群の中には、樹頂点付近が、帯状に示されている第１階層にまで達している樹木が存在していることがわかる。
【００７５】
なお、図５において、第１階層樹冠高点群領域Ａ１，Ａ２，Ａ３，・・・は、実際には、多数存在するが、代表的な第１階層樹冠高点群領域のみに符号が付されている。これは、後述する他の階層において抽出される樹冠高点群領域についても同様である。
【００７６】
また、図５はモノクロ画像であるため、図５においては、帯状に示されている第１階層及び第１階層樹冠高点群領域Ａ１，Ａ２，Ａ３，・・・は、それぞれ黒色で示されているが、後述する第２階層以下において得られる各階層との区別を明確にするために、図５の元となるカラー画像上では、帯状に示されている第１階層及び樹冠高点群領域Ａ１，Ａ２，Ａ３，・・・は、それぞれ赤色で示されている。
【００７７】
このようにして、第１階層樹冠高点群領域Ａ１，Ａ２，Ａ３，・・・が抽出されると、当該第１階層樹冠高点群領域Ａ１，Ａ２，Ａ３，・・・に対して、局所最大値フィルタ（例えば、３メッシュ×３メッシュ）を用いて、個々の第１階層樹冠高点群領域Ａ１，Ａ２，Ａ３，・・・のそれぞれにおける最大値を抽出して、抽出した最大値を個々の樹冠の樹頂点とする樹頂点算定処理（ステップＳ５４）を行う。なお、樹冠の樹頂点は、当該樹冠を有する樹木の樹頂点であるため、樹冠の樹頂点を得ることによって当該樹冠を有する樹木の樹高を得ることができる。
【００７８】
図６は、図５に示されている第１階層樹冠高点群領域Ａ１，Ａ２，Ａ３，・・・において抽出された樹頂点ａ１，ａ２，ａ３，・・・を示す図である。なお、現時点の処理対象階層（この場合、第１階層）において得られた当樹頂点は、黒丸の点で示されている。
【００７９】
このようにして、第１階層樹冠高点群領域Ａ１，Ａ２，Ａ３，・・・において樹頂点ａ１，ａ２，ａ３，・・・を抽出する処理が終了したら、Ｈｍａｘ（調査対象森林域内の最大樹冠高）−Ｄ（階層幅）×Ｎｉ（現時点の処理対象階層）＝Ｈｉ（次の処理対象階層内の最大樹冠高）を求める処理（ステップＳ５５）を行う。
【００８０】
なお、次の処理対象階層内の最大樹冠高Ｈｉは、調査対象森林域内の最大樹冠高（Ｈｍａｘ）＝２５ｍ、階層幅（Ｄ）＝１ｍであり、また、この時点での処理対象階層Ｎｉ＝１であるため、２５ｍ−１ｍ×１＝２４ｍとなり、次の処理対象階層内の最大樹冠高Ｈｉは２４ｍとなる。そして、Ｈｉ＞０か否かを判定し（ステップＳ）、この場合、Ｈｉ＝２４であり、Ｈｉ＞０であるため、Ｎｉ＝Ｎｉ＋１を行って（ステップＳ５７）、Ｎｉ＝２として、処理対象階層を第２階層に設定して、階層別樹冠高点群データ作成処理（ステップＳ５３）に戻る。
【００８１】
ここで、第２階層は、樹冠高が２３ｍ〜２４ｍの階層であり、当該第２階層において、階層別樹冠高点群データ作成処理（ステップＳ５３）を行う。すなわち、各メッシュの樹冠高のうち、当該第２階層に属する樹冠高点群データ（第２階層樹冠高点群データとする。）を作成する。そして、この第２階層樹冠高点群データに対応する第２階層樹冠高点群領域をディスプレイ上に表示する。
【００８２】
図７は、ディスプレイ上に表示された第２階層樹冠高点群領域Ｂ１，Ｂ２，Ｂ３，・・・を示す図である。なお、図７においては、第２階層樹冠高点群領域Ｂ１，Ｂ２，Ｂ３，・・・だけではなく、当該第２階層樹冠高点群領域Ｂ１，Ｂ２，Ｂ３，・・・において抽出された樹頂点も示されている。また、図７において、灰色で塗りつぶされている領域が第２階層樹冠高点群領域Ｂ１，Ｂ２，Ｂ３，・・・であり、当該第２階層樹冠高点群領域Ｂ１，Ｂ２，Ｂ３，・・・は、樹頂点付近が第２階層（２３ｍ〜２４ｍ）に達している樹木の存在を示している。例えば、図７において、黒枠内に含まれる樹木群の様子（横断面図）を見ると、これらの樹木群の中には、樹頂点付近が、帯状に示されている第２階層にまで達している樹木が存在していることがわかる。
【００８３】
なお、図７はモノクロ画像であるため、図７においては、帯状に示されている第２階層及び第２階層樹冠高点群領域Ｂ１，Ｂ２，Ｂ３，・・・は、それぞれ灰色で示されているが、他の階層において得られる各階層の樹冠高点群領域との区別を明確にするために、図７の元となるカラー画像上では、帯状に示されている第２階層及び第２階層樹冠高点群領域Ｂ１，Ｂ２，Ｂ３，・・・は、それぞれ橙色で示されている。
【００８４】
なお、図７においては、新たに抽出された第２階層樹冠高点群領域と、既に抽出されている第１階層樹冠高点群領域Ａ１，Ａ２，Ａ３，・・・との区別をし易くするために、既に抽出されている第１階層樹冠高点群領域Ａ１，Ａ２，Ａ３，・・・は白抜きで示されている。
【００８５】
図７に示すように、第２階層樹冠高点群領域Ｂ１，Ｂ２，Ｂ３，・・・のうち、既に抽出されている樹冠高点群領域（この場合、第１階層樹冠高点群領域Ａ１，Ａ２，Ａ３、・・・）と同じ樹冠において抽出された第２階層樹冠高点群領域（例えば、第２階層樹冠高点群領域Ｂ４，Ｂ５など）は、白抜きで示されている領域（第１階層樹冠高点群領域Ａ１，Ａ２）を取り囲むような環状をなす樹冠高点群領域として表示されている。ここで、既に抽出されている樹冠高点群領域と同じ樹冠において抽出された樹冠高点群領域を「既抽出樹冠高点群領域」とし、当該「既抽出冠高点群領域」以外の新たに抽出された樹冠高点群領域を「新規抽出樹冠高点群領域」とする。
【００８６】
このようにして、第２階層樹冠高点群領域Ｂ１，Ｂ２，Ｂ３，・・・が抽出されると、当該第２階層樹冠高点群領域Ｂ１，Ｂ２，Ｂ３，・・・に対して、局所最大値フィルタ（例えば、３メッシュ×３メッシュ）を用いて、個々の第２階層樹冠高点群領域Ｂ１，Ｂ２，Ｂ３，・・・のそれぞれにおける最大値を算定して、算定した最大値を個々の樹冠の樹頂点とする樹頂点算定処理（ステップＳ５４）を行う。
【００８７】
このとき、最大値を算定すべき階層（この場合、第２階層）において作成された個々の樹冠に対応する階層別樹冠高点群データのうち、階層別樹冠高点群データが既に作成されている樹冠と同じ樹冠に対応して作成された階層別樹冠高点群データは除外して、最大値を算定する。図７を例に取って説明すると、第２階層樹冠高点群領域のうち、例えば、樹冠高点群領域Ｂ４，Ｂ５などは、既抽出樹冠高点群領域であるとして、当該既抽出樹冠高点群領域Ｂ４，Ｂ５などは樹頂点算定対象とはしない。
【００８８】
このようにして、第２階層樹冠高点群領域Ｂ１，Ｂ２，Ｂ３，・・・において抽出された樹頂点ｂ１，ｂ２，ｂ３，・・・は、図７において黒点で示されている。なお、図７はモノクロ画像であるため、図７においては、第２階層樹冠高点群領域Ｂ１，Ｂ２，Ｂ３，・・・において抽出された樹頂点ｂ１，ｂ２，ｂ３，・・・及び既に抽出された第１階層樹冠高点群領域における樹頂点それぞれが黒丸の点で示されているが、両者の区別を明確にするために、図７の元となるカラー画像上では、既に抽出された階層（この場合、第１階層）の樹冠高点群領域における樹頂点は赤丸の点で示されている。
【００８９】
第２階層樹冠高点群領域において樹頂点ｂ１，ｂ２，ｂ３，・・・を抽出する処理が終了したら、Ｈｍａｘ（調査対象森林域内の最大樹冠高）−Ｄ（階層幅）×Ｎｉ（現時点の処理対象階層）＝Ｈｉ（次の処理対象階層内の最大樹冠高）を求める処理（ステップＳ５５）を行う。なお、次の処理対象階層内の最大樹冠高Ｈｉは、調査対象森林域内の最大樹冠高（Ｈｍａｘ）＝２５ｍ、階層幅（Ｄ）＝１ｍであり、また、この時点での処理対象階層Ｎｉ＝２であるため、２５ｍ−１ｍ×２＝２３ｍとなり、次の処理対象階層内の最大樹冠高Ｈｉは２３ｍとなる。そして、Ｈｉ＞０か否かを判定し（ステップＳ）、この場合、Ｈｉ＝２３であり、Ｈｉ＞０であるため、Ｎｉ＝Ｎｉ＋１＝２＋１を行って（ステップＳ５７）、Ｎｉ＝３として、処理対象階層を第２階層に設定して、階層別樹冠高点群データ作成処理Ｓ５３に戻る。
【００９０】
ここで、第３階層は、樹冠高が２２ｍ〜２３ｍの階層であり、当該第３階層において、階層別樹冠高点群データ作成処理（ステップＳ５３）を行う。すなわち、各メッシュの樹冠高のうち、当該第３階層に属する樹冠高点群データ（第３階層樹冠高点群データとする。）を作成する。そして、この第３階層樹冠高点群データに対応する第３階層樹冠高点群領域をディスプレイ上に表示する。
【００９１】
図８は、ディスプレイ上に表示された第３階層樹冠高点群領域Ｃ１，Ｃ２，Ｃ３，・・・を示す図である。なお、図８においては、第３階層樹冠高点群領域Ｃ１，Ｃ２，Ｃ３，・・・だけではなく、当該第３階層樹冠高点群領域Ｃ１，Ｃ２，Ｃ３，・・・において抽出された樹頂点も示されている。また、図８において、灰色で塗りつぶされている領域が第３階層樹冠高点群領域Ｃ１，Ｃ２，Ｃ３，・・・であり、当該第３階層樹冠高点群領域Ｃ１，Ｃ２，Ｃ３，・・・は、樹頂点付近が第３階層（２２ｍ〜２３ｍ）に達している樹木の存在を示している。例えば、図８において、黒枠内に含まれる樹木群の様子（横断面図）を見ると、これらの樹木群の中には、樹頂点付近が、帯状に示されている第３階層にまで達している樹木が存在していることがわかる。
【００９２】
なお、図８はモノクロ画像であるため、図８においては、帯状に示されている第３階層及び第３階層樹冠高点群領域Ｃ１，Ｃ２，Ｃ３，・・・は、それぞれ灰色で示されているが、他の階層において得られる各階層の樹冠高点群領域との区別を明確にするために、図８の元となるカラー画像上では、帯状に示されている第３階層及び第３階層樹冠高点群領域Ｃ１，Ｃ２，Ｃ３，・・・は、黄色で示されている。
【００９３】
なお、図８においては、新たに抽出された第２階層樹冠高点群領域と、既に抽出されている第１階層樹冠高点群領域Ａ１，Ａ２，Ａ３，・・・及び第２階層樹冠高点群領域Ｂ１．Ｂ２，Ｂ３，・・・との区別をし易くするために、既に抽出されている第１階層樹冠高点群領域Ａ１，Ａ２，Ａ３，・・・及び第２階層樹冠高点群領域Ｂ１．Ｂ２，Ｂ３，・・・は白抜きで示されている。
【００９４】
このようにして、第３階層樹冠高点群領域Ｃ１，Ｃ２，Ｃ３，・・・が抽出されると、当該第３階層樹冠高点群領域Ｃ１，Ｃ２，Ｃ３，・・・に対して、局所最大値フィルタ（例えば、３メッシュ×３メッシュ）を用いて、個々の第３階層樹冠高点群領域Ｃ１，Ｃ２，Ｃ３，・・・のそれぞれにおける最大値を抽出して、抽出した最大値を当該樹木の樹頂点とする樹頂点算定処理（ステップＳ５４）を行う。
【００９５】
この場合も、最大値を算定すべき階層（この場合、第３階層）において作成された個々の樹冠に対応する階層別樹冠高点群データのうち、階層別樹冠高点群データが既に作成されている樹冠と同じ樹冠に対応する階層別樹冠高点群データは除外して、最大値を算定する。第８図を例に取って説明すると、第３階層樹冠高点群領域のうち、例えば、樹冠高点群領域Ｃ４，Ｃ５などは、既抽出樹冠高点群領域であるとして、当該既抽出樹冠高点群領域Ｃ４，Ｃ５などは樹頂点算定対象とはしない。
【００９６】
このようにして、第３階層樹冠高点群領域Ｂ１，Ｂ２，Ｂ３，・・・において抽出された樹頂点ｃ１，ｃ２，ｃ３，・・・は、図８において黒丸の点で示されている。なお、図８はモノクロ画像であるため、図８においては、これまでのすべての階層の樹冠高点群領域において抽出された樹頂点それぞれが黒丸の点で示されているが、既に抽出された第１階層樹冠高点群領域における樹頂点と、第２階層及び第１階層樹冠高点群領域において抽出された樹頂点との区別を明確にするために、図８の元となるカラー画像上では、第２階層及び第１階層樹冠高点群領域において抽出された樹頂点は赤丸の点で示されている。
【００９７】
以上のようにして、第３階層樹冠高点群領域において樹頂点ｃ１，ｃ２，ｃ３，・・・を抽出する処理が終了したら、上述同様の処理（ステップＳ５５，Ｓ５６，Ｓ５７を行って、階層別樹冠高点群データ作成処理（ステップＳ５３）と、樹木頂点算定処理（ステップＳ５４）とを順次繰り返し行う。このような処理は、Ｈｉ＝０となるまで行ってもよいが、必ずしも、Ｈｉ＝０となるまで行う必要はなく、調査対象森林域の状況などに応じて適宜設定可能であり、例えば、Ｈｉ＝５（次の処理対処階層内の最大樹冠高が５ｍ）程度で終了してもよい。
【００９８】
図９は、第１階層から第５層までの各階層の階層別樹冠高点群領域と当該階層別樹冠高点群領域における樹頂点を示す図である。図９においては、各階層において抽出された階層別樹冠高点群領域（第１階層〜第５階層樹冠高点群領域）は、各階層別に異なった色で示されている。なお、図９はモノクロ画像であるため、各階層を色別として読み取ることは困難であるが、実際には、図９の元となるカラー画像上では、第１階層は赤色、第２階層は橙色、第３階層は黄色、第４階層は黄緑色、第５階層は緑色というように色別で示されている。そして、各階層において抽出された新規樹冠高点群領域における最大値（樹頂点）が黒丸の点で示されている。このような処理をさらに下位の階層まで順次行うことによって、階層ごとに個々の樹冠の樹冠高点群画像を得ることができる。
【００９９】
図１０は、階層別樹冠高点群データ作成処理（ステップＳ５３）及び樹頂点算定処理（ステップＳ５４）を第１階層（最上層）から最下層まで順次行うことによって得られた全階層の階層別樹冠高点群領域と樹頂点を示す図である。階層別樹冠高点群データ作成処理（ステップＳ５３）及び樹頂点算定処理（ステップＳ５４）を第１階層（最上層）から最下層まで順次行うことによって、調査対象森林域において第１階層（最上層）から最下層までの階層ごとの階層別樹冠高点群領域と調査対象森林域に存在するすべての樹冠の樹頂点を抽出することができる。図１０において、黒丸の点が階層ごとの個々の樹冠の樹頂点である。ここで、個々の樹冠の樹頂点は、当該樹冠を有する各樹木の樹頂点であるため、樹冠の樹頂点を得ることによって当該樹冠を有する樹木の樹高を得ることができる。各樹木の樹高は、平均化処理（ステップＳ５１）によってノイズが除去されたノイズ除去済みの樹冠高画像データから得ることができる。
【０１００】
図１０において、各階層において抽出された階層別樹冠高点群領域は、各階層別に異なった色で示されている。なお、図１０はモノクロ画像であるため、各階層を色別として読み取ることは困難であるが、実際には、図１０の元となるカラー画像上では、第１階層は赤色、第２階層は橙色、第３階層は黄色、第４階層は黄緑色、第５階層は緑色というように色別で示されている。
【０１０１】
このようにして、調査対象森林域に存在するすべての樹頂点が抽出されたら、図１のフローチャートにおける精密樹冠画像データ作成処理（ステップＳ６０）を行う。この精密樹冠画像データ作成処理は、各樹木に対応した樹冠を高精度に抽出可能とする精密樹冠画像データを作成する処理である。
【０１０２】
精密樹冠画像データ作成処理は、階層ごとにそれぞれ抽出された樹冠の樹頂点を基にして、当該樹冠の樹冠高から、隣接して存在する樹冠を領域分割する領域分割アルゴリズムを用いて、個々の樹冠が区分された精密樹冠画像データを作成する。具体的には、階層ごとにそれぞれ抽出された樹冠の樹頂点を種（Seed）にして、ｚ値（ステップＳ４０の樹冠高画像データ作成処理で得られた樹冠高画像データ）から、公知の領域分割法であるWatershedアルゴリズムを用いて精密な樹冠画像を作成するものである。なお、実施形態に係る森林資源情報算定方法においては、ｚ値としては、図４に示す平均化処理（ステップＳ５１）によってノイズが除去されたノイズ除去済みの樹冠高画像データを用いている。
【０１０３】
前述したように、ノイズ除去済みの樹冠高画像データは、個々の樹冠に存在する極端な凹凸（枝葉が周囲に比べて疎になっている部分や徒長枝葉など）が除去されて、個々の樹冠が平滑化されたものとなっている。このため、Watershedアルゴリズによる領域分割を行う際に、このようなノイズ除去済みの樹冠高画像データを用いることによって、Watershedアルゴリズによる領域分割を高精度に行うことができる。
【０１０４】
Watershedアルゴリズムは、隣接して密集する多数の物体同士が繋がった状態で存在する場合において、個々の物体の境界を高精度に抽出して個々の物体を区分させることができるアルゴリズムである。当該Watershedアルゴリズムについては、下記に示す公知文献２に記載されている。
【０１０５】
公知文献２：Salman N. and Liu C. Q., “Image Segmentation and Edge Detection Based on Watershed Techniques,” International Journal of Computers and Applications, vol. 25, no. 4, pp. 258-263,2003年
【０１０６】
精密樹冠画像データ作成処理（ステップＳ６０）は、Watershedアルゴリズムを用いて精密樹冠画像データを作成し、当該精密樹冠画像データ作成処理（ステップＳ６０）によって作成された精密樹冠画像データに対応する精密樹冠画像をディスプレイ上に表示する。
【０１０７】
図１１は、ディスプレイ上に表示された精密樹冠画像を示す図である。図１１に示すように、Watershedアルゴリズムを用いて作成された精密樹冠画像は、個々の樹木の樹頂点を取り囲むように個々の樹冠が領域分割されたものとなっている。図１１において、枠線Ｌで囲まれた領域が個々の樹冠である。
【０１０８】
図１１において、各階層において抽出された階層別樹冠高点群領域は、各階層別に異なった色で示されている。なお、図１１はモノクロ画像であるため、各階層を色別として読み取ることは困難であるが、実際には、図１１の元となるカラー画像上では、図１０と同様に、第１階層は赤色、第２階層は橙色、第３階層は黄色、第４階層は黄緑色、第５階層は緑色というように階層ごとに色別で示されている。
【０１０９】
図１１に示す精密樹冠画像からもわかるように、Watershedアルゴリズムによって領域分割されて得られた各樹木に対応する樹冠は、階層ごとに、それぞれの階層に属する樹木の樹冠が高精度に表わされたものとなる。上述したように、Watershedアルゴリズムは、隣接して密集する多数の物体同士が繋がった状態で存在する場合において、個々の物体の境界を高精度に抽出して個々の物体を区分させることができるアルゴリズムである。従って、Watershedアルゴリズムを用いることによって、密集して存在する樹冠を、樹冠の大きさ（樹冠直径や高さ）に係わらず個々の樹冠として高精度に区分することができる。
【０１１０】
このようにして、精密樹冠画像データが作成されると、当該精密樹冠画像データに基づいて精密樹冠情報作成処理（ステップＳ７０）を行い、各樹木の樹冠に関する精密な情報（精密樹冠情報とする。）を作成する。精密樹冠情報作成処理（ステップＳ７０）において作成される精密樹冠情報としては、少なくとも、各樹木に対応する樹冠のラベル番号（通し番号）、各樹木に対応する樹冠の樹冠位置（樹冠の中心位置）、各樹木に対応する樹冠の樹冠直径、各樹木に対応する樹冠の樹冠面積などを例示することができ、この他に、各樹木に対応する樹冠の樹冠形状なども精密樹冠情報として得ることができる。
【０１１１】
続いて、森林資源情報算定処理（ステップＳ８０）を行う。当該森林資源情報算定処理は、精密樹冠情報作成処理（ステップＳ７０）で作成された精密樹冠情報に基づいて、森林資源に関する情報を算定する。当該森林資源情報算定処理（ステップＳ８０）において算定する森林資源に関する情報としては、各樹木の胸高直径、各樹木の材積を例示できる。なお、森林資源に関する情報としては、上述の各樹木の胸高直径、各樹木の材積の他に、各樹木の樹高を挙げることができるが、当該樹高は、階層別樹頂点抽出処理（ステップＳ５０）の樹頂点算定処理（ステップＳ５４）で算定された樹頂点に基づいて、各樹木の樹高を得ることができる。
【０１１２】
これら森林資源に関する情報を算定するために、森林資源情報算定処理（ステップＳ８０）においては、各樹木の胸高直径を算定する胸高直径算定処理（ステップＳ８１）と、各樹木の資源としての材積を算定する材積算定処理（ステップＳ８２）と、森林資源情報を集計する森林資源情報集計処理（ステップＳ８３）とを行う。
【０１１３】
胸高直径算定処理（ステップＳ８１）は、胸高直径と強い相関のある樹冠面積及び樹冠直径の少なくとも一方と、樹高とから重回帰式で胸高直径を算定する。また、内挿法で未測定木の胸高直径を算定することができる。これによって、調査対象森林域に存在する調査対象樹木すべての胸高直径を算定することができる。なお、重回帰式を求める場合、樹種ごとに現地調査において標準木を１０数本程度選び、選んだ標準木の樹冠面積、樹冠直径及び樹高を測定することから重回帰式の変数を求めることが好ましい。
【０１１４】
材積算定処理（ステップＳ８２）は、胸高直径と樹高の２変数材積式から材積を求める。なお、胸高直径と樹高とから材積が求められている早見表が公知であり、このような早見表から木資源（材積）を求めることができる。当該早見表については、下記文献（公知文献３及び公知文献４）に記載されている。
【０１１５】
公知文献３：林野庁計画課、立木幹材積表  東日本編、日本林業調査会、３３４ページ、２００３年
公知文献４：林野庁計画課、立木幹材積表  西日本編、日本林業調査会、３２０ページ、１９７０年
【０１１６】
森林資源情報集計処理（ステップＳ８３）は、調査対象樹木の本数、各樹木についての樹冠直径、樹冠面積、胸高直径、樹高及び材積などの各樹木の資源情報、すなわち、単木資源情報を集計する。そして、集計した単木資源情報を属性データベースに登録する。属性データベースに登録されている単木資源情報は、一覧表としてディスプレイ上に表示することができる。なお、一覧表として表示される単木資源情報には、樹冠の中心位置を表すＸＹ座標も含まれている。また、調査対象樹木の本数は、精密樹冠情報作成処理（ステップＳ７０）によって得られたラベル番号に基づいて、当該調査対象森林域における調査対象樹木の本数を求めることができる。
【０１１７】
ところで、上述した「樹冠直径」は、樹冠を平面視したときの樹冠の水平方向の広がりの大きさを表すものであり、また、「胸高直径」は、樹木の幹の太さ（各樹木において人間のほぼ胸の高さの直径）を表すものである。なお、樹冠及び樹木の幹は実際には真円ではないため、厳密には「直径」とは言えないが、この明細書においては、樹冠の水平方向の広がりの大きさを表すものとして「樹冠直径」と表記し、幹の太さを表すものとして「胸高直径」と表記するものとする。
【０１１８】
図１２は、ディスプレイ上に表示された単木資源情報の一覧表を示す図である。図１２に示すように、調査対象森樹木の単木資源情報として、各樹木に対応する樹冠のナンバー（ラベル番号）と、個々の樹冠の中心位置を表すＸＹ座標、樹冠直径、樹冠面積、胸高直径（図１２においては「ＤＢＨ」と表記されている。）、樹高、材積が、樹冠のナンバー（ラベル番号）に対応付けられて表示される。
【０１１９】
また、森林資源情報集計処理（ステップＳ８３）は、調査対象樹木の本数とともに、各樹木に対して算定された樹冠直径の平均値、各樹木に対して算定された樹冠直径のうちの最大値及び最小値、各樹木に対して算定された樹冠面積の平均値、各樹木に対して算定された樹冠面積のうちの最大値及び最小値、各樹木に対して算定された胸高直径の平均値、各樹木に対して算定された胸高直径の平均値、各樹木に対して算定された胸高直径のうちの最大値及び最小値、各樹木に対して算定された樹高の平均値、各樹木に対して算定された樹高のうちの最大値及び最小値と、各樹木に対して算定された材積の平均値、各樹木に対して算定された材積のうちの最大値及び最小値を算定するとともに、材積については合計値をも算定し、さらに、森林経営で使用する１ヘクタール（ｈａ）当りの森林資源量（ここでは、樹木の本数及び材積とする。）を算定し、これら算定した結果を、森林資源の概要を示す情報（森林資源概要情報）として属性データベースに登録する処理も行う。属性データベースに登録されている森林資源概要情報は、一覧表としてディスプレイ上に表示することができる。
【０１２０】
図１３は、ディスプレイ上に表示された森林資源概要情報の一覧表を示す図である。図１３に示すように、調査対象樹木の本数、樹冠直径の平均値、樹冠直径の最大値及び最小値、樹冠面積の平均値、樹冠面積の最大値及び最小値、胸高直径の平均値、胸高直径の平均値、胸高直径の最大値及び最小値、樹高の平均値、樹高の最大値及び最小値と、材積の平均値、材積の最大値及び最小値、材積の合計値が表示され、さらに、森林経営で使用する１ヘクタール（ｈａ）当りの森林資源量（ここでは、樹木の本数及び材積とする。）が表示されている。ディスプレイ上にこのような森林資源概要情報を表示することによって、調査対処森林域の森林資源の概要を把握することができる。
【０１２１】
なお、この明細書においては、図１２に示す単木資源情報と図１３に示す森林資源概要情報とを合わせて森林資源情報とする。このような森林資源情報と既存の森林調査簿データベース（例えば、都道府県などで整備されている森林の管理台帳などのデータベース）とを組み合わせることで、森林現況の把握及び間伐などの森林施業立案などの森林管理に活用できる。例えば、森林の間伐を行うおうとする際に、間伐計画のシミュレーションを立てる場合に役立つものとなる。
【０１２２】
以上説明したように、実施形態に係る森林資源情報算定方法においては、調査対象森林域を含む地域をドローンなどによって撮影した空撮画像データから３次元点群データを作成したのち、当該３次元点群データと地理的情報とに基づいて調査対象森林域画像データを作成し、作成された調査対象森林域画像データに基づいて、調査対象森林域に存在する樹木の樹冠高を表す樹冠高画像データを作成する。そして、平均化処理によって樹冠高画像データからノイズ除去したノイズ除去済みの樹冠高画像データを作成したのち、当該ノイズ除去済みの樹冠高画像データに基づいて、第１階層を含む各階層に存在する個々の樹冠に対応する階層別樹冠高点群データを階層ごとに作成し、階層ごとに作成された個々の樹冠に対応する階層別樹冠高点群データから、個々の樹冠に対応する階層別樹冠高点群データの最大値を算定するようにしている。
【０１２３】
ここで、階層別樹冠高点群データを階層ごとに作成し、階層ごとに作成された個々の樹冠に対応する階層別樹冠高点群データから、個々の樹冠に対応する階層別樹冠高点群データの最大値を算定する処理は、例えて言えば、個々の樹冠をＣＴ（Computed Tomography）スキャンのように、階層ごとに輪切りにして階層別樹冠高点群データを作成して、階層ごとに求められた階層別樹冠高点群データから最大値を算定するといった処理である。このような処理を行うことによって、各階層に存在する個々の樹冠の最大値を、階層ごとに算定できる。なお、階層ごとに算定された「樹冠の樹頂点」は、当該樹冠を有する樹木の樹頂点であり、当該算定された樹頂点から当該樹木の樹高を得ることができる。
【０１２４】
このようにして、各階層に存在する個々の樹冠の最大値（樹頂点）を、階層ごとに算定した後に、階層ごとにそれぞれ抽出された樹冠の樹頂点を基（種）にして、各樹冠の樹冠高から、領域分割アルゴリズム（watershedアルゴリズム）を用いて、隣接する樹冠が個々の樹冠に区分された精密樹冠画像データを作成するようにしている。
【０１２５】
このようにして作成された精密樹冠画像データは、調査対象樹木の各樹木に対応した個々の樹冠を高精度に抽出可能な精密樹冠画像データとなる。このため、当該精密樹冠画像データに基づいて作成した精密樹冠情報は、調査対象樹木の各樹木に対応した高精度な樹冠情報となり、当該精密樹冠情報に基づいて、森林資源に関する情報（森林資源情報）を算定することにより、高精度な森林資源情報を得ることができる。これにより、計測誤差が少なく客観性の高い高精度な森林資源情報を、調査対象森林域の全域、小班ごと、任意の範囲において提供できることから、間伐などの森林施業立案などの森林管理に活用できる。
【０１２６】
また、実施形態においては、空撮画像データを得るための撮影機材としては一般的なカメラと、普及型のドローンを使用することができるため、森林資源算定を行うためのシステムを安価でかつ容易に実現することができる。
【０１２７】
また、調査対象森林域をカメラで撮影して得られた空撮画像データを用いているため、樹冠の下側（樹冠と地面との間）に存在する下層植生（調査対処としない低木や草本など）が画像データとして取り込まれにくいため、下層植生をノイズとして除去するノイズ処理を削減できるといった効果も得られる。
【０１２８】
以上、実施形態に係る森林資源情報算定方法１について説明したが、続いて、実施形態に係る森林資源情報算定装置について説明する。
【０１２９】
図１４は、実施形態に係る森林資源情報算定装置１を説明するために示す図である。実施形態に係る森林資源情報算定装置１は、図１及び図４に示した各ステップの処理を行うための装置であり、図１におけるステップＳ１０の処理を行う機能を有する空撮データ入力部１０と、図１におけるステップＳ２０の処理を行う機能を有する３次元点群データ作成部２０と、図１におけるステップＳ３０の処理を行う機能を有する調査対象森林域画像データ作成部３０と、図１におけるステップＳ４０の処理を行う機能を有する樹冠高画像データ作成部４０と、図１におけるステップＳ５０の処理を行う機能を有する階層別樹頂点抽出部５０と、図１におけるステップＳ６０の処理を行う機能を有する精密樹冠画像データ作成部６０と、図１におけるステップＳ７０の処理を行う機能を有する精密樹冠情報作成部７０と、図１におけるステップＳ８０の処理を行う機能を有する森林資源情報算定部８０と、を有する。
【０１３０】
また、階層別樹頂点抽出部５０には、複数の処理を行うための構成要素が含まれており、これについては図１５により後述する。また、森林資源情報算定部８０には、図１におけるステップＳ８１の処理を行う機能を有する胸高直径算定部８１と、図１におけるステップＳ８２の処理を行う機能を有する材積算定部８２と、図１におけるステップＳ８３の処理を行う機能を有する森林資源情報集計部８３と、が含まれている。
【０１３１】
図１５は、図１４に示されている階層別樹頂点抽出部５０の構成を示す図である。階層別樹頂点抽出部５０には、図１５に示すように、図４におけるステップＳ５１の平均化処理を行うための平均化処理部５１と、図４におけるステップＳ５２の初期値設定処理を行うための初期値設定部５２と、図４におけるステップＳ５３の階層別樹冠高点群データ作成処理を行うための階層別樹冠高点群データ作成部５３と、図４におけるステップＳ５４の樹木頂点算定処理を行うための樹木頂点算定部５４と、図４におけるステップＳ５５のＨｍａｘ−Ｄ×Ｎｉ＝Ｈｉ、ステップＳ５６のＨｉ＞０、ステップＳ５７のＮｉ＝Ｎｉ＋１といった演算を行う演算処理部５８とが含まれている。
【０１３２】
このように、図１４及び図１５に示す実施形態に係る森林資源情報算定装置１の各構成要素は、前述した実施形態に係る森林資源情報算定方法の各ステップの処理（図１及び図４参照。）を行う機能を有しているため、実施形態に係る森林資源情報算定装置１は、前述した実施形態に係る森林資源情報算定方法で得られる効果と同様の効果が得られる。なお、実施形態に係る森林資源情報算定装置１の処理の内容についての説明は省略する。
【０１３３】
また、実施形態に係る森林資源情報算定装置１は、当該当該森林資源情報算定装置１に含まれる上記各構成要素（図１４及び図５参照。）それぞれが有する機能がコンピュータのプログラムとしてインストールされており、上記各構成要素に所定のデータを与えることによって、各構成要素それぞれが持つ機能がコンピュータのソフトウエア上で実行されるものである。
【０１３４】
このようなソフトウエアは、調査対象森林域を含む地域に対応する３次元点群データと地理的情報とに基づいて調査対象森林域画像データを作成する調査対象森林域画像データ作成処理ステップと、調査対象森林域画像データに基づいて、調査対象森林域に存在する各樹木の樹冠の高さを樹冠高として表す樹冠高画像データの作成を行う樹冠高画像データ作成処理ステップと、樹冠高画像データから得られる樹冠高のうち、最大樹冠高が存在する階層を第１階層として、前記最大樹冠高を起点に所定長さの階層幅を有する階層を前記第１階層から地表面に向かって複数設定したときに、前記第１階層を含む各階層に樹頂点が存在する樹冠の樹頂点を、前記階層ごとにそれぞれ抽出する階層別樹頂点抽出処理ステップと、を有するように構成してもよい。
【０１３５】
なお、本発明は上述の実施形態に限られるものではなく、本発明の要旨を逸脱しない範囲で種々変形実施可能となるものである。例えば、下記に示すような変形実施も可能である。
【０１３６】
（１）上記実施形態においては、調査対象樹木は、例えば、カラマツ、アカマツ、スギ、ヒノキなどの植栽木とした場合を例示したが、調査対象樹木としては、上記した植栽木に限られるものではなく、例えば、ブナやナラなどが所定の広さの範囲にまとまって群生している林が存在する場合においては、当該ブナやナラを調査対象樹木とすることも可能である。
【０１３７】
（２）上記実施形態においては、ドローンを用いて空撮画像データを得るようにしたが、ドローンに限られるものではなく、例えば、航空機を使用して空撮画像データを得るようにしてもよい。
【０１３８】
（３）上記実施形態においては、３次元形状復元技術であるＳＦＭ(Structure from Motion)によって、空撮画像データから歪みのない地図補正されたオルソ画像を作成して、作成したオルソ画像から３次元点群データを作成するようにしたが、上空からレーザー光を照射して得られたレーザー計測データに基づいて３次元点群データを得るようにしてもよい。
【符号の説明】
【０１３９】
１・・・森林資源情報算定装置、１０・・・空撮データ入力部、２０・・・３次元点群作成部、３０・・・調査対象森林域画像データ作成部、４０・・・樹冠高画像データ作成部、５０・・・階層別樹頂点抽出部、５１・・・平均化処理部、５２…初期値設定部、５３・・・階層別樹冠高点群データ作成部、５４・・・樹頂点算定部、５８・・・演算処理部、６０・・・精密樹冠画像データ作成部、７０・・・精密樹冠情報作成部、８０・・・森林資源情報算定部、８１・・・胸高直径算定部、８２・・・材積算定部、８３・・・森林資源情報集計部、Ｓ１０・・・空撮データ入力処理、Ｓ２０・・・３次元点群作成処理、Ｓ３０・・・調査対象森林域画像データ作成処理、Ｓ４０・・・樹冠高画像データ作成処理部、Ｓ５０・・・階層別樹頂点抽出処理、Ｓ５１・・・平均化処理、Ｓ５２…初期値設定処理、Ｓ５３・・・階層別樹冠高点群データ作成処理、Ｓ５４・・・樹頂点算定処理、Ｓ５５・・・Ｈｍａｘ−Ｄ×Ｎｉ＝Ｈｉの演算処理、Ｓ５６・・・Ｈｉ＞０の演算処理、Ｓ５７・・・Ｎｉ＝Ｎｉ＋１の演算処理、Ｓ６０・・・精密樹冠画像データ作成処理、Ｓ７０・・・精密樹冠情報作成処理、Ｓ８０・・・森林資源情報算定処理、Ｓ８１・・・胸高直径算定処理、Ｓ８２・・・材積算定処理、Ｓ８３・・・森林資源情報集計処理、Ａ１，Ａ２，Ａ３，・・・第１階層樹冠高点群領域、ａ１，ａ２，ａ３・・・第１階層樹冠高点群領域の最大値、Ｂ１，Ｂ２，Ｂ３，・・・第２階層樹冠高点群領域、ｂ１，ｂ２，ｂ３・・・第２階層樹冠高点群領域の最大値、Ｃ１，Ｃ２，Ｃ３，・・・第３階層樹冠高点群領域、ｃ１，ｃ２，ｃ３・・・第３階層樹冠高点群領域の最大値、Ｈｍｚｘ・・・調査対象森林域内の最大樹冠高、Ｄ・・・階層幅、Ｎｉ・・・現時点の処理対象階層、Ｈｉ・・・次の処理対象階層内の最大樹冠高
【図１】
2019185449.tif 000003
【図２】
2019185449.tif 000004
【図３】
2019185449.tif 000005
【図４】
2019185449.tif 000006
【図５】
2019185449.tif 000007
【図６】
2019185449.tif 000008
【図７】
2019185449.tif 000009
【図８】
2019185449.tif 000010
【図９】
2019185449.tif 000011
【図１０】
2019185449.tif 000012
【図１１】
2019185449.tif 000013
【図１２】
2019185449.tif 000014
【図１３】
2019185449.tif 000015
【図１４】
2019185449.tif 000016
【図１５】
2019185449.tif 000017