JP4177017B2 - オルソフォト生成方法、オルソフォト生成システム及びオルソフォト生成プログラム - Google Patents
オルソフォト生成方法、オルソフォト生成システム及びオルソフォト生成プログラム Download PDFInfo
- Publication number
- JP4177017B2 JP4177017B2 JP2002118033A JP2002118033A JP4177017B2 JP 4177017 B2 JP4177017 B2 JP 4177017B2 JP 2002118033 A JP2002118033 A JP 2002118033A JP 2002118033 A JP2002118033 A JP 2002118033A JP 4177017 B2 JP4177017 B2 JP 4177017B2
- Authority
- JP
- Japan
- Prior art keywords
- triangle
- orthophoto
- image
- aerial
- basic
- 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 - Fee Related
Links
- 238000000034 method Methods 0.000 title claims description 41
- 238000003860 storage Methods 0.000 claims description 23
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 claims description 9
- 230000002194 synthesizing effect Effects 0.000 claims description 2
- 230000008707 rearrangement Effects 0.000 claims 3
- 238000010586 diagram Methods 0.000 description 21
- 238000004364 calculation method Methods 0.000 description 13
- 230000006870 function Effects 0.000 description 13
- 238000006243 chemical reaction Methods 0.000 description 10
- 230000014509 gene expression Effects 0.000 description 10
- 239000000872 buffer Substances 0.000 description 9
- 230000009466 transformation Effects 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 4
- 238000013500 data storage Methods 0.000 description 4
- 238000002360 preparation method Methods 0.000 description 4
- 238000003786 synthesis reaction Methods 0.000 description 4
- 241000839309 Thesea Species 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 2
- 208000006440 Open Bite Diseases 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000010365 information processing Effects 0.000 description 2
- 241000772991 Aira Species 0.000 description 1
- 101150064138 MAP1 gene Proteins 0.000 description 1
- XSTXAVWGXDQKEL-UHFFFAOYSA-N Trichloroethylene Chemical compound ClC=C(Cl)Cl XSTXAVWGXDQKEL-UHFFFAOYSA-N 0.000 description 1
- 230000003139 buffering effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 101150077939 mapA gene Proteins 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000012536 storage buffer Substances 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Images
Landscapes
- Instructional Devices (AREA)
- Processing Or Creating Images (AREA)
- Image Processing (AREA)
Description
【発明の属する技術分野】
本発明は、空中撮影された地表画像から等縮尺で平行投影したオルソフォトを生成するオルソフォト生成方法、オルソフォト生成システム及びオルソフォト生成プログラムに関する。
【0002】
【従来の技術】
航空機等による空中撮影に用いられる画像センサは一般に中心投影であり、画像センサに近い物ほど大きく写り、速いものほど小さく写る。このため、このような画像センサで撮影された地表の画像を、地図のように等縮尺で描かれている画像と単純に重ね合わせることはできない。
【0003】
このような撮影時の縮尺比の影響を、地表面モデルを用いて補正し、図24に示すように、鈷直方向に向かって平行投影によって投影した場合の画像を人工的に再現したものがオルソフォトである。このようなオルソフォトは、カメラの標定要素と地表面の形状を示すモデルがあれば生成する事ができる。
【0004】
空気の密度や温度差、レンズの特性等による画像の補正を行い、いわゆる共線条件を満たす場合には、カメラの標定要素(内部標定要素・外部標定要素)が既知であるとすると、3次元空間から画像平面への写像、つまり地上座標系で座標(X,Y,Z)の点Pが写るデジタル写真画像I上のピクセル座標(xI,yI)を求める関数は次式で与えられる。
【0005】
xI=Fx(X,Y,Z) (1)
yI=Fy(X,Y,Z) (2)
具体的には、図25に示すように、カメラの投影中心の座標を(X0,Y0,Z0)、カメラの地上座標に対する回転角を(κ,φ,ω)、カメラの焦点距離をcとすると、地上座標(X,Y,Z)はカメラの投影中心と傾きに相対的な座標系(カメラ座標系)(U,V,W)に次式で変換される。
【0006】
【数1】
また、この座標が写る写真上の座標(写真座標系、投影中心の像(主点)を中心とする2次元座標系)(x,y)は、次式で表される。
【0007】
x=−c・U/W (4)
y=−c・V/W (5)
更に写真をデジタル画像Iに変換するとき、非線形な歪が入らなければ、写真座標(x,y)は一般に次式に示すアフィン変換で画像座標(xI,yI)に変換することができる。
【0008】
【数2】
以上より、投影中心座標・回転角・焦点距離・アフィン変換係数がわかれば、地上座標(X,Y,Z)を画像座標(xI,yI)に変換することができる。この手続きをFx,Fyとすると、上述の式(1)及び式(2)のようになる。
【0009】
また、地表面の平面座標(X,Y)における標高値Zを与える関数(地表面モデル)が、Z=S(X,Y)で与えられるとすると、平面座標(X,Y)に対応する地表面上の座標は(X,Y,S(X,Y))である。このため、平面座標(X,Y)に対応する写真画像上の座標(xI,yI)は、次式のようになる。
【0010】
xI=Fx(X,Y,S(X,Y)) (7)
yI=Fy(X,Y,S(X,Y)) (8)
なお、図26に示すように、生成したいオルソフォトの原点を(X0,Y0),オルソフォトのカラム方向の単位方向ベクトルをec=(Dxx,Dxy)、オルソフォトのライン方向の単位方向ベクトルをel=(Dyx,Dyy)、カラム軸方向の解像度をdc,ライン紬方向の分解能をdlとすれば、オルソフォト画像の座標(i,j)に対応する平面座標(X,Y)の関係は、次式で与えられる。
【0011】
【数3】
【数4】
式(7),式(8)と式(9)を組み合わせれば、オルソフォト画像Iorthの座標(i,j)での写真画像上の座標を得ることができる。最終的に写真画像Iを用いると、
Iorth(i,j)=I(Fx(X,Y,S(X,Y)),Fy(X,Y,S(X,Y))
と表現することができる。なお、座標(Fx(X,Y,S(X,Y)),Fy(X,Y,S(X,Y)))の座標値が整数ではない場合は、周囲の画素値を適当な方法(例えばニアレストネイバー・バイリニア・キュービックコンボルーション等)で内挿計算する。
【0012】
【発明が解決しようとする課題】
ところで、上述のようなアルゴリズムは、オルソフォト生成のために一般的に用いられているが、このアルゴリズムでは、例えば地面や建物が別の建物の影になって写真画像に写っていない「オクルージョン」が起こらないことを前提としている。このオクルージョンとは、例えば地面や建物が別の建物の影になって写真画像に写っていないような状況のことである。
【0013】
このようなオクルージョンが起こった状況では、図27に示すように、一つの写真画像上の座標(xI,yI)に複数の地表座標(X1,Y1,Z1),(X2,Y2,Z2),…,(XN,YN,ZN)が対応してしまう(同図はN=3の場合を示している)。これらの地表座標のうち、実際に画像に写っているものに対応するのは、カメラに最も近い点である。すなわち、この図27の場合では、写真画像のP点には、3次元空間上のP1,P2,P3の各点が投影されるが、実際にはP2,P3は点P1によって隠されて写真画像には写らない。
【0014】
このような状況を考慮しないで、上述の全ての座標{(X1,Y1,Z1),(X2,Y2,Z2),…,(XN,YN,ZN)}に対して上述の式(11)を適用してしまうと、図28に示すように、オルソフォト上で同じ像が何度も写ったり、建物の壁面がオルソフォトに混入したりしてしまうことになる。
【0015】
このような問題を解決するためには、写真画像上の座標(xI,yI)に対応する地表面モデル上の点(X1,Y1,Z1),(X2,Y2,Z2),…,(XN,YN,ZN)を求め、これらの投影中心に対する前後関係を調べ、最も投影中心に近いもののみをオルソフォト化するという手続が必要である。この手続はコンピュータグラフィックスにおけるZバッファリングの逆手法ともいえるが、各画素について地表面モデル上の候補点の探索と前後関係の計算を行う必要があり、膨大な計算量を要する問題があった。
【0016】
本発明は、上述の課題に鑑みてなされたものであり、処理負荷の増加を抑制しつつ正確なオルソフォトを生成することができるオルソフォト生成方法等を提供することを目的とする。
【0017】
【課題を解決するための手段】
本発明のオルソフォト生成方法は、コンピュータが、
地面オブジェクト及び空中オブジェクトの三次元地図を第1の記憶手段に記憶するステップと、
空中撮影された空中写真画像を第2の記憶手段に記憶するステップと、
前記空中写真画像の内部標定要素及び外部標定要素を含む標定要素を第3の記憶手段に記憶するステップと、
オルソフォト画像を得るための範囲が入力されると、この範囲に対応する前記空中写真画像の標定要素に基づいて、オルソフォト座標とその位置の標高を、前記空中写真画像上の位置に関係を定めるステップと、
前記三次元地図の地面上オブジェクト及び空中オブジェクトの境界線を拘束線としてTINを発生させて三角形を生成し、これらの三角形を基礎三角形群とし、該基礎三角形群の各々の基礎三角形に前記三次元地図に基づいて高さを与えた三角柱モデルを生成するステップと、
前記各々の基礎三角形に、前記空中写真画像の投影中心を投影した二次元視点と前記各々の基礎三角形の各頂点の関係から非隠蔽三角形及び隠蔽三角形を判断し、前記二次元視点から見るとき、前記各々の基礎三角形を、順番が前の基礎三角形が後ろの基礎三角形に隠蔽されないということを保証するように並べるステップと、
前記並べ替えが不可能な場合、これを可能であるように保障するために、前記二次元視点から前記基礎三角形に引いた半直線を境界線とし、この境界線で前記基礎三角形および三角柱モデルを分割して新たな前記基礎三角形を得るステップと、
前記並べられた基礎三角形に対応する三角柱モデル内の三角柱を順に取り出して、この三角柱について前記オルソフォト座標とその位置の標高を、前記空中写真画像上の位置に関係を定めることによって、前記範囲の前記オルソフォト画像を得るステップと
を行うことを要旨とする。
【0018】
本発明のオルソフォト生成システムは、地面オブジェクト及び空中オブジェクトの三次元地図を記憶した第1の記憶手段と、
空中撮影された空中写真画像を記憶した第2の記憶手段と、
前記空中写真画像の内部標定要素及び外部標定要素を含む標定要素を記憶した第3の記憶手段と、
オルソフォト画像を得るための範囲が入力されると、この範囲に対応する前記空中写真画像の標定要素に基づいて、オルソフォト座標とその位置の標高を、前記空中写真画像上の位置に関係を定める手段と、
前記三次元地図の地面上オブジェクト及び空中オブジェクトの境界線を拘束線としてTINを発生させて三角形を生成し、これらの三角形を基礎三角形群とし、該基礎三角形群の各々の基礎三角形に前記三次元地図に基づいて高さを与えた三角柱モデルを生成する手段と、
前記各々の基礎三角形に、前記空中写真画像の投影中心を投影した二次元視点と前記各々の基礎三角形の各頂点の関係から非隠蔽三角形及び隠蔽三角形を判断し、前記二次元視点から見るとき、前記各々の基礎三角形を、順番が前の基礎三角形が後ろの基礎三角形に隠蔽されないということを保証するように並べる手段と、
前記並べ替えが不可能な場合、これを可能であるように保障するために、前記二次元視点から前記基礎三角形に引いた半直線を境界線とし、この境界線で前記基礎三角形および三角柱モデルを分割して新たな前記基礎三角形を得る手段と、
前記並べられた基礎三角形に対応する三角柱モデル内の三角柱を順に取り出して、この三角柱について前記オルソフォト座標とその位置の標高を、前記空中写真画像上の位置に関係を定めることによって、前記範囲の前記オルソフォト画像を得る手段とを備えたことを要旨とする。
【0019】
さらに、オルソフォト生成プログラムは、コンピュータに、
地面オブジェクト及び空中オブジェクトの三次元地図を第1の記憶手段に記憶する手段、
空中撮影された空中写真画像を第2の記憶手段に記憶する手段、
前記空中写真画像の内部標定要素及び外部標定要素を含む標定要素を第3の記憶手段に記憶する手段、
オルソフォト画像を得るための範囲が入力されると、この範囲に対応する前記空中写真画像の標定要素に基づいて、オルソフォト座標とその位置の標高を、前記空中写真画像上の位置に関係を定める手段、
前記三次元地図の地面上オブジェクト及び空中オブジェクトの境界線を拘束線としてTINを発生させて三角形を生成し、これらの三角形を基礎三角形群とし、該基礎三角形群の各々の基礎三角形に前記三次元地図に基づいて高さを与えた三角柱モデルを生成する手段、
前記各々の基礎三角形に、前記空中写真画像の投影中心を投影した二次元視点と前記各々の基礎三角形の各頂点の関係から非隠蔽三角形及び隠蔽三角形を判断し、前記二次元視点から見るとき、前記各々の基礎三角形を、順番が前の基礎三角形が後ろの基礎三角形に隠蔽されないということを保証するように並べる手段、
前記並べ替えが不可能な場合、これを可能であるように保障するために、前記二次元視点から前記基礎三角形に引いた半直線を境界線とし、この境界線で前記基礎三角形および三角柱モデルを分割して新たな前記基礎三角形を得る手段、
前記並べられた基礎三角形に対応する三角柱モデル内の三角柱を順に取り出して、この三角柱について前記オルソフォト座標とその位置の標高を、前記空中写真画像上の位置に関係を定めることによって、前記範囲の前記オルソフォト画像を得る手段
としての機能を実行させる。
【0020】
【発明の実施の形態】
本発明は、例えば航空写真からオルソフォトを生成するオルソフォト生成システムに適用することができる。
【0021】
(構成)
本発明の一実施形態に係るオルソフォト生成システムは、例えばプロセッサ(MPU)、メモリ、ハードディスクドライブ(HDD),光ディスクドライブ等の補助記憶装置、キーボード,ネットワークインターフェース(NIC)等の入力装置、ディスプレイ装置,画像信号出力インターフェース等の出力装置を備えるパーソナルコンピュータ,ワークステーション等の情報処理装置から構成することができる。
【0022】
このような情報処理装置は、例えばMPUが本発明に係るオルソフォト生成プログラムを実行することにより、図1に示すように、3次元地図を保持する3次元地図DB10と、空中撮影された地表画像を保持する空中写真画像DB20と、各々の地表画像毎に標定要素を保持する標定要素DB30と、3次元地図DB10に格納された3次元地図から三角柱モデルを生成する三角柱モデル生成部40と、生成された三角柱モデルと標定要素DB30に格納されている標定要素に応じて空中写真画像DB20に格納されている個々の写真画像をオルソフォト化するオルソフォト化部50と、オルソフォト化された個々の写真画像を格納するオルソフォト化データ格納部60と、オルソフォト化された個々の写真画像を合成する合成処理部70と、合成されたオルソフォト画像を格納するオルソフォト格納部80として機能する。
【0023】
なお、この図1では、本発明のオルソフォト生成装置を1つの装置として構成した例について示しているが、各DB10〜30を各々別の装置として構成する等、適宜構成を変更することができる。また、この図1と同様の構成をハードウェアで構成することもできる。
【0024】
このオルソフォト生成装置には、入力データとして3次元地図,空中写真画像及び標定要素が入力される。これらのデータは、例えば上述の入力装置を介して外部から供給され、各々3次元地図DB10、空中写真画像DB20及び標定要素DB30に格納される。
【0025】
(3次元地図)
3次元地図DB10に格納される3次元地図は、全ての座標が標高付で保存されている地図である。具体的には、3次元の点、もしくは3次元の点を結んだ線分・連続線分・閉曲線を示すデータからなる。この3次元地図は更に、地面上オブジェクトと空中オブジェクトの2つに分けられているものとする。
【0026】
地面上オブジェクトは、道路・等高線・河川等の地面上にあると考えられるオブジェクト(物体)である。これに対し、空中オブジェクトは、建物・橋梁等、図化された部分が地面より上にあるオブジェクトであり、多角形の閉領域で記述されている。
【0027】
(空中写真画像)
空中写真画像DB20には、例えばデジタル形式のni個の空中写真画像I1−Iniが格納される。なお、これらの空中写真画像では、レンズディストーションのような非線形な歪は無視できるか、予め取り除かれているものとする。
【0028】
(標定要素)
標定要素DB30に格納される標定要素は、各写真画像の内部標定要素(焦点距離・画像座標と写真座標の変換行列)及び外部標定要素(投影中心座標・カメラの地上座標系に対する回転角)である。
【0029】
また、このオルソフォト生成装置には、オルソフォトを生成するためのパラメータ(オルソフォトパラメータ)が供給される。このオルソフォトパラメータは、少なくとも次のパラメータを含んでいる。
【0030】
▲1▼オルソフォトの原点(オルソフォト画像原点の中心座標)(X0,Y0)
▲2▼カラム軸及びライン軸方向の単位方向ベクトルec(Dxx,Dxy)・el(Dyx,Dyy)
▲3▼カラム軸及びライン軸方向の解像度dc・dl
▲4▼オルソフォトのカラム数及びライン数Sc・Sl
なお、これらのパラメータによって、オルソフォトの範囲は、例えば図2に示すように、次の4点(X0,Y0),(X1,Y1),(X2,Y2),(X3,Y3)で囲まれた矩形範囲となる。
【0031】
また、このオルソフォト生成装置は、出力データとしてオルソフォトパラメータで指定される範囲のオルソフォト画像を出力する。このオルソフォト画像は、例えばオルソフォト格納部80に格納され、外部からの要求に応じて出力される。
【0032】
(動作)
上述のように構成されたオルソフォト生成装置では、例えば図3に示すような三角柱モデルを用い、オブジェクトの隠蔽関係を考慮してオルソフォトの計算を行う。
【0033】
この三角柱モデルは、同図中に示すように、建物等の構造物を含む地表面を、すべて三角柱で表現したものであり、以下の条件を満たす有限個の三角柱の集合(三角柱群)である。
【0034】
[条件1]
三角柱の側面は鉛直方向に直立し、底面と上面を水平面に鉛直に投影した像は同じ三角形に一致する。この三角形を基礎三角柱と呼ぶ。
【0035】
[条件2]
全ての三角柱の基礎三角形は、境界を除いて互いに疎、つまり共通部分を持たない。これらの三角形を合せて基礎三角柱群と呼ぶ。
【0036】
ここで、三角柱の上面と底面は必ずしも平行でなくとも良い。また、例えば図4に示すように、上面と底面が一致していても良い。つまり、一般の三角メッシュを使った地形モデルも概念的には三角柱モデルに含まれる。
【0037】
また、上述の条件2から、三角柱モデルを構成する三角柱同士も側面を除いて共通部分を持たないことが帰結できる。
【0038】
オルソフォト計算の全体は、例えば図5に示すように、大きく▲1▼座標変換及び地面標高計算準備(S1),▲2▼三角柱モデルの生成(S2),▲3▼各写真のオルソフォト化(S3),▲4▼オルソフォトの合成(S4)の4つの部分からなっている。
【0039】
▲1▼座標変換及び地面標高計算準備(S1)
座標変換及び地面標高計算準備では、例えば図6に示すように、座標変換のパラメータの設定(S11)と、地面標高計算のパラメータの設定(S12)を行う。
【0040】
座標変像のためのパラメータの設定(S11)では、オルソフォト化部50は、次の3つの座標変換を実行するために、オルソフォトパラメータ及び標定要素を読み込む。
【0041】
(i)上述の式(1)及び式(2)による3次元座標(X,Y,Z)から画像座標(xI,yI)への変換
(ii)上述の式(9)によるオルソフォト画像座標(i,j)から平面座標(X,Y)への変換
(iii)上述の式(10)による平面座標(X,Y)からオルソフォト座標(i,j)への変換
また、地面標高計算のためのパラメータの設定(S12)では、後の処理で必要となる地面標高計算のための準備を行う。地面標高計算では、任意の平面座標(X,Y)に対して、その位置における地面の標高を求める。
【0042】
三角柱モデル生成部40に設けられた地形標高計算部は、3次元地図の内、地面上オブジェクトを読み込み、地面上オブジェクトを拘束線としてTINを発生し、地形モデルとする。与えられた(X,Y)のZ座標は、平面座標で(X,Y)を含む三角形の頂点の3次元座標(X1,Y1,Z1),(X2,Y2,Z2),(X3,Y3,Z3)から線形内挿で求める。地面標高計算を、平面座標(X,Y)の関数として、Z=T(X,Y) (13)
と表現する。
【0043】
▲2▼三角柱モデルの生成(S2)
三角柱モデルの生成では、例えば図7に示すように、オルソフォト生成に利用する三角柱モデルを生成する。
【0044】
まず、三角柱モデル生成部40は、3次元地図の読み込みを行う(S21)。具体的には、出力範囲の矩形(X0,Y0)−(X3,Y3)及びその周りL(Lは矩形の長辺の1/20〜1/5程度)の範囲にある地面上オブジェクト及び空中オブジェクトを読み込む。なお、多角形等は、範囲にかかるものはすべて読み出しておく。
【0045】
次に、空中オブジェクト直下の地面標高の計算を行う(S22)。具体的には、空中オブジェクトの各点について、その直下の地面標高を、上述の式(13)による地面標高計算により計算する。
【0046】
この後、TINの生成を行う(S23)。具体的には、空中オブジェクト及び地面上オブジェクトの両方の平面座標より、新たにTINを発生させる。この際、地上オブジェクトの中の線(道路縁)と共に空中オブジェクトの多角形境界も拘束線(ブレークライン)とする。
【0047】
次に、三角柱モデルの生成(S24)を行う。具体的には、S23で生成したTINを基礎三角形とする三角柱を生成し、上述の図3に示す三角柱モデルを構成する。同図中に示すように、空中オブジェクト内部の三角柱をタイプA、それ以外をタイプBと呼ぶ。各々の三角柱は各々、以下の上面三角形,下面三角形及び側面四角形によって定義される。
【0048】
・上面三角形
タイプA:空中オブジェクトの各ノードの標高を持つ三角形とする。
【0049】
タイプB:地面上オブジェクトの標高もしくは、S22で計算した空中オブジェクトノード直下の地上位置標高を持つ三角形とする。
【0050】
・底面三角形
空中オブジェクト及び地面上オブジェクトの最下点標高より低い一定の値Hbottomを持つ三角形とする。
【0051】
・側面四角形
上記上面と側面を結ぶ3つの四角形とする。
【0052】
ここで生成した三角柱群をTpri(i=1,・・・,ntpr)とする。また、各三角柱の上面三角形をTtop(Tpri)、底面三角形をTbottom(Tpri)、側面四角形をTsidej(Tpri)(j=1,2,3)とする。
【0053】
このような三角柱モデルでは、上述の図3に示すように、タイプA・タイプB共に、標高Hbottomに底面を持つ三角柱で表すことができる。このような三角柱モデルを用いると、後述のように、図形の隠蔽関係を単純化できるので、画素毎に前後関係を求める必要がなく、高速にオルソフォトを生成することができる。
【0054】
最後に、後述する隠蔽関係ソート(S43)を保証するために、隠蔽関係ソート可能性の保証のための三角柱モデルの半直線分割(S25)を行う。この三角柱モデルの半直線分割は、隠蔽関係ソート可能性を保証するために行うものであり、具体的には以下の手順で実行される。
【0055】
(半直線分割を利用した隠蔽関係ソート可能性保証)
三角形群の中には、例えば図8の左端の例に示すような、隠蔽関係ソート不可能な例も含まれる。
【0056】
このような場合は、視点から引いた半直線によって三角形群を多角形に分割し(同図中央)、更にその多角形を三角形に分割する(同図右)によって、隠蔽関係ソート可能な三角形群にすることができる。これを、三角形群の半直線分割と呼ぶ。
【0057】
この半直線分割を利用して、三角柱モデルの隠蔽関係ソートの可能性を保証するよう三角柱を分割する方法について述べる。
【0058】
半直線分割を利用すると、例えば図9に示す方法で、三角柱モデルの隠蔽関係ソート可能性が保証された三角柱モデルに分割することができる。
【0059】
・基礎三角形群の取得(S31)
三角柱モデルTpr:={Tpr1,Tpr2,・‥,Tprn}から、基礎三角形群Tri:={Tri1,Tri2,・‥,Trin}を取得する。また、写真画像の投影中心を2次元平面に投影した点(2次元視点)をもとめ、これをOとする。
【0060】
・基礎三角形群の隠蔽関係ソート可能性の判定(S32)
Oに対するTriの隠蔽関係ソートを行い、隠蔽関係ソート可能性を判定する。隠蔽関係ソート可能なら、三角柱モデルTprも隠蔽ソート可能である。
【0061】
・基礎三角形群の半直線分割(S33)
隠蔽関係ソート不可能なら、2次元視点Oから発する1本の半直線によりTriに半直線分割を施す。半直録分割の結果をTrinewとする。
【0062】
・三角柱モデルの分割(S34)
Trinewを基礎三角形群とするよう、三角柱モデルTprの三角柱を細分化した三角柱モデルTprnewを生成する。上述のようにTrinewは隠蔽関係ソート可能であり、Tprnewも隠蔽関係ソート可能であることが保証される。
【0063】
(3)各写真のオルソフォト化(S3)
半直線分割が終了すると、オルソフォト化部50は、例えば図10に示すように、入力された全ての空中写真画像I1〜Ini(個数ni)について、オルソフォトパラメータで指定される範囲のオルソフォトを作成する。この処理においては、作業領域として、サイズがSc×Slのオルソフォト出力用バッファ配列Iorth(i,j)と、写真画像と同じサイズの写真画像格納バッファI(i,j)及びオクルージョンコードバッファ配列Iocl(i,j)を使用する。
【0064】
・範囲のチェック(S41)
オルソフォトの出力範囲と各写真画像が重なるかどうかをチェックする。具体的には、まずオルソフォト出力範囲の4隅(Xi,Yi)(i=0,…,3)の写る画像座標(xIi,yIi)を次式によって求める。
【0065】
xIi=Fx(Xi,Yi,T(Xi,Yi)) (14)
yIi=Fy(Xi,Yi,T(Xi,Yi)) (15)
この(xIi,yIi)で形成する四角形と画像の4隅が構成する矩形が重なり合うかどうかを調べることによって、出力範囲が画像と重なるかどうかをチェックすることができる。
【0066】
チェックの結果、オルソ出力範囲が写真に重ならない場合は、次の写真の処理(S40)に移る。重なる場合は次のS42に進む。
【0067】
・画像の読み込みとバッファの初期化(S42)
写真画像I(i,j)を読み込む。また、オルソフォト出力バッファIorth(i,j)を初期コードOinitに初期化する。Oinitは、画素値としてあり得ない値にしておく。また、オクルージョンコードバッファ配列Iocl(i,j)を非オクルージョンコードMinit(典型的には0)に初期化する。なお、Iocl(i,j)は、非オクルージョンコードMinitとオクルージョンコードMocl(典型的には1)の2値を取るものとする。
【0068】
・三角柱モデルの隠蔽関係ソート(S43)
三角柱群Tpri(i=1,…,ntpr)を、写真の投影中心座標を視点とした隠蔽関係によって予めソートする。
【0069】
隠蔽関係によってソートするとは、三角柱(若しくは三角形)をある3次元的視点(三角形のときは2次元的視点)から見るとき、順番が前の物が後ろの物に隠蔽される事がないということを保証するように並べることである。このように並べられた三角柱(若しくは三角形)を隠蔽関係ソート列と呼ぶ。
【0070】
また、与えられた三角柱群(若しくは三角系群)の全てを隠蔽関係ソート列に構成する事ができる場合、隠蔽関係ソート可能と呼ぶ。なお、ここで使う隠蔽という言葉を数学的に定義しておく。
【0071】
例えば図11(A)及び同図(B)に示す状態において、視点Oからみて三角柱(三角形)Aが他の三角柱(三角形)Bに隠蔽されるとは、Aの内部の点PAとOを結ぶ線分PAO上に、Bの内部の点PBが存在する事を言う。逆に隠蔽されないとは、そのようなBの内部の点が存在しない事を言う。
【0072】
三角柱(三角形)Aが視点Oに対して三角柱(三角形)Bに隠蔽されることを、次の式で表現することとする。
【0073】
【数5】
また、Aが三角柱(三角形)Bに隠蔽されない
【数6】
ことを、次の式で表現することとする。
【0074】
【数7】
三角柱モデルと基礎三角形群については、次の重要な定理が存在する。
【0075】
[定理]
三角柱モデルTprの基礎三角形群Triが隠蔽関係ソート可能ならば、Tprも隠蔽関係ソート可能である。また、STriがTriの隠蔽関係ソート列ならば、STriの各三角形に対応する三角柱をSTriと同じ順に並べた三角柱列STprは、Tprの隠蔽関係ソート列である。つまり、三角柱モデルの隠蔽関係ソートを行うには、その基礎三角形群の隠蔽関係ソート列を行えばよい。以下、基礎三角形群の隠蔽関係ソートを利用した三角柱モデルの隠蔽関係ソート方法について説明する。
【0076】
上述の三角柱の隠蔽関係ソートは、例えば図12に示す処理によって行う。
【0077】
・基礎三角形群の取得(S51)
三角柱モデルTpr:={Tpr1,Tpr2,…,Tprn}から、基礎三角形群Tri:={Tri1,Tri2,・・・,Trin}を取得する。また、画像の投影中心を水平面に投影したものを2次元視点Oとする。
【0078】
・三角形分割の取得(S52)
例えば図13に示すように、与えられた基礎三角形群Triの頂点より、三角形の辺をブレークラインとして、全体が凸包で囲まれるよう三角形分割を行う(同図(1))。これは普通TINを発生することによって得られる。なお、視点が三角形群の外部にある場合は視点もTINを構成する点に追加しておく(同図(2))。ここで得られた三角形分割の三角形群をTri′={Tri1′,Tri2′,・・・,Trin′′}とする。
【0079】
・最初の三角形の取り出し(S53)
まず、図14(1)に示すように、三角形分割領域内で2次元視点Oを含む三角形を隠蔽関係ソート列の最初の三角形S1とする。視点Oが三角形境界上にあるときは、その境界を構成するいずれかの三角形をS1とする。
【0080】
・前線多角形の初期化(S54)
視点の周囲で、既に取り出した三角形内部を除いた多角形境界を前線多角形と呼ぶ(図14(2)の太線部)。最初の前線多角形は、上記STri1′それ自身である。
【0081】
・非隠蔽三角形のサーチ(S55)
前線多角形の辺を1辺とする三角形の中で、非隠蔽三角形であるかどうかを調べる。非隠蔽三角形である条件は、図15に示すように、次の2つの条件のいずれかである。
【0082】
[条件1]
前線多角形の辺をAB、この辺を1辺とする三角形の頂点Cとすると、線分ABと線分OCが交点を持つ(図15(1))。
【0083】
[条件2]
前線多角形の辺をAB、この辺を1辺とする三角形の頂点Cとすると、Cが前線多角形のAもしくはBの隣の点である(図15(2))。
【0084】
・隠蔽関係ソート列の追加(S57)
S56において、非隠蔽三角形が見つかったときは、その三角形を隠蔽関係ソート列Siに追加する。
【0085】
・前線多角形の更新(S58)
ソート列に追加した三角形を三角形群から取り除く。すなわち、前線多角形を結合した図形の外周を新たに前線多角形とする(上述の図14中の(3)〜(11))。ソート対象となる三角形が残っている場合はS55に戻る。
【0086】
・処理の終了(S59)
全ての前線多角形の辺について非隠蔽三角形が見つからない、もしくはソート対象となる三角形がなくなったときは、処理を終了する(図14(11))。得られた隠蔽関係ソート列をS:={S1,S2,…,Sn′}とする。
【0087】
・隠蔽関係ソート列の抽出
上述の三角形分割の隠蔽関係ソート(S51〜S59)で得られた隠蔽関係ソート列Sから、元の三角形群に含まれるもののみを抽出し、これをSTri:={STri1,STri2,…,STrin}とする。
【0088】
・三角柱の隠蔽関係ソート列の取得(S60)
STriの各三角形を基礎三角形とする三角柱をSTriと同じ順に取り出したものをSTpr:={STpr1,STpr2,・‥,STprn}とすると、これが求めるTprの隠蔽関係ソート列である。
【0089】
上述の処理によって隠蔽関係ソートの順に並べられた三角柱群を改めてTpri(i=1,…,ntpr)とし、図10中のS44に進む。
【0090】
・三角柱のオルソフォト化(S44)
三角柱を隠蔽関係ソートソート順に順次オルソフォト化する。三角柱Tpriを順次取り出し、以下に示すS45〜S48の処理を実行する。
【0091】
S45では、順次、次の三角柱を選択する。
【0092】
(i)オルソフォトの生成(S46)
オルソフォト化部50に設けられた座標変換モジュールにより三角柱の基礎三角形Tbase(Tpri)をオルソフォト画像座標の三角形Torth(Tpri)に変換する。また、オルソフォト化部50に設けられたオルソフォト出力バッファで、上述の三角形Torth(Tpri)の内部及び周上の各画素について、上述の式(11)によって、Iorth(i,j)を求める。なお、式(11)に用いる地表面モデルZ=S(X,Y)には、上面三角形Ttop(Tpri)を通る平面上の標高を与える関数を用いる。
【0093】
具体的には、次の手順で処理を実行する。
【0094】
▲1▼式(9)により、Torth(Tpri)内の点(i,j)について、対応する平面座標(X,Y)を算出する。
【0095】
▲2▼上述の式(1)及び式(2)により、平面座標(X,Y)に対応する写真画像上の座標(xI,yI)を求める。この座標を四捨五入して得られる整数座標を(ixI,iyI)としたとき、オクルージョンコードバッファIoci(ixI,iyI)がMinitならば、上述の(11)式によってIorth(i,j)を計算する。また、Moclならば、Iorth(i,j)は初期コードOinitのままとする。
【0096】
(ii)オクルージョンの記録(S47)
三角柱Tpriによってオクルージョンとなる部分を記録しておく。写真画像上でこの構造物によって隠される部分は、この構造物の上面三角形Ttop(Tpri))、底面三角形Tbottom(Tpri)及び3個の側面四角形Tsidej(Tpri)(j=1,2,3)が投影される部分である。
【0097】
オクルージョンの記録にあたっては、Ttop(Tpri)・Tbottom(Tpri)・Tsidej(Tpri)を順次上述の式(1)及び式(2)によってデジタル写真画像座標に変換する。変換後の2次元画像座標多角形の内部にある全ての整数座標(ii,jj)についてIocl(ii,jj)をオクルージョンコードMoclに設定し、オクルージョンとなっていることを記録する。
【0098】
以上の処理を実行した後、S48において三角柱が残っているか否かを判定し、三角柱が残っていなければS49に進み、三角柱が残っていれば、S45に戻る。
【0099】
・オルソフォトの格納(S49)
全ての三角柱についてオルソフォトを生成した後、最終的にオルソフォトIorthをオルソフォト化データ格納部60にIorthjとして格納する。
【0100】
上述のような処理によるオルソフォトの生成過程を示したものが、図16,図17,図18及び図19である。
【0101】
▲1▼三角柱モデルを構成する三角柱の中で、他の三角柱に隠されることがない三角柱(非隠蔽三角柱)を取り出す。例えば視点の直下にある三角柱(例えば図18(1)の番号1の三角柱)は、他の三角柱に隠されることはない。
【0102】
▲2▼取り出した三角柱を水平面に鉛直に投影した三角形(これを基礎三角形と呼ぶ。)の周上及び内部について上述の式(11)によってオルソフォトヘの出力値を計算する(図18(2)のIorth)。この際、地表面モデルZ=S(X,Y)には、三角柱の上面を使用する(図17)。
【0103】
▲3▼この三角柱によって隠される写真画像上の領域を記録しておく(図18(2)のIocl)。
【0104】
▲4▼上記三角柱を三角柱モデルから取り除く(図18(2)の三角柱モデル)。
【0105】
▲5▼三角柱モデルを構成する三角柱(既に処理したもの以外)の中で、他の三角柱に隠されることがない三角柱を取り出す(例えば図18(1)の三角柱2,3,4,…の順)。
【0106】
▲6▼▲2▼と同様に取り出した三角柱の基確三角形の図上及び内部についてオルソフォトを出力する(図18(3),図19(1)〜(2)のIorth)。ただし、既にオルソフォトを作成した三角柱によって隠されていると記録されている写真画像上の部分に対応した場合は、オルソフォト出力バッファに記録しない(図19(1)のIorth参照)。
【0107】
▲7▼この三角柱によって隠される写真画像上の領域を記録しておく(図18(3),図19(1)〜(2)のIocl)。
【0108】
▲8▼上記三角柱を三角柱モデルから取り除く(図18(3),図19(1)〜(2)の三角柱モデル)。
【0109】
▲9▼三角柱モデルに三角柱が残っていれば▲5▼に戻る。残っていなければ処理を終了する。
【0110】
以上の手順によってできたオルソフォトは、オクルージョンによる影部分が欠損したものになる(図20左)。この図20左の図は、遮蔽関係を考慮したオルソフォトであり、黒い部分はオクルージョンによる欠損を示している。
【0111】
しかし、ある写真画像から生成したオルソフォトで欠損している部分は、別の視点から撮影した写真から生成したオルソフォトでは欠損していない可能性がある。そこで複数の写真で作成したオルソフォトを更に合成することにより、オクルージョンを補い合う。こうして、光学的に正しく、かつオクルージョンによる欠損の少ない精密なオルソフォト画像を作成することができる(図20右)。
【0112】
(4)オルソフォトの合成(S4)
オルソフォトの合成では、合成処理部70は、上述の各写真のオルソフォト化(S3)で生成され、オルソフォト化データ格納部60に格納されたオルソフォトをメモリ上に読み出す。ここで生成されたオルソフォトをIorthk(k=1,…,nort)とする。
【0113】
▲1▼投影中心のオルソ画像座標系への変換
Iorthkの原画像の投影中心座標(これは標定要素から既知である)を、上述の式(10)によりオルソ画像座標系に変換したものをOk=(Xok,Yok)とする。
【0114】
▲2▼オルソフォトの合成
最終出力オルソフォトバッファ(オルソフォト格納部80)を
【数8】
とする。合成処理部70は、
【数9】
上の全ての座標(i,j)について、オルソ画像出力値
【数10】
を次の重み付平均によって計算する。
【0115】
【数11】
ここで、w(i,j,k)は画像座標(i,j)に対するk番目の画像の重み関数である。実際の各画像の重みW(i,j,k)は、重み関数を使って次式で表される。
【0116】
【数12】
w(i,j,k)の取り方により、様々なオルソフォトの接合を行うことができるが、いずれの場合も、
Iorthk(i,j)=Oinitのとき、w(i,j,k)=0 (18)とする。つまり、各オルソフォト内で、初期値のままの部分については重み0である。また、式(16)よりわかるように、全てのオルソフォトについてIorthk(i,j)=Oinitなら、出力は背景色Ibackになる。また、一つの写真でオクルージョンのためIorthk(i,j)=Oinitとなっていても、他の写真でIorthm(i,j)≠Oinitならば、上述の図20に示すように、オクルージョン部分を補うことができる。
【0117】
一般に投影中心が近い写真を優先して画像を合成するので、重み関数w(k,i,j)としてOkと(i,j)の距離D(i,j,k)を利用する最近傍法,最近傍法の境界付近で画像の合成を滑らかに切り替える方法等によって合成を行うことが多い。
【0118】
最近傍法は、最も主点位置が近い画像を優先する方法で、重み関数w(k,i,j)として、例えば次式のような関数を用いる。
【0119】
【数13】
最近傍法の境界付近で画像の合成を滑らかに切り替える方法では、w(k)が距離に応じて急激に小さくなる関数、例えば次のような関数を用いる。
【0120】
【数14】
ここでpは適当な整数とする。pが大きいほど、境界部分での画像の重みの切り替えが急激になる(図21・図22・図23)。
【0121】
(効果)
以上説明したように、この実施形態に係るオルソフォト生成装置では、三角柱モデルを用いることにより、隠蔽関係を判定することが容易になり、詳細なオルソフォトを高速に生成することができる。
【0122】
なお、本発明の技術的思想は、上述の実施形態に限定されず、適宜変更することができる。
【0123】
【発明の効果】
本発明に係るオルソフォト生成システムでは、三角柱モデルを用いることにより、隠蔽関係を判定することが容易になり、処理負荷の増加を抑制しつつ正確なオルソフォトを生成することができる。これにより、従来はあまり利用されなかった都市部の精密オルソフォトの生成が可能になる。
【図面の簡単な説明】
【図1】本発明の一実施形態に係るオルソフォト生成システムの構成例を示すブロック図である。
【図2】オルソフォトの範囲の例を示す図である。
【図3】前記オルソフォト生成システムで用いる三角柱モデルを示す図である。
【図4】三角柱モデルと基礎三角形群を示す図である。
【図5】前記オルソフォト生成システムによるオルソフォトの計算処理の例を示すフローチャートである。
【図6】前記オルソフォト生成システムを構成するオルソフォト化部によるパラメータの設定処理の例を示すフローチャートである。
【図7】前記オルソフォト生成システムを構成する三角柱モデル生成部による三角柱モデルの生成処理の例を示すフローチャートである。
【図8】前記三角柱モデル生成部による半直線分割処理の例を示すフローチャートである。
【図9】前記三角柱モデル生成部による三角柱モデルの隠蔽関係ソート可能性が保証された三角柱モデルに分割する処理の例を示すフローチャートである。
【図10】前記オルソフォト化部によるオルソフォト化処理の例を示すフローチャートである。
【図11】三角柱と三角形の隠蔽を示す図である。
【図12】前記オルソフォト化部による三角柱の隠蔽関係ソート処理の例を示すフローチャートである。
【図13】前記オルソフォト化部による三角柱の分割処理の例を示す図である。
【図14】前記オルソフォト化部による三角柱分割の隠蔽関係ソート手順を示す図である。
【図15】非隠蔽三角形の条件を示す図である。
【図16】三角柱モデルを用いたオルソフォト生成過程を示す図である。
【図17】三角柱モデルを用いたオルソフォト生成過程を示す図である。
【図18】三角柱モデルを用いたオルソフォト生成過程を示す図である。
【図19】三角柱モデルを用いたオルソフォト作生成過程を示す図である。
【図20】オルソフォト生成方法を示す図である。
【図21】合成する2枚の画像の位置関係を示す図である。
【図22】画像の重み関数の例を示す図である。
【図23】画像の重み関数の例を示す図である。
【図24】オルソフォトの原理を示す図である。
【図25】各座標系の関係を示す図である。
【図26】平面座標系とオルソフォト座標系の関係を示す図である。
【図27】従来のオルソフォトの問題点を示す図である。
【図28】従来のオルソフォトの問題点を示す図である。
【符号の説明】
10…3次元地図DB、
20…空中写真画像DB、
30…標定要素DB、
40…三角柱モデル生成部、
50…オルソフォト化部、
60…オルソフォト化データ格納部、
70…合成処理部、
80…オルソフォト格納部
Claims (4)
- コンピュータが、
地面オブジェクト及び空中オブジェクトの三次元地図を第1の記憶手段に記憶するステップと、
空中撮影された空中写真画像を第2の記憶手段に記憶するステップと、
前記空中写真画像の内部標定要素及び外部標定要素を含む標定要素を第3の記憶手段に記憶するステップと、
オルソフォト画像を得るための範囲が入力されると、この範囲に対応する前記空中写真画像の標定要素に基づいて、オルソフォト座標とその位置の標高を、前記空中写真画像上の位置に関係を定めるステップと、
前記三次元地図の地面上オブジェクト及び空中オブジェクトの境界線を拘束線としてTINを発生させて三角形を生成し、これらの三角形を基礎三角形群とし、該基礎三角形群の各々の基礎三角形に前記三次元地図に基づいて高さを与えた三角柱モデルを生成するステップと、
前記各々の基礎三角形に、前記空中写真画像の投影中心を投影した二次元視点と前記各々の基礎三角形の各頂点の関係から非隠蔽三角形及び隠蔽三角形を判断し、前記二次元視点から見るとき、前記各々の基礎三角形を、順番が前の基礎三角形が後ろの基礎三角形に隠蔽されないということを保証するように並べるステップと、
前記並べ替えが不可能な場合、これを可能であるように保障するために、前記二次元視点から前記基礎三角形に引いた半直線を境界線とし、この境界線で前記基礎三角形および三角柱モデルを分割して新たな前記基礎三角形を得るステップと、
前記並べられた基礎三角形に対応する三角柱モデル内の三角柱を順に取り出して、この三角柱について前記オルソフォト座標とその位置の標高を、前記空中写真画像上の位置に関係を定めることによって、前記範囲の前記オルソフォト画像を得るステップと
を行うことを特徴とするオルソフォト生成方法。 - 前記範囲を更新し、該更新毎の範囲で求められた各々の前記オルソフォト画像を合成するステップと
を行うことを特徴とする請求項1記載のオルソフォト生成方法。 - 地面オブジェクト及び空中オブジェクトの三次元地図を記憶した第1の記憶手段と、
空中撮影された空中写真画像を記憶した第2の記憶手段と、
前記空中写真画像の内部標定要素及び外部標定要素を含む標定要素を記憶した第3の記憶手段と、
オルソフォト画像を得るための範囲が入力されると、この範囲に対応する前記空中写真画像の標定要素に基づいて、オルソフォト座標とその位置の標高を、前記空中写真画像上の位置に関係を定める手段と、
前記三次元地図の地面上オブジェクト及び空中オブジェクトの境界線を拘束線としてTINを発生させて三角形を生成し、これらの三角形を基礎三角形群とし、該基礎三角形群の各々の基礎三角形に前記三次元地図に基づいて高さを与えた三角柱モデルを生成する手段と、
前記各々の基礎三角形に、前記空中写真画像の投影中心を投影した二次元視点と前記各々の基礎三角形の各頂点の関係から非隠蔽三角形及び隠蔽三角形を判断し、前記二次元視点から見るとき、前記各々の基礎三角形を、順番が前の基礎三角形が後ろの基礎三角形に隠蔽されないということを保証するように並べる手段と、
前記並べ替えが不可能な場合、これを可能であるように保障するために、前記二次元視点から前記基礎三角形に引いた半直線を境界線とし、この境界線で前記基礎三角形および三角柱モデルを分割して新たな前記基礎三角形を得る手段と、
前記並べられた基礎三角形に対応する三角柱モデル内の三角柱を順に取り出して、この 三角柱について前記オルソフォト座標とその位置の標高を、前記空中写真画像上の位置に関係を定めることによって、前記範囲の前記オルソフォト画像を得る手段と
を有することを特徴とするオルソフォト生成システム。 - コンピュータに、
地面オブジェクト及び空中オブジェクトの三次元地図を第1の記憶手段に記憶する手段、
空中撮影された空中写真画像を第2の記憶手段に記憶する手段、
前記空中写真画像の内部標定要素及び外部標定要素を含む標定要素を第3の記憶手段に記憶する手段、
オルソフォト画像を得るための範囲が入力されると、この範囲に対応する前記空中写真画像の標定要素に基づいて、オルソフォト座標とその位置の標高を、前記空中写真画像上の位置に関係を定める手段、
前記三次元地図の地面上オブジェクト及び空中オブジェクトの境界線を拘束線としてTINを発生させて三角形を生成し、これらの三角形を基礎三角形群とし、該基礎三角形群の各々の基礎三角形に前記三次元地図に基づいて高さを与えた三角柱モデルを生成する手段、
前記各々の基礎三角形に、前記空中写真画像の投影中心を投影した二次元視点と前記各々の基礎三角形の各頂点の関係から非隠蔽三角形及び隠蔽三角形を判断し、前記二次元視点から見るとき、前記各々の基礎三角形を、順番が前の基礎三角形が後ろの基礎三角形に隠蔽されないということを保証するように並べる手段、
前記並べ替えが不可能な場合、これを可能であるように保障するために、前記二次元視点から前記基礎三角形に引いた半直線を境界線とし、この境界線で前記基礎三角形および三角柱モデルを分割して新たな前記基礎三角形を得る手段、
前記並べられた基礎三角形に対応する三角柱モデル内の三角柱を順に取り出して、この三角柱について前記オルソフォト座標とその位置の標高を、前記空中写真画像上の位置に関係を定めることによって、前記範囲の前記オルソフォト画像を得る手段
としての機能を実行させるためのオルソフォト生成プログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002118033A JP4177017B2 (ja) | 2002-04-19 | 2002-04-19 | オルソフォト生成方法、オルソフォト生成システム及びオルソフォト生成プログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002118033A JP4177017B2 (ja) | 2002-04-19 | 2002-04-19 | オルソフォト生成方法、オルソフォト生成システム及びオルソフォト生成プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2003316256A JP2003316256A (ja) | 2003-11-07 |
JP4177017B2 true JP4177017B2 (ja) | 2008-11-05 |
Family
ID=29535058
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002118033A Expired - Fee Related JP4177017B2 (ja) | 2002-04-19 | 2002-04-19 | オルソフォト生成方法、オルソフォト生成システム及びオルソフォト生成プログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4177017B2 (ja) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4642431B2 (ja) * | 2003-11-13 | 2011-03-02 | パナソニック株式会社 | 地図表示装置、地図表示システム、地図表示方法およびプログラム |
JP5137448B2 (ja) * | 2007-04-17 | 2013-02-06 | 飛島建設株式会社 | トンネル切羽面の測量システム |
JP5088627B2 (ja) * | 2008-04-30 | 2012-12-05 | 国際航業株式会社 | マルチマッチング方法 |
JP6121746B2 (ja) * | 2013-02-20 | 2017-04-26 | Necソリューションイノベータ株式会社 | 撮影位置決定装置、プログラム、および方法 |
JP6486069B2 (ja) * | 2014-10-31 | 2019-03-20 | 株式会社東芝 | 画像処理装置、検査装置、検査方法及び画像処理プログラム |
JP7344488B2 (ja) * | 2021-12-13 | 2023-09-14 | アイサンテクノロジー株式会社 | 三次元モデルデータ生成装置、コンピュータプログラム、及び三次元モデルデータ生成方法 |
-
2002
- 2002-04-19 JP JP2002118033A patent/JP4177017B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2003316256A (ja) | 2003-11-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7162933B2 (ja) | オブジェクトの内部空間モデルを確立するための方法、装置及びシステム、並びに、コンピュータ装置及びコンピュータ可読記憶媒体 | |
US6466892B2 (en) | Method of composing three-dimensional multi-viewpoints data | |
CN105026997B (zh) | 投影系统、半导体集成电路及图像修正方法 | |
JP4065488B2 (ja) | 3次元画像生成装置、3次元画像生成方法及び記憶媒体 | |
ES2351961T3 (es) | Procedimiento basado en imágenes de representación y reproducción de objetos tridimensionales. | |
JP6201476B2 (ja) | 自由視点画像撮像装置およびその方法 | |
US20020113865A1 (en) | Image processing method and apparatus | |
JP4025442B2 (ja) | 三次元モデル変換装置及び方法 | |
US7129943B2 (en) | System and method for feature-based light field morphing and texture transfer | |
KR20090031113A (ko) | 다수의 2차원 실사 영상들을 이용하여 3차원 영상을생성하는 방법 및 장치 | |
CN110246161B (zh) | 一种360度全景图像无缝拼接的方法 | |
JP2000268179A (ja) | 三次元形状情報取得方法及び装置,二次元画像取得方法及び装置並びに記録媒体 | |
KR20220053332A (ko) | 파노라마 영상으로부터 공간 모델을 생성하는 서버, 방법 및 컴퓨터 프로그램 | |
US20060066613A1 (en) | Method and system for partitioning the surface of a three dimentional digital object model in order to map a texture | |
JP2003323640A (ja) | レーザスキャナデータと空中写真画像を用いた高精度都市モデルの生成方法及び高精度都市モデルの生成システム並びに高精度都市モデルの生成のプログラム | |
JPH11504452A (ja) | 2次元投影図に基づいて3次元の対象物を再現し、取り扱うための装置と方法 | |
JP4177017B2 (ja) | オルソフォト生成方法、オルソフォト生成システム及びオルソフォト生成プログラム | |
JP3971714B2 (ja) | 仮想視点画像生成方法及び仮想視点画像生成装置、ならびに仮想視点画像生成プログラム及び記録媒体 | |
JP2001243497A (ja) | 3次元画像処理方法,3次元モデリング方法および3次元画像処理プログラムを記録した記録媒体 | |
JP6802928B2 (ja) | 画像生成装置および画像生成方法 | |
JP2002109518A (ja) | 三次元形状復元方法及びシステム | |
JP2003323603A (ja) | ステレオマッチング方法、3次元計測方法及び3次元計測装置並びにステレオマッチング方法のプログラム及び3次元計測のプログラム | |
JP2002520969A (ja) | 動き画像からの自動化された3次元シーン走査方法 | |
JP2002157603A (ja) | 画像処理装置及びその方法並びに記憶媒体 | |
Fusiello | Image-based rendering |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050413 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070828 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20071029 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080401 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080602 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20080805 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080821 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110829 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110829 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120829 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120829 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130829 Year of fee payment: 5 |
|
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 |
|
LAPS | Cancellation because of no payment of annual fees |