(57)【要約】【課題】より安定して精度よく樹木の頂点を特定するための樹木幅設定装置及びプログラムを提供する。【解決手段】樹木幅設定装置は、数値表層モデル及び数値標高モデルを用いて抽出された樹木の頂点候補の各々における樹木の特性に基づいて、当該樹木の領域幅を設定する幅設定部を備える。【選択図】図４
2020098442.tif 000002
【特許請求の範囲】
【請求項１】
数値表層モデル及び数値標高モデルを用いて抽出された樹木の頂点候補の各々における樹木の特性に基づいて、当該樹木の領域幅を設定する幅設定部を備えることを特徴とする樹木幅設定装置。
【請求項２】
前記特性には、前記頂点候補の位置における前記樹木の高さが含まれることを特徴とする請求項１記載の樹木幅設定装置。
【請求項３】
前記頂点候補の位置のそれぞれにおける樹木種を特定する種別特定部を備え、
前記種別特定部は、樹木種の分布情報を取得し、当該分布情報における前記頂点候補の位置を含む領域の樹木種の情報により当該樹木種を特定し、
前記特性には、特定された前記樹木種が含まれる
ことを特徴とする請求項１又は２記載の樹木幅設定装置。
【請求項４】
前記頂点候補の位置の環境特性情報を特定する特性特定部を備え、
前記環境特性情報には、前記頂点候補の位置の地理的特性が含まれ、
前記幅設定部は、前記環境特性情報に基づいて前記領域幅を設定する
ことを特徴とする請求項１〜３のいずれか一項に記載の樹木幅設定装置。
【請求項５】
コンピュータを、
数値表層モデル及び数値標高モデルを用いて抽出された樹木の頂点候補の各々における樹木の特性に基づいて、当該樹木の領域幅を設定する幅設定手段
として機能させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、樹木幅設定装置及びプログラムに関する。
【背景技術】
【０００２】
従来、センシングデータを用いて樹木の頂点を特定する技術がある。特定された頂点は、例えば、樹木の本数の計数に用いられたり、各樹木の伐採計画の立案に設けられたりといった森林資源の把握や利用に係る各種用途で用いられる。
【０００３】
樹木の頂点位置を特定する手法の一つとして、樹冠の上縁に係る地面からの高さの分布を取得し、当該分布における高さの空間的な極大点を抽出することで樹木の頂点を特定する技術がある。高さの分布を示す樹冠高データ（Digital Canopy Height Model、ＤＣＨＭ）は、樹木表面を含む高さ分布を示す数値表層モデル（Digital Surface Model、ＤＳＭ）と、地表面の高さ分布を示す数値標高モデル（Digital Elevation Model、ＤＥＭ）との差分をとることで得られる。これらＤＳＭモデル及びＤＥＭモデルを得るための計測手法の一つとして、上空から照射したレーザ光の反射を用いたレーザプロファイラが知られている。
【０００４】
このとき、樹木の枝ぶりなどによっては、頂点が適切に検出できず、過小又は過大に特定される場合がある。特許文献１では、ＤＣＨＭデータから地上開度及び地下開度の分布を求め、これらの値に基づいて予め定められた樹木間隔の範囲内でそれぞれ頂点を検出する技術が開示されている。樹木間隔は、例えば、ユーザにより別途設定される。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００９−０２２２７８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、森林における立木密度や樹木間の間隔は、各樹木について一様ではない。したがって、従来の取り扱いでは、全体として樹木の頂点の特定精度が低くなりやすいという問題がある。
【０００７】
この発明の目的は、より安定して精度よく樹木の頂点を特定するための樹木幅設定装置及びプログラムを提供することにある。
【課題を解決するための手段】
【０００８】
上記目的を達成するため、本発明は、
数値表層モデル及び数値標高モデルを用いて抽出された樹木の頂点候補の各々における樹木の特性に基づいて、当該樹木の領域幅を設定する幅設定部を備えることを特徴とする樹木幅設定装置である。
【発明の効果】
【０００９】
本発明に従うと、より安定して精度よく樹木の頂点を特定することが可能になるという効果がある。
【図面の簡単な説明】
【００１０】
【図１】処理装置の機能構成を示すブロック図である。
【図２】ＤＳＭ、ＤＥＭ及びＤＣＨＭを説明する模式図である。
【図３】樹木分布情報及び地理情報を模式的に示す図である。
【図４】樹木特定処理の流れの概略を説明する図である。
【図５】林分密度テーブルの一部の例を示す図表である。
【図６】頂点候補のリストの例を示す図表と、分布の例を示す図である。
【図７】樹木頂点特定処理の制御手順を示すフローチャートである。
【図８】樹木の本数を特定した結果の例を示す図表である。
【発明を実施するための形態】
【００１１】
以下、本発明の実施の形態を図面に基づいて説明する。
図１は、樹木幅設定装置の実施形態である処理装置１の機能構成を示すブロック図である。
【００１２】
処理装置１は、例えば、通常のコンピュータ（コンピュータ端末又はサーバ）であり、制御部１１と、記憶部１２と、入出力インターフェイス１３と、表示部１４と、操作受付部１５などを備える。
【００１３】
制御部１１は、各種演算処理を行うプロセッサであり、例えば、ＣＰＵ（Central Processing Unit）などを備える。ＣＰＵは、記憶部１２からプログラムなどを読み込んで実行することで各種制御処理を行う。
【００１４】
記憶部１２は、ＲＡＭと、不揮発性メモリなどを備え、各種データを記憶する。ＲＡＭは、ＣＰＵに作業用のメモリ空間を提供し、一時データを記憶する。不揮発性メモリは、プログラム１２１や設定データなどを記憶保持する。不揮発性メモリとしては、フラッシュメモリが含まれ、あるいは、これに代えて又は加えてＨＤＤ（Hard Disk Drive）などが含まれてもよい。また、初期制御プログラムなどは、ＲＯＭ（Read Only Memory）に格納されていてもよい。プログラム１２１には、後述の樹木特定処理プログラムが含まれる。
【００１５】
入出力インターフェイス１３は、周辺機器と接続するための接続端子１３１や外部機器との間で信号の送受信を行うための通信部１３２などを備える。周辺機器としては、大型の補助記憶装置、例えば、データベース装置２１、並びにＣＤＲＯＭ、ＤＶＤ及びＢｌｕ−ｒａｙ（登録商標）などの可搬型記憶媒体（光ディスク）を読み取る光学読取装置２２などが挙げられる。また、可搬型記憶媒体に磁気テープが含まれ、この磁気テープを読み取る読取装置が周辺機器に含まれてもよい。接続端子としては、ＵＳＢ（Universal Serial Bus）など各種規格に応じたものが利用されてよい。あるいは、ＬＡＮ（Local Area Network）などのネットワーク回線を介して通信部が外部から取得するものであってもよい。データベース装置２１及び／又は可搬型記憶媒体に記憶されたデータには、ＤＳＭデータ２０１、ＤＥＭデータ２０２、樹木分布情報２０３（分布情報）、地理情報２０４（環境特性情報）及び林分密度テーブル２０５が含まれる。
【００１６】
ＤＳＭデータ２０１及びＤＥＭデータ２０２は、水平面に平行な２次元面（２次元地図面）について、２次元地図を区切った複数のメッシュごとに中心位置の水平座標と高度（標高値）とを対応付けたものである。
【００１７】
図２は、ＤＳＭ、ＤＥＭ及びＤＣＨＭを説明する模式図である。
図２（ａ）に示すように、ＤＳＭの高度は、地上の構造物や樹林などを含めたものであり、樹林帯では、樹冠の上縁に係る高度分布を示す。これらの値は、例えば、レーザプロファイラによって計測、取得されたものである。この高度には、後述のように、送電線Ｗなども含まれ得る。また、葉や枝の隙間を抜けて地面や下草の高度が取得されたものが混じっている場合がある。後述のように、樹木の高さとして適切でない高さのデータを除外するために、上限基準高さＨｍａｘや下限基準高さＨｍｉｎが定められている。水平座標は、緯度経度の値であってもよいし、所定の原点位置からの相対距離を示すものであってもよい。
【００１８】
また、ＤＥＭでは、水平面に平行な２次元面について、２次元地図を区切った複数のメッシュごとに地表面の高度を水平座標と対応付けたものである。すなわち、ＤＥＭデータ２０２では、レーザプロファイラなどで取得されたＤＳＭデータから地上の構造物や樹林などの高さを取り除いた高度値（標高）の分布を示す。水平座標は、ＤＳＭデータ２０１と同様に表される。隣り合う樹木の間の距離が樹木間距離Ｄである。
【００１９】
図２（ｂ）に示すように、樹林帯では、同一水平座標におけるＤＳＭデータ２０１の値とＤＥＭデータ２０２の値との差が樹冠上縁の高さＨ（すなわち、各位置における樹木の高さ）となり、その分布がＤＣＨＭである。樹木又は構造物がなく、直接地表が観測される部分では、ＤＳＭデータ２０１の値とＤＥＭデータ２０２の値は等しくなる。
【００２０】
樹木分布情報２０３は、樹木種の分布に係る情報である。樹木分布情報２０３では、例えば、単独の樹木種又は優勢な樹木種の情報が単位位置範囲ごとにメッシュ状に規定されている。樹木分布情報２０３を参照することで、指定された緯度経度の情報に対して、対応する位置における樹木種Ｓが得られる。樹木分布情報２０３には、予め実測などにより各樹木種の生育範囲が特定されている林相図が用いられる。あるいは、樹木分布情報２０３は、後述の頂点候補の抽出処理と並行で樹木種を特定する処理がなされてもよい。例えば、ＤＳＭデータ２０１及びＤＥＭデータ２０２に係るレーザ計測と同時になされた写真撮影データに基づく樹木種の特定処理がなされてもよい。
【００２１】
地理情報２０４は、上記樹木種Ｓの分布よりも大きいエリアであるエリアＧの範囲情報（地理的特性）を保持する。エリアＧは、少なくとも、気候や積雪などの環境条件が異なる範囲、例えば、北海道、東北の日本海側、東北の太平洋側といった範囲で定められる。一般的には、エリアＧは、例えば、都道府県単位や市町村単位で区分されていればよい。
【００２２】
図３は、樹木分布情報２０３及び地理情報２０４を模式的に示す図である。
図３（ａ）の樹木分布情報２０３において、単位位置範囲Ｂは、例えば、矩形（正方形）形状であり、大きさは、上記のＤＳＭデータ２０１などにおける単位面積と同一であってもよいし、当該単位面積の所定倍などであってもよい。林相図では、樹木種の分布範囲ごとにエリアＳ１、Ｓ２、Ｓ３などが定められている。空白部分は、森林以外の部分である。これらのうち、特に単一樹種のエリアについては、具体的な樹種が明示されている。なお、メッシュ単位では樹木種がそれぞれＩＤ値などで示され、別途ＩＤ値などと樹木種の名称とが対応付けられていてもよい。また、最終的に樹木種に応じた領域幅の算出が可能であればよいので、処理装置１内で理解可能な対応付けがなされていれば（ＩＤ値のみで全ての処理が可能であれば）、樹木種の名称が明示されている必要はない。樹木種としては、具体的には、ヒノキ、スギ、アカマツなどが挙げられる。複数種類の雑木林などについては、より広く広葉樹林又は針葉樹林などで定められてもよい。例えば、ＤＳＭ、ＤＥＭ、ＤＣＨＭなどにおける水平位置の緯度経度情報に応じて指定された地点Ｐ１は、エリアＳ２に属し、当該エリアＳ２に対応する樹木種の樹木であると特定される。
【００２３】
図３（ｂ）の地理情報２０４において、エリア境界ＢＲが定められており、当該エリア境界ＢＲを境に両側がそれぞれエリアＧ１、Ｇ２として設定されている。ＤＳＭ、ＤＥＭ、ＤＣＨＭなどにおける水平位置情報に基づいて、当該位置に対応する単位位置範囲が所属するエリア情報が得られる。単位位置範囲は、樹木分布情報２０３の単位位置範囲Ｂと同一であってもよいし、より細かくてもよい。あるいは、境界に沿った部分だけ他の部分よりも細かく定められるなどであってもよい。例えば、上記地点Ｐ１は、エリアＧ１に属すると特定される。各エリアＧは、後述の最多立木密度に応じた値を算出するときに対応付けが可能に区分されていればよく、直接都道府県名が指定されていてもよいし、エリアごとにＩＤが付されていてもよい。
【００２４】
林分密度テーブル２０５は、林分の樹木種Ｓ、エリアＧ及び上層樹高の平均高さに応じた最多立木密度を求めるための算出パラメータが、樹木種Ｓ及びエリアＧの組み合わせに対応付けられたテーブルである。林分密度テーブル２０５については、後述する。
【００２５】
表示部１４は、制御部１１の制御により、表示画面にメニュー、ステータスや処理結果などの表示を行わせる。表示画面としては、特には限られないが、例えば、液晶表示画面を有する。操作受付部１５としてタッチパネルを有する場合には、表示画面は、タッチパネルと併用され、タッチセンサと重ねて設けられる。
【００２６】
操作受付部１５は、例えば、キーボードやマウスなどの入力デバイスを備え、ユーザなど外部からの入力操作を受け付ける。入力デバイスには、タッチパネル（タッチセンサ）などが含まれていてもよい。
【００２７】
次に、樹木特定処理の原理について説明する。樹木特定処理は、定められたエリア内における樹木の頂点位置を特定する処理である。
【００２８】
図４は、樹木特定処理の流れの概略を説明する図である。
ここでは、まず、数値表層モデル（ＤＳＭ）と、地表面の高さ分布を示す数値標高モデル（ＤＥＭ）との差分をとることで樹冠高データ（ＤＣＨＭ）を取得する（ＰＲ１）。当該ＤＣＨＭデータから高さの極大点を樹木の頂点候補として抽出する（ＰＲ２；候補抽出部）。抽出方法としては、例えば、従来の局所的な最大値（極大値）を取得する方法などが用いられてよい。
【００２９】
このとき、極大点は、枝ぶりなどにより樹頂点以外の場所でも生じ得る。後述の処理では、抽出された頂点候補から頂点ではないと判断された候補点が除去されていくので、抽出は、検出漏れがある過小検出よりも誤検出がある過剰検出の方がよい。したがって、局所的最大値の決定範囲は、実際に想定され得る樹頂点間の最小間隔に比して狭く定められるとよい。例えば、人工林では、多くの場合、植林時の植林密度には上限（例えば、４０００本／ｈａなど）があり、以後は伐採されたり枯れたりして単調減少する。したがって、この植林密度に基づいて、すなわち、植林密度での樹木を分離可能な値に最小間隔が設定される。ここでは、例えば、１２０ｃｍ四方又は直径１２０ｃｍ（半径６０ｃｍ）の円などに設定される。
【００３０】
樹木分布情報２０３が参照（取得）されて、頂点候補として抽出された位置を含むエリア（領域）の樹木種Ｓが特定され（ＰＲ３；種別特定部）、また、地理情報２０４が参照されて、当該位置のエリアＧが特定される（ＰＲ４；特性特定部）。
【００３１】
各頂点候補の位置における樹頂点間の最小間隔である樹木間距離Ｄが求められ（ＰＲ５）、これにより、頂点候補の領域幅Ｒ（候補領域）が定められる。樹木の頂点を中心として領域幅Ｒを直径とする円形エリア内は、他の樹木の円形エリアと重複しない排他エリアとして設定される。ここでは、対象とされる頂点候補の位置における樹木の高さＨと、当該樹木が含まれるエリアの樹木種Ｓと、地理的なエリアＧとに基づいて、最多立木密度を求める式を用いて当該最多立木密度に当たる数値Ｎが求められ、数値Ｎから樹木間距離Ｄが求められる。上述のように、最多立木密度は、同様の高さの樹木が並んでいる樹林（林分）の上層樹高の平均的な高さに基づいて当該林分に対して求められる値である。ここでは、各頂点候補について、高さＨにより各々数値Ｎが算出される。すなわち、求める樹木間距離Ｄは、頂点候補ごとに各々別個に定められる変数である。領域幅Ｒは、樹木間距離Ｄと等しいものとして設定される。
【００３２】
数値Ｎの算出情報は、林分密度テーブル２０５から取得される。最多立木密度は、樹木種Ｓと、地理的なエリアＧと、上層樹高の平均値Ｈａとに依存する。ここでは、日本林業技術協会による人工林林分密度管理図などが用いられる。最多立木密度Ｎｒは、以下の（数式１）により定められている。
ｌｏｇ（Ｎｒ）＝ａ−ｂ×ｌｏｇ（Ｈａ）  …  （数式１）
係数ａ、ｂは、樹木種Ｓ及びエリアＧにより定まる数値パラメータである。すなわち、林分密度テーブル２０５には、樹木種Ｓ及びエリアＧの組み合わせに応じた係数ａ、ｂの値がそれぞれ設定、記憶されている。
【００３３】
図５は、林分密度テーブル２０５の一部の例を示す図表である。林分密度テーブル２０５では、樹木種Ｓ、エリアＧ、及びこれらに対応する係数ａ、ｂが対応付けられてそれぞれリスト記憶されている。これらの値は、上記の人工林林分密度管理図の他に、特定の地域や諸外国で定められた独自の又は第三者作成のデータが用いられてもよい。この場合には、樹木種Ｓ及びエリアＧ以外の要素が数値パラメータ（ここでは、係数ａ、ｂ）の特定に用いられてもよい。
【００３４】
特定されている樹木種Ｓ及びエリアＧにより、林分密度テーブル２０５からこれらに応じた係数ａ、ｂが取得される。これらの係数ａ、ｂとともに、上層樹高の平均値Ｈａの代わりに高さＨを用いて、数値Ｎが算出される。数値Ｎの平方根の逆数（Ｎ）−１／２が最小の樹木間距離Ｄ、すなわち、領域幅Ｒとなる。なお、係数ａ、ｂ及び高さＨから直接領域幅Ｒが求められるように計算式が予め定められていてもよい。（数式１）で分かるように、高さＨが大きいほど数値Ｎが低下し、これに応じて領域幅Ｒが広くなる。
【００３５】
上記のように抽出された頂点候補の位置情報（すなわち、エリアＧ）及び高さＨに基づいて、当該樹木の周囲に他の樹木の候補領域と重複しない範囲（候補領域）を示す領域幅Ｒが得られる。この領域幅Ｒは、同一エリアの同一樹木種であっても、高さＨに依存することから、本実施形態の処理装置１では、各頂点候補に対して領域幅Ｒが各々異なる値となる場合がある。
【００３６】
図６は、このようにして得られた頂点候補のリストの例を示す図表（ａ）と、分布の例を示す図（ｂ）である。
図６（ａ）に示すように、緯度経度で示される頂点候補に対し、それぞれ、その高さＨ及び樹木種Ｓ（樹木の特性）、並びにエリアＧ（地理的特性）が定められている。これらに基づいて得られた樹木間距離Ｄに応じて、頂点候補が頂点である場合の領域幅Ｒが設定されて、頂点候補に係る領域（候補領域）として対応付けられる（ＰＲ６；幅設定部（ＰＲ５を含んでよい））。
【００３７】
この候補領域の設定データに基づき、頂点が特定される（ＰＲ７；樹頂点特定部）。
図６（ｂ）に示すように、各頂点候補（Ｐ１１〜Ｐ１４）を中心として領域幅（Ｒ１１〜Ｒ１４）をそれぞれ直径とする範囲（候補領域）を設定した場合に、当該候補領域が他の頂点候補について定められた候補領域と重複するか否かを検出する。ここでは、頂点候補Ｐ１１の候補領域と頂点候補Ｐ１２の候補領域とが重複している。重複する場合には、いずれかの頂点候補が誤抽出されたものであると判断される。したがって、重複範囲がなくなるように、いずれかの頂点候補が除去されて、残った頂点候補が頂点として特定される。
【００３８】
重複がある場合に、いずれの頂点候補が除去されるかは、機械的に行われてもよいし、ユーザが判断して手動で決定されてもよい。機械的に行われる場合には、例えば、重複している複数点のうち高さＨが最大の点が選択される（ここでは、頂点候補Ｐ１２の方が頂点候補Ｐ１１よりも高い）。ユーザが判断する場合には、重複状況を表示部１４に表示させ、ユーザが表示内容を見て操作受付部１５により一つを選択する又は一つ以外を除去する操作を行えばよい。
【００３９】
必要に応じて確定操作がなされることで、頂点位置が特定される。特定された頂点位置のデータは、直接、及び／又は統計処理などがなされて、表示部１４の表示画面への表示がなされたり、データベース装置２１及び／又は外部機器に送信されたり、光学読取装置２２にセットされた光学ディスクに書き込みが可能とされたりしてよい。
【００４０】
図７は、本実施形態の処理装置１で実行される樹木頂点特定処理の制御部１１による制御手順を示すフローチャートである。この処理は、例えば、ユーザが操作受付部１５から所定の位置範囲を指定して開始命令を行うことで開始される。
【００４１】
樹木頂点特定処理が開始されると、制御部１１は、指定されたエリアに対応するＤＳＭデータ２０１及びＤＥＭデータ２０２を取得する（ステップＳ１０１）。制御部１１は、ＤＳＭデータ２０１の各位置における標高値と、ＤＥＭデータ２０２の各位置における標高値とを取得して、差分値を算出し、水平座標と対応付けてＤＣＨＭデータを生成する（ステップＳ１０２）。
【００４２】
制御部１１は、予め定められた局所範囲に基づいて最大値（極大値）を検出し、頂点候補として抽出する（ステップＳ１０３）。制御部１１は、抽出された頂点候補のうち基準外のデータを除外する（ステップＳ１０４）。基準外のデータは、例えば、基準高さ範囲から外れたものである。例えば、上述のように、送電線が上限基準高さ（例えば、３０ｍ）より高い高さＨで検出され、下草などが下限基準高さ（例えば、５ｍ）未満の高さＨで検出される。なお、ステップＳ１０４の処理がステップＳ１０３の処理の前に行われて、予め基準外のデータが除外されてもよい。
【００４３】
制御部１１は、抽出された頂点候補から未選択のものを順番に選択する（ステップＳ１０５）。制御部１１は、ＤＣＨＭを用いて、選択した頂点候補の位置における高さＨを取得する（ステップＳ１０６）。
【００４４】
制御部１１は、樹木分布情報２０３を参照して、頂点候補の水平座標に対応する樹木種Ｓを特定する（ステップＳ１０７）。制御部１１は、地理情報２０４を参照して、頂点候補の水平座標に対応するエリアＧを特定する（ステップＳ１０８）。制御部１１は、樹木種Ｓ及びエリアＧに対応する係数ａ、ｂを特定し、これら係数ａ、ｂ及び高さＨを用いて領域幅Ｒを算出する（ステップＳ１０９）。
【００４５】
制御部１１は、全ての頂点候補を選択したか否かを判別する（ステップＳ１１０）。選択していない頂点候補があると判別された場合には（ステップＳ１１０で“ＮＯ”）、制御部１１の処理は、ステップＳ１０５に戻る。全ての頂点候補が選択されたと判別された場合には（ステップＳ１１０で“ＹＥＳ”）、制御部１１の処理は、ステップＳ１１１に移行する。
【００４６】
ステップＳ１１１の処理に移行すると、制御部１１は、各頂点候補の位置（水平座標）を中心として、当該頂点候補に対して求められた領域幅Ｒを半径とするエリアを候補領域として各々設定する（ステップＳ１１１；幅設定ステップ、幅設定手段（ステップＳ１０９を含んでよい））。制御部１１は、設定された各候補領域の間で重複部分があるか否かを判別する（ステップＳ１１２）。重複部分がないと判別された場合には（ステップＳ１１２で“ＮＯ”）、制御部１１の処理は、ステップＳ１１６に移行する。
【００４７】
重複部分があると判別された場合には（ステップＳ１１２で“ＹＥＳ”）、制御部１１は、重複する頂点候補を表示部１４により表示させる（ステップＳ１１３）。表示は、処理がなされた全体についてなされてもよいし（制御部１１は、重複部分のみ色を異ならせるなどでユーザが容易に知得可能とさせることができる）、重複部分のみを拡大してなされてもよい。
【００４８】
制御部１１は、重複部分について、いずれか一つが選択される操作が検出されたか否かを判別する（ステップＳ１１４）。検出されていないと判別された場合には（ステップＳ１１４で“ＮＯ”）、制御部１１は、ステップＳ１１４の処理を繰り返す。
【００４９】
いずれか一つが選択される操作が検出されたと判別された場合には（ステップＳ１１４で“ＹＥＳ”）、制御部１１は、選択されていない頂点候補を全て除去する（ステップＳ１１５）。それから、制御部１１の処理は、ステップＳ１１６に移行する。
【００５０】
ステップＳ１１２、Ｓ１１５の処理からステップＳ１１６の処理に移行すると、制御部１１は、残っている頂点候補を樹木の頂点として確定する（ステップＳ１１６）。そして、制御部１１は、樹木頂点特定処理を終了する。
【００５１】
上述のように、手動選択を用いずに機械的に自動で重複エリアから一つの頂点候補を選択する場合には、ステップＳ１１３、Ｓ１１４の処理の代わりに適宜上述のような選択処理がなされて、選択されていない頂点候補が除去されればよい。
【００５２】
以上のように、本実施形態の処理装置１は、制御部１１を備え、制御部１１は、幅設定部として、ＤＳＭモデル及びＤＥＭモデルを用いて抽出された樹木の頂点候補の各々における樹木の特性、すなわち、樹木の高さＨ及び樹木種Ｓに基づいて、当該樹木の領域幅Ｒを設定する。
このように、樹木ごとに領域幅Ｒを定量的に評価して設定するので、容易かつより安定して精度のよい領域幅Ｒを定めることができる。この領域幅Ｒを用いることによって、より安定して精度よく樹木の頂点を特定することが可能になる。
また、樹木の頂点の特定時に、最小樹木幅の最適値を試行錯誤したり、ユーザによる立木密度の読み取りに経験を要したりする必要がない。特に、頂点位置の特定範囲における立木密度の読み取りを予め行って設定を行ってから処理を開始したり、一度頂点候補を定めてからその周囲の立木密度を読み取って再度処理をしなおしたりするような二度手間を必要としないので、より効率よく正確な樹木の頂点位置の特定に寄与することができる。
【００５３】
また、樹木の特性には、頂点候補の位置における樹木の高さＨが含まれる。同一の林分内の同時期に植えられた樹木であっても、地形などの生育条件に応じて高さにはばらつきを生じる場合がある。また、異なる時期に植えられた樹木であれば、更に大きな高さの差が生じ得る。これらに応じて、頂点候補の位置における樹木の高さＨを用いて当該樹木の領域幅Ｒを定めるので、個々の木に対してより適切に領域幅Ｒが定められ、当該樹木の周囲における過剰な特定や必要以上の樹木の除外を避けることができる。
【００５４】
また、制御部１１は、種別特定部として、頂点候補の位置のそれぞれにおける樹木種を特定する。領域幅Ｒの算出に係る樹木の特性には、特定された樹木種が含まれる。
これにより、対象となる樹木の樹冠の形状の違いに適切に対応し、より正確な領域幅Ｒを定めることができる。これにより、より適切に各樹頂点を特定し、過剰な特定や除外を避けることができる。
【００５５】
また、制御部１１は、種別特定部として、樹木分布情報２０３を取得し、当該樹木分布情報２０３における頂点候補の位置を含む領域の樹木種Ｓの情報により当該樹木種を特定する。すなわち、予め定められている樹木分布情報２０３を参照して対応する位置の樹木種を取得するので、容易に樹木種を特定して適切な領域幅Ｒを定めることができる。
【００５６】
また、制御部１１は、特性特定部として、地理情報２０４を取得し、当該地理情報２０４における頂点候補の位置を含む領域のエリアＧの情報により、環境特性情報を特定する。制御部１１は、幅設定部として、環境特性情報に基づいて領域幅Ｒを設定する。すなわち、制御部１１は、樹木の特性だけではなく、当該樹木の生育環境も考慮して領域幅Ｒを定める。これにより、より精度よく各樹木の領域幅Ｒを定めることができる。
【００５７】
また、環境特性情報には、頂点候補の位置の地理的特性が含まれる。すなわち、地理的なエリアに応じた気候条件を反映して、領域幅Ｒをより適切に設定することができる。また、各地方や都道府県自治体などの地理的なエリアの情報は、比較的容易に用いることができるので、少ない手間でより確実に適切な領域幅Ｒに基づき、樹木の頂点を特定することができる。
【００５８】
また、制御部１１は、幅設定部として、最多立木密度を求める式を用いて樹木の特性に基づく値を算出し、当該値に応じて領域幅Ｒを設定する。最多立木密度は、最小の樹木間距離に対応する。樹木間距離は、当該樹木と、この樹木に隣り合う樹木との勢力範囲を示す半径ｒの和となるので、結局この樹木間距離から領域幅Ｒが得られる。したがって、樹木ごとにその高さＨを用いて最多立木密度に当たる値を求めることで、よりきめ細かく各樹木の領域幅Ｒを想定される最低範囲で定めることができる。これにより、樹木の頂点の特定時に、より正確に重複候補の除外を行い、また、不要な排除を避けることができる。
【００５９】
また、ＤＳＭモデル及びＤＥＭモデルは、レーザプロファイラを用いて算出されたものである。樹林帯では、隙間のある樹木を通してＤＳＭモデルとＤＥＭモデルがよい精度で同時に取得されるので、効率よくＤＣＨＭモデルを取得することができ、したがって、適切な頂点候補の抽出及び当該頂点候補の高さＨの特定が可能となる。よって、処理装置１では、より精度のよい領域幅Ｒの算出が可能となり、これに応じて頂点の特定をより正確に行うことが可能になる。
【００６０】
また、制御部１１は、候補抽出部として、ＤＳＭモデル及びＤＥＭモデルを用いて樹木の頂点候補を抽出する。これにより、処理装置１では、頂点候補を抽出して、抽出した頂点候補に対してより適切に頂点であるか否かを判別するための適切な領域幅Ｒを定めることができる。
【００６１】
また、樹木頂点特定処理の実施形態としての処理装置１は、制御部１１が、上述の樹木幅設定装置としての各機能動作と、樹頂点特定部として、設定された領域幅Ｒで定まる各樹木の領域が重複しないように頂点候補から頂点を特定する処理とを行う。すなわち、上述のように領域幅Ｒが定められた各頂点候補における領域（頂点候補の位置を頂点とした直径が領域幅Ｒの円形領域）が重複している場合には、いずれかの頂点候補が誤検出であると判断されて、重複領域について１つのみが選択される。この重複領域を定める領域が上述の領域幅Ｒの設定によりより適切になされているので、処理装置１では、重複の判定及び頂点の特定をより正確に行うことができる。
【００６２】
また、上記のような処理により、上述の処理装置１と同様に容易かつより安定して精度のよい領域幅Ｒを定めることができる。この領域幅Ｒを用いることによって、より安定して精度よく樹木の頂点を特定することが可能になる。
【００６３】
また、本実施形態のプログラム１２１は、コンピュータ（処理装置１）を、ＤＳＭモデル及びＤＥＭモデルを用いて抽出された樹木の頂点候補の各々における樹木の特性に基づいて、当該樹木の領域幅Ｒを設定する幅設定手段として機能させる。
このような処理を行うプログラム１２１をインストールしてコンピュータに処理を行わせることで、ソフトウェア的に容易かつ適切に、より精度のよい領域幅Ｒを安定して定めることができる。これにより、より安定して精度良く樹木の頂点を特定することが可能になる。
【００６４】
なお、本発明は、上記実施の形態に限られるものではなく、様々な変更が可能である。
例えば、上記実施の形態では、高さＨ、樹木種Ｓ及びエリアＧを用いて最多立木密度Ｎｒを求める式に基づいて領域幅Ｒを算出し、また、独自データや第三者データがもちいられてもよい旨示したが、例えば、特定の樹木種Ｓに係る全地域の算出パラメータの代表値（平均値）などを取得して、エリアＧの情報なしに算出が可能に変形してもよい。データの読み込みや処理の負荷を若干低減し、処理を簡略化することができる。同様に、算出における高さＨ及び樹木種Ｓのうちいずれかを省いても（固定値などとしても）、従来に比して若干の精度の向上を図ることができる。
【００６５】
また、本実施形態では、ＤＳＭデータ２０１、ＤＥＭデータ２０２、樹木分布情報２０３、地理情報２０４及び林分密度テーブル２０５の各データは、外部のデータベース装置２１又は可搬型記憶媒体から取得されるものとして説明したが、処理装置１がこれらのうち一部又は全部を予め内部に保持していてもよい。
【００６６】
また、上記実施の形態では、頂点候補の抽出にＤＣＨＭの局所的な最大値の検出を利用したが、これに限られるものではない。その他の機械的に抽出可能な方法を用いて又は併用してもよいし、あるいは、ＤＣＨＭデータからノイズが除去されたデータをそのまま用いてもよい。
【００６７】
また、上記実施の形態では、環境特性情報としての地理情報２０４は、環境条件に応じたエリア分けの情報として、都道府県や、地方の太平洋側／日本海側などの区分を例に挙げて説明したが、都道府県内でも積雪状況、例えば、北側斜面や南側斜面などによってより細かく区分を行ってもよい。あるいは、上記の区分に対して、積雪状況に基づいて補正がなされるように、積雪状況の情報やその数値補正用のデータを保持していてもよい。
【００６８】
また、上記実施の形態では、ＤＳＭデータ２０１及びＤＥＭデータ２０２がレーザプロファイラによって取得されたものであるとして説明したが、これに限られない。例えば、異なる位置で撮影された航空写真を用いたステレオマッチングなどにより、高度分布が取得されてもよい。
【００６９】
また、上記実施の形態では、候補の抽出、領域幅の設定及び頂点の特定に係る各処理を単一の処理プログラム（樹木特定処理プログラム）により行われることとして説明したが、これらの各処理、及びこれら各処理の一部は、異なるプログラムに分割されて作成され、別個に呼び出されて実行されてもよい。
【００７０】
また、樹木特定処理プログラムに係るソフトウェア的な処理の一部が専用の論理回路などのハードウェアによってなされてもよい。
【００７１】
また、これらの各処理は、異なる処理装置（コンピュータ）に分散されてなされてもよい。
【００７２】
その他、上記実施の形態で示した装置構成、処理内容や処理手順などの具体的な細部は、本発明の趣旨を逸脱しない範囲において適宜変更可能である。
【００７３】
また、上記実施の形態では、プログラム１２１が記憶部１２の不揮発性メモリ又はＨＤＤなどに記憶されるものとして説明したが、これに限られない。プログラム１２１は、ＣＤＲＯＭ、ＤＶＤ、又はＢｌｕ−ｒａｙなどの光学ディスクといった可搬型記憶媒体に記憶されてよい。また、本発明に係るプログラムのデータを通信回線を介して提供する媒体として、キャリアウェーブ（搬送波）も本発明に適用される。
【００７４】
［実施例］
北近畿・中国地方エリアにおけるスギ及びヒノキの樹林について、スギ及び／又はヒノキを樹木種とし（一部のエリアではスギとヒノキが共存しており、スギが優勢である）、立木密度の異なる半径１１．２８ｍの円形エリアを設定して目視により実際に計測した樹木の本数と、機械的に特定した樹木の本数とを比較した。
【００７５】
これらの各円形エリアに対し、本実施形態を適用した場合として、まず半径６０ｃｍのエリア設定を用いて局所的最大値を求めて樹木の頂点候補を抽出した。次いで、上記エリア及び主な樹木種を適用して得られた最多立木密度の係数ａ、ｂと、頂点候補における高さＨを用いて、最小の樹木間距離Ｄを算出した。そして、頂点候補を中心として、得られた樹木間距離Ｄを領域幅Ｒとして、当該領域幅Ｒを直径とする円領域を候補領域として定め、頂点候補間で候補領域の重複があった場合には、重複がなくなるまで頂点候補を除去した。
【００７６】
一方、比較例としては、局所的最大値を求めるエリアの半径ｒを６０ｃｍ〜１８０ｃｍの固定値として定めて抽出された頂点候補を、そのまま樹頂点として特定した。
【００７７】
図８は、樹木の本数を特定した結果の例を示す図表である。
設定された２０箇所の円形エリアにおいて、目視によって得られた実際の本数は、１７〜６７本の間に分布している。すなわち、樹木間距離には、大きな開きがある。これに対し、半径ｒ＝６０ｃｍの場合には、特定数が目視数より多く、過剰に特定されている場合が多い。このときの平均平方二乗誤差（Root Mean Square Error；ＲＭＳＥ）は、１２．２である。
【００７８】
半径ｒが６０ｃｍよりも大きくなると、特定本数は減少していく。半径Ｄが１００ｃｍでは、特定数が目視数よりも少ない場合が多く、既に過小検出となっている。ＲＭＳＥは、１０．２と、半径ｒが６０ｃｍの場合より若干小さい。半径ｒが１４０ｃｍ、１８０ｃｍでは、更に検出数が減少し、ＲＭＳＥがそれぞれ１７．８、２２．１となっている。
【００７９】
これに対し、本実施例では、目視数に対して過剰検出又は過小検出のいずれにも偏らず、また、全体としてずれが小さく抑えられている。ＲＭＳＥは、４．６であり、比較例における半径ｒ＝６０ｃｍ、１００ｃｍの場合よりも十分に小さい。事例ごとに見ると、比較例のいずれかにおける特定数が、実施例による特定数よりも目視数に近い場合もあるが、本実施例では、全体としてのずれ量が安定して小さく抑えられていることが示されている。
【００８０】
さらに、樹木種がスギのエリア（８エリア）とヒノキのエリア（１０エリア）とで区別した場合、目視での本数がスギで１７〜６７本、ヒノキで２６〜４９本に対し、本実施例では、スギ及びヒノキに係るＲＭＳＥは、それぞれ、５．０、４．７と安定して小さい値が維持されている。これに対し、比較例のｒ＝６０ｃｍでは、スギに係るＲＭＳＥの値（４．１）に比してヒノキに係るＲＭＳＥの値（１６．７）が大きい。反対に、ｒ＝１００ｃｍでは、スギに係るＲＭＳＥの値（１２．８）がヒノキに係るＲＭＳＥの値（８．５）より大きくなっている。すなわち、本実施例では、樹木種に応じた違いに対して柔軟に対応して、より安定して精度よく特定されるという結果が得られた。
【符号の説明】
【００８１】
１     処理装置
１１   制御部
１２   記憶部
１２１ プログラム
１３   入出力インターフェイス
１３１ 接続端子
１３２ 通信部
１４   表示部
１５   操作受付部
２１   データベース装置
２２   光学読取装置
２０１ ＤＳＭデータ
２０２ ＤＥＭデータ
２０３ 樹木分布情報
２０４ 地理情報
２０５ 林分密度テーブル
Ｄ     樹木間距離
Ｇ     エリア
Ｒ     領域幅
Ｓ     樹木種
【図１】
2020098442.tif 000003
【図２】
2020098442.tif 000004
【図３】
2020098442.tif 000005
【図４】
2020098442.tif 000006
【図５】
2020098442.tif 000007
【図６】
2020098442.tif 000008
【図７】
2020098442.tif 000009
【図８】
2020098442.tif 000010