(57)【要約】【課題】簡易な方式で特徴的な地形を検出することが可能な地形状態検出装置を提供する。【解決手段】地形状態検出装置は、３次元計測情報を取得する取得部と、取得部で取得した３次元計測情報に基づいて、比較対象となる第１範囲の第１傾斜量を算出し、第１範囲よりも狭い第２範囲の第２傾斜量を算出し、第１傾斜量と第２傾斜量との差分に基づいて、所定の地形状態を検出する地形状態検出部とを備える。【選択図】図１
2019158530.tif 000002
【特許請求の範囲】
【請求項１】
３次元計測情報を取得する取得部と、
前記取得部で取得した３次元計測情報に基づいて、比較対象となる第１範囲の第１傾斜量を算出し、前記第１範囲よりも狭い第２範囲の第２傾斜量を算出し、前記第１傾斜量と前記第２傾斜量との差分に基づいて、所定の地形状態を検出する地形状態検出部とを備える、地形状態検出装置。
【請求項２】
前記３次元計測情報に基づいて複数地点のグリッドデータを生成するグリッドデータ生成部をさらに備え、
前記地形状態検出部は、
前記第１範囲に含まれる隣接する複数地点のグリッドデータに基づいて第１傾斜量を算出し、
前記第２範囲に含まれる隣接する複数地点のグリッドデータに基づいて第２傾斜量を算出する、請求項１記載の地形状態検出装置。
【請求項３】
前記地形状態検出部は、
前記第１傾斜量と前記第２傾斜量との差分が所定の閾値を超えるか否かを判断し、
前記差分が所定の閾値を超える場合には、亀裂を含む地形があると検出する、請求項１記載の地形状態検出装置。
【請求項４】
３次元計測情報を取得するステップと、
取得した３次元計測情報に基づいて、比較対象となる第１範囲の第１傾斜量を算出するステップと、
取得した３次元計測情報に基づいて、前記第１範囲よりも狭い第２範囲の第２傾斜量を算出するステップと、
前記第１傾斜量と前記第２傾斜量との差分に基づいて、所定の地形状態を検出するステップとを備える、地形状態検出方法。
【請求項５】
コンピュータを、地形状態検出装置として機能させるためのプログラムであって、
前記プログラムは、前記コンピュータに、
３次元計測情報を取得する取得部と、
前記取得部で取得した３次元計測情報に基づいて、比較対象となる第１範囲の第１傾斜量を算出し、前記第１範囲よりも狭い第２範囲の第２傾斜量を算出し、前記第１傾斜量と前記第２傾斜量との差分に基づいて、所定の地形状態を検出する地形状態検出部として、機能させる地形状態検出プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、地形状態検出装置に関し、斜面の亀裂を含む地形等の特徴的な地形状態を検出する方式に関する。
【背景技術】
【０００２】
豪雨などの自然災害によって、地滑りなどの斜面崩壊が発生すると、交通路が土砂で塞がれたり、道路や鉄道盛土が変形したりすることによって、交通網が麻痺する被害を蒙る。さらには、崩壊した土砂が土石流となって、下流側の住民などに大きな被害をもたらす場合もある。このような被害を未然に防ぐため、あるいは、被害を最小限に抑えるためには、斜面崩壊のおそれがある場所を検出することが重要である。
【０００３】
この点で、空中写真等を用いて画像処理により特徴的な地形を検出する方式が提案されている（特許文献１）。
【０００４】
この場合、撮影条件により得られる空中写真が異なるため検出結果が異なる可能性があり、精度の高い検出が難しいという可能性がある。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００３−１４０５４４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本開示は、上記の問題点を解決するためになされたものであり、簡易な方式で特徴的な地形を検出することが可能な地形状態検出装置、地形状態検出方法および地形状態検出プログラムを提供することを目的とする。
【課題を解決するための手段】
【０００７】
ある実施形態に地形状態検出装置は、３次元計測情報を取得する取得部と、取得部で取得した３次元計測情報に基づいて、比較対象となる第１範囲の第１傾斜量を算出し、第１範囲よりも狭い第２範囲の第２傾斜量を算出し、第１傾斜量と第２傾斜量との差分に基づいて、所定の地形状態を検出する地形状態検出部とを備える。
【０００８】
好ましくは、３次元計測情報に基づいて複数地点のグリッドデータを生成するグリッドデータ生成部をさらに備える。地形状態検出部は、第１範囲に含まれる隣接する複数地点のグリッドデータに基づいて第１傾斜量を算出し、第２範囲に含まれる隣接する複数地点のグリッドデータに基づいて第２傾斜量を算出する。
【０００９】
好ましくは、地形状態検出部は、第１傾斜量と第２傾斜量との差分が所定の閾値を超えるか否かを判断し、差分が所定の閾値を超える場合には、亀裂を含む地形があると検出する。
【００１０】
ある実施形態に従う地形状態検出方法は、３次元計測情報を取得するステップと、取得した３次元計測情報に基づいて、比較対象となる第１範囲の第１傾斜量を算出するステップと、取得した３次元計測情報に基づいて、第１範囲よりも狭い第２範囲の第２傾斜量を算出するステップと、第１傾斜量と第２傾斜量との差分に基づいて、所定の地形状態を検出するステップとを備える。
【００１１】
ある実施形態に従う地形状態検出プログラムは、コンピュータを、地形状態検出装置として機能させるためのプログラムであって、プログラムは、コンピュータに、３次元計測情報を取得する取得部と、取得部で取得した３次元計測情報に基づいて、比較対象となる第１範囲の第１傾斜量を算出し、第１範囲よりも狭い第２範囲の第２傾斜量を算出し、第１傾斜量と第２傾斜量との差分に基づいて、所定の地形状態を検出する地形状態検出部として、機能させる。
【発明の効果】
【００１２】
本開示の地形状態検出装置、地形状態検出方法および地形状態検出プログラムは、簡易な方式で特徴的な地形を検出することが可能である。
【図面の簡単な説明】
【００１３】
【図１】実施形態に従う地形状態検出装置１０について説明する図である。
【図２】実施形態に基づくグリッドデータ生成部２４の処理について説明する図である。
【図３】実施形態に基づくグリッドデータに基づいて地形の標高を陰影表現で描画した場合を説明する図である。
【図４】実施形態に基づく地形状態検出装置１０の処理について説明するフロー図である。
【図５】実施形態に基づく傾斜量を算出する範囲を説明する図である。
【図６】実施形態に基づく任意のセル座標と隣接するセルとの位置関係を説明する図である。
【図７】実施形態に基づいて所定エリアにおける第１範囲に基づいて算出した第１傾斜量の傾斜区分図である。
【図８】実施形態に基づいて所定エリアにおける第２範囲に基づいて算出した第２傾斜量の傾斜区分図である。
【図９】実施形態に基づく亀裂を含む地形の検出の方式を説明する概念図である。
【図１０】実施形態に基づく平面モデルで傾斜量の差分を説明する概念図である。
【図１１】実施形態に基づく差分のイメージを説明する概念図である。
【図１２】実施形態に基づく所定エリアの亀裂を含む地形を描画した図である。
【発明を実施するための形態】
【００１４】
実施形態について図面を参照しながら詳細に説明する。なお、図中同一または相当部分には同一符号を付し、その説明は繰り返さない。
【００１５】
図１は、実施形態に従う地形状態検出装置１０について説明する図である。
図１を参照して、地形状態検出装置１０は、ＣＰＵ（Central Processing Unit）２０と、表示部１４と、メモリ１６と、操作部１８とを含む。
【００１６】
ＣＰＵ２０は、地形状態検出装置１０全体を制御する。
ＣＰＵ２０は、メモリ１６に格納されているプログラムを実行することにより種々の機能ブロックを実現する。
【００１７】
具体的には、ＣＰＵ２０は、データ取得部２２と、グリッドデータ生成部２４と、地形状態算出部２５とを含む。
【００１８】
地形状態算出部２５は、セル設定部２７と、第１算出部２６と、第２算出部２８と、検出部３０とを含む。
【００１９】
メモリ１６は、コンピュータに内蔵されるＲＯＭ(Read Only Memory)、ＲＡＭ(Random Access Memory)、ハードディスクなどの記憶装置である。メモリ１６は、データ取得部２２、グリッドデータ生成部２４、および地形状態算出部２５の処理プログラムを含む各種プログラムや各種データを記憶し、ＣＰＵ２０との間でこれらの情報を入出力する。
【００２０】
表示部１４は液晶モニタ等の画像表示装置であり、操作部１８はキーボードやマウスなどで構成される。
【００２１】
地形状態検出装置１０は、取得した３次元計測情報に基づいて特徴的な地形状態を検出する。
【００２２】
本例においては、航空レーザ測量で計測した３次元計測データを用いる場合について説明する。
【００２３】
データ取得部２２は、航空レーザ測量の計測結果である３次元計測データに基づいてグラウンドデータを取得する。
【００２４】
なお、データ取得部２２は、航空レーザ測量で計測される３次元計測データのオリジナルデータである数値表層モデル、ＤＳＭ(Digital Surface Model)データを取得して、建物や樹木の高さの情報を除くフィルタリング処理してグラウンドデータを取得するようにしても良い。
【００２５】
図２は、実施形態に基づくグリッドデータ生成部２４の処理について説明する図である。
【００２６】
図２に示されるように、グリッドデータ生成部２４は、グラウンドデータから内挿補間により格子状の標高データ（グリッドデータ）を作成する。
【００２７】
なお、内挿補間の方式としては、ＴＩＮ(Triangulated irregular Network)、最近隣法等を用いることが可能である。あるいは、Ｋｒｉｎｇｉｎｇ法を用いることも可能である。
【００２８】
図３は、実施形態に基づくグリッドデータに基づいて地形の標高を陰影表現で描画した場合を説明する図である。
【００２９】
図３に示されるように、一例として、地形の標高が高いところが色が濃くなり、標高が低いところが色が薄く示されている。本例においては、一例として当該所定エリアにおける特徴的な地形を検出する方式について説明する。特徴的な地形として亀裂を含む地形を検出する場合について説明する。
【００３０】
図４は、実施形態に基づく地形状態検出装置１０の処理について説明するフロー図である。
【００３１】
図４を参照して、ＣＰＵ２０は、３次元計測データを取得する（ステップＳ２）。具体的には、データ取得部２２は、グラウンドデータを取得する。当該グラウンドデータは、図示しない外部装置から取得する。たとえば、航空測量計測装置から取得するようにしても良い。あるいは、外部のサーバからグラウンドデータをダウンロードして取得するようにしても良い。
【００３２】
次に、ＣＰＵ２０は、グリッドデータを生成する（ステップＳ４）。グリッドデータ生成部２４は、取得したグランドデータに基づいて図２で説明したグリッドデータを生成する。
【００３３】
次に、ＣＰＵ２０は、地形量算出セルを設定する（ステップＳ５）。
セル設定部２７は、地形状態を検出する所定の領域を複数のセル状に分割する。
【００３４】
次に、ＣＰＵ２０は、第１傾斜量を算出する（ステップＳ６）。
第１算出部２６は、第１範囲に含まれる複数地点のグリッドデータに基づいて第１傾斜量を算出する。
【００３５】
次に、ＣＰＵ２０は、第２傾斜量を算出する（ステップＳ８）。
第２算出部２８は、第２範囲に含まれる複数地点のグリッドデータに基づいて第２傾斜量を算出する。
【００３６】
次に、ＣＰＵ２０は、差分量を算出する（ステップＳ１０）。検出部３０は、第１算出部２６で算出した第１傾斜量と、第２算出部２８で算出した第２傾斜量との差分量を算出する。
【００３７】
次に、ＣＰＵ２０は、差分量が所定の閾値以上か否かを判断する（ステップＳ１２）。検出部３０は、算出した差分量が所定の閾値以上か否かを判断する。
【００３８】
次に、ＣＰＵ２０は、ステップＳ１２において、算出した差分量が所定の閾値以上であると判断した場合（ステップＳ１２においてＹＥＳ）には、亀裂を含む地形有りと判定する（ステップＳ１４）。検出部３０は、第１算出部２６で算出した第１傾斜量と、第２算出部２８で算出した第２傾斜量との差分量が所定の閾値以上であると判断した場合には、亀裂を含む地形有りと判定する。
【００３９】
そして、処理を終了する（エンド）。
一方、ＣＰＵ２０は、ステップＳ１２において、算出した差分量は所定の閾値以上でないと判断した場合（ステップＳ１２においてＮＯ）には、亀裂を含む地形無しと判定する（ステップＳ１６）。検出部３０は、第１算出部２６で算出した第１傾斜量と、第２算出部２８で算出した第２傾斜量との差分量が所定の閾値以上でないと判断した場合には、亀裂を含む地形無しと判定する。
【００４０】
そして、処理を終了する（エンド）。
図５は、実施形態に基づく地形量を算出するセルおよび範囲について説明する図である。
【００４１】
図５を参照して、本例においては、セル設定部２７により所定の領域が複数のセル状に分割された場合が示されている。セル間隔は、一例として２ｍ×２ｍに設定されている。なお、当該セル間隔は、任意の値に設定することが可能である。
【００４２】
各セルには、それぞれグリッドデータが対応付けられている。なお、１つのセルの中に複数のグリッドデータが含まれる場合には、その代表的なグリッドデータを用いても良いし、平均化したグリッドデータを用いるようにしても良い。
【００４３】
実施形態においては、複数地点（本例においては複数セル）のグリッドデータに基づいて傾斜量を算出する。
【００４４】
具体的には、地形量を算出するセル（地形量算出セル）に対して隣接するセルのグリッドデータを用いて傾斜量を算出する。
【００４５】
本例においては、第１範囲として広範囲の領域を設定する。第２範囲として狭範囲の領域を設定する。具体的には、第１範囲として、１０ｍ離れた隣接するセルのグリッドデータを用いて第１傾斜量を算出する。第２範囲として、２ｍ離れた隣接するセルのグリッドデータを用いて第２傾斜量を算出する。なお、当該数値は一例であり当該値に限定されるものではない。
【００４６】
図６は、実施形態に基づく任意のセル座標と隣接するセルとの位置関係を説明する図である。
【００４７】
図６には、地形量を算出するセルの座標（xi，yi）が示されている。また、地形量を算出するセルに対して隣接するセルの座標が示されている。
【００４８】
本例においては、地形量を算出するセルに対して、隣接する上下左右のセルのグリッドデータに基づいて傾斜量を算出する。
【００４９】
具体的には、傾斜量は、地形量を算出するセルのグリッドデータｕ（xi，yi）に隣接するセルのグリッドデータｕ（xi−１，yi）、ｕ（xi＋１，yi）、ｕ（xi，yi−１）、ｕ（xi，yi＋１）を用いて次式に基づいて算出することが可能である。
【００５０】
【数１】
2019158530.tif 000003
【００５１】
図７は、実施形態に基づいて所定エリアにおける第１範囲に基づいて算出した第１傾斜量の傾斜区分図である。
【００５２】
図８は、実施形態に基づいて所定エリアにおける第２範囲に基づいて算出した第２傾斜量の傾斜区分図である。
【００５３】
図７および図８に示されるように傾斜量が大きい箇所が濃くなり、傾斜量が小さい箇所が薄く示されている。
【００５４】
また、図７の方が粗い傾斜量の変化が示されており、図８の方が密な傾斜量の変化が示されている。
【００５５】
図９は、実施形態に基づく亀裂を含む地形の検出の方式を説明する概念図である。
図９に示されるように、第１範囲として広範囲の領域を設定して第１傾斜量を算出する場合には地形量を算出するセルの周辺に隣接するグリッドデータに基づいて平滑化した基準平面に対する傾斜量が算出される。
【００５６】
一方、第２範囲として狭範囲の領域を設定して第２傾斜量を算出する場合には、局所的な地形に対する傾斜量が算出される。
【００５７】
したがって、第１傾斜量と第２傾斜量とを比較すると、亀裂を含む地形は、第１傾斜量と第２傾斜量との差分の値が大きくなる。
【００５８】
一方で、第１傾斜量と第２傾斜量とを比較すると、亀裂無しの地形は、第１傾斜量と第２傾斜量との差分の値は小さくなる。
【００５９】
本例の場合には、差分の値が所定の閾値以上である場合に亀裂を含む地形有りと判定する。
【００６０】
図１０は、実施形態に基づく平面モデルで傾斜量の差分を説明する概念図である。
図１０に示されるように、広範囲に基づく数値地形モデルは、地形を全体的に平滑化した基準平面の傾斜量として算出される。
【００６１】
狭範囲に基づく数値地形モデルは、局所平面の精度の高い傾斜量が算出される。
図１１は、実施形態に基づく差分のイメージを説明する概念図である。
【００６２】
図１１に示されるように、全て０の数値モデルは、広範囲に基づく基準平面の傾斜量である。５が含まれる数値モデルは、狭範囲に基づく局所平面の傾斜量である。
【００６３】
その差分により局所的な傾斜量の変化が表面化する。
これにより亀裂を含む地形を判定することが可能である。
【００６４】
図１２は、実施形態に基づく所定エリアの亀裂を含む地形を描画した図である。
図１２に示されるように、亀裂を含む地形のみ抽出することが可能である。
【００６５】
これにより、簡易な方式で斜面の亀裂を含む地形を判定することが可能であり、亀裂を含む地形を精度よく判定することが可能である。
【００６６】
当該亀裂を含む地形の判定により深層崩壊が生じ易い地形形状を簡易に把握することが可能である。
【００６７】
また、本実施形態におけるプログラムとして、パーソナルコンピュータで実行可能なアプリケーションを提供してもよい。このとき、本実施の形態に係るプログラムは、パーソナルコンピュータ上で実行される各種アプリケーションの一部の機能として組み込まれてもよい。
【００６８】
今回開示された実施の形態はすべての点で例示であって制限的なものではないと考えられるべきである。本開示の範囲は、上記した説明ではなく、特許請求の範囲によって示され、特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれることが意図される。
【符号の説明】
【００６９】
１０  地形状態検出装置、１４  表示部、１６  メモリ、１８  操作部、２２  データ取得部、２４  グリッドデータ生成部、２５  地形状態算出部、２６  第１算出部、２８  第２算出部、３０  検出部。
【図１】
2019158530.tif 000004
【図２】
2019158530.tif 000005
【図３】
2019158530.tif 000006
【図４】
2019158530.tif 000007
【図５】
2019158530.tif 000008
【図６】
2019158530.tif 000009
【図７】
2019158530.tif 000010
【図８】
2019158530.tif 000011
【図９】
2019158530.tif 000012
【図１０】
2019158530.tif 000013
【図１１】
2019158530.tif 000014
【図１２】
2019158530.tif 000015