(57)【特許請求の範囲】
【請求項１】
第１画像パターンと第２画像パターンとの対応関係をピクセルピッチの単位で求め、前記第１画像パターンにおける特定位置に対する前記第２画像パターンにおける対応位置の相対変位量（ａｘ，ａｙ）を求める第１変位量演算部と、
前記第１画像パターンを表す第１画像データのうち前記特定位置を中心とする画像データｐに対して離散的フーリエ変換を施して、その離散的フーリエ変換後のデータＰを求める第１フーリエ変換部と、
前記第２画像パターンを表す第２画像データのうち前記対応位置を中心とする画像データｑに対して離散的フーリエ変換を施して、その離散的フーリエ変換後のデータＱを求める第２フーリエ変換部と、
前記データＰと前記データＱの複素共役との積を表すデータＲ１を求める複素積演算部と、
ピクセルピッチ未満の相対変位量に相当する傾斜位相を前記データＲ１に付与して、その付与後のデータＲ２を求める傾斜位相付与部と、
前記データＲ２に対して離散的逆フーリエ変換または離散的フーリエ変換を施して、その変換後のデータｒ２を求める変換部と、
前記データｒ２のピーク位置を中心とする前記データｒ２の点対称の度合を求める点対称度合演算部と、
前記傾斜位相付与部において付与される傾斜位相を変化させながら前記傾斜位相付与部，前記変換部および前記点対称度合演算部の各処理を繰り返させて、前記点対称度合演算部において求められる前記データｒ２の点対称の度合が最も大きくなる傾斜位相を求め、その傾斜位相に基づいて前記画像データｐに対する前記画像データｑの相対変位量（ｓｘ，ｓｙ）を求めて、前記第１画像パターンにおける特定位置に対する前記第２画像パターンにおける対応位置の相対変位量（ａｘ＋ｓｘ，ａｙ＋ｓｙ）を求める第２変位量演算部と、
を備えることを特徴とする画像パターン照合装置。
【請求項２】
第１画像パターンと第２画像パターンとの対応関係をピクセルピッチの単位で求め、前記第１画像パターンにおける特定位置に対する前記第２画像パターンにおける対応位置の相対変位量（ａｘ，ａｙ）を求める第１変位量演算ステップと、
前記第１画像パターンを表す第１画像データのうち前記特定位置を中心とする画像データｐに対して離散的フーリエ変換を施して、その離散的フーリエ変換後のデータＰを求める第１フーリエ変換ステップと、
前記第２画像パターンを表す第２画像データのうち前記対応位置を中心とする画像データｑに対して離散的フーリエ変換を施して、その離散的フーリエ変換後のデータＱを求める第２フーリエ変換ステップと、
前記データＰと前記データＱの複素共役との積を表すデータＲ１を求める複素積演算ステップと、
ピクセルピッチ未満の相対変位量に相当する傾斜位相を前記データＲ１に付与して、その付与後のデータＲ２を求める傾斜位相付与ステップと、
前記データＲ２に対して離散的逆フーリエ変換または離散的フーリエ変換を施して、その変換後のデータｒ２を求める変換ステップと、
前記データｒ２のピーク位置を中心とする前記データｒ２の点対称の度合を求める点対称度合演算ステップと、
を備え、
前記傾斜位相付与ステップにおいて付与される傾斜位相を変化させながら前記傾斜位相付与ステップ，前記変換ステップおよび前記点対称度合演算ステップの各処理を繰り返させて、前記点対称度合演算ステップにおいて求められる前記データｒ２の点対称の度合が最も大きくなる傾斜位相を求め、その傾斜位相に基づいて前記画像データｐに対する前記画像データｑの相対変位量（ｓｘ，ｓｙ）を求めて、前記第１画像パターンにおける特定位置に対する前記第２画像パターンにおける対応位置の相対変位量（ａｘ＋ｓｘ，ａｙ＋ｓｙ）を求める、
ことを特徴とする画像パターン照合方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像パターン照合装置および方法に関するものである。
【背景技術】
【０００２】
第１画像パターンと第２画像パターンとを互いに照合する画像パターン照合は、画像の解析・認識などの様々な分野で重要な基本的な処理である。例えば、被写体を第１カメラにより撮像することで第１画像パターンを得るとともに、該被写体を第２カメラにより撮像することで第２画像パターンを得て、これら第１画像パターンと第２画像パターンとを互いに照合することで、被写体の三次元形状を求めることができる。また、例えば、互いに異なる時刻に共通のカメラにより被写体を撮像することで第１画像パターンおよび第２画像パターンを得て、これら第１画像パターンと第２画像パターンとを互いに照合することで、被写体の動きを測定することができる。このような画像パターン照合技術として幾つかの技術が知られている。
【０００３】
第１の画像パターン照合技術では、第１画像パターンを表す第１画像データ（または、その一部）をｘ方向へａｘピクセル分だけ移動させるとともにｙ方向へａｙピクセル分だけ移動させたものと、第２画像パターンを表す第２画像データ（または、その一部）とを対比して、両者の類似度を求める。様々な移動量（ａｘ，ａｙ）に対する類似度を求めて、二次元平面（ａｘ，ａｙ）における類似度分布を作成する。そして、この類似度分布において類似度がピークとなる移動量（ａｘ，ａｙ）を、第１画像パターンと第２画像パターンとの間の相対変位量を表すものとして求める。
【０００４】
この第１の画像パターン照合技術において、第１画像パターンと第２画像パターンとの間の相対変位量をピクセルピッチの単位で求める方法だけでなく、ピクセルピッチ未満の分解能で求める方法も知られている（特許文献１参照）。後者の場合、二次元平面（ａｘ，ａｙ）における類似度分布が所定の形状を有するものと仮定して、その仮定の下に相対変位量をピクセルピッチ未満の分解能で求める。
【０００５】
第２の画像パターン照合技術（特許文献２〜４および非特許文献１を参照）では、第１画像パターンを表す第１画像データ（または、その一部）ｐを離散的フーリエ変換してデータＰ（＝Ｆ[ｐ]）を得るとともに、第２画像パターンを表す第２画像データ（または、その一部）ｑを離散的フーリエ変換してデータＱ（＝Ｆ[ｑ]）を得る。さらに、データＰおよびデータＱそれぞれの振幅を一定値とした上で、データＰとデータＱの複素共役Ｑ＊との積を表すデータＲ（＝ＰＱ＊）を求め、このデータＲを離散的逆フーリエ変換して相関データｒ（＝Ｆ-1[Ｒ]）を求める。そして、この相関データｒにおける中心位置からピーク位置への変位量を、第１画像パターンと第２画像パターンとの間の相対変位量を表すものとして求める。Ｆは離散的フーリエ変換を表す。Ｆ-1は離散的逆フーリエ変換を表す。
【０００６】
この第２の画像パターン照合技術においても、第１画像パターンと第２画像パターンとの間の相対変位量をピクセルピッチの単位で求める方法だけでなく、ピクセルピッチ未満の分解能で求める方法も知られている。後者の場合、相関データｒの分布が所定の形状を有するものと仮定して、その仮定の下に相対変位量をピクセルピッチ未満の分解能で求める。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００８−１１７４１６号公報
【特許文献２】特開平５−１５９０５６号公報
【特許文献３】特開平１０−１３２５３４号公報
【特許文献４】特許第３９３００６７号公報
【非特許文献】
【０００８】
【非特許文献１】青木孝文、他、「位相限定相関法に基づく高精度マシンビジョン − ピクセル分解能の壁を越える画像センシング技術を目指して −」、Fundamentals Review, Vol.1, No.1, pp.30-40(2007).
【発明の概要】
【発明が解決しようとする課題】
【０００９】
しかしながら、上記従来の第１および第２の画像パターン照合技術は、類似度分布または相関データが或る形状を有するとの仮定の下に相対変位量をピクセルピッチ未満の分解能で求めるものであることから、その仮定が妥当でない場合には相対変位量を正確に求めることができない。
【００１０】
また、上記従来の第２の画像パターン照合技術は、フーリエ変換データＰ，Ｑそれぞれの振幅を一定値とした上で相関データｒを求めるものであることから、本来は振幅が少ないため相関データのピーク位置に影響を与えることの少ない位相値について振幅を一定値に持ち上げることになって、相関データｒのピーク位置に影響を与える。特に、フーリエ変換の面積が少ない場合、それらの位相のランダム性による打ち消しの効果が少なくなるので、相関データのピーク位置は曖昧となる。
【００１１】
本発明は、上記問題点を解消する為になされたものであり、第１画像パターンと第２画像パターンとの間の相対変位量をピクセルピッチ未満の分解能で正確に求めることができる画像パターン照合装置および方法を提供することを目的とする。
【課題を解決するための手段】
【００１２】
本発明に係る画像パターン照合装置は、(1) 第１画像パターンと第２画像パターンとの対応関係をピクセルピッチの単位で求め、第１画像パターンにおける特定位置に対する第２画像パターンにおける対応位置の相対変位量（ａｘ，ａｙ）を求める第１変位量演算部と、(2)第１画像パターンを表す第１画像データのうち特定位置を中心とする画像データｐに対して離散的フーリエ変換を施して、その離散的フーリエ変換後のデータＰを求める第１フーリエ変換部と、(3)第２画像パターンを表す第２画像データのうち対応位置を中心とする画像データｑに対して離散的フーリエ変換を施して、その離散的フーリエ変換後のデータＱを求める第２フーリエ変換部と、(4)データＰとデータＱの複素共役との積を表すデータＲ１を求める複素積演算部と、(5) ピクセルピッチ未満の相対変位量に相当する傾斜位相をデータＲ１に付与して、その付与後のデータＲ２を求める傾斜位相付与部と、(6)データＲ２に対して離散的逆フーリエ変換または離散的フーリエ変換を施して、その変換後のデータｒ２を求める変換部と、(7)データｒ２のピーク位置を中心とするデータｒ２の点対称の度合を求める点対称度合演算部と、(8) 傾斜位相付与部において付与される傾斜位相を変化させながら傾斜位相付与部，変換部および点対称度合演算部の各処理を繰り返させて、点対称度合演算部において求められるデータｒ２の点対称の度合が最も大きくなる傾斜位相を求め、その傾斜位相に基づいて画像データｐに対する画像データｑの相対変位量（ｓｘ，ｓｙ）を求めて、第１画像パターンにおける特定位置に対する第２画像パターンにおける対応位置の相対変位量（ａｘ＋ｓｘ，ａｙ＋ｓｙ）を求める第２変位量演算部と、を備えることを特徴とする。
【００１３】
本発明に係る画像パターン照合方法は、(1) 第１画像パターンと第２画像パターンとの対応関係をピクセルピッチの単位で求め、第１画像パターンにおける特定位置に対する第２画像パターンにおける対応位置の相対変位量（ａｘ，ａｙ）を求める第１変位量演算ステップと、(2)第１画像パターンを表す第１画像データのうち特定位置を中心とする画像データｐに対して離散的フーリエ変換を施して、その離散的フーリエ変換後のデータＰを求める第１フーリエ変換ステップと、(3)第２画像パターンを表す第２画像データのうち対応位置を中心とする画像データｑに対して離散的フーリエ変換を施して、その離散的フーリエ変換後のデータＱを求める第２フーリエ変換ステップと、(4)データＰとデータＱの複素共役との積を表すデータＲ１を求める複素積演算ステップと、(5)ピクセルピッチ未満の相対変位量に相当する傾斜位相をデータＲ１に付与して、その付与後のデータＲ２を求める傾斜位相付与ステップと、(6)データＲ２に対して離散的逆フーリエ変換または離散的フーリエ変換を施して、その変換後のデータｒ２を求める変換ステップと、(7)データｒ２のピーク位置を中心とするデータｒ２の点対称の度合を求める点対称度合演算ステップと、を備えることを特徴とする。さらに、本発明に係る画像パターン照合方法は、傾斜位相付与ステップにおいて付与される傾斜位相を変化させながら傾斜位相付与ステップ，変換ステップおよび点対称度合演算ステップの各処理を繰り返させて、点対称度合演算ステップにおいて求められるデータｒ２の点対称の度合が最も大きくなる傾斜位相を求め、その傾斜位相に基づいて画像データｐに対する画像データｑの相対変位量（ｓｘ，ｓｙ）を求めて、第１画像パターンにおける特定位置に対する第２画像パターンにおける対応位置の相対変位量（ａｘ＋ｓｘ，ａｙ＋ｓｙ）を求めることを特徴とする。
【発明の効果】
【００１４】
本発明によれば、第１画像パターンと第２画像パターンとの間の相対変位量をピクセルピッチ未満の分解能で正確に求めることができる。
【図面の簡単な説明】
【００１５】
【図１】本実施形態に係る画像パターン照合装置１の構成を示す図である。
【図２】本実施形態に係る画像パターン照合方法を示すフローチャートである。
【図３】本実施形態に係る画像パターン照合装置１の第１変位量演算部１１または本実施形態に係る画像パターン照合方法の第１変位量演算ステップＳ１１における処理を説明する図である。
【図４】本実施形態に係る画像パターン照合方法の各ステップで求められるデータを示す図である。
【図５】比較例の画像パターン照合方法による実験結果を示す図である。
【図６】本実施形態に係る画像パターン照合方法による実験結果を示す図である。
【発明を実施するための形態】
【００１６】
以下、添付図面を参照して、本発明を実施するための形態を詳細に説明する。なお、図面の説明において同一の要素には同一の符号を付し、重複する説明を省略する。
【００１７】
図１は、本実施形態に係る画像パターン照合装置１の構成を示す図である。本実施形態に係る画像パターン照合装置１は、第１変位量演算部１１、第１フーリエ変換部１２、第２フーリエ変換部１３、複素積演算部１４、傾斜位相付与部１５、変換部１６、点対称度合演算部１７および第２変位量演算部１８を備える。なお、第１フーリエ変換部１２と第２フーリエ変換部１３とは共通であってもよい。
【００１８】
第１変位量演算部１１は、第１画像パターンと第２画像パターンとを入力し、両者の対応関係をピクセルピッチの単位で求めて、第１画像パターンにおける特定位置に対する第２画像パターンにおける対応位置の相対変位量（ａｘ，ａｙ）を求める。この第１変位量演算部１１においてピクセルピッチの単位で相対変位量（ａｘ，ａｙ）を求めるには、第１画像パターンを表す第１画像データ（または、その一部）をｘ方向へａｘピクセル分だけ移動させるとともにｙ方向へａｙピクセル分だけ移動させたものと、第２画像パターンを表す第２画像データ（または、その一部）とを対比して、両者の類似度を求める。様々な移動量（ａｘ，ａｙ）に対する類似度を求めて、二次元平面（ａｘ，ａｙ）における類似度分布を作成する。そして、この類似度分布において類似度がピークとなる移動量（ａｘ，ａｙ）を、第１画像パターンと第２画像パターンとの間の相対変位量を表すものとして求める。類似度としては、対応する画素毎の両画像データの差の絶対値和や相関係数が用いられる。
【００１９】
第１フーリエ変換部１２は、第１画像パターンを表す第１画像データのうち特定位置を中心とする画像データｐに対して離散的フーリエ変換を施して、その離散的フーリエ変換後のデータＰを求める。すなわち、第１フーリエ変換部１２は以下の演算を行う。Ｆは離散的フーリエ変換を表す。
【００２０】
Ｐ＝Ｆ[ｐ]       …(1)
【００２１】
第２フーリエ変換部１３は、第２画像パターンを表す第２画像データのうち特定位置を中心とする画像データｑに対して離散的フーリエ変換を施して、その離散的フーリエ変換後のデータＱを求める。すなわち、第２フーリエ変換部１３は以下の演算を行う。
【００２２】
Ｑ＝Ｆ[ｑ]       …(2)
【００２３】
複素積演算部１４は、第１フーリエ変換部１２により求められたデータＰと第２フーリエ変換部１３により求められたデータＱの複素共役Ｑ＊との積を表すデータＲ１を求める。すなわち、複素積演算部１４は以下の演算を行う。
【００２４】
Ｒ１＝ＰＱ＊      …(3)
【００２５】
傾斜位相付与部１５は、複素積演算部１４により求められたデータＲ１に対して、ピクセルピッチ未満の相対変位量（ｓｘ，ｓｙ）に相当する傾斜位相（φｘ，φｙ）を付与して、その付与後のデータＲ２を求める。すなわち、傾斜位相付与部１５は以下の演算を行う。Ｎｘ×ＮｙはデータＲ１のサイズである。ｋｘ，ｋｙはフーリエ空間での座標値である。
【００２６】
Ｒ２＝Ｒ１・exp(ｉ(φｘｋｘ＋φｙｋｙ))      …(4)
φｘ＝(２π／Ｎｘ)・ｓｘ                  …(5)
φｙ＝(２π／Ｎｙ)・ｓｙ                  …(6)
【００２７】
変換部１６は、傾斜位相付与部１５により求められたデータＲ２に対して離散的逆フーリエ変換を施して、その離散的逆フーリエ変換後のデータｒ２を求める。すなわち、変換部１６は以下の演算を行う。Ｆ-1は離散的逆フーリエ変換を表す。なお、傾斜位相量（φｘｋｘ＋φｙｋｙ）が０〜２πの範囲を超える場合には、傾斜位相量（φｘｋｘ＋φｙｋｙ）を２πで除算したときの剰余を用いる。
【００２８】
ｒ２＝Ｆ-1[Ｒ２]      …(7)
【００２９】
この変換部１６により求められたデータｒ２は、データＲ１を離散的逆フーリエ変換して得られるデータｒ１（＝Ｆ-1[Ｒ１]）対して相対変位量（ｓｘ，ｓｙ）だけ平行移動したものとなる。ｓｘ，ｓｙはピクセルピッチ未満の量である。
【００３０】
点対称度合演算部１７は、変換部１６により求められたデータｒ２のピーク位置を中心とするデータｒ２の点対称の度合を求める。点対称の度合を求めるに際しては、データｒ２の絶対値のピーク位置に隣接する位置であって該ピーク位置を中心とする点対称関係にある２つの位置それぞれのデータｒ２の絶対値の差を求めることで得られ、この差が小さいほど点対称の度合が大きいことを表す。
【００３１】
第２変位量演算部１８は、傾斜位相付与部１５において付与される傾斜位相を変化させながら傾斜位相付与部１５，変換部１６および点対称度合演算部１７の各処理を繰り返させる。そして、第２変位量演算部１８は、点対称度合演算部１７において求められるデータｒ２の点対称の度合が最も大きくなる傾斜位相を求め、その傾斜位相に基づいて画像データｐに対する画像データｑの相対変位量（ｓｘ，ｓｙ）を求めて、第１画像パターンにおける特定位置に対する第２画像パターンにおける対応位置の相対変位量（ａｘ＋ｓｘ，ａｙ＋ｓｙ）を求める。
【００３２】
なお、変換部１６は、上述したようにデータＲ２に対して離散的逆フーリエ変換を施してもよいし、データＲ２に対して離散的フーリエ変換を施してもよい。データＲ２に対する離散的逆フーリエ変換の結果と、データＲ２に対する離散的フーリエ変換の結果とは、画像中心に対して点対称の関係にあるから、何れの場合にも点対称度合演算部１７において変換部１６の結果から同様にして点対称度合を求めることができる。
【００３３】
次に、本実施形態に係る画像パターン照合装置１の動作について説明するとともに、本実施形態に係る画像パターン照合方法について説明する。図２は、本実施形態に係る画像パターン照合方法を示すフローチャートである。図３は、本実施形態に係る画像パターン照合装置１の第１変位量演算部１１または本実施形態に係る画像パターン照合方法の第１変位量演算ステップＳ１１における処理を説明する図である。また、図４は、本実施形態に係る画像パターン照合方法の各ステップで求められるデータを示す図である。
【００３４】
本実施形態に係る画像パターン照合方法は、図２に示されるように、第１変位量演算ステップＳ１１、第１フーリエ変換ステップＳ１２、第２フーリエ変換ステップＳ１３、複素積演算ステップＳ１４、傾斜位相付与ステップＳ１５、変換ステップＳ１６、点対称度合演算ステップＳ１７、最大点対称度合更新ステップＳ１８および終了判断ステップＳ１９を備える。
【００３５】
第１変位量演算ステップＳ１１では、第１変位量演算部１１により、図３に示されるように、第１画像パターン２１と第２画像パターン２２との対応関係がピクセルピッチの単位で求められて、第１画像パターン２１における特定位置３１に対する第２画像パターン２２における対応位置３２の相対変位量（ａｘ，ａｙ）が求められる。
【００３６】
第１フーリエ変換ステップＳ１２では、第１フーリエ変換部１２により、第１画像パターン２１を表す第１画像データのうち特定位置３１を中心とする所定範囲４１の画像データｐに対して離散的フーリエ変換が施されて、その離散的フーリエ変換後のデータＰが求められる（上記(1)式）。なお、第１画像データが強度データ（図４（ａ））である場合には、その強度データの平方根の値が画像データｐの振幅とされ（図４（ｂ））、画像データｐの位相については値０とされる（図４（ｃ））。そして、この画像データｐの振幅および位相が用いられて離散的フーリエ変換が行われて、その離散的フーリエ変換後のデータＰが求められる。
【００３７】
第２フーリエ変換ステップＳ１３では、第２フーリエ変換部１３により、第２画像パターン２２を表す第２画像データのうち対応位置３２を中心とする所定範囲４２の画像データｑに対して離散的フーリエ変換が施されて、その離散的フーリエ変換後のデータＱが求められる（上記(2)式）。なお、第２画像データが強度データ（図４（ｄ））である場合には、その強度データの平方根の値が画像データｑの振幅とされ（図４（ｅ））、画像データｑの位相については値０とされる（図４（ｆ））。そして、この画像データｑの振幅および位相が用いられて離散的フーリエ変換が行われて、その離散的フーリエ変換後のデータＱが求められる。
【００３８】
複素積演算ステップＳ１４では、複素積演算部１４により、第１フーリエ変換ステップＳ１２で求められたデータＰと第２フーリエ変換ステップＳ１３で求められたデータＱの複素共役Ｑ＊との積を表すデータＲ１の振幅（図４（ｇ））および位相（図４（ｈ））が求められる（上記(3)式）。
【００３９】
傾斜位相付与ステップＳ１５では、傾斜位相付与部１５により、複素積演算ステップＳ１４で求められたデータＲ１に対して、ピクセルピッチ未満の相対変位量（ｓｘ，ｓｙ）に相当する傾斜位相（φｘ，φｙ）が付与されて、その付与後のデータＲ２が求められる（上記(4)〜(6)式）。すなわち、データＲ１の位相（図４（ｈ））に傾斜位相（図４（ｉ））が加えられて、データＲ２の位相が得られる。データＲ２の振幅はデータＲ１の振幅と同じである。
【００４０】
変換ステップＳ１６では、変換部１６により、傾斜位相付与ステップＳ１５で求められたデータＲ２に対して離散的逆フーリエ変換が施されて、その離散的逆フーリエ変換後のデータｒ２の振幅（図４（ｊ））および位相（図４（ｋ））が求められる（上記(7)式）。なお、変換ステップＳ１６では、データＲ２に対して離散的フーリエ変換が施されもよい。
【００４１】
点対称度合演算ステップＳ１７では、点対称度合演算部１７により、変換ステップＳ１６で求められたデータｒ２のピーク位置を中心とするデータｒ２の点対称の度合が求められる。
【００４２】
最大点対称度合更新ステップＳ１８では、直前の点対称度合演算ステップＳ１７で求められた点対称度合と、これまで記憶していた点対称度合とのうち、大きい値の点対称度合が更新記憶される。すなわち、これまでに求められた幾つかの点対称度合のうち最大値が常に記憶される。
【００４３】
終了判断ステップＳ１９では、傾斜位相付与ステップＳ１５，変換ステップＳ１６，点対称度合演算ステップＳ１７および最大点対称度合更新ステップＳ１８を含むループ処理を終了してよいか否かが判断されて、さらにループ処理を続ける必要があると判断されたときに傾斜位相付与ステップＳ１５へ戻る。このループ処理は、傾斜位相付与ステップＳ１５において付与される傾斜位相（φｘ，φｙ）を変化させながら行われる。すなわち、傾斜位相付与ステップＳ１５において付与される傾斜位相を変化させながら傾斜位相付与ステップＳ１５，変換ステップＳ１６，点対称度合演算ステップＳ１７および最大点対称度合更新ステップＳ１８の各処理が繰り返される。
【００４４】
終了判断ステップＳ１９におけるループ処理終了可否の判断は、傾斜位相（φｘ，φｙ）について、すなわち、相対変位量（ｓｘ，ｓｙ）について、変化させるべき範囲の全体についてループ処理が終了したか否かに基づいて為され、或いは、所望のサブピクセル精度が得られたか否かに基づいて為される。
【００４５】
以上の処理が終了すると、第２変位量演算部１８により、最終的に得られた最大点対称度合に対応する傾斜位相（φｘ，φｙ）に基づいて、画像データｐに対する画像データｑの相対変位量（ｓｘ，ｓｙ）が求められて、第１画像パターン２１における特定位置３１に対する第２画像パターン２２における対応位置３２の相対変位量（ａｘ＋ｓｘ，ａｙ＋ｓｙ）が求められる。
【００４６】
本実施形態に係る画像パターン照合装置または画像パターン照合方法は、類似度分布または相関データが或る形状を有するとの仮定を用いる必要がないので、第１画像パターンと第２画像パターンとの間の相対変位量をピクセルピッチ未満の分解能で正確に求めることができる。
【００４７】
本実施形態に係る画像パターン照合装置または画像パターン照合方法において、第１画像データが第１カメラにより被写体が撮像されて得られたものであって、第２画像データが第２カメラにより被写体が撮像されて得られたものである場合、相対変位量（ａｘ＋ｓｘ，ａｙ＋ｓｙ）に基づいて特定位置に関する第１カメラと第２カメラとの視差を求めて、この視差に基づいて特定位置に対応する被写体上の位置までの距離を求めることができる。また、複数の特定位置それぞれについて距離を求めることができ、さらに、被写体の三次元形状を求めることができる。
【００４８】
また、本実施形態に係る画像パターン照合装置または画像パターン照合方法において、第１画像データおよび第２画像データが互いに異なる時刻に共通のカメラにより被写体が撮像されて得られたものである場合、相対変位量（ａｘ＋ｓｘ，ａｙ＋ｓｙ）に基づいて特定位置の移動方向および移動量を求めることができ、また、複数の特定位置それぞれについて移動方向および移動量を求めることができる。
【００４９】
図５は、比較例の画像パターン照合方法による実験結果を示す図である。また、図６は、本実施形態に係る画像パターン照合方法による実験結果を示す図である。所定パターンをＸ軸方向およびＹ軸方向それぞれに単位移動量ずつステージにより移動させて、基準位置にある所定パターンを撮像して得られた画像パターンと、各位置にある所定パターンを撮像して得られた画像パターンとを用いて、ピーク位置の移動量を求めた。画像データｐおよび画像データｑの各サイズを３２×３２ピクセルとした。比較例は上記の第２の画像パターン照合技術に拠る。両図から判るように、比較例と対比して本実施形態ではピーク位置がピクセルピッチ未満の分解能で高精度に得られている。なお、画像データｐおよび画像データｑの各サイズを６４×６４ピクセルおよび１２８×１２８ピクセルそれぞれとした場合にも、比較例と対比して本実施形態ではピーク位置がピクセルピッチ未満の分解能で高精度に得られた。
【符号の説明】
【００５０】
１…画像パターン照合装置、１１…第１変位量演算部、１２…第１フーリエ変換部、１３…第２フーリエ変換部、１４…複素積演算部、１５…傾斜位相付与部、１６…変換部、１７…点対称度合演算部、１８…第２変位量演算部。
【図１】
0005491162.tif 000002
【図２】
0005491162.tif 000003
【図３】
0005491162.tif 000004
【図５】
0005491162.tif 000005
【図６】
0005491162.tif 000006
【図４】
0005491162.tif 000007