JP2623449B2 - 地形モデル作成方法 - Google Patents
地形モデル作成方法Info
- Publication number
- JP2623449B2 JP2623449B2 JP6289295A JP28929594A JP2623449B2 JP 2623449 B2 JP2623449 B2 JP 2623449B2 JP 6289295 A JP6289295 A JP 6289295A JP 28929594 A JP28929594 A JP 28929594A JP 2623449 B2 JP2623449 B2 JP 2623449B2
- Authority
- JP
- Japan
- Prior art keywords
- terrain
- terrain model
- model
- value
- operator
- 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.)
- Expired - Lifetime
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Processing Or Creating Images (AREA)
- Instructional Devices (AREA)
- Image Generation (AREA)
Description
を再現するための地形モデル作成装置および方法、特
に、微細な襞等を含む自然な形状特徴を持つ数値化地形
モデルを等高線データから再現することができる装置お
よび方法に関する。
モデルは、等高線に囲まれた領域の標高値は一定であ
り、領域を取り囲む等高線の平均値であるとしたモデル
である。これはスタッキングモデルと呼ばれている。ス
タッキングモデルの持つ情報は等高線と全く等価であ
る。スタッキングモデルでは、等高線の密度が低い領域
において、各等高線ごとの段差や等高線間の平坦性が顕
著なものになってしまう。そこで、等高線に囲まれた領
域の標高値を変化させ、等高線と等高線の間を連続的な
面で接続することが問題になる。これが等高線からの地
形再現問題である。標高値を変化させるために必要な情
報は、等高線から直接与えられるわけではない。従っ
て、この問題は、等高線デ−タに適切な情報を付加して
妥当な標高面デ−タを得るという、情報付加の問題であ
り、情報の付加内容および情報の付加形式の2つの問題
に分解できる。
される。 (1)幾つかの方向に関して、スプライン曲線の補間等
により地形断面を計算し、 これらの断面に適当に
重みをつけ平均をとる方法。 (2)等高線と等高線の間に三角パッチを張り巡らし、
その三角パッチの上で補間を行なう方法。 (3)スッタキングモデル等の容易に生成できる初期モ
デルに対し、2次元的なロ−パスフィルタ−等を作用さ
せることにより、初期モデルの角ばった部分を滑らかに
する方法。 (4)等高線を点集合とみなし、適当な評価関数を作っ
てこれを最小化することにより、面の当てはめを行なう
方法。
上記(1)の方法の場合、段差や溝および星状のノイズ
が表れ、(2)については、図5に示すように、三角パ
ッチ501が地形の癖として残存してしまう。また上記
(3)、(4)の方法では、図6に示すように、地表面
が細かい襞や褶曲のない異様な丸みをおびたものにな
る、という問題があった。
るのは、情報の付加内容として、地形の地形学的な特質
ではなく、単に幾何学的で技巧的な条件だけを与えるか
らである。また情報の付加形式の点から見ても、従来技
術には困難がある。上記(1)〜(3)の方法は、その
原理上情報の付加形式に従いほぼ情報の付加内容が決定
されるという、一体化処理であり、必要に応じて付加内
容の一部分を追加したり削除することは極めて困難であ
る。また(4)の方法では、評価関数の変更により情報
の付加内容を選択できるが、いずれにせよこの付加内容
は大域的な最適化問題として定式化しておかねばなら
ず、局所的な要請を取り入れることは極めて困難であ
り、選択の余地も狭い。
ので、地形モデル作成装置および方法は、従来の、単に
幾何学的で技巧的な方法が持つ、幾何学的ノイズや不自
然な形状特徴等が結果に表れるという問題点を解消し、
地形学的に妥当で、微細な襞等を含む自然な形状特徴を
持つ数値化地形モデルを、等高線から再現することを目
的とする。
め、本発明においては、等高線を境界値として設定し、
2次元面上で定義された関数空間上の作用素の境界値問
題として地形再現問題を定式化する。情報の付加形式を
境界値問題としたので、作用素の変更により情報の付加
内容を選択でき、局所的な要請を取り入れることも容易
であり、広い範囲の物理学的、地形学的特徴を表現する
ことができる。特に、雨による侵食が主に地形形状を支
配している地域では、情報の付加内容として、地形標高
値は流水線方向に沿って滑らかであるという地形学的特
質を与えることは最も自然であり、このことが地形再現
問題にとって基本的な課題となる。そこで、本発明にお
いては流水線またはその近傍に沿って地形モデルの標高
値を平滑化する作用素を導入することにより、少なくと
もこの地形学的特質を地形モデルに与えることを可能に
した。また、他の作用を付加した作用素を構成すること
により、種々の特質を地形モデルに与えることができ
る。本明細書において、流水線またはその近傍に沿って
標高値を平滑化する作用素とは、この平滑化作用を、少
なくとも作用の一部として含む作用素を意味する。
よび方法は、従来技術が持つ、幾何学的ノイズや不自然
な形状特徴等が結果に表れるという問題を解消し、さら
に、地形学的に妥当で、微細な襞等を含む自然な形状特
徴を持つ数値化地形モデルを等高線から再現することを
可能にした。
る。まず記号法を説明し、これをもとに本発明の実施例
における情報の付加形式を表す。
Aに代入する操作を表す。地形を再現する地域を Dとお
く。簡単化のためこれは正方形であるものとしておく。
地域Dを縦横とも正方形のブロックに N分割して量子化
し、それぞれのブロックを2個の整数の組すなわち点
(x,y) ( x,y = 0,1,2、・・・,N-1)で表す。しかしこの量
子化は説明の簡単化のために導入したものに過ぎず、不
均等な大きさを持つ三角形によるもの等、他の量子化を
使うこともできる。等高線は D内の1次元部分集合であ
り、各等高線を C1, C2,・・・ とし、対応する標高値をそ
れぞれ v1, v2,・・・ とする。
点 (x,y) 上に値(標高値)をとる関数、すなわち2次
元配列 h(x,y) 等により表される。2次元面上で定義さ
れた関数の空間上の作用素とは、2次元配列に2次元配
列を対応させる変換であり、2次元配列 hを h'に変換
する作用素 Tの作用を h'= T(h) と表記する。また、2
個の作用素 S、T の合成変換 h'= S(T(h))、を S・T と
表記する。作用素に関する境界値問題とは、例えば h =
T(h) というような作用素を含む h に関する方程式
を、 D の1次元部分集合 C の任意の点 (x,y) で h(x,
y) = v が成立しなくてはならない、といった制約条件
の下で解き、 h を求める問題である。
の付加形式の1つは、次の境界値問題として設定でき
る。
滑化する作用素を T とする。地形モデルが十分に平滑
化されれば、もうそれ以上平滑化が進行しない平衡状態
に達すると考えられる。このとき地形モデル h は、 h
= T(h) という方程式を、各等高線 Ci ( i = 1,2、・・・)
の任意の点 (x,y) において h(x,y) = vi ( i = 1,2、
・・・) という条件の下で満足する、境界値問題の解とし
て表される。方程式をこの形においたのは求解が比較的
簡単であるからであり、例えば別の考察により、 T'(h)
= 0 とする定式化から出発してもよい。
えられ、作用素 T が初期値の近傍で安定であれば、等
高線の情報を代入し境界値を設定する作用素 B を使
い、反復法により解くことができる。
変換規則により定義される。すなわち、点 (x,y)が等高
線 Ci に含まれる場合は h'(x,y) = vi であり、含ま
れない場合は h'(x,y) = h(x,y) である。
き ho に S を反復作用させ、 h =S・S・ ・・・ ・S(ho)、を
計算すると、この h が境界値問題の近似解になること
が知られている。また解が収束しない場合にも、漸近的
な近似解が得られる。方程式そのものが、現象論的なモ
デルに過ぎないから、解の厳密性をあまり要求しても意
味がない。
を用いてもよい。
形状が非相関的な場合には、得られた近似解 h に対し
て T だけを何回か反復作用させる修正フィルタリング
を行なうことにより、自然な結果が得られるときがあ
る。これは、等高線の情報よりも、作用素のもっている
情報の付加内容を優先する、地形形状決定の方法である
(請求項6に記載の発明)。
図1のブロック図を参照して説明する。実施例の装置1
は、記憶装置2、演算装置3、入力装置4、出力装置
5、制御装置6を有しており、これらはバスライン7,
8によって相互に接続されている。
データを格納する等高線データ格納部、初期地形データ
を作成するプログラムが格納される初期地形データ作成
プログラム部、流水線またはその近傍に沿って入力地形
モデルの標高値を平滑化する作用素の作用プログラムが
格納される平滑化プログラム部、境界条件を設定するプ
ログラムが格納される境界条件設定プログラム部、本装
置により生成された数値化地形データを格納する数値化
地形データ格納部、OS等の制御プログラムを格納する
制御プログラム部を備えている。
装置4はたとえばデジタイザ、マウス、キーボード、数
値ファイルあるいはライトペン等により構成されてい
る。出力装置5はたとえば数値ファイル、ネットワーク
中継装置等である。制御装置6はプログラムを実行する
ための各装置を制御している。
のフローチャートを用いて説明する。まず、初期地形モ
デル hoを作成し(101)、ho を初期値として設定す
る(102)。次に、流水線またはその近傍に沿って標
高値を平滑化する作用素 T を地形モデル h に作用させ
(103)、さらに等高線の情報を代入し境界値を設定
する作用素 B を地形モデル h に作用させる(10
4)。ステップ103,104の演算回数 nをチェック
して(105)、演算回数 n が規定回数 nmax 未満な
らステップ103へ戻り、規定回数 nmaxに達したらス
テップ106へ進んで数値化地形 モデルを hを出力す
る。
ら初期地形モデルを作成する第1のステップはステップ
101とステップ102より成り、流水線またはその近
傍に沿って標高値を平滑化する第2のステップはステッ
プ103からステップ106までより成る。
や、作用素の選択により、本発明の各種の実現形態が生
じる。しかし第1のステップは共通であるから、ここで
説明し、作用素の具体的な定義方法や装置の拡張構成に
ついては、続く各項の実施例の中で説明する。
タッキングモデルを作成し初期地形とすればよい。これ
は、画面領域の塗りつぶし技術等により容易に作成でき
る。等高線に囲まれた領域の1点を取りだし、この点を
種として、この領域を等高線から決まる標高値で塗りつ
ぶす等の手法布つかえばよい。スタッキングモデル以外
にも、三角パッチを利用する補間方法等、著しいノイズ
が発生しない方法であればなんでも初期値作成手段に利
用できる。
点近傍において、流水線の決定、標高値の平滑化等全て
の演算を行なう。作用が局所的であるので流水線の積分
器部分が必要なく、インプリメントが簡単という特徴を
持つ。
る、点(x,y)における作用を次のように定義する
(以下適宜図3を用いて説明する)。 (1)点(x,y)と8点近傍の標高値を基に、この点
における勾配ベクトルgradh(x,y)を決定す
る。このときには、8点近傍全ての点の標高値を計算に
入れることが望ましい。 (2)この勾配ベクトルが0ベクトルではない場合(こ
の点を図3においてP1で示す)平滑化作用素Toを作
用させる。この作用h’=To(h)を次のように定め
る。 1)勾配ベクトル(図3に301で示す)を、水平方
向、垂直方向、右上がり45度方向、左上がり45度方
拘のうち、最も近い方向に量子化する。図3の例では右
上がり45度方向に決定される。 2)量子化方向の直線と8点近傍(図3にU1で示す)
との3つの交点を取り出し、これらの集合を局所流水線
(図3に302で示す)とする。 3)局所流水線の3点の標高値に、総和が1である適当
な重みをつけて平均をとることにより、局所流水線に沿
った平滑化を行ない、点P1:(x,y)における新た
な標高値h’(x,y)を与える。この重みとしては、
図3で示せば、例えば、P1の重みを0とし、他の2点
P2,P3の重みを1/2ずつとしたものを用いること
ができる。 (3)勾配ベクトルが0ベクトルである場合(この点を
図3においてP4で示す)流水線(図3にcで示す)を
延長するために、8点近傍(図3にU2で示す)の点の
標高値に、総和が1である適当な重みをつけて平均をと
ることにより、2次元的な平滑化を行ない、新たな標高
値を与える。この重みとしては、例えば、図3に点P4
の周りに斜線で示す4点(x−1,y),(x+1,
y),(x,y−1),(x,y+1)の重みを1/4
ずつとし、他の点の重みを0としたものを用いることが
できる。
は、前記手段に限られるわけではなく、勾配ベクトルに
直接依存するフィルタ−や、非線型フィルタ−等であっ
てもよい。また、上記(3)の平滑化手段も、非線型フ
ィルタ−等であってもよい。たとえば、hx=h(x+1,y)-h
(x-1,y), hy=h(x,y+1)-h(x,y-1)としたとき、 h'(x,y)=(h(x-1,y)+h(x+1,y))・|hx|/4(|hx|+|hy|) +(h(x,y-1)+h(x,y+1))・|hy|/4(|hx|+|hy|) で定まる非線形作用素を利用することができる。初期値
として滑らかな地形モデルを選べば、上記(3)の場合
はほとんど生じることはなく、(1)、(2)だけから
構成される作用素を使うこともできる。
組み込んで使用すれば、数値化地形モデルが得られる。
この実施例により、従来技術の持つ、幾何学的ノイズや
不自然な形状特徴等が結果に表れるという問題は解消で
きる。
は、作用素の作用の範囲を広げ、積分器を用いて必要な
長さの流水線を求め、ある程度広い範囲で平滑化を行な
う。また、尾根、谷領域に対してはこの平滑化の度合を
低減することによりよりシャープな尾根、谷の地形再現
を行う。
に対する作用を次のように定義する。 (1)地形モデル h の勾配ベクトル場 grad h を決定
する。
を通る整数格子点上の流水線を求める離散的な積分演
算を行う。流水線 c(t)は次のとおり定義される。 c(t) =(c1(t), c2(t)) (c(0) =(x,y), c(t) ≠ c(t+1)) 流水線 c(t)は、上記条件に加えて、|c(t)-c(t-1)|≠
0 かつ |c(t+1)-c(t)|≠ 0 のとき、内積 (c(t)-c
(t-1), c(t+1)-c(t)) が正であるという条件を満足する
ことが望ましい。 この積分は、流水線の始点から
終点までの全域に渡ってではなく、平滑化 手段の
作用範囲程度の局所的な領域上だけで行なってもよい。
上で存在する場合。平滑化作用素 T1 を作用させる。こ
の作用、h' = T1(h) を、 h'(x,y) = Σ w(c(t)-(x,y)) h(c(t)) / Σ w(c(t)-(x,y)) で定める。ここで和は両者とも t についてとり、この
範囲は、流水線と、 重み w の非零定義領域
(台)とを交差させる t の値全体とする。ここで重み
w は例えば、w(x,y) = exp( -(x2+y2)/(2σ2)) ( -3σ
< x,y < 3σ,σ:定数)w(x,y) = 0
( x,y が上の範囲以外のとき)とすればよい。これ
は重み付けをガウス分布的にした例である。ここに示し
た作用素は最も簡単なものであり、非線型フィルタ−等
を利用してもよい。
合。 1)点 (x,y)が尾根、谷領域の点に該当するか否かの判
定を行なう。この判定は、例えば、 a+ = h(x+1,y)-h(x,y), a- = h(x-1,y)-h(x,y), b+ = h(x,y+1)-h(x,y), b- = h(x,y-1)-h(x,y), とおき、a+ ・ a- > 0 または b+ ・ b- > 0 のとき該
当、これ以外のとき非該当、とすることにより行なうこ
とができる。この方法は一例に過ぎず、より大域的な方
法等、これ以外の方法で判定を行なってもよい。2)点
(x,y)が尾根、谷領域の点である場合。 平滑化作用素 T2 を作用させる。この作用、h' = T2(h)
を例えば、 h'(x,y) = ( h(x,y) + wd(x,y) ・ hL(x,y)) /( 1 + wd(x,y)) で定める。ここで hL は、 hL(x,y) = ( h(x-1,y) + h(x+1,y) + h(x,y-1) + h(x,y+1))/ 4 とした。wd は点 (x,y) の尾根度、谷度により変化する
重みであり、例えば、h のラプラシアン(物理的には水
の流入出量を表わす)等の関数として与えてやればよ
い。特に、 T2 の定義方法は色々な手法を考えることが
できる。例えば、 wd を、勾配ベクトルが零である点か
らの距離により増加する関数にとることや、hL の代わ
りに、非線型フィルタ−等の作用結果を利用することも
できる。 3)2)以外の場合。 前記 hL の定める平滑化作用素 h'(x,y) = hL(x,y) を
作用させる。
域に関する判定を行なわずに、常に上記3)の平滑化を
行なう、やや簡単な構成を考えることもできる。この場
合は、実施例1の作用素 To を単に T1 に置き換えた構
成と見なすことができる。
組み込んで使用すれば、数値化地形モデルが得られる。
本実施例によれば、等高線密度の低い地域まで処理が及
びにくい、尾根、谷領域の地形が平坦になりやすい、お
よび地形斜面がややシャ−プさに欠ける、といった問題
を解消し、地形学的に妥当で、微細な襞等を含む自然な
形状特徴を持つ数値化地形モデルを、等高線から再現す
ることができる。
めの、拡張的構成を示す。本境界値問題の解法に反復法
を用いる場合、平滑化の平衡状態に達するまでに必要な
作用素の反復回数は、一般にほぼ等高線同士の間隔に比
例する。例えば、量子化単位の大きさを、基本単位の4
倍にして、地域Dの各辺の分割数をN/4とした量子化
では、等高線同士のブロック間隔は1/4になり、した
がって反復回数も1/4になる。また計算すべき格子点
の数も1/4×1/4すなわち1/16になる。よって
基本の量子化に比べ、1/64の計算量で結果が得られ
る。量子化単位を基本単位の2倍に変更し、各辺の分割
数がN/2である格子面上に、この結果を補間拡大して
得られる地形モデルは、スタッキングモデルに比べては
るかに良い平衡状態の近似になっている。したがって、
この量子化の段階での反復回数はごくわずかで済む。以
下、結果を補間拡大する同様な手順を再帰的に実行する
ことにより、基本の量子化での地形モデルが得られる。
る。 (1)入力等高線データのスケールを縦横それぞれ1/
4に縮小し、実施例2の装置により1/4スケールの地
形モデルを作成する。 (2)前記作成結果を2倍に補間拡大する。 (3)前記拡大結果を初期値とし、また入力等高線デー
タのスケールを縦横それぞれ1/2に縮小し、これを境
界値として、実施例2の第2のステップにより、1/2
スケールの地形モデルを作成する。 (4)前記作成結果を2倍に補間拡大する。 (5)前記拡大結果を初期値とし、入力等高線データを
境界値として、実施例2の第2のステップにより、もと
の1/1スケールの地形モデルを作成する。
としたが、別の比をとってもよく、各段階ごとに変化さ
せてもよい。開始スケ−ルも 1/4 である必要はない。
補間拡大の際に、平面の当てはめ等の1次補間を利用し
て、補間点の標高値をその最隣接点の標高値から求める
ことは望ましくない。このような補間では、尾根や谷は
無視されて平坦な拡大が行なわれるからである。一旦失
われた凹凸を第2のステップにより復元することは困難
であり、尾根や谷は潰れたままになってしまう。
ュの補間法等による3次曲面の当てはめ等を利用すれ
ば、地形の凹凸まで補間できる。図4はそのことを説明
する地形断面図であり、図の横軸はx軸を、縦軸は標高
値を表わしている。点P1〜P4は標高値が既知の点であ
り、Pは標高値を補間したい点である。実線cは1次補
間した平面であり、破線c’は高次補間した曲面であ
る。hは1次補間した場合の標高値であり、h’は高次
補間した場合の標高値であって、両者の間に図に示すよ
うな差が生じる。このように地形モデルの補間拡大に
は、上記のような高次補間が望ましい。
式等の高速解法として知られていたものである。従来、
この方法はラプラス方程式等の境界値問題に適用され、
通常、数倍程度の高速化効果が得られていた。この方法
は、その構成から明らかなように、解の周波数成分の大
きさが、周波数の増大につれてゆっくり減少するという
構造をもつ場合に、解の概形を求める際、特に有効であ
る。しかし、従来の問題では、解がフラクタル(自己相
似)構造とよばれるこのような構造を持たなかったうえ
に、結果の厳密性が必要とされることが多かった。
因は有利に働く。というのも、近年の地形に対するフラ
クタル科学の成果により、地形モデルがもともとフラク
タル構造を持つことが解明されており、また、結果の厳
密性も本質的には必要ないからである。実在の幾つかの
地形の等高線を入力とするテストによれば、数十倍から
百倍程度の高速化が可能である。従来の通常の問題と高
速化の度合いを比較すると、十倍以上の効果が得られて
いる。従って、公知の方法を新しい問題に、単に適用し
たという以上の顕著な効果が認められる。
化作用素の、反復作用の回数を著しく低減できるので、
尾根、谷領域での平滑化抑制効果も得られ、尾根、谷筋
の通った地形モデルが得られる。
辺図)は、実施例3の装置により作成された地形モデル
を、実施例2の作用素 T を用いた請求項6に記載の修
正フィルタによって修正して作成した地形モデルを利用
して作成したものである。この地形モデルに対して仮想
的に左上方向から光を当て、この光線に対する地形表面
の各点における傾きより陰影を計算し、さらに地形モデ
ルの標高値により変化する色調を計算し、この陰影と色
調を組み合わせた画像データを作成した。この画像デー
タに、実施例3の装置へ入力した等高線データおよびそ
の他のデータを重ねて印刷したものが添付参考図であ
る。同参考図においては、シャープな尾根線、谷線ある
いは細かい山の襞などが表現され、全体として自然な地
形再現に成功している。
ル作成装置および方法によれば、従来の、単に幾何学的
で技巧的な方法が持つ、幾何学的ノイズや不自然な形状
特徴等が結果に表れるという問題を解消し、地形学的に
妥当で、微細な襞等を含む自然な形状特徴を持つ数値化
地形モデルを、等高線から再現することができる。
ブロック図である。
る。
る。
形モデルの例を示す図である。
形モデルの他の例を示す図である。
Claims (11)
- 【請求項1】 計算機によって、地形標高値は流水線方
向に沿って滑らかであるという地形学的特質をシミュレ
ートするために、等高線データから、スタッキングモデ
ル等の2次元面上で定義された関数として地形モデルを
計算し、このモデルを初期地形として設定する第1のス
テップと、等高線を境界値とし、流水線またはその近傍
に沿って入力地形モデルの標高値を平滑化する作用素に
関する境界値問題を、前記初期地形を初期値とする反復
法により近似的に解いて、地形モデルを求める第2のス
テップとを備えたことを特徴とする地形モデル作成方
法。 - 【請求項2】 前記作用素が、入力地形モデルの勾配ベ
クトル場を計算し、これを積分して流水線を少なくとも
局所的に求め、前記入力地形モデルにおける各流水線上
の標高値の凹凸をそれぞれの流水線またはその近傍に沿
って平滑化する請求項1に記載の地形モデル作成方法。 - 【請求項3】 前記作用素が、入力地形モデルの尾根領
域および谷領域を検出し、前記勾配ベクトルが零になる
点の近傍で2次元的に平滑化し、前記検出された領域お
よびその近傍で前記2次元的平滑化の度合を低減させる
請求項2に記載の地形モデル作成方法。 - 【請求項4】 地形モデルの定義領域の量子化単位を必
要な精度より大きく設定し、請求項1に記載の第1のス
テップおよび第2のステップによって低精度の地形モデ
ルを作成する第3のステップと、前記量子化単位を小さ
く設定し直し、これに応じて前記地形モデルを補間拡大
して初期値として設定し、この初期値に対して再び請求
項1に記載の第2のステップを行う第4のステップと、
必要な精度に達するまで前記第4のステップを再帰的に
反復して行う第5のステップとを備えたことを特徴とす
る地形モデル作成方法。 - 【請求項5】 前記補間拡大は、補間点の最隣接近傍よ
り広い範囲の標高値を利用する高次補間である請求項4
に記載の地形モデル作成方法。 - 【請求項6】 計算機によって、地形標高値は流水線方
向に沿って滑らかであるという地形学的特質をシミュレ
ートするために、等高線データから計算された2次元面
上で定義された関数である地形モデル作成結果を初期値
とし、境界条件をつけることなく流水線またはその近傍
に沿って入力地形モデルの標高値を平滑化する作用素
を、前記初期値に対し一回または再帰的に作用させるこ
とを特徴とする地形モデル修正方法。 - 【請求項7】 計算機によって、地形標高値は流水線方
向に沿って滑らかであるという地形学的特質をシミュレ
ートするために、等高線データから、2次元面上で定義
された関数としての地形モデルを作成する方法であっ
て、前記等高線データを境界値とし、流水線またはその
近傍に沿って地形モデルの標高値を平滑化する作用素に
関する境界値問題を近似的に解くことにより地形モデル
を求めることを特徴とする地形モデル作成方法。 - 【請求項8】 請求項1または7に記載の方法を使用し
て地形モデルを作成し、この地形モデルに対して仮想的
に光を当て、この光線に対する地形表面の各点における
傾きより陰影を計算して画像データを作成することを特
徴とする画像データ作成方法。 - 【請求項9】 請求項1または7に記載の方法を使用し
て地形モデルを作成し、この地形モデルに対して仮想的
に光を当て、この光線に対する地形表面の各点における
傾きより陰影を計算し、さらに地形モデルの標高値によ
り変化する色調を計算し、前記陰影と前記色調とを組み
合わせた画像データを作成することを特徴とする画像デ
ータ作成方法。 - 【請求項10】 請求項1または7に記載の方法を使用
して地形モデルを作成し、この地形モデルに対して仮想
的に光を当て、この光線に対する地形表面の各点におけ
る傾きより陰影を計算して画像データを作成し、この画
像データを印刷することを特徴とする地形モデル図作成
方法。 - 【請求項11】 請求項1または7に記載の方法を使用
して地形モデルを作成し、この地形モデルに対して仮想
的に光を当て、この光線に対する地形表面の各点におけ
る傾きより陰影を計算し、さらに地形モデルの標高値に
より変化する色調を計算し、前記陰影と前記色調とを組
み合わせた画像データを作成し、この画像データを印刷
することを特徴とする地形モデル図作成方法。
Priority Applications (7)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP6289295A JP2623449B2 (ja) | 1994-10-29 | 1994-10-29 | 地形モデル作成方法 |
DE69521436T DE69521436T2 (de) | 1994-10-29 | 1995-10-26 | Vorrichtung und verfahren für die wiedergabe von terrainmodellen |
PCT/JP1995/002202 WO1996013812A1 (en) | 1994-10-29 | 1995-10-26 | Apparatus and method for terrain model reproduction |
AU37539/95A AU702316B2 (en) | 1994-10-29 | 1995-10-26 | Apparatus and method for terrain model reproduction |
US08/817,483 US6272448B1 (en) | 1994-10-29 | 1995-10-26 | Apparatus and method for terrain model reproduction |
EP95935572A EP0788639B1 (en) | 1994-10-29 | 1995-10-26 | Apparatus and method for terrain model reproduction |
CA002202655A CA2202655C (en) | 1994-10-29 | 1995-10-26 | Apparatus and method for terrain model reproduction |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP6289295A JP2623449B2 (ja) | 1994-10-29 | 1994-10-29 | 地形モデル作成方法 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP8312679A Division JPH09106451A (ja) | 1996-11-08 | 1996-11-08 | 地形モデル作成装置および地形モデル作成プログラムを記憶した媒体 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH08129636A JPH08129636A (ja) | 1996-05-21 |
JP2623449B2 true JP2623449B2 (ja) | 1997-06-25 |
Family
ID=17741339
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP6289295A Expired - Lifetime JP2623449B2 (ja) | 1994-10-29 | 1994-10-29 | 地形モデル作成方法 |
Country Status (7)
Country | Link |
---|---|
US (1) | US6272448B1 (ja) |
EP (1) | EP0788639B1 (ja) |
JP (1) | JP2623449B2 (ja) |
AU (1) | AU702316B2 (ja) |
CA (1) | CA2202655C (ja) |
DE (1) | DE69521436T2 (ja) |
WO (1) | WO1996013812A1 (ja) |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7216003B2 (en) * | 2000-04-06 | 2007-05-08 | Solid Terrain Modeling | Hi-accuracy three-dimensional topographic modeling apparatus |
US8958654B1 (en) * | 2001-04-25 | 2015-02-17 | Lockheed Martin Corporation | Method and apparatus for enhancing three-dimensional imagery data |
KR100491723B1 (ko) * | 2001-12-20 | 2005-05-27 | 한국전자통신연구원 | 지리정보시스템에서의 디지털지형모델 가시화방법 |
FR2837949A1 (fr) * | 2002-03-29 | 2003-10-03 | Emg | Procede pour la visualisation d'une image geographique semi-reelle |
JP3670274B2 (ja) | 2002-11-05 | 2005-07-13 | アジア航測株式会社 | 視覚化処理システム、視覚化処理方法、及び視覚化処理プログラム |
US7006106B2 (en) * | 2004-03-09 | 2006-02-28 | Honeywell International Inc. | Real time contour line generation |
US20050251332A1 (en) * | 2004-03-22 | 2005-11-10 | Daimlerchrysler Ag | Method for correlating altitude and/or grade information with route points of a digital map |
US8334079B2 (en) * | 2004-04-30 | 2012-12-18 | NanoCell Systems, Inc. | Metastable ceramic fuel cell and method of making the same |
US8207966B2 (en) * | 2004-05-26 | 2012-06-26 | Sony Online Entertainment Llc | Terrain editor tool for rule-based procedural terrain generation |
US8115765B2 (en) * | 2004-05-26 | 2012-02-14 | Sony Online Entertainment Llc | Rule-based procedural terrain generation |
KR100916474B1 (ko) * | 2009-04-28 | 2009-09-08 | 최준수 | 등고선이 표현하는 지역을 다각형 세부 영역들로 분해하고 수치표고모델 데이터와 지형 표면 구성 데이터를 산출하는 방법과, 그 방법을 실행하기 위한 프로그램이 기록된 기록 매체 |
CN102117500B (zh) * | 2009-12-30 | 2015-03-25 | 新奥特(北京)视频技术有限公司 | 一种三维模型构建方法及*** |
KR101764341B1 (ko) * | 2012-01-03 | 2017-08-02 | 한국전자통신연구원 | 절차적 방법에 의해 생성된 지형 데이터를 편집하는 방법 |
CN102831639A (zh) * | 2012-07-13 | 2012-12-19 | 南信大影像技术工程(苏州)有限公司 | 一种应用于智能电网三维运行管理***的层次细节模型 |
CN106097450B (zh) * | 2016-06-21 | 2019-01-22 | 武汉瑞得信息工程有限责任公司 | 一种等高线生成方法及装置 |
US10535119B2 (en) * | 2017-08-11 | 2020-01-14 | Intermap Technologies Inc. | Method and apparatus for enhancing 3D model resolution |
US11532093B2 (en) | 2019-10-10 | 2022-12-20 | Intermap Technologies, Inc. | First floor height estimation from optical images |
US11551366B2 (en) | 2021-03-05 | 2023-01-10 | Intermap Technologies, Inc. | System and methods for correcting terrain elevations under forest canopy |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2091526B (en) * | 1981-01-13 | 1985-10-02 | Harris Corp | Digital map generator and display system |
WO1988002156A2 (en) * | 1986-09-11 | 1988-03-24 | Hughes Aircraft Company | Digital simulation system for generating realistic scenes |
US4823287A (en) * | 1987-05-28 | 1989-04-18 | Honeywell Inc. | Digital contour line generator |
US5053778A (en) * | 1989-08-10 | 1991-10-01 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Generation of topographic terrain models utilizing synthetic aperture radar and surface level data |
US5168161A (en) * | 1990-04-18 | 1992-12-01 | Texas Instruments Incorporated | System and method of determining surface characteristics using infrared imaging |
DE69204241T2 (de) * | 1991-03-27 | 1996-02-29 | Exxon Production Research Co | Darstellung von n-dimensionalen seismischen Daten in einem n-1 dimensionalen Format. |
WO1993004437A1 (en) * | 1991-08-15 | 1993-03-04 | Hughes Aircraft Company | Efficient storage of geographic data for visual displays |
US5815596A (en) * | 1994-04-14 | 1998-09-29 | Narendra Ahuja | Multiscale image edge and region detection method and apparatus |
FR2731819B1 (fr) * | 1995-03-17 | 1997-04-11 | Alsthom Cge Alcatel | Methode d'extraction de contours utilisant une analyse multi-fractale |
JP3522408B2 (ja) * | 1995-09-18 | 2004-04-26 | 富士通株式会社 | 数値流体解析結果の誤差見積方法、数値流体解析結果の誤差見積装置、数値流体解析方法、及び数値流体解析装置 |
US5798923A (en) * | 1995-10-18 | 1998-08-25 | Intergraph Corporation | Optimal projection design and analysis |
US5835099A (en) * | 1996-06-26 | 1998-11-10 | Xerox Corporation | Representing a region of a color image using a space-color separable model |
-
1994
- 1994-10-29 JP JP6289295A patent/JP2623449B2/ja not_active Expired - Lifetime
-
1995
- 1995-10-26 CA CA002202655A patent/CA2202655C/en not_active Expired - Fee Related
- 1995-10-26 WO PCT/JP1995/002202 patent/WO1996013812A1/en active IP Right Grant
- 1995-10-26 DE DE69521436T patent/DE69521436T2/de not_active Expired - Fee Related
- 1995-10-26 EP EP95935572A patent/EP0788639B1/en not_active Expired - Lifetime
- 1995-10-26 US US08/817,483 patent/US6272448B1/en not_active Expired - Fee Related
- 1995-10-26 AU AU37539/95A patent/AU702316B2/en not_active Ceased
Also Published As
Publication number | Publication date |
---|---|
EP0788639B1 (en) | 2001-06-20 |
US6272448B1 (en) | 2001-08-07 |
DE69521436D1 (de) | 2001-07-26 |
JPH08129636A (ja) | 1996-05-21 |
DE69521436T2 (de) | 2002-05-02 |
CA2202655A1 (en) | 1996-05-09 |
CA2202655C (en) | 2005-08-16 |
AU702316B2 (en) | 1999-02-18 |
AU3753995A (en) | 1996-05-23 |
EP0788639A1 (en) | 1997-08-13 |
WO1996013812A1 (en) | 1996-05-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2623449B2 (ja) | 地形モデル作成方法 | |
JP4012258B2 (ja) | プログレッシブメッシュの適応細分方法および装置 | |
Desmet | Effects of interpolation errors on the analysis of DEMs | |
US7170516B2 (en) | Approximation of Catmull-Clark subdivision surfaces by Bezier patches | |
US6639592B1 (en) | Curve network modeling | |
JP2006523891A (ja) | 局所洗練を用いてt−スプライン及びt−nurcc表面を定義するためのシステム及び方法 | |
US20100246973A1 (en) | Image processing apparatus and image processing method | |
Pfeifer | A subdivision algorithm for smooth 3D terrain models | |
Shen et al. | A PDE patch-based spectral method for progressive mesh compression and mesh denoising | |
US20010030646A1 (en) | System and method for multi-resolution fairing of non-manifold models | |
KR100414058B1 (ko) | 3차원 다각형 표면의 리메슁 최적화 방법 및 장치 | |
JP2837584B2 (ja) | 地形データの作成方法 | |
US7453457B2 (en) | Computer graphics using coarse level meshes | |
JPH09106451A (ja) | 地形モデル作成装置および地形モデル作成プログラムを記憶した媒体 | |
Nelson et al. | Adaptive tessellation method for creating TINs from GIS data | |
Dyn et al. | Boundary correction for piecewise linear interpolation defined over data-dependent triangulations | |
WO2000075865A1 (fr) | Procede de traitement d'image | |
Pakdel et al. | Incremental catmull-clark subdivision | |
JP3103790B2 (ja) | 三次元オブジェクトのモデリング方法、装置および記録媒体 | |
JPH11339071A (ja) | 点群からのポリゴン自動生成システム | |
JP2008152511A (ja) | 等高線編集装置、等高線編集方法、及び等高線編集プログラム | |
US7308391B1 (en) | Universal hydraulic solver with techniques for improving the representation of flow data | |
Shen | Building interpolating and approximating implicit surfaces using moving least squares | |
Brooks et al. | Integrating procedural textures with replicated image editing | |
JPH07220118A (ja) | 地形表示装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080411 Year of fee payment: 11 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090411 Year of fee payment: 12 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090411 Year of fee payment: 12 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100411 Year of fee payment: 13 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110411 Year of fee payment: 14 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120411 Year of fee payment: 15 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120411 Year of fee payment: 15 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130411 Year of fee payment: 16 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130411 Year of fee payment: 16 |
|
S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130411 Year of fee payment: 16 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130411 Year of fee payment: 16 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150411 Year of fee payment: 18 |
|
EXPY | Cancellation because of completion of term |