JP7248419B2 - 樹木幅設定装置及びプログラム - Google Patents
樹木幅設定装置及びプログラム Download PDFInfo
- Publication number
- JP7248419B2 JP7248419B2 JP2018235989A JP2018235989A JP7248419B2 JP 7248419 B2 JP7248419 B2 JP 7248419B2 JP 2018235989 A JP2018235989 A JP 2018235989A JP 2018235989 A JP2018235989 A JP 2018235989A JP 7248419 B2 JP7248419 B2 JP 7248419B2
- Authority
- JP
- Japan
- Prior art keywords
- tree
- vertex
- area
- height
- width
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
樹木の特性と、当該樹木の領域幅に係る情報との対応関係を記憶する記憶部と、
樹木表面を含む高さ分布を示す数値表層モデル及び地表面の高さ分布を示す数値標高モデルとの差分により得た樹木の高さの極大点を樹木の頂点候補として抽出し、抽出された前記樹木の頂点候補の各々における樹木の特性に応じた領域幅を前記記憶部から取得することで、当該樹木のそれぞれの領域幅を設定する幅設定部と、
を備えることを特徴とする樹木幅設定装置である。
図1は、樹木幅設定装置の実施形態である処理装置1の機能構成を示すブロック図である。
図2(a)に示すように、DSMの高度は、地上の構造物や樹林などを含めたものであり、樹林帯では、樹冠の上縁に係る高度分布を示す。これらの値は、例えば、レーザプロファイラによって計測、取得されたものである。この高度には、後述のように、送電線Wなども含まれ得る。また、葉や枝の隙間を抜けて地面や下草の高度が取得されたものが混じっている場合がある。後述のように、樹木の高さとして適切でない高さのデータを除外するために、上限基準高さHmaxや下限基準高さHminが定められている。水平座標は、緯度経度の値であってもよいし、所定の原点位置からの相対距離を示すものであってもよい。
図3(a)の樹木分布情報203において、単位位置範囲Bは、例えば、矩形(正方形)形状であり、大きさは、上記のDSMデータ201などにおける単位面積と同一であってもよいし、当該単位面積の所定倍などであってもよい。林相図では、樹木種の分布範囲ごとにエリアS1、S2、S3などが定められている。空白部分は、森林以外の部分である。これらのうち、特に単一樹種のエリアについては、具体的な樹種が明示されている。なお、メッシュ単位では樹木種がそれぞれID値などで示され、別途ID値などと樹木種の名称とが対応付けられていてもよい。また、最終的に樹木種に応じた領域幅の算出が可能であればよいので、処理装置1内で理解可能な対応付けがなされていれば(ID値のみで全ての処理が可能であれば)、樹木種の名称が明示されている必要はない。樹木種としては、具体的には、ヒノキ、スギ、アカマツなどが挙げられる。複数種類の雑木林などについては、より広く広葉樹林又は針葉樹林などで定められてもよい。例えば、DSM、DEM、DCHMなどにおける水平位置の緯度経度情報に応じて指定された地点P1は、エリアS2に属し、当該エリアS2に対応する樹木種の樹木であると特定される。
ここでは、まず、数値表層モデル(DSM)と、地表面の高さ分布を示す数値標高モデル(DEM)との差分をとることで樹冠高データ(DCHM)を取得する(PR1)。当該DCHMデータから高さの極大点を樹木の頂点候補として抽出する(PR2;候補抽出部)。抽出方法としては、例えば、従来の局所的な最大値(極大値)を取得する方法などが用いられてよい。
log(Nr)=a-b×log(Ha) … (数式1)
係数a、bは、樹木種S及びエリアGにより定まる数値パラメータである。すなわち、林分密度テーブル205には、樹木種S及びエリアGの組み合わせに応じた係数a、bの値がそれぞれ設定、記憶されている。
図6(a)に示すように、緯度経度で示される頂点候補に対し、それぞれ、その高さH及び樹木種S(樹木の特性)、並びにエリアG(地理的特性)が定められている。これらに基づいて得られた樹木間距離Dに応じて、頂点候補が頂点である場合の領域幅Rが設定されて、頂点候補に係る領域(候補領域)として対応付けられる(PR6;幅設定部(PR5を含んでよい))。
図6(b)に示すように、各頂点候補(P11~P14)を中心として領域幅(R11~R14)をそれぞれ直径とする範囲(候補領域)を設定した場合に、当該候補領域が他の頂点候補について定められた候補領域と重複するか否かを検出する。ここでは、頂点候補P11の候補領域と頂点候補P12の候補領域とが重複している。重複する場合には、いずれかの頂点候補が誤抽出されたものであると判断される。したがって、重複範囲がなくなるように、いずれかの頂点候補が除去されて、残った頂点候補が頂点として特定される。
このように、樹木ごとに領域幅Rを定量的に評価して設定するので、容易かつより安定して精度のよい領域幅Rを定めることができる。この領域幅Rを用いることによって、より安定して精度よく樹木の頂点を特定することが可能になる。
また、樹木の頂点の特定時に、最小樹木幅の最適値を試行錯誤したり、ユーザによる立木密度の読み取りに経験を要したりする必要がない。特に、頂点位置の特定範囲における立木密度の読み取りを予め行って設定を行ってから処理を開始したり、一度頂点候補を定めてからその周囲の立木密度を読み取って再度処理をしなおしたりするような二度手間を必要としないので、より効率よく正確な樹木の頂点位置の特定に寄与することができる。
これにより、対象となる樹木の樹冠の形状の違いに適切に対応し、より正確な領域幅Rを定めることができる。これにより、より適切に各樹頂点を特定し、過剰な特定や除外を避けることができる。
このような処理を行うプログラム121をインストールしてコンピュータに処理を行わせることで、ソフトウェア的に容易かつ適切に、より精度のよい領域幅Rを安定して定めることができる。これにより、より安定して精度良く樹木の頂点を特定することが可能になる。
例えば、上記実施の形態では、高さH、樹木種S及びエリアGを用いて最多立木密度Nrを求める式に基づいて領域幅Rを算出し、また、独自データや第三者データがもちいられてもよい旨示したが、例えば、特定の樹木種Sに係る全地域の算出パラメータの代表値(平均値)などを取得して、エリアGの情報なしに算出が可能に変形してもよい。データの読み込みや処理の負荷を若干低減し、処理を簡略化することができる。同様に、算出における高さH及び樹木種Sのうちいずれかを省いても(固定値などとしても)、従来に比して若干の精度の向上を図ることができる。
北近畿・中国地方エリアにおけるスギ及びヒノキの樹林について、スギ及び/又はヒノキを樹木種とし(一部のエリアではスギとヒノキが共存しており、スギが優勢である)、立木密度の異なる半径11.28mの円形エリアを設定して目視により実際に計測した樹木の本数と、機械的に特定した樹木の本数とを比較した。
設定された20箇所の円形エリアにおいて、目視によって得られた実際の本数は、17~67本の間に分布している。すなわち、樹木間距離には、大きな開きがある。これに対し、半径r=60cmの場合には、特定数が目視数より多く、過剰に特定されている場合が多い。このときの平均平方二乗誤差(Root Mean Square Error;RMSE)は、12.2である。
11 制御部
12 記憶部
121 プログラム
13 入出力インターフェイス
131 接続端子
132 通信部
14 表示部
15 操作受付部
21 データベース装置
22 光学読取装置
201 DSMデータ
202 DEMデータ
203 樹木分布情報
204 地理情報
205 林分密度テーブル
D 樹木間距離
G エリア
R 領域幅
S 樹木種
Claims (5)
- 樹木の特性と、当該樹木の領域幅に係る情報との対応関係を記憶する記憶部と、
樹木表面を含む高さ分布を示す数値表層モデル及び地表面の高さ分布を示す数値標高モデルとの差分により得た樹木の高さの極大点を樹木の頂点候補として抽出し、抽出された前記樹木の頂点候補の各々における樹木の特性に応じた領域幅を前記記憶部から取得することで、当該樹木のそれぞれの領域幅を設定する幅設定部と、
を備えることを特徴とする樹木幅設定装置。 - 前記特性には、前記頂点候補の位置における前記樹木の高さが含まれることを特徴とする請求項1記載の樹木幅設定装置。
- 前記頂点候補の位置のそれぞれにおける樹木種を特定する種別特定部を備え、
前記種別特定部は、樹木種の分布情報を取得し、当該分布情報における前記頂点候補の位置を含む領域の樹木種の情報により当該樹木種を特定し、
前記特性には、特定された前記樹木種が含まれる
ことを特徴とする請求項1又は2記載の樹木幅設定装置。 - 前記頂点候補の位置の環境特性情報を特定する特性特定部を備え、
前記環境特性情報には、前記頂点候補の位置の地理的特性が含まれ、
前記幅設定部は、前記環境特性情報に基づいて前記領域幅を設定する
ことを特徴とする請求項1~3のいずれか一項に記載の樹木幅設定装置。 - コンピュータを、
樹木の特性と、当該樹木の領域幅に係る情報との対応関係を記憶する記憶手段、
樹木表面を含む高さ分布を示す数値表層モデル及び地表面の高さ分布を示す数値標高モデルとの差分により得た樹木の高さの極大点を樹木の頂点候補として抽出し、抽出された前記樹木の頂点候補の各々における樹木の特性に応じた領域幅を前記記憶手段の記憶内容から取得することで、当該樹木のそれぞれの領域幅を設定する幅設定手段
として機能させることを特徴とするプログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018235989A JP7248419B2 (ja) | 2018-12-18 | 2018-12-18 | 樹木幅設定装置及びプログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018235989A JP7248419B2 (ja) | 2018-12-18 | 2018-12-18 | 樹木幅設定装置及びプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2020098442A JP2020098442A (ja) | 2020-06-25 |
JP7248419B2 true JP7248419B2 (ja) | 2023-03-29 |
Family
ID=71106552
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2018235989A Active JP7248419B2 (ja) | 2018-12-18 | 2018-12-18 | 樹木幅設定装置及びプログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP7248419B2 (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7323689B1 (ja) | 2022-08-25 | 2023-08-08 | 朝日航洋株式会社 | 表示装置、表示方法、及び、表示プログラム |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008111724A (ja) | 2006-10-30 | 2008-05-15 | Tokyo Electric Power Co Inc:The | 樹木頂点認識方法及び樹木頂点認識のプログラム並びに樹木頂点認識装置 |
US20080319673A1 (en) | 2007-06-22 | 2008-12-25 | Weyerhaeuser Co. | Identifying vegetation attributes from LiDAR data |
JP2009022278A (ja) | 2007-06-19 | 2009-02-05 | Tokyo Electric Power Co Inc:The | 樹木頂点認識方法及び樹木頂点認識装置並びに樹木頂点認識のプログラム |
JP2010086276A (ja) | 2008-09-30 | 2010-04-15 | Shinshu Univ | 樹種分類方法及び樹種分類システム、森林現況情報の作成方法及び森林現況情報の作成システム、間伐対象区域の選定方法及び間伐対象区域の選定システム |
US20100250482A1 (en) | 2009-03-24 | 2010-09-30 | Weyerhaeuser Nr Company | System and method for identifying trees using lidar tree models |
JP2012098247A (ja) | 2010-11-05 | 2012-05-24 | Pasuko:Kk | 樹木位置検出装置、樹木位置検出方法、及びプログラム |
JP2012103155A (ja) | 2010-11-11 | 2012-05-31 | Pasuko:Kk | レーザ密度分布推定装置、レーザ密度分布推定方法、及びプログラム |
JP2014100099A (ja) | 2012-11-20 | 2014-06-05 | Kodaira Associates Kk | 森林情報管理装置 |
JP2018084472A (ja) | 2016-11-22 | 2018-05-31 | 国立大学法人信州大学 | 森林資源情報算定方法及び森林資源情報算定装置 |
-
2018
- 2018-12-18 JP JP2018235989A patent/JP7248419B2/ja active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008111724A (ja) | 2006-10-30 | 2008-05-15 | Tokyo Electric Power Co Inc:The | 樹木頂点認識方法及び樹木頂点認識のプログラム並びに樹木頂点認識装置 |
JP2009022278A (ja) | 2007-06-19 | 2009-02-05 | Tokyo Electric Power Co Inc:The | 樹木頂点認識方法及び樹木頂点認識装置並びに樹木頂点認識のプログラム |
US20080319673A1 (en) | 2007-06-22 | 2008-12-25 | Weyerhaeuser Co. | Identifying vegetation attributes from LiDAR data |
JP2010086276A (ja) | 2008-09-30 | 2010-04-15 | Shinshu Univ | 樹種分類方法及び樹種分類システム、森林現況情報の作成方法及び森林現況情報の作成システム、間伐対象区域の選定方法及び間伐対象区域の選定システム |
US20100250482A1 (en) | 2009-03-24 | 2010-09-30 | Weyerhaeuser Nr Company | System and method for identifying trees using lidar tree models |
JP2012098247A (ja) | 2010-11-05 | 2012-05-24 | Pasuko:Kk | 樹木位置検出装置、樹木位置検出方法、及びプログラム |
JP2012103155A (ja) | 2010-11-11 | 2012-05-31 | Pasuko:Kk | レーザ密度分布推定装置、レーザ密度分布推定方法、及びプログラム |
JP2014100099A (ja) | 2012-11-20 | 2014-06-05 | Kodaira Associates Kk | 森林情報管理装置 |
JP2018084472A (ja) | 2016-11-22 | 2018-05-31 | 国立大学法人信州大学 | 森林資源情報算定方法及び森林資源情報算定装置 |
Also Published As
Publication number | Publication date |
---|---|
JP2020098442A (ja) | 2020-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Leitold et al. | Airborne lidar-based estimates of tropical forest structure in complex terrain: opportunities and trade-offs for REDD+ | |
JP5360989B2 (ja) | 地理情報生成システム及び地理情報生成方法 | |
US8897483B2 (en) | System and method for inventorying vegetal substance | |
US20090185741A1 (en) | Apparatus and method for automatic airborne LiDAR data processing and mapping using data obtained thereby | |
US8582808B2 (en) | Methods for identifying rooftops using elevation data sets | |
CN111340012B (zh) | 一种地质灾害解译方法、装置、终端设备 | |
US9098745B2 (en) | Sampling position-fixing system | |
US20160093056A1 (en) | Multi-spectral image labeling with radiometric attribute vectors of image space representation components | |
Xie et al. | Accuracy assessment and error analysis for diameter at breast height measurement of trees obtained using a novel backpack LiDAR system | |
Yadav et al. | Identification of trees and their trunks from mobile laser scanning data of roadway scenes | |
Arkin et al. | Integrated fire severity–land cover mapping using very-high-spatial-resolution aerial imagery and point clouds | |
CN111222539A (zh) | 基于多源多时相遥感影像优化和扩充监督分类样本的方法 | |
JP7248419B2 (ja) | 樹木幅設定装置及びプログラム | |
Potůčková et al. | Comparison of quality measures for building outline extraction | |
CN115457408A (zh) | 一种土地监测方法、装置、电子设备及介质 | |
JP4946072B2 (ja) | 樹頂探索方法、装置及びプログラム | |
Alvites et al. | Unsupervised algorithms to detect single trees in a mixed-species and multilayered Mediterranean forest using LiDAR data | |
Norman et al. | Fusion of multispectral imagery and LiDAR data for roofing materials and roofing surface conditions assessment | |
CN114519712A (zh) | 一种点云数据处理方法、装置、终端设备及存储介质 | |
Xiao et al. | Individual tree detection from multi-view satellite images | |
Yin et al. | Object-based larch tree-crown delineation using high-resolution satellite imagery | |
Sarıtaş et al. | Enhancing Ground Point Extraction in Airborne LiDAR Point Cloud Data Using the CSF Filter Algorithm | |
Wężyk et al. | Use of airborne laser scanning data for a revision and update of a digital forest map and its descriptive database: a case study from the Tatra National Park | |
US20220277418A1 (en) | Method and electronic portal system for building façade inspection | |
CN113111793A (zh) | 树木识别方法与装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20211021 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20220728 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20220906 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20221101 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20230307 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20230316 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7248419 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |