(57)【要約】【課題】処理の効率化を図ることのできるレーザデータのフィルタリング方法を提供することを目的とする。【解決手段】上空を飛行する飛行体から地上をレーザ掃射して取得したレーザデータ群に平面座標に基づくメッシュを設定した後、各メッシュ内の最低標高点データをメッシュ代表点１として抽出し、次いで、隣接するメッシュのメッシュ代表点１を近傍として標高要素に基づいてクラスタリングし、この後、構成メッシュ数が最大建物想定面積相当数以上のクラス内に属するレーザデータを地表面反射成分として抽出する。【選択図】  図１
2004272688.tif 000002
【特許請求の範囲】
【請求項１】
上空を飛行する飛行体から地上をレーザ掃射して取得したレーザデータ群に平面座標に基づくメッシュを設定した後、
各メッシュ内の最低標高点データをメッシュ代表点として抽出し、
次いで、隣接するメッシュのメッシュ代表点を近傍として標高要素に基づいてクラスタリングし、
この後、構成メッシュ数が最大建物想定面積相当数以上のクラス内に属するレーザデータを地表面反射成分として抽出するレーザデータのフィルタリング方法。
【請求項２】
上空を飛行する飛行体から地上をレーザ掃射して取得したレーザデータ群に平面座標に基づくメッシュを設定した後、
各メッシュ内の最低標高点データをメッシュ代表点として抽出し、
次いで、隣接するメッシュのメッシュ代表点を近傍として標高要素に基づいてクラスタリングし、
この後、構成メッシュ数が最大建物想定面積相当数以下、最小建物想定面積以上のクラス内のメッシュ代表点の平面位置及び標高値の統計的代表値を平面位置及び標高値とするクラス代表点を設定し、
次いで、前記クラス代表点との間の傾斜角度が所定値以下である前記クラスの近傍における地表面反射レーザデータの割合を求め、
前記割合が所定値以上である場合に前記対象クラス内のレーザデータを地表面反射成分として抽出するレーザデータのフィルタリング方法。
【請求項３】
上空を飛行する飛行体から地上をレーザ掃射して取得したレーザデータ群に平面座標に基づくメッシュを設定した後、
各メッシュ内の最低標高点データをメッシュ代表点として抽出し、
次いで、隣接するメッシュのメッシュ代表点を近傍として標高要素に基づいてクラスタリングし、
この後、構成メッシュ数が最小建物想定面積相当数以下のクラス内で、かつ、セカンドパルス以降のレーザデータからなるメッシュ代表点との間の傾斜角度が所定値以下である前記メッシュ代表点の近傍における地表面反射レーザデータ及び他のセカンドパルス以降のレーザデータの割合を求め、
前記割合が所定値以上である場合に前記対象クラス内のレーザデータを地表面反射成分として抽出するレーザデータのフィルタリング方法。
【請求項４】
前記傾斜角度算定に際して使用される地表面反射成分は、請求項１記載のレーザデータのフィルタリング方法により地表面反射成分として抽出されたクラス内のメッシュ代表点である請求項２または３記載のレーザデータのフィルタリング方法。
【請求項５】
前記クラスタリングが、注目メッシュのメッシュ代表点との間の傾斜角度が所定値以下であるメッシュ代表点の数を基準に行われる請求項１ないし４のいずれかに記載のレーザデータのフィルタリング方法。
【請求項６】
請求項１ないし５のいずれかに記載のレーザデータのフィルタリング方法により地表面反射成分として抽出されたメッシュ代表点の属するメッシュ間を補間して形成した地表面モデルから等高線図と、適宜の断面取得線で切断した断面図とを作成し、
次いで、前記断面図上に前記断面取得線が通過するメッシュを中心として断面幅方向メッシュ数を変化させながら対象メッシュ内の全レーザデータを表示するフィルタリング結果の表示方法。
【請求項７】
上空を飛行する飛行体から地上をレーザ掃射して取得したレーザデータ群に平面座標に基づくメッシュを設定するステップと、
各メッシュ内の最低標高点データをメッシュ代表点として抽出するステップと、
隣接するメッシュのメッシュ代表点を近傍として標高要素に基づいてクラスタリングするステップと、
構成メッシュ数が最大建物想定面積相当数以上のクラス内に属するレーザデータを地表面反射成分として抽出するステップとをコンピュータに実行させるためのレーザデータのフィルタリングプログラム。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、レーザデータのフィルタリング方法に関するものである。
【０００２】
【従来の技術】
航空機レーザスキャナ等により取得したレーザデータから地表面データ、あるいはそれ以外の地物データを抽出することは、例えば、ディジタル標高モデル（ＤＴＭ：Ｄｉｇｉｔａｌ  Ｔｅｒｒａｉｎ  Ｍｏｄｅｌ）を作成する場合に必要であり、このような方法としては、例えば、特許文献１に記載されたものが知られている。
【０００３】
この従来例において、取得されたレーザデータ群はクラスタリングにより類似した特徴を持つグループ（クラスタ）に分類される。クラスタリングのために設定される特徴空間には、各レーザデータの平面位置、標高を要素とする３次元空間がそのまま利用される。
【０００４】
また、クラスタリングに際して定義される近傍にはボロノイ近傍が使用され、レーザデータ群が形成する特徴空間は、先ず、ボロノイ分割により各レーザデータの近傍を定義し、探索される。
【０００５】
【特許文献１】
特開２００２−３６６９７７号公報
【０００６】
【発明が解決しようとする課題】
しかし、上述した従来例において、クラスタリングに際する近傍点の探索は、全レーザデータに対して行う必要があるために、膨大な作業量を必要とするという欠点がある。
【０００７】
本発明は、以上の欠点を解消すべくなされたものであって、処理の効率化を図ることのできるレーザデータのフィルタリング方法及びプログラムを提供することを目的とする。
【０００８】
また、本発明の他の目的は、このようにしてフィルタリングした結果の表示方法の提供にある。
【０００９】
【課題を解決するための手段】
本発明によれば、上記目的は、
上空を飛行する飛行体から地上をレーザ掃射して取得したレーザデータ群に平面座標に基づくメッシュを設定した後、
各メッシュ内の最低標高点データをメッシュ代表点１として抽出し、
次いで、隣接するメッシュのメッシュ代表点１を近傍として標高要素に基づいてクラスタリングし、
この後、構成メッシュ数が最大建物想定面積相当数以上のクラス内に属するレーザデータを地表面反射成分として抽出するレーザデータのフィルタリング方法を提供することにより達成される。
【００１０】
航空機レーザ測量で上空からのレーザ掃射により取得されるレーザデータは、平面位置（Ｘ，Ｙ）情報及び標高（Ｚ）情報を有するポイントデータであり、各レーザデータ単位では地表面等の識別を行うことはできない。このため、地表面等の識別には、隣接するレーザデータとの位置関係を利用することが必要であり、レーザデータ群に対する種々の処理が必要になる。
【００１１】
そして、レーザデータ群に対する処理は、上述した従来例をはじめとして、一般に各レーザデータを平等に扱い、これらに対する統計的処理として行われていた。
【００１２】
一方、レーザデータを種々の目的によって分類、抽出するレーザデータのフィルタリングに当たって、先ず、地表面データの抽出を優先することによって、地表面であることの特殊性をレーザデータ群の処理対象量の軽減に利用する可能性が生じる。すなわち、レーザデータ群に設定された好適なサイズのメッシュ内における最低標高点は、地表面である確率が高いために、これら最低標高点に処理対象を限定して処理対象量を軽減した方が、地表同士を検索の対象とするので、地表抽出精度は高くなる。
【００１３】
本発明は、以上の認識の下になされたものであって、メッシュ化及びメッシュ内でのメッシュ代表値の設定を行うことによって処理対象量を大幅に軽減させるとともに、クラスタリングに伴う近接点の検索を隣接メッシュに限定することによって処理効率を高め、かつ、正確な地表面抽出を可能にするものである。
【００１４】
すなわち、本発明において、先ず、レーザデータ群に対して平面座標に基づくメッシュを設定するステップ（Ｓ２）が実行される。メッシュサイズは、レーザデータが含まれない空メッシュが発生することのないように十分な大きさで、かつ、過度に周囲の条件に平均化されてメッシュの個性が失われない程度の小ささであることが必要であり、レーザデータの密度、正確にはポイントデータ間の間隔を考慮して、例えば、ポイントデータの間隔の２倍程度に決定される。
【００１５】
次いで、上記各メッシュ内のレーザデータのうち、最低標高を有するものをメッシュ代表点１として抽出するステップ（Ｓ３）が実行された後、隣接するメッシュのメッシュ代表点１を近傍として標高要素に基づいてクラスタリングするステップ（Ｓ４）が実行される。メッシュ代表点１を最低標高点データとすることにより、車、樹木等を予め除去することも可能になる。
【００１６】
クラスタリング操作は、標高要素を判断基準として隣接するメッシュとの間で行われる。近傍は、注目メッシュを包囲するように隣接する８個であっても、注目メッシュに対して縦横に隣接する４個であってもよいが、メッシュ代表点１が最低標高点により定義されているために、縦横の４個をとれば十分な精度を維持可能であり、かつ、処理量も少なくてすむ。
【００１７】
また、標高要素による判断は、例えば、隣接するメッシュ代表点１間の比高差をとることも可能であるが、注目メッシュのメッシュ代表点１との間の傾斜角度が所定値以下であるメッシュ代表点１の数を基準に行うことができる。
【００１８】
メッシュ代表点１は、メッシュ内のレーザデータから抽出されたものであり、Ｘ，Ｙ値及びＺ値をそのまま保持しているために、クラスタリングの基準として傾斜角度を利用すると、メッシュ全体を同一標高として考えることとなる比高差を利用する場合に比して、より現実の状態に即した判断が可能になる。
【００１９】
閾値として利用される傾斜角度は、経験則上、水平距離５（ｍ）に対する標高差が０．５（ｍ）、すなわち、５．７°程度に設定するのが望ましい。
【００２０】
以上のクラスタリングステップの後に、構成メッシュ数が最大建物想定面積相当数以上のクラス内に属するレーザデータを地表面反射成分として抽出するステップ、すなわち、クラス内のメッシュ数による地表判定ステップ（Ｓ５）が実行される。
【００２１】
最大建物想定面積は、通常想定される建物の最大面積であり、通常、５０００ｍ２程度に設定しておけば足りる。最大建物想定面積に相当するメッシュ数、すなわち、最大建物想定面積をＳ、１メッシュの面積をｓ、クラス内メッシュ数をＭとした場合、
Ｍ＞Ｓ／ｓ
の条件を満たすクラスは、全体に比較的なだらかな傾斜が建物と想定される場合の面積を超えて連続していることとなるために、非常に高い精度で地表部分であることを推定できる。
【００２２】
また、以上のフィルタリング方法は、なだらかで広い平野部の抽出に適しているが、
上空を飛行する飛行体から地上をレーザ掃射して取得したレーザデータ群に平面座標に基づくメッシュを設定した後、
各メッシュ内の最低標高点データをメッシュ代表点１として抽出し、
次いで、隣接するメッシュのメッシュ代表点１を近傍として標高要素に基づいてクラスタリングし、
この後、構成メッシュ数が最大建物想定面積相当数以下、最小建物想定面積以上のクラス内のメッシュ代表点１の平面位置及び標高値の統計的代表値を平面位置及び標高値とするクラス代表点２を設定し、
次いで、前記クラス代表点２との間の傾斜角度が所定値以下である前記クラスの近傍における地表面反射レーザデータの割合を求め、
前記割合が所定値以上である場合に前記対象クラス内のレーザデータを地表面反射成分として抽出するレーザデータのフィルタリング方法により、クラスの面積により直ちに地表面と判断できない場合の地表抽出方法が提供される。
【００２３】
この方法において、クラスタリングまでのステップ（Ｓ２〜Ｓ４）は上述した方法と全く同様に行われ、この後、構成メッシュ数が最大建物想定面積相当数以下、最小建物想定面積以上のクラスを検索するステップ（Ｓ６１）が実行され、次いで、該当クラス内にクラス代表点２を設定するステップ（Ｓ６２）が実行される。
【００２４】
クラス代表点２は、該当クラス内の複数のメッシュ代表点１を代表する平面位置及び標高値を有する点であり、メッシュ代表点１に対して適宜の統計的処理を施した統計的代表値として与えられる。統計的代表値を求めるためには、平均、あるいはメディアン値等を使用できる。
【００２５】
以上のようにして求められたクラス代表点２は、該当クラスを代表する平面位置及び標高情報を有しており、次いで、これらクラス代表点２とクラス近傍における地表面反射レーザデータとの間の傾斜角度を求めるステップ（Ｓ６３）が実行される。
【００２６】
傾斜角度を求めるに当たって、先ず、クラス代表点２から、適宜の手段によって地表面と判断されたレーザデータとの間で形成されるＴＩＮ（Ｔｒｉａｇｕｌａｔｅｄ  Ｉｒｒｅｇｕｌａｒ  Ｎｅｔｗｏｒｋ）ネットワーク上の各点を近接点として定義することができる。また、この方法は、上述したメッシュ数による地表面抽出方法に続いて行うことが可能であり、この場合、傾斜角度は、メッシュ代表点１を対象として算出される。
【００２７】
該当クラスの地表面判定は、以上のようにしてクラス代表点２と近傍の地表面反射レーザデータの間の傾斜角度を求めた後、傾斜角度が所定値以下の地表面反射データが全測定対象の地表面反射データに占める割合を求めるステップ（Ｓ６４）を実行することにより行われる。
【００２８】
この方法における目的は、建物が地表面からの判断が曖昧なクラスからの地表面の判定であるために、判断の閾値となる傾斜角度は、地表面の想定される最大傾斜角度を基準に設定することができ、例えば、３０°に設定される。
【００２９】
また、割合の閾値は、判定する地形等を考慮して判定精度が高くなるように適宜決定可能であるが、１／３程度にするのが経験則上望ましく、閾値以下の傾斜角度の割合が上記閾値以上である場合には、当該クラスは地表面と判定される。
【００３０】
このように、判定に当たってクラス内の各点を対象とすることなく、クラス代表点２を基準にクラス全体を判定することにより、計算量を少なくすることができる。
【００３１】
さらに、以上のフィルタリング方法に加え、
上空を飛行する飛行体から地上をレーザ掃射して取得したレーザデータ群に平面座標に基づくメッシュを設定した後、
各メッシュ内の最低標高点データをメッシュ代表点１として抽出し、
次いで、隣接するメッシュのメッシュ代表点１を近傍として標高要素に基づいてクラスタリングし、
この後、構成メッシュ数が最小建物想定面積相当数以下のクラス内で、かつ、セカンドパルス以降のレーザデータからなるメッシュ代表点１との間の傾斜角度が所定値以下である前記メッシュ代表点１の近傍における地表面反射レーザデータ及び他のセカンドパルス以降のレーザデータの割合を求め、
前記割合が所定値以上である場合に前記対象クラス内のレーザデータを地表面反射成分として抽出するレーザデータのフィルタリング方法を採用することにより、樹木等の多い山間部での地表面の抽出精度を高めることができる。
【００３２】
この方法は、上述した方法において、構成メッシュ数が最小建物想定面積相当数以下となったクラスに対する有効な処理を提供するもので、レーザスキャンに際して掃射されたレーザは複数の反射点から反射され、各々がリターン情報とともに受領される点が利用される。すなわち、対象物に到達したレーザは、所定のスポット径を有しているために、対象物に孔、隙間等がある場合には、最高点において先ず反射し（ファーストリターン）、この後、木漏れ日のように徐々に標高方向に進行し、セカンド、サードリターンとして反射する。このようなセカンドリターン以降のレーザデータは、性質上樹木が植生している山間部等において取得されることが多く、深度の深い、すなわち、リターン数の多いものほど地表面からの反射レーザデータである確率が高いこととなる。
【００３３】
一方、ファーストパルスが樹木の頂部で反射したとも限らないし、ラストパルスが地表面で反射したとも限らない上に、ラストパルスは反射強度が弱いために誤差成分が多く含まれやすい。したがって、これらセカンドリターン以降のレーザデータを扱うには、統計的解析が必要となる。
【００３４】
そこで、本方法において、先ず、当該クラス内でのセカンドパルス以降のレーザデータの抽出ステップ（Ｓ７１）を実行した後、抽出された各レーザデータと、近傍の地表反射レーザデータ及び他のセカンドパルス以降のレーザデータとの間で、上述した方法と同様に、傾斜角度算出ステップ（Ｓ７２）と割合による判定ステップ（Ｓ７３）を実施し、所定の傾斜角度以下の割合が所定の閾値以上である場合に地表面と判定する。
【００３５】
また、以上のようにして行われたフィルタリングの結果を確認するためには、
地表面反射成分として抽出されたメッシュ代表点１の属するメッシュ間を補間して形成した地表面モデルから等高線図３と、適宜の断面取得線４で切断した断面図とを作成し、
次いで、前記断面図上に前記断面取得線４が通過するメッシュを中心として断面幅方向メッシュ数を変化させながら対象メッシュ内の全レーザポイントデータを表示するフィルタリング結果の表示方法を使用することができる。
【００３６】
【発明の実施の形態】
図１、２にレーザデータからＤＴＭを作成する方法に適用された本発明の実施の形態を示す。レーザデータは、上空からのレーザ掃射により得られた地上物（地表、地物）の平面位置及び標高値を有するポイントデータであり、地表面判定に当たり、予め建物部分のデータをポリゴンデータ等として有している場合には、前処理として建物部分を除去し（ステップ（Ｓ１））、次いで、平面視において矩形のメッシュを設定する（ステップ（Ｓ２））。レーザデータ間の平均間隔がほぼ２．５（ｍ）であることを考慮して、メッシュサイズは５（ｍ）とする。
【００３７】
次いで、ステップＳ３において、各メッシュ内で標高値が最低のレーザデータを検索し、メッシュ代表点１として指定する。図３（ａ）にレーザデータ群に設定された各メッシュ内に１個のメッシュ代表点１を指定した状態を示す。メッシュ代表点１は、生データから抽出されたデータであり、抽出前の平面位置および標高値をそのまま保持する。
【００３８】
この後、ステップ（Ｓ４）において、クラスタリングを行う。クラスタリングは、高さ方向傾斜量が所定の閾値内であれば同一クラスとする処理を行うことによってなされ、先ず、図３（ｂ）に示すように、注目するメッシュ代表点１と、このメッシュ代表点１が属するメッシュに対して縦方向及び横方向に隣接する４個のメッシュ内のメッシュ代表点１との間の高さ方向の傾斜角を求め、所定の閾値と比較する。閾値には、５．７°が採用される。図３（ｃ）において縦横隣接するメッシュ内のメッシュ代表点１間の傾斜角が閾値以下のもののみを線分で結んで示す。
【００３９】
また、クラスタリングに際して、各メッシュ代表点１からは、隣接するメッシュが空メッシュでない限り、隣接するメッシュ内のメッシュ代表点１に向けて４本の傾斜線が引かれるが、このうち、閾値以下の傾斜角のものが所定の本数以上（この実施の形態においては１本）以上あるものが同一クラスタとして取り扱われる。
【００４０】
クラスタの境界は、直接隣接しているメッシュ間が閾値以上の傾斜であっても、回り込んで隣接する場合には同一のグループとして分類される。すなわち、図４（ａ）において、点ＡとＤとの間には閾値を超える傾斜があるが、点Ｂ、Ｃ、Ｄを結ぶ閾値以下の傾斜ルートが存在するために、Ａ〜Ｄが属するメッシュは同一クラスに分類される。これに対し、点Ｅ、Ｆが属するメッシュは、Ａ〜Ｄが属するメッシュとは別のクラスに分類される。なお、図４（ａ）においてメッシュ代表点１を黒点で、メッシュ代表点１間の傾斜が閾値以下の場合には、メッシュ代表点１間を実線で、閾値以上の場合には、破線で示される。
【００４１】
また、図４（ｂ）に、図３（ｃ）の結果から、クラスタリングした状態を、同一クラスタを連続する模様を施して示す。
【００４２】
この後、ステップ（Ｓ５）において、メッシュ数を基準とした地表判定を行う。上述したように、地表判定基準は、メッシュ数＝面積を示すことに着目し、最大建物想定面積より大きなクラスを地表面として、最小建物想定面積より小さなクラスを樹木等の建物以外の地物として抽出する。図に最大建物想定面積を５０００ｍ２＝２５（単位メッシュ面積）×２００（メッシュ）、最小建物想定面積を２００ｍ２＝２５×８（メッシュ）と設定した場合の地表面、建物及び地物を示す。
【００４３】
図４（ｂ）の結果から、メッシュ数２００以上で、地表と判定されたクラスを図５（ａ）において右上がりのハッチングを施して示す。
【００４４】
以上のステップ（Ｓ１）からステップ（Ｓ５）まではコンピュータにより実行可能であり、これらのステップを記述したプログラムとして独立して使用することも可能である。
【００４５】
この実施の形態において、上記ステップ（Ｓ５）に引き続き、ステップ（Ｓ６）において傾斜量とメッシュ数による地表判定が行われる。このステップは、上記ステップ（Ｓ５）において最大建物想定面積と最小建物想定面積の中間にあるクラスから地表面を抽出するもので、このステップにおける処理対象を図５（ａ）において右下がりのハッチングを施して示す。
【００４６】
傾斜量とメッシュ数による地表判定を行うに際して、先ず、ステップ（Ｓ６１）において上記条件を満足するクラスを検索し、選択し、次いで、対象クラス内のメッシュ代表点１を相加平均してクラス代表点２を設定した（ステップ（Ｓ６２））後、ステップ（Ｓ６３）でクラス代表点２と、ステップ（Ｓ５）において地表面と判定されたメッシュ内のメッシュ代表点１との間の傾斜角度を算出する。図５（ｂ）に対象クラスと、傾斜角度の対象候補である地表面クラス内のメッシュ代表点１を示し、図５（ｃ）にクラス代表点２を白丸で、傾斜角度算出部を線分で示す。
【００４７】
クラス代表点２からの傾斜角度は、図６（ａ）に示すように、クラス内のメッシュ数が１個で、隣接するメッシュの中心との傾斜角がαであっても、図６（ｂ）に示すように、メッシュ数が４（幅方向２メッシュ）になった場合には、中心点から隣接メッシュの中心までの傾斜角（β）は、αに対して、
β＝ｔａｎ−１（２×ｔａｎα／３）、
メッシュ数が９（幅方向３メッシュ）になった場合（図６（ｃ））には、
β＝ｔａｎ−１（２×ｔａｎα／４）、
メッシュ数が１６（幅方向４メッシュ）の場合（図６（ｄ））には、
β＝ｔａｎ−１（２×ｔａｎα／５）というように、
低下する。
【００４８】
これらαとβとの間には、一般に
メッシュサイズをＭ、要素数をＮとして、
β＝ｔａｎ−１（２×ｔａｎα／（ｓｑｒｔ（Ｎ）＋１））・・・（１）
の関係があるため、式（１）において適宜の管理角度（α）を代入して閾値とする。なお、図６において縦線はメッシュ境界を示す。
【００４９】
このようにして傾斜角度の算出が終了すると、次にステップ（Ｓ６４）で地表面の判定が行われる。地表判定は、クラス代表点２とのなす傾斜量が上記閾角度以内である傾斜角度算出対象点の全体に対する割合が所定値以上か否かにより行われる。
【００５０】
図５（ｃ）は、ステップ（Ｓ６４）における判定方法を示すもので、図５（ｃ）においてクラス代表点２から上方に線分が引かれていることからも明らかなように、図の上方にも地表面クラスが存在し、角度測定対象点はクラス代表点２から６本の線分が延びていることから、全部で６点あることが分かる。
【００５１】
上述したように、クラス代表点２は、平面位置及び標高値を有しており、この点と角度測定対象点との標高傾斜を例えば、α＝３０°の条件で仕分け、それ以下の傾斜角度を有する角度測定対象点が図５（ｃ）において実線で示すように、４点であったとすると、上記閾角度以内の隣接点数の割合は、
４／６＝０．６６６となり、
地表判定の閾値を０．３３３とすると、
０．６６６＞０．３３３
より、このクラスは地表面であると判定される。
【００５２】
なお、このステップ（Ｓ６）は、ステップ（Ｓ１）において建物部分が完全に除去されている場合には、省略することができる。また、以上のステップ（Ｓ１〜Ｓ４及びＳ６）を独立にコンピュータに実行させるようなプログラムを作ることも可能であり、さらに、ステップ（Ｓ１）からステップ（Ｓ６）までを連続してコンピュータに実行させるプログラムも作ることができる。
【００５３】
以上のようにして一応の地表面抽出を行った後、ステップ（Ｓ７）でリターン情報によるパルス解析を行う。
【００５４】
パルス解析は、セカンドパルス以下のパルスを利用して上記ステップ（Ｓ５）で建物以外の地物と判定されたクラスから地表成分を分離するために行われる。パルス解析に際して、先ず、該当クラス内でセカンドリターン以降のレーザデータを抽出する（ステップ（Ｓ７１））。図７（ａ）においてセカンドリターン以降のレーザデータを黒丸で、ステップ（Ｓ６）までで地表面と判定されたメッシュ代表点１を白丸で示す。
【００５５】
次いで、ステップ（Ｓ７２）でセカンドリターン以降のレーザデータを中心とする傾斜角度を個別に算出する。傾斜角度の算出に先立って、図７（ｂ）において実線で示すように、セカンドリターン以降のレーザデータ間及び上記地表面反射レーザデータ間でＴＩＮネットワークを形成し、傾斜角度の算出は、ネットワーク対象点間を対象にして行われる。
【００５６】
地表面か否かの判定は、ステップ（Ｓ６）と同様の基準で行われ、処理方法がステップ（Ｓ６）とほぼ同一なので、このステップ（Ｓ７）はステップ（Ｓ６）と同時に行うことができる。
【００５７】
図７（ｂ）に地表面判定の条件を充たさなかったレーザデータを二重丸で、この結果、地表面として新たに指定された領域を図７（ｃ）においてハッチングを施して示す。
【００５８】
なお、以上のステップ（Ｓ１〜Ｓ４及びＳ７）を独立にコンピュータに実行させるようなプログラムを作ることも可能であり、さらに、ステップ（Ｓ１）からステップ（Ｓ７）までを連続してコンピュータに実行させるプログラム、あるいは、ステップ（Ｓ１〜Ｓ５及びＳ７）をコンピュータに実行させるプログラムも作ることができる。
【００５９】
以上のようにして地表面の抽出を行った後、ステップ（Ｓ８）で補間を行って地表面を求める。補間には、スプライン、クリギング等の種々の手法が使用できるが、トポグリッド手法を使用することによって、尾根・谷にめりはりをつけることができる。
【００６０】
さらに、この実施の形態は、以上のようにして求めた地表面の確認（ステップ（Ｓ１０））を含む。確認は、上述したフィルタリング結果を表示し（ステップ（Ｓ９））、この結果を目視することにより不自然な点を発見、修正して行われる。
【００６１】
確認時の手助けとなる結果表示に際し、先ず、ステップ（Ｓ８）により形成された地表面データから等高線図３を作成し（ステップ（９１））、さらに、ユーザの指定した断面取得線４に沿った断面図を作成する（ステップ（９２））。
【００６２】
この後、この断面図上に断面取得線４が通過するメッシュを中心として断面幅方向の所定数のメッシュ内の全レーザデータを表示する（ステップ（９３））。断面図上に表示するレーザデータ数、すなわち、断面取得線４に直交する方向のメッシュ数は適宜変更可能であり、当該メッシュ数を増減させた状態の変化から、地表面データの確からしさを確認できる。
【００６３】
図８はこの状態を示すもので、図８（ａ）に断面取得線４を一点鎖線で示した等高線図を、図８（ｂ）に断面取得線４に沿う断面図と、断面取得線４が通過する１列のメッシュ列内のレーザデータを表示した状態、図８（ｃ）に断面取得線４を中心に上下計３列のメッシュ列内の全レーザデータを表示した状態、図８（ｄ）に上下計５列のメッシュ列内のレーザデータを表示した状態を示す。
【００６４】
このレーザフィルタリング結果の表示方法は、上記ステップ（Ｓ５、Ｓ６、Ｓ７）のいずれかに続いて、ステップ（Ｓ８）を実行し、さらに、与えられた断面取得線４に沿って、与えられた断面幅で、ステップ（Ｓ９）をコンピュータに実行させるフィルタリング結果表示用プログラムを走行させることにより達成できる。
【００６５】
以上のようにしてフィルタリング結果の表示によりＤＴＭに問題がある場合には、閾値を変えてステップ（Ｓ３）からやり直すか、あるいはステップ（Ｓ１１）で手動修正が行われる。
【００６６】
これに対し、確認結果が良好である場合には、必要に応じて図化編集（ステップ（Ｓ１２））、および図面間接合（ステップ（１３））が実行されて作業が終了する。
【００６７】
【発明の効果】
以上の説明から明らかなように、本発明によれば、処理の効率化を図り、かつ、正確に地表面の抽出等のフィルタリングを行うことができる。
【図面の簡単な説明】
【図１】本発明を示すフローチャートである。
【図２】本発明を示すフローチャートである。
【図３】クラスタリング方法を示す説明図である。
【図４】クラスタリング結果を示す説明図である。
【図５】傾斜量とメッシュ数による地表判定手順を示す説明図である。
【図６】傾斜角の換算を示す説明図である。
【図７】リターン情報によるパルス解析の手順を示す説明図である。
【図８】フィルタリングの結果表示を説明する図である。
【符号の説明】
１        メッシュ代表点
２        クラス代表点
３        等高線図
４        断面取得線
【図１】
2004272688.tif 000003
【図２】
2004272688.tif 000004
【図３】
2004272688.tif 000005
【図４】
2004272688.tif 000006
【図５】
2004272688.tif 000007
【図６】
2004272688.tif 000008
【図７】
2004272688.tif 000009
【図８】
2004272688.tif 000010