JP2008530595A - 地形モデル化のためのデジタル標高モデル(dem)を向上させる方法および装置 - Google Patents

地形モデル化のためのデジタル標高モデル(dem)を向上させる方法および装置 Download PDF

Info

Publication number
JP2008530595A
JP2008530595A JP2007554275A JP2007554275A JP2008530595A JP 2008530595 A JP2008530595 A JP 2008530595A JP 2007554275 A JP2007554275 A JP 2007554275A JP 2007554275 A JP2007554275 A JP 2007554275A JP 2008530595 A JP2008530595 A JP 2008530595A
Authority
JP
Japan
Prior art keywords
dem
refined
original
null
computer screen
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.)
Withdrawn
Application number
JP2007554275A
Other languages
English (en)
Inventor
ラームズ,マーク
カープ,ジョン
スミス,アンソニー
コネッティー,スティーヴン,ジュニア
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.)
Harris Corp
Original Assignee
Harris 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 Harris Corp filed Critical Harris Corp
Publication of JP2008530595A publication Critical patent/JP2008530595A/ja
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05Geographic models
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C11/00Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/88Lidar systems specially adapted for specific applications
    • G01S17/89Lidar systems specially adapted for specific applications for mapping or imaging
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/30Noise filtering

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Software Systems (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Multimedia (AREA)
  • Computer Graphics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Electromagnetism (AREA)
  • Processing Or Creating Images (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Instructional Devices (AREA)

Abstract

土地(52)およびその上の複数のオブジェクト(54、56)についてのデータを含む独自のデジタル標高モデル(DEM)を処理するためのコンピュータ実装される方法。本方法は、元のDEM(102)から低解像度DEM(104)を生成し、低解像度DEM(104)と元のDEM(102)との比較に基づいてオブジェクトのみのDEM(108)を生成することを含む。本方法はさらに、オブジェクトのみDEM(114、116)におけるノイズを削減することを含む。

Description

地理的領域の地形モデルは、フライト・シミュレータおよび氾濫原解析を含め多くの用途に使用されうる。さらに、人工構造(たとえば都市)の地形モデルは、たとえば携帯電話のアンテナ配置、都市計画、災害の備えおよび解析ならびに地図作成といった用途においてきわめて役に立つことがありうる。
現在、地形モデルを作成するためのさまざまな種類および方法が使われている。一つの一般的な地形モデルはデジタル標高モデル(DEM: digital elevation model)である。DEMは、ある地理的領域のサンプリングされたマトリクス表現であり、コンピュータによって自動化された方法で生成されうる。DEMでは、座標点に高さの値が対応付けられる。DEMは典型的には、ある標高から次の標高にかけて、異なる標高(たとえば谷、山など)の間の遷移が概してなめらかである土地をモデル化するために使われる。すなわち、DEMは典型的には土地を複数の曲面としてモデル化し、その間に不連続があったとしても「ならされ」てしまう。この理由で、DEMは一般に、都市中心域の高層ビルのような人工構造を、上記の用途の多くのために十分な精度でモデル化するには好適でない。
Rahmes et al.に対する米国特許第6,654,690号は地形学における著しい進歩を開示している。この'690特許は、土地およびその上の建物を含む領域の地形モデルを、標高対位置のランダムに離間したデータに基づいて作成する自動化された方法を開示している。'690特許は、本発明の被譲渡者に譲渡されており、ここに参照によってその全体において組み込まれる。その方法は、前記のランダムに離間したデータを処理して所定の位置グリッドに合う標高対位置のグリッド化データを生成し、該グリッド化データを処理して土地データから建物データを区別し、建物データについてポリゴン抽出を実行してその領域の土地およびその上の建物を含む地形モデルを作成することを含む。
さらに、DEMの品質は土地のモデル化ならびに該土地上の建物および植生のモデル化に影響する。上述のように、'690特許は、初期DEMが、土地のみのDEMと建物のみのDEMに分離されることを開示する。これはモデル化前の各個別DEMの品質向上の役に立つが、表示される地形モデルがよりリアルになりうるよう、個別DEMの向上をさらに改善する必要がある。
以上の背景に鑑み、デジタル標高モデル(DEM)を、該DEM内の建物、植生および土地をモデル化するのに先立ってさらに向上させることが本発明の一つの目的である。
本発明に基づくこの目的およびその他の目的、特徴および効果は、土地およびその上の複数のオブジェクトを含む元のデジタル標高モデル(DEM)を処理するためのコンピュータ実装される方法によって提供される。該方法は、元のDEMからより低い解像度の低解像度DEMを生成し、該低解像度DEMと元のDEMとの比較に基づいてオブジェクトのみのDEMを生成し、該オブジェクトのみのDEMにおけるノイズを削減することを含みうる。
オブジェクトのみのDEMにおけるノイズ削減は、第一のウィンドウ・サイズに基づいてオブジェクトのみのDEMにおけるヌル(null)を同定し、オブジェクトのみのDEM中に同定されたヌルを拡張し、オブジェクトのみのDEMの拡張されたヌルを埋めて第一の洗練されたオブジェクトのみDEMを生成することを含みうる。
本方法はさらに、前記第一のウィンドウ・サイズに基づいて前記第一の洗練されたオブジェクトのみのDEMにおけるヌルを同定し、前記第一の洗練されたオブジェクトのみのDEM中に同定されたヌルを拡張し、前記第一の洗練されたオブジェクトのみのDEMの拡張されたヌルを埋めて第二の洗練されたオブジェクトのみDEMを生成することを含みうる。
前記第一の洗練されたオブジェクトのみDEMと前記第二の洗練されたオブジェクトのみDEMとの比較に基づいて、ノイズのみのDEMが生成されうる。該ノイズのみのDEMが、前記第二の洗練されたオブジェクトのみDEMから引き算されて、最終的な洗練されたオブジェクトのみDEMが生成される。最終的な洗練されたオブジェクトのみDEMは、従来技術で地形モデル化のために典型的に使われる当初のオブジェクトのみDEMよりも著しく少ないノイズを有しうる。これは、よりリアルな地形モデル化につながる。
オブジェクトのみDEMを向上させるもう一つの実施形態は次のようなものである。低解像度DEMが元のDEMから生成されることができ、該低解像度DEM中のヌルが第一のウィンドウ・サイズに基づいて同定され、該ヌルが埋められる。次いで、そのヌルが埋められた低解像度DEMと元のDEMとの比較に基づいて、オブジェクトのみDEMが生成されうる。該オブジェクトのみDEMと元のDEMとの比較に基づいて、該オブジェクトのない元のDEMが生成されうる。オブジェクトのない元のDEMにおけるヌルが第一のウィンドウ・サイズに基づいて同定され、ヌルが拡張される。
オブジェクトのない元のDEMからオブジェクトのない低解像度DEMが生成され、該オブジェクトのない低解像度DEM中のヌルが第一のウィンドウ・サイズに基づいて同定され、該ヌルが埋められて土地のみのDEMが生成される。すると、該土地のみDEMと元のDEMとの比較に基づいて、洗練されたオブジェクトのみDEMが生成されうる。該洗練されたオブジェクトのみDEMは、従来技術で地形モデル化のために典型的に使われる当初のオブジェクトのみDEMより、著しくノイズが少ないことがありうる。
洗練された土地のみDEMも生成されうる。これは、前記の洗練されたオブジェクトのみDEMと元のDEMとの比較に基づいて洗練された該オブジェクトのないDEMを生成し、該洗練されたオブジェクトのないDEMから低解像度のオブジェクトのないDEMを生成し、前記第一のウィンドウ・サイズに基づいて前記低解像度のオブジェクトのないDEM中にヌルを同定することを含む。すると該ヌルが埋められて、前記した洗練された土地のみDEMが生成される。さらに、前記の洗練された土地のみDEMと元のDEMとの比較に基づいて、最終的な洗練されたオブジェクトのみDEMが生成されうる。
本発明に基づくもう一つの側面は、土地およびその上の複数のオブジェクトについてのデータを含むデジタル標高モデル(DEM)を処理するためのデータ構造が記憶されているコンピュータ可読媒体に向けられる。該コンピュータ可読媒体の第一の実施形態は、元のDEMからの低解像度DEMを生成するためのデータを含む第一のデータ・フィールドと、前記低解像度DEMと元のDEMとの比較に基づいてオブジェクトのみDEMを生成するためのデータを含む第二のデータ・フィールドと、前記オブジェクトのみDEMのノイズを削減するためのデータを含む第三のデータ・フィールドとを含みうる。
前記コンピュータ可読媒体の第二の実施形態は、元のDEMからの低解像度DEMを生成するためのデータを含む第一のデータ・フィールドと、第一のウィンドウ・サイズに基づいて前記低解像度DEM中のヌルを同定し、該ヌルを埋めるためのデータを含む第二のデータ・フィールドと、ヌルを埋めた前記低解像度DEMと元のDEMとの比較に基づいてオブジェクトのみDEMを生成するためのデータを含む第三のデータ・フィールドとを含みうる。第四のデータ・フィールドは、前記オブジェクトのみDEMと元のDEMとの比較に基づいて該オブジェクトのない元のDEMを生成するためのデータを含みうる。第五のデータ・フィールドは、前記第一のウィンドウ・サイズに基づくオブジェクトのない元のDEMにおけるヌルを同定し、該ヌルを拡張するためのデータを含みうる。第六のデータ・フィールドは、オブジェクトのない元のDEMからオブジェクトのない低解像度DEMを生成するためのデータを含みうる。第七のデータ・フィールドは、前記第一のウィンドウ・サイズに基づくオブジェクトなしの低解像度DEMにおけるヌルを同定し、土地のみのDEMを生成するために該ヌルを埋めるためのデータを含みうる。第八のデータ・フィールドは、前記土地のみのDEMと元のDEMとの比較に基づいて洗練されたオブジェクトのみDEMを生成するためのデータを含みうる。
本発明に基づくさらにもう一つの側面は、土地およびその上の複数のオブジェクトについてのデータを含む元のデジタル標高モデル(DEM)を処理するプロセッサを有する、地形モデル化のためのコンピュータ・システムに向けられる。前記プロセッサは、上で論じたコンピュータ可読媒体の第一および第二の実施形態を処理しうる。
本発明について、これから付属の図面を参照しつつより完全に記述する。その際、本発明の好ましい実施形態を示すが、本発明は多くの異なる形で具現されうるものであり、ここに述べられる実施形態に限定されると解釈されるべきではない。むしろ、これらの実施形態は本開示を十全かつ完備なものとし、当業者に本発明の範囲を完全に伝達するために与えられるものである。全体を通じて同様の参照符号は同様の要素を指す。
地形データを収集する収集機50および収集された地形データからデジタル標高モデル(DEM)を生成するためのシステム60についてこれから図1を参照しつつ説明する。該DEMは土地52およびその土地の上のオブジェクトを含むある地域のものであり、ここでオブジェクトは建物54および植生56でありうる。植生56は主として樹木を含み、したがって植生および樹木は交換可能である。土地52、建物54および植生56のモデル化は、当該地域でのランダムかつ任意に離間した標高対位置のデータに基づく。
光検出および測距(LIDAR: light detection and ranging)のような収集機50はそのランダムに離間したデータを収集するために使用されうる。ランダムに離間したデータは、名目上は、非一様に離間した位置および高度の測定値の組でありうる。LIDAR収集機50は、都市のような関心のある地域の上を航空機70によって運ばれうる。当該地域は、たとえば道路58のような、建物54および樹木56に比べて比較的小さな特徴をも含んでいてもよい。
当業者は、LIDAR源は、単一の画像からの標高対位置情報を含むデータを提供することを認識するであろう。標高対位置データを与えるには、一般には、種々の視点から撮影された当該地域の複数の光学画像が必要とされるが、その同じ情報を単一のLIDAR画像から得ることができる。もちろん、当業者が認識するであろうように、本発明はLIDAR収集機に加えて、たとえば光学的(たとえば写真)、電気光学的および赤外線の源といった種々の源からの標高対位置データを使うことができる。LIDARデータによって与えられる位置情報は、たとえば緯度および経度の情報を含みうるが、他の好適な位置指標も使われうる。
ひとたびランダムに離間したデータが収集されたら、該データは、コンピュータ62への転送のために、たとえば磁気ディスクまたは光ディスクといった記憶媒体80に記憶されうる。もちろん、当業者がすぐ認識するであろうように、データ転送のための他の好適な方法を使ってもよい。ランダムに離間されたデータは、コンピュータ62によって見るためのDEMを生成するために使われる。
DEMを見るためにディスプレイ64がコンピュータ62に接続されている。キーボード66およびマウス68のような入力装置もコンピュータ62に接続されている。本発明によれば、コンピュータ62は、1)土地のみDEMおよびオブジェクトのみDEMを生成し、次いでオブジェクトのみDEMからノイズを除去することによってDEMを向上させ、2)オブジェクトのみDEMを建物DEMおよび植生DEMに分離するためのプロセッサ69を含む。
元のDEMすなわち当初のDEMの生成および向上について、これから図2の流れ図および図3〜25に示されるコンピュータ表示画面を参照しつつ論じる。流れ図では、ステップ(1)〜(19)についてまず論じる。ここで、ステップ(2)〜(19)は、のちにより詳細に論じるように、バッチ・プロセスの一部と考えられる。バッチ・プロセスでは、流れ図に示されるブロックのいくつかは二度以上論じられる。そのそれぞれの機能がループ・プロセスに基づいて繰り返されるからである。元のDEMを向上させる本方法をよりよく解説するため、論じられる各ステップの番号を、対応するブロック内に括弧に入れて与えてある。
開始(ブロック100)から、初期ステップ(1)はブロック102でコンピュータ62を使って、記憶媒体80を介して提供されるランダムに離間されたデータから元のDEMを生成する。図3に示される初期コンピュータ画面200を参照すると、ユーザーはフィールド202で「点からDEMを生成」を選択する。すると図4に示されるように「点からDEMの設定」コンピュータ画面204が生成される。
「点からDEMの設定」コンピュータ画面204では、収集データを記憶しているファイルの名前がフィールド206に入力される。点の書式がフィールド208で選択される。この場合、点はユニバーサル横メルカトル(UTM)グリッドに基づく。点の測定単位がフィールド210で選択され、たとえばメートルである。UTMグリッドは、60個の南北方向のゾーンを含んでおり、各ゾーンは緯度において6度の幅である。関心のあるUTMゾーンがフィールド212で選択される。たとえばフィールド212でゾーン15が選択される。データの解像度がフィールド214で選択され、ウィンドウ・フィルタ・サイズがフィールド216で選択される。これは当業者はすぐ認識するであろう。生成された元のDEMは、図5に示されるようにコンピュータ画面300に与えられる。
元のDEMを向上させるためのステップ(2)〜(19)は、図3に示される初期コンピュータ画面200からフィールド220で「バッチ・プロセス実行」を選択することによって開始される。バッチ・プロセスが実行されると、フィールド222、224および226によって、ユーザーは、バッチ・プロセスに関連するある種のパラメータを設定することが許容される。これらのパラメータについては下記で述べる。
ブロック104でのステップ(2)は、元のDEMの再サンプリングである。再サンプリングに関連する設定は、図6に示されるようにコンピュータ画面230で与えられる。元のDEMの再サンプリングはたとえば1メートルの解像度を有しているが、ここで、より低い解像度で再サンプリングされる。その解像度はフィールド232で設定され、たとえば30メートルである。ウィンドウ・フィルタ・サイズも、フィールド234において選択される。結果は図7に示されるようなコンピュータ画面302で与えられる。この結果は、土地52の上のオブジェクト54、56をならしたものである。
ブロック106では、再サンプリングされたDEMに対してヌル充填(null filling)が実行される。該ヌル充填は、図3に示される初期コンピュータ画面200からのフィールド226で与えられるヌル操作と関連している。ヌル操作は、図8のコンピュータ画面240で与えられるようにヌル拡張(null expansion)またはヌル充填に分けられうる。フィールド242がヌル拡張に対応し、フィールド244がヌル充填に対応する。ヌル充填が実行されるので、図9に示されるようなコンピュータ画面250が表示される。ヌル充填に関連する設定は、充填方法についてのフィールド252、実行されるべき充填パス数についてのフィールド254および充填到達距離(fill reach)についてのフィールド256を含む。ヌル充填後に結果として得られる再サンプリングされたDEMは、図10に示されるようにコンピュータ画面304で与えられる。
ブロック108では、DEM引き算(subtraction)が実行される。コンピュータ画面260は、図11に示されるDEM引き算に関連している。DEM引き算で使われる閾値はフィールド262で選択される。ステップ(3)でのヌル充填後の再サンプリングされたDEMが、ステップ(1)での元のDEMから引き算されて、オブジェクトのみDEMが生成される。オブジェクトのみDEMは図12に示されるようなコンピュータ画面306で与えられる。
ブロック110で、ステップ(5)は別のDEM除去(subtraction)である。ステップ(4)からのオブジェクトのみDEMはステップ(1)での元のDEMから除去されて、該オブジェクトのないDEMが生成される。このDEMは図13に示されるようにコンピュータ画面308で与えられる。
ブロック112では、ステップ(6)に対応して、前記のオブジェクトのないDEMに対してヌル拡張が実行される。図14に示されるように、コンピュータ画面270はヌル拡張に関連付けられている。ヌルはフィールド272で選択された値に対応して拡張される。ヌル拡張は、すべてのオブジェクトが除去され、コンピュータ画面310で与えられ、図15に示されるように結果が1メートルの解像度でオブジェクトがないDEMとなることを確実にする。
実行バッチ・プロセスは、ステップ(7)のためにブロック104にループで戻る。ブロック104は、図15で与えられるようなオブジェクトのないDEMに対して再サンプリングを実行する。この再サンプリングはより低い解像度、すなわち1メートルないし30メートルの解像度で実行される。その結果は、図16に示されるようなコンピュータ画面312で与えられる。
ブロック106では、ステップ(8)に対応して二度目のヌル充填が実行される。この二度目のヌル充填は、ブロック112で与えられるような再サンプリングされたオブジェクトのないDEMに対して実行される。このプロセスは、図17のコンピュータ画面314で与えられるような土地のみDEMを生成する。ブロック108では、ステップ(9)に対応して二度目のDEM引き算が実行される。ステップ(8)からの土地のみDEMがここでステップ(1)での元のDEMから引き算されて、向上されたオブジェクトのみDEMが生成される。向上されたオブジェクトのみDEMは図18に示されるようなコンピュータ画面316で与えられる。ブロック110では、ステップ(10)は別のDEM除去ステップである。ステップ(9)からの向上されたオブジェクトのみDEMがステップ(1)での元のDEMから除去されて、向上された、オブジェクトのないDEMが生成される。この向上された、オブジェクトのないDEMは図19に示されるようにコンピュータ画面318で与えられる。
実行バッチ・プロセスは、ステップ(11)のためにブロック104にループで戻る。ブロック104は、ブロック110で与えられるような向上された、オブジェクトのないDEMに対して、今一度再サンプリングを実行する。前回同様、この再サンプリングはより低い解像度、すなわち1メートルないし30メートルの解像度で実行される。その結果は、図20に示されるようなコンピュータ画面320で与えられる。
ブロック106では、ステップ(12)に対応して今一度、ヌル充填が実行される。この三度目のヌル充填は、ブロック104によって与えられるような再サンプリングされたオブジェクトのないDEMに対して実行され、図21のコンピュータ画面322で与えられるような向上された土地のみDEMを生成する。このDEMは、最終的な土地のみDEMとも称される。
ブロック108では、ステップ(13)に対応して三度目のDEM引き算が実行される。ステップ(12)からの向上された土地のみDEMがステップ(1)での元のDEMから引き算されて、一層向上されたオブジェクトのみDEMが生成される。一層向上されたオブジェクトのみDEMは図22に示されるようなコンピュータ画面324で与えられる。
ブロック114では、ステップ(14)に対応して、前記の一層向上されたオブジェクトのみDEMに対してヌル拡張が実行される。ブロック116では、ステップ(15)に対応してヌル充填が実行される。ステップ(14)および(15)は、オブジェクトのまわりからノイズを除去して、一層向上されたオブジェクトのみDEMを生成するために実行される。プロセスはブロック114にループで戻り、ステップ(16)および(17)が実行される。すなわち、今一度ヌル拡張およびヌル充填が実行されて、図23のコンピュータ画面326で与えられるような第二の一層向上されたオブジェクトのみDEMが生成される。
ブロック118では、DEM引き算が実行される。ステップ(17)からの第二の一層向上されたオブジェクトのみDEMが、ステップ(15)からの一層向上されたオブジェクトのみDEMから引き算されて、ノイズのみDEMが生成される。ノイズのみDEMは図24に示されるようなコンピュータ画面328で与えられる。
ブロック120では、別のDEM除去が実行される(ステップ19)。ステップ(18)からのノイズのみDEMはステップ(15)からの第二の一層向上されたオブジェクトのみDEMから除去されて、最終的なオブジェクトのみDEMが生成される。この最終的なオブジェクトのみDEMは図25に示されるようにコンピュータ画面330で与えられる。
ステップ(2)〜(19)について上で論じたようにして、最終的な土地のみDEMおよび最終的なオブジェクトのみDEMが生成された。従来技術のDEMに比べ、これらのDEMは、ステップ(2)〜(19)において実行されたループ反復の結果として向上されている。これから論じる本発明に基づくさらにもう一つの側面は、最終的なオブジェクトのみDEMの建物DEMおよび植生DEMへの分離である。換言すれば、最終的なオブジェクトのみDEMは二つの別個のDEMに分離され、各DEMは別個に処理されうるようになる。
最終的なオブジェクトのみDEMの建物DEMおよび植生DEMへの分離についてこれから図26に示される流れ図および図27〜28を参照しつつ論じる。最初に図3に示されていたコンピュータ画面200を参照するに、ユーザーはフィールド227で「建物と樹木を分離」を選択する。すると、図27に示されるような「建物と樹木を分離」のコンピュータ画面410が表示される。
引き続き図27を参照するに、ユーザーはいくつかの閾値パラメータを選択するオプションを有している。分離は、各オブジェクトについて面積当たりの周および各オブジェクトの高さを計算することに基づいて実行されるので、対応する比較閾値がコンピュータ画面410を介して設定される。たとえば、面積当たりの周の閾値がフィールド412で設定される。考慮すべき各オブジェクトの最小サイズがフィールド414で設定される。このフィールドは最小ポスト(minimum posts)と表示されている。フィールド416でコード残差(chord residue)が選択されるが、これは考慮されているオブジェクトの幅または一側の長さである。最大樹木高さがフィールド418で選択される。第二の面積当たりの周テストに関連した第二の閾値がフィールド420で選択される。フィールド420で選択されたこの第二の閾値は、フィールド412で選択された第一の閾値とは異なっていてもよい。
最終的なオブジェクトのみDEMの建物DEMおよび植生DEMへの分離を開始するため、最終的なオブジェクトのみDEMにおける各オブジェクトについて周対面積パラメータがブロック502で決定される。本発明を説明する目的のため、図28も参照すると、最終的なオブジェクトのみDEMがどのように建物DEMと植生DEMに分離されるかが与えられている。たとえば、最初に、図25で示されたような最終的なオブジェクトのみDEMを表す単純化された表現が図28のフレーム600に与えられる。
フレーム600内に示されるように、オブジェクトは、同じDEM内にまとめられた建物54および樹木56を含む。各オブジェクトについての周対面積パラメータは、ブロック504で選択された閾値と比較されて、DEM中の各オブジェクトが建物54であるか植生56であるかが識別される。閾値との比較に基づいて、フレーム602および604に示されるように、オブジェクトは建物DEMと植生DEMに分離される。
各オブジェクトについてのデータは高さの値を含んでおり、フレーム604中の植生56として識別された各オブジェクトについての高さの値が、ブロック506で高さ閾値と比較される。ブロック508では、フレーム604中の各植生56は、それに関連する高さ値が高さ閾値より大きければ、建物54と識別し直される。フレーム606および608に示されるように、フレーム604で植生56と識別された建物54がフレーム606では建物として識別し直されている。しかしながら、フレーム606に示されるように、フレーム604中の高い樹木56が今度は高さ閾値との比較に基づいて建物と識別されてしまった。
ブロック510では、フレーム606で建物と識別されている各オブジェクトについて、第二の周対面積パラメータが決定される。第二の周対面積パラメータのそれぞれは、ブロック512で第二の閾値と比較される。各建物54は、第二の周対面積パラメータが第二の閾値より大きければ、ブロック514で植生56として識別し直される。
ブロック516では、建物54として識別されたオブジェクトは建物DEMに分離され、植生として識別されたオブジェクトは植生DEMに分離される。別個の建物DEMおよび植生DEMは有利なことに、著しく少ないユーザー介入で、よりリアルな地形モデルが生成されることを許容する。建物DEMはフレーム610によって表され、植生DEMはフレーム612によって表されている。最終的なオブジェクトのみDEMを建物DEMと植生DEMに分離する本方法は、ブロック518で終了する。最終的なオブジェクトのみDEMを二つの別個のDEMに分離する上記の諸ステップは、図2のステップ(20)に対応する。
ここで、残りのステップ(21)〜(26)について論じる。建物54および樹木56はそれぞれ異なってモデル化される。ブロック124はステップ(21)に対応し、任意的であるが、自動プロセスが各オブジェクトを正しく識別しそこなった場合に建物54と樹木56の分離をユーザーが手動で仕上げまたは編集することを許容する。
ブロック126では、ステップ(22)に対応して、植生DEM612中の各植生が3D点としてモデル化される。点をx(緯度)、y(経度)およびz(高度)のリストに変換するために、ユーザーは図3に示されるようなフィールド228で「DEMから点を生成」を選択する。
ブロック128は、ステップ(23)に対応して、建物がベクトルとしてモデル化される。建物をベクトルとしてモデル化することは、背景技術の部分で論じた米国特許第6,654,690号に開示されている。ブロック130で、ステップ(24)に対応して、建物を表すポリゴンにテクスチャーが配置される。換言すれば、ポリゴン上に画像が配置されて、地形モデルにリアルな見かけが与えられる。RealSite(商標)はこのタスクを実行する一つの市販のツールである。RealSite(商標)は本発明の被譲渡者であるハリス・コーポレーションによって開発されたものである。
ブロック132では、コンピュータ・システム60上での表示のために生成された幾何学図形およびテクスチャーのすべてを整形するためにSceneBuilder(商標)が使用される。SceneBuilder(商標)も市販のツールである。ブロック134でInReality(商標)を使って、図29に示されるような表示のための最終的な地形モデルが提供される。InReality(商標)は、ハリス・コーポレーションによって開発されたもう一つの市販のツールであり、ユーザーが仮想シーンを行き来し、さまざまな解析を遂行することを許容する。InReality(商標)は、RealSite(商標)ソフトウェアに対する相棒として設計されている。プロセスはブロック136で終了する。
本発明に基づく、地形データを収集し、収集された地形データから地形モデルを生成することのブロック概略図である。 本発明に基づく地形モデルを生成するための流れ図である。 本発明に基づく元のDEMの生成に対応するコンピュータ画面表示の図である。 本発明に基づく元のDEMの生成に対応するコンピュータ画面表示の図である。 本発明に基づく元のDEMの生成に対応するコンピュータ画面表示の図である。 本発明に基づく元のDEMの再サンプリングに対応するコンピュータ画面表示の図である。 本発明に基づく元のDEMの再サンプリングに対応するコンピュータ画面表示の図である。 本発明に基づく再サンプリングされたDEMに対して実行されるヌル充填プロセスに対応するコンピュータ画面表示の図である。 本発明に基づく再サンプリングされたDEMに対して実行されるヌル充填プロセスに対応するコンピュータ画面表示の図である。 本発明に基づく再サンプリングされたDEMに対して実行されるヌル充填プロセスに対応するコンピュータ画面表示の図である。 本発明に基づく、オブジェクトのみDEMと該オブジェクトのないDEMを生成するためのDEM引き算に対応するコンピュータ画面表示の図である。 本発明に基づく、オブジェクトのみDEMと該オブジェクトのないDEMを生成するためのDEM引き算に対応するコンピュータ画面表示の図である。 本発明に基づく、オブジェクトのみDEMと該オブジェクトのないDEMを生成するためのDEM引き算に対応するコンピュータ画面表示の図である。 図13で与えられるオブジェクトのないDEMに対して実行されるヌル拡張プロセスに対応するコンピュータ画面表示の図である。 図13で与えられるオブジェクトのないDEMに対して実行されるヌル拡張プロセスに対応するコンピュータ画面表示の図である。 本発明に基づく土地のみDEMの生成に対応するコンピュータ画面表示の図である。 本発明に基づく土地のみDEMの生成に対応するコンピュータ画面表示の図である。 本発明に基づく、向上されたオブジェクトのみDEMと向上された、該オブジェクトのないDEMを生成するためのDEM引き算に対応するコンピュータ画面表示の図である。 本発明に基づく、向上されたオブジェクトのみDEMと向上された、該オブジェクトのないDEMを生成するためのDEM引き算に対応するコンピュータ画面表示の図である。 本発明に基づく最終的な土地のみDEMおよび一層向上されたオブジェクトのみDEMの生成に対応するコンピュータ画面表示の図である。 本発明に基づく最終的な土地のみDEMおよび一層向上されたオブジェクトのみDEMの生成に対応するコンピュータ画面表示の図である。 本発明に基づく最終的な土地のみDEMおよび一層向上されたオブジェクトのみDEMの生成に対応するコンピュータ画面表示の図である。 本発明に基づく、ヌル充填/ヌル拡張に基づく第二の一層向上されたオブジェクトのみDEMの生成に対応するコンピュータ画面表示の図である。 本発明に基づく、ノイズのみDEMの生成に対応するコンピュータ画面表示の図である。 本発明に基づく最終的なオブジェクトのみDEMの生成に対応するコンピュータ画面表示の図である。 図25で与えられる最終的なオブジェクトのみDEMを建物DEMおよび植生DEMに分離するための流れ図である。 本発明に基づく、最終的なオブジェクトのみDEMの建物DEMおよび植生DEMへの分離に関連するパラメータを設定するためのコンピュータ画面表示の図である。 本発明に基づく、最終的なオブジェクトのみDEMの建物DEMおよび植生DEMへの分離のためのステップの図的な表現である。 本発明に基づいて生成された地形モデルのコンピュータ画面表示の図である。

Claims (10)

  1. 土地およびその上の複数のオブジェクトについてのデータを含む元のデジタル標高モデル(DEM)を処理するためのコンピュータ実装される方法であって:
    元のDEMからより低い解像度の低解像度DEMを生成する段階と;
    該低解像度DEMと元のDEMとの比較に基づいてオブジェクトのみのDEMを生成する段階と;
    該オブジェクトのみのDEMにおけるノイズを削減する段階とを有する方法。
  2. 請求項1記載のコンピュータ実装される方法であって:
    第一のウィンドウ・サイズに基づいて前記オブジェクトのみのDEMにおけるヌルを同定する段階と;
    前記オブジェクトのみのDEM中に同定されたヌルを拡張する段階と;
    前記オブジェクトのみのDEM中の拡張されたヌルを埋めて第一の洗練されたオブジェクトのみDEMを生成する段階とを有する方法。
  3. 請求項2記載のコンピュータ実装される方法であって:
    前記第一のウィンドウ・サイズに基づいて前記第一の洗練されたオブジェクトのみのDEMにおけるヌルを同定する段階と;
    前記第一の洗練されたオブジェクトのみDEM中に同定されたヌルを拡張する段階と;
    前記第一の洗練されたオブジェクトのみDEM中の拡張されたヌルを埋めて第二の洗練されたオブジェクトのみDEMを生成する段階とをさらに有する方法。
  4. 前記第一の洗練されたオブジェクトのみDEMと前記第二の洗練されたオブジェクトのみDEMとの比較に基づいて、ノイズのみのDEMを生成する段階をさらに有する、請求項3記載のコンピュータ実装される方法。
  5. 前記ノイズのみのDEMを、前記第二の洗練されたオブジェクトのみDEMから引き算して、最終的な洗練されたオブジェクトのみDEMを生成する段階をさらに有する、請求項4記載の方法。
  6. 地形モデル化のためのコンピュータ・システムであって:
    土地およびその上の複数のオブジェクトについてのデータを含む元のデジタル標高モデル(DEM)を処理するためのプロセッサを有しており、前記処理は:
    元のDEMからより低い解像度の低解像度DEMを生成する段階と;
    該低解像度DEMと元のDEMとの比較に基づいてオブジェクトのみのDEMを生成する段階と;
    該オブジェクトのみのDEMにおけるノイズを削減する段階とを有しており、
    当該コンピュータ・システムがさらに、前記プロセッサに結合され、前記処理に基づく地形モデルを表示するディスプレイを有する、コンピュータ・システム。
  7. 請求項6記載のコンピュータ・システムであって、前記プロセッサが、前記オブジェクトのみのDEMにおけるノイズを削減するために:
    第一のウィンドウ・サイズに基づいて前記オブジェクトのみのDEMにおけるヌルを同定する段階と;
    前記オブジェクトのみのDEM中に同定されたヌルを拡張する段階と;
    前記オブジェクトのみのDEM中の拡張されたヌルを埋めて第一の洗練されたオブジェクトのみDEMを生成する段階とを実行する、コンピュータ・システム。
  8. 請求項7記載のコンピュータ・システムであって、前記プロセッサがさらに:
    前記第一のウィンドウ・サイズに基づいて前記第一の洗練されたオブジェクトのみのDEMにおけるヌルを同定する段階と;
    前記第一の洗練されたオブジェクトのみDEM中に同定されたヌルを拡張する段階と;
    前記第一の洗練されたオブジェクトのみDEM中の拡張されたヌルを埋めて第二の洗練されたオブジェクトのみDEMを生成する段階とを実行する、コンピュータ・システム。
  9. 前記プロセッサが、前記第一の洗練されたオブジェクトのみDEMと前記第二の洗練されたオブジェクトのみDEMとの比較に基づいて、ノイズのみのDEMを生成する、請求項8記載のコンピュータ・システム。
  10. 前記プロセッサが、前記ノイズのみのDEMを、前記第二の洗練されたオブジェクトのみDEMから引き算して、最終的な洗練されたオブジェクトのみDEMを生成し;前記ディスプレイが前記最終的な洗練されたオブジェクトのみDEMに基づく地形モデルを表示する、請求項9記載のコンピュータ・システム。
JP2007554275A 2005-02-08 2006-02-06 地形モデル化のためのデジタル標高モデル(dem)を向上させる方法および装置 Withdrawn JP2008530595A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US11/053,495 US7142984B2 (en) 2005-02-08 2005-02-08 Method and apparatus for enhancing a digital elevation model (DEM) for topographical modeling
PCT/US2006/003979 WO2006086253A1 (en) 2005-02-08 2006-02-06 Method and apparatus for enhancing a digital elevation model (dem) for topographical modeling

Publications (1)

Publication Number Publication Date
JP2008530595A true JP2008530595A (ja) 2008-08-07

Family

ID=36779413

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007554275A Withdrawn JP2008530595A (ja) 2005-02-08 2006-02-06 地形モデル化のためのデジタル標高モデル(dem)を向上させる方法および装置

Country Status (6)

Country Link
US (1) US7142984B2 (ja)
EP (1) EP1854049A4 (ja)
JP (1) JP2008530595A (ja)
CA (1) CA2597057C (ja)
TW (1) TWI306229B (ja)
WO (1) WO2006086253A1 (ja)

Families Citing this family (42)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7356774B2 (en) * 2002-08-13 2008-04-08 National Instruments Corporation Grouping components of a measurement system
US7424133B2 (en) 2002-11-08 2008-09-09 Pictometry International Corporation Method and apparatus for capturing, geolocating and measuring oblique images
US7424335B2 (en) * 2005-07-13 2008-09-09 Swift Lawrence W Identification of terrestrial foliage location, type and height for scaled physical models
US8433457B2 (en) * 2006-01-10 2013-04-30 Harris Corporation Environmental condition detecting system using geospatial images and associated methods
US7630797B2 (en) * 2006-01-10 2009-12-08 Harris Corporation Accuracy enhancing system for geospatial collection value of an image sensor aboard an airborne platform and associated methods
US7603208B2 (en) * 2006-01-10 2009-10-13 Harris Corporation Geospatial image change detecting system with environmental enhancement and associated methods
US7750902B2 (en) * 2006-07-20 2010-07-06 Harris Corporation Geospatial modeling system providing non-linear inpainting for voids in geospatial model cultural feature data and related methods
US7376513B2 (en) * 2006-08-09 2008-05-20 Harris Corporation Geospatial modeling system for separating foliage data from building data based upon noise filtering operations and related methods
US7292941B1 (en) * 2006-08-09 2007-11-06 Harris Corporation Geospatial modeling system for separating foliage data from building data based upon loose and strict tolerance noise filtering operations and related methods
US7873238B2 (en) 2006-08-30 2011-01-18 Pictometry International Corporation Mosaic oblique images and methods of making and using same
US8593518B2 (en) 2007-02-01 2013-11-26 Pictometry International Corp. Computer system for continuous oblique panning
US7881913B2 (en) * 2007-02-12 2011-02-01 Harris Corporation Exemplar/PDE-based technique to fill null regions and corresponding accuracy assessment
US8520079B2 (en) 2007-02-15 2013-08-27 Pictometry International Corp. Event multiplexer for managing the capture of images
US8385672B2 (en) 2007-05-01 2013-02-26 Pictometry International Corp. System for detecting image abnormalities
US9262818B2 (en) 2007-05-01 2016-02-16 Pictometry International Corp. System for detecting image abnormalities
US7991226B2 (en) 2007-10-12 2011-08-02 Pictometry International Corporation System and process for color-balancing a series of oblique images
US7983474B2 (en) * 2007-10-17 2011-07-19 Harris Corporation Geospatial modeling system and related method using multiple sources of geographic information
US8531472B2 (en) 2007-12-03 2013-09-10 Pictometry International Corp. Systems and methods for rapid three-dimensional modeling with real façade texture
US8588547B2 (en) 2008-08-05 2013-11-19 Pictometry International Corp. Cut-line steering methods for forming a mosaic image of a geographical area
US20100142814A1 (en) * 2008-12-04 2010-06-10 Harris Corporation Image processing device for tonal balancing of mosaic images and related methods
US8401222B2 (en) 2009-05-22 2013-03-19 Pictometry International Corp. System and process for roof measurement using aerial imagery
US9330494B2 (en) 2009-10-26 2016-05-03 Pictometry International Corp. Method for the automatic material classification and texture simulation for 3D models
US8503761B2 (en) * 2009-11-12 2013-08-06 Harris Corporation Geospatial modeling system for classifying building and vegetation in a DSM and related methods
CN101763659B (zh) * 2010-01-08 2012-03-21 南京师范大学 基于图像集的大规模数字化高程数据模型自动生成方法
US8477190B2 (en) 2010-07-07 2013-07-02 Pictometry International Corp. Real-time moving platform management system
US8823732B2 (en) 2010-12-17 2014-09-02 Pictometry International Corp. Systems and methods for processing images with edge detection and snap-to feature
CA2835290C (en) 2011-06-10 2020-09-08 Pictometry International Corp. System and method for forming a video stream containing gis data in real-time
US10382285B2 (en) * 2011-08-25 2019-08-13 Siemens Industry, Inc. Smart grid communication assessment and co-simulation tool
US9317966B1 (en) * 2012-02-15 2016-04-19 Google Inc. Determine heights/shapes of buildings from images with specific types of metadata
US9183538B2 (en) 2012-03-19 2015-11-10 Pictometry International Corp. Method and system for quick square roof reporting
US9244272B2 (en) 2013-03-12 2016-01-26 Pictometry International Corp. Lidar system producing multiple scan paths and method of making and using same
US9881163B2 (en) 2013-03-12 2018-01-30 Pictometry International Corp. System and method for performing sensitive geo-spatial processing in non-sensitive operator environments
US9275080B2 (en) 2013-03-15 2016-03-01 Pictometry International Corp. System and method for early access to captured images
US9753950B2 (en) 2013-03-15 2017-09-05 Pictometry International Corp. Virtual property reporting for automatic structure detection
WO2015105886A1 (en) 2014-01-10 2015-07-16 Pictometry International Corp. Unmanned aircraft structure evaluation system and method
US9292913B2 (en) 2014-01-31 2016-03-22 Pictometry International Corp. Augmented three dimensional point collection of vertical structures
CA2938973A1 (en) 2014-02-08 2015-08-13 Pictometry International Corp. Method and system for displaying room interiors on a floor plan
US9576373B2 (en) 2015-04-30 2017-02-21 Harris Corporation Geospatial imaging system providing segmentation and classification features and related methods
US10402676B2 (en) 2016-02-15 2019-09-03 Pictometry International Corp. Automated system and methodology for feature extraction
US10671648B2 (en) 2016-02-22 2020-06-02 Eagle View Technologies, Inc. Integrated centralized property database systems and methods
EP3654064B1 (en) * 2018-11-16 2021-01-13 Bayerische Motoren Werke Aktiengesellschaft Apparatus and method for characterizing an object based on measurement samples from one or more location sensors
CN110059385B (zh) * 2019-04-09 2022-10-18 集美大学 一种耦合异速生长的网格动力学情景模拟方法及终端设备

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2759803B1 (fr) * 1997-02-20 1999-03-26 Alsthom Cge Alcatel Procede d'aide a la detection de structures humaines dans un modele numerique de terrain (mnt)
US6249613B1 (en) * 1997-03-31 2001-06-19 Sharp Laboratories Of America, Inc. Mosaic generation and sprite-based coding with automatic foreground and background separation
US6338027B1 (en) * 1999-05-27 2002-01-08 Arborcom Technologies Inc. Canopy modification using computer modelling
US6664529B2 (en) * 2000-07-19 2003-12-16 Utah State University 3D multispectral lidar
AUPR301401A0 (en) * 2001-02-09 2001-03-08 Commonwealth Scientific And Industrial Research Organisation Lidar system and method
US6654690B2 (en) * 2001-04-05 2003-11-25 Harris Corporation Automated method for making a topographical model and related system

Also Published As

Publication number Publication date
CA2597057C (en) 2012-06-05
TWI306229B (en) 2009-02-11
WO2006086253A1 (en) 2006-08-17
EP1854049A4 (en) 2012-02-08
US7142984B2 (en) 2006-11-28
EP1854049A1 (en) 2007-11-14
US20060176213A1 (en) 2006-08-10
TW200636614A (en) 2006-10-16
CA2597057A1 (en) 2006-08-17

Similar Documents

Publication Publication Date Title
US7142984B2 (en) Method and apparatus for enhancing a digital elevation model (DEM) for topographical modeling
US7191066B1 (en) Method and apparatus for distinguishing foliage from buildings for topographical modeling
US11454500B2 (en) Map feature extraction system for computer map visualizations
US6654690B2 (en) Automated method for making a topographical model and related system
Chen et al. A methodology for automated segmentation and reconstruction of urban 3-D buildings from ALS point clouds
Yu et al. Automated derivation of urban building density information using airborne LiDAR data and object-based method
WO2022227910A1 (zh) 虚拟场景生成方法、装置、计算机设备和存储介质
JP2010518458A (ja) 事例/偏微分式ベースの空白領域への埋め込み技術および関連する精度評価
KR102243437B1 (ko) 3d 지형정보 영상 매핑을 위한 타일링 처리를 이용하는 영상 표출 방법 및 시스템
JP2011501301A (ja) 地理的情報の複数のソースを用いた地理空間モデリングシステム及び関連する方法
Truong-Hong et al. Octree-based, automatic building facade generation from LiDAR data
CN113112603A (zh) 三维模型优化的方法和装置
BRPI0714264A2 (pt) sistema de modelaÇço geoespacial, mÉtodo de modelaÇço geoespacial e mÍdia que pode ser lida por computador que possui màdulos executÁveis por computador
JP2004093632A (ja) 地形形状抽出方法、地形形状抽出システム及びプログラム
BRPI0714260A2 (pt) sistema de modelaÇço geoespacial e mÉtodos modelaÇço geoespacial
CN116912817A (zh) 三维场景模型拆分方法、装置、电子设备和存储介质
CN114529689B (zh) 基于对抗神经网络的陶瓷杯缺陷样本扩增方法及系统
US10140758B1 (en) Texturing models of structures
Juckette et al. Using Virtual Reality and Photogrammetry to Enrich 3D Object Identity
KR100782152B1 (ko) 3차원 수치지도를 제작하기 위하여 항공사진 db로부터건물의 3차원 데이터를 획득하는 방법
JP7305919B2 (ja) 建物設計支援装置及び建物設計支援プログラム
JP7232552B1 (ja) 情報処理装置、方法、およびプログラム
Hasegawa et al. Application of stochastic point-based rendering to laser-scanned point clouds of various cultural heritage objects
Otero et al. An application of computer graphics for landscape impact assessment
CN111105488B (zh) 成像仿真方法、装置、电子设备和存储介质

Legal Events

Date Code Title Description
A761 Written withdrawal of application

Free format text: JAPANESE INTERMEDIATE CODE: A761

Effective date: 20101206