特許 7366227 地表点抽出装置、地表点抽出方法およびプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B1)

(11)【特許番号】

(24)【登録日】2023-10-12

(45)【発行日】2023-10-20

(54)【発明の名称】地表点抽出装置、地表点抽出方法およびプログラム

(51)【国際特許分類】

   G06T 7/00 20170101AFI20231013BHJP

【ＦＩ】

G06T7/00 640

【請求項の数】 8

(21)【出願番号】P 2022185766

(22)【出願日】2022-11-21

【審査請求日】2022-11-21

(73)【特許権者】

【識別番号】000135771

【氏名又は名称】株式会社パスコ

(74)【代理人】

【識別番号】100099759

【弁理士】

【氏名又は名称】青木 篤

(74)【代理人】

【識別番号】100123582

【弁理士】

【氏名又は名称】三橋 真二

(74)【代理人】

【識別番号】100114018

【弁理士】

【氏名又は名称】南山 知広

(74)【代理人】

【識別番号】100180806

【弁理士】

【氏名又は名称】三浦 剛

(74)【代理人】

【識別番号】100207778

【弁理士】

【氏名又は名称】阿形 直起

(72)【発明者】

【氏名】李 勇鶴

(72)【発明者】

【氏名】坂元 光輝

【審査官】佐藤 実

(56)【参考文献】

【文献】特開２０２１－１７９３６５（ＪＰ，Ａ）

【文献】特開２０１１－１５８２７８（ＪＰ，Ａ）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｇ０６Ｔ ７／００

Ｇ０９Ｂ ２９／００

(57)【特許請求の範囲】

【請求項1】

飛行体から対象領域を計測して得られた３次元点群に含まれる複数の構成点を取得するとともに、前記取得された複数の構成点を候補点に設定する取得手段と、
複数の段階のそれぞれにおいて、前記対象領域を、直前の段階の単位領域よりも小さい複数の単位領域に区分する区分手段と、
前記複数の候補点のそれぞれについて、当該候補点の高度の低さに応じた指標を算出し、前記複数の段階のそれぞれにおいて、前記指標に基づいて、前記区分された各単位領域に含まれる候補点のうちから各単位領域の代表点を設定する設定手段と、
前記複数の段階のそれぞれにおいて、各単位領域の代表点から近似曲面を算出する算出手段と、
前記複数の段階のうちの最終段階以外の段階において、前記取得された複数の構成点のうち、前記近似曲面に対する相対高度が所定範囲である構成点を候補点として抽出する候補点抽出手段と、
前記複数の段階のうちの最終段階において、前記取得された複数の構成点のうち、前記近似曲面に対する相対高度が所定範囲である構成点を地表を示す地表点として抽出する地表点抽出手段と、
を有することを特徴とする地表点抽出装置。

【請求項2】

前記設定手段は、前記複数の単位領域のうち、少なくとも一つの候補点を含む単位領域のみにおいて前記代表点を設定する、
請求項１に記載の地表点抽出装置。

【請求項3】

前記設定手段は、各単位領域に直前の段階において代表点に設定された候補点が含まれる場合に当該候補点を当該単位領域の代表点に設定するとともに、直前の段階において代表点に設定されていない候補点のうちから新たな代表点を設定する、
請求項１または２に記載の地表点抽出装置。

【請求項4】

前記算出手段は、
前記代表点のうち、直前の段階において代表点に設定されていた代表点の重みが、直前の段階において代表点に設定されていなかった代表点の重みよりも大きくなるようにして仮の近似曲面を算出し、
前記仮の近似曲面に対する相対高度が高い代表点ほど重みが小さくなるようにして前記近似曲面を算出する、
請求項３に記載の地表点抽出装置。

【請求項5】

前記設定手段は、各候補点と当該候補点を含む単位領域の中心位置との間の距離に基づく指標を算出する、
請求項１に記載の地表点抽出装置。

【請求項6】

前記算出手段は、各段階において同一の数の前記単位領域からなるグループごとに前記近似曲面を算出する、
請求項１に記載の地表点抽出装置。

【請求項7】

飛行体から対象領域を計測して得られた３次元点群に含まれる複数の構成点を取得するとともに、前記取得された複数の構成点を候補点に設定し、
複数の段階のそれぞれにおいて、前記対象領域を、直前の段階における単位領域よりも小さい複数の単位領域に区分し、
前記複数の候補点のそれぞれについて、当該候補点の高度の低さに応じた指標を算出し、前記複数の段階のそれぞれにおいて、前記指標に基づいて、前記区分された各単位領域に含まれる候補点のうちから各単位領域の代表点を設定し、
前記複数の段階のそれぞれにおいて、各単位領域の代表点から近似曲面を算出し、
前記複数の段階のうちの最終段階以外の段階において、前記取得された複数の構成点のうち、前記近似曲面に対する相対高度が所定範囲である構成点を候補点として抽出し、
前記複数の段階のうちの最終段階において、前記取得された複数の構成点のうち、前記近似曲面に対する相対高度が所定範囲である構成点を地表を示す地表点として抽出する、
ことを含むことを特徴とする地表点抽出方法。

【請求項8】

飛行体から対象領域を計測して得られた３次元点群に含まれる複数の構成点を取得するとともに、前記取得された複数の構成点を候補点に設定し、
複数の段階のそれぞれにおいて、前記対象領域を、直前の段階における単位領域よりも小さい複数の単位領域に区分し、
前記複数の候補点のそれぞれについて、当該候補点の高度の低さに応じた指標を算出し、前記複数の段階のそれぞれにおいて、前記指標に基づいて、前記区分された各単位領域に含まれる候補点のうちから各単位領域の代表点を設定し、
前記複数の段階のそれぞれにおいて、各単位領域の代表点から近似曲面を算出し、
前記複数の段階のうちの最終段階以外の段階において、前記取得された複数の構成点のうち、前記近似曲面に対する相対高度が所定範囲である構成点を候補点として抽出し、
前記複数の段階のうちの最終段階において、前記取得された複数の構成点のうち、前記近似曲面に対する相対高度が所定範囲である構成点を地表を示す地表点として抽出する、
ことをコンピュータに実行させることを特徴とするプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、地表点抽出装置、地表点抽出方法およびプログラムに関する。

【背景技術】

【0002】

ＤＥＭ（Digital Elevation Model）の生成等を目的として航空レーザ測量等により得られた３次元点群から地表点を抽出する、フィルタリングと呼ばれる技術が知られている。

【0003】

非特許文献１には、ＭＣＣ（Multiscale Curvature Classification）と呼ばれるフィルタリング技術が記載されている。ＭＣＣでは、処理対象の地理的領域に設定されるウィンドウ内の３次元点群の近似曲面が算出され、近似曲面よりも高度が一定以上高い点が非地表点として除外される。ウィンドウを段階的に大きくしながら近似曲面の算出および非地表点の除外が繰り返され、最後まで除外されずに残った点が地表点として抽出される。

【先行技術文献】

【非特許文献】

【0004】

【文献】A Multiscale Curvature Algorithm for Classifying Discrete Return LiDAR in Forested Environments, S. Evans and T. Hudak, IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, VOL. 45, NO. 4，2007

【発明の概要】

【発明が解決しようとする課題】

【0005】

ＭＣＣでは、ウィンドウ内に非地表点が密集している場合に、その非地表点に基づいて近似曲面が算出されるため、非地表点が地表点として誤抽出（過剰抽出）される場合がある。

【0006】

本発明は、上述の課題を解決するためになされたものであり、３次元点群から地表点を高精度に抽出することを可能とする地表点抽出装置、地表点抽出方法およびプログラムを提供することを目的とする。

【課題を解決するための手段】

【0007】

本発明の実施形態に係る地表点抽出装置は、飛行体から対象領域を計測して得られた３次元点群に含まれる複数の構成点を取得するとともに、取得された複数の構成点を候補点に設定する取得手段と、複数の段階のそれぞれにおいて、対象領域を、直前の段階の単位領域よりも小さい複数の単位領域に区分する区分手段と、複数の候補点のそれぞれについて、候補点の高度の低さに応じた指標を算出し、複数の段階のそれぞれにおいて、指標に基づいて、区分された各単位領域に含まれる候補点のうちから各単位領域の代表点を設定する設定手段と、複数の段階のそれぞれにおいて、各単位領域の代表点から近似曲面を算出する算出手段と、複数の段階のうちの最終段階以外の段階において、取得された複数の構成点のうち、近似曲面に対する相対高度が所定範囲である構成点を候補点として抽出する候補点抽出手段と、複数の段階のうちの最終段階において、取得された複数の構成点のうち、近似曲面に対する相対高度が所定範囲である構成点を地表を示す地表点として抽出する地表点抽出手段と、を有することを特徴とする。

【0008】

また、設定手段は、複数の単位領域のうち、少なくとも一つの候補点を含む単位領域のみにおいて代表点を設定することが好ましい。

【0009】

また、設定手段は、各単位領域に直前の段階において代表点に設定された候補点が含まれる場合に候補点を単位領域の代表点に設定するとともに、直前の段階において代表点に設定されていない候補点のうちから新たな代表点を設定する、ことが好ましい。

【0010】

また、算出手段は、代表点のうち、直前の段階において代表点に設定されていた代表点の重みが、直前の段階において代表点に設定されていなかった代表点の重みよりも大きくなるようにして仮の近似曲面を算出し、仮の近似曲面に対する相対高度が高い代表点ほど重みが小さくなるようにして近似曲面を算出することが好ましい。

【0011】

また、設定手段は、各候補点と当該候補点を含む単位領域の中心位置との間の距離に基づく指標を算出することが好ましい。

【0012】

また、算出手段は、各段階において同一の数の前記単位領域からなるグループごとに近似曲面を算出することが好ましい。

【0013】

本発明の実施形態に係る地表点抽出方法は、飛行体から対象領域を計測して得られた３次元点群に含まれる複数の構成点を取得するとともに、取得された複数の構成点を候補点に設定し、複数の段階のそれぞれにおいて、対象領域を、直前の段階における単位領域よりも小さい複数の単位領域に区分し、複数の候補点のそれぞれについて、候補点の高度の低さに応じた指標を算出し、複数の段階のそれぞれにおいて、指標に基づいて、区分された各単位領域に含まれる候補点のうちから各単位領域の代表点を設定し、複数の段階のそれぞれにおいて、各単位領域の代表点から近似曲面を算出し、複数の段階のうちの最終段階以外の段階において、取得された複数の構成点のうち、近似曲面に対する相対高度が所定範囲である構成点を候補点として抽出し、複数の段階のうちの最終段階において、取得された複数の構成点のうち、近似曲面に対する相対高度が所定範囲である構成点を地表を示す地表点として抽出する、ことを含むことを特徴とする。

【0014】

本発明の実施形態に係るプログラムは、飛行体から対象領域を計測して得られた３次元点群に含まれる複数の構成点を取得するとともに、取得された複数の構成点を候補点に設定し、複数の段階のそれぞれにおいて、対象領域を、直前の段階における単位領域よりも小さい複数の単位領域に区分し、複数の候補点のそれぞれについて、候補点の高度の低さに応じた指標を算出し、複数の段階のそれぞれにおいて、指標に基づいて、区分された各単位領域に含まれる候補点のうちから各単位領域の代表点を設定し、複数の段階のそれぞれにおいて、各単位領域の代表点から近似曲面を算出し、複数の段階のうちの最終段階以外の段階において、取得された複数の構成点のうち、近似曲面に対する相対高度が所定範囲である構成点を候補点として抽出し、複数の段階のうちの最終段階において、取得された複数の構成点のうち、近似曲面に対する相対高度が所定範囲である構成点を地表を示す地表点として抽出する、ことをコンピュータに実行させることを特徴とする。

【発明の効果】

【0015】

本発明に係る地表点抽出装置、地表点抽出方法およびプログラムは、３次元点群から地表点を高精度に抽出することを可能とする。

【図面の簡単な説明】

【0016】

【図1】地表点抽出装置１の機能ブロック図である。

【図2】構成点テーブルＴ１のデータ構造の例を示す図である。

【図3】地表点抽出処理の流れの例を示すフロー図である。

【図4】単位領域およびグループの設定について説明するための模式図である。

【図5】代表点設定処理の流れの例を示すフロー図である。

【図6】指標の算出および代表点の設定について説明するための模式図である。

【図7】検出処理の流れの例を示すフロー図である。

【図8】近似平面の傾斜角度に関する条件について説明するための模式図である。

【図9】近似曲面算出処理の流れの例を示すフロー図である。

【図10】（Ａ）および（Ｂ）は、近似曲面の算出について説明するための模式図である。

【発明を実施するための形態】

【0017】

以下、図面を参照しつつ、本発明の様々な実施形態について説明する。本発明の技術的範囲はそれらの実施形態に限定されず、特許請求の範囲に記載された発明及びその均等物に及ぶ点に留意されたい。

【0018】

図１は、本発明の実施形態に係る地表点抽出装置１の機能ブロック図である。地表点抽出装置１は、飛行体から対象領域を計測して得られた３次元点群を構成する点のうちから地表を示す地表点を抽出する。３次元点群は、飛行体から対象領域の地表面に向けてレーザを照射したときの反射点の位置を算出することにより得られる。レーザは対象領域の樹木や建造物のような地物でも反射されるため、３次元点群には、地表面を示す地表点と地表面以外の地物を示す非地表点とが含まれる。すなわち、地表点抽出装置１は、地表点と非地表点とを含む３次元点群から地表点のみを抽出する。地表点抽出装置１は、記憶部１１、通信部１２、表示部１３、操作部１４および処理部１５を有する。

【0019】

記憶部１１は、プログラムおよびデータを記憶するための構成であり、例えば半導体メモリを備える。記憶部１１は、プログラムとして、処理部１５による処理に用いられるオペレーティングシステムプログラム、ドライバプログラム、アプリケーションプログラム等を記憶する。プログラムは、ＣＤ－ＲＯＭ（Compact Disc Read Only Memory）、ＤＶＤ－ＲＯＭ（Digital Versatile Disc Read Only Memory）等のコンピュータ読取可能かつ非一時的な可搬型記憶媒体から記憶部１１にインストールされる。

【0020】

通信部１２は、地表点抽出装置１を他の装置と通信可能にする構成であり、通信インタフェース回路を備える。通信部１２が備える通信インタフェース回路は、有線ＬＡＮ（Local Area Network）又は無線ＬＡＮ等の通信インタフェース回路である。通信部１２は、データを他の装置から受信して処理部１５に供給するとともに、処理部１５から供給されたデータを他の装置に送信する。

【0021】

表示部１３は、画像を表示するための構成であり、例えば液晶ディスプレイまたは有機ＥＬ（Electro-Luminescence）ディスプレイを備える。表示部１３は、処理部１５から供給された表示データに基づいて画像を表示する。

【0022】

操作部１４は、地表点抽出装置１に対する操作を受け付けるための構成であり、例えばキーボード、キーパッド、マウス等を備える。操作部１４は、表示部１３と一体化されたタッチパネルでもよい。操作部１４は、受け付けた操作に応じた操作信号を生成して処理部１５に供給する。

【0023】

処理部１５は、地表点抽出装置１の動作を統括的に制御する構成であり、一つまたは複数のプロセッサおよびその周辺回路を備える。処理部１５は、例えば、ＣＰＵ（Central Processing Unit）を備える。処理部１５は、ＧＰＵ（Graphics Processing Unit）、ＤＳＰ（Digital Signal Processor）、ＬＳＩ（Large Scale Integration）、ＡＳＩＣ（Application Specific Integrated Circuit）、ＦＰＧＡ（Field Programmable Gate Array）等を備えてもよい。処理部１５は、記憶部１１に記憶されているプログラムに基づいて地表点抽出装置１の各種処理が適切な手順で実行されるように、各構成の動作を制御するとともに、各種の処理を実行する。

【0024】

処理部１５は、取得部１５１、区分部１５２、選択部１５３、設定部１５４、検出部１５５、算出部１５６、候補点抽出部１５７、地表点抽出部１５８および出力部１５９を機能ブロックとして備える。これらの各部は、処理部１５が実行するプログラムによって実現される機能モジュールである。これらの各部は、ファームウェアとして地表点抽出装置１に実装されてもよい。取得部１５１、区分部１５２、選択部１５３、設定部１５４、検出部１５５、算出部１５６、候補点抽出部１５７、地表点抽出部１５８および出力部１５９は、それぞれ取得手段、区分手段、選択手段、設定手段、検出手段、算出手段、候補点抽出手段、地表点抽出手段および出力手段の一例である。

【0025】

図２は、記憶部１１に記憶される構成点テーブルＴ１のデータ構造の例を示す図である。構成点は、３次元点群を構成する点である。構成点テーブルＴ１は、構成点ＩＤ、座標およびパルス種別を相互に関連付けて記憶する。

【0026】

構成点ＩＤは、３次元点群を構成する構成点のそれぞれを識別する情報である。座標は、構成点の位置を示す３次元座標であり、例えば平面直角座標系で規定される。なお、以降では、３次元座標のうち、ｘ座標またはｙ座標によって規定される位置を水平位置と称し、ｚ座標によって規定される位置を高度と称することがある。パルス種別として、ファーストパルス、中間パルスおよびラストパルスのうちのいずれか１つが記憶される。

【0027】

構成点テーブルＴ１の各データは、飛行体から対象領域を計測して得られた計測データに基づいてあらかじめ生成される。飛行体には、ＧＮＳＳ（Global Navigation Satellite System）、ＩＭＵ（Inertial Measurement Unit）およびレーザ計測器が搭載される。レーザ計測器は、対象領域に向けてレーザを照射してその反射光を検出することにより反射点までの距離を計測する。ＧＮＳＳおよびＩＭＵによる計測結果に基づく飛行体の位置および姿勢と、レーザ計測器によって計測された反射点までの距離とに基づいて、反射点である構成点の位置が算出される。

【0028】

レーザ計測器から照射されたレーザが複数の反射点で反射され、複数の反射光が検出されることがある。この場合、最初に検出された反射光がファーストパルスであり、最後に検出された反射光がラストパルスであり、ファーストパルスとラストパルスとの間に検出された反射光が中間パルスである。一般に、ファーストパルスおよび中間パルスの反射点はラストパルスの反射点よりも上方であるため、ファーストパルス、中間パルスおよびラストパルスのうち、地表点で反射した可能性が最も高いのはラストパルスである。なお、一回のレーザの照射に対して１つの反射光のみが検出された場合、その反射光はラストパルスであるものとする。

【0029】

パルス種別に代えて、一回のレーザの照射に対して検出された複数の反射光のうち何番目の反射光であるかを示す反射順が記憶されてもよい。この場合、反射順の値が１である反射光がファーストパルスであり、反射順の値が最大である反射光がラストパルスであり、他の反射光が中間パルスである。

【0030】

図３は、地表点抽出装置１によって実行される地表点抽出処理の流れの例を示すフロー図である。地表点抽出処理は、記憶部１１に記憶されたプログラムに基づいて、処理部１５が地表点抽出装置１の他の構成と協働することにより実現される。

【0031】

以下では、地表点抽出処理の流れの概略について説明する。

【0032】

最初に、対象領域が複数の単位領域に区分される。また、対象領域に含まれる構成点の少なくとも一部が候補点に設定される。例えば、構成点テーブルＴ１に記憶された構成点のうちラストパルスに対応する構成点が候補点に設定される。次に、候補点の高度に基づいて、地表点である可能性が高い候補点が各単位領域の代表点に設定される。次に、代表点の高度に基づいて、地表面の傾斜が急変する傾斜変化点を含む単位領域が検出される。次に、傾斜変化点を含む単位領域の代表点の重みが他の代表点の重みよりも大きくなるようにして、代表点の近似曲面が算出される。次に、近似曲面に対する相対高度が所定範囲である構成点が候補点として抽出される。所定範囲は、例えば、所定の閾値未満の範囲である。相対高度が高い候補点は、地表面よりも上方の樹木や建造物で反射された反射点に対応し、非地表点である可能性が高いため、相対高度が低い構成点が候補点として抽出される。処理対象の領域を段階的に縮小するとともに、相対高度の所定範囲の閾値を減少させながら上述の処理が繰り返される。構成点テーブルＴ１に記憶された構成点のうち、最終段階において近似曲面に対する相対高度が所定範囲である構成点が地表点として抽出される。

【0033】

単位領域の候補点の多くが非地表点である場合には、非地表点が代表点に設定される可能性が高くなり、地表面の形状が適切に反映された近似曲面を算出することができない場合がある。しかし、単位領域を最初に広く設定することにより、最初の段階における単位領域のそれぞれに地表点が含まれる可能性が高くなる。そして、最初の段階において非地表点を候補点から除外することにより、後段階では地表点である候補点の割合が多くなり、地表面の形状が適切に反映された近似曲面が算出される。したがって、単位領域が段階的に縮小されることにより、非地表点が密集している領域を含む３次元点群から地表点を高精度に抽出することが可能となる。

【0034】

また、一般に、近似曲面は傾斜が連続的に変化する面として算出されるため、地表面の傾斜が急変する領域においては地表面の形状を適切に反映することができない場合がある。しかし、傾斜変化点を含む単位領域の代表点に大きい重みを設定することにより、地表面の傾斜が急変する領域においても地表面の形状が近似曲面に適切に反映されるようになる。したがって、傾斜変化点を含む単位領域の代表点の重みが大きく設定されることにより、傾斜が急変する領域を含む３次元点群から地表点を高精度に抽出することが可能となる。

【0035】

以下では、地表点抽出処理の流れの詳細について説明する。

【0036】

最初に、取得部１５１は、地表点を抽出する対象となる対象領域および複数の構成点を取得するとともに、複数の構成点を候補点に設定する（ステップＳ１０１）。例えば、取得部１５１は、操作部１４に対する利用者の操作によって指定された地理的領域を対象領域として取得する。対象領域は任意の形状の領域であってよいが、以下では矩形領域であるものとする。また、取得部１５１は、構成点テーブルＴ１に記憶された構成点のうち、対象領域に含まれる構成点を取得する。対象領域に含まれる構成点は、ｘ座標およびｙ座標によって規定される候補点の水平位置が対象領域に含まれる構成点である。

【0037】

取得部１５１は、操作部１４に対する利用者の操作によって指定された計測データに基づく全ての構成点の水平位置を包含する地理的領域を対象領域として取得してもよい。その場合、対象領域は、後述する単位領域の辺の長さの倍数である長さの辺を有する矩形領域であるのが好ましい。また、取得部１５１は、通信部１２を介して他の装置から対象領域および構成点を取得してもよい。

【0038】

取得部１５１は、複数の構成点のうちの少なくとも一部の構成点を候補点に設定する。例えば、取得部１５１は、対象領域に含まれる構成点のうち、ラストパルスに対応する構成点を候補点に設定する。

【0039】

次に、区分部１５２は、対象領域を複数の単位領域に区分する（ステップＳ１０２）。単位領域は任意の形状の領域であってよいが、以下では矩形領域であるものとする。また、区分部１５２は、相互に隣接する所定個数の単位領域からなる複数のグループを設定する。

【0040】

図４（Ａ）および（Ｂ）は、単位領域およびグループの設定について説明するための模式図である。図４（Ａ）は、最初の段階における単位領域およびグループについて示している。図４（Ａ）に示す例では、正方形である対象領域Ｔが１６個（４×４個）の単位領域Ｕ１に区分されており、各単位領域の境界が破線により図示されている。また、相互に隣接する４個（２×２個）の単位領域Ｕ１からなるグループＧ１が設定されており、各グループの境界が実線により図示されている。

【0041】

ステップＳ１０２－Ｓ１０９の処理は、単位領域の大きさを段階的に縮小させながら繰り返し実行される。図４（Ｂ）は、図４（Ａ）の次の段階における単位領域およびグループについて示している。図４（Ｂ）に示す例では、２段階目の単位領域Ｕ２（２回目に実行されるステップＳ１０２において設定される単位領域）は１段階目の単位領域Ｕ１（１回目に実行されるステップＳ１０２において設定される単位領域）の半分の大きさ（４分の１の面積）に設定されている。これにより、対象領域Ｔは６４個（８×８個）の単位領域Ｕ２に区分されている。すなわち、対象領域Ｔは、二段階目以降において、直前の段階における単位領域Ｕ１よりも小さい複数の単位領域Ｕ２に区分される。また、相互に隣接する４個（２×２個）の単位領域Ｕ２からなるグループＧ２が設定されている。このように、単位領域は後段階の単位領域が前段階の単位領域よりも小さくなるように設定され、グループは各段階において同一の数の単位領域からなるように設定される。したがって、グループが段階的に縮小される。

【0042】

段階数、各段階における単位領域の大きさ、グループを構成する単位領域の数は、任意に設定されてよい。段階数は、例えば６段階に設定される。この場合、各段階における単位領域は、例えば四辺がそれぞれ３２ｍ、１６ｍ、８ｍ、４ｍ、２ｍ、１ｍの矩形領域に設定される。また、グループを構成する単位領域の数は、例えば１００個（１０×１０個）に設定される。このような例に限られず、段階数は、非地表点が地表点として誤抽出（過剰抽出）される可能性が小さくなるように事前の実験を通じて適宜に設定されてよい。単位領域の大きさ、特に最初の段階における単位領域の大きさは、非地表点が密集し得る建造物の大きさや樹木の密集部の広さを考慮して適宜に設定されてよい。グループを構成する単位領域の数は、近似曲面の精度と近似曲面の算出処理の負荷とのトレードオフを考慮して適宜に設定されてよい。

【0043】

なお、図４（Ａ）および（Ｂ）に示す例では、２段階目の単位領域Ｕ２は、１段階目の単位領域Ｕ１の半分の大きさを有し、１段階目の単位領域Ｕ１を４つに区分したものとなっている。すなわち、後段階の単位領域は前段階の単位領域を区分することにより設定されている。このような例に限られず、後段階の単位領域は前段階の単位領域にかかわらずに設定されてよい。例えば、単位領域が直前の単位領域の０．９倍の大きさに設定される場合には、後段階の単位領域は前段階の単位領域を区分したものではなくなる。

【0044】

図３に戻り、次に、選択部１５３は、処理対象となるグループを選択する（ステップＳ１０３）。選択部１５３は、ステップＳ１０２で設定された複数のグループのうち、まだ処理対象として選択されていないグループを処理対象のグループとして選択する。

【0045】

次に、設定部１５４は、処理対象のグループの各単位領域に含まれる候補点のうちから代表点を設定する代表点設定処理を実行する（ステップＳ１０４）。代表点設定処理の詳細は後述する。

【0046】

次に、検出部１５５は、傾斜変化点を含む単位領域を検出する検出処理を実行する（ステップＳ１０５）。検出処理の詳細は後述する。

【0047】

次に、算出部１５６は、各単位領域の代表点から近似曲面を算出する近似曲面算出処理を実行する（ステップＳ１０６）。近似曲面算出処理の詳細は後述する。

【0048】

次に、選択部１５３は、ステップＳ１０２で設定された全てのグループが処理対象のグループとして選択されたか否かを判定する（ステップＳ１０７）。

【0049】

全てのグループが処理対象のグループとして選択されていない場合（ステップＳ１０７－Ｎｏ）、地表点抽出処理はステップＳ１０３に戻り、選択部１５３は、まだ処理対象として選択されていないグループを処理対象のグループとして選択する。

【0050】

全てのグループが処理対象のグループとして選択された場合（ステップＳ１０７－Ｙｅｓ）、候補点抽出部１５７は、最終段階の処理が実行されたか否かを判定する（ステップＳ１０８）。候補点抽出部１５７は、ステップＳ１０２－Ｓ１０７の処理が予め設定された段階数だけ実行された場合に、最終段階の処理が実行されたと判定する。

【0051】

まだ最終段階の処理が実行されていない場合（ステップＳ１０８－Ｎｏ）、候補点抽出部１５７は、ステップＳ１０１で取得された構成点のうち、現段階において算出された近似曲面に対する相対高度が所定範囲である構成点を候補点として抽出する（ステップＳ１０９）。相対高度は、構成点の高度から、構成点と同一の水平位置における近似曲面の高度を減じた値である。所定範囲は、例えば、所定の閾値未満の範囲である。この場合、閾値は、段階的に小さくなるように設定される。すなわち、閾値は、二段階目以降において、直前の段階の閾値よりも小さくなるように設定される。例えば、閾値は、直前の段階の閾値の半分となるように設定される。

【0052】

ステップＳ１０８の後に、地表点抽出処理はステップＳ１０２に戻り、区分部１５２は、対象領域を次段階の大きさの単位領域に区分する。また、区分部１５２は、相互に隣接する所定個数の単位領域からなる複数のグループを設定する。

【0053】

最終段階の処理が実行された場合（ステップＳ１０８－Ｙｅｓ）、地表点抽出部１５８は、ステップＳ１０１で取得された複数の構成点のうち、最終段階において各グループについて算出された近似曲面に対する相対高度が所定範囲である構成点を地表点として抽出する（ステップＳ１１０）。所定範囲は、例えば、所定の閾値未満の範囲である。閾値は、直前の段階のステップＳ１０９における閾値よりも小さくなるように設定されてもよく、ステップＳ１０９における閾値にかかわらずあらかじめ設定されてもよい。

【0054】

次に、出力部１５９は、抽出された地表点を出力する（ステップＳ１１１）。例えば、出力部１５９は、通信部１２を介して、抽出された地表点のデータを他の装置に送信する。出力部１５９は、抽出された地表点を表示部１３に表示してもよい。以上で、地表点抽出処理が終了する。

【0055】

図５は、代表点設定処理の流れの例を示すフロー図である。代表点設定処理は、処理対象のグループの各単位領域に含まれる候補点のうちから代表点を設定する処理である。代表点設定処理は、地表点抽出処理のステップＳ１０４において実行される。

【0056】

最初に、設定部１５４は、代表点を設定する対象となる処理対象の単位領域を選択する（ステップＳ２０１）。設定部１５４は、処理対象のグループの単位領域のうち、まだ処理対象として選択されていない単位領域を処理対象の単位領域として選択する。

【0057】

次に、設定部１５４は、処理対象の単位領域に候補点が含まれるか否かを判定する（ステップＳ２０２）。設定部１５４は、処理対象の単位領域に構成点テーブルＴ１に記憶されているいずれかの候補点の水平位置が含まれるか否かにより、処理対象の単位領域に候補点が含まれるか否かを判定する。

【0058】

処理対象の単位領域の候補点が含まれない場合（ステップＳ２０２－Ｎｏ）、設定部１５４は、処理対象の単位領域に代表点を設定しない（ステップＳ２０３）。

【0059】

処理対象の単位領域に候補点が含まれる場合（ステップＳ２０２－Ｙｅｓ）、設定部１５４は、処理対象の単位領域の候補点に、直前の段階において代表点に設定された候補点が含まれるか否かを判定する（ステップＳ２０４）。現段階が最初の段階である場合には、設定部１５４は、処理対象の単位領域の候補点に、直前の段階において代表点に設定された候補点は含まれないと判定する。

【0060】

処理対象の単位領域の候補点に直前段階に代表点として設定された候補点が含まれる場合（ステップＳ２０４－Ｙｅｓ）、設定部１５４は、直前段階の代表点であった候補点を代表点に設定する（ステップＳ２０５）。すなわち、直前の段階において代表点に設定された候補点は、現段階でも引き続き代表点に設定される。

【0061】

処理対象の単位領域の候補点に直前段階に代表点として設定された候補点が含まれない場合（ステップＳ２０４－Ｎｏ）、代表点設定処理はステップＳ２０６に進む。

【0062】

次に、設定部１５４は、処理対象の単位領域の候補点のそれぞれについて、候補点の高度の低さに応じた指標を算出する（ステップＳ２０６）。例えば、設定部１５４は、候補点と候補点を含む単位領域の中心位置との間の距離に基づく指標を算出する。指標は、候補点の高度が低いほど小さくなり、候補点と単位領域の中心位置との間の距離が短いほど小さくなる。例えば、指標Ｉは、候補点の高度と、候補点と単位領域の中心位置との間の水平距離との重み付け和であり、次の式に従う。

【数1】

ここで、Ｈは候補点の高度であり、Ｄは候補点の位置と単位領域の中心位置との間の水平距離である。また、ａは重みであり、事前の実験を通じて０より大きい任意の値に設定される。

【0063】

次に、設定部１５４は、指標に基づいて、処理対象の単位領域に含まれる、代表点に設定されていない候補点のうちのいずれかを処理対象の単位領域の新たな代表点に設定する（ステップＳ２０７）。例えば、設定部１５４は、代表点に設定されていない候補点のうち最小の指標を有する候補点を代表点に設定する。すなわち、直前の段階に代表点に設定された候補点が処理対象の単位領域に含まれる場合には、直前の段階に代表点に設定された候補点と、直前の段階で代表点に設定されていない候補点のうち最小の指標を有する候補点とが代表点に設定される。

【0064】

図６は、指標の算出および代表点の設定について説明するための模式図である。上述したように、代表点は、候補点の高度Ｈの低さに応じた指標が最小となる候補点である。非地表点である候補点は地物で反射されたレーザの反射点であるため、地表面よりも高い高度を有する。代表点が候補点の高度Ｈの低さに応じた指標が最小となる候補点に設定されることにより、地表点である可能性が高い候補点が代表点に設定される。

【0065】

また、代表点は、各候補点と当該候補点を含む単位領域の中心位置との間の距離に基づく指標が最小となる候補点である。すなわち、候補点が処理対象の単位領域の中心に近いほど、その候補点が代表点に設定される可能性が高くなる。代表点が候補点の水平位置に基づかずに設定される場合、図６の領域Ａのように、単位領域の境界付近に地表面が窪んだ領域が存在するときに、多くの代表点がその領域に含まれる候補点に設定される可能性がある。この場合、地表面の形状が代表点に対する近似曲面に適切に反映されなくなる。処理対象の単位領域の中心に近い候補点を代表点に設定されやすくすることにより、単位領域の境界付近の候補点が代表点に設定されにくくなる。したがって、代表点が分散され、地表面の形状が適切に反映された近似曲面が算出される可能性が高くなる。

【0066】

図５に戻り、ステップＳ２０３またはＳ２０７の次に、設定部１５４は、処理対象のグループの全ての単位領域が処理対象の単位領域として選択されたか否かを判定する（ステップＳ２０８）。

【0067】

全ての単位領域が処理対象の単位領域として選択されていない場合（ステップＳ２０７－Ｎｏ）、代表点設定処理はステップＳ２０１に戻り、設定部１５４は、まだ処理対象として選択されていない単位領域を処理対象の単位領域として選択する。

【0068】

全ての単位領域が処理対象の単位領域として選択された場合（ステップＳ２０７－Ｙｅｓ）、代表点設定処理が終了する。

【0069】

図７は、検出処理の流れを示すフロー図である。検出処理は、各単位領域が傾斜変化点を含むか否かを判定することにより傾斜変化点を含む単位領域を検出する処理である。検出処理は、地表点抽出処理のステップＳ１０５で実行される。

【0070】

最初に、検出部１５５は、傾斜変化点を含むか否かの判定の対象となる処理対象の単位領域を選択する（ステップＳ３０１）。検出部１５５は、処理対象のグループの単位領域のうち、まだ処理対象として選択されていない単位領域を処理対象の単位領域として選択する。

【0071】

次に、検出部１５５は、処理対象の単位領域の代表点と処理対象の単位領域の周囲の単位領域の代表点とに対する二つの近似平面を算出する（ステップＳ３０２）。例えば、検出部１５５は、ＲＡＮＳＡＣ（Random Sample Consensus）法により、境界が処理対象の単位領域に含まれるような二つの近似平面を算出する。

【0072】

例えば、検出部１５５は、処理対象の単位領域の周囲の８個の単位領域の代表点を取得する。検出部１５５は、取得した８個の代表点のうちから２個の代表点を無作為に抽出し、処理対象の単位領域の代表点と、抽出された２個の代表点とを含む第１の平面を算出する。また、検出部１５５は、同様にして第１の平面と異なる第２の平面を算出する。検出部１５５は、各代表点と第１の平面および第２の平面との間の残差を算出する。残差は、例えば、各代表点から第１の平面までの距離および各代表点から第２の平面までの距離のうちの小さいものの二乗平均平方根である。検出部１５５は、上述の処理を繰り返して算出された第１の平面と第２の平面との組合せのうち残差が最も小さいものを二つの近似平面とする。

【0073】

次に、検出部１５５は、二つの近似平面の傾斜角度に関する第１の条件が満たされるか否かを判定する（ステップＳ３０３）。第１の条件は、二つの近似平面がなす角の角度に関する条件である。

【0074】

第１の条件が満たされる場合（ステップＳ３０３－Ｙｅｓ）、検出部１５５は、処理対象の単位領域が傾斜変化点を含むと判定する（ステップＳ３０４）。

【0075】

第１の条件が満たされない場合（ステップＳ３０３－Ｎｏ）、二つの近似平面の傾斜角度に関する第２の条件が満たされるか否かを判定する（ステップＳ３０５）。第２の条件は、二つの近似平面のそれぞれの傾斜角度に関する条件である。近似平面の傾斜角度は、近似平面と水平面がなす角の角度である。

【0076】

第２の条件が満たされる場合（ステップＳ３０５－Ｙｅｓ）、検出部１５５は、処理対象の単位領域が傾斜変化点を含むと判定する（ステップＳ３０４）。

【0077】

第２の条件が満たされない場合（ステップＳ３０５－Ｎｏ）、検出部１５５は、処理対象の単位領域が傾斜変化点を含まないと判定する（ステップＳ３０６）。

【0078】

図８は、近似平面の傾斜角度に関する第１の条件および第２の条件について説明するための模式図である。図８に示す例では、処理対象の単位領域の代表点と処理対象の単位領域の周囲の単位領域の代表点とが二つの近似平面Ｐ１およびＰ２により近似されている。また、図８では、水平面Ｈの方向が図示されている。

【0079】

第１の条件は、例えば、二つの近似平面Ｐ１およびＰ２がなす２つの角の角度のうち小さいもの（図８に示す例では、角度α）が所定角度以上であることである。第２の条件は、例えば、二つの近似平面のうちの一方の近似平面Ｐ１の傾斜角度（図８に示す例では、角度β）が０度を含む所定角度範囲内であり他方の近似平面Ｐ２の傾斜角度（図８に示す例では、角度γ）が所定傾斜角度以上であることである。所定角度範囲は、例えば－１０度以上＋１０度以下の範囲であり、所定傾斜角度は、例えば４０度である。検出部１５５は、第１の条件および第２の条件のうちの少なくとも一つが満たされる場合に、処理対象の単位領域が傾斜変化点を含むと判定する。

【0080】

なお、検出部１５５は、第１の条件および第２の条件の両方が満たされる場合に処理対象の単位領域が傾斜変化点を含むと判定してもよい。また、検出部１５５は、第１の条件および第２の条件のうちのいずれか一方のみに基づいて処理対象の単位領域が傾斜変化点を含むか否かを判定してもよい。

【0081】

図７に戻り、ステップＳ３０４またはＳ３０６の後に、検出部１５５は、処理対象のグループの全ての単位領域が処理対象の単位領域として選択されたか否かを判定する（ステップＳ３０６）。

【0082】

全ての単位領域が処理対象として選択されていない場合（ステップＳ３０６－Ｎｏ）、推定処理はステップＳ３０１に戻り、検出部１５５は、まだ処理対象として選択されていない単位領域を処理対象の単位領域として選択する。

【0083】

全ての単位領域が処理対象として選択された場合（ステップＳ３０６－Ｙｅｓ）、推定処理が終了する。

【0084】

図９は、近似曲面算出処理の流れの例を示すフロー図である。近似曲面算出処理は、代表点に対する近似曲面を算出する処理である。近似曲面算出処理は、地表点抽出処理のステップＳ１０６において実行される。

【0085】

最初に、算出部１５６は、各単位領域の代表点に重みを設定する（ステップＳ４０１）。算出部１５６は、傾斜変化点を含むと推定された単位領域の代表点に、他の単位領域の代表点よりも大きい重みを設定する。また、算出部１５６は、傾斜変化点を含まないと推定された単位領域の代表点のうち、直前の段階において代表点に設定された代表点に、直前の段階において代表点に設定されていない代表点よりも大きい重みを設定する。例えば、算出部１５６は、傾斜変化点を含むと推定された単位領域の代表点に重み８．０を設定し、直前の段階において代表点に設定された代表点に重み２．０を設定し、他の代表点に重み１．０を設定する。

【0086】

次に、算出部１５６は、設定された重みに基づいて、代表点から仮の近似曲面を算出する（ステップＳ４０２）。すなわち、算出部１５６は、直前の段階において代表点に設定されていた代表点の重みが、直前の段階において代表点に設定されていなかった代表点の重みよりも大きくなり、かつ傾斜変化点を含む単位領域の代表点の重みが前記傾斜変化点を含まない単位領域の代表点の重みよりも大きくなるようにして仮の近似曲面を算出する。例えば、算出部１５６は、ＴＰＳ（Thin Plate Spline）法により仮の近似曲面を算出する。

【0087】

次に、算出部１５６は、各代表点と仮の近似曲面との相対高度の標準偏差を算出する（ステップＳ４０３）。算出部１５６は、次の式により標準偏差σを算出する。

【数2】

ここで、ｗｉはｉ番目の代表点の重みであり、ｒｉはｉ番目の代表点と仮の近似曲面との間の相対高度である。相対高度は、代表点の高度から、代表点と同一の水平位置における仮の近似曲面の高度を減じた値である。

【0088】

次に、算出部１５６は、相対高度に基づいて代表点の重みを更新する（ステップＳ４０４）。例えば、算出部１５６は、傾斜変化点を含まない単位領域の代表点のうち、直前の段階において代表点に設定されていない代表点の重みを、相対高度が低いほど大きくなるように更新する。例えば、算出部１５６は、相対高度ｒｉに基づいて、次の式に従い代表点の重みｗｉを更新する。

【数3】

【0089】

次に、算出部１５６は、更新された重みに基づいて、代表点から近似曲面を算出する（ステップＳ４０５）。すなわち、算出部１５６は、仮の近似曲面に対する相対高度が高い代表点ほど重みが小さくなるようにして近似曲面を算出する。例えば、算出部１５６は、更新された重みを用いて、ステップＳ４０２と同様のＴＰＳ法により近似曲面を算出する。以上で、近似曲面算出処理を終了する。

【0090】

図１０（Ａ）および（Ｂ）は、近似曲面の算出について説明するための模式図である。図１０（Ａ）は、ステップＳ４０２で算出された仮の近似曲面を模式的に示している。図１０（Ａ）は、黒塗りの丸で示される地表点である代表点よりも、白抜きの丸で示される非地表点である代表点が多い例を示している。この場合、仮の近似曲面は、地表点である代表点よりも非地表点である代表点に近くなるように算出される。すなわち、仮の近似曲面は、地表面ではなく地物の表面の形状を反映したものとなる。

【0091】

図１０（Ｂ）は、ステップＳ４０５で算出された近似曲面を模式的に示している。ステップＳ４０４において仮の近似曲面に対する代表点の相対高度が低いほど重みが大きくなるように代表点の重みが更新されるため、図１０（Ｂ）に示す例では、低い高度を有する地表点である代表点の重みが大きくなっている。したがって、近似曲面は、非地表点である代表点よりも地表点である代表点に近くなるように算出される。すなわち、近似曲面は地表面の形状を反映したものとなる。

【0092】

このように、近似曲面算出処理は、処理対象のグループの代表点から近似曲面を算出する。すなわち、近似曲面は、処理対象のグループの全ての構成点から算出されるのではなく、構成点の一部である代表点のみから算出される。これにより、近似曲面を算出する処理の負荷が低減される。また、近似曲面算出処理は、処理対象のグループのそれぞれについて実行される。すなわち、近似曲面は、対象領域に含まれる全ての代表点から算出されるのではなく、処理対象のグループに含まれる代表点のみから算出される。これにより、近似曲面を算出する処理の負荷がより低減される。

【0093】

地表点抽出装置１には、次に述べるような変形例が適用されてもよい。

【0094】

上述した説明では、地表点抽出処理のステップＳ１０１において、取得部１５１がラストパルスに対応する構成点を候補点に設定するものとしたが、このような例に限られない。取得部１５１は、取得した構成点の全てを候補点に設定してもよい。すなわち、取得部１５１は、ラストパルスに対応する構成点に加えて、ファーストパルスおよび中間パルスに対応する構成点を候補点に設定してもよい。このようにしても、地表点抽出装置１は、３次元点群から地表点を高精度に抽出することができる。この場合、ステップＳ１０９において、候補点抽出部１５７は、ファーストパルスに対応する候補点、中間パルスに対応する候補点および近似曲面に対する相対高度が閾値以上であるラストパルスに対応する候補点を除外してもよい。また、取得部１５１は、構成点テーブルＴ１に記憶された構成点のうちラストパルスに対応する構成点のみを取得し、ファーストパルスまたは中間パルスに対応する構成点を取得しなくてもよい。

【0095】

上述した説明では、地表点抽出処理のステップＳ１０２において、区分部１５２が各段階に共通する個数の単位領域からなるグループを設定するものとしたが、このような例に限られない。区分部１５２は、段階ごとに異なる個数の単位領域からなるグループを設定してもよい。このようにしても、地表点抽出装置１は地表面の傾斜が急変する領域を含む３次元点群から地表点を高精度に抽出することができる。また、このようにしても、グループが段階的に縮小するように設定されていれば、地表点抽出装置１は非地表点が密集している領域を含む３次元点群から地表点を高精度に抽出することができる。

【0096】

上述した説明では、地表点抽出処理のステップＳ１０７において、候補点抽出部１５７が近似曲面に対する相対高度が所定範囲である構成点として、相対高度が閾値未満である構成点を候補点として抽出するものとしたが、このような例に限られない。所定範囲は、第１閾値未満であり、かつ第２閾値より大きい範囲でもよい。すなわち、候補点抽出部１５７は、近似曲面に対する相対高度が所定範囲内である構成点として、相対高度が第１閾値未満でありかつ第２閾値より大きい構成点を候補点として抽出してもよい。例えば、第１閾値は正値であり第２閾値は負値である。この場合、所定範囲は±０ｍを含む範囲であり、例えば、－０．５ｍ以上＋０．３ｍ以下の範囲である。このようにすることで、近似曲面に対する相対高度が著しく低く、ノイズであると疑われる候補点が除外されるため、地表点抽出装置１は３次元点群から地表点をより高精度に抽出することができる。

【0097】

上述した説明では、地表点抽出処理のステップＳ１０７において、候補点抽出部１５７が、ステップＳ１０１で取得された構成点のうちから候補点を抽出するものとしたが、このような例に限られない。例えば、現段階における候補点のうち、近似曲面に対する相対高度が閾値未満である候補点または近似曲面に対する相対高度が所定範囲内である候補点を次段階における候補点として抽出してもよい。すなわち、候補点抽出部１５７は、複数の段階のうちの最終段階以外の段階において、相対高度が閾値以上である候補点を次段階の候補点から除外してもよい。

【0098】

上述した説明では、代表点設定処理のステップＳ２０５において、設定部１５４が直前の段階で代表点に設定されていた候補点を代表点に設定するものとしたが、設定部１５４は、そのような代表点を候補点に設定しなくてもよい。この場合、代表点設定処理のステップＳ２０３およびＳ２０４が省略される。また、設定部１５４は、直前の段階で代表点に設定されていた候補点が代表点に設定された場合には、指標に基づいて代表点を新たに設定しなくてもよい。

【0099】

上述した説明では、代表点設定処理のステップＳ２０６において、設定部１５４が候補点の高度が低いほど小さくなり、かつ候補点とその候補点を含む単位領域の中心との間の水平距離が短いほど小さくなる指標を算出するものとしたが、このような例に限られない。設定部１５４は、候補点の水平位置には基づかない指標を算出してもよい。例えば、設定部１５４は、候補点の高度をそのまま指標としてもよい。このようにしても、地表点抽出装置１は、３次元点群から地表点を高精度に抽出することができる。また、指標は、候補点の高度が高いほど大きくなるものでもよく、候補点と単位領域の中心位置との間の距離が長いほど大きくなるものでもよい。この場合、ステップＳ２０７において、設定部１５４は最大の指標を有する候補点を代表点に設定する。

【0100】

上述した説明では、代表点設定処理のステップＳ２０７において、設定部１５４が代表点に設定されていない候補点のうち最小の指標を有する候補点を処理対象の単位領域の代表点に設定するものとしたが、このような例に限られない。設定部１５４は、代表点に設定されていない候補点のうち指標が小さい複数の候補点を代表点に設定してもよい。

【0101】

上述した説明では、推定処理のステップＳ３０２において、検出部１５５は境界が処理対象の単位領域に含まれるような二つの近似平面を算出するものとしたが、このような例に限られない。検出部１５５は、境界が処理対象の単位領域に含まれるか否かを考慮せずに二つの近似平面を算出してもよい。例えば、検出部１５５は、最小二乗法を用いて代表点から第１の近似平面を算出し、代表点のうち第１の近似平面からの距離が閾値以下であるものを除外し、最小二乗法を用いて除外されていない代表点から第２の近似平面を算出してもよい。この場合、ステップＳ３０３において、検出部１５５は二つの近似平面の境界が処理対象の単位領域に含まれかつ第１の条件が満たされるか否かを判定する。また、ステップＳ３０５において、検出部１５５は二つの近似平面の境界が処理対象の単位領域に含まれかつ第２の条件が満たされるか否かを判定する。

【0102】

上述した説明では、近似曲面算出処理のステップＳ４０１において、算出部１５６が直前の段階において代表点に設定された代表点に、直前の段階において代表点に設定されていない代表点よりも大きい重みを設定するものとしたが、このような例に限られない。算出部１５６は、直前の段階において代表点に設定された代表点に他の代表点と同じ重みを設定してもよく、他の代表点よりも小さい重みを設定してもよい。

【0103】

同様に、算出部１５６は傾斜変化点を含むと推定された単位領域の代表点に他の代表点と同じ重みを設定してもよい。この場合、地表点抽出処理のステップＳ１０５の推定処理が省略されてもよい。このようにしても、地表点抽出装置１は、非地表点が密集している領域を含む３次元点群から地表点を高精度に抽出することができる。また、算出部１５６は傾斜変化点を含むと推定された単位領域の代表点に他の代表点よりも小さい重みを設定してもよい。

【0104】

上述した説明では、近似曲面算出処理のステップＳ４０４において、算出部１５６は、傾斜変化点を含まない単位領域の代表点のうち、直前の段階において代表点に設定されていない代表点のみの重みを更新するものとしたが、このような例に限られない。算出部１５６は、傾斜変化点を含む単位領域の代表点または直前の段階において代表点に設定されていた代表点の重みを更新してもよい。ステップＳ４０１において、傾斜変化点を含む単位領域の代表点または直前の段階において代表点に設定されていた代表点には他の代表点よりも大きい重みが設定される。したがって、ステップＳ４０２で算出される仮の近似曲面に対するこれらの代表点の相対高度は、他の代表点の相対高度よりも小さくなり、ステップＳ４０４において、これらの代表点の重みは他の代表点の重みよりも大きくなるように更新される。その結果、ステップＳ４０５において、地表面の形状がより適切に反映された地表面が算出されるようになる。

【0105】

上述した近似曲面算出処理のステップＳ４０３－Ｓ４０５が省略されてもよい。すなわち、算出部１５６は、仮の近似曲面の算出および代表点の重みの更新をすることなく、ステップＳ４０１で設定された重みに基づいて近似曲面を算出してもよい。

【0106】

上述した説明では、地表点抽出処理において、ステップＳ１０２－Ｓ１０８の処理は複数の段階において実行されるものとしたが、一度だけ実行されるものとしてもよい。このようにしても、地表点抽出装置１は、地表面の傾斜が急変する領域を含む３次元点群から地表点を高精度に抽出することができる。

【0107】

上述した説明では、３次元点群は飛行体から対象領域を計測して得られるものとしたが、航空写真等の画像をステレオ解析または多視点画像解析して得られるものであってもよい。この場合、構成点テーブルＴ１はパルス種別を記憶しなくてもよい。

【0108】

以上説明したように、地表点抽出装置１は、対象領域を、後段階の単位領域が前段階の単位領域よりも小さくなるように設定された複数の単位領域に区分し、複数の段階のそれぞれにおいて、各候補点が示す高度の低さに応じた指標に基づいて代表点を設定し、代表点に対する近似曲面を算出する。また、地表点抽出装置１は、複数の段階のうちの最終段階以外の段階において、近似曲面に対する相対高度が所定範囲である構成点を候補点として抽出し、複数の段階のうちの最終段階において、近似曲面に対する相対高度が所定範囲である構成点を地表点として抽出する。このようにすることで、地表点抽出装置１は、非地表点が密集している領域を含む３次元点群から地表点を高精度に抽出することを可能とする。

【0109】

また、地表点抽出装置１は、少なくとも１つの候補点を含む単位領域のみにおいて代表点を設定する。このようにすることで、非地表点が密集しているために候補点が設定または抽出されなかった単位領域においては代表点が設定されなくなるため、地表点抽出装置１は地表面の形状を適切に反映した近似曲面を算出し、地表点をより高精度に抽出することを可能とする。

【0110】

また、地表点抽出装置１は、直前の段階において代表点に設定された候補点が存在する場合にその候補点を代表点に設定するとともに、直前の段階において代表点に設定されていない候補点のうち最小の指標を有する候補点とを代表点に設定する。直前の段階において代表点に設定された候補点は、処理対象の単位領域が広い段階において代表点に設定されているため、地表点である可能性が特に高い。地表点抽出装置１は、各段階において、地表点である可能性が特に高い候補点を代表点に設定することにより、地表面の形状をより適切に反映した近似曲面を算出し、地表点をより高精度に抽出することを可能とする。

【0111】

また、地表点抽出装置１は、直前の段階において代表点に設定された代表点に、直前の段階において代表点に設定されていない代表点よりも大きい重みを設定する。地表点抽出装置１は、地表点である可能性が特に高い代表点に他の代表点よりも大きい重みを設定することにより、地表面の形状をより適切に反映した近似曲面を算出し、地表点をより高精度に抽出することを可能とする。

【0112】

また、地表点抽出装置１は、近似曲面に対する相対高度が高い代表点ほど重みが小さくなるように代表点の重みを更新し、更新された重みに基づいて近似曲面を更新する。これにより、地表点抽出装置１は、地表面の形状をより適切に反映した近似曲面を算出し、地表点をより高精度に抽出することを可能とする。

【0113】

また、地表点抽出装置１は、各候補点の高度が低いほど小さくなり、かつ各候補点の水平位置が候補点を含む単位領域の中心に近いほど小さくなる指標を算出する。これにより、単位領域の境界付近の候補点が代表点に設定されにくくなるため、地表点抽出装置１は、地表面の形状をより適切に反映した近似曲面を算出し、地表点をより高精度に抽出することを可能とする。

【0114】

また、地表点抽出装置１は、対象領域を複数の単位領域に区分し、傾斜変化点を含む単位領域の代表点の重みを他の単位領域の代表点の重みよりも大きくなるようにして近似曲面を算出し、近似曲面に対する相対高度が所定範囲である構成点を抽出する。このようにすることで、地表点抽出装置１は、地表面の傾斜が急変する領域を含む３次元点群から地表点を高精度に抽出することを可能とする。

【0115】

また、地表点抽出装置１は、傾斜変化点を含む単位領域の代表点の重みが他の単位領域の代表点の重みよりも大きくなるようにして仮の近似曲面を算出し、仮の近似曲面に対する相対高度が高いほど重みが小さくなるように重みを更新し、更新された重みに基づいて近似曲面を算出する。このようにすることで、地表点抽出装置１は、地表面の形状を適切に反映した近似曲面を算出し、地表点をより高精度に抽出することを可能とする。

【0116】

当業者は、本発明の範囲から外れることなく、様々な変更、置換及び修正をこれに加えることが可能であることを理解されたい。例えば、上述した実施形態及び変形例は、本発明の範囲において、適宜に組み合わせて実施されてもよい。

【符号の説明】

【0117】

１ 地表点抽出装置
１５１ 取得部
１５２ 区分部
１５３ 選択部
１５４ 設定部
１５５ 検出部
１５６ 算出部
１５７ 候補点抽出部
１５８ 地表点抽出部
１５９ 出力部

【要約】

【課題】３次元点群から地表点を高精度に抽出することを可能とする地表点抽出装置等を提供する。
【解決手段】地表点抽出装置は、飛行体から対象領域を計測して得られた３次元点群に含まれる複数の構成点を取得するとともに、取得された複数の構成点を候補点に設定し、各段階において、対象領域を、直前の段階の単位領域よりも小さい複数の単位領域に区分し、各候補点について、候補点の高度の低さに応じた指標を算出し、各段階において、指標に基づいて、区分された各単位領域に含まれる候補点のうちから各単位領域の代表点を設定し、各段階において、各単位領域の代表点から近似曲面を算出する算出手段と、最終段階以外の段階において、取得された複数の構成点のうち、近似曲面に対する相対高度が閾値未満である構成点を候補点として抽出し、最終段階において、取得された複数の構成点のうち、相対高度が所定範囲である構成点を地表点として抽出する。
【選択図】図３

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】