Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP3915664B2 - Image processing apparatus, image processing method used therefor, and program therefor - Google Patents
[go: Go Back, main page]

JP3915664B2 - Image processing apparatus, image processing method used therefor, and program therefor - Google Patents

Image processing apparatus, image processing method used therefor, and program therefor Download PDF

Info

Publication number
JP3915664B2
JP3915664B2 JP2002323250A JP2002323250A JP3915664B2 JP 3915664 B2 JP3915664 B2 JP 3915664B2 JP 2002323250 A JP2002323250 A JP 2002323250A JP 2002323250 A JP2002323250 A JP 2002323250A JP 3915664 B2 JP3915664 B2 JP 3915664B2
Authority
JP
Japan
Prior art keywords
image
ortho
dimensional data
images
stereo
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
Application number
JP2002323250A
Other languages
Japanese (ja)
Other versions
JP2004046777A (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.)
NEC Corp
Original Assignee
NEC 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 NEC Corp filed Critical NEC Corp
Priority to JP2002323250A priority Critical patent/JP3915664B2/en
Publication of JP2004046777A publication Critical patent/JP2004046777A/en
Application granted granted Critical
Publication of JP3915664B2 publication Critical patent/JP3915664B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Processing Or Creating Images (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Description

【0001】
【発明の属する技術分野】
本発明は画像処理装置及びそれに用いる画像処理方法並びにそのプログラムに関し、特にステレオ画像からの3 次元データ生成及び正射影変換に用いる画像処理方法に関する。
【0002】
【従来の技術】
従来、ディジタル写真測量や地図作成作業等においては、地上に設置したカメラから得られるステレオ画像、あるいは衛星や航空機から撮影される画像を基に3次元データ[以下、DEM(Digital Elevation Map)データとする]をステレオマッチングによって生成し、それを用いて画像を正射変換してオルソ画像を生成し、図面化等に応用する方法が広く行われている。
【0003】
ここで、オルソ画像とは、通常のカメラの場合、図11(a)に示すように、中心投影のために傾いて撮影される対象を、図11(b)に示すように、所定の面に対して正射影して得られる画像のことであり、複数の画像を連結して広範囲の対象を把握し、図面等を作成するためには不可欠な画像である。
【0004】
上記の手法を応用した例としては、地上でのディジタル画像計測を行う際に上記の手法を用いる技術がある(例えば、特許文献1,2参照)。また、上記の技術を航空写真や衛星画像に応用したシステム例としは、航空写真や衛星画像をステレオマッチングして正射変換する技術がある(例えば、特許文献3参照)。
【0005】
図12は従来の画像処理装置の構成例を示すブロック図である。この図12に示すように、従来の画像処理装置5は複数の画像に対してステレオマッチング処理を行ってDEMデータを抽出するステレオ処理手段11と、DEMデータを用いて画像を正射変換するオルソ化手段13とを備えている。
【0006】
上記のような構成を持つ画像処理装置5では、図13に示すような動作が行われている。すなわち、複数の写真から得られる画像データがステレオ処理手段11に与えられ、従来のシステムと同様に、ステレオマッチングを行うことによって自動的にDEMデータを抽出する(図13ステップS41)。
【0007】
これによって得られるDEMデータは画像データとともに、オルソ化手段13に与えられ、DEM情報と標定された撮影角度情報とに基づいて画像内の画素を移動させることによってオルソ画像が生成される(図13ステップS42)。
【0008】
【特許文献1】
特開2000−182059号公報(第3〜7頁、図1,2)
【特許文献2】
特開平11−351865号公報(第4,5頁、図1)
【特許文献3】
特開2000−348160号公報(第2,3頁、図1)
【0009】
【発明が解決しようとする課題】
しかしながら、上述した従来のオルソ画像生成手法では、ステレオマッチングによって得られたDEMデータに、対応点のミスマッチングによって生じた誤差が含まれており、正射変換の際に誤ったDEMデータを持つ画素が誤った位置に移されてしまうため、上記の方法によって生成されたオルソ画像の画質が低いという問題がある。
【0010】
また、従来のオルソ画像生成手法では、市街地を撮影した航空画像等、画像中の撮影対象が不連続な形状を多く含む場合、ビル等によって隠れた道路等の隠蔽領域で画素が撮像されていないため、オルソ画像中に対応する画素を移すことができず、図14に示すように、隠蔽領域に画素の欠損した抜け領域が生じるので、上記の方法によって生成されたオルソ画像に画素が欠損した領域を多く含むこととなり、不自然な画像になるという問題がある。
【0011】
従来のオルソ画像生成手法では、上述したような問題があるため、作図やCG(Computer Graphics)生成に活用することができ、より自然なオルソ画像が得られる処理方法が望まれている。
【0012】
そこで、本発明の目的は上記の問題点を解消し、DEMデータに含まれるノイズの影響を受けにくく、また欠損領域に起因する違和感の少ない、より自然なオルソ画像を得ることができる画像処理装置及びそれに用いる画像処理方法並びにそのプログラム画像処理装置及びそれに用いるオルソ画像自動整形方法並びにそのプログラムを提供することにある。
【0013】
【課題を解決するための手段】
本発明による画像処理装置は、複数の画像に対してステレオマッチング処理を行って3次元データを生成する画像処理装置であって、前記ステレオマッチング処理によって前記3次元データを抽出するステレオ処理手段と、前記ステレオ処理手段で抽出された前記3次元データに平滑化フィルタを適用して前記3次元データのノイズを除去して平滑化するDEM平滑化手段と、前記DEM平滑化手段で平滑化された3次元データを用いて前記画像に含まれる全画素個別に正射変換するオルソ化手段とを備えている。
【0014】
本発明による他の画像処理装置は、複数の画像に対してステレオマッチング処理を行って3次元データを生成する画像処理装置であって、前記ステレオマッチング処理によって前記3次元データを抽出するステレオ処理手段と、前記3次元データを用いて前記画像を正射変換するオルソ化手段と、前記画像を正射投影面に対する標定角にしたがって射影変換する射影変換手段と、前記オルソ化手段で得られたオルソ画像と前記射影変換手段で得られた射影変換画像とを合成して欠損した画素が補間された補間オルソ画像を生成する画像合成手段とを備えている。
【0015】
本発明による画像処理方法は、画像処理装置において複数の画像に対してステレオマッチング処理を行って3次元データを生成する画像処理方法であって、前記画像処理装置に、前記ステレオマッチング処理によって前記3次元データを抽出するステップと、その抽出された前記3次元データに平滑化フィルタを適用して前記3次元データのノイズを除去して平滑化するステップと、平滑化された3次元データを用いて前記画像に含まれる全画素個別に正射変換するステップとを備えている。
【0016】
本発明による他の画像処理方法は、複数の画像に対してステレオマッチング処理を行って3次元データを生成する画像処理方法であって、前記ステレオマッチング処理によって前記3次元データを抽出するステップと、前記3次元データを用いて前記画像を正射変換するステップと、前記画像を正射投影面に対する標定角にしたがって射影変換するステップと、前記画像の正射変換で得られたオルソ画像と前記射影変換で得られた射影変換画像とを合成して欠損した画素が補間された補間オルソ画像を生成するステップとを備えている。
【0017】
本発明による画像処理方法のプログラムは、複数の画像に対してステレオマッチング処理を行って3次元データを生成する画像処理方法のプログラムであって、コンピュータに、前記ステレオマッチング処理によって前記3次元データを抽出する処理と、その抽出された前記3次元データに平滑化フィルタを適用して前記3次元データのノイズを除去して平滑化する処理と、平滑化された3次元データを用いて前記画像に含まれる全画素個別にを正射変換する処理とを実行させている。
【0018】
本発明による他の画像処理方法のプログラムは、複数の画像に対してステレオマッチング処理を行って3次元データを生成する画像処理方法のプログラムであって、コンピュータに、前記ステレオマッチング処理によって前記3次元データを抽出する処理と、前記3次元データを用いて前記画像を正射変換する処理と、前記画像を正射投影面に対する標定角にしたがって射影変換する処理と、前記画像の正射変換で得られたオルソ画像と前記射影変換で得られた射影変換画像とを合成して欠損した画素が補間された補間オルソ画像を生成する処理とを実行させている。
【0019】
すなわち、本発明の画像処理装置は、複数の画像に対してステレオマッチング処理を行って3次元データ[以下、DEM(Digital Elevation Map)データとする]を抽出するステレオ処理手段と、DEMデータに平滑化フィルタを適用してノイズを除去するDEM平滑化手段と、DEMデータを用いて画像を正射変換するオルソ化手段と、画像を正射投影面に対する標定角にしたがって射影変換する射影変換手段と、射影変換画像とオルソ画像とを合成して補間オルソ画像を生成する画像合成手段とを備えている。
【0020】
本発明の画像処理装置では上記のような構成とすることで、複数の画像に対してステレオマッチング処理を自動的に行い、得られた視差情報による画像中の各点における高さ情報、すなわちDEMデータを抽出し、平滑化フィルタを適用してDEMデータ中のノイズを除去した後、画像を平滑化したDEMデータを用いて正射変換してオルソ画像を生成し、さらに画像を正射投影面と平行な面に対して射影変換した射影変換画像と合成し、欠損した画素を補間した補間オルソ画像を生成するよう動作している。
【0021】
これによって、本発明の画像処理装置では、DEMデータ中のノイズを低減してオルソ化に利用し、また隠蔽による画素の欠損領域を射影変換画像で補間することで、ノイズによる画質低下や欠損領域による違和感の少ない、より自然なオルソ画像を生成することが可能となる。
【0022】
【発明の実施の形態】
次に、本発明の実施例について図面を参照して説明する。図1は本発明の第1の実施例による画像処理装置の構成を示すブロック図である。図1において、画像処理装置1は複数の画像(航空画像データa)に対してステレオマッチング処理を行って3次元データ[以下、DEM(Digital ElevationMap)データとする]bを抽出するステレオ処理手段11と、抽出されたDEMデータbに平滑化フィルタ(図示せず)を適用してノイズを除去するDEM平滑化手段12と、平滑化DEMデータcを用いて画像(航空画像データa)を正射変換するオルソ化手段13と、画像処理装置1の各手段を実現するためのプログラムを格納する記録媒体14とから構成されている。
【0023】
上記の画像処理装置1の各手段はそれぞれ概略、次のように動作する。つまり、画像処理装置1を主に構成するコンピュータが記録媒体14のプログラムを実行することで、以下の処理が行われる。
【0024】
まず、ステレオ処理手段11は複数の画像データ(航空画像データa)に対してステレオマッチング処理を行い、DEMデータbを抽出する。DEM平滑化手段12はステレオ処理手段11によって得られたDEMデータbに平滑化フィルタを適用して平滑化DEMデータcを生成する。
【0025】
オルソ化手段13はDEM平滑化手段12によって平滑化された平滑化DEMデータcを用いて画像(航空画像データa)の正射変換を行い、オルソ画像dを生成する。ここで、オルソ画像とは、上述したように、通常のカメラの場合、図11(a)に示すように、中心投影のため傾いて撮影される対象を、図11(b)に示すように、所定の面に対して正射影して得られる画像のことであり、複数の画像を連結して広範囲の対象を把握し、図面等を作成するためには不可欠な画像である。
【0026】
図2は本発明の第1の実施例による画像処理装置1の動作を示すフローチャートである。これら図1及び図2を参照して本発明の第1の実施例による画像処理装置1の動作について説明する。尚、図2に示す処理は画像処理装置1を主に構成するコンピュータが記録媒体14のプログラムを実行することで実現される。
【0027】
まず、ステレオ処理手段11は複数の画像(航空画像データa)に対して自動的にステレオマッチング処理を行い、DEMデータを抽出する(図2ステップS1)。ここで用いるステレオマッチング処理には特に制限がなく、例えば、特開平3−167678号公報に開示された方法で実現することができる。
【0028】
次に、DEM平滑化手段12は、上記のステレオ処理手段11によって得られたDEMデータbに対して平滑化フィルタを適用し、平滑化DEMデータcを生成する(図2ステップS2)。ここで用いる平滑化フィルタはノイズを軽減する効果を持つフィルタであれば、特に制限はなく、例えば、メディアンフィルタ、多値モルフォロジカルフィルタ等を用いてもよい。また、それらのフィルタを複数組合せて適用することも可能である。
【0029】
オルソ化手段13は上記のDEM平滑化手段12によって平滑化された平滑化DEMデータcを用い、画像(航空画像データa)をある基準面に対して正射変換し、オルソ画像dを作成する(図2中のステップS3)。この時の正射投影面は、例えば、衛星画像や航空写真を用いる場合であれば、地上面に平行な面となるが、用途によって任意の面を定めて正射変換することができる。
【0030】
上述したように、本実施例では、平滑化したDEMデータをオルソ化処理に用いる構成とすることで、元々のDEMデータbに含まれるノイズの影響を軽減したオルソ画像dを生成することができる。
【0031】
図3は本発明の第2の実施例による画像処理装置の構成を示すブロック図である。図3において、本発明の第2の実施例はDEM平滑化手段12を省き、射影変換手段21と画像合成手段22と追加した以外は図1に示す本発明の第1の実施例と同様の構成となっており、同一構成要素には同一符号を付してある。また、同一構成要素の動作は本発明の第1の実施例と同様である。
【0032】
すなわち、画像処理装置2は複数の画像(航空画像データa)に対してステレオマッチング処理を行ってDEMデータbを抽出するステレオ処理手段11と、その抽出されたDEMデータbを用いて画像(航空画像データa)を正射変換するオルソ化手段13と、画像(航空画像データa)を正射投影面に対する回転角にしたがって射影変換する射影変換手段21と、射影変換画像(射影変換航空画像e)とオルソ画像dとを合成し、欠損した画素が補間された補間オルソ画像fを生成する画像合成手段22と、画像処理装置2の各手段を実現するためのプログラムを格納する記録媒体14とから構成されている。
【0033】
上記の画像処理装置2の各手段はそれぞれ概略、次のように動作する。つまり、画像処理装置2を主に構成するコンピュータが記録媒体14のプログラムを実行することで、以下の処理が行われる。
【0034】
ステレオ処理手段11は、上述したように、複数の画像データ(航空画像データa)に対してステレオマッチング処理を行い、DEMデータbを抽出する。オルソ化手段13は上記のステレオ処理手段11によって得られたDEMデータaを用いて画像(航空画像データa)を正射変換する。
【0035】
射影変換手段21は画像(航空画像データa)を標定計算によって得られた正射投影面への回転角にしたがって射影変換する。画像合成手段22は上記のオルソ化手段13によって得られたオルソ画像dと、上記の射影変換手段21によって変換された画像(射影変換航空画像e)とを重ね合わせて補間し、欠損した画素が補間された補間オルソ画像fを生成する。
【0036】
図4は本発明の第2の実施例による画像処理装置2の動作を示すフローチャートである。これら図3及び図4を参照して本発明の第2の実施例による画像処理装置2の動作について説明する。尚、図4に示す処理は画像処理装置2を主に構成するコンピュータが記録媒体14のプログラムを実行することで実現される。
【0037】
まず、ステレオ処理手段11は複数の画像(航空画像データa)に対して自動的にステレオマッチング処理を行い、DEMデータbを抽出する(図4ステップS11)。ここで用いるステレオマッチング処理には特に制限はなく、例えば、特開平3−167678号公報に開示された方法で実現することができる。
【0038】
次に、オルソ化手段13は上記のステレオ処理手段11によって得られたDEMデータbを用いて、画像(航空画像データa)をある基準面に対して正射変換し、オルソ画像dを作成する(図4ステップS12)。この時の正射投影面は、例えば、衛星画像や航空写真を用いる場合であれば、地上面に平行な面となるが、用途によって任意の面を定めて正射変換することができる。
【0039】
射影変換手段21は上記のステレオ処理手段11で用いた複数の画像(航空画像データa)のうちのいずれかの画像を、正射投影面への回転角にしたがって射影変換した画像(射影変換航空画像e)を生成する(図4ステップS13)。これは画像を平面で近似し、視点を固定したまま正射投影面に再投影する操作である。その際に用いる画像はステレオ処理手段11で用いた画像(航空画像データa)のうちの任意の画像を用いて良いが、正射投影面に対する回転角が小さい画像の方が望ましい。
【0040】
最後に、画像合成手段22はオルソ化手段13において得られたオルソ画像d及び射影変換変換手段21において得られた射影変換航空画像eを重ね合わせて補間し、補間オルソ画像fを生成する(図4ステップS14)。具体的には、オルソ画像d中に値が得られない抜け領域の画素に対し、対応する位置の射影変換航空画像eの画素を用いて補間する処理を行う。
【0041】
上述したように、本実施例では、画像(航空画像データa)を正射投影面に対する回転角に合せて射影変換し、その画像(射影変換航空画像e)の対応領域の画素をオルソ画像dの欠損領域において代用する構成としているため、オルソ画像dの抜け領域を補間することができ、その結果、オルソ画像dをCG等に応用する際、抜け領域によって発生していた違和感を減少させることができる。
【0042】
図5は本発明の第3の実施例による画像処理装置の構成を示すブロック図である。図5において、本発明の第3の実施例はDEM平滑化手段12を追加した以外は図3に示す本発明の第2の実施例と同様の構成となっており、同一構成要素には同一符号を付してある。また、同一構成要素の動作は本発明の第2の実施例と同様である。
【0043】
すなわち、画像処理装置3は複数の画像(航空画像データa)に対してステレオマッチング処理を行い、DEMデータbを抽出するステレオ処理手段11と、DEMデータbに平滑化フィルタを適用してノイズを除去するDEM平滑化手段12と、平滑化DEMデータcを用いて画像(航空画像データa)を正射変換するオルソ化手段13と、画像(航空画像データa)を正射投影面に対する標定角にしたがって射影変換する射影変換手段21と、射影変換航空画像eとオルソ画像dとを合成し、欠損した画素が補間された補間オルソ画像fを生成する画像合成手段22と、画像処理装置2の各手段を実現するためのプログラムを格納する記録媒体14とから構成されている。
【0044】
上記の画像処理装置3の各手段はそれぞれ概略、次のように動作する。つまり、画像処理装置3を主に構成するコンピュータが記録媒体14のプログラムを実行することで、以下の処理が行われる。
【0045】
ステレオ処理手段11は複数の画像データ(航空画像データa)に対してステレオマッチング処理を行い、DEMデータbを抽出する。DEM平滑化手段12は上記のステレオ処理手段11によって得られたDEMデータbに平滑化フィルタを適用し、平滑化DEMデータcを生成する。
【0046】
オルソ化手段13は上記のDEM平滑化手段12によって平滑化された平滑化DEMデータcを用いて画像(航空画像データa)を正射変換し、オルソ画像dを作成する。
【0047】
射影変換手段21は画像(航空画像データa)を、標定された正射投影面への回転角にしたがって射影変換する。画像合成手段22は上記のオルソ化手段13によって得られたオルソ画像dと、上記の射影変換手段21によって変換された画像(射影変換航空画像e)とを重ね合わせて補間し、欠損した画素が補間された補間オルソ画像fを生成する。
【0048】
図6は本発明の第3の実施例による画像処理装置3の動作を示すフローチャートである。これら図5及び図6を参照して本発明の第3の実施例による画像処理装置3の動作について説明する。尚、図6に示す処理は画像処理装置3を主に構成するコンピュータが記録媒体14のプログラムを実行することで実現される。
【0049】
まず、ステレオ処理手段11は複数の画像(航空画像データa)に対して自動的にステレオマッチング処理を行い、DEMデータbを抽出する(図6ステップS21)。ここで用いるステレオマッチング処理には、特に制限はなく、例えば、特開平3−167678号公報に開示された方法で実現することができる。
【0050】
次に、DEM平滑化手段12は上記のステレオ処理手段11によって得られたDEMデータbに対し、平滑化フィルタを適用して平滑化DEMデータcを作成する(図6ステップS22)。ここで用いる平滑化フィルタはノイズを軽減する効果を持つフィルタであれば、特に制限はなく、例えば、メディアンフィルタ、多値モルフォロジカルフィルタ等を用いてよい。また、それらのフィルタを複数組合せて適用することも可能である。
【0051】
オルソ化手段13は上記のDEM平滑化手段12によって平滑化された平滑化DEMデータcを用いて、画像(航空画像データa)をある基準面に対して正射変換し、オルソ画像dを作成する(図6ステップS23)。この時の正射投影面は、例えば、衛星画像や航空写真を用いる場合であれば、地上面に平行な面となるが、用途によって任意の面を定めて正射変換することができる。
【0052】
続いて、射影変換手段21は上記のステレオ処理手段11で用られた複数の画像(航空画像データa)のうちのいずれかの画像を、正射投影面への回転角にしたがって射影変換した画像(射影変換航空画像e)を生成する(図6ステップS24)。これは画像を平面で近似し、視点を固定したまま正射投影面に再投影する操作である。この際に用いる画像は、ステレオ処理手段11で用いた画像(航空画像データa)のうちの任意の画像を用いて良いが、正射投影面に対する回転角が小さい画像の方が望ましい。
【0053】
最後に、画像合成手段22は上記のオルソ化手段13において得られたオルソ画像d及び上記の射影変換変換手段21において得られた射影変換航空画像eを重ね合わせて補間し、欠損した画素が補間された補間オルソ画像fを生成する(図6ステップS25)。具体的には、オルソ画像d中に値が得られない抜け領域の画素に対し、対応する位置の射影変換航空画像eの画素を用いて補間する処理を行う。
【0054】
上述したように、本実施例では平滑化された平滑化DEMデータcをオルソ化処理に用いる構成としているので、元々のDEMデータbに含まれるノイズの影響を軽減したオルソ画像dを生成することができる。
【0055】
また、本実施例では、画像(航空画像データa)を正射投影面に対する回転角に合せて射影変換し、その画像(射影変換航空画像e)の対応領域の画素をオルソ画像dの欠損領域において代用する構成としているので、オルソ画像dの抜け領域を補間することができ、その結果、オルソ画像dをCG等に応用する際、抜け領域によって発生していた違和感を減少させることができる。
【0056】
図7は本発明の第4の実施例による画像処理装置の構成を示すブロック図である。図7において、本発明の第4の実施例はステレオペア設定手段42を追加した以外は図5に示す本発明の第3の実施例と同様の構成となっており、同一構成要素には同一符号を付してある。また、同一構成要素の動作は本発明の第3の実施例と同様である。
【0057】
すなわち、画像処理装置4は複数の画像(航空画像データa)の中から対になる画像を選択するステレオペア設定手段42と、オルソ画像dを生成するオルソ画像生成手段41と、画像(航空画像データa)を正射投影面に対する標定角にしたがって射影変換する射影変換手段21と、射影変換航空画像eとオルソ画像dとを合成し、欠損した画素が補間された補間オルソ画像fを生成する画像合成手段22と、画像処理装置2の各手段を実現するためのプログラムを格納する記録媒体14とから構成されている。
【0058】
ここで、オルソ画像生成手段41は複数の画像(航空画像データa)に対してステレオマッチング処理を行い、DEMデータbを抽出するステレオ処理手段11と、その抽出されたDEMデータbに平滑化フィルタを適用してノイズを除去するDEM平滑化手段12と、平滑化された平滑化DEMデータcを用いて画像(航空画像データa)を正射変換するオルソ化手段13とから構成されている。
【0059】
上記の画像処理装置4の各手段はそれぞれ概略、次のように動作する。つまり、画像処理装置4を主に構成するコンピュータが記録媒体14のプログラムを実行することで、以下の処理が行われる。
【0060】
ステレオペア設定手段42は同一地点を撮影した3枚以上の画像データ群から、うち2枚を選択して得られるステレオ画像対を複数決定し、処理単位とする。ステレオ処理手段11は上記のステレオペア設定手段42によって与えられる各ステレオ画像対に対してステレオマッチング処理を行い、DEMデータbを抽出する。
【0061】
DEM平滑化手段12は上記のステレオ処理手段11によって得られたDEMデータbに平滑化フィルタを適用して平滑化する。オルソ化手段13は上記のDEM平滑化手段12によって平滑化された平滑化DEMデータcを用いて画像(航空画像データa)を正射変換する。
【0062】
射影変換手段21は画像(航空画像データa)を、標定された正射投影面への回転角にしたがって射影変換する。画像合成手段22は上記のオルソ化手段13によって得られた各オルソ画像dと、上記の射影変換手段21によって変換された画像(射影変換航空画像e)とを重ね合わせて補間し、欠損した画素が補間された補間オルソ画像fを生成する。
【0063】
図8は本発明の第4の実施例による画像処理装置4の動作を示すフローチャートである。これら図7及び図8を参照して本発明の第4の実施例による画像処理装置4の動作について説明する。尚、図8に示す処理は画像処理装置4を主に構成するコンピュータが記録媒体14のプログラムを実行することで実現される。
【0064】
まず、ステレオペア設定手段42は同一地点を撮影した3枚以上の画像データ群から、うち2枚を選択して得られるステレオ画像対を複数決定し、処理単位とする(図8ステップS31)。例えば、4枚の画像を得た場合、これらのうちある2枚で一組、残りの2枚でもう一組のステレオ画像対を設定する。選べる組合せは全部で6通り存在するが、必ずしもこれらを全てステレオ画像対として設定する必要はなく、隠蔽領域を減らすために全ての画像が使われるよう、ステレオ画像対を複数定めればよい。
【0065】
次に、ステレオ処理手段11は上記のステレオペア設定手段42によって設定された各ステレオ画像対についてステレオマッチング処理を行い、DEMデータbを抽出する(図8ステップS32)。ここで用いるステレオマッチング処理には、特に制限はなく、例えば、特開平3−167678号公報に開示された方法で実現することができる。
【0066】
続いて、DEM平滑化手段12は上記のステレオ処理手段11によって得られた各DEMデータbに対し、平滑化フィルタを適用して平滑化処理を行う(図8ステップS33)。ここで用いる平滑化フィルタはノイズを軽減する効果を持つフィルタであれば、特に制限はなく、例えば、メディアンフィルタ、多値モルフォロジカルフィルタ等を用いてよい。また、それらのフィルタを複数組合せて適用することも可能である。
【0067】
オルソ化手段13は上記のDEM平滑化手段12によって平滑化された各DEMデータ(平滑化DEMデータc)を用いて、画像(航空画像データa)をある基準面に対して正射変換し、オルソ画像dを作成する(図8ステップS34)。この時の正射投影面は、例えば、衛星画像や航空写真を用いる場合であれば、地上面に平行な面となるが、用途によって任意の面を定めて正射変換することができる。
【0068】
射影変換手段21は上記のステレオ処理手段11で用いた複数の画像(航空画像データa)のうちのいずれかの画像を正射投影面への回転角にしたがって射影変換した画像(射影変換航空画像e)を生成する(図8ステップS35)。これは画像(航空画像データa)を平面で近似し、視点を固定したまま正射投影面に再投影する操作である。その際に用いる画像は、ステレオ処理手段11で用いた画像のうちの任意の画像を用いて良いが、正射投影面に対する回転角が小さい画像の方が望ましい。
【0069】
最後に、画像合成手段22は上記のオルソ化手段13において得られた各オルソ画像d及び上記の射影変換変換手段21において得られた射影変換航空画像eを重ね合わせて補間し、欠損した画素が補間された補間オルソ画像fを生成する(図8ステップS36)。具体的には、オルソ画像d中に値が得られない抜け領域の画素に対し、他のオルソ画像中の対応画素を用いて補間する。この処理で残った抜け領域の画素に対しては、射影変換画像中の対応する画素を用いて補間する処理を行う。
【0070】
尚、上述した本発明の第4の実施例では、複数のステレオペアに対する処理をステレオ処理手段11、DEM平滑化手段12、オルソ化手段13において繰り返し行っているが、例えば上記のオルソ画像生成手段41を複数設け、複数のステレオペアについての処理を並列に行わせることも可能である。
【0071】
上述したように、本実施例では、平滑化した平滑化DEMデータcをオルソ化処理に用いているため、元々のDEMデータに含まれるノイズの影響を軽減したオルソ画像を生成することができる。
【0072】
また、本実施例では、視線方向の異なる複数のステレオペアからオルソ画像を相互補完して生成し、加えて画像を正射投影面に対する回転角に合せて射影変換し、その画像の対応領域の画素をオルソ画像の欠損領域において代用しているので、オルソ画像の抜け領域を補間することができ、その結果、オルソ画像をCG等に応用する際、抜け領域によって発生していた違和感を減少させることが可能となる。
【0073】
尚、本発明は、上述した本発明の第1〜第4の実施例各々に限定されるものではなく、様々な条件での実施が可能である。例えば、対象となる画像は航空写真や衛星画像に限定されるものではなく、地上あるいは室内においてステレオ撮影を行った画像に対しても適用が可能である。また画像は2台のカメラを使用して撮影しても、単眼カメラを2度用いて撮影したものでもよく、画像データについても、ディジタルカメラで直接撮影したもの、あるいはアナログ写真をスキャンしたデータであっても実施可能である。
【0074】
続いて、具体例を用いて本発明の実施例の動作について説明する。この説明ではステレオ航空写真から地上面に平行なオルソ画像を、上記の本発明の第3の実施例にしたがって生成する場合について述べる。
【0075】
ステレオ航空写真は連続して航空機から撮影された画像の重複部分をステレオ画像と見なして扱われ、まずスキャナを用いてディジタルデータ化される。続いて、撮影時のカメラパラメータを決定するため、標定作業が行われる。通常、鉛直に近い方向で撮影した航空写真でも、航空機の姿勢変化等の影響によって、撮影時にカメラパラメータを把握するのは難しいため、通常左右の画像中で共通に映っている対象の、それぞれの画像中における座標値を読取り、この対応点の座標対を用いて左右画像間のカメラ回転角等のパラメータを標定する。ここで用いる標定方法には、特に制限はなく、例えば、図化ソフトによって算出される値を使用しても良い。
【0076】
カメラパラメータの標定が行われた後、図9に示すように、左右の画像を必要に応じて共通平行面に再投影する。これは平行化と呼ばれる操作であり、平行化を行った後の左右画像中の同一対象は同一走査線上にくることになる。
【0077】
また、平行化を行った左右の画像に対し、ステレオ処理手段11において、自動的にステレオマッチング処理を行い、左右の画像間の対応を求めて密な視差画像を生成する。ここで用いるステレオマッチング処理には、特に制限はなく、例えば、特開平3−167678号公報に開示された方法で実現することができる。
【0078】
視差画像が得られたならば、標定されたカメラパラメータと組合せて各画素に対応する3次元座標系での位置を計算し、DEMデータを抽出する。得られたDEMデータはステレオマッチング処理時に対応付けを誤った画素等の影響によって、結果として不連続な3次元座標を持つスパイク状の孤立点をノイズとして多く含むこととなる。
【0079】
これら孤立点の影響を除くため、DEM平滑化手段12において、上記のステレオ処理手段11によって得られたDEMデータに対し、平滑化フィルタを適用して平滑化処理を行う。ここで用いる平滑化フィルタはノイズを軽減する効果を持つフィルタであれば、特に制限はないが、例えば、対象画素を中心とする小領域における画素値のメディアンを代表値として置換するメディアンフィルタを用いれば、実際の形状に依存するエッジを維持しつつ、上記のような孤立した雑音を効率的に除去することができる。
【0080】
その他、小領域に対する最大値フィルタと最小値フィルタとを組合せて構成される多値モルフォロジカルフィルタ等も、上記と同様の効果が期待できる。このようなフィルタを複数組合せて適用することも可能である。
【0081】
続いて、上記のDEM平滑化手段12によって平滑化され、ノイズが除去された平滑化DEMデータ及び標定されているカメラパラメータを用いて、オルソ化手段13によって元となる航空画像の地上面に対する正射影変換を行い、オルソ画像を作成する。
【0082】
具体的には、まずDEMデータ中の全ての点における3次元座標を地上面に投影した場合の地上座標系(例えば、平面直角座標系)での平面座標と、元の画像中の座標との対応を求める。次に、求めた対応にしたがって、元となる航空画像の画素を地表面に対応する画像に補間しながら写像する。
【0083】
さらに、射影変換手段21において、航空画像データを地表面に平行な面に射影変換によって再投影する作業が行われる。上述したように、航空写真撮影時のカメラの光軸は地上面に対して完全に鉛直ではなく、地表を斜めから撮影しているため、オルソ画像との重ね合わせは難しい。このため、航空画像に射影変換を適用して、カメラ中心位置を変えず、地上に平行な投影面に再投影する操作を適用する。この操作によって完全な鉛直写真に相当する画像を生成することができる。
【0084】
射影変換の具体的な方法は、撮影時のカメラパラメータのうち、地上面に対するカメラ座標系の回転角と、焦点距離とを用いて、例えば、特開平8−567号公報に開示された方法によって行うことができる。
【0085】
最後に、画像合成手段22は上記のオルソ化手段13によって生成されたオルソ画像と、上記の射影変換手段21によって生成された射影変換画像とを合成し、補間オルソ画像を生成する処理を行う。具体的には、オルソ画像と射影変換画像とを重ね合わせた後、隠蔽の影響によって値が得られていないオルソ画像中の抜け領域の画素について、対応する位置の射影変換画像の画素を引用し、補間する処理を行うことで補間オルソ画像を生成する。
【0086】
この際に画素を引用する元となる画像は、左右の射影変換画像のうち、対象となる位置に主点が近い位置に存在する画像、すなわち対象点になるべく正対する方を用いる。こうすることによって、平面近似のため倒れこんで写る建物壁面が隠蔽領域に使用されてしまう場合を減少させ、画質を向上させることができる。
【0087】
また、引用の際に境界部あるいは引用部をぼかすことによって、発生するエッジの影響を抑え、より違和感の少ない補間オルソ画像を生成することもできる。
【0088】
次に、上記第4の実施の形態の実施例に従って、4枚のステレオ航空写真から地上面に平行なオルソ画像を生成する場合の具体例を説明する。
【0089】
4枚のステレオ航空写真は、例えば、図10に示すように、2つの隣接する撮影コース#1,#2間の重なりを大きくとって撮影を行い、同一地点を撮影している画像対をそれぞれのコースから1組ずつ選択することによって得ることができる。しかしながら、撮影方法はこれに限られるものではなく、個別に撮影して得られた4枚の画像#1〜#4を用いても良い。得られた航空写真群は、まずスキャナを用いてディジタルデータ化される。
【0090】
次に、ステレオペア設定手段42において、ディジタルデータ化された4枚のステレオ航空写真からステレオペアを設定する。画像に付加されるメタデータ等を参照し、同じ撮影コース内の画像2枚をそれぞれ組合せて2組のステレオペアを得る。しかしながら、ペアの設定方法はこれに限るものではなく、4枚の画像#1〜#4全てを用いるならば、異なるコース間で画像対を設定しても良い。また、各画像#1〜#4を複数回使用して、最大6組のステレオペアを生成しても構わない。
【0091】
設定された各ステレオペアについて、本発明の第3の実施例と同様に、撮影時のカメラパラメータを決定する標定作業、平行化、ステレオマッチング処理、DEMデータの抽出、DEMデータの平滑化が行われる。
【0092】
続いて、上記のDEM平滑化手段12によって平滑化し、ノイズを除去された平滑化DEMデータ及び標定されているカメラパラメータを用いて、オルソ化手段13によって各航空写真毎に元となる航空画像の地上面に対する正射影変換を行い、オルソ画像を作成する。
【0093】
具体的には、まずDEMデータ中の全ての点における3次元座標を地上面に投影した場合の、地上座標系(例えば、平面直角座標系)での平面座標と、元の画像中の座標の対応を求める。次に、求めた対応にしたがって、元となる航空画像の画素を地表面に対応する画像に補間しながら写像する。
【0094】
さらに、射影変換手段21において、各ステレオペアの航空画像データを地表面に平行な面に射影変換によって再投影する作業が行われる。上述したように、航空写真撮影時のカメラの光軸は地上面に対して完全に鉛直ではなく、地表を斜めから撮影しているため、オルソ画像との重ね合わせは難しい。このため、航空画像に射影変換を適用し、カメラ中心位置を変えず、地上に平行な投影面に再投影する操作を適用する。この操作によって完全な鉛直写真に相当する画像を生成することができる。
【0095】
射影変換の具体的な方法は、撮影時のカメラパラメータのうちの地上面に対するカメラ座標系の回転角と焦点距離とを用いて、例えば特開平8−567号公報に記載されているような方法によってこれを行える。
【0096】
最後に、画像合成手段22において、上記のオルソ化手段13によって生成された各オルソ画像と、上記射影変換手段21によって生成された各射影変換画像とを合成し、補間オルソ画像を生成する処理を行う。
【0097】
具体的には、各オルソ画像と各射影変換画像を重ね合わせた後、まず複数のオルソ画像を画素を相互補間する形で合成する。すなわち、あるオルソ画像についてのみ画素が得られている点についてはその値を用い、複数のオルソ画像で値が得られている場合には、そのうちの対象となる位置に主点が近い位置に存在する画像、つまりその画素をなるべく正対して撮影している画像中の画素値を用いて合成する。
【0098】
さらに、隠蔽の影響によって値が得られていないオルソ画像中の抜け領域の画素については、対応する位置の射影変換画像の画素を引用し、補間する処理を行うことで補間オルソ画像を生成する。この際に画素を引用する元となる画像は、複数の射影変換画像のうち、対象となる位置に主点が近い位置に存在する画像、すなわち対象点になるべく正対する方を用いる。これによって、平面近似のために倒れこんで写る建物壁面が隠蔽領域に使用されてしまう場合を減少させ、画質を向上させることができる。また、引用の際に境界部あるいは引用部をぼかすことによって、発生するエッジの影響を抑え、より違和感の少ない補間オルソ画像を生成することもできる。
【0099】
このように、本発明は、平滑化したDEMデータをオルソ化処理に用いることによって、元のDEMデータに含まれるノイズを軽減することができ、ノイズのためにオルソ画像の画質が低下するのを防ぐことができるので、ノイズの少ない自然なオルソ画像を生成することができる。
【0100】
また、本発明は、画像を正射投影面に対する回転角に合せて射影変換し、その画像の対応領域の画素をオルソ画像の欠損領域において代用することによって、オルソ画像の抜け領域を補間することができるので、抜け領域による違和感の少ないオルソ画像を生成することができる。
【0101】
【発明の効果】
以上説明したように本発明の画像処理装置は、複数の画像に対してステレオマッチング処理を行って3次元データを生成する画像処理装置において、ステレオマッチング処理によって3次元データを抽出し、その抽出された3次元データに平滑化フィルタを適用して3次元データのノイズを除去して平滑化し、この平滑化された3次元データを用いて画像を正射変換することによって、3次元データに含まれるノイズの影響を受けにくく、また欠損領域に起因する違和感の少ない、より自然なオルソ画像を得ることができるという効果が得られる。
【0102】
また、本発明の他の画像処理装置は、複数の画像に対してステレオマッチング処理を行って3次元データを生成する画像処理装置において、ステレオマッチング処理によって3次元データを抽出し、この3次元データを用いて画像を正射変換し、画像を正射投影面に対する標定角にしたがって射影変換し、画像の正射変換で得られたオルソ画像と射影変換で得られた射影変換画像とを合成して欠損した画素が補間された補間オルソ画像を生成することによって、3次元データに含まれるノイズの影響を受けにくく、また欠損領域に起因する違和感の少ない、より自然なオルソ画像を得ることができるという効果が得られる。
【図面の簡単な説明】
【図1】本発明の第1の実施例による画像処理装置の構成を示すブロック図である。
【図2】本発明の第1の実施例による画像処理装置の動作を示すフローチャートである。
【図3】本発明の第2の実施例による画像処理装置の構成を示すブロック図である。
【図4】本発明の第2の実施例による画像処理装置の動作を示すフローチャートである。
【図5】本発明の第3の実施例による画像処理装置の構成を示すブロック図である。
【図6】本発明の第3の実施例による画像処理装置の動作を示すフローチャートである。
【図7】本発明の第4の実施例による画像処理装置の構成を示すブロック図である。
【図8】本発明の第4の実施例による画像処理装置の動作を示すフローチャートである。
【図9】本発明のステレオ撮影幾何を示す模式図である。
【図10】本発明の実施例における4方向からの航空写真撮影例を示す模式図である。
【図11】(a),(b)は中心投影画像と正射画像との違いを示す模式図である。
【図12】従来例による画像処理装置の構成を示すブロック図である。
【図13】従来例による画像処理装置の動作を示すフローチャートである。
【図14】オルソ画像における抜け領域の発生を示す模式図である。
【符号の説明】
1,2,3,4 画像処理装置
11 ステレオ処理手段
12 DEM平滑化手段
13 オルソ化手段
14 記録媒体
21 射影変換手段
22 画像合成手段
41 オルソ画像生成手段
42 ステレオペア設定手段
[0001]
BACKGROUND OF THE INVENTION
The present invention relates to an image processing apparatus, an image processing method used therefor, and a program therefor, and more particularly to an image processing method used for three-dimensional data generation from stereo images and orthographic transformation.
[0002]
[Prior art]
Conventionally, in digital photogrammetry and cartographic work, three-dimensional data [hereinafter referred to as DEM (Digital Evelopment Map) data, based on stereo images obtained from cameras installed on the ground or images taken from satellites or aircraft, Is generated by stereo matching, and the orthorectified image is generated by using the stereo matching to generate an ortho image and applied to drawing or the like.
[0003]
Here, in the case of a normal camera, an ortho image is an object that is photographed with an inclination for central projection, as shown in FIG. 11 (a), as shown in FIG. 11 (b). Is an image obtained by orthogonally projecting a plurality of images, grasping a wide range of objects, and creating a drawing or the like.
[0004]
As an example of applying the above method, there is a technique using the above method when performing digital image measurement on the ground (for example, see Patent Documents 1 and 2). In addition, as an example of a system in which the above technique is applied to aerial photographs and satellite images, there is a technique for stereo-matching aerial photographs and satellite images to perform orthogonal conversion (see, for example, Patent Document 3).
[0005]
FIG. 12 is a block diagram illustrating a configuration example of a conventional image processing apparatus. As shown in FIG. 12, a conventional image processing apparatus 5 performs stereo matching processing on a plurality of images to extract DEM data, and an orthorectifier that orthorectifies images using the DEM data. And a conversion means 13.
[0006]
In the image processing apparatus 5 having the above configuration, an operation as shown in FIG. 13 is performed. That is, image data obtained from a plurality of photographs is given to the stereo processing means 11, and DEM data is automatically extracted by performing stereo matching as in the conventional system (step S41 in FIG. 13).
[0007]
The DEM data obtained in this way is given to the orthorizing means 13 together with the image data, and an ortho image is generated by moving the pixels in the image based on the DEM information and the standardized shooting angle information (FIG. 13). Step S42).
[0008]
[Patent Document 1]
JP 2000-182059 A (3rd to 7th pages, FIGS. 1 and 2)
[Patent Document 2]
Japanese Patent Laid-Open No. 11-351865 (pages 4, 5 and 1)
[Patent Document 3]
Japanese Unexamined Patent Publication No. 2000-348160 (pages 2, 3 and 1)
[0009]
[Problems to be solved by the invention]
However, in the conventional ortho image generation method described above, an error caused by mismatching of corresponding points is included in the DEM data obtained by stereo matching, and pixels having erroneous DEM data at the time of orthonormal transformation are included. Is moved to an incorrect position, and there is a problem that the image quality of the ortho image generated by the above method is low.
[0010]
In addition, in the conventional ortho image generation method, when an imaging target in an image includes many discontinuous shapes such as an aerial image capturing a city area, pixels are not captured in a concealed area such as a road hidden by a building or the like. Therefore, the corresponding pixel cannot be moved in the ortho image, and as shown in FIG. 14, a missing region in which the pixel is lost is generated in the concealed region, so that the pixel is lost in the ortho image generated by the above method. There is a problem in that many regions are included, resulting in an unnatural image.
[0011]
Since the conventional ortho image generation method has the above-described problems, there is a demand for a processing method that can be used for drawing and CG (Computer Graphics) generation and that can obtain a more natural ortho image.
[0012]
Accordingly, an object of the present invention is to solve the above-described problems, and to be able to obtain a more natural ortho image that is less affected by noise included in DEM data and less uncomfortable due to a defective area. And an image processing method used therefor, a program image processing apparatus thereof, an ortho image automatic shaping method used therefor, and a program thereof.
[0013]
[Means for Solving the Problems]
An image processing apparatus according to the present invention is an image processing apparatus that generates a three-dimensional data by performing a stereo matching process on a plurality of images, and a stereo processing unit that extracts the three-dimensional data by the stereo matching process; A DEM smoothing unit that applies a smoothing filter to the three-dimensional data extracted by the stereo processing unit to remove noise from the three-dimensional data and smoothes the three-dimensional data, and is smoothed by the DEM smoothing unit. The image using dimensional data All pixels included in The Individually And orthorectifying means for orthogonal transformation.
[0014]
Another image processing apparatus according to the present invention is an image processing apparatus that generates a three-dimensional data by performing a stereo matching process on a plurality of images, and a stereo processing unit that extracts the three-dimensional data by the stereo matching process An orthorectifying means for orthorectifying the image using the three-dimensional data, a projective transforming means for projectively transforming the image in accordance with an orientation angle with respect to an orthographic projection plane, and an orthorectification obtained by the orthorectifying means. Image synthesizing means for synthesizing the image and the projective transformation image obtained by the projective transformation means to generate an interpolated ortho image in which the missing pixels are interpolated.
[0015]
An image processing method according to the present invention is an image processing method for generating three-dimensional data by performing stereo matching processing on a plurality of images in an image processing device. A step of extracting dimensional data, a step of applying a smoothing filter to the extracted three-dimensional data to remove noise of the three-dimensional data and smoothing, and using the smoothed three-dimensional data The image All pixels included in The Individually And an orthophoto conversion step.
[0016]
Another image processing method according to the present invention is an image processing method for generating three-dimensional data by performing stereo matching processing on a plurality of images, and extracting the three-dimensional data by the stereo matching processing; Orthorectifying the image using the three-dimensional data, projecting transforming the image in accordance with an orientation angle with respect to an orthographic projection plane, an ortho image obtained by orthographic transformation of the image, and the projection And a step of generating an interpolated ortho image in which the missing pixels are interpolated by synthesizing the projected conversion image obtained by the conversion.
[0017]
A program for an image processing method according to the present invention is a program for an image processing method for generating three-dimensional data by performing stereo matching processing on a plurality of images. The computer program stores the three-dimensional data by the stereo matching processing. A process of extracting, a process of applying a smoothing filter to the extracted three-dimensional data to remove noise of the three-dimensional data and smoothing, and the image using the smoothed three-dimensional data All pixels included in The Individually And orthogonal transformation processing is executed.
[0018]
Another image processing method program according to the present invention is a program of an image processing method for generating three-dimensional data by performing stereo matching processing on a plurality of images. Obtained by processing for extracting data, processing for orthorectifying the image using the three-dimensional data, processing for projective transformation of the image according to the orientation angle with respect to the orthographic projection plane, and orthographic transformation of the image The generated ortho image and the projective transformation image obtained by the projective transformation are combined to generate an interpolated ortho image in which the missing pixels are interpolated.
[0019]
That is, the image processing apparatus of the present invention performs stereo matching processing on a plurality of images to extract three-dimensional data [hereinafter referred to as DEM (Digital Elevation Map) data], and smoothes the DEM data. A DEM smoothing means for removing noise by applying a normalization filter, an orthorectifying means for orthogonally transforming an image using DEM data, and a projective transforming means for projectively transforming an image according to an orientation angle with respect to an orthographic projection plane And an image synthesizing means for synthesizing the projected transformation image and the ortho image to generate an interpolated ortho image.
[0020]
In the image processing apparatus according to the present invention, the stereo matching process is automatically performed on a plurality of images with the above-described configuration, and height information at each point in the image based on the obtained parallax information, that is, DEM. After extracting the data and applying a smoothing filter to remove noise in the DEM data, orthorectification is generated by using the DEM data obtained by smoothing the image, and the orthographic projection plane is generated. Is combined with a projective transformation image obtained by projective transformation with respect to a plane parallel to the axis, and an interpolated ortho image is generated by interpolating the missing pixels.
[0021]
Thus, in the image processing apparatus of the present invention, noise in the DEM data is reduced and used for orthorectification, and the pixel defect area due to concealment is interpolated with the projective transformation image, thereby reducing the image quality due to noise and the defect area. It is possible to generate a more natural ortho image with less sense of incongruity.
[0022]
DETAILED DESCRIPTION OF THE INVENTION
Next, embodiments of the present invention will be described with reference to the drawings. FIG. 1 is a block diagram showing the configuration of an image processing apparatus according to the first embodiment of the present invention. In FIG. 1, an image processing apparatus 1 performs stereo matching processing on a plurality of images (aerial image data a) to extract three-dimensional data [hereinafter referred to as DEM (Digital Elevation Map) data] b. A DEM smoothing means 12 for removing noise by applying a smoothing filter (not shown) to the extracted DEM data b, and an image (aerial image data a) using the smoothed DEM data c. An orthorectifying means 13 for conversion and a recording medium 14 for storing a program for realizing each means of the image processing apparatus 1 are configured.
[0023]
Each means of the image processing apparatus 1 generally operates as follows. That is, the following processing is performed by the computer that mainly configures the image processing apparatus 1 executing the program of the recording medium 14.
[0024]
First, the stereo processing means 11 performs stereo matching processing on a plurality of image data (aerial image data a), and extracts DEM data b. The DEM smoothing means 12 generates a smoothed DEM data c by applying a smoothing filter to the DEM data b obtained by the stereo processing means 11.
[0025]
The orthorectifying means 13 performs an orthogonal transformation of the image (aerial image data a) using the smoothed DEM data c smoothed by the DEM smoothing means 12 to generate an ortho image d. Here, as described above, in the case of an ordinary camera, an ortho image is an object that is photographed with an inclination for central projection, as shown in FIG. 11 (a), as shown in FIG. 11 (b). This is an image obtained by orthogonal projection with respect to a predetermined plane, and is an indispensable image for connecting a plurality of images to grasp a wide range of objects and creating a drawing or the like.
[0026]
FIG. 2 is a flowchart showing the operation of the image processing apparatus 1 according to the first embodiment of the present invention. The operation of the image processing apparatus 1 according to the first embodiment of the present invention will be described with reference to FIGS. Note that the processing shown in FIG. 2 is realized by a computer that mainly configures the image processing apparatus 1 executing the program of the recording medium 14.
[0027]
First, the stereo processing means 11 automatically performs a stereo matching process on a plurality of images (aerial image data a) to extract DEM data (step S1 in FIG. 2). There is no restriction | limiting in particular in the stereo matching process used here, For example, it can implement | achieve by the method disclosed by Unexamined-Japanese-Patent No. 3-167678.
[0028]
Next, the DEM smoothing means 12 applies a smoothing filter to the DEM data b obtained by the stereo processing means 11 to generate smoothed DEM data c (step S2 in FIG. 2). The smoothing filter used here is not particularly limited as long as it has an effect of reducing noise. For example, a median filter, a multi-value morphological filter, or the like may be used. It is also possible to apply a combination of these filters.
[0029]
The orthorectifying means 13 uses the smoothed DEM data c smoothed by the DEM smoothing means 12 and orthorectifies the image (aerial image data a) with respect to a certain reference plane to create an ortho image d. (Step S3 in FIG. 2). The orthographic projection plane at this time is, for example, a plane parallel to the ground surface in the case of using a satellite image or an aerial photograph.
[0030]
As described above, in the present embodiment, by using the smoothed DEM data for the orthorectification process, it is possible to generate an ortho image d in which the influence of noise included in the original DEM data b is reduced. .
[0031]
FIG. 3 is a block diagram showing the configuration of the image processing apparatus according to the second embodiment of the present invention. In FIG. 3, the second embodiment of the present invention is the same as the first embodiment of the present invention shown in FIG. 1 except that the DEM smoothing means 12 is omitted and a projective transformation means 21 and an image composition means 22 are added. The same components are denoted by the same reference numerals. The operation of the same component is the same as that of the first embodiment of the present invention.
[0032]
That is, the image processing apparatus 2 performs stereo matching processing on a plurality of images (aerial image data a) to extract DEM data b, and uses the extracted DEM data b to generate an image (aerial image). Orthorectification means 13 for orthographic transformation of the image data a), projective transformation means 21 for projective transformation of the image (aerial image data a) according to the rotation angle with respect to the orthographic projection plane, and a projective transformation image (projection transformation aerial image e ) And the ortho image d to generate an interpolated ortho image f in which the missing pixels are interpolated, and a recording medium 14 for storing a program for realizing each means of the image processing apparatus 2 It is composed of
[0033]
Each means of the above-described image processing apparatus 2 generally operates as follows. That is, the following processing is performed by the computer that mainly configures the image processing apparatus 2 executing the program of the recording medium 14.
[0034]
As described above, the stereo processing means 11 performs stereo matching processing on a plurality of image data (aerial image data a), and extracts DEM data b. The orthorectifying means 13 orthorectifies the image (aerial image data a) using the DEM data a obtained by the stereo processing means 11 described above.
[0035]
The projective conversion means 21 performs projective conversion on the image (aerial image data a) according to the rotation angle to the orthographic projection plane obtained by the orientation calculation. The image synthesizing means 22 interpolates the ortho image d obtained by the orthorectifying means 13 and the image (projective transformation aerial image e) transformed by the projective transform means 21 so that the missing pixel is replaced. An interpolated ortho ortho image f is generated.
[0036]
FIG. 4 is a flowchart showing the operation of the image processing apparatus 2 according to the second embodiment of the present invention. The operation of the image processing apparatus 2 according to the second embodiment of the present invention will be described with reference to FIGS. Note that the processing shown in FIG. 4 is realized by a computer that mainly configures the image processing apparatus 2 executing the program of the recording medium 14.
[0037]
First, the stereo processing means 11 automatically performs stereo matching processing on a plurality of images (aerial image data a) to extract DEM data b (step S11 in FIG. 4). There is no restriction | limiting in particular in the stereo matching process used here, For example, it can implement | achieve by the method disclosed by Unexamined-Japanese-Patent No. 3-167678.
[0038]
Next, using the DEM data b obtained by the stereo processing means 11, the orthorectifying means 13 orthorectifies the image (aerial image data a) with respect to a certain reference plane to create an ortho image d. (FIG. 4, step S12). The orthographic projection plane at this time is, for example, a plane parallel to the ground surface in the case of using a satellite image or an aerial photograph.
[0039]
Projection conversion means 21 is an image (projection conversion aviation) obtained by projective conversion of any one of the plurality of images (aerial image data a) used in the stereo processing means 11 according to the rotation angle to the orthographic projection plane. An image e) is generated (step S13 in FIG. 4). This is an operation of approximating an image on a plane and reprojecting it on an orthographic projection plane with the viewpoint fixed. The image used at that time may be any image of the images (aerial image data a) used in the stereo processing means 11, but an image with a small rotation angle with respect to the orthographic projection plane is desirable.
[0040]
Finally, the image synthesizing means 22 superimposes the orthoimage d obtained by the orthorectifying means 13 and the projective transformation aerial image e obtained by the projective transformation conversion means 21 to interpolate to generate an interpolated orthoimage f (FIG. 4 step S14). Specifically, a process of interpolating using pixels of the projective transformation aerial image e at the corresponding position is performed on the pixels of the missing region where no value is obtained in the ortho image d.
[0041]
As described above, in this embodiment, the image (aerial image data a) is projectively converted in accordance with the rotation angle with respect to the orthographic projection plane, and the pixels in the corresponding region of the image (projection-converted aerial image e) are orthoimage d. Therefore, the missing area of the ortho image d can be interpolated. As a result, when the ortho image d is applied to CG or the like, the uncomfortable feeling caused by the missing area is reduced. Can do.
[0042]
FIG. 5 is a block diagram showing the configuration of the image processing apparatus according to the third embodiment of the present invention. In FIG. 5, the third embodiment of the present invention has the same configuration as that of the second embodiment of the present invention shown in FIG. 3 except that the DEM smoothing means 12 is added, and the same components are the same. The code | symbol is attached | subjected. The operation of the same component is the same as that of the second embodiment of the present invention.
[0043]
That is, the image processing device 3 performs stereo matching processing on a plurality of images (aerial image data a), and extracts stereo noise by applying a smoothing filter to the DEM data b and stereo processing means 11 that extracts DEM data b. DEM smoothing means 12 to be removed, orthorectifying means 13 for orthorectifying the image (aerial image data a) using the smoothed DEM data c, and an orientation angle for the image (aerial image data a) with respect to the orthographic projection plane The projection conversion means 21 for performing projective conversion according to the above, the image synthesis means 22 for synthesizing the projection-transformed aerial image e and the ortho image d, and generating an interpolated ortho image f in which the missing pixels are interpolated; The recording medium 14 stores a program for realizing each means.
[0044]
Each means of the image processing apparatus 3 operates roughly as follows. That is, the following processing is performed by the computer that mainly configures the image processing apparatus 3 executing the program of the recording medium 14.
[0045]
The stereo processing means 11 performs stereo matching processing on a plurality of image data (aerial image data a), and extracts DEM data b. The DEM smoothing means 12 applies a smoothing filter to the DEM data b obtained by the stereo processing means 11 to generate smoothed DEM data c.
[0046]
The orthorectifying means 13 orthorectifies the image (aerial image data a) using the smoothed DEM data c smoothed by the DEM smoothing means 12 to create an ortho image d.
[0047]
The projective conversion means 21 performs projective conversion of the image (aerial image data a) according to the rotation angle to the orthographic projection plane that has been determined. The image synthesizing means 22 interpolates the ortho image d obtained by the orthorectifying means 13 and the image (projective transformation aerial image e) transformed by the projective transform means 21 so that the missing pixel is replaced. An interpolated ortho ortho image f is generated.
[0048]
FIG. 6 is a flowchart showing the operation of the image processing apparatus 3 according to the third embodiment of the present invention. The operation of the image processing apparatus 3 according to the third embodiment of the present invention will be described with reference to FIGS. Note that the processing shown in FIG. 6 is realized by a computer that mainly configures the image processing apparatus 3 executing the program of the recording medium 14.
[0049]
First, the stereo processing means 11 automatically performs stereo matching processing on a plurality of images (aerial image data a) to extract DEM data b (step S21 in FIG. 6). There is no restriction | limiting in particular in the stereo matching process used here, For example, it can implement | achieve by the method disclosed by Unexamined-Japanese-Patent No. 3-167678.
[0050]
Next, the DEM smoothing means 12 applies a smoothing filter to the DEM data b obtained by the stereo processing means 11 to create smoothed DEM data c (step S22 in FIG. 6). The smoothing filter used here is not particularly limited as long as it has an effect of reducing noise. For example, a median filter, a multi-value morphological filter, or the like may be used. It is also possible to apply a combination of these filters.
[0051]
The orthorectifying means 13 orthorectifies the image (aerial image data a) with respect to a certain reference plane using the smoothed DEM data c smoothed by the DEM smoothing means 12 to create an ortho image d. (Step S23 in FIG. 6). The orthographic projection plane at this time is, for example, a plane parallel to the ground surface in the case of using a satellite image or an aerial photograph.
[0052]
Subsequently, the projective transformation means 21 is an image obtained by projective transformation of any one of the plurality of images (aerial image data a) used in the stereo processing means 11 according to the rotation angle to the orthographic projection plane. (Projection conversion aerial image e) is generated (step S24 in FIG. 6). This is an operation of approximating an image on a plane and reprojecting it on an orthographic projection plane with the viewpoint fixed. As an image used at this time, an arbitrary image of the images (aerial image data a) used in the stereo processing unit 11 may be used, but an image having a small rotation angle with respect to the orthographic projection plane is desirable.
[0053]
Finally, the image synthesizing means 22 interpolates the orthoimage d obtained by the orthorectifying means 13 and the projected transformation aerial image e obtained by the projective transformation converting means 21 so as to interpolate the missing pixels. The interpolated ortho image f is generated (step S25 in FIG. 6). Specifically, a process of interpolating using pixels of the projective transformation aerial image e at the corresponding position is performed on the pixels of the missing region where no value is obtained in the ortho image d.
[0054]
As described above, since the smoothed smoothed DEM data c is used in the orthorectification process in the present embodiment, the ortho image d in which the influence of noise included in the original DEM data b is reduced is generated. Can do.
[0055]
In this embodiment, the image (aerial image data a) is projectively converted in accordance with the rotation angle with respect to the orthographic projection plane, and pixels in the corresponding region of the image (projection-transformed aerial image e) are deleted regions of the ortho image d. Therefore, the missing area of the ortho image d can be interpolated. As a result, when the ortho image d is applied to CG or the like, the uncomfortable feeling caused by the missing area can be reduced.
[0056]
FIG. 7 is a block diagram showing a configuration of an image processing apparatus according to the fourth embodiment of the present invention. 7, the fourth embodiment of the present invention has the same configuration as that of the third embodiment of the present invention shown in FIG. 5 except that a stereo pair setting means 42 is added, and the same components are the same. The code | symbol is attached | subjected. The operation of the same components is the same as that of the third embodiment of the present invention.
[0057]
That is, the image processing apparatus 4 includes a stereo pair setting unit 42 that selects a pair of images from a plurality of images (aerial image data a), an ortho image generation unit 41 that generates an ortho image d, and an image (aerial image). The projection transformation means 21 for projecting transformation of the data a) according to the orientation angle with respect to the orthographic projection plane, the projection transformation aerial image e and the ortho image d are synthesized to generate an interpolated ortho image f in which the missing pixels are interpolated. The image synthesizing unit 22 and a recording medium 14 for storing a program for realizing each unit of the image processing apparatus 2 are configured.
[0058]
Here, the ortho image generation means 41 performs a stereo matching process on a plurality of images (aerial image data a) and extracts the DEM data b, and a smoothing filter on the extracted DEM data b. Is applied to the DEM smoothing means 12 for removing noise and the orthorectifying means 13 for orthogonally transforming the image (aerial image data a) using the smoothed smoothed DEM data c.
[0059]
Each means of the image processing apparatus 4 operates generally as follows. That is, the following processing is performed by the computer that mainly configures the image processing apparatus 4 executing the program of the recording medium 14.
[0060]
The stereo pair setting means 42 determines a plurality of stereo image pairs obtained by selecting two of the three or more image data groups obtained by photographing the same spot, and sets it as a processing unit. The stereo processing means 11 performs a stereo matching process on each stereo image pair given by the stereo pair setting means 42 to extract DEM data b.
[0061]
The DEM smoothing means 12 smoothes the DEM data b obtained by the stereo processing means 11 by applying a smoothing filter. The orthorectifying means 13 orthorectifies the image (aerial image data a) using the smoothed DEM data c smoothed by the DEM smoothing means 12 described above.
[0062]
The projective conversion means 21 performs projective conversion of the image (aerial image data a) according to the rotation angle to the orthographic projection plane that has been determined. The image synthesizing means 22 interpolates the orthoimages d obtained by the orthorectifying means 13 and the image (projective transformation aerial image e) transformed by the projective transformation means 21 by superimposing the missing pixels. Is generated as an interpolated ortho image f.
[0063]
FIG. 8 is a flowchart showing the operation of the image processing apparatus 4 according to the fourth embodiment of the present invention. The operation of the image processing apparatus 4 according to the fourth embodiment of the present invention will be described with reference to FIGS. Note that the processing shown in FIG. 8 is realized by a computer that mainly configures the image processing apparatus 4 executing the program of the recording medium 14.
[0064]
First, the stereo pair setting means 42 determines a plurality of stereo image pairs obtained by selecting two of three or more image data groups obtained by photographing the same spot, and sets them as a processing unit (step S31 in FIG. 8). For example, when four images are obtained, one set of two of these sets and another set of stereo images are set with the remaining two. Although there are six combinations that can be selected, it is not always necessary to set all of them as stereo image pairs, and a plurality of stereo image pairs may be determined so that all images are used to reduce the concealment area.
[0065]
Next, the stereo processing means 11 performs a stereo matching process for each stereo image pair set by the stereo pair setting means 42, and extracts DEM data b (step S32 in FIG. 8). There is no restriction | limiting in particular in the stereo matching process used here, For example, it can implement | achieve by the method disclosed by Unexamined-Japanese-Patent No. 3-167678.
[0066]
Subsequently, the DEM smoothing unit 12 applies a smoothing filter to each DEM data b obtained by the stereo processing unit 11 to perform a smoothing process (step S33 in FIG. 8). The smoothing filter used here is not particularly limited as long as it has an effect of reducing noise. For example, a median filter, a multi-value morphological filter, or the like may be used. It is also possible to apply a combination of these filters.
[0067]
The orthorectifying means 13 orthorectifies the image (aerial image data a) with respect to a reference plane using each DEM data (smoothed DEM data c) smoothed by the DEM smoothing means 12, An ortho image d is created (step S34 in FIG. 8). The orthographic projection plane at this time is, for example, a plane parallel to the ground surface in the case of using a satellite image or an aerial photograph.
[0068]
The projective transformation means 21 is an image (projection-transformed aerial image) obtained by projective transformation of any one of the plurality of images (aerial image data a) used in the stereo processing means 11 according to the rotation angle to the orthographic projection plane. e) is generated (step S35 in FIG. 8). This is an operation of approximating an image (aerial image data a) with a plane and reprojecting it on an orthographic projection plane with the viewpoint fixed. As an image used at that time, an arbitrary image among the images used in the stereo processing unit 11 may be used, but an image having a small rotation angle with respect to the orthographic projection plane is desirable.
[0069]
Finally, the image synthesizing means 22 superimposes each ortho image d obtained by the orthorectifying means 13 and the projective transformation aerial image e obtained by the projective transformation converting means 21 so as to interpolate. An interpolated ortho ortho image f is generated (step S36 in FIG. 8). Specifically, interpolation is performed using the corresponding pixels in the other orthoimages for the pixels in the omission area where no value is obtained in the orthoimage d. Interpolation processing is performed on the pixels in the missing area remaining in this processing using the corresponding pixels in the projective transformation image.
[0070]
In the above-described fourth embodiment of the present invention, the processing for a plurality of stereo pairs is repeatedly performed in the stereo processing means 11, the DEM smoothing means 12, and the orthorectifying means 13. For example, the ortho image generating means described above is used. It is also possible to provide a plurality of 41 and perform processing for a plurality of stereo pairs in parallel.
[0071]
As described above, in the present embodiment, since the smoothed smoothed DEM data c is used for the orthorectification process, an ortho image in which the influence of noise included in the original DEM data is reduced can be generated.
[0072]
Further, in this embodiment, ortho images are generated by mutually complementing from a plurality of stereo pairs with different line-of-sight directions, and in addition, the image is projectively converted according to the rotation angle with respect to the orthographic projection plane, and the corresponding region of the image is displayed. Since pixels are substituted in the missing area of the ortho image, the omission area of the ortho image can be interpolated, and as a result, when the ortho image is applied to CG, the uncomfortable feeling caused by the omission area is reduced. It becomes possible.
[0073]
The present invention is not limited to each of the first to fourth embodiments of the present invention described above, and can be implemented under various conditions. For example, the target image is not limited to an aerial photograph or a satellite image, and can be applied to an image obtained by stereo shooting on the ground or indoors. The image may be taken using two cameras or taken twice using a monocular camera, and the image data may be taken directly with a digital camera or scanned from an analog photo. Even if it is, it can be implemented.
[0074]
Subsequently, the operation of the embodiment of the present invention will be described using a specific example. In this description, a case will be described in which an ortho image parallel to the ground surface is generated from a stereo aerial photograph according to the third embodiment of the present invention.
[0075]
Stereo aerial photographs are treated by treating overlapping portions of images continuously taken from an aircraft as stereo images, and are first converted into digital data using a scanner. Subsequently, orientation work is performed to determine camera parameters at the time of shooting. Normally, even in aerial photographs taken in a direction close to vertical, it is difficult to grasp the camera parameters at the time of shooting due to the effects of changes in the attitude of the aircraft, etc. The coordinate values in the image are read, and parameters such as the camera rotation angle between the left and right images are determined using the coordinate pairs of the corresponding points. There is no restriction | limiting in particular in the orientation method used here, For example, you may use the value calculated by plotting software.
[0076]
After the camera parameters are determined, as shown in FIG. 9, the left and right images are re-projected onto a common parallel plane as necessary. This is an operation called parallelization, and the same object in the left and right images after parallelization comes on the same scanning line.
[0077]
Further, the stereo processing unit 11 automatically performs stereo matching processing on the left and right images that have been made parallel, and obtains a correspondence between the left and right images to generate dense parallax images. There is no restriction | limiting in particular in the stereo matching process used here, For example, it can implement | achieve by the method disclosed by Unexamined-Japanese-Patent No. 3-167678.
[0078]
If the parallax image is obtained, the position in the three-dimensional coordinate system corresponding to each pixel is calculated in combination with the determined camera parameters, and DEM data is extracted. As a result, the obtained DEM data includes many spike-like isolated points having discontinuous three-dimensional coordinates as noise due to the influence of pixels or the like that are mismatched during stereo matching processing.
[0079]
In order to remove the influence of these isolated points, the DEM smoothing means 12 applies a smoothing filter to the DEM data obtained by the stereo processing means 11 to perform a smoothing process. The smoothing filter used here is not particularly limited as long as it has an effect of reducing noise. For example, a median filter that replaces the median of pixel values in a small area centered on the target pixel as a representative value is used. For example, the above isolated noise can be efficiently removed while maintaining an edge depending on the actual shape.
[0080]
In addition, a multi-value morphological filter configured by combining a maximum value filter and a minimum value filter for a small region can be expected to have the same effect as described above. It is also possible to apply a combination of a plurality of such filters.
[0081]
Subsequently, using the smoothed DEM data that has been smoothed by the DEM smoothing means 12 and from which the noise has been removed, and the camera parameters that have been standardized, the orthorectifying means 13 corrects the original aerial image with respect to the ground surface. Projective transformation is performed to create an ortho image.
[0082]
Specifically, first, the plane coordinates in the ground coordinate system (for example, the plane rectangular coordinate system) when the three-dimensional coordinates at all points in the DEM data are projected on the ground surface, and the coordinates in the original image Ask for action. Then, according to the obtained correspondence, the original aerial image pixels are mapped while being interpolated into an image corresponding to the ground surface.
[0083]
Further, the projective conversion means 21 performs an operation of reprojecting the aerial image data on a plane parallel to the ground surface by projective conversion. As described above, the optical axis of the camera at the time of aerial photography is not completely perpendicular to the ground surface, and since the ground surface is photographed from an oblique direction, it is difficult to superimpose with the ortho image. For this reason, the projective transformation is applied to the aerial image, and the operation of reprojecting onto the projection plane parallel to the ground without changing the camera center position is applied. By this operation, an image corresponding to a complete vertical photograph can be generated.
[0084]
A specific method of projective transformation is, for example, a method disclosed in Japanese Patent Laid-Open No. 8-567 using the rotation angle of the camera coordinate system with respect to the ground surface and the focal length among the camera parameters at the time of shooting. It can be carried out.
[0085]
Finally, the image synthesizing unit 22 synthesizes the ortho image generated by the orthorizing unit 13 and the projective transformation image generated by the projective transformation unit 21 to generate an interpolated ortho image. Specifically, after superposing the ortho image and the projective transformation image, the pixel of the projective transformation image at the corresponding position is cited for the pixel of the omission area in the ortho image where the value is not obtained due to the influence of concealment. Then, an interpolation ortho image is generated by performing an interpolation process.
[0086]
In this case, the image from which the pixel is cited is an image of the left and right projection-transformed images that has a principal point close to the target position, that is, the one that faces the target point as much as possible. By doing so, it is possible to reduce the case where the building wall surface that falls down due to plane approximation is used for the concealment region, and to improve the image quality.
[0087]
Further, by blurring the boundary portion or the quoting portion at the time of quoting, it is possible to suppress the influence of the generated edge, and to generate an interpolated ortho image with less sense of incongruity.
[0088]
Next, a specific example in the case where an ortho image parallel to the ground plane is generated from four stereo aerial photographs according to the example of the fourth embodiment will be described.
[0089]
For example, as shown in FIG. 10, four stereo aerial photographs are taken with a large overlap between two adjacent shooting courses # 1 and # 2, and each pair of images shooting the same point is taken. Can be obtained by selecting one set from each course. However, the photographing method is not limited to this, and four images # 1 to # 4 obtained by individually photographing may be used. The obtained aerial photograph group is first converted into digital data using a scanner.
[0090]
Next, the stereo pair setting means 42 sets a stereo pair from the four stereo aerial photographs converted into digital data. By referring to metadata or the like added to the image, two stereo pairs are obtained by combining two images in the same shooting course. However, the pair setting method is not limited to this, and if all four images # 1 to # 4 are used, image pairs may be set between different courses. Further, each image # 1 to # 4 may be used a plurality of times to generate a maximum of 6 stereo pairs.
[0091]
For each set stereo pair, as in the third embodiment of the present invention, orientation work for determining camera parameters at the time of photographing, parallelization, stereo matching processing, extraction of DEM data, and smoothing of DEM data are performed. Is called.
[0092]
Subsequently, using the DEM smoothing means 12 smoothed by the above DEM smoothing means 12 and noise-removed smoothed DEM data and the camera parameters that have been standardized, the orthorectifying means 13 generates the original aerial image for each aerial photograph. Orthographic transformation is performed on the ground surface to create an ortho image.
[0093]
Specifically, first, when the three-dimensional coordinates at all points in the DEM data are projected onto the ground surface, the plane coordinates in the ground coordinate system (for example, the plane rectangular coordinate system) and the coordinates in the original image Ask for action. Then, according to the obtained correspondence, the original aerial image pixels are mapped while being interpolated into an image corresponding to the ground surface.
[0094]
Further, the projective transformation means 21 performs an operation of reprojecting the aerial image data of each stereo pair onto a plane parallel to the ground surface by projective transformation. As described above, the optical axis of the camera at the time of aerial photography is not completely perpendicular to the ground surface, and since the ground surface is photographed from an oblique direction, it is difficult to superimpose with the ortho image. For this reason, projective transformation is applied to the aerial image, and the operation of reprojecting onto a projection plane parallel to the ground without changing the camera center position is applied. By this operation, an image corresponding to a complete vertical photograph can be generated.
[0095]
A specific method of projective transformation is, for example, a method described in Japanese Patent Application Laid-Open No. 8-567 using the rotation angle and focal length of the camera coordinate system with respect to the ground surface among the camera parameters at the time of shooting. Can do this.
[0096]
Finally, the image synthesizing unit 22 synthesizes each ortho image generated by the above orthorizing unit 13 and each projective transformation image generated by the projective transformation unit 21 to generate an interpolated ortho image. Do.
[0097]
Specifically, after each ortho image and each projective transformation image are overlapped, first, a plurality of ortho images are synthesized by interpolating pixels. In other words, the value is used for the point where a pixel is obtained only for a certain ortho image, and if the value is obtained for multiple ortho images, the main point is located close to the target position. Are synthesized using pixel values in an image to be photographed, that is, an image in which the pixel is photographed with the pixel facing as much as possible.
[0098]
Further, with respect to pixels in the omission region in the ortho image in which a value is not obtained due to the influence of concealment, an interpolation ortho image is generated by quoting the pixel of the projective transformation image at the corresponding position and performing an interpolation process. At this time, the image from which the pixel is cited is an image of the plurality of projective transformation images having a principal point close to the target position, that is, the one facing the target point as much as possible. As a result, it is possible to reduce the case where the building wall surface that falls down due to the plane approximation is used for the concealment region, and to improve the image quality. Further, by blurring the boundary portion or the quoting portion at the time of quoting, it is possible to suppress the influence of the generated edge, and to generate an interpolated ortho image with less sense of incongruity.
[0099]
As described above, according to the present invention, noise included in the original DEM data can be reduced by using the smoothed DEM data for the orthorectification process, and the image quality of the ortho image is reduced due to the noise. Therefore, a natural ortho image with less noise can be generated.
[0100]
In addition, the present invention interpolates a missing region of an ortho image by projective transforming the image in accordance with a rotation angle with respect to the orthographic projection plane, and substituting a corresponding region pixel of the image in a missing region of the ortho image. Therefore, it is possible to generate an ortho image with less sense of incongruity due to the missing area.
[0101]
【The invention's effect】
As described above, the image processing apparatus of the present invention extracts three-dimensional data by stereo matching processing in an image processing apparatus that performs stereo matching processing on a plurality of images to generate three-dimensional data. The three-dimensional data is included in the three-dimensional data by applying a smoothing filter to remove and smooth the noise of the three-dimensional data, and orthorectifying the image using the smoothed three-dimensional data. There is an effect that it is possible to obtain a more natural ortho image that is not easily affected by noise and that is less uncomfortable due to a defective region.
[0102]
According to another image processing apparatus of the present invention, in an image processing apparatus that performs stereo matching processing on a plurality of images to generate three-dimensional data, the three-dimensional data is extracted by stereo matching processing. The image is orthographically transformed, and the image is projectively transformed according to the orientation angle with respect to the orthographic projection plane, and the ortho image obtained by the orthographic transformation of the image and the projective transformation image obtained by the projective transformation are synthesized. By generating an interpolated ortho image in which the missing pixels are interpolated, it is possible to obtain a more natural ortho image that is less affected by noise included in the three-dimensional data and that has less uncomfortable feeling caused by the defective area. The effect is obtained.
[Brief description of the drawings]
FIG. 1 is a block diagram showing a configuration of an image processing apparatus according to a first embodiment of the present invention.
FIG. 2 is a flowchart showing the operation of the image processing apparatus according to the first embodiment of the present invention.
FIG. 3 is a block diagram showing a configuration of an image processing apparatus according to a second embodiment of the present invention.
FIG. 4 is a flowchart showing the operation of the image processing apparatus according to the second embodiment of the present invention.
FIG. 5 is a block diagram showing a configuration of an image processing apparatus according to a third embodiment of the present invention.
FIG. 6 is a flowchart showing the operation of the image processing apparatus according to the third embodiment of the present invention.
FIG. 7 is a block diagram showing a configuration of an image processing apparatus according to a fourth embodiment of the present invention.
FIG. 8 is a flowchart showing an operation of the image processing apparatus according to the fourth embodiment of the present invention.
FIG. 9 is a schematic diagram showing stereo shooting geometry of the present invention.
FIG. 10 is a schematic diagram showing an example of aerial photography taken from four directions in the embodiment of the present invention.
FIGS. 11A and 11B are schematic diagrams showing the difference between a central projection image and an orthogonal image. FIGS.
FIG. 12 is a block diagram illustrating a configuration of an image processing apparatus according to a conventional example.
FIG. 13 is a flowchart illustrating an operation of an image processing apparatus according to a conventional example.
FIG. 14 is a schematic diagram showing the occurrence of a missing area in an ortho image.
[Explanation of symbols]
1, 2, 3, 4 image processing apparatus
11 Stereo processing means
12 DEM smoothing means
13 Orthorectification means
14 Recording media
21 Projective transformation means
22 Image composition means
41 Ortho-image generation means
42 Stereo pair setting means

Claims (13)

複数の画像に対してステレオマッチング処理を行って3次元データを生成する画像処理装置であって、前記ステレオマッチング処理によって前記3次元データを抽出するステレオ処理手段と、前記ステレオ処理手段で抽出された前記3次元データに平滑化フィルタを適用して前記3次元データのノイズを除去して平滑化するDEM平滑化手段と、前記DEM平滑化手段で平滑化された3次元データを用いて前記画像に含まれる全画素個別に正射変換するオルソ化手段とを有することを特徴とする画像処理装置。An image processing apparatus for generating three-dimensional data by performing stereo matching processing on a plurality of images, wherein the stereo processing means extracts the three-dimensional data by the stereo matching processing, and is extracted by the stereo processing means and DEM smoothing means for smoothing to remove noise of the three-dimensional data by applying a smoothing filter to the three-dimensional data, the image using the three-dimensional data smoothed by the DEM smoothing means An image processing apparatus comprising orthorectifying means for orthogonally transforming all contained pixels individually . 前記画像を正射投影面に対する標定角にしたがって射影変換する射影変換手段と、前記オルソ化手段で得られたオルソ画像と前記射影変換手段で得られた射影変換画像とを合成して欠損した画素が補間された補間オルソ画像を生成する画像合成手段とを含むことを特徴とする請求項1記載の画像処理装置。  Projection conversion means for projective conversion of the image according to the orientation angle with respect to the orthographic projection plane, and a missing pixel by combining the ortho image obtained by the orthorectification means and the projection conversion image obtained by the projection conversion means The image processing apparatus according to claim 1, further comprising: an image synthesizing unit that generates an interpolated ortho image obtained by interpolating. 前記複数の画像の中から対になる画像を選択するステレオペア設定手段を含み、前記ステレオ処理手段と前記オルソ化手段と前記射影変換手段とが前記ステレオペア設定手段で選択された画像対に対して各々の処理を行うことを特徴とする請求項2記載の画像処理装置。  Stereo pair setting means for selecting a pair of images from the plurality of images, wherein the stereo processing means, the orthorectifying means, and the projective transformation means are for the image pair selected by the stereo pair setting means. The image processing apparatus according to claim 2, wherein each processing is performed. 複数の画像に対してステレオマッチング処理を行って3次元データを生成する画像処理装置であって、前記ステレオマッチング処理によって前記3次元データを抽出するステレオ処理手段と、前記3次元データを用いて前記画像を正射変換するオルソ化手段と、前記画像を正射投影面に対する標定角にしたがって射影変換する射影変換手段と、前記オルソ化手段で得られたオルソ画像と前記射影変換手段で得られた射影変換画像とを合成して欠損した画素が補間された補間オルソ画像を生成する画像合成手段とを有することを特徴とする画像処理装置。  An image processing apparatus that performs stereo matching processing on a plurality of images to generate three-dimensional data, wherein stereo processing means for extracting the three-dimensional data by the stereo matching processing, and using the three-dimensional data, the stereo processing means Orthorectification means for orthographic transformation of the image, projection transformation means for projective transformation of the image according to the orientation angle with respect to the orthographic projection plane, orthoimage obtained by the orthorectification means, and the projection transformation means An image processing apparatus, comprising: an image synthesizing unit that generates an interpolated ortho image in which a missing pixel is interpolated by synthesizing a projective conversion image. 前記画像合成手段は、複数のオルソ画像において値の得られていない画素を相互に補間するとともに、重複して値が得られる画素について値を持つ画像のうちの当該点に最も近い撮影主軸を持つオルソ画像を選択して合成し、複数生成されたオルソ画像のいずれにおいても値が得られない画素について当該点に最も近い撮影主軸を持つ射影変換画像から対応する画素を用いて補間することで前記補間オルソ画像を合成することを特徴とする請求項4記載の画像処理装置。  The image synthesizing means interpolates pixels for which values are not obtained in a plurality of ortho-images, and has a photographing main axis closest to the point in an image having values for pixels with which values are obtained in duplicate. By selecting and synthesizing an ortho image and interpolating using a corresponding pixel from a projective transformation image having a photographing principal axis closest to the point for a pixel for which no value is obtained in any of the plurality of generated ortho images The image processing apparatus according to claim 4, wherein the interpolated ortho image is synthesized. 画像処理装置において複数の画像に対してステレオマッチング処理を行って3次元データを生成する画像処理方法であって、前記画像処理装置に、前記ステレオマッチング処理によって前記3次元データを抽出するステップと、その抽出された前記3次元データに平滑化フィルタを適用して前記3次元データのノイズを除去して平滑化するステップと、平滑化された3次元データを用いて前記画像に含まれる全画素個別に正射変換するステップとを有することを特徴とする画像処理方法。An image processing method for generating three-dimensional data by performing stereo matching processing on a plurality of images in an image processing device, wherein the image processing device extracts the three-dimensional data by the stereo matching processing; Applying a smoothing filter to the extracted three-dimensional data to remove and smooth the noise of the three-dimensional data; and using the smoothed three-dimensional data, all pixels included in the image An image processing method comprising the step of performing orthogonal transformation individually . 前記画像処理装置に、前記画像を正射投影面に対する標定角にしたがって射影変換するステップと、前記画像の正射変換で得られたオルソ画像と前記射影変換で得られた射影変換画像とを合成して欠損した画素が補間された補間オルソ画像を生成するステップとを含むことを特徴とする請求項6記載の画像処理方法。  Projection conversion of the image according to the orientation angle with respect to the orthographic projection plane is performed on the image processing apparatus, and the ortho image obtained by orthographic conversion of the image and the projection conversion image obtained by the projective transformation are combined. 7. An image processing method according to claim 6, further comprising: generating an interpolated ortho image in which the missing pixel is interpolated. 前記画像処理装置に、前記複数の画像の中から対になる画像を選択するステップを含み、この選択された画像対に対して前記3次元データの抽出と前記画像の正射変換と前記射影変換とを行うことを特徴とする請求項7記載の画像処理方法。  The image processing apparatus includes a step of selecting a pair of images from the plurality of images, and extracting the three-dimensional data, orthonormal transformation of the image, and projective transformation for the selected image pair. The image processing method according to claim 7, wherein: 画像処理装置において複数の画像に対してステレオマッチング処理を行って3次元データを生成する画像処理方法であって、前記画像処理装置に、前記ステレオマッチング処理によって前記3次元データを抽出するステップと、前記3次元データを用いて前記画像を正射変換するステップと、前記画像を正射投影面に対する標定角にしたがって射影変換するステップと、前記画像の正射変換で得られたオルソ画像と前記射影変換で得られた射影変換画像とを合成して欠損した画素が補間された補間オルソ画像を生成するステップとを有することを特徴とする画像処理方法。  An image processing method for generating three-dimensional data by performing stereo matching processing on a plurality of images in an image processing device, wherein the image processing device extracts the three-dimensional data by the stereo matching processing; Orthorectifying the image using the three-dimensional data, projecting transforming the image in accordance with an orientation angle with respect to an orthographic projection plane, an ortho image obtained by orthographic transformation of the image, and the projection And a step of generating an interpolated ortho image in which the missing pixels are interpolated by synthesizing with the projective transformation image obtained by the conversion. 前記補間オルソ画像を生成するステップは、複数のオルソ画像において値の得られていない画素を相互に補間するとともに、重複して値が得られる画素について値を持つ画像のうちの当該点に最も近い撮影主軸を持つオルソ画像を選択して合成し、複数生成されたオルソ画像のいずれにおいても値が得られない画素について当該点に最も近い撮影主軸を持つ射影変換画像から対応する画素を用いて補間することによって前記補間オルソ画像を合成することを特徴とする請求項9記載の画像処理方法。  The step of generating the interpolated ortho image interpolates pixels for which values are not obtained in a plurality of ortho images, and is closest to the point in the image having values for pixels that are obtained with duplicate values. Select and synthesize an ortho image with a shooting principal axis, and interpolate using the corresponding pixel from the projective transformation image with the shooting principal axis closest to that point for pixels for which no value is obtained in any of the multiple generated ortho images The image processing method according to claim 9, wherein the interpolated ortho image is synthesized. 複数の画像に対してステレオマッチング処理を行って3次元データを生成する画像処理方法のプログラムであって、コンピュータに、前記ステレオマッチング処理によって前記3次元データを抽出する処理と、その抽出された前記3次元データに平滑化フィルタを適用して前記3次元データのノイズを除去して平滑化する処理と、平滑化された3次元データを用いて前記画像に含まれる全画素個別にを正射変換する処理とを実行させるためのプログラム。A program of an image processing method for generating three-dimensional data by performing stereo matching processing on a plurality of images, the computer extracting the three-dimensional data by the stereo matching processing, and the extracted 3-dimensional data to apply a smoothing filter processing for smoothing and removing noise of the three-dimensional data, orthographic and individually all the pixels included in the image using the three-dimensional data smoothed A program for executing conversion processing. 複数の画像に対してステレオマッチング処理を行って3次元データを生成する画像処理方法のプログラムであって、コンピュータに、前記ステレオマッチング処理によって前記3次元データを抽出する処理と、前記3次元データを用いて前記画像を正射変換する処理と、前記画像を正射投影面に対する標定角にしたがって射影変換する処理と、前記画像の正射変換で得られたオルソ画像と前記射影変換で得られた射影変換画像とを合成して欠損した画素が補間された補間オルソ画像を生成する処理とを実行させるためのプログラム。  A program of an image processing method for generating a three-dimensional data by performing a stereo matching process on a plurality of images, the computer extracting the three-dimensional data by the stereo matching process; Obtained by orthographic conversion of the image, orthographic image obtained by orthographic transformation of the image, and processing of projecting transformation of the image according to the orientation angle with respect to the orthographic projection plane. A program for executing a process of generating an interpolated ortho image in which a missing pixel is interpolated by synthesizing with a projective transformation image. 前記コンピュータに、前記補間オルソ画像を生成する処理を実行する際に、複数のオルソ画像において値の得られていない画素を相互に補間するととも、重複して値が得られる画素について値を持つ画像のうちの当該点に最も近い撮影主軸を持つオルソ画像を選択して合成し、複数生成されたオルソ画像のいずれにおいても値が得られない画素について当該点に最も近い撮影主軸を持つ射影変換画像から対応する画素を用いて補間することによって前記補間オルソ画像を合成することを特徴とする請求項12記載のプログラム。  When executing the process of generating the interpolated ortho image on the computer, an image having a value with respect to a pixel whose value is obtained redundantly while interpolating pixels whose values are not obtained in a plurality of ortho images Of these, an orthoimage having the photographing principal axis closest to the point is selected and synthesized, and a projected transformation image having the photographing principal axis closest to the point is obtained for a pixel for which no value is obtained in any of a plurality of generated orthoimages. 13. The program according to claim 12, wherein the interpolated ortho image is synthesized by performing interpolation using corresponding pixels.
JP2002323250A 2002-05-15 2002-11-07 Image processing apparatus, image processing method used therefor, and program therefor Expired - Fee Related JP3915664B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2002323250A JP3915664B2 (en) 2002-05-15 2002-11-07 Image processing apparatus, image processing method used therefor, and program therefor

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2002139378 2002-05-15
JP2002323250A JP3915664B2 (en) 2002-05-15 2002-11-07 Image processing apparatus, image processing method used therefor, and program therefor

Publications (2)

Publication Number Publication Date
JP2004046777A JP2004046777A (en) 2004-02-12
JP3915664B2 true JP3915664B2 (en) 2007-05-16

Family

ID=31719354

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2002323250A Expired - Fee Related JP3915664B2 (en) 2002-05-15 2002-11-07 Image processing apparatus, image processing method used therefor, and program therefor

Country Status (1)

Country Link
JP (1) JP3915664B2 (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4631606B2 (en) * 2005-03-07 2011-02-16 日本電信電話株式会社 Image depth distribution estimation method, apparatus, and program
KR20090052889A (en) * 2006-09-04 2009-05-26 코닌클리케 필립스 일렉트로닉스 엔.브이. Method for determining depth map from images and device for determining depth map
KR100912936B1 (en) 2008-11-27 2009-08-20 삼아항업(주) Digital Elevation Model Modification System and Method for Image Mapping
US8189943B2 (en) * 2009-03-17 2012-05-29 Mitsubishi Electric Research Laboratories, Inc. Method for up-sampling depth images
CN103714543A (en) * 2013-12-26 2014-04-09 南京理工大学 Simple tree dynamic programming binocular and stereo matching method based on invariant moment spatial information
JP6471428B2 (en) * 2014-08-13 2019-02-20 株式会社ソシオネクスト Image correction method, image correction apparatus, and image correction program

Also Published As

Publication number Publication date
JP2004046777A (en) 2004-02-12

Similar Documents

Publication Publication Date Title
US9959653B2 (en) Mosaic oblique images and methods of making and using same
US8855441B2 (en) Method and apparatus for transforming a non-linear lens-distorted image
JP4942221B2 (en) High resolution virtual focal plane image generation method
JP4658223B2 (en) Image generating method, apparatus, program thereof, and recording medium recording program
WO2015127847A1 (en) Super resolution processing method for depth image
JP2015022510A (en) Free viewpoint image capturing apparatus and method
JP3524147B2 (en) 3D image display device
WO2016125371A1 (en) Image generation device, coordinate conversion table creation device and creation method
JPH05303629A (en) Shape synthesis method
JP3915664B2 (en) Image processing apparatus, image processing method used therefor, and program therefor
CN115176284A (en) Image processing apparatus, image processing method, and image processing program
JP2011160354A (en) Image processing device, image processing method, and program therefor
Marí et al. L1B+: A perfect sensor localization model for simple satellite stereo reconstruction from push-frame image strips
JP7514125B2 (en) Image processing device and program thereof
JP3384764B2 (en) Image corresponding area extracting method, image synthesizing method, image synthesizing apparatus, computer-readable recording medium storing image synthesizing program
JPH0969170A (en) Method and apparatus for reconstructing shape and surface pattern of three-dimensional scene
JP2008282338A (en) Image processing apparatus, image processing program and method
KR20220092032A (en) Stitching system for panoramic image contets and stiching method usiting therewith
WO2002007093A2 (en) Method and apparatus for modeling via a three-dimensional image mosaic system
JP2000270206A (en) Device for compositing image of wide visual field and its method

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20061006

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20061017

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20061218

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: 20070116

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20070129

R150 Certificate of patent or registration of utility model

Ref document number: 3915664

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100216

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110216

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110216

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120216

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120216

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130216

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130216

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140216

Year of fee payment: 7

LAPS Cancellation because of no payment of annual fees