JP4004332B2 - 森林情報処理システム - Google Patents
森林情報処理システム Download PDFInfo
- Publication number
- JP4004332B2 JP4004332B2 JP2002147774A JP2002147774A JP4004332B2 JP 4004332 B2 JP4004332 B2 JP 4004332B2 JP 2002147774 A JP2002147774 A JP 2002147774A JP 2002147774 A JP2002147774 A JP 2002147774A JP 4004332 B2 JP4004332 B2 JP 4004332B2
- Authority
- JP
- Japan
- Prior art keywords
- tree
- forest
- information
- height
- crown
- 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.)
- Expired - Lifetime
Links
- 230000010365 information processing Effects 0.000 title claims description 17
- 238000005259 measurement Methods 0.000 claims description 25
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000000605 extraction Methods 0.000 claims description 5
- 239000000463 material Substances 0.000 claims description 4
- 241000894007 species Species 0.000 description 28
- 238000000034 method Methods 0.000 description 14
- 238000012545 processing Methods 0.000 description 11
- 238000001228 spectrum Methods 0.000 description 10
- 241000218691 Cupressaceae Species 0.000 description 5
- 241000219492 Quercus Species 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000010792 warming Methods 0.000 description 4
- 239000002023 wood Substances 0.000 description 4
- 241000218645 Cedrus Species 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 239000003550 marker Substances 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 241000218692 Cryptomeria Species 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 241000039951 Lithocarpus glaber Species 0.000 description 1
- 240000004050 Pentaglottis sempervirens Species 0.000 description 1
- 235000004522 Pentaglottis sempervirens Nutrition 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000009532 heart rate measurement Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Image Processing (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Image Analysis (AREA)
Description
【発明の属する技術分野】
本発明は、森林地域に関して上空から取得した樹木情報に基づいて森林を評価する森林情報処理システムに関し、特に、森林地域を航空機等で上空から取得した森林画像情報と、上空から距離計測手段により計測した該地域の樹高情報とにより、該地域の森林を評価するための森林材積を求める森林情報処理システムに関する。
【0002】
【従来の技術】
従来における森林の評価は、広大な森林地域の比較的狭い地域について予め調査区域として定めて、森林地域に実際に赴き、異なる樹木の分布状態、樹木の埴生の疎密状態、樹木の幹の直径及び樹高などを各種計測器又は目視により測定し、この調査区域内での実地調査結果に基づいて、当該調査区域の森林を評価していた。狭い森林地域を広大な森林地域を代表させ、このような部分的な調査結果に基づいた推測により森林地域全体を評価している。
【0003】
従来の森林地域を評価するための手順は、
(a)自治体が有する縮尺5000分の1の森林計画図、或いは、航空写真を参考にして森林地域内の一部分を調査区域として定める。
(b)適当な交通手段を用いて調査区域の近くまで行って、徒歩により調査区域に入る。
(C)GPS(Global Positioning System)が利用可能であれば、GPSを用いて調査地点の位置を求める。GPSが利用不可能であれば周辺の地形状況などから調査地点の位置を特定する。
(d)調査結果を地図及び調査ノートなどに記入する。
というようなものであった。
【0004】
さらに、近年では、世界各国において、地球温暖化対策が検討されている。その温暖化防止策の例として、世界各地における植林活動が、より経済的に炭酸ガス削減目標を達成するものとして行われている。どのような対策が温暖化防止策として採られるにしても、森林の炭素ストック変化を正確に計測することが必要である。そこで、森林による炭素ストック変化を算定するためには、材積・樹高・樹冠疎密度・樹冠径区分などの情報をモニタリングすることは重要である。
【0005】
【発明が解決しようとする課題】
このモニタリングの中で採用されている従来の材積計測方法においては、森林地域を上空から撮影した航空写真を目視で判読するというやり方で行われ、樹種・樹高・立木密度情報を推定するものであった。さらに、それらの情報の精度を上げるために、上述した現地調査による森林地域評価手順に従った調査員による現地調査情報を加えて推定していた。
【0006】
特に、航空写真の判読精度は、判読する者の技術レベルに大きく依存するものであり、その判読結果は、科学的な、或いは、客観的な情報とは言いがたい場合もあった。また、現地調査する場合には、調査地が危険な急傾斜であったり、崩落しやすい地形であったりと、種々の危険性があり、森林地域を現地調査することは、厳しいものであった。
【0007】
上述のように、調査区域は、車両などが使用できない山地となることが多いため、調査区域内での移動は徒歩に頼らざるを得ず、現地調査には多大な費用と労力がかかるという問題があった。更に、調査区域内の樹種及び樹齢が、他の地域の樹種及び樹齢と同一と仮定しても、樹木の立地条件の違いによって樹木の生育状態が異なるため、限定された調査区域での調査結果に基づいて森林地域全体を正確に評価することは極めて困難であった。
【0008】
更に、樹木が繁茂している地域では、GPSの人工衛星からの位置情報データを受信できない場合が多く、このため、調査地点の位置を正確に求めるのは困難な場合が多かった。
【0009】
従って、本発明は、所定の森林地域を評価する際に重要な要素となっている樹木材積を科学的かつ客観的に計測するために、上空から撮影された森林画像情報と、上空から取得した樹木高さ情報とから当該森林地域の樹木材積を自動計測できる森林情報処理システムを提供することを目的とし、この森林情報処理システムによって、広い地域の森林について樹木材積の調査が必要である場合、各地の森林特性や判読者の技量といった計測誤差の要因を回避でき、現地調査を省略できるようにすることを目的とする。
【0010】
【課題を解決するための手段】
そこで、上述の問題点を解決するため、本発明では、飛行体搭載の画像センサーと距離センサーとで取得される情報に基づいた森林情報取得システムにおいて、森林地域を上空から取得した森林画像情報にウォーターシェッド・アルゴリズムを適用して樹冠形状情報を抽出する樹冠形状抽出手段と、前記森林地域の上空から計測した樹高情報と前記森林画像情報を重ね合せ、前記樹冠形状の範囲内に含まれる前記樹高情報のうち最大値を抽出して当該樹冠の樹高とする樹高計測手段と、前記森林画像内で生成された前記樹冠形状情報に基づいて、前記森林地域内に形成された前記樹冠形状の重心点を求めて該重心点を計数し、該計数値を前記森林地域内の立木本数として取得する立木本数計測手段と、前記森林画像情報に含まれる色情報を前記樹冠形状に基づいて抽出し、色統計情報を参照して当該樹冠の樹種を取得する樹種分類手段と、前記樹高と前記立木本数とに基づいて、前記樹種に応じて樹木の材積を計算する材積計算手段とを備えた。
【0013】
さらに、前記樹高と前記立木本数とに基づいて、前記樹種に応じて樹木の材積を計算する計算手段が備えられ、該計算手段は、前記樹木毎の材積を、材積推定式
V=(B1・HB2+B3・HB4/N)-1
〔但し、Vは材積(m3/ha)、Hは樹高(m)、Nは1haあたりの立木本数であり、B1乃至B4は、樹種や森林地域に従って決定される係数である〕に従って計算することとした。
【0014】
そして、前記森林画像情報は、当該取得位置情報と共に取得され、前記樹高情報は、当該計測位置情報と共に計測され、前記森林画像情報と前記樹高情報との重ね合せは、該取得位置情報と該計測位置情報とに基づいて行われ、或いは、前記森林情報処理システムに使用される前記森林画像情報と前記樹高情報は、同一飛行体に搭載された画像取得手段と距離計測手段により取得されることとした。
【0016】
【発明の実施の形態】
次に、本発明による森林情報処理システムの実施形態について、図を参照しながら説明する。
【0017】
図1(a)に、これから樹木材積を計測しようとしている森林の一部について、その概略を鳥瞰図で示した。同図には、種類の異なる或いは大きさの異なる樹木T1乃至T6が地表面に生えている様子を示している。図中では、代表例を示したものであり、実際には、地表面全域を樹木が埋め尽くしている場合もあれば、単独樹が分散している場合もある。ここで、「樹冠」とは、森林を構成する各樹木の上部において枝と葉で構成された丸みを持った部分を指し、「樹冠形状」とは、森林の各樹木を真上(或いは、略真上)から見た樹冠の形状である。
【0018】
図1(a)に示される各樹木が含まれる森林地域の上方から、航空機或いは人工衛星等に搭載されたカメラによって撮影された画像を図1(b)に示した。図1においては、X軸とY軸とで、水平面を表現しており、Z軸が樹木の高さ方向を表している。ここでは、上方から見た樹木T1乃至T6の樹冠形状を太線で示している。
【0019】
本実施形態では、森林情報処理システムにおける樹木の材積を計測するための原理を説明する。この材積を算出する手法として、人工林の生育に伴う森林密度管理の設計や検討に用いられている地域の材積推定方法に着目した。次式に示す材積推定式は、収量密度効果線の逆算式である。
【0020】
V=(B1・HB2+B3・HB4/N)-1
ここで、Vは、材積(m3/ha)を、Hは、樹高(m)を、Nは、1haあたりの立木本数をそれぞれ表している。また、式中のB1乃至B4は、係数を示しており、樹種や地域によって異なった値をとる。これらの係数は、現地調査結果などを踏まえて、樹種毎、地域毎に決定されている。
【0021】
上記の式によって、各樹木の材積Vを計算するには、樹高Hと、立木本数N、更には、樹種が取得できれば良いことが分かる。
【0022】
商用高解像度衛星IKONOSは、これまでの人工衛星よりも高精度に森林情報を計測することが可能であり、樹種などの情報を調査することができる。この様な高解像度衛星で取得できる画像の他に、航空機搭載型画像センサーにより観測される高解像度画像によっても、より詳細な森林情報を計測できる。
【0023】
先ず、立木本数の計測から説明する。高解像度衛星又は航空機によって取得された森林画像情報から立木本数を計測するために、ウォーターシェッド(Watershed)アルゴリズムという領域分割手法により樹冠形状を抽出する。
【0024】
ここで、樹冠が画像に表示された場合に、一般には、樹木の高い部分ほど、輝度が明るくなり、樹木の周囲部分になるほど、暗くなる傾向がある。そこで、この傾向を利用して、この樹冠における輝度の変化を計測し、樹冠形状を抽出する。この計測に、ウォーターシェッド・アルゴリズムが利用される。先ず、一つの樹冠形状の中において、画像上で輝度の変化率が極小の点にマーカとよばれる領域の核を設定する。そして、輝度の変化方向に領域を近接画素へと拡大していき、隣の樹冠の領域と接するか、輝度の変化率があるしきい値以下になったら領域の拡大を終了し、領域を閉じる。
【0025】
この様にして、樹冠形状が抽出される。さらに、抽出した樹冠形状から樹冠形状ポリゴンを作成し、その重心点を求め、この点を樹木の幹の位置と仮定する。そうすると、その重心点が当該森林画像の地域範囲内に生えている各樹木に該当することになるので、これらをカウントすることにより、その値が立木本数Nとして計測される。
【0026】
図3(a)には、図1(b)に示した例による樹冠形状に基づいて生成された樹冠形状ポリゴンを示した。この樹冠形状ポリゴンの生成手順を説明する。
【0027】
先ず、飛行体に搭載された画像センサーで取得したカラーのデジタル画像を多値の白黒デジタル画像(モノクロデジタル画像)に変換し、ガウスフィルタ等により画像の平滑化を行ない後続の画像処理に備える。この平滑化処理の目的は、画像中に細かい画素値の変動が存在すると、樹冠の形状を求める際、単一の樹冠形状が複数個に分割されるという問題を除去することにある。
【0028】
次に、樹冠形状を求めやすくするため、例えば、ハイパスフィルタ処理を行って、隣接する樹冠の境界線を強調し、樹冠形状T1乃至T6を求める。これらの樹冠形状生成原理は、上空から樹冠を撮影すると、樹冠の盛り上がりに応じて輝度が高くなることに基づいている。一般に、樹冠の頂点部分では、輝度が最も高く、且つこの頂点部分では、輝度変化が最も小さくなっており、樹冠の周辺部に向かうにしたがって輝度が低下し、隣接する樹冠の境界部分で最も輝度が低くなる。本実施形態では、このような樹冠の丸みに伴う輝度の変化を利用し、ウォーターシェッド・アルゴリズムにより樹冠形状を求めている。
【0029】
ウォーターシェッド・アルゴリズムに従って、画像中で輝度の変化が最も小さい部分を有する領域の核にマーカM1乃至M6を設定する。即ち、各マーカは、夫々の樹冠の頂点或いは頂点付近に対応して設定されることになる。図1(b)では、そのマーカについて「×」印で示した。
【0030】
次に、輝度の変化方向に領域を近接画素へと拡大していき、隣の樹冠の領域と接するか、或いは、輝度の変化率がある閾値以下になったとき、ここで領域の拡大を終了し、領域を閉じる。この領域の外殻を樹冠の境界線とすることにより、樹冠形状T1乃至T6を抽出している。抽出された樹冠形状T1乃至T6をポリゴン化して、図3(a)に示されるように、樹冠形状ポリゴンP1乃至P6が生成される。
【0031】
次に、樹種の取得について説明する。樹木毎の樹種データを取得するには、個々の樹木に対応して生成された樹冠形状ポリゴン内の波長スペクトルを分析すればよいので、先ず、取得された森林画像に、生成された各樹冠形状ポリゴンを重ね合わせ、そして、各樹冠形状ポリゴン内においてスペクトル分析を行う。
【0032】
一般に、森林地域には、常緑樹木、落葉樹木等が混在し、取得したカラー画像中の各樹木に関するスペクトルは、例えば、季節によっても異なる。しかし、このような森林地域のカラー画像に現れる「季節による変動データ」に関しては、適切な期間にわたって蓄積して検証することにより、樹冠形状と、樹冠形状を求めたカラー画像のスペクトル分析とに基づいて、正確な樹木特定が可能である。更に、画像取得時の太陽光線の相違によってもデータ変動が生じるが、これについても撮影データを蓄積して検証することにより解決される。
【0033】
この様にして得られた樹木に関するスペクトル分布図の一例を、図2に示した。同図では、横軸が光波長を、縦軸が輝度を表し、樹種に対応したスペクトル分布が見られる。例えば、スギ(杉)、ヒノキ(櫓)、コナラ(小楢)、カシ(樫)等では、スギの樹冠部分の輝度は、他の樹種に比較して青から近赤外線にいたるスペクトル全体にわたって総じて低く、ヒノキの樹冠部分は、青、緑、赤の可視光線ではスギと同様に輝度が低いが、近赤外線領域では輝度が高くなっている。一方、コナラはスギ及びヒノキに比べて全体的に輝度が高く、カシは、コナラと同様に全体的に輝度が高いが、コナラに比較して緑及び赤のスペクトルが低くなっている。
【0034】
この様に、樹種によって、スペクトル分布が異なっているので、これらの特性を参照して、作成された樹冠形状ポリゴン内のスペクトル分布を分析することにより、当該樹冠形状ポリゴンに係る樹木の樹種を特定することができる。また、森林画像情報が、マルチ/ハイパースペクトルセンサーによって取得されたものであれば、赤外領域のスペクトルも含め詳細なスペクトル分布を取得でき、より正確な樹種の特定が可能となる。更には、同種類の樹木である場合には、当該樹木の樹冠の大小を参照すると、幹の直径、樹高、樹齢を推定することも可能であり、高木に限定されることなく、低木の分布評価にも適用可能である。
【0035】
これまで述べてきたように、高解像度人工衛星画像あるいは航空機による森林画像を用いて立木本数を推定し、樹木毎の樹種を決定し、森林簿や現地調査で計測した樹高情報を組み合わせることにより、材積を精度良く推定できる。しかし、画像センサーだけでなく、対象地域の樹高を計測することのできる、航空機搭載型のレーザー距離計(LiDAR)や合成開口レーダ(SAR)等の距離計測センサーと組み合わせ、マルチセンシングプラットフォームとすることで、森林情報を総合的に取得することができる。
【0036】
ここで、材積を計算するために必要な樹木毎の樹高Hを取得する処理手順について説明する。
【0037】
衛星や航空機といった飛行体によるセンシングでは、GPSと慣性計測装置(IMU)を併用した高精度位置姿勢計測技術(GPS/IMU)が用いられる。GPS/IMUによれば、センサーのイベントタグ(データ取得タイミング)と同期してセンサーのX、Y、Z軸情報と同時にセンサーの3軸姿勢情報が取得できるので、センシングデータの幾何学的な歪補正が簡単に行える。GPS/IMUを備えたマルチセンサー対応のセンシングプラットフォームとすることにより、森林画像だけでなく、様々な情報を一度に計測することができる。
【0038】
なお、森林画像情報と樹高情報とが、同一飛行体に搭載された各センサーで同時に取得されなくとも、GPS/IMUを用いることによって、取得点とセンサーの傾きとが分かり、それらの情報から、取得対象である樹木の地上位置座標を特定できるので、森林画像情報と樹高情報とが別々に取得されても、各位置情報を知ることができ、森林画像情報と樹高情報とを重ね合せるときの各々の位置合せが、簡単に実施できる。
【0039】
樹高を測定するために、距離計測センサーとしてレーザー距離計を用いた場合では、例えば、同一地点からのレーザー反射を5パルス観測することができるものを使用する。5パルス計測で取得された5種類の高さ情報に従って、樹高を計測することができる。
【0040】
5つのパルスのうち、1番目のパルスは、プラットフォームの樹冠頂点からの高さを示している。そして、最後に取得した5番目のパルスは、地表面からのものであり、プラットフォームの地表面からの高さを示している。従って、最初と最後のパルスで計測された高さ情報による差が、樹高Hを表わしていることになる。
【0041】
この様に計測された樹高データは、撮影した森林画像に映し出されている樹木に係る樹高点群を表しているので、この樹高点群から各樹木の樹高データを抽出する必要がある。そこで、前述した、森林画像から抽出された樹冠形状を樹高データに重ねあわせることにより、当該樹冠形状に対応して樹木毎の樹高データを求めることができる。なお、合成開口レーダを使用した場合は、各画素が距離を表す画像(樹高画像)が得られる。
【0042】
図3及び図4に示した(a)乃至(e)を参照しながら、飛行体搭載の画像センサーから撮影された森林画像と、同飛行体搭載の距離計測センサーで取得された樹高点群又は樹高画像による樹高情報とから、樹木毎の樹高データを生成する本実施形態における具体的処理について説明する。
【0043】
飛行体搭載の距離計測センサーによって、森林画像と同時に、該森林画像と同じ位置に関する樹高点群又は樹高画像による樹高情報が取得されている。図3(b)には、前述のレーザー距離計で測定された樹高点群の例を示した。同図では、X−Z平面で、つまり、図3(a)におけるX軸方向で1ラインをスキャンしたときの高さ情報をドットで表している。通常、飛行体搭載のレーザー距離計では、プラットフォームの測定対象からの高さ を、測定対象のジオイド(一般に、平均海水面)からの高さに変換した値を出力する。Hsは、最初のパルスで測定されたジオイドからの樹冠の高さを示しており、樹冠の高さを表している。そして、Hgは、最後のパルスで計測されたジオイドからの地表面の高さを表している。ここから、各ドットに係る樹高をHtとすれば、Ht=Hs−Hgの関係があり、各ドットの樹高Htが求められる。
【0044】
次いで、図3(c)に示されるように、求められた各ドットに係る樹高Htを図3(a)の樹冠形状ポリゴンP1乃至P6とXY座標に合せて重ねる。そこで、樹高Htの樹高データ群から、樹冠形状ポリゴンP1乃至P6の内側に含まれる樹高データhxyの群を抽出する。その様子が、模式的に図4(d)に示される。同図では、樹冠形状ポリゴンP4の場合を示した。なお、当該樹冠形状に係る樹高Htについて、図3(c)の段階で求めた場合を説明したが、図4(d)に示されるように、樹冠形状ポリゴンP4で絞られたドットのHs及びHgに対する樹高データからHtを求めても良い。
【0045】
樹冠形状ポリゴンP1乃至P6の内側に含まれる樹高データhxyの群が抽出できたならば、各樹冠形状ポリゴンP1乃至P6内において、含まれる樹高データhxyのうちで最大値のものを選択する。この最大値を有するドットが樹冠の頂点に対応していると評価でき、これを樹高とする。その最大値が選択された状態を、図4(e)に示した。図中では、その最大値を有するドットの位置に「・」印を付して、樹高点を示した。樹冠形状ポリゴンP4の場合には、樹高データh44が示されている。
【0046】
このようにして、飛行体搭載の画像センサーによる森林画像から生成した樹冠形状ポリゴンと、飛行体搭載の距離計測センサーから取得された樹高情報から、樹木毎の樹高データを取得できた。
【0047】
以上によって、飛行体に搭載された画像センサーから取得された森林画像情報と、距離計測センサーによって取得された樹高情報とから、当該森林地域における立木本数、該森林地域内にある樹木毎の樹高データを計測し、そして、該樹木の樹種を特定して該当係数を選定することができた。従って、立木本数N、樹高データH及び樹種に係る係数Bによって、上述の材積推定式から材積を求めることができる。
【0048】
次に、所定の森林地域に関する材積を求める全体的な処理手順を説明するフローチャートを、図5に示した。図5のフローチャートは、大きく分けて、森林観測システムと地上解析システムとの処理を含んでいる。
【0049】
森林観測システムの処理は、上述したセンシングプラットフォーム内で実行される。本実施形態の森林情報処理システムでは、この森林観測システムには、飛行体に搭載された画像センサー、例えば、マルチ/ハイパースペクトルセンサーによって、所定高度から、該当森林地域についての森林画像情報が取得されるステップS1と、飛行体に搭載されたレーザー距離計又は合成開口レーダによって、該森林地域に関する樹高点群又は樹高画像による樹高情報が取得されるステップS2とが含まれている。
【0050】
地上解析システムでは、森林観測システムで取得された森林画像情報と樹高情報に基づいて、森林材積の計算が行われる。この森林材積を計算する処理のために、地上解析システムは、樹冠形状抽出手段、樹種分類手段、樹高計測手段、材積計算手段(いずれも図示せず)から構成されている。
【0051】
樹冠形状抽出手段1では、先ず、ステップS1において取得された森林画像情報に、ウォーターシェッド・アルゴリズムを適用して、該画像に映し出されている樹冠に係る領域を抽出し、該領域の外殻線を求めることにより、図1(b)に示されるように、樹冠境界線の強調処理を行う(ステップS3)。
【0052】
次いで、強調処理された樹冠境界線に基づいて樹冠形状ポリゴンP1乃至P6を抽出する(ステップS4)。そして、これらの樹冠形状ポリゴンの各重心点を求め、当該森林画像中において、各重心点をカウントすることにより、該森林地域に含まれる立木本数Nを取得する(ステップS5)。
【0053】
次に、樹種分類手段2では、ステップS1において取得された森林画像に、ステップS4において抽出された樹冠形状ポリゴンP1乃至P6を重ね合わせる(ステップS6)。
【0054】
そこで、重ね合わされた各樹冠形状ポリゴン内に含まれる森林画像の色情報について分析処理が行われる(ステップS7)。そして、分析された樹冠形状ポリゴン毎の色情報に基づいて、樹種に関する色統計情報を参照し、当該樹冠形状ポリゴンに対応する樹木の樹種が特定され、樹木毎の樹種データが取得される(ステップS8)。
【0055】
一方、樹高計測手段3では、ステップS2において取得された樹高点群又は樹高画像の樹高情報と、ステップS4において取得された樹冠形状ポリゴンP1乃至P6とを重ね合せ(ステップS9)、各樹冠形状ポリゴン内に含まれる樹高情報の樹木反射点データを抽出する(ステップS10)。ここで、上述した5パルスによるレーザー計測の例であれば、樹木反射点は、最初のパルスによるデータであってもよい。
【0056】
そして、抽出された各樹冠形状ポリゴン内の反射点データのうちで、最も高い樹高値を示す樹木反射点データを、樹冠頂点を示すものとして抽出する(ステップS11)。この抽出された樹木最高値の樹木反射点データが、樹冠部頂点を示すものとして取得される。次に、樹高情報Hsと地表面高情報Hgとから樹木毎の樹高Hが求められる(ステップS12)。
【0057】
以上のようにして、ステップS5において立木本数N、ステップS8において樹木毎の樹種データ、そして、ステップS12において樹木毎の樹高Hがそれぞれ取得されたので、材積計算手段4では、上述した材積推定式を利用して、各樹種における樹木毎の材積Vを計算する(ステップS13)。このとき、同式における係数B1乃至B4については、現地調査結果を踏まえ、各森林地域に対応して決定されているので、材積計算時に樹種毎に選定して使用される。
【0058】
この様にして、飛行体に搭載された計測機器を用いて、森林地域の上空から撮影された森林画像情報と、上空から取得した樹高情報とから、当該森林地域の樹木材積を自動計測できる森林情報処理システムを構築できる。
【0059】
【発明の効果】
以上のように、本発明の森林情報処理システムによれば、該当森林地域の上空から取得された森林画像情報と、同じく上空から計測された樹高情報とに基づいて、立木本数、樹種、そして、樹木毎の樹高を把握し、樹木毎の材積を計算できるようにしたので、森林地域の材積計算を自動化でき、しかも、科学的、客観的に判定することが可能となった。
【0060】
さらに、地球温暖化対策として、全国的に広い範囲について樹木材積の調査を行わなければならないが、本発明の森林情報処理システムは、この広い範囲に対しても、調査対象とすることができる。そして、各地域の森林特性や、判読者の技術レベルの差といった計測誤差の要因を回避することが可能となった。
【0061】
また、現地調査を省略することができるため、森林計測コストの削減を図ることができ、現地調査者の安全性を担保できる効果もある。
【図面の簡単な説明】
【図1】 森林地域における各樹木の樹冠と樹冠形状を説明する図である。
【図2】 樹木の樹種に対応したスペクトル分布の一例を表した図である。
【図3】 樹高計測を行う処理手順を説明する図である。
【図4】 図3に続く樹高計測を行う処理手順を説明する図である。
【図5】 本実施形態による森林情報処理システムに係る処理手順を示すフロー図である。
【符号の説明】
T1〜T6…樹木
P1〜P6…樹冠形状ポリゴン
Hs…樹高情報
Hg…地表面高情報
hxy…樹高データ
Claims (4)
- 森林地域を上空から取得した森林画像情報にウォーターシェッド・アルゴリズムを適用して樹冠形状情報を抽出する樹冠形状抽出手段と、
前記森林地域の上空から計測した樹高情報と前記森林画像情報を重ね合せ、前記樹冠形状の範囲内に含まれる前記樹高情報のうち最大値を抽出して当該樹冠の樹高とする樹高計測手段と、
前記森林画像内で生成された前記樹冠形状情報に基づいて、前記森林地域内に形成された前記樹冠形状の重心点を求めて該重心点を計数し、該計数値を前記森林地域内の立木本数として取得する立木本数計測手段と、
前記森林画像情報に含まれる色情報を前記樹冠形状に基づいて抽出し、色統計情報を参照して当該樹冠の樹種を取得する樹種分類手段と、
前記樹高と前記立木本数とに基づいて、前記樹種に応じて樹木の材積を計算する材積計算手段とを有する森林情報処理システム。 - 前記材積計算手段は、前記樹木毎の材積を、材積推定式
V=(B1・HB2+B3・HB4/N)-1
〔但し、Vは材積(m3/ha)、Hは樹高(m)、Nは1haあたりの立木本数であり、B1乃至B4は樹種及び森林地域に従って決定される係数である〕に従って計算することを特徴とする請求項1に記載の森林情報処理システム。 - 前記森林画像情報は、当該取得位置情報と共に取得され、前記樹高情報は、当該計測位置情報と共に計測され、前記森林画像情報と前記樹高情報との重ね合せは、該取得位置情報と該計測位置情報とに基づいて行われることを特徴とする請求項1又は2に記載の森林情報処理システム。
- 前記森林画像情報と前記樹高情報は、同一飛行体に搭載された画像取得手段と距離計測手段によりそれぞれ取得されることを特徴とする請求項1乃至3のいずれか一項に記載の森林情報処理システム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002147774A JP4004332B2 (ja) | 2002-05-22 | 2002-05-22 | 森林情報処理システム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002147774A JP4004332B2 (ja) | 2002-05-22 | 2002-05-22 | 森林情報処理システム |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2003344048A JP2003344048A (ja) | 2003-12-03 |
JP2003344048A5 JP2003344048A5 (ja) | 2005-10-06 |
JP4004332B2 true JP4004332B2 (ja) | 2007-11-07 |
Family
ID=29766675
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002147774A Expired - Lifetime JP4004332B2 (ja) | 2002-05-22 | 2002-05-22 | 森林情報処理システム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4004332B2 (ja) |
Families Citing this family (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006285310A (ja) * | 2005-03-31 | 2006-10-19 | Kanazawa Univ | 森林の樹冠の評価方法及びその樹冠評価プログラム |
US7411196B2 (en) * | 2005-08-18 | 2008-08-12 | Itt Manufacturing Enterprises, Inc. | Multi-sensors and differential absorption LIDAR data fusion |
JP2007107962A (ja) * | 2005-10-12 | 2007-04-26 | Yamamoto Denko Kk | 測定対象物計測システム、及び送電線下離隔計測システム |
JP5034052B2 (ja) * | 2005-12-15 | 2012-09-26 | 国立大学法人山口大学 | 樹体生産能力を評価する方法及び樹体生産能力を評価するための撮像装置 |
JP4946072B2 (ja) * | 2006-01-24 | 2012-06-06 | 朝日航洋株式会社 | 樹頂探索方法、装置及びプログラム |
JP2008210165A (ja) * | 2007-02-27 | 2008-09-11 | Tohoku Kensetsu Kyokai | 植生判定方法及び画像処理装置 |
CN101978395B (zh) * | 2008-04-23 | 2012-10-03 | 株式会社博思科 | 建筑物屋顶轮廓识别装置及建筑物屋顶轮廓识别方法 |
JP4900356B2 (ja) * | 2008-09-30 | 2012-03-21 | 国立大学法人信州大学 | 樹種分類方法及び樹種分類システム、森林現況情報の作成方法及び森林現況情報の作成システム、間伐対象区域の選定方法及び間伐対象区域の選定システム |
KR100935857B1 (ko) | 2009-08-25 | 2010-01-07 | 홍종옥 | 항공라이다와 디지털항공사진을 이용한 3차원 산림지리정보 생성 시스템 및 그 방법 |
JP4858793B2 (ja) * | 2009-11-12 | 2012-01-18 | 国立大学法人信州大学 | 樹木本数算定方法及び樹木本数算定装置 |
JP5520786B2 (ja) * | 2010-11-11 | 2014-06-11 | 株式会社パスコ | レーザ密度分布推定装置、レーザ密度分布推定方法、及びプログラム |
JP2013088188A (ja) * | 2011-10-14 | 2013-05-13 | Fuji Architect Co Ltd | 三次元測定対象物の形態調査方法 |
JP6273473B2 (ja) * | 2012-11-20 | 2018-02-07 | 小平アソシエイツ株式会社 | 森林情報管理装置 |
JP6200821B2 (ja) * | 2014-01-29 | 2017-09-20 | 株式会社パスコ | 林相解析装置、林相解析方法及びプログラム |
JP6200826B2 (ja) * | 2014-02-12 | 2017-09-20 | 株式会社パスコ | 林相解析装置、林相解析方法及びプログラム |
JP6303764B2 (ja) * | 2014-04-23 | 2018-04-04 | 日本電気株式会社 | データ融合装置、土地被覆分類システム、方法およびプログラム |
KR101545479B1 (ko) | 2015-03-31 | 2015-08-24 | 김명준 | 무인항공기 정사영상을 이용한 산림 병해충 피해목 검출 방법 |
US10546195B2 (en) * | 2016-12-02 | 2020-01-28 | Geostat Aerospace & Technology Inc. | Methods and systems for automatic object detection from aerial imagery |
JP7123366B2 (ja) * | 2018-02-26 | 2022-08-23 | 株式会社 ジツタ | 森林地域での立木の評価方法、及びこの評価方法における評価対象エリアの特定に適した境界線測量方法 |
CN109064556B (zh) * | 2018-08-10 | 2022-07-08 | 四川大学 | 一种面向isr的地貌高精度仿真建模*** |
CN109343093A (zh) * | 2018-10-10 | 2019-02-15 | 北华大学 | 森林样地内上层树木树高精确测量方法 |
JP7169640B2 (ja) * | 2018-11-16 | 2022-11-11 | 国立大学法人信州大学 | 間伐木の選定方法及び間伐木の選定装置 |
JP2020091640A (ja) * | 2018-12-05 | 2020-06-11 | 国立大学法人京都大学 | 物体分類システム、学習システム、学習データ生成方法、学習済モデル生成方法、学習済モデル、判別装置、判別方法、およびコンピュータプログラム |
JP7298823B2 (ja) * | 2019-08-22 | 2023-06-27 | 株式会社 ジツタ | 森林地域での立木の評価方法、及びこの評価方法における評価対象エリアの特定に適した境界線測量方法 |
CN114279431A (zh) * | 2021-11-22 | 2022-04-05 | 北京林业大学 | 一种用于密郁闭度森林的林木位置测量方法及其装置 |
JP7374977B2 (ja) * | 2021-12-17 | 2023-11-07 | ヤマハ発動機株式会社 | 樹木情報推定システム、樹木情報推定方法およびコンピュータプログラム |
CN116485822B (zh) * | 2023-04-28 | 2023-12-08 | 中建一局集团第三建筑有限公司 | 一种高郁闭度人工林单木树冠分层分割方法 |
-
2002
- 2002-05-22 JP JP2002147774A patent/JP4004332B2/ja not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
JP2003344048A (ja) | 2003-12-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4004332B2 (ja) | 森林情報処理システム | |
Lim et al. | Calculation of tree height and canopy crown from drone images using segmentation | |
Popescu et al. | Measuring individual tree crown diameter with lidar and assessing its influence on estimating forest volume and biomass | |
US8208689B2 (en) | Method for determination of stand attributes and a computer program for performing the method | |
JP4900356B2 (ja) | 樹種分類方法及び樹種分類システム、森林現況情報の作成方法及び森林現況情報の作成システム、間伐対象区域の選定方法及び間伐対象区域の選定システム | |
Gini et al. | Aerial images from an UAV system: 3D modeling and tree species classification in a park area | |
Ni et al. | Mapping three-dimensional structures of forest canopy using UAV stereo imagery: Evaluating impacts of forward overlaps and image resolutions with LiDAR data as reference | |
Zhang et al. | Estimation of forest leaf area index using height and canopy cover information extracted from unmanned aerial vehicle stereo imagery | |
Valbuena et al. | Comparing airborne laser scanning-imagery fusion methods based on geometric accuracy in forested areas | |
Coluzzi et al. | On the LiDAR contribution for landscape archaeology and palaeoenvironmental studies: the case study of Bosco dell'Incoronata (Southern Italy) | |
Meneses et al. | Modelling heights of sparse aquatic reed (Phragmites australis) using Structure from Motion point clouds derived from Rotary-and Fixed-Wing Unmanned Aerial Vehicle (UAV) data | |
Raup et al. | Remote sensing of glaciers | |
Durrieu et al. | Optical remote sensing of tree and stand heights | |
Yin et al. | Individual tree parameters estimation for chinese fir (cunninghamia lanceolate (lamb.) hook) plantations of south china using UAV Oblique Photography: Possibilities and Challenges | |
Yilmaz et al. | Determination of tree crown diameters with segmentation of a UAS-Based canopy height model | |
CN115760885B (zh) | 基于消费级无人机影像的高郁闭度湿地森林参数提取方法 | |
Hosingholizade et al. | Height estimation of pine (Pinus eldarica) single trees using slope corrected shadow length on unmanned aerial vehicle (UAV) imagery in a plantation forest | |
Adhikari et al. | Does topographic normalization of landsat images improve fractional tree cover mapping in tropical mountains? | |
Viau et al. | The use of airborne LIDAR and multispectral sensors for orchard trees inventory and characterization | |
Jiménez-Bello et al. | Use of remote sensing and geographic information tools for irrigation management of citrus trees | |
Pepe | Use of digital aerial photogrammetry sensors for land cover classification | |
Li et al. | High-throughput plant height estimation from RGB Images acquired with Aerial platforms: a 3D point cloud based approach | |
Huang et al. | Experiment on extracting forest canopy height from Worldview-2 | |
Kim et al. | Identification of tree species from high-resolution satellite imagery by using crown parameters | |
Pang et al. | An airborne multi-angle hyperspectral experiment in a boreal forest of Northeast China |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050519 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050519 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070206 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070406 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070508 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070703 |
|
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: 20070724 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070821 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4004332 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100831 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100831 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110831 Year of fee payment: 4 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120831 Year of fee payment: 5 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130831 Year of fee payment: 6 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
EXPY | Cancellation because of completion of term |