JP2014191796A - 地物表面検出装置、地物表面検出方法、及びプログラム - Google Patents

地物表面検出装置、地物表面検出方法、及びプログラム Download PDF

Info

Publication number
JP2014191796A
JP2014191796A JP2013069520A JP2013069520A JP2014191796A JP 2014191796 A JP2014191796 A JP 2014191796A JP 2013069520 A JP2013069520 A JP 2013069520A JP 2013069520 A JP2013069520 A JP 2013069520A JP 2014191796 A JP2014191796 A JP 2014191796A
Authority
JP
Japan
Prior art keywords
point
feature
cutting plane
boundary
point group
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.)
Granted
Application number
JP2013069520A
Other languages
English (en)
Other versions
JP6013962B2 (ja
Inventor
Shizuo Manabe
鎮男 真鍋
Ikuo Kitagawa
育夫 北川
Hirotsugu Yoshioka
裕嗣 吉岡
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Pasco Corp
Original Assignee
Pasco Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Pasco Corp filed Critical Pasco Corp
Priority to JP2013069520A priority Critical patent/JP6013962B2/ja
Publication of JP2014191796A publication Critical patent/JP2014191796A/ja
Application granted granted Critical
Publication of JP6013962B2 publication Critical patent/JP6013962B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

【課題】地物表面から抽出された点群の3次元座標データに基づいて植生が繁茂した地物表面の形状を精度良く自動的に検出する。
【解決手段】点群抽出手段20は、鉛直な切断平面から近傍距離内に位置する注目点群を抽出する。境界探索手段22は注目点群の切断平面への射影像である射影点群に基づいて、切断平面内にて地物表面により形成される境界線を探索する。単位領域設定手段30は切断平面を格子状に分割して複数の単位領域を設定する。下限単位領域抽出手段32は単位領域の各鉛直列を下から上へ順番に調べて、射影点群が最初に見つかった単位領域を下限単位領域として抽出する。特徴点抽出手段34は下限単位領域内の射影点群のうち高さが最も低い点を特徴点として抽出する。境界決定手段38は単位領域の鉛直列ごとの特徴点を節点に有する折線を境界線として求める。
【選択図】図1

Description

本発明は、レーザスキャナによって地物表面から抽出された点群の3次元座標データに基づいて植生が繁茂した堤防や道路沿線の法面の表面を検出する地物表面検出装置、地物表面検出方法、及びプログラムに関する。
河川堤防は主に土で構成されている。一般にその天端は河川管理用通路として利用され、必要に応じて敷砂利やアスファルト舗装が施されている。一方、法面は表面の崩落を防ぐために通常は芝を生やしている。ちなみに、雑草が繁茂し法面が弱体化するのを抑制するために、あるいは堤防の定期的な巡視点検のために年に1,2回程度の除草管理が行われている。
また、河川堤防は流水による浸食や水の浸透による決壊や、地盤沈下などによる越水現象が生じないように管理する必要がある。当該管理において、堤防形状の変化の監視が有効である。
ここで、地物の形状を計測する技術として、特許文献1には、レーザスキャナを用いて、地物の形状を表す3次元点群データを取得する技術が示されている。例えば、モービルマッピングシステムでは、車両に搭載したレーザスキャナを用い道路に沿って地物の形状を表す3次元点群データを取得する。当該システムでは、自動車に搭載されたレーザスキャナは車体の上部から斜め下方向や斜め上方向にレーザを照射する。レーザの光軸は横方向に走査され、走査角度範囲内にて微小角度ごとにレーザパルスが発射される。レーザの発射から反射光の受信までの時間に基づいて距離が計測され、またその際、レーザの発射方向、時刻、及び車体の位置・姿勢などが計測される。それら計測データから、レーザパルスを反射した点の3次元座標を表す点群データが求められる。
特開2009−204615号公報
従来、点群データをもとに地物を判読するためには人手を要しており、3次元CADで編集ツール等を利用して手作業で地物の抽出作業が行われていた。そのため、堤防から点群データを取得し、それに基づいて堤防形状を抽出するのに手間がかかるという問題があった。特に、法面等の植生によるレーザ反射は、堤防形状を抽出する上でノイズとなる点群を生じるという問題がある。
本発明は、地物表面から抽出された点群の3次元座標データに基づいて植生が繁茂した地物表面の形状を精度良く自動的に検出する地物表面検出装置、地物表面検出方法、及びプログラムを提供することを目的とする。
本発明に係る地物表面検出装置は、レーザスキャナによって地物表面から抽出された点群の3次元座標データに基づき、対象空間に存在する植生が繁茂した地物表面を検出するものであって、前記対象空間に鉛直面である切断平面を仮想的に設定し、前記点群のうち当該切断平面から予め設定した近傍距離内に位置する注目点群を抽出する点群抽出手段と、前記注目点群の前記切断平面への射影像である射影点群に基づいて、前記切断平面内にて前記地物表面により形成される境界線を探索する境界探索手段と、を有し、前記境界探索手段は、前記切断平面を水平方向及び鉛直方向に沿って格子状に分割して複数の単位領域を設定する単位領域設定手段と、前記切断平面における前記単位領域の各鉛直列を下から上へ順番に調べて、前記射影点群が最初に見つかった前記単位領域を下限単位領域として抽出する下限単位領域抽出手段と、前記下限単位領域内の前記射影点群のうち高さが最も低い点を特徴点として抽出する特徴点抽出手段と、前記単位領域の鉛直列ごとの前記特徴点を節点に有する折線を前記境界線として求める境界決定手段と、を有する。
他の本発明に係る地物表面検出装置はさらに、前記射影点群のうち、隣接する前記特徴点同士を結ぶ線分からなる暫定境界線より高さが低い点を追加特徴点として抽出する追加特徴点抽出手段を有し、前記境界決定手段は、前記特徴点及び前記追加特徴点を節点に有する折線を前記境界線として求める。
さらに他の本発明に係る地物表面検出装置においては、前記追加特徴点抽出手段は、隣接する前記特徴点間にて前記暫定境界線を複数区間に分割し、当該区間ごとに、矩形領域でありその上辺を当該区間内における前記暫定境界線の最高点に応じた高さに設定した探索領域内にて前記追加特徴点を探索する。
本発明に係る地物表面検出方法は、レーザスキャナによって地物表面から抽出された点群の3次元座標データに基づき、対象空間に存在する植生が繁茂した地物表面を検出する方法であって、前記対象空間に鉛直面である切断平面を仮想的に設定し、前記点群のうち当該切断平面から予め設定した近傍距離内に位置する注目点群を抽出する点群抽出ステップと、前記注目点群の前記切断平面への射影像である射影点群に基づいて、前記切断平面内にて前記地物表面により形成される境界を探索する境界探索ステップと、を有し、前記境界探索ステップは、前記切断平面を水平方向及び鉛直方向に沿って格子状に分割して複数の単位領域を設定する単位領域設定ステップと、前記切断平面における前記単位領域の各鉛直列を下から上へ順番に調べて、前記射影点群が最初に見つかった前記単位領域を下限単位領域として抽出する下限単位領域抽出ステップと、前記下限単位領域内の前記射影点群のうち高さが最も低い点を特徴点として抽出する特徴点抽出ステップと、前記単位領域の鉛直列ごとの前記特徴点を節点に有する折線を前記境界として求める境界決定ステップと、を有する。
本発明に係るプログラムは、コンピュータに、レーザスキャナによって地物表面から抽出された点群の3次元座標データに基づき、対象空間に存在する植生が繁茂した地物表面を検出するデータ解析を行わせるためのプログラムであって、当該コンピュータを、前記対象空間に鉛直面である切断平面を仮想的に設定し、前記点群のうち当該切断平面から予め設定した近傍距離内に位置する注目点群を抽出する点群抽出手段、及び、前記注目点群の前記切断平面への射影像である射影点群に基づいて、前記切断平面内にて前記地物表面により形成される境界を探索する境界探索手段、として機能させ、前記境界探索手段は、前記切断平面を水平方向及び鉛直方向に沿って格子状に分割して複数の単位領域を設定する単位領域設定手段と、前記切断平面における前記単位領域の各鉛直列を下から上へ順番に調べて、前記射影点群が最初に見つかった前記単位領域を下限単位領域として抽出する下限単位領域抽出手段と、前記下限単位領域内の前記射影点群のうち高さが最も低い点を特徴点として抽出する特徴点抽出手段と、前記単位領域の鉛直列ごとの前記特徴点を節点に有する折線を前記境界として求める境界決定手段と、を備える。
本発明によれば、地物表面から抽出された点群の3次元座標データに基づいて地物表面の形状を精度良く自動的に検出することができる。
本発明の実施形態に係る地物表面検出システムの概略の構成を示すブロック図である。 堤防の横断方向に設定した切断平面を示す模式的な斜視図である。 本発明の実施形態に係る地物表面検出システムによる堤防表面形状抽出処理の概略のフロー図である。 切断平面における堤防の形状及び単位領域の配置の一例を示す模式図である。 堤防の一部分における射影点群の分布の一例を示す模式図である。 本発明の実施形態に係る地物表面検出システムによる堤防表面の境界線検出の処理フロー図である。 下限単位領域の例を示す切断平面の一部分の模式図である。 特徴点の例を示す切断平面の一部分の模式図である。 探索領域の設定例を示す切断平面の一部分の模式図である。 境界線の例を示す切断平面の一部分の模式図である。
以下、本発明の実施の形態(以下実施形態という)である地物表面検出システム2について、図面に基づいて説明する。本システムは、地物表面の3次元形状を表す点群データに基づいて例えば、河川堤防の表面形状を検出する。点群データは例えば、上述のモービルマッピングシステムのように地上を走行する車両に搭載されたレーザスキャナにより取得される。また、レーザスキャナを地上に設置して計測を行っても良い。点群データが地物表面の3次元形状を表すには、地物表面の凹凸、段差のスケールに応じた密度でレーザスキャンが行われる必要がある。この点、車両や三脚等の高さから行うレーザスキャンの走査密度、距離精度及びフットプリントの大きさは、例えば、地物の形状を数センチメートル程度の精度で捉えることができる諸元を有し、河川堤防の法面の変形を航空レーザや衛星レーザなどに比べて高精度に検出することが可能である。
図1は、地物表面検出システム2の概略の構成を示すブロック図である。本システムは、演算処理装置4、記憶装置6、入力装置8及び出力装置10を含んで構成される。演算処理装置4として、本システムの各種演算処理を行う専用のハードウェアを作ることも可能であるが、本実施形態では演算処理装置4は、コンピュータ及び、当該コンピュータ上で実行されるプログラムを用いて構築される。
当該コンピュータのCPU(Central Processing Unit)が演算処理装置4を構成し、後述する点群抽出手段20及び境界探索手段22として機能する。境界探索手段22は単位領域設定手段30、下限単位領域抽出手段32、特徴点抽出手段34、追加特徴点抽出手段36及び境界決定手段38を含んで構成される。
記憶装置6はコンピュータに内蔵されるハードディスクなどで構成される。記憶装置6は演算処理装置4を点群抽出手段20及び境界探索手段22として機能させるためのプログラム及びその他のプログラムや、本システムの処理に必要な各種データを記憶する。例えば、記憶装置6は、処理対象データとして解析の対象空間の点群データを格納される。地物表面検出システム2は基本的には堤防が存在する空間を解析の対象空間に設定され利用される。特に、地物表面検出システム2は例えば、河川堤防の法面などのように植生が存在する部分の表面検出に有効である。
入力装置8は、キーボード、マウスなどであり、ユーザが本システムへの操作を行うために用いる。
出力装置10は、ディスプレイ、プリンタなどであり、本システムにより求められた堤防の3次元形状や断面形状を画面表示、印刷等によりユーザに示す等に用いられる。
後述するように本実施形態において堤防の表面形状を抽出する際の基本処理は、対象空間に仮想的に設定した切断平面に現れる地物表面の境界線(断面形状)を検出する処理である。堤防の3次元形状は、切断平面を複数設定してそれぞれについて解析を行うことによって把握することができる。
地物表面検出システム2において切断平面の設定の仕方は任意であるが、本実施形態では切断平面を堤防の横断方向に設定する場合を例に説明する。図2はこの場合を示す模式的な斜視図である。例えば、ユーザが基準線50を堤防52の横断方向に沿って設定すると、切断平面54は基準線50を通る鉛直面に設定される。基準線50は地物表面検出システム2により自動的に設定することもできる。また、例えば、ユーザが基準線50の最初の位置、間隔、本数や堤防の縦断方向における解析範囲などのパラメータを指定すると、地物表面検出システム2が指定された条件に基づいて複数の基準線を順次、自動設定して解析を行う構成とすることもできる。例えば、モービルマッピングシステムは天端の道路を走行しながらレーザ点群を取得することができる。そこで、モービルマッピングシステムにより取得される車両の走行軌跡に基づいて、天端に縦断方向に沿ったセンター線を設定し、これをある距離(ピッチ)で刻んで横断方向の基準線50を発生させて連続的にあるいは単独で断面形状を検出する処理を行うことができる。ここで、所定距離で刻む以外に、距離標(キロポスト)の値に基づいて刻んでも良い。また、車両の走行軌跡に基づくセンター線の代わりに既存の縦断センター線を用いてもよい。ここで、基準線に沿った水平方向をX軸、鉛直方向をY軸とすると切断平面はXY平面であり、Z軸は切断平面に垂直な方向に設定される。なお、Y軸の正の向きは上向きに設定する。
図3は、地物表面検出システム2による堤防表面形状抽出処理の概略のフロー図である。図3を参照しながら、演算処理装置4の各手段を説明する。
点群抽出手段20は対象空間に切断平面54を仮想的に設定し、点群のうち切断平面54から予め設定した近傍距離内に位置する注目点群を抽出する。具体的には、点群抽出手段20は切断平面54を中心としてZ軸方向に寸法(奥行き)Dを有する空間(以下、注目部分空間と称する。)を設定し、注目部分空間から注目点群を抽出する(S10)。ちなみに図2では、注目部分空間の水平断面56を矩形で模式的に表している。
境界探索手段22は、注目点群の切断平面54への射影像である射影点群に基づいて、切断平面54内にて堤防52の表面により形成される境界線58を探索する。
単位領域設定手段30は切断平面54を水平方向及び鉛直方向に区画して、X,Y各軸に沿った辺を有する矩形の領域(単位領域)を設定する(S12)。図4は単位領域を説明する模式図であり、切断平面54(XY平面)における堤防52の形状(境界線58)及び単位領域70の配置の一例を示している。また、図5は堤防52の一部分における射影点群60の分布の一例を示す模式図である。切断平面54は直交格子(メッシュ)状に区画され、単位領域70はX方向及びY方向を配列方向とする2次元配列をなす。本実施形態では単位領域のX軸方向の寸法(幅W)及びY軸方向の寸法(高さH)を等しくし、単位領域を正方形に設定する。なお、点群抽出手段20及び境界探索手段22を構成するプログラムでは、幅W、高さH及び奥行きDはパラメータ化されており、例えば、ユーザが入力装置8を用いて変更することができる。
ちなみに、図4におけるメッシュは図示の都合上、実際より粗く描いている。一方、図5は実際の射影点群及びメッシュの例を示している。図5の例では単位領域70は幅W、高さHがそれぞれ30cmの正方形である。
下限単位領域抽出手段32は切断平面54における単位領域70の鉛直列ごとに下から上へ順番に単位領域70を選択して当該単位領域70内における射影点群の有無を調べる。そして、下限単位領域抽出手段32は、単位領域70の鉛直列ごとに、射影点群が最初に見つかった単位領域を下限単位領域として抽出する(S14)。
特徴点抽出手段34は各下限単位領域内の射影点群のうち高さが最も低い点を特徴点として抽出する(S16)。
追加特徴点抽出手段36は、隣接する特徴点同士を結ぶ線分からなる折線を暫定的な境界線として求め(S18)、射影点群のうち、暫定の境界線より高さが低い点を追加特徴点として抽出する(S20)。
境界決定手段38は特徴点抽出手段34により抽出された特徴点及び追加特徴点抽出手段36により抽出された追加特徴点を節点に有する折線を境界線として求める(S22)。
以下、地物表面検出システム2の処理をより詳しく説明する。図6は地物表面検出システム2による堤防表面検出の概略の処理フロー図であり、切断平面に現れる堤防表面の境界線を検出する処理を示している。
点群抽出手段20は対象空間に切断平面54を仮想的に設定し、切断平面54を中心としてZ軸方向に寸法Dを有する注目部分空間内の点群を注目点群として、例えばハードディスク等の記憶装置6から演算処理装置4のRAM(Random Access Memory)等の作業領域に取り込む(S30)。
単位領域設定手段30は切断平面54をメッシュ状に分割して、単位領域70に区画する。例えば、注目部分空間のサイズが設定されている場合には、単位領域のサイズW,Hに応じて、X,Y各方向の単位領域の配列個数が定まり、例えば、各方向の単位領域の位置を示すインデックスの範囲が決定される。ここで、X,Y方向のインデックスをそれぞれI,I、またIの範囲を1〜N、Iの範囲を1〜Nと表す。ここで、IはY軸の正の向き、つまり下から上に順に増加するように設定される。
下限単位領域抽出手段32は単位領域の鉛直列を指定するインデックスIの値を0に初期設定し(S32)、Iをインクリメントし、単位領域の鉛直列を1列ずつ処理対象として選択する(S36)。
下限単位領域抽出手段32は、各列の処理開始に際しインデックスIの値を0に初期設定する(S38)。そして、下限単位領域抽出手段32は、Iをインクリメントして単位領域を下から順番に1つずつ選択し(S42)、注目点群をXY平面上へ射影した射影点群を含む単位領域のうち最も下に位置する下限単位領域を探索する。具体的には、選択した単位領域内に射影点があるか否かを調べ(S44)、射影点が存在しない場合には(S46にて「No」の場合)、順次、1つ上の単位領域が選択される(S42)。一方、選択した単位領域内に射影点が存在した場合(S46にて「Yes」の場合)、当該単位領域が下限単位領域となる。
図7は下限単位領域の例を示す切断平面54の一部分の模式図である。図7には、第k列〜第(k+3)列の単位領域70及び注目点群の射影像(射影点群80)の例が示され、また当該射影点群80に対し各単位領域列にて抽出された下限単位領域82が太線の矩形で示されている。
特徴点抽出手段34は下限単位領域内にて最も下に位置する射影点を探索し、当該射影点を処理対象としている第I列の単位領域列の特徴点P(I)として抽出する(S48)。特徴点P(I)が見つかると、下限単位領域抽出手段32は次の単位領域列の処理を始める(S36)。なお、I=Nの単位領域、つまり単位領域列の上端まで射影点が見つからなかった場合も(S40にて「Yes」の場合)、次の単位領域列の処理に移る(S36)。
図8は図7の例に対応して抽出される特徴点P(k)〜P(k+3)を示す切断平面54の一部分の模式図である。全ての単位領域列について特徴点の探索・抽出処理が完了すると(S34にて「Yes」の場合)、当該特徴点を節点とする折線が暫定境界線84として設定される(S50)。
追加特徴点抽出手段36は暫定境界線より下に位置する射影点を追加特徴点として抽出する(S52,S54)。具体的には、追加特徴点抽出手段36は、暫定境界線を構成する線分ごとに当該処理S50〜S54を行うことができる。下限単位領域抽出手段32は暫定境界線を構成する線分の両端となる特徴点対、すなわちX方向に関して隣接する特徴点の対を順次選択し(S50)、当該選択した2つの特徴点の間をX方向に関し複数区間に分割し、当該区間ごとに追加特徴点の探索領域を設定する(S52)。探索領域はX方向に沿った辺とY方向に沿った辺とからなる矩形領域に設定され、その上辺は当該区間内における暫定境界線の最高点に応じた高さに設定される。
図9は探索領域の設定例を示す切断平面54の一部分の模式図であり、図7及び図8の例に対応している。例えば、追加特徴点抽出手段36は特徴点間をX座標が小さい側(図において左側)から所定間隔で分割して区間を設定する。なお、この場合、特徴点間にて最も右側の区間は他の区間より小さくなり得る。例えば、特徴点P(k)とP(k+1)との間は暫定境界線上の点Q〜Qで分割され5つの区間[P(k),Q],[Q,Q],[Q,Q],[Q,Q],[Q,P(k+1)]が設定される。これら各区間に対応して探索領域90が設定される。例えば、区間[Q,Q]ではQが最高点となるので、当該区間の探索領域90の上辺の高さは基本的にQの高さに設定することができる。この設定の趣旨は、追加特徴点を抽出する暫定境界線より下側の領域は排除せずに、堤防表面の上側に多数存在する植生に起因する点群をできるだけ含まないように探索領域を設定することで、処理効率の向上を図ることにある。
なお、上辺の高さは当該趣旨からは区間の最高点に一致させることが好適であるが、当該趣旨に与える影響が許容される範囲内にて最高点からずれた高さに設定することもできる。また、上記趣旨からは、探索領域の上辺は必ずしも水平でなくてもよく、暫定境界線に沿って設定することが好適であり、そのように構成しても良い。一方、上述のように区間内で一定のY座標を有する上辺とする構成は、X方向の位置に応じて上辺の高さを計算する必要がない分、処理が簡素化される。
探索領域90の下辺の高さに関しては、例えば、特徴点P(k)とP(k+1)との間の各区間について、特徴点P(k)とP(k+1)とのうち低い方であるP(k+1)の高さ以下で任意に設定することができる。これは、特徴点が下限単位領域内の最下点であることから、特徴点P(k)とP(k+1)との間には、P(k+1)より低い点は存在せず、追加特徴点はP(k+1)の高さ以上にしか存在しないからである。
本実施形態では、探索領域90のY方向の寸法を暫定境界線の各線分及び当該線分毎の各区間にて画一的に設定しており、これにより処理の簡素化を図れる。そのため、図9に示す探索領域90の下辺の高さは暫定境界線の各線分の低い方の高さよりも十分低く設定している。この設定における探索領域90のY方向の必要寸法は堤防表面の想定される勾配と単位領域の幅Wとに基づいて定めることができる。例えば、河川堤防の法面の傾斜は原則として50%以下と定められているので、探索領域90のY方向の必要寸法は通常、単位領域の幅Wの半分程度にすれば足りる。
なお、特徴点P(k)とP(k+1)との間での追加特徴点の探索において、両特徴点のうち高い方であるP(k)が抽出された単位領域列に対応したX座標の範囲では追加特徴点の探索を省略してもよい。これは当該範囲では、P(k)が最下点であり暫定境界線より低い点は存在しないからである。
追加特徴点抽出手段36は上述のように設定した探索領域内にて、暫定境界線より高さが低い射影点を探索し、得られた射影点を追加特徴点とする(S54)。例えば、図9の例では点92が追加特徴点となる。
境界決定手段38は点群抽出手段20が設定した切断平面54での堤防表面の境界線として、特徴点抽出手段34及び追加特徴点抽出手段36により抽出された特徴点及び追加特徴点をX軸方向の並び順に従って接続した折線を生成する(S56)。図10は境界線58の例を示す切断平面54の一部分の模式図であり、図7〜図9の例に対応している。
境界決定手段38は例えば、求めた境界線のデータを基準線50(又は切断平面54)の位置・方向と関連づけて記憶装置6に登録する(S58)。
演算処理装置4は複数の切断平面54での境界線を求め、記憶装置6に蓄積する。そして、蓄積された複数の境界線のデータに基づいて堤防の3次元形状や断面形状を表す画像データを生成し、出力装置10であるディスプレイ、プリンタなどに出力する。
なお、上述の実施形態では、特徴点及び追加特徴点を節点に含む折線を境界線として求めたが、特徴点だけを節点に含む折線を境界線として求めても良い。
地物表面検出システム2は基準線50に沿った方向に関し単位領域の幅W程度の分解能で地物断面形状を抽出する。単位領域列は、図5に示した例のように、多数の点群を含むように設定され、単位領域列ごとに多数の点群のうちの最下点を特徴点として抽出する。ここで、レーザスキャナはX軸方向の任意の位置で堤防表面から反射点を得られるわけではなく、或るX座標での点群の最下点は植生からの反射点である可能性がある。すなわち、基準線に沿って微細に最下点を抽出し点群の下側の輪郭を求めても、当該輪郭は植生の反射点をノイズとして含んでおり、堤防表面の形状を好適に表しておらず、堤防表面の形状の変化の観察には適していない。これに対し、本システムの単位領域列ごとに抽出した最下点は堤防表面の高さを好適に示していることが期待できる。つまり、本システムが生成する境界線は植生によるノイズの影響を受けにくい。また、地上からのレーザスキャンにより取得される点群は航空レーザや衛星レーザなどに比べて高密度であるので、単位領域の幅Wは比較的小さくできる。よって、本システムによれば、ノイズの影響を受けにくく高精度に高さが表現された堤防表面の形状を高分解能で求めることが可能であり、堤防の変形を高精度に検出することができる。
また、上述の実施形態では、河川堤防の法面を対象に説明したが、これに限定されるものではなく、たとえば植生が繁茂した道路沿線の法面にも適用することができる。
2 地物表面検出システム、4 演算処理装置、6 記憶装置、8 入力装置、10 出力装置、20 点群抽出手段、22 境界探索手段、30 単位領域設定手段、32 下限単位領域抽出手段、34 特徴点抽出手段、36 追加特徴点抽出手段、38 境界決定手段、50 基準線、52 堤防、54 切断平面、58 境界線。

Claims (5)

  1. レーザスキャナによって地物表面から抽出された点群の3次元座標データに基づき、対象空間に存在する植生が繁茂した地物表面を検出する地物表面検出装置であって、
    前記対象空間に鉛直面である切断平面を仮想的に設定し、前記点群のうち当該切断平面から予め設定した近傍距離内に位置する注目点群を抽出する点群抽出手段と、
    前記注目点群の前記切断平面への射影像である射影点群に基づいて、前記切断平面内にて前記地物表面により形成される境界線を探索する境界探索手段と、を有し、
    前記境界探索手段は、
    前記切断平面を水平方向及び鉛直方向に沿って格子状に分割して複数の単位領域を設定する単位領域設定手段と、
    前記切断平面における前記単位領域の各鉛直列を下から上へ順番に調べて、前記射影点群が最初に見つかった前記単位領域を下限単位領域として抽出する下限単位領域抽出手段と、
    前記下限単位領域内の前記射影点群のうち高さが最も低い点を特徴点として抽出する特徴点抽出手段と、
    前記単位領域の鉛直列ごとの前記特徴点を節点に有する折線を前記境界線として求める境界決定手段と、
    を有することを特徴とする地物表面検出装置。
  2. 請求項1に記載の地物表面検出装置において、
    前記射影点群のうち、隣接する前記特徴点同士を結ぶ線分からなる暫定境界線より高さが低い点を追加特徴点として抽出する追加特徴点抽出手段を有し、
    前記境界決定手段は、前記特徴点及び前記追加特徴点を節点に有する折線を前記境界線として求めること、
    を特徴とする地物表面検出装置。
  3. 請求項2に記載の地物表面検出装置において、
    前記追加特徴点抽出手段は、隣接する前記特徴点間にて前記暫定境界線を複数区間に分割し、当該区間ごとに、矩形領域でありその上辺を当該区間内における前記暫定境界線の最高点に応じた高さに設定した探索領域内にて前記追加特徴点を探索すること、
    を特徴とする地物表面検出装置。
  4. レーザスキャナによって地物表面から抽出された点群の3次元座標データに基づき、対象空間に存在する植生が繁茂した地物表面を検出する地物表面検出方法であって、
    前記対象空間に鉛直面である切断平面を仮想的に設定し、前記点群のうち当該切断平面から予め設定した近傍距離内に位置する注目点群を抽出する点群抽出ステップと、
    前記注目点群の前記切断平面への射影像である射影点群に基づいて、前記切断平面内にて前記地物表面により形成される境界を探索する境界探索ステップと、を有し、
    前記境界探索ステップは、
    前記切断平面を水平方向及び鉛直方向に沿って格子状に分割して複数の単位領域を設定する単位領域設定ステップと、
    前記切断平面における前記単位領域の各鉛直列を下から上へ順番に調べて、前記射影点群が最初に見つかった前記単位領域を下限単位領域として抽出する下限単位領域抽出ステップと、
    前記下限単位領域内の前記射影点群のうち高さが最も低い点を特徴点として抽出する特徴点抽出ステップと、
    前記単位領域の鉛直列ごとの前記特徴点を節点に有する折線を前記境界として求める境界決定ステップと、
    を有することを特徴とする地物表面検出方法。
  5. コンピュータに、レーザスキャナによって地物表面から抽出された点群の3次元座標データに基づき、対象空間に存在する植生が繁茂した地物表面を検出するデータ解析を行わせるためのプログラムであって、当該コンピュータを、
    前記対象空間に鉛直面である切断平面を仮想的に設定し、前記点群のうち当該切断平面から予め設定した近傍距離内に位置する注目点群を抽出する点群抽出手段、及び、
    前記注目点群の前記切断平面への射影像である射影点群に基づいて、前記切断平面内にて前記地物表面により形成される境界を探索する境界探索手段、として機能させ、
    前記境界探索手段は、
    前記切断平面を水平方向及び鉛直方向に沿って格子状に分割して複数の単位領域を設定する単位領域設定手段と、
    前記切断平面における前記単位領域の各鉛直列を下から上へ順番に調べて、前記射影点群が最初に見つかった前記単位領域を下限単位領域として抽出する下限単位領域抽出手段と、
    前記下限単位領域内の前記射影点群のうち高さが最も低い点を特徴点として抽出する特徴点抽出手段と、
    前記単位領域の鉛直列ごとの前記特徴点を節点に有する折線を前記境界として求める境界決定手段と、
    を備えることを特徴とするプログラム。
JP2013069520A 2013-03-28 2013-03-28 地物表面検出装置、地物表面検出方法、及びプログラム Active JP6013962B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2013069520A JP6013962B2 (ja) 2013-03-28 2013-03-28 地物表面検出装置、地物表面検出方法、及びプログラム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2013069520A JP6013962B2 (ja) 2013-03-28 2013-03-28 地物表面検出装置、地物表面検出方法、及びプログラム

Publications (2)

Publication Number Publication Date
JP2014191796A true JP2014191796A (ja) 2014-10-06
JP6013962B2 JP6013962B2 (ja) 2016-10-25

Family

ID=51837931

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013069520A Active JP6013962B2 (ja) 2013-03-28 2013-03-28 地物表面検出装置、地物表面検出方法、及びプログラム

Country Status (1)

Country Link
JP (1) JP6013962B2 (ja)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016142535A (ja) * 2015-01-29 2016-08-08 田中 成典 測量データ処理装置
WO2019111323A1 (ja) * 2017-12-05 2019-06-13 三菱電機株式会社 折れ線作画装置及び折れ線作画方法
WO2020080088A1 (ja) * 2018-10-15 2020-04-23 三菱電機株式会社 情報処理装置
WO2020101424A1 (ko) * 2018-11-16 2020-05-22 한국과학기술원 지상 로봇의 움직임을 고려한 강인한 레이저 스캐너 매칭 방법 및 그 장치
CN112560747A (zh) * 2020-12-23 2021-03-26 苏州工业园区测绘地理信息有限公司 基于车载点云数据的车道边界交互式提取方法
JP7449506B2 (ja) 2020-08-06 2024-03-14 株式会社イクシス 計測システム
JP7452797B2 (ja) 2020-08-06 2024-03-19 株式会社イクシス 計測システム

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005202510A (ja) * 2004-01-13 2005-07-28 Asahi Koyo Kk 3次元デ−タ処理方法及び装置
JP2006003332A (ja) * 2004-06-16 2006-01-05 Kokusai Kogyo Co Ltd 地盤高推定方法
JP2011158278A (ja) * 2010-01-29 2011-08-18 Pasuko:Kk レーザデータのフィルタリング方法及び装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005202510A (ja) * 2004-01-13 2005-07-28 Asahi Koyo Kk 3次元デ−タ処理方法及び装置
JP2006003332A (ja) * 2004-06-16 2006-01-05 Kokusai Kogyo Co Ltd 地盤高推定方法
JP2011158278A (ja) * 2010-01-29 2011-08-18 Pasuko:Kk レーザデータのフィルタリング方法及び装置

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016142535A (ja) * 2015-01-29 2016-08-08 田中 成典 測量データ処理装置
WO2019111323A1 (ja) * 2017-12-05 2019-06-13 三菱電機株式会社 折れ線作画装置及び折れ線作画方法
JPWO2019111323A1 (ja) * 2017-12-05 2020-05-28 三菱電機株式会社 折れ線作画装置
WO2020080088A1 (ja) * 2018-10-15 2020-04-23 三菱電機株式会社 情報処理装置
JPWO2020080088A1 (ja) * 2018-10-15 2021-06-10 三菱電機株式会社 情報処理装置、情報処理方法及びプログラム
JP7046218B2 (ja) 2018-10-15 2022-04-01 三菱電機株式会社 情報処理装置、情報処理方法及びプログラム
JP2022089828A (ja) * 2018-10-15 2022-06-16 三菱電機株式会社 情報処理装置
JP7232946B2 (ja) 2018-10-15 2023-03-03 三菱電機株式会社 情報処理装置、情報処理方法及びプログラム
WO2020101424A1 (ko) * 2018-11-16 2020-05-22 한국과학기술원 지상 로봇의 움직임을 고려한 강인한 레이저 스캐너 매칭 방법 및 그 장치
JP7449506B2 (ja) 2020-08-06 2024-03-14 株式会社イクシス 計測システム
JP7452797B2 (ja) 2020-08-06 2024-03-19 株式会社イクシス 計測システム
CN112560747A (zh) * 2020-12-23 2021-03-26 苏州工业园区测绘地理信息有限公司 基于车载点云数据的车道边界交互式提取方法

Also Published As

Publication number Publication date
JP6013962B2 (ja) 2016-10-25

Similar Documents

Publication Publication Date Title
JP6013962B2 (ja) 地物表面検出装置、地物表面検出方法、及びプログラム
JP6097182B2 (ja) データ解析装置、データ解析方法、及びプログラム
JP5161936B2 (ja) データ解析装置、データ解析方法、及びプログラム
Chen et al. A multiresolution hierarchical classification algorithm for filtering airborne LiDAR data
JP5566353B2 (ja) データ解析装置、データ解析方法、及びプログラム
CN106408604A (zh) 一种点云数据的滤波方法及装置
JP6840050B2 (ja) 鉄塔の傾き検出方法および傾き検出装置
JP2014190962A (ja) データ解析装置、データ解析方法、及びプログラム
JP6006179B2 (ja) データ解析装置、データ解析方法、及びプログラム
Liu et al. High-resolution DEM generated from LiDAR data for water resource management
CN108427741A (zh) 一种基于大量高精度控制点的dem相对误差评价方法
JP6032678B2 (ja) データ解析装置、データ解析方法、及びプログラム
KR101236882B1 (ko) 항공 라이다 자료를 이용한 토석류 발생지역의 지형복원 방법
JP6514901B2 (ja) 測量データ処理装置
JP7144993B2 (ja) 三次元地図生成装置、及びプログラム
JP5991492B2 (ja) データ解析装置、データ解析方法、及びプログラム
JP6018959B2 (ja) データ解析装置、データ解析方法、及びプログラム
KR100741155B1 (ko) 레이저 펄스를 이용한 수목의 높이 데이터 추출장치 및 그방법
KR101139796B1 (ko) 라이다 측량자료를 이용한 수치지도 제작에서의 도로 주변의 가로수 추출 시스템 및 그 방법
JP6157990B2 (ja) データ解析装置
JP2008009723A (ja) 地盤抽出方法、装置及びプログラム
Gao et al. PyVF: A python program for extracting vertical features from LiDAR-DEMs
Makinen et al. Creating visually consistent DEMs over water bodies using laser scanning data
JP7333163B2 (ja) 地形状態検出装置、地形状態検出方法および地形状態検出プログラム
JP2022076188A (ja) 小領域生成装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20151007

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20160905

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: 20160913

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160923

R150 Certificate of patent or registration of utility model

Ref document number: 6013962

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

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