JP5491162B2 - 画像パターン照合装置および方法 - Google Patents

Description

本発明は、画像パターン照合装置および方法に関するものである。

第１画像パターンと第２画像パターンとを互いに照合する画像パターン照合は、画像の解析・認識などの様々な分野で重要な基本的な処理である。例えば、被写体を第１カメラにより撮像することで第１画像パターンを得るとともに、該被写体を第２カメラにより撮像することで第２画像パターンを得て、これら第１画像パターンと第２画像パターンとを互いに照合することで、被写体の三次元形状を求めることができる。また、例えば、互いに異なる時刻に共通のカメラにより被写体を撮像することで第１画像パターンおよび第２画像パターンを得て、これら第１画像パターンと第２画像パターンとを互いに照合することで、被写体の動きを測定することができる。このような画像パターン照合技術として幾つかの技術が知られている。

第１の画像パターン照合技術では、第１画像パターンを表す第１画像データ（または、その一部）をｘ方向へａ_ｘピクセル分だけ移動させるとともにｙ方向へａ_ｙピクセル分だけ移動させたものと、第２画像パターンを表す第２画像データ（または、その一部）とを対比して、両者の類似度を求める。様々な移動量（ａ_ｘ，ａ_ｙ）に対する類似度を求めて、二次元平面（ａ_ｘ，ａ_ｙ）における類似度分布を作成する。そして、この類似度分布において類似度がピークとなる移動量（ａ_ｘ，ａ_ｙ）を、第１画像パターンと第２画像パターンとの間の相対変位量を表すものとして求める。

この第１の画像パターン照合技術において、第１画像パターンと第２画像パターンとの間の相対変位量をピクセルピッチの単位で求める方法だけでなく、ピクセルピッチ未満の分解能で求める方法も知られている（特許文献１参照）。後者の場合、二次元平面（ａ_ｘ，ａ_ｙ）における類似度分布が所定の形状を有するものと仮定して、その仮定の下に相対変位量をピクセルピッチ未満の分解能で求める。

第２の画像パターン照合技術（特許文献２〜４および非特許文献１を参照）では、第１画像パターンを表す第１画像データ（または、その一部）ｐを離散的フーリエ変換してデータＰ（＝Ｆ[ｐ]）を得るとともに、第２画像パターンを表す第２画像データ（または、その一部）ｑを離散的フーリエ変換してデータＱ（＝Ｆ[ｑ]）を得る。さらに、データＰおよびデータＱそれぞれの振幅を一定値とした上で、データＰとデータＱの複素共役Ｑ^＊との積を表すデータＲ（＝ＰＱ^＊）を求め、このデータＲを離散的逆フーリエ変換して相関データｒ（＝Ｆ^-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)式）。なお、変換ステップＳ１６では、データＲ_２に対して離散的フーリエ変換が施されもよい。

点対称度合演算ステップＳ１７では、点対称度合演算部１７により、変換ステップＳ１６で求められたデータｒ_２のピーク位置を中心とするデータｒ_２の点対称の度合が求められる。

最大点対称度合更新ステップＳ１８では、直前の点対称度合演算ステップＳ１７で求められた点対称度合と、これまで記憶していた点対称度合とのうち、大きい値の点対称度合が更新記憶される。すなわち、これまでに求められた幾つかの点対称度合のうち最大値が常に記憶される。

終了判断ステップＳ１９では、傾斜位相付与ステップＳ１５，変換ステップＳ１６，点対称度合演算ステップＳ１７および最大点対称度合更新ステップＳ１８を含むループ処理を終了してよいか否かが判断されて、さらにループ処理を続ける必要があると判断されたときに傾斜位相付与ステップＳ１５へ戻る。このループ処理は、傾斜位相付与ステップＳ１５において付与される傾斜位相（φ_ｘ，φ_ｙ）を変化させながら行われる。すなわち、傾斜位相付与ステップＳ１５において付与される傾斜位相を変化させながら傾斜位相付与ステップＳ１５，変換ステップＳ１６，点対称度合演算ステップＳ１７および最大点対称度合更新ステップＳ１８の各処理が繰り返される。

終了判断ステップＳ１９におけるループ処理終了可否の判断は、傾斜位相（φ_ｘ，φ_ｙ）について、すなわち、相対変位量（ｓ_ｘ，ｓ_ｙ）について、変化させるべき範囲の全体についてループ処理が終了したか否かに基づいて為され、或いは、所望のサブピクセル精度が得られたか否かに基づいて為される。

以上の処理が終了すると、第２変位量演算部１８により、最終的に得られた最大点対称度合に対応する傾斜位相（φ_ｘ，φ_ｙ）に基づいて、画像データｐに対する画像データｑの相対変位量（ｓ_ｘ，ｓ_ｙ）が求められて、第１画像パターン２１における特定位置３１に対する第２画像パターン２２における対応位置３２の相対変位量（ａ_ｘ＋ｓ_ｘ，ａ_ｙ＋ｓ_ｙ）が求められる。

本実施形態に係る画像パターン照合装置または画像パターン照合方法は、類似度分布または相関データが或る形状を有するとの仮定を用いる必要がないので、第１画像パターンと第２画像パターンとの間の相対変位量をピクセルピッチ未満の分解能で正確に求めることができる。

本実施形態に係る画像パターン照合装置または画像パターン照合方法において、第１画像データが第１カメラにより被写体が撮像されて得られたものであって、第２画像データが第２カメラにより被写体が撮像されて得られたものである場合、相対変位量（ａ_ｘ＋ｓ_ｘ，ａ_ｙ＋ｓ_ｙ）に基づいて特定位置に関する第１カメラと第２カメラとの視差を求めて、この視差に基づいて特定位置に対応する被写体上の位置までの距離を求めることができる。また、複数の特定位置それぞれについて距離を求めることができ、さらに、被写体の三次元形状を求めることができる。

また、本実施形態に係る画像パターン照合装置または画像パターン照合方法において、第１画像データおよび第２画像データが互いに異なる時刻に共通のカメラにより被写体が撮像されて得られたものである場合、相対変位量（ａ_ｘ＋ｓ_ｘ，ａ_ｙ＋ｓ_ｙ）に基づいて特定位置の移動方向および移動量を求めることができ、また、複数の特定位置それぞれについて移動方向および移動量を求めることができる。

図５は、比較例の画像パターン照合方法による実験結果を示す図である。また、図６は、本実施形態に係る画像パターン照合方法による実験結果を示す図である。所定パターンをＸ軸方向およびＹ軸方向それぞれに単位移動量ずつステージにより移動させて、基準位置にある所定パターンを撮像して得られた画像パターンと、各位置にある所定パターンを撮像して得られた画像パターンとを用いて、ピーク位置の移動量を求めた。画像データｐおよび画像データｑの各サイズを３２×３２ピクセルとした。比較例は上記の第２の画像パターン照合技術に拠る。両図から判るように、比較例と対比して本実施形態ではピーク位置がピクセルピッチ未満の分解能で高精度に得られている。なお、画像データｐおよび画像データｑの各サイズを６４×６４ピクセルおよび１２８×１２８ピクセルそれぞれとした場合にも、比較例と対比して本実施形態ではピーク位置がピクセルピッチ未満の分解能で高精度に得られた。

１…画像パターン照合装置、１１…第１変位量演算部、１２…第１フーリエ変換部、１３…第２フーリエ変換部、１４…複素積演算部、１５…傾斜位相付与部、１６…変換部、１７…点対称度合演算部、１８…第２変位量演算部。

Claims (2)

1. 第１画像パターンと第２画像パターンとの対応関係をピクセルピッチの単位で求め、前記第１画像パターンにおける特定位置に対する前記第２画像パターンにおける対応位置の相対変位量（ａ_ｘ，ａ_ｙ）を求める第１変位量演算部と、
   前記第１画像パターンを表す第１画像データのうち前記特定位置を中心とする画像データｐに対して離散的フーリエ変換を施して、その離散的フーリエ変換後のデータＰを求める第１フーリエ変換部と、
   前記第２画像パターンを表す第２画像データのうち前記対応位置を中心とする画像データｑに対して離散的フーリエ変換を施して、その離散的フーリエ変換後のデータＱを求める第２フーリエ変換部と、
   前記データＰと前記データＱの複素共役との積を表すデータＲ_１を求める複素積演算部と、
   ピクセルピッチ未満の相対変位量に相当する傾斜位相を前記データＲ_１に付与して、その付与後のデータＲ_２を求める傾斜位相付与部と、
   前記データＲ_２に対して離散的逆フーリエ変換または離散的フーリエ変換を施して、その変換後のデータｒ_２を求める変換部と、
   前記データｒ_２のピーク位置を中心とする前記データｒ_２の点対称の度合を求める点対称度合演算部と、
   前記傾斜位相付与部において付与される傾斜位相を変化させながら前記傾斜位相付与部，前記変換部および前記点対称度合演算部の各処理を繰り返させて、前記点対称度合演算部において求められる前記データｒ_２の点対称の度合が最も大きくなる傾斜位相を求め、その傾斜位相に基づいて前記画像データｐに対する前記画像データｑの相対変位量（ｓ_ｘ，ｓ_ｙ）を求めて、前記第１画像パターンにおける特定位置に対する前記第２画像パターンにおける対応位置の相対変位量（ａ_ｘ＋ｓ_ｘ，ａ_ｙ＋ｓ_ｙ）を求める第２変位量演算部と、
   を備えることを特徴とする画像パターン照合装置。
2. 第１画像パターンと第２画像パターンとの対応関係をピクセルピッチの単位で求め、前記第１画像パターンにおける特定位置に対する前記第２画像パターンにおける対応位置の相対変位量（ａ_ｘ，ａ_ｙ）を求める第１変位量演算ステップと、
   前記第１画像パターンを表す第１画像データのうち前記特定位置を中心とする画像データｐに対して離散的フーリエ変換を施して、その離散的フーリエ変換後のデータＰを求める第１フーリエ変換ステップと、
   前記第２画像パターンを表す第２画像データのうち前記対応位置を中心とする画像データｑに対して離散的フーリエ変換を施して、その離散的フーリエ変換後のデータＱを求める第２フーリエ変換ステップと、
   前記データＰと前記データＱの複素共役との積を表すデータＲ_１を求める複素積演算ステップと、
   ピクセルピッチ未満の相対変位量に相当する傾斜位相を前記データＲ_１に付与して、その付与後のデータＲ_２を求める傾斜位相付与ステップと、
   前記データＲ_２に対して離散的逆フーリエ変換または離散的フーリエ変換を施して、その変換後のデータｒ_２を求める変換ステップと、
   前記データｒ_２のピーク位置を中心とする前記データｒ_２の点対称の度合を求める点対称度合演算ステップと、
   を備え、
   前記傾斜位相付与ステップにおいて付与される傾斜位相を変化させながら前記傾斜位相付与ステップ，前記変換ステップおよび前記点対称度合演算ステップの各処理を繰り返させて、前記点対称度合演算ステップにおいて求められる前記データｒ_２の点対称の度合が最も大きくなる傾斜位相を求め、その傾斜位相に基づいて前記画像データｐに対する前記画像データｑの相対変位量（ｓ_ｘ，ｓ_ｙ）を求めて、前記第１画像パターンにおける特定位置に対する前記第２画像パターンにおける対応位置の相対変位量（ａ_ｘ＋ｓ_ｘ，ａ_ｙ＋ｓ_ｙ）を求める、
   ことを特徴とする画像パターン照合方法。