[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

JP2013254097A - 画像処理装置及びその制御方法、並びにプログラム - Google Patents

画像処理装置及びその制御方法、並びにプログラム Download PDF

Info

Publication number
JP2013254097A
JP2013254097A JP2012129865A JP2012129865A JP2013254097A JP 2013254097 A JP2013254097 A JP 2013254097A JP 2012129865 A JP2012129865 A JP 2012129865A JP 2012129865 A JP2012129865 A JP 2012129865A JP 2013254097 A JP2013254097 A JP 2013254097A
Authority
JP
Japan
Prior art keywords
image
distance image
distance
subject
images
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.)
Granted
Application number
JP2012129865A
Other languages
English (en)
Other versions
JP2013254097A5 (ja
JP6071257B2 (ja
Inventor
Masaaki Matsuoka
正明 松岡
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.)
Canon Inc
Original Assignee
Canon Inc
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 Canon Inc filed Critical Canon Inc
Priority to JP2012129865A priority Critical patent/JP6071257B2/ja
Priority to US13/910,644 priority patent/US10205928B2/en
Publication of JP2013254097A publication Critical patent/JP2013254097A/ja
Publication of JP2013254097A5 publication Critical patent/JP2013254097A5/ja
Application granted granted Critical
Publication of JP6071257B2 publication Critical patent/JP6071257B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N13/00Stereoscopic video systems; Multi-view video systems; Details thereof
    • H04N13/20Image signal generators
    • H04N13/204Image signal generators using stereoscopic image cameras
    • H04N13/207Image signal generators using stereoscopic image cameras using a single 2D image sensor
    • H04N13/218Image signal generators using stereoscopic image cameras using a single 2D image sensor using spatial multiplexing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/50Depth or shape recovery
    • G06T7/55Depth or shape recovery from multiple images
    • G06T7/571Depth or shape recovery from multiple images from focus
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/50Depth or shape recovery
    • G06T7/55Depth or shape recovery from multiple images
    • G06T7/593Depth or shape recovery from multiple images from stereo images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image
    • G06T2207/10012Stereo images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10141Special mode during image acquisition
    • G06T2207/10148Varying focus
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20221Image fusion; Image merging

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Studio Devices (AREA)
  • Automatic Focus Adjustment (AREA)

Abstract

【課題】結像光学系の絞り値や測定距離範囲の影響で、撮像する焦点位置の数が増大したり、距離精度が低下したりすることがない画像処理装置及びその制御方法、並びにプログラムを提供する。
【解決手段】画像処理装置100は、撮像手段により撮像された異なる視差を有する複数の画像から第1距離画像を生成し、撮像手段により撮像された画像内の各位置で合焦度が異なる複数の画像から第2距離画像を生成し、被写体像に示される複数の被写体から選択された主要被写体の被写界深度に応じて合成係数を算出し、第1距離画像と第2距離画像とを、算出手段により算出された合成係数を用いて合成することで合成距離画像を生成する。
【選択図】図9

Description

本発明は、画像処理装置及びその制御方法、並びにプログラムに関する。
従来より、コンピュータビジョン等の分野では、撮像した画像から距離情報を得ることが求められている。
例えば、視差を有した複数の画像から相関を用いて対応点探索を行って距離情報を推定する技術が開示されている(例えば、特許文献1、非特許文献1参照)。
また、焦点位置の異なる複数の画像から合焦度の変化量を計算して距離情報を推定する技術が開示されている(例えば、特許文献2参照)。
特許文献1ではさらに、結像光学系の被写界深度から外れた対象物に対しても精度よく距離情報を推定するために、複数の焦点位置で画像を撮像することでそのような対象物が被写界深度から外れないようにして距離情報を推定する技術が開示されている。
図21は、複数の焦点位置を説明するための図である。図21において、距離測定装置9117は、視差を有した複数の画像を撮像する。距離測定装置9117の結像光学系の被写界深度は、以下の式で計算される。
Figure 2013254097
ここで、Dは前方被写界深度、Dは後方被写界深度、Dは被写界深度、rは許容錯乱円径、Avは絞り値、Zは被写体距離、fは焦点距離である。
例えば、許容錯乱円径を0.03328mm、絞り値を2.8、焦点距離を50mmとし、1m〜20mの距離範囲で被写界深度から外れる被写体が存在しないように複数の焦点位置を決定したとする。このとき、式1,2より、各々の焦点位置に対応する被写界深度の範囲は図21に示される矢印の範囲となる。図21における矢印9101〜9116は各被写界深度範囲の境界を示している。
特開平6−74762号公報 特開2010−183174号公報 Gary Bradski, Adrian Kaebler(著)、松田晃一(訳)、「詳解OpenCV」、オライリー・ジャパン、オーム社、2009年8月、p.377―466
しかしながら、上述の特許文献1に開示された従来技術では、距離測定装置の結像光学系の絞り値や距離測定する距離範囲によっては距離測定する焦点位置の数が増大してしまい、処理時間の増大を招いてしまう。
例えば図21に示されるような、距離測定する焦点位置が15か所の場合、デジタルカメラなど即時性が要求される装置に対し、従来技術は不向きである。
また、非特許文献1に開示された従来技術では、結像光学系の被写界深度から外れた対象物に対して精度よく距離情報を推定することができない。
さらに、特許文献2に開示された従来技術では、距離測定装置の結像光学系の絞り値や距離測定する距離範囲によっては被写界深度が深くなり、焦点位置の違いによる合焦度の変化量が小さくなってしまい、距離精度が低下してしまう。
そこで、本発明の目的は、結像光学系の絞り値や測定距離範囲の影響で、撮像する焦点位置の数が増大したり、距離精度が低下したりすることがない画像処理装置及びその制御方法、並びにプログラムを提供することにある。
上記目的を達成するために、請求項1の画像処理装置は、被写体像を結像する光学結像手段と、前記光学結像手段により結像された被写体像を撮像する撮像手段と、前記撮像手段により撮像された異なる視差を有する複数の画像から第1距離画像を生成する第1距離画像生成手段と、前記撮像手段により撮像された画像内の各位置で合焦度が異なる複数の画像から第2距離画像を生成する第2距離画像生成手段と、前記被写体像に示される複数の被写体から選択された主要被写体の被写界深度に応じて合成係数を算出する算出手段と、前記第1距離画像生成手段により生成された第1距離画像と前記第2距離画像生成手段により生成された第2距離画像とを、前記算出手段により算出された合成係数を用いて合成することで合成距離画像を生成する生成手段とを備えたことを特徴とする。
本発明によれば、結像光学系の絞り値や測定距離範囲の影響で、撮像する焦点位置の数が増大したり、距離精度が低下したりすることがない画像処理装置及びその制御方法、並びにプログラムを提供することができる。
本発明の実施の形態に係る画像処理装置の概略構成を示す図である。 図1における信号処理部の概略構成を示す図である。 被写界深度Dに対する合成係数を決定するグラフを示す図であり、(A)は連続的に変化するように合成係数を決定する場合のグラフを示し、(B)合成係数を2値のいずれかに決定する場合のグラフを示している。 図2におけるステレオ方式距離画像生成部の概略構成を示す図である。 左右の視点が形成する撮像系をピンホールモデルで表現した図である。 ホモグラフィHおよびHの写像手順を説明するための図である。 撮像された画像を、ステレオ平行化された左右画像へ写像する手順を説明するための図である。 ステレオ対応点探索で用いるブロックマッチング手法を説明するための図である。 図1におけるシステム制御部により実行される合成距離画像生成処理の手順を示す図である。 画像上におけるAF測距枠の位置を説明するための図である。 図9におけるステップS105のステレオ方式距離画像生成処理の手順を示すフローチャートである。 図9におけるステップS106のDFD方式距離画像生成処理の手順を示すフローチャートである。 図12におけるステップS304の合焦評価値演算処理の手順を示すフローチャートである。 図13の合焦評価値演算処理の処理概要を説明するための図である。 被写界深度Dに対するK1を決定するグラフを示す図である。 K2を決定するグラフを示す図である。 画面上におけるAF測距枠の位置を説明するための図である。 図1におけるシステム制御部により実行される合成距離画像生成処理の手順を示す図である。 K2_C、K2_L、K2_Rを決定するグラフを示す図である。 K1_C、K1_LRを決定するグラフを示す図である。 複数の焦点位置を説明するための図である。
以下、本発明の実施の形態について図面を参照しながら詳述する。
[第1の実施の形態]
図1は、本発明の実施の形態に係る画像処理装置100の概略構成を示す図である。
図1において、結像光学系101はレンズや絞りなどからなり、フォーカス調節や焦点距離調整を行う。撮像素子102は光学像を電気信号に変換する。この撮像素子102には例えばCCDが用いられる。上記結像光学系101は、被写体像を結像する光学結像手段に対応する。また、撮像素子102は、結像光学系101により結像された被写体像を撮像する撮像手段に対応する。
A/D変換部103は撮像素子102から出力されたアナログ画像信号をデジタル画像信号に変換する。信号処理部104はA/D変換部103から出力されたデジタル画像信号に対して信号処理を施す。
表示部109は撮像素子102で撮像された画像などを表示する。AF評価値取得部108は結像光学系101のオートフォーカス(AF)調整用のAF評価値を取得する。メモリ制御部105はメモリ(例えばDRAM)106との間で画像信号データなどの記憶または読み出しを行う。
操作部110はユーザからの指示を受け付ける。システム制御部107は画像処理装置100全体の動作を制御する。
このように構成された画像処理装置100において、結像光学系101のフォーカス調整をオートフォーカスで実現するためには、一般によく知られている山登りコントラストAFの技術を利用すればよい。
つまり、結像光学系101を通過した被写体光束を撮像素子102で受光し、AF測距枠として指定した領域内で撮像出力のコントラスト(AF評価値)が極大となるようにフォーカス調整することでオートフォーカスを実現することができる。
また、フォーカス調整された結像光学系101の焦点位置から、AF測距枠中の合焦被写体までの距離Zを求めることができる。
図2は、図1における信号処理部104の概略構成を示す図である。
図2において、基準画像201はステレオ方式距離画像生成部210とDFD(Depth From Defocus)方式距離画像生成部224が参照する画像である。また、ステレオ方式参照画像202はステレオ方式距離画像生成部210が参照する画像であり、この画像は基準画像201に対して視差を有している。
DFD方式参照画像221,222は、DFD方式距離画像生成部224が参照する画像であり、これらの画像は基準画像201とは異なる焦点位置で撮像されている。
ステレオ方式距離画像生成部210は、視差を有する複数の画像から相関を用いて対応点探索を行うことで距離情報を推定し、ステレオ方式距離画像209を生成する。
DFD方式距離画像生成部224は、焦点位置の異なる複数の画像から合焦度の変化量を計算して距離情報を推定し、DFD方式距離画像223を生成する。
距離画像合成部243は、ステレオ方式距離画像(第1距離画像)209とDFD方式距離画像(第2距離画像)223を合成することで、合成距離画像241を生成する。
この距離画像合成部243において、画像の合成は、以下の式4を用いて行われる。
Figure 2013254097
ここで、(x,y)は画素位置、ZMIX(x,y)は合成距離画像の(x,y)での画素値、ZSTR(x,y)はステレオ方式距離画像の(x,y)での画素値、ZDFD(x,y)はDFD方式距離画像の(x,y)での画素値、Kは合成係数である。
距離画像合成部243は、コントラストAFでフォーカスが合った被写体を主要被写体として選択し、主要被写体の被写界深度に応じて合成係数Kを決定する。
例えば、許容錯乱円径を0.03328mm、絞り値を2.8、焦点距離を50mmとし、主要被写体の被写体距離を9mとすると、被写界深度は式1〜式3により、6.80mとなる。同様に、主要被写体の被写体距離を1mとすると、被写界深度は0.075mとなる。このように、まず被写体像に示される複数の被写体から選択された主要被写体の被写界深度に応じて合成係数が算出される。そして、第1距離画像と第2距離画像とを、算出された合成係数を用いて合成することで合成距離画像が生成される。
さらに、合成係数kは、主要被写体の被写界深度が増加するにしたがって単調増加する係数である。式4に示されるように、ステレオ方式距離画像及びDFD方式距離画像における位置(x,y)に応じて対応するステレオ方式距離画像の画素及びDFD方式距離画像の画素の、各々の画素値の加重平均を新たな画素値とした画像が合成距離画像として生成される。そして、当該加重平均において合成係数kはステレオ方式距離画像への画素値に対する重み係数である。
また、DFD方式距離画像への画素値に対する重み係数(1−k)は、合成係数kが増加するにしたがって単調減少する。
なお、上述した単調増加は、書物によっては広義の単調増加または単調非減少と表現されることがある。同様に単調減少は、書物によっては広義の単調減少または単調非増加と表現されることがある。
図3は、被写界深度Dに対する合成係数を決定するグラフを示す図であり、(A)は連続的に変化するように合成係数を決定する場合のグラフを示し、(B)合成係数を2値のいずれかに決定する場合のグラフを示している。
図3(A)(B)において、例えば主要被写体の被写界深度が6.80mであれば合成係数は1.0となり、式4から合成距離画像はステレオ方式距離画像優先となる。また、主要被写体の被写界深度が0.075mであれば合成係数は0.0となり、式4から合成距離画像はDFD方式距離画像優先となる。
さらに、図3(A)において、被写界深度が4.0m〜5.0mの場合は、DFD方式距離画像優先からステレオ方式距離画像優先へと徐々に変化する。
図4は、図2におけるステレオ方式距離画像生成部210の概略構成を示す図である。
図4におけるステレオ方式距離画像生成部210の説明に先立ち、ステレオ方式距離画像の生成について非特許文献1では、以下のステップでステレオ方式距離画像を生成することができることが示されている。
まず、左右の異なる視点で撮像された2つの画像に対して、結像光学系101により生じる歪み補正を行う。次いで、左右の異なる視点間でカメラの角度と距離を調整するために、左右画像を写像する(ステレオ平行化)。そして、左右のカメラのビューで同じ特徴を求めることで、視差画像を生成する(ステレオ対応点探索)。左右のカメラの幾何学的な位置情報を用いて、視差画像を距離画像に変換する(3D再投影)。
図4において、左右のカメラで撮像された画像は、基準画像201とステレオ方式参照画像202である。
また、ホモグラフィ推定部203はステレオ平行化のためのホモグラフィを推定する。ここでのホモグラフィ推定処理においては、結像光学系101による歪みが補正された座標系を用いる必要がある。
その座標変換式は以下の式5である。
Figure 2013254097
ここで、(x,y)は歪まされた点の元の位置であり、(xcorrected,ycorrected)は補正の結果得られる新しい位置である。さらに、rは歪まされた点の元の位置におけるレンズ中心からの距離であり、k1、k2、k3はレンズの形状から決定される歪み補正係数である。
ここでの歪み補正は半径方向歪みのみを補正するものとしたが、それだけに限定されるものではなく、例えば円周方向歪みも補正するようにしてもよい。
図5は、左右の視点が形成する撮像系をピンホールモデルで表現した図である。
図5においてO及びOは左右の視点の投影中心位置、Pは被写体位置、p及びpは左右の視点の投影平面に投影される被写体の結像位置、e及びeは左右の視点におけるエピポール、p、e及びp、eは左右の視点におけるエピポーラ線である。
エピポーラ拘束により、1つの画像内の特徴点が与えられると、もう1つの画像内で対応する点は対応するエピポーラ線上に存在する。
このため、ステレオ対応点探索はエピポーラ線に沿って対応点探索することになるが、2つの画像でエピポーラ線を平行にして対応点の高さを合わせると水平方向の1次元での対応点探索となる。この場合、2次元で探索するよりも計算量が大幅に軽減されるし、結果の信頼度も高い。
2つの画像でエピポーラ線を平行にして対応点の高さを合わせる座標変換が、ステレオ平行化である。ステレオ平行は以下の順に処理が実行される。
(1−1)2つの画像で対応する8点以上の特徴点を求める。
(1−2)対応する8点以上の特徴点から基礎行列Fを求める。
(1−3)基礎行列Fからエピポールeおよびeを求める。
(1−4)エピポールeから右画像のステレオ平行化のためのホモグラフィHを推定する。
(1−5)対応する特徴点、ホモグラフィH、およびエピポールeから左画像のステレオ平行化のためのホモグラフィHを推定する。
上記(1−1)の処理では、2つの画像間でブロックマッチング手法による2次元の対応点探索を行うことで、必要な対応する特徴点を求める。画素単位で対応点探索するステレオ対応点探索とは異なり、画像一面内で対応する特徴点を最低8点だけ求めればよいので、2次元であっても計算量は少なく、結果の信頼度にも大きな影響はない。
ブロックマッチング手法は、後述するステレオ対応点探索の手法と同様である。ただしステレオ対応点探索のブロックマッチング手法は1次元方向のみの探索なので、ここではその探索範囲を2次元方向に拡張したブロックマッチング手法が適用される。
次の(1−2)の処理では、8点の特徴点から3×3行列の基礎行列Fを求める。また、特徴点の座標は式5で歪み補正されたものを用いる。以下の式6は一般にエピポーラ方程式と呼ばれている方程式である。
Figure 2013254097
ここで、Eは、3×3行列の基本行列である。また3×3行列のカメラ内部行列Mを用いると、一般に以下の式7で3次元物理座標pを投影平面の2次元ピクセル座標qに変換できることが知られている。
Figure 2013254097
またカメラ内部行列Mを用いると、一般に以下の式8で基本行列Eを基礎行列Fに変換できることが知られている。
Figure 2013254097
上述した式6から式8を用いると、基礎行列Fに関する以下の式9を導くことができる。
Figure 2013254097
基礎行列Fは3×3行列で、かつスカラの自由度を持つので、9個の行列要素のうちの8個が変数となる。このため式9に左右画像において対応する特徴点のq及びqの8点分を代入すれば、単なる連立1次方程式となり、この連立1次方程式を解くことで基礎行列Fを求めることができる。
このように、本実施の形態では、基礎行列Fを求めるのに左右画像において対応する8店の特徴点のq及びqを用いているが、これに限定されるものではない。例えば非特許文献1で開示されているように8点以上の特徴点を用いた最小二乗法によって基礎行列Fを推定してもよい。
次の(1−3)の処理では、基礎行列Fからエピポールe及びeを求める。基礎行列Fとエピポールe及びeの関係は一般に以下の式10で表されることが知られており、この式10に求めた基礎行列Fを代入することでエピポールe及びeを求めることができる。
Figure 2013254097
(1−4)の処理では、エピポールeからホモグラフィHを推定する。
図6は、ホモグラフィHおよびHの写像手順を説明するための図である。
図6において、網掛けされた三角形は、右画像の投影平面におけるオブジェクトの結像イメージを示している。
そして、図6(A)では、エピポールeのピクセル座標(x,y)が示されている。まず、エピポールeをx軸上の(k,0)に写像する回転行列Rを求める。つまり、以下の式11のようになる。
Figure 2013254097
式11をエピポールeに適用すると、以下の式12のようになる。
Figure 2013254097
図6(B)は、(A)に対して回転行列Rを適用したピクセル座標系を示している。次に、x軸上に写像されたエピポール(k,0)を無限遠点(±∞,0)に写像するホモグラフィGを求める。つまり、ホモグラフィGは以下の式13のようになる。
Figure 2013254097
このようなホモグラフィは、すべてのエピポーラ線をx軸に対して平行な直線に写像する。
図6(C)は、(B)に対してホモグラフィGを適用したピクセル座標系を示している。(C)におけるエピポーラ線はx軸に対して平行な直線に写像されている。
式11から式13により、ホモグラフィHは以下の式14のようになる。
Figure 2013254097
(1−5)の処理では、特徴点の対応、ホモグラフィH、およびエピポールeから左画像のステレオ平行化のためのホモグラフィHを推定する。このようなホモグラフィHは、左画像のすべてのエピポーラ線をx軸に対して平行な直線に写像し、左右画像の対応するエピポーラ線の高さを揃える。
左画像のすべてのエピポーラ線をx軸に対して平行な直線に写像する手順は、(1−4)の処理と同様の手順で実現することができる。また同様の手順で左画像に対して求められた行列を、G’及びR’とする。
左右画像の対応するエピポーラ線の高さを揃えることは、左右画像間に存在し、対応する特徴点において、特徴点間距離の総和が最小になるようなホモグラフィKを推定することで実現することができる。つまり、推定したいホモグラフィHは以下の式15ようになる。
Figure 2013254097
さらに、ホモグラフィKを推定するために、以下の式16を用いる。
Figure 2013254097
なお式16のdd(Q,R)は、2つのピクセル座標QR間の距離自乗和を表している。また対応する特徴点の組は(q ,q )は、(1−1)の処理で求められたものを用いればよい。
図4の説明に戻り、逆写像処理部204,205は、ステレオ平行化のためのホモグラフィを用いて、左画像(基準画像201)および右画像(ステレオ方式参照画像202)をステレオ平行化する。このステレオ平行化のための逆写像処理の手法は非特許文献1に示される手法と同様である。
図7は、撮像された画像(図では生画像)を、ステレオ平行化された左右画像へ写像する手順を説明するための図である。
図7において、逆写像処理は、(C)に示されるステレオ平行化された画像内の各整数座標に対して、(B)に示される歪み補正された画像内の座標を求め、それらを用いて(A)に示される生画像内の小数点座標を参照することで実現できる。
つまり、この小数点座標を用いて生画像内で画素補間処理を行い、(C)に示されるステレオ平行化画像内での画素値を求める。
なお、図中のDIST()は、式5の(xcorrected,ycorrected)から(x,y)を求める式(2)の逆変換処理である。また、(D)は切り出された視差画像を示している。
図4の説明に戻り、ステレオ対応点探索部206はステレオ平行化された左画像と右画像から視差画像を生成する。ステレオ対応点探索の手法は非特許文献1に示される手法と同様である。
図8は、ステレオ対応点探索で用いるブロックマッチング手法を説明するための図である。
ブロックマッチングは差分絶対値和(SAD)ウィンドウを行方向にスライドすることで行われる。具体的にブロックマッチングは、左の画像内でそれぞれの特徴に対して、右の画像で対応する行から最もよくマッチするものを探す。
左右画像はステレオ平行化されていて各行がエピポーラ線になっており、右の画像内でマッチングする箇所は左の画像内の同じ行にあるため、探索方向は図8に示されるように、例えば(x,y)から行方向だけに限定することができる。
左の画像内のある特徴点のx座標xに対して、対応する右の画像内の対応点のx座標xまでの距離d=(x−x)が、x座標xに対応する視差量となる。
図8では、基準画像201が左の画像の場合を示しているが、基準画像が右の画像の場合は、右の画像内の特徴点を基準に視差量を計算する。つまり、右の画像内のある特徴点のx座標xに対して、対応する左の画像内の対応点のx座標xまでの距離d=(x−x)を、x座標xに対応する視差量とする。
図4の説明に戻り、3D再投影部207は視差画像からステレオ方式距離画像209を生成する。
3D再投影の手法は非特許文献1の手法と同様である。三角測量の手法により、距離Zと視差量dの関係は一般に以下の式17で表されることが知られている。
Figure 2013254097
ここで、fは撮像系の焦点距離、Tは左右2つの視点の投影中心の距離である。この式17を用いて3D再投影部207により生成されたステレオ方式距離画像209は、メモリ制御部105によってメモリ106に記憶される。
以上説明した構成により実行される処理について、以下フローチャートを用いて説明する。
図9は、図1におけるシステム制御部107により実行される合成距離画像生成処理の手順を示す図である。
図9において、まずステレオ方式参照画像を撮像し、更新する(ステップS101)。焦点位置は、ステップS101の撮像前に予めAF評価値取得部108によるコントラストAFで調整しておく。図10は画像上におけるAF測距枠の位置を説明するための図であり、AF測距枠301が示されている。
このステップS101の処理は、ユーザによって撮像指示が操作部110で操作されるまで行われる。
撮像指示が操作部110で操作されると(ステップS102でYES)、基準画像を撮像し、更新する(ステップS103)。このようにユーザにより撮像指示が操作されるまでの間、ユーザはフレーミング動作で画像処理装置100の位置を移動させている。そして、撮像指示が操作された後は画像処理装置100の位置を確定させるため、撮像指示が操作された前後で撮像したステレオ方式参照画像と基準画像との間にはユーザのフレーミング動作に応じた視差が発生している。
続いて、焦点位置を変えて、DFD方式参照画像を撮像する(ステップS104)。撮像時の焦点位置は、被写体距離が1mと20mに対応する焦点位置としている。
次いで、ステレオ方式距離画像生成部210によるステレオ方式距離画像生成処理によって、異なる視差を有する複数の画像からステレオ方式距離画像が生成される(ステップS105)。このステップS105は、撮像手段により撮像された異なる視差を有する複数の画像からステレオ方式距離画像を生成する第1距離画像生成手段に対応する。
さらに、DFD方式距離画像生成部224によるDFD方式距離画像生成処理によって、撮像手段により撮像された画像内の各位置で合焦度が異なる複数の画像からDFD方式距離画像が生成される(ステップS106)。このステップS106は、撮像手段により撮像された画像内の各位置で合焦度が異なる複数の画像からDFD方式距離画像を生成する第2距離画像生成手段に対応する。
そして、生成されたステレオ方式距離画像とDFD方式距離画像は、距離画像合成部243により合成され、合成距離画像が生成され(ステップS107)、本処理を終了する。ここでは、上述したように合成係数が算出され、合成距離画像を生成するので、このステップS107は、算出手段及び生成手段に対応する。
図11は、図9におけるステップS105のステレオ方式距離画像生成処理の手順を示すフローチャートである。
図11において、まず視差画像を生成する(ステップS201)。上述したように、視差画像は、逆写像処理部204,205でステレオ平行化された左右画像を、ステレオ対応点探索部206で処理することで生成される。
次いで、AF測距枠301でオートフォーカスされた被写体の被写体距離Zと視差量dから、投影中心の並進移動量を取得する(ステップS202)。上記式17から、並進移動量T=(Z/f)となる。なお、被写体距離Zは、ステップS101の撮像前に予めAF評価値取得部108によるコントラストAFで算出されている。また、視差量dは、ステップS201の視差画像から抽出される。
そして、視差画像の画素ごとに式17を適用することで距離画像を生成し(ステップS203)、本処理を終了する。
図12は、図9におけるステップS106のDFD方式距離画像生成処理の手順を示すフローチャートである。
図12において、DFD方式距離画像生成部224は、基準画像及び参照画像のうちから、画像の合焦度を示す合焦評価値を求めるための処理に用いる処理対象画像を指定する(ステップS303)。なお、ここではDFD方式距離画像生成部224は、基準画像及び参照画像のすべてを結像光学系101の焦点位置及び絞り値の一方が異なる画像からなる処理対象画像として順次指定してゆくものとする。
次いで、DFD方式距離画像生成部224は、合焦評価値演算処理を実行し、処理対象画像について画像内の各領域または位置での合焦評価値をそれぞれ求める(ステップS304)。
ここで、画像内の各領域とは、撮像画像を複数の領域に(XY平面上で)2次元的に分割した領域を指すものとする。各領域の最小単位は1画素となる。なお、画像内の位置とは、その処理対象画像を撮像した撮像素子内に含まれる画素エリア(または単一の画素)の、撮像素子内でのXY位置を示すものとする。合焦評価値を求める領域または位置は、全ての処理対象画像に渡り共通している。以下の説明では、領域または位置を、単に領域と表現する。
またステップS304では、DFD方式距離画像生成部224は、処理対象画像での各領域または位置と合焦評価値との対応関係を示す合焦評価値マップを生成する。各々の合焦評価値は、画像の合焦度が高いほど高い値を示すものとする。
次いで、DFD方式距離画像生成部224は、全ての処理対象画像で合焦評価値演算処理が終了したか否かを判別する(ステップS305)。
ステップS305の判別の結果、全ての処理対象画像で合焦評価値演算処理が終了していないときは(ステップS305でNO)、上記ステップS303に戻る。
一方、ステップS305の判別の結果、全ての処理対象画像で合焦評価値演算処理が終了したときは(ステップS305でYES)、以下の処理を行う。
すなわち、DFD方式距離画像生成部224は、基準画像に対応する第1の合焦評価値マップと、参照画像に対応する第2の合焦評価値マップと、各々の画像の撮像時の焦点位置情報とを用いる。これらを用いることで、DFD方式距離画像生成部224は、基準画像内の各領域に対応する被写体距離をそれぞれ推測し(ステップS306)、本処理を終了する。
具体的にDFD方式距離画像生成部224は、第1の合焦評価値マップと第2の合焦評価値マップとの間で対応する画像上の領域ごとに合焦評価値の差分をとって、各々の領域で合焦評価値の変化量を求める。なお、参照画像が複数存在する場合には、DFD方式距離画像生成部224は各々の参照画像と基準画像との間で合焦評価値の変化量を求めるものとする。
ここで、基準画像内のある領域において、注目被写体が合焦状態であって、参照画像がぼけているケースを考える。このケースでは、第1の合焦評価値マップで注目被写体に対応する領域の合焦評価値は高い値を示すが、第2の合焦評価値マップで注目被写体に対応する領域の合焦評価値は第1の合焦評価値マップと比べて非常に低い値を示す。
そのため、このケースでは、注目被写体が対応する領域の合焦評価値の変化量は大きな値を示すこととなる。
一方、基準画像及び参照画像の両方で注目被写体がぼけているケースでは、第1の合焦評価値マップ及び第2の合焦評価値マップの両方で、注目被写体に対応する領域の合焦評価値の値がいずれも小さくなる。
そのため、このケースでは、注目被写体が対応する領域の合焦評価値の変化量は相対的に小さな値を示す。
よって、基準画像での合焦評価値が高く、参照画像では合焦評価値が大きく減少する領域は、焦点位置の移動により合焦状態から非合焦状態となった箇所と推測できる。
逆に、基準画像での合焦評価値が低く、参照画像では合焦評価値が大きくなる領域は、参照画像の焦点位置で合焦状態となった箇所と推測できる。
このとき、基準画像と参照画像との間で焦点位置の移動方向とその移動量とを参照すれば、注目する被写体の存在する領域が、基準画像の焦点位置を基準として前ピンまたは後ピンのいずれかにあるかが判別できる。また基準画像の焦点面とどの程度離れているのかを、DFD方式距離画像生成部224が判別できる。
同様に、基準画像での合焦評価値が低く、参照画像との間で合焦評価値の変化量も小さい領域は、基準画像での焦点位置と参照画像での焦点位置とのいずれもが合焦面から遠い箇所と推測できる。
また、基準画像での合焦評価値が高く、参照画像との間で合焦評価値の変化量が小さい領域は、基準画像での焦点位置および参照画像での焦点位置の両方でほぼ合焦状態にある箇所と推測できる。
上述したステップS306でのDFD方式距離画像生成部224は、上記の推測結果に基づいて、基準画像内の各領域に対応する被写体距離の情報をそれぞれ生成する。
例えば、DFD方式距離画像生成部224は、合焦評価値の高さと変化量とに基づいて、基準画像の焦点位置で合焦した被写体との相対距離またはその度数で各々の被写体距離を表現してもよい。
または、DFD方式距離画像生成部224は、画像の各位置で合焦状態となる画像を特定し、合焦状態となる画像に対応する被写体距離の値をそのまま被写体距離を表現してもよい。
図13は、図12におけるステップS304の合焦評価値演算処理の手順を示すフローチャートである。
図13においては、DFD方式距離画像生成部224が、処理対象画像に各々の周波数特性が異なる一対のエッジ抽出処理を施すとともに、これらのエッジ抽出処理で求まる一対のエッジ量の差分を用いて合焦評価値を求める。
まず、DFD方式距離画像生成部224は、処理対象画像に第1のエッジ抽出処理を施し、抽出したエッジ成分の強度と画像での位置との対応関係を示す第1エッジ画像を生成する(ステップS401)。
次いで、DFD方式距離画像生成部224は、ステップS401での抽出結果よりも低域側の周波数成分が含まれるように、処理対象の画像に第2のエッジ抽出処理を施す。そして、DFD方式距離画像生成部224は、抽出したエッジ成分の強度と画像での位置との対応関係を示す第2エッジ画像を生成する(ステップS402)。
そして、DFD方式距離画像生成部224は、第1エッジ画像及び第2エッジ画像を、それぞれ例えば20画素×20画素程度のサイズに領域分割するとともに、各領域ごとにそれぞれ特徴量を求める(ステップS403)。
次いで、DFD方式距離画像生成部224は、第1エッジ画像の各領域の特徴量と、第2エッジ画像の各領域の特徴量とを、対応する領域ごとに差分を求めることで(ステップS404)、その変位を求める。
DFD方式距離画像生成部224は、ステップS404で求めた各領域の特徴量の差分値を、第1エッジ画像または第2エッジ画像から求めたエッジ量を用いて正規化する(ステップS405)。
DFD方式距離画像生成部224は、ステップS405で求めた全ての正規化された特徴量の平均値を求め、この平均値を正規化された特徴量の基準値として設定する(ステップS406)。ただし、既に別の入力画像を用いて、この基準値が求まっている場合には、DFD方式距離画像生成部224は、その値を基準値としてもよい。
DFD方式距離画像生成部224は、各々の領域の正規化された特徴量(ステップS405)と正規化された特徴量の基準値(ステップS406)とを比較し、各領域の合焦評価値を求め(ステップS407)、本処理を終了する。
以上説明した各ステップでの処理に関して、図14を用いて説明する。
図14は、図13の合焦評価値演算処理の処理概要を説明するための図である。
まず、ステップS401について、DFD方式距離画像生成部224は、以下の(A1)、(A2)の各処理を実行する。
(A1)DFD方式距離画像生成部224は、処理対象画像に対して、ガウシアン型平滑化フィルタによるフィルタ演算を行う。一例として、ステップS401でのDFD方式距離画像生成部224は、P(r)∝exp(−r/2σ)の平滑化フィルタを用いるものとする。
なお、ステップS401での平滑化フィルタの強度とフィルタ半径とに関係するパラメータσの値は任意に設定できる。この処理により、図14に示されるように、元画像から平滑化された第1ボケ画像が生成される。以下の説明では、パラメータσを「フィルタ半径」という。
(A2)DFD方式距離画像生成部224は、元画像の各画素値から第1ボケ画像で対応する位置にある画素値をそれぞれ減算するとともに、減算後の画素値の絶対値をとる。これにより、DFD方式距離画像生成部224は、元画像を平滑化したときに失われるエッジ成分を抽出することで、図14に示されるような第1エッジ画像を生成する。
ステップS402で生成される第2エッジ画像も同様に生成されるが、ステップS401におけるA1の処理と比較して、平滑化フィルタのフィルタ半径σが大きい平滑化が行なわれる。
従って、ステップS402では、A1の処理と比べて平滑化の度合いが強くなる。すなわち、第2ボケ画像は、第1ボケ画像と比べて高域側の周波数成分がより少なく、また低域側の周波数成分がより多く含まれるので、第1ボケ画像よりもエッジがなまった状態となる。
ここで、画像のピントが合っていない領域では、フィルタ半径σの値を少し変えてもエッジ量はほとんど変化しない。
一方、画像のピントが合っている領域では、フィルタ半径σの値を少し変えただけでもエッジ量は大きく変化する。
そこで、ステップS401,402で、フィルタ半径σの値を互いに近い値に設定すると、後述の特徴量の差分をとったときに合焦領域と非合焦領域とでエッジ量の変位の振る舞いが大きく異なり、合焦度の測定が比較的容易となる。
一例として、ステップS401ではフィルタ半径σを0.5とし、ステップS402ではフィルタ半径σの値を0.6としてもよい。
次に、ステップS403の処理の一例について説明する。ステップS403において、合焦度算出部27は、各エッジ画像の領域毎にその領域に含まれる全画素のエッジ量を積算し、この積算値を領域の画素数で割った値を特徴量として算出する。すなわち、ここでの特徴量は領域内のエッジ量の平均値である。これにより、各領域の特徴をそれぞれ1つの値で表すことができる。
次に、ステップS404の処理について説明する。上述したように、ピントが合っていない領域では、フィルタ半径σの値を少しだけ変えてもエッジ量はほとんど変化しないが、画像のピントが合っている領域では、フィルタ半径σの値を少し変えただけでもエッジ量は大きく変化する。
そのため、平滑化フィルタのフィルタ半径を小さい範囲で変化させたときの各領域におけるエッジ量の変位をみれば、領域ごとの合焦度の強弱を判定できることが分かる。
次に、ステップS405の処理について説明する。例えば、DFD方式距離画像生成部224は、領域ごとの差分値を、第1エッジ画像で対応する領域の特徴量で割ることにより、差分値を正規化してもよい。
なお、ここでの特徴量として、第2エッジ画像で対応する領域の特徴量や、同じ領域における第1エッジ画像の特徴量と第2エッジ画像の特徴量との平均値を用いてもよい。
そして、DFD方式距離画像生成部224は、S405の正規化後の値を各々の領域における合焦量とする。
上記のステップS404の段階では、画像に含まれるエッジ量の値が元々大きい箇所でのエッジの変位と、画像に含まれるエッジ量の値が小さい箇所でのエッジの変位をそのまま比べると、前者の方の変位量が大きくなる可能性が高い。なお、エッジ量の値が大きい箇所は、高輝度なので画素値が大きい箇所であり、エッジ量の値が小さい箇所は、低輝度なので画素値が小さい箇所となっている。
すなわち、画像はボケていても輝度差が大きな箇所では合焦度を高く評価してしまう可能性がある。そのため、ステップS405でのDFD方式距離画像生成部224は、エッジ画像のエッジ量で差分値を正規化することで、画像の輝度差が合焦度に及ぼす影響を抑制している。
次に、ステップS407の処理について説明する。ステップS407では、各々の領域の正規化された特徴量と正規化された特徴量の基準値とを比較するが、この比較には、例えば、合焦量と基準値との減算値または除算値を用いてもよい。
このようにして合焦評価値を求めることにより、DFD方式距離画像生成部224は合焦評価値マップを取得できる。なお、異なる画像から求めた合焦評価値マップを比較する場合、DFD方式距離画像生成部224は、画像間の輝度差を考慮して必要に応じて合焦評価値マップの正規化を行うものとする。
以上説明した第1の実施の形態によれば、結像光学系の絞り値や測定距離範囲の影響で、撮像する焦点位置の数が増大したり、距離精度が低下したりすることがない。
[第2の実施の形態]
第2の実施の形態では、合成係数Kの決定方法が第1の実施の形態と異なっている。
第2の実施の形態における距離画像合成部243は、式4を用いて距離画像を合成することは第1の形態と同じであり、合成係数Kは以下の式18で決定する。
Figure 2013254097
距離画像合成部243は、コントラストAFでフォーカスが合った被写体を主要被写体として選択し、主要被写体の被写界深度に応じて合成係数K1を決定する。
例えば、許容錯乱円径を0.03328mm、絞り値を2.8、焦点距離を50mmとし、主要被写体の被写体距離を9mとすると、被写界深度は式1〜式3により、6.80mとなる。
図15は、被写界深度Dに対するK1を決定するグラフを示す図である。
図15において、例えば主要被写体の被写界深度が6.80mであれば合成係数K1は1.0となる。また、被写界深度が4.0m〜5.0mの場合は、K1の値が0から1へと単調増加する。
また、距離画像合成部243は、コントラストAFでフォーカスが合った被写体を主要被写体として選択し、主要被写体の被写界深度の範囲に応じて合成係数K2を決定する。
図16は、K2を決定するグラフを示す図である。
図16において、主要被写体402はコントラストAFでフォーカスが合った被写体である。
例えば、許容錯乱円径を0.03328mm、絞り値を2.8、焦点距離を50mmとし、主要被写体402の被写体距離を9mとすると、式1〜式3から、範囲403が主要被写体402の被写界深度の範囲となる。
合成係数K2を決定するために必要となる距離Zについては、ステレオ方式距離画像209を参照する。このように、被写界深度の範囲でとり得る値が、この範囲外でとり得る値と比較して大きくなる値を用いて合成係数Kが算出される。
以上説明した第2の実施の形態によれば、結像光学系の絞り値や測定距離範囲の影響で、撮像する焦点位置の数が増大したり、距離精度が低下したりすることがない。さらに、ステレオ方式距離画像について、主要被写体の被写界深度から外れた距離範囲がDFD方式距離画像に切り替わるので、そのような距離範囲の距離精度を向上させることができる。
[第3の実施の形態]
第3の実施の形態における画像処理装置100は、AF測距枠の数と位置、システム制御部107の動作を説明するためのフローチャート、距離画像合成部243が第1の実施の形態と異なっている。
図17は画面上におけるAF測距枠の位置を説明するための図である。
図17において、AF測距枠として、主要被写体用の中央測距枠501、非主要被写体用の左測距枠502、及び非主要被写体用の右測距枠503が示されている。第3の実施の形態では、これらの各AF測距枠が用いられる。
図18は、図1におけるシステム制御部107により実行される合成距離画像生成処理の手順を示す図である。
図18において、まずDFD方式参照画像を撮像する(ステップS501)。ここでの撮像での焦点位置は、被写体距離が1mと20mに対応する位置である。
次いで、主要被写体用の中央測距枠501でコントラストAFされた焦点位置で基準画像を撮像する(ステップS502)。
続けて、非主要被写体用の左測距枠502でコントラストAFされた焦点位置で基準画像を撮像する(ステップS503)。さらに非主要被写体用の右測距枠503でコントラストAFされた焦点位置で基準画像を撮像する(ステップS504)。
次いで、表示部109にユーザに視点の移動を促す内容を表示する(ステップS505)。そして、主要被写体用の中央測距枠501でコントラストAFされた焦点位置でステレオ方式参照画像を撮像する(ステップS506)。
続けて、非主要被写体用の左測距枠502でコントラストAFされた焦点位置でステレオ方式参照画像を撮像する(ステップS507)。さらに非主要被写体用の右測距枠503でコントラストAFされた焦点位置でステレオ方式参照画像を撮像する(ステップS508)。
そして、中央測距枠501の焦点位置の基準画像とDFD方式参照画像を用いて、DFD方式距離画像を生成する(ステップS509)。
さらに、主要被写体用の中央測距枠501におけるステレオ方式距離画像を生成する(ステップS510)。
続けて、非主要被写体用の左測距枠502におけるステレオ方式距離画像を生成する(ステップS511)。さらに非主要被写体用の右測距枠503におけるステレオ方式距離画像を生成する(ステップS512)。
生成された3つの測距枠の各焦点位置におけるステレオ方式距離画像とDFD方式距離画像を、距離画像合成部243を用いて合成して、合成距離画像を生成し(ステップS513)、本処理を終了する。
距離画像の合成は、以下の式19〜式24を用いて行われる。
Figure 2013254097
式19〜式24での記号の意味は以下の通りである。
(x,y):画素位置
MIX(x,y):合成距離画像
MIX_L(x,y)、ZMIX_R(x,y):左右の中間合成距離画像
STR_C(x,y):中央測距枠501のステレオ方式距離画像
STR_L(x,y):左測距枠502のステレオ方式距離画像
STR_R(x,y):右測距枠503のステレオ方式距離画像
DFD(x,y):DFD方式距離画像
K_C:中央測距枠501の合成係数
K_L:左測距枠502の合成係数
K_R:右測距枠503の合成係数
図19は、K2_C、K2_L、K2_Rを決定するグラフを示す図である。
図19において、主要被写体602は中央測距枠501におけるコントラストAFでフォーカスが合った被写体である。非主要被写体605は左測距枠、非主要被写体608は右測距枠におけるコントラストAFでフォーカスが合った被写体である。
例えば、許容錯乱円径を0.03328mm、絞り値を2.8、焦点距離を50mmとし、主要被写体602の被写体距離を9m、非主要被写体605の被写体距離を15m、非主要被写体608の被写体距離を1mとする。このとき、式1〜式3から、各被写体の被写界深度の範囲は範囲603,606,609となる。
図20は、K1_C、K1_LRを決定するグラフを示す図である。
非主要被写体の場合は主要被写体の場合に比べて、ステレオ方式優先となる被写界深度が主要被写体より深くなっている。
例えば、許容錯乱円径を0.03328mm、絞り値を2.8、焦点距離を50mmとし、主要被写体602の被写体距離を9mとすると、式1〜式3から、主要被写体602の被写界深度は、6.80mとなる。この場合、対応する合成係数K1_Cは、図20に示されるように1.0となる。
また、非主要被写体605の被写体距離を15m、非主要被写体608の被写体距離を1mとすると、式1〜式3から、各非主要被写体の被写界深度は、24.4m、0.075mとなる。この場合、非主要被写体605に対応する合成係数K1_LRは、1.0であり、非主要被写体608に対応する合成係数K1_LRは、0.0となる。このように、主要被写体とは異なる非主要被写体の被写界深度をさらに用いて合成係数が算出される。
また、第3の実施の形態では、非主要被写体に対応する他の第1距離画像として、左測距枠502のステレオ方式距離画像、及び右測距枠503のステレオ方式距離画像がさらに生成される。また、非主要被写体の被写界深度に応じて他の合成係数K_L、K_Rがさらに算出される。
そして、左測距枠502のステレオ方式距離画像、右測距枠503のステレオ方式距離画像とDFD方式距離画像を他の合成係数を用いて合成することで、左右の中間合成距離画像としてZMIX_L(x,y)、ZMIX_R(x,y)が生成される。
この左右の中間合成距離画像と中央測距枠501のステレオ方式距離画像とを合成係数K_Cを用いてさらに合成することで合成距離画像が生成される。
以上説明した第3の実施の形態によれば、結像光学系の絞り値や測定距離範囲の影響で、撮像する焦点位置の数が増大したり、距離精度が低下したりすることがない。
さらに、主要被写体の被写界深度から外れていて、かつ被写界深度が深くなる距離範囲において、DFD方式ではなくステレオ方式を優先するため、撮像する焦点位置の数を大きく増大させることなく、そのような距離範囲の距離精度を向上させることができる。
以上、本発明の好ましい実施の形態について説明したが、本発明はこれらの実施の形態に限定されず、その要旨の範囲内で種々の変形及び変更が可能である。
(他の実施の形態)
本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア(プログラム)をネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ(又はCPUやMPU等)がプログラムコードを読み出して実行する処理である。この場合、そのプログラム、及び該プログラムを記憶した記憶媒体は本発明を構成することになる。
100 画像処理装置
101 結像光学系
102 撮像素子
103 A/D変換部
104 信号処理部
105 メモリ制御部
106 メモリ
107 システム制御部
108 AF評価値取得部
109 表示部
110 操作部
201 基準画像
202 ステレオ方式参照画像
203 ホモグラフィ推定部
204,205 逆画像処理部
206 ステレオ対応点探索部
207 3D再投影部
208 視差画像
209 ステレオ方式距離画像
210 ステレオ方式距離画像生成部
221,222 DFD方式参照画像
223 DFD方式距離画像
224 DFD方式距離画像生成部
241 合成距離画像
243 距離画像合成部

Claims (10)

  1. 被写体像を結像する光学結像手段と、
    前記光学結像手段により結像された被写体像を撮像する撮像手段と、
    前記撮像手段により撮像された異なる視差を有する複数の画像から第1距離画像を生成する第1距離画像生成手段と、
    前記撮像手段により撮像された画像内の各位置で合焦度が異なる複数の画像から第2距離画像を生成する第2距離画像生成手段と、
    前記被写体像に示される複数の被写体から選択された主要被写体の被写界深度に応じて合成係数を算出する算出手段と、
    前記第1距離画像生成手段により生成された第1距離画像と前記第2距離画像生成手段により生成された第2距離画像とを、前記算出手段により算出された合成係数を用いて合成することで合成距離画像を生成する生成手段と
    を備えたことを特徴とする画像処理装置。
  2. 前記合成係数は、前記主要被写体の被写界深度が増加するにしたがって単調増加する係数であり、
    前記生成手段は、前記第1距離画像及び前記第2距離画像における位置に応じて対応する前記第1距離画像の画素及び前記第2距離画像の画素の、各々の画素値の加重平均を新たな画素値とした画像を前記合成距離画像として生成し、当該加重平均において前記合成係数は前記第1距離画像への画素値に対する重み係数であることを特徴とする請求項1記載の画像処理装置。
  3. 前記第2距離画像への画素値に対する重み係数は、前記合成係数が増加するにしたがって単調減少することを特徴とする請求項2記載の画像処理装置。
  4. 前記算出手段は、前記被写界深度の範囲でとり得る値が、その範囲外でとり得る値と比較して大きくなる値を用いて前記合成係数を算出することを特徴とする請求項2または請求項3記載の画像処理装置。
  5. 前記算出手段は、前記主要被写体とは異なる非主要被写体の被写界深度をさらに用いて前記合成係数を算出することを特徴とする請求項2乃至請求項4のいずれか1項に記載の画像処理装置。
  6. 前記第1距離画像生成手段は、前記非主要被写体に対応する他の第1距離画像をさらに生成し、
    前記算出手段は前記非主要被写体の被写界深度に応じて他の合成係数をさらに算出し、
    前記生成手段は、前記他の第1距離画像と前記第2距離画像とを前記他の合成係数を用いて合成することで生成した画像と、前記第1距離画像とを前記合成係数を用いてさらに合成することで前記合成距離画像を生成することを特徴とする請求項5記載の画像処理装置。
  7. 前記被写界深度は、被写体距離、前記光学結像手段と前記撮像手段で決定される焦点距離、絞り値、及び許容錯乱円径により定まることを特徴とする請求項1乃至請求項5のいずれか1項に記載の画像処理装置。
  8. 前記撮像手段で撮像された画像内の各位置で合焦度が異なる複数の画像は、前記光学結像手段の焦点位置及び絞り値の一方が異なる画像からなることを特徴とする請求項1乃至請求項6のいずれか1項に記載の画像処理装置。
  9. 被写体像を結像する光学結像手段と、前記光学結像手段により結像された被写体像を撮像する撮像手段とを備えた画像処理装置の制御方法であって、
    前記撮像手段により撮像された異なる視差を有する複数の画像から第1距離画像を生成する第1距離画像生成ステップと、
    前記撮像手段により撮像された画像内の各位置で合焦度が異なる複数の画像から第2距離画像を生成する第2距離画像生成ステップと、
    前記被写体像に示される複数の被写体から選択された主要被写体の被写界深度に応じて合成係数を算出する算出ステップと、
    前記第1距離画像生成ステップにより生成された第1距離画像と前記第2距離画像生成ステップにより生成された第2距離画像とを前記算出ステップにより算出された合成係数を用いて合成することで合成距離画像を生成する生成ステップと
    を備えたことを特徴とする制御方法。
  10. 被写体像を結像する光学結像手段と、前記光学結像手段により結像された被写体像を撮像する撮像手段とを備えた画像処理装置の制御方法をコンピュータに実行させるためのプログラムであって、
    前記制御方法は、
    前記撮像手段により撮像された異なる視差を有する複数の画像から第1距離画像を生成する第1距離画像生成ステップと、
    前記撮像手段により撮像された画像内の各位置で合焦度が異なる複数の画像から第2距離画像を生成する第2距離画像生成ステップと、
    前記被写体像に示される複数の被写体から選択された主要被写体の被写界深度に応じて合成係数を算出する算出ステップと、
    前記第1距離画像生成ステップにより生成された第1距離画像と前記第2距離画像生成ステップにより生成された第2距離画像とを前記算出ステップにより算出された合成係数を用いて合成することで合成距離画像を生成する生成ステップと
    を備えたことを特徴とするプログラム。
JP2012129865A 2012-06-07 2012-06-07 画像処理装置及びその制御方法、並びにプログラム Expired - Fee Related JP6071257B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2012129865A JP6071257B2 (ja) 2012-06-07 2012-06-07 画像処理装置及びその制御方法、並びにプログラム
US13/910,644 US10205928B2 (en) 2012-06-07 2013-06-05 Image processing apparatus that estimates distance information, method of controlling the same, and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2012129865A JP6071257B2 (ja) 2012-06-07 2012-06-07 画像処理装置及びその制御方法、並びにプログラム

Publications (3)

Publication Number Publication Date
JP2013254097A true JP2013254097A (ja) 2013-12-19
JP2013254097A5 JP2013254097A5 (ja) 2015-07-23
JP6071257B2 JP6071257B2 (ja) 2017-02-01

Family

ID=49714989

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012129865A Expired - Fee Related JP6071257B2 (ja) 2012-06-07 2012-06-07 画像処理装置及びその制御方法、並びにプログラム

Country Status (2)

Country Link
US (1) US10205928B2 (ja)
JP (1) JP6071257B2 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016197857A (ja) * 2015-04-03 2016-11-24 コグネックス・コーポレーション ホモグラフィの修正
US10275863B2 (en) 2015-04-03 2019-04-30 Cognex Corporation Homography rectification

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2513112B (en) * 2013-04-08 2020-01-08 Snell Advanced Media Ltd Video sequence processing
KR20150074348A (ko) * 2013-12-24 2015-07-02 삼성전자주식회사 디스플레이 장치 및 이를 이용한 영상표시방법
JP2015156938A (ja) * 2014-02-24 2015-09-03 ソニー株式会社 画像処理装置、画像処理方法
KR101932547B1 (ko) * 2014-10-23 2018-12-27 한화테크윈 주식회사 카메라 시스템 및 이의 영상 정합 방법
JP6723079B2 (ja) * 2016-06-08 2020-07-15 日立オートモティブシステムズ株式会社 物体距離検出装置
JP6636963B2 (ja) * 2017-01-13 2020-01-29 株式会社東芝 画像処理装置及び画像処理方法
JP7431527B2 (ja) * 2019-08-07 2024-02-15 キヤノン株式会社 深度情報生成装置、撮像装置、深度情報生成方法、画像処理装置、画像処理方法、及びプログラム

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07181024A (ja) * 1993-12-24 1995-07-18 Canon Inc 3次元形状計測方法および装置
JP2008116308A (ja) * 2006-11-02 2008-05-22 Fujifilm Corp 距離画像生成方法及びその装置
JP2010016743A (ja) * 2008-07-07 2010-01-21 Olympus Corp 測距装置、測距方法、測距プログラム又は撮像装置
US20110298898A1 (en) * 2010-05-11 2011-12-08 Samsung Electronics Co., Ltd. Three dimensional image generating system and method accomodating multi-view imaging
WO2011158508A1 (ja) * 2010-06-17 2011-12-22 パナソニック株式会社 画像処理装置および画像処理方法
WO2011158515A1 (ja) * 2010-06-17 2011-12-22 パナソニック株式会社 距離推定装置、距離推定方法、集積回路、コンピュータプログラム
JP2013509104A (ja) * 2009-10-23 2013-03-07 クゥアルコム・インコーポレイテッド 2dビデオデータの3dビデオデータへの変換のための深度マップ生成技法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2594454B2 (ja) * 1988-03-16 1997-03-26 旭光学工業株式会社 カメラの被写界深度の表示方法
JPH0674762A (ja) 1992-08-31 1994-03-18 Olympus Optical Co Ltd 距離測定装置
US6148120A (en) * 1997-10-30 2000-11-14 Cognex Corporation Warping of focal images to correct correspondence error
EP1237363B1 (en) * 1999-11-22 2006-01-18 Matsushita Electric Industrial Co., Ltd. Solid-state imaging device
JP4524717B2 (ja) * 2008-06-13 2010-08-18 富士フイルム株式会社 画像処理装置、撮像装置、画像処理方法及びプログラム
JP5369729B2 (ja) 2009-02-03 2013-12-18 株式会社ニコン 画像処理装置、撮像装置およびプログラム
US20130090415A1 (en) * 2010-06-18 2013-04-11 Yukio Kusano Additive for reclamation of asphalt, reclaimed asphalt pavement material containing same, modified asphalt, and asphalt pavement material containing same
US9270970B2 (en) * 2010-10-27 2016-02-23 Dolby International Ab Device apparatus and method for 3D image interpolation based on a degree of similarity between a motion vector and a range motion vector
US8937644B2 (en) * 2012-03-21 2015-01-20 Canon Kabushiki Kaisha Stereoscopic image capture
JP6312487B2 (ja) * 2014-03-26 2018-04-18 キヤノン株式会社 画像処理装置及びその制御方法、並びに、プログラム

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07181024A (ja) * 1993-12-24 1995-07-18 Canon Inc 3次元形状計測方法および装置
JP2008116308A (ja) * 2006-11-02 2008-05-22 Fujifilm Corp 距離画像生成方法及びその装置
JP2010016743A (ja) * 2008-07-07 2010-01-21 Olympus Corp 測距装置、測距方法、測距プログラム又は撮像装置
JP2013509104A (ja) * 2009-10-23 2013-03-07 クゥアルコム・インコーポレイテッド 2dビデオデータの3dビデオデータへの変換のための深度マップ生成技法
US20110298898A1 (en) * 2010-05-11 2011-12-08 Samsung Electronics Co., Ltd. Three dimensional image generating system and method accomodating multi-view imaging
WO2011158508A1 (ja) * 2010-06-17 2011-12-22 パナソニック株式会社 画像処理装置および画像処理方法
WO2011158515A1 (ja) * 2010-06-17 2011-12-22 パナソニック株式会社 距離推定装置、距離推定方法、集積回路、コンピュータプログラム

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016197857A (ja) * 2015-04-03 2016-11-24 コグネックス・コーポレーション ホモグラフィの修正
US10275863B2 (en) 2015-04-03 2019-04-30 Cognex Corporation Homography rectification

Also Published As

Publication number Publication date
US10205928B2 (en) 2019-02-12
JP6071257B2 (ja) 2017-02-01
US20130329019A1 (en) 2013-12-12

Similar Documents

Publication Publication Date Title
JP6071257B2 (ja) 画像処理装置及びその制御方法、並びにプログラム
JP6112824B2 (ja) 画像処理方法および装置、プログラム。
JP6245885B2 (ja) 撮像装置およびその制御方法
JP5968107B2 (ja) 画像処理方法、画像処理装置およびプログラム
CN105453136B (zh) 使用自动聚焦反馈进行立体侧倾校正的系统、方法及设备
JP4216824B2 (ja) 3次元モデル生成装置、3次元モデル生成方法および3次元モデル生成プログラム
US8144974B2 (en) Image processing apparatus, method, and program
CN110689581A (zh) 结构光模组标定方法、电子设备、计算机可读存储介质
JPWO2011158498A1 (ja) 撮像装置及び撮像方法
JP2011232330A (ja) 撮像装置、長さ測定方法、及びプログラム
JP6452360B2 (ja) 画像処理装置、撮像装置、画像処理方法およびプログラム
JP2016129289A (ja) 画像処理装置、撮像装置、画像処理方法、プログラム、および、記憶媒体
JPWO2011125937A1 (ja) キャリブレーションデータ選択装置、選択方法、及び選択プログラム、並びに三次元位置測定装置
JP2010117593A (ja) 距離情報取得装置、撮像装置、及びプログラム
JP7378219B2 (ja) 撮像装置、画像処理装置、制御方法、及びプログラム
JPWO2012176556A1 (ja) 対応点探索装置、および距離測定装置
JP2022024688A (ja) デプスマップ生成装置及びそのプログラム、並びに、デプスマップ生成システム
JP2013044597A (ja) 画像処理装置および方法、プログラム
JP6732440B2 (ja) 画像処理装置、画像処理方法、及びそのプログラム
JP2016213578A (ja) 画像処理装置、撮像装置、画像処理方法、プログラム
JP5925109B2 (ja) 画像処理装置、その制御方法、および制御プログラム
JP2013120435A (ja) 画像処理装置および画像処理方法、プログラム
JP2016156702A (ja) 撮像装置および撮像方法
JP6755737B2 (ja) 距離測定装置、撮像装置、および距離測定方法
JP6292785B2 (ja) 画像処理装置、画像処理方法およびプログラム

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150604

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20150604

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20160316

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20160322

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160520

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20160809

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20161011

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20161227

R151 Written notification of patent or utility model registration

Ref document number: 6071257

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

LAPS Cancellation because of no payment of annual fees