(57)【特許請求の範囲】
【請求項１】
測量データに基づいて地表面を抽出する測量データ処理装置であって、
対象物の計測点について計測した、三次元計測データの集まりである点群データを取得する点群データ取得手段と、
対象物の延長方向に所定の間隔をおいて複数の横断面を設定し、当該横断面に近接する計測点を横断面上に投影する設定投影手段と、
前記横断面上の計測点に基づいて、平坦な地表面を抽出する平坦地表面抽出手段と、
を備えた測量データ処理装置において、
前記平坦地表面抽出手段は、前記横断面を高さ方向に複数のエリアに分割し、標高の低いエリアから順に、当該エリアに含まれる計測点を取得し、取得した計測点の横断面における幅方向の距離が所定値以下であれば、当該取得した計測点を選択計測点とし、当該幅に含まれる他の全ての計測点を削除し、
当該選択計測点に基づいて平坦地表面を抽出することを特徴とする測量データ処理装置。
【請求項２】
測量データに基づいて地表面を抽出する測量データ処理装置をコンピュータによって実現するための測量データ処理プログラムであって、コンピュータを、
対象物の計測点について計測した、三次元計測データの集まりである点群データを取得する点群データ取得手段と、
対象物の延長方向に所定の間隔をおいて複数の横断面を設定し、当該横断面に近接する計測点を横断面上に投影する設定投影手段と、
前記横断面上の計測点に基づいて、平坦な地表面を抽出する平坦地表面抽出手段として機能させるための測量データ処理プログラムにおいて、
前記平坦地表面抽出手段は、前記横断面を高さ方向に複数のエリアに分割し、標高の低いエリアから順に、当該エリアに含まれる計測点を取得し、取得した計測点の横断面における幅方向の距離が所定値以下であれば、当該取得した計測点を選択計測点とし、当該幅に含まれる他の全ての計測点を削除し、
当該選択計測点に基づいて平坦地表面を抽出することを特徴とする測量データ処理プログラム。
【請求項３】
請求項１の装置または請求項２のプログラムにおいて、
前記設定投影手段は、断面形状データによるブレークラインの方向に所定の間隔をおいて複数の横断面を設定することを特徴とする装置またはプログラム。
【請求項４】
請求項１〜３のいずれかの装置またはプログラムにおいて、
前記平坦地表面抽出手段は、隣接する選択計測点の角度に基づいて平坦地表線を抽出し、各横断面の平坦地表線をつないで平坦地表面とすることを特徴とする装置またはプログラム。
【請求項５】
請求項１〜４のいずれかの装置またはプログラムにおいて、
前記平坦地表面抽出手段は、所定以上の幅があり、最も標高の高い平坦地表線より上部の計測点を削除して平坦にすることを特徴とする装置またはプログラム。
【請求項６】
請求項１、３〜５のいずれかの装置において、さらに、
各横断面の平坦地表線の端部の変化点をつないだブレークラインを抽出するブレークライン抽出手段を備えていることを特徴とする装置。
【請求項７】
請求項２〜５のいずれかのプログラムにおいて、コンピュータを、さらに、
各横断面の平坦地表線の端部の変化点をつないだブレークラインを抽出するブレークライン抽出手段として機能させるためのプログラム。
【請求項８】
測量データに基づく処理を行うために、測量データを選択する測量データ処理装置であって、
対象物の計測点について計測した、三次元計測データの集まりである点群データを取得する点群データ取得手段と、
対象物の延長方向に所定の間隔をおいて複数の横断面を設定し、当該横断面に近接する計測点を横断面上に投影する設定投影手段と、
前記横断面を高さ方向に複数のエリアに分割し、標高の低いエリアから順に、当該エリアに含まれる計測点を取得し、取得した計測点の横断面における幅方向の距離が所定値以下であれば、当該計測点を選択計測点とし、当該幅に含まれる他の全ての計測点を削除し、当該選択計測点に基づいて平坦な地表面を抽出する平坦地表面を抽出する選択する計測点選択手段と、
を備えた測量データ処理装置。
【請求項９】
測量データに基づく処理を行うために、測量データを選択する測量データ処理装置をコンピュータによって実現するための測量データ処理プログラムであって、コンピュータを、
対象物の計測点について計測した、三次元計測データの集まりである点群データを取得する点群データ取得手段と、
対象物の延長方向に所定の間隔をおいて複数の横断面を設定し、当該横断面に近接する計測点を横断面上に投影する設定投影手段と、
前記横断面を高さ方向に複数のエリアに分割し、標高の低いエリアから順に、当該エリアに含まれる計測点を取得し、取得した計測点の横断面における幅方向の距離が所定値以下であれば、当該計測点を選択計測点とし、当該幅に含まれる他の全ての計測点を削除し、当該選択計測点に基づいて平坦な地表面を抽出する平坦地表面を抽出する選択する計測点選択手段として機能させるための測量データ処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は測量データに基づいて、水平地表面やブレークラインなどの特徴形状を抽出するための装置に関するものである。
【背景技術】
【０００２】
河川などの管理において、航空写真などから形状の特徴として水平地表面やブレークライン（地形線）を導くことが行われている。たとえば、作業者が写真を見ながら目視で、水平地表面やブレークラインを描いている。
【０００３】
また、航空レーザ測量（ＬＰ）などの点群データ（三次元地形データ）に基づいて形状を把握することも行われている。この場合、測量によって得られた点群データに基づいて、コンピュータにより水平地表面やブレークラインを生成することが行われている。
【０００４】
たとえば、特許文献１には、点群データに基づいてブレークライン（地形線）を算出することが開示されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００３−３２９４４０
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、作業者が目視によって水平地表面やブレークラインを見いだす方法は、効率が悪く、また、作業者によって抽出の精度がまちまちであった。
【０００７】
一方、特許文献１の技術では、コンピュータによって自動的にブレークラインなどを見いだすことができるという点において、客観性や効率の面から好ましいものである。
【０００８】
しかし、特許文献１の技術では、あらゆる方向におけるブレークラインを考慮しているため、本来必要な方向のブレークラインを抽出するためにはノイズが多く、処理が複雑になるという問題があった。このため、水平地表面も性格に抽出できないという問題があった。加えて、三次元空間から直接的に地表面やブレークラインを抽出しようとしているので、この点においても処理が複雑であるという問題があった。
【０００９】
さらに、水平地表面やブレークラインを生成する際に、測定誤差を含むデータも用いているので、処理が煩雑になるだけでなく、抽出の精度が十分であるとはいえなかった。
【００１０】
この発明は、上記のような問題点の少なくとも一つを解決して、効率的で精度のよい測量データ処理を行うことのできる測量データ処理装置を提供することを目的とする。
【課題を解決するための手段】
【００１１】
独立して適用可能な、本発明のいくつかの特徴を以下に示す。
【００１２】
(1)(2)この発明にかかる測量データ処理装置は、測量データに基づいて地表面を抽出する測量データ処理装置であって、対象物の計測点について計測した、三次元計測データの集まりである点群データを取得する点群データ取得手段と、対象物の延長方向に所定の間隔をおいて複数の横断面を設定し、当該横断面に近接する計測点を横断面上に投影する設定投影手段と、前記横断面上の計測点に基づいて、平坦な地表面を抽出する平坦地表面抽出手段とを備えている。
【００１３】
したがって、横断面における形状に基づいて地表面を抽出しているので、処理効率がよく精度のよい抽出を行うことができる。
【００１４】
(3)この発明にかかる測量データ処理装置は、設定投影手段が、断面形状データによるブレークラインの方向に所定の間隔をおいて複数の横断面を設定することを特徴としている。
【００１５】
したがって、平坦地表面を抽出するための横断面を適切に設定することができる。
【００１６】
(4)この発明にかかる測量データ処理装置は、平坦地表面抽出手段が、前記横断面の幅方向に近接する複数の計測点について、標高の低いものを選択計測点とし、当該選択計測点に基づいて平坦地表面を抽出することを特徴としている。
【００１７】
したがって、計測誤差を低減して正確に平坦地表面を抽出することができる。
【００１８】
(5)この発明にかかる測量データ処理装置は、平坦地表面抽出手段が、前記横断面を高さ方向に複数のエリアに分割し、標高の低いエリアから順に、当該エリアに含まれる計測点を取得し、取得した計測点の横断面における幅方向の距離が所定値以下であれば、全てのエリアに関して当該幅に含まれる計測点のうち前記取得した計測点を選択計測点とすることを特徴としている。
【００１９】
したがって、より詳細に誤差データを排除することができる。
【００２０】
(6)この発明にかかる測量データ処理装置は、平坦地表面抽出手段が、隣接する選択計測点の角度に基づいて平坦地表線を抽出し、各横断面の平坦地表線をつないで平坦地表面とすることを特徴としている。
【００２１】
したがって、精度よく平坦地表面を抽出することができる。
【００２２】
(7)この発明にかかる測量データ処理装置は、平坦地表面抽出手段が、所定以上の幅があり、最も標高の高い平坦地表線より上部の計測点を削除して平坦にすることを特徴としている。
【００２３】
したがって、平坦面の上にあるものに起因したノイズデータを排除することができる。
【００２４】
(8)(9)この発明にかかる測量データ処理装置は、各横断面の平坦地表線の端部の変化点をつないだブレークラインを抽出するブレークライン抽出手段を備えている。
【００２５】
したがって、精度よくブレークラインを抽出することができる。
【００２６】
(10)(11)この発明にかかる測量データ処理装置は、測量データに基づく処理を行うために、測量データを選択する測量データ処理装置であって、対象物の計測点について計測した、三次元計測データの集まりである点群データを取得する点群データ取得手段と、対象物の延長方向に所定の間隔をおいて複数の横断面を設定し、当該横断面に近接する計測点を横断面上に投影する設定投影手段と、前記横断面を高さ方向に複数のエリアに分割し、標高の低いエリアから順に、当該エリアに含まれる計測点を取得し、取得した計測点の横断面における幅方向の距離が所定値以下であれば、全てのエリアに関して当該幅に含まれる計測点のうち前記取得した計測点を選択する計測点選択手段とを備えている。
【００２７】
したがって、ノイズデータを排除して測量データ処理のために必要なデータを選択することができる。
【００２８】
「点群データ取得手段」は、実施形態においては、ステップＳ１がこれに対応する。
【００２９】
「断面設定手段」は、実施形態においては、ステップＳ３がこれに対応する。
【００３０】
「計測点投影手段」は、実施形態においては、ステップＳ５がこれに対応する。
【００３１】
「計測点選択手段」は、実施形態においては、ステップＳ６がこれに対応する。
【００３２】
「角度変化抽出手段」は、実施形態においては、ステップＳ７がこれに対応する。
【００３３】
「ブレークライン生成手段」は、実施形態においては、ステップＳ１０がこれに対応する。
【００３４】
この発明において、「プログラム」とは、ＣＰＵにより直接実行可能なプログラムだけでなく、ソース形式のプログラム、圧縮処理がされたプログラム、暗号化されたプログラム等を含む概念である。
【図面の簡単な説明】
【００３５】
【図１】この発明の実施形態による測量データ処理装置の機能ブロック図である。
【図２】測量データ処理装置のハードウエア構成である。
【図３】測量データ処理プログラムのフローチャートである。
【図４】横断面設定の処理を説明するための図である。
【図５】点群データ抽出のフローチャートである。
【図６】点群データ抽出のフローチャートである。
【図７】点群データの例である。
【図８】点群データ抽出処理を説明するための図である。
【図９】点群データ抽出処理を説明するための図である。
【図１０】横断面に投影された計測点を示す図である。
【図１１】計測点選択のフローチャートである。
【図１２】計測点選択処理を説明するための図である。
【図１３】計測点選択処理を説明するための図である。
【図１４】天端面抽出のフローチャートである。
【図１５】ノイズ排除処理を説明するための図である。
【図１６】天端面生成を説明するための図である。
【図１７】格子データ生成を説明するための図である。
【発明を実施するための形態】
【００３６】
１．全体構成
図１に、この発明の一実施形態による測量データ処理装置の機能ブロック図を示す。点群データ取得手段２は、河川などの対象物について、航空レーザ測量（ＬＰ）などの点群データ（三次元地形データ）１８を取得する。点群データ１８は、各計測点の三次元計測データの集まりである。点群データ１８は、記録部２０に記録される。
【００３７】
この実施形態において、設定投影手段４は、断面設定手段１０と計測点投影手段１２を有している。断面設定手段１０は、河川などの対象物の延長方向に所定間隔を置いて複数の横断面Ｓを設定する。この実施形態では、対象物に沿った候補ブレークラインに所定間隔を置いて複数の横断面Ｓを設定するようにしている。
【００３８】
計測点投影手段１２は、各横断面Ｓについて近接する計測点を選択する。この実施形態では、横断面Ｓに沿って所定区間ΔＬを設定し、横断面Ｓに垂直な方向に選択幅ΔＷを設定している。各所定区間ΔＬにおいて選択される計測点の数は、選択幅ΔＷの大きさによって変わる。つまり、選択幅ΔＷが大きいほど、多くの計測点が選択されることになる。この実施形態では、所定区間ΔＬのそれぞれにおいて、所定基準によって（たとえば、選択される計測点の数所定数を上回るように、所定区間ΔＬ内における計測点のＸ値の標準偏差が所定値以上（所定以上のばらつきとなるよう）となるように）所定区間ΔＬごとに選択幅ΔＷを変化させて、計測点を選択する様にしている。
【００３９】
さらに、計測点投影手段１２は、選択した計測点を、横断面Ｓ上に来るように平行移動する。つまり、高さ方向の計測値を維持しつつ、横断面Ｓ上まで移動する。
【００４０】
したがって、横断面Ｓの近傍に多くの計測点があれば、当該区間では、これら計測点が抽出される。横断面の近傍にある計測点が少ない場合には、当該区間では、より遠くにある計測点も含めて抽出される。
【００４１】
このようにして抽出された計測点に基づいて、平坦地表面抽出手段６が、平坦な地表面を抽出する。この実施形態において、平坦地表面抽出手段６は、計測点選択手段１４と角度変化抽出手段１６を有している。
【００４２】
計測点選択手段１４は、横断面Ｓの計測点Ｐのうち標高の低いもの（図において丸印を付した計測点）を選択する。角度変化抽出手段１６は、このようにして選択した計測点Ｐを並べたときの角度変化をみて、平坦地表面ＦＳを見いだす。
【００４３】
ブレークライン生成手段８は、各横断面Ｓの変化点Ｖを接続してブレークライン２２を生成する。
【００４４】
２．ハードウエア構成
図２に、この発明の一実施形態による測量データ処理装置のハードウエア構成を示す。ＣＰＵ３０に、メモリ３２、ディスプレイ３４、ハードディスク３６、ＤＶＤ−ＲＯＭドライブ３８、マウス／キーボード４０が接続されている。ハードディスク３６には、オペレーティングシステム４２、測量データ処理プログラム４４、横断面形状データ４８などが記録されている。
【００４５】
測量データ処理プログラム４４は、オペレーティングシステム４２と協働してその機能を発揮するものである。これらプログラムは、ＤＶＤ−ＲＯＭ４６に記録されていたものを、ＤＶＤ−ＲＯＭドライブ３８を介して、ハードディスク３６にインストールしたものである。なお、インターネットを介してダウンロードしたものであってもよい。
【００４６】
３．測量データ処理プログラム４４の処理
図３に、測量データ処理プログラム４４のフローチャートを示す。ＣＰＵ３０は、河川や道路などの計測点について計測した三次元計測データの集まりである点群データ１８を取得し、ハードディスク３６に記録する（ステップＳ１）。点群データ１８は、各測量点についてのＸ、Ｙ、Ｚ座標の数値データの集合である。点群データ１８は、ＤＶＤ−ＲＯＭなどの記録媒体から読み込んでもいいし、ネットワークなどを介して取り込むようにしてもよい。
【００４７】
この実施形態では、レーザプロファイラ測量（ＬＰ測量）によって得た点群データを用いている。ＬＰ測量とは、航空機に搭載したレーザスキャナから地上にレーザ光を照射し、反射したレーザ光との時間差により得られる地上までの距離と、ＧＰＳ測量機、ＩＭＵ（慣性計測装置）から得られる位置とに基づいて、地形の三次元データを測量するものである。なお、車などを用いて測量した三次元データを用いることもできる。
【００４８】
3.1横断面の設定と計測点の投影
次に、ＣＰＵ３０は、候補となるブレークラインを設定する（ステップＳ２）。この実施形態では、既に生成されている過去の断面形状データ４８（河川などの所定間隔に生成された断面形状データ）の変化点を結ぶ線（スプライン曲線）によって候補ブレークラインを生成するようにしている。たとえば、図４Ａに示すように、断面形状データ４８の法肩を結ぶ線を候補ブレークライン５０として生成する。
【００４９】
ＣＰＵ３０は、この候補ブレークライン５０を点群データ１８に重ねる。両者の座標位置は、基準となる座標（たとえば距離標の座標）に基づいて合わせるようにすればよい。
【００５０】
なお、過去の断面形状データ４８がない場合などは、航空写真などの地図をディスプレイ３４に表示し、ユーザがマウス４０を操作して、候補ブレークライン５０を入力するようにしてもよい。この場合も、基準となる座標に基づいて両者の座標位置を合わせるようにすればよい。
【００５１】
続いて、ＣＰＵ３０は、図４Ｂに示すように、この候補ブレークライン５０に垂直な横断面Ｓを、候補ブレークライン５０に沿って生成する（ステップＳ３）。なお、この実施形態では、候補ブレークライン５０に対して垂直に横断面Ｓを生成しているが、垂直でなく任意の角度でよい。
【００５２】
次に、ＣＰＵ３０は、生成した横断面Ｓに対し、近傍の計測点Ｐを抽出して投影する（ステップＳ５）。この処理の詳細を、図５、図６に示す。
【００５３】
まず、ＣＰＵ３０は、対象横断面を含む矩形を生成する（ステップＳ５０１）。図７に、点群データ１８を模式的に示す。多数の計測点Ｐが示されている。各計測点Ｐは、Ｘ値、Ｙ値、Ｚ値を有している。なお、図７において、高さを表すＺ値は紙面に垂直な方向であるため表現されていない。
【００５４】
ＣＰＵ３０は、対象横断面Ｓを含む矩形５２を生成する。生成された矩形５２と、これに含まれる点群データ１８の計測点Ｐを、図８Ａに示す。
【００５５】
ＣＰＵ３０は、この矩形５２を、対象横断面Ｓに沿って複数の区間ΔＬに分割する（ステップＳ５０２）。この実施形態では、等間隔の区間ΔＬに分割しているが、異なる大きさの区分ΔＬに分割するようにしてもよい。
【００５６】
次に、ＣＰＵ３０は、対象横断面Ｓを中心とした対象幅ΔＷ1の範囲にある計測点Ｐの数を計数する（ステップＳ５０４）。計測点Ｐの数がしきい値を超えていれば、当該対象幅ΔＷ1の範囲内にある計測点Ｐを選択し（ステップＳ５０８）、これらを用いて横断面形状データの生成を行う。
【００５７】
計測点Ｐの数がしきい値を超えていなければ、対象幅をΔＷ1からΔＷ2に広げる（ステップＳ５０６）。つまり、対象横断面Ｓに垂直な方向に広げる。そして、対象幅ΔＷ2の範囲にある計測点Ｐの数を計数する（ステップＳ５０４）。計測点Ｐの数がしきい値を超えていれば、当該対象幅ΔＷ2の範囲内にある計測点Ｐを選択し（ステップＳ５０８）、これらを用いて横断面形状データの生成を行う。
【００５８】
以上の処理を繰り返し、対象幅に含まれる計測点Ｐが所定数を超えるまで、対象幅を広げていく。このようにすることで、できるだけ対象横断面Ｓに近い計測点Ｐを用いつつ、形状を把握するために必要な数の計測点Ｐを選択することができる。
【００５９】
なお、対象幅を広くした結果、対象幅が矩形５２の幅を超えた場合には、計測点Ｐの数が条件を満たしていなくとも、これを以後の処理のために選択する（ステップＳ５０８）。
【００６０】
ＣＰＵ３０は、以上のようにして選択した計測点Ｐを対象横断面Ｓに投影する。つまり、各計測点Ｐを対象横断面Ｓに垂直に移動し、対象横断面Ｓ上に位置させる。対象横断面Ｓ上に、多数の計測点Ｐがないと、ブレークラインを得るための変化点を得ることができないからである。
【００６１】
この投影処理の詳細は以下のとおりである。各計測点や左基準位置（横断面Ｓの左端など）５０Ｌ、右基準位置（横断面Ｓの右端など）５０Ｒについて、高さ方向を考慮せず、図９Ａに示すような平面を想定する。ＣＰＵ３０は、左基準位置５０Ｌと右基準位置５０Ｒを結ぶ横断側線ＳＬ（対象横断面Ｓの高さがないもの）に対し、対象とする計測点Ｐから垂線５４を引く。垂線５４と横断側線ＳＬとの交点Ｋを求める（ステップＳ５１０）。
【００６２】
ＣＰＵ３０は、左基準位置５０Ｌから交点Ｋまでの距離を算出し、これをｘ値とする（ステップＳ５１１）。また、計測点Ｐの高さ方向の値をそのままｚ値とする（ステップＳ３１２）。このようにして、計測点Ｐは、図９に示すように、対象横断面Ｓ（水平方向にｘ軸、高さ方向にｚ軸）の上に投影される（ステップＳ３１３）
ＣＰＵ３０は、この投影処理をステップＳ３０８において取得した全ての計測点Ｐについて行う（ステップＳ５０９、Ｓ５１４）。このようにして、当該区間ΔＬにおける対象横断面Ｓ上の点群データを得ることができる。得られた必要点群データの例を、図１０Ａに示す。
【００６３】
ＣＰＵ３０は、以上の点群データを抽出する処理を、全ての区間について実行する（ステップＳ３０３、Ｓ３１４）。したがって、図１０Ｂに示すように、各区間における点群データを得ることができる。
【００６４】
3.2計測点の選択
次に、ＣＰＵ３０は、横断面に投影した計測点Ｐについて、変化点を見いだすために用いる計測点を選択する（図３、ステップＳ６）。計測点選択のフローチャートを、図１１に示す。
【００６５】
まず、ＣＰＵ３０は、当該横断面に投影した計測点Ｐのうち、最大のｚ値（標高）をもつものと、最小のｚ値（標高）を持つものを選択する（ステップＳ６０１）。次に、この最大計測点と最小計測点とを含むように、Ｚ方向に所定区分に分割する。この状態を示したのが、図１２Ａである。図においては、区分Ｅ1から区分Ｅ5までの５つの区分に分割されている。
【００６６】
次に、ＣＰＵ３０は、最も下の区分Ｅ1に含まれる計測点Ｐを選択する（ステップＳ６０４）。図１２Ｂに示すように、計測点Ｐ1とＰ2が選択される。ＣＰＵ３０は、選択した計測点Ｐ1、Ｐ2のＸ方向距離Ｆが、しきい値（たとえば５０cm）より小さいかどうかを判断する（ステップＳ６０５）。しきい値より小さければ、計測点Ｐ1と計測点Ｐ2のＸ座標範囲にある計測点を対象外とする（ステップＳ６０６）。図１２Ｂの場合であれば、計測点Ｐ3、Ｐ4、Ｐ5を対象外とする。しきい値より大きければ、ステップＳ６０６は実行しない。
【００６７】
近接したＸ距離に計測点がある場合、そのＸ範囲であって、それより上の区分にある計測点は、計測誤差である可能性が高いからである。つまり、短い距離の間に、標高が乱高下することは考えにくく、レーザ測距であれ写真測距であれ、地表を計測する場合、地表上の草木などによって、上部方向に誤差を生じる可能性が高いからである。
【００６８】
一方、２点間の距離が離れている場合には、両計測点のＸ範囲における上の区分の計測点は、実際に地表の標高が高くなっている可能性もあることから、対象外としないようにしている。
【００６９】
次に、ＣＰＵ３０は、一つ上の区分Ｅ2について同様の処理を行う。図１３Ａに示すように、対象外とした計測点Ｐ3、Ｐ4、Ｐ5を除外して処理を行う。したがって、図１３Ｂに示すように、計測点Ｐ6が新たに選択され、計測点Ｐ7が対象外とされる。
【００７０】
ＣＰＵ３０は、上記の処理を全ての区分について実行する。これにより、当該横断面について、所望の計測点が選択（選択計測点）されることになる。
【００７１】
次に、ＣＰＵ３０は、選択された計測点のそれぞれにつきメディアンフィルタ処理を行う（ステップＳ６０８）。対象とする計測点を含む前後数個（この実施形態では２個）のｚ値の平均をとり、これを対象とする計測点のｚ値とする。これにより、急激な変化（ノイズ）が取り除かれる。
【００７２】
3.3角度変化に基づく平坦面の抽出
次に、ＣＰＵ３０は、当該横断面について選択された計測点に基づき、計測点を接続する直線の角度変化に基づき、水平面を抽出する（ステップＳ７）。水平面抽出のフローチャートを図１４に示す。ここでは、水平に対して所定角度にある面を水平面と呼んでいる。
【００７３】
ＣＰＵ３０は、まず、横断面の左端の選択計測点を対象選択計測点とし、右隣の選択計測点との角度を算出する。ＣＰＵ３０は、この角度が所定の範囲内（たとえば、−５度から５度）に入っていれば、当該対象選択計測点と隣接する選択計測点は、水平部分の一部であると判断する。次に、隣接する選択計測点を対象選択計測点として、上記の処理を繰り返す。この処理を、右端の選択計測点まで行えば、水平部分を抽出することができる（ステップＳ７０１）。
【００７４】
なお、この水平部分の抽出処理は、各計測点を、所定の標高の範囲（下標高から上標高までの範囲（たとえば３０ｃｍの範囲））に区分して行う。したがって、各区分において、平坦部分があればこれが選択される。
【００７５】
図１５Ａに示すような場合には、異なる区分において水平部分Ｈ1とＨ2が抽出されることになる。なお、一つの区分において、水平部分が複数見いだされた場合には、誤差である可能性が高いので、低い方の水平部分を、高い方の水平部分にあわせる。つまり、低い方の水平部分の標高を上げて、一つの水平部分とする。
【００７６】
次に、ＣＰＵ３０は、抽出された水平部分のうち、所定以上の幅（たとえば、１ｍ以上）があって、最も高い区分にあるものを見いだし、これを天端面とする（ステップＳ７０２）。図１５Ａの場合であれば、水平部分Ｈ2が天端面として選択されることになる。ここで、天端面とは、河川の両岸にある堤防の一番高い水平部分をいうものである。なお、所定の幅以上のものとしたのは、天端面より上に木などがある場合、この部分を天端面とする誤判断を避けるためである。
【００７７】
次に、ＣＰＵ３０は、見いだした天端面より高い部分を削除する（ステップＳ７０３）。図１５Ｂの場合であれば、選択計測点ＤＰ（×印を付した点）を削除する。これらは、天端面に植えられた植物や木などである可能性が高いからである。
【００７８】
次に、ＣＰＵ３０は、図１５Ｂに示すように、天端面Ｈ2と隣接する斜線ＣＬを延長し、両者の交点を求め、これを変化点Ｖとする（ステップＳ７０４）。このようにして、左端の変化点Ｖと右端の変化点Ｖとの間が天端面であるとして抽出することができる。
【００７９】
以上のように、図３のステップＳ５〜Ｓ７によって、横断面における天端面（正確には天端面を構成する点端線）を抽出することができる。ＣＰＵ３０は、この処理を、すべての横断面について繰り返し、すべての横断面について天端線を抽出する（ステップＳ４、ステップＳ８）。
【００８０】
続いて、ＣＰＵ３０は、図１６に示すように、各横断面の天端線Ｕを結んで、破線にて示す天端面ＵＦを生成する（ステップＳ９）。さらに、ＣＰＵ３０は、変化点Ｖを接続して、ブレークライン５１を生成する（ステップＳ１０）。これらは、対象物の形状の特徴を表すものとして用いられる。
【００８１】
以上のようにして、平坦面（水平面）Ｈ1、Ｈ2を抽出することができる。
【００８２】
４．その他
(1)上記実施形態では、対象物の形状の特徴を表すものとして、天端面やブレークラインを抽出している。しかし、いずれか一方のみを抽出するようにしてもよいし、他の特徴を抽出するようにしてもよい。
【００８３】
(2)上記実施形態では、河川を対象として説明した。しかし、道路など他の対象物にも適用することができる。道路の場合には、天端面に代えて、道路面を抽出することができる。
【００８４】
(3)上記実施形態では、最も高い水辺部分を天端面として抽出している。しかし、候補ブレークラインの高さに近い水平部分を天端面として抽出するようにしてもよい。
【００８５】
(4)上記実施形態において算出した選択計測点を、格子状の水平位置を有する測量データに変換するようにしてもよい。選択計測点の水平位置（Ｘ、Ｙ座標値）は、ランダムである。このため、これら選択計測点に基づいて任意の位置で断面形状を得ようとすると、整った形状を得ることが難しい。そこで、図１７Ａに示すように、ＸＹ座標に関し等間隔の格子となる交点において、Ｚ値を有する点ＣＰ、ＣＰ・・・のデータに変換する。
【００８６】
この場合、図１７Ｂに示すように、点ＣＰのＺ値は、近接する複数の計測点ＰのＺ値を平均したものとすればよい。なお、この平均値の算出において、距離ｄが小さい計測点ＰのＺ値ほど、ウエイトが高くなるように加重平均を求めるようにすることが好ましい。
【００８７】
さらに、このような格子状の点データに加えて、図１７Ａに示すように、ブレークライン５１上の点ＢＰ、ＢＰ・・・についても、点データを生成してもよい。点ＢＰのＺ値の算出方法は、図１７Ｂと同じ方法とすることができる。なお、この実施形態では、格子状の点ＣＰの間隔よりも、ブレークライン５１上の点ＢＰの間隔を密にしている。これにより、これら点データに基づいて任意位置での断面図などを生成する場合に、ブレークライン５１が明瞭に現れるようになる。
【００８８】
(5)上記実施形態では、必要点群データ抽出処理において、計測点の数によって計測点を選択するようにしている（ステップＳ５０４、Ｓ５０５）。しかし、計測点の数に代えて、あるいは、計測点の数に加えて、次の条件を選択のための指標としてもよい。
【００８９】
たとえば、図８Ｂの所定区間ΔＬに含まれる計測点のＸ値のばらつき（標準偏差）が所定値以上あることを選択条件としてもよい。標準偏差が所定値に達していなければ、対象幅を広げ、選択条件を満足するようにする。
【００９０】
あるいは、所定区間ΔＬをさらに細かい微小区間に区分し、各微小区間において計測点があるか否かを判断し、計測点がない微小区間が所定数以上連続していないことを選択条件としてもよい。計測点がない微小区間が所定数以上連続していれば、対象幅を広げ、選択条件を満足するようにする。
【００９１】
(6)上記実施形態では、ＰＣを測量データ処理装置としているが、サーバ装置によってこれを実現するよいうにしてもよい。  
【図１】
0006514901.tif 000002
【図２】
0006514901.tif 000003
【図３】
0006514901.tif 000004
【図４】
0006514901.tif 000005
【図５】
0006514901.tif 000006
【図６】
0006514901.tif 000007
【図７】
0006514901.tif 000008
【図８】
0006514901.tif 000009
【図９】
0006514901.tif 000010
【図１０】
0006514901.tif 000011
【図１１】
0006514901.tif 000012
【図１２】
0006514901.tif 000013
【図１３】
0006514901.tif 000014
【図１４】
0006514901.tif 000015
【図１５】
0006514901.tif 000016
【図１６】
0006514901.tif 000017
【図１７】
0006514901.tif 000018