JP2014186659A - 画像照合装置、画像照合方法、画像照合プログラム - Google Patents

Abstract

【課題】 画像を高精度に照合する画像照合装置、画像照合方法、画像照合プログラムを提供すること。
【解決手段】 登録紙片の画像を登録する紙片登録部２１と、照合紙片の画像を登録紙片の画像と照合する紙片識別部２２とを備える。画像から特徴点を抽出し、それらの特徴点同士を用いて画像変換パラメータを算出して照合をする画像照合処理において、画像変換パラメータによって推定された画像の妥当性を評価することで、推定するパラメータを順次入れ替えながら画像の照合を行う。また、局所的な回転や伸縮に対応するため、一組の画像変換パラメータで対応づけられない対応点同士を別のパラメータで対応づける処理を再帰的に行うことで、画像全体を照合する。
【選択図】 図２

Description

本発明は、画像照合装置、画像照合方法、画像照合プログラムに係り、特に、紙片やシート等の物体の真贋を判定するための画像照合装置、画像照合方法、画像照合プログラムに関する。

近年、複写技術の向上により紙幣や有価証券、切手等の精巧な偽造券が、比較的簡便に作成可能となり、社会問題化する傾向にあろう。一般的に、紙片の偽造の防止には、漉き込み技術や特殊なインクが必要で、印刷コストを上げる原因となっている。例えば、特許文献１には、証券の文字を磁気インクで印刷し、複写機による複製を防止する方法が開示されている。特許文献２には、磁性体を漉き込んで真贋判定を行う方法が開示されている。
また、画像を用いた照合システムとして、画像から複数の特徴点を抽出し、抽出した特徴点から得られる特徴量や座標情報を用いたシステムが提案されている。例えば、特許文献３には、画像間の特徴量、距離、角度を算出して画像同士をマッチングする手法が開示されている。また、特許文献４には、画像から複数の特徴点を抽出し、特徴量との距離により画像マッチングを行う方法が開示されている。さらに、特許文献５には、特徴点とユークリッド距離に基づいて画像同士の類似度を求める方法が開示されている。また、特許文献６、７には、アフィン変換を用いて登録画像と照合画像の対応関係を求め、画像の回転や拡大・縮小に対して頑健な画像照合方式が開示されている。

特開２０１０−９４８２号公報 特開２００７−１６４２９３号公報 特開２００７−３３４７９５号公報 特開２０１２−１６００４７号公報 特開２０１２−１９０２２０号公報 特開２０１０−２７２０９１号公報 特開２００４−３２６６９３号公報 特開２００７−３３４７９５号公報 特開２００６−３３８３３０号公報

Ｄａｖｉｄ Ｇ． Ｌｏｗｅ、"Ｄｉｓｔｉｎｃｔｉｖｅ Ｉｍａｇｅ Ｆｅａｔｕｒｅｓ ｆｒｏｍ Ｓｃａｌｅ−Ｉｎｖａｒｉａｎｔ Ｋｅｙｐｏｉｎｔｓ"、Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｏｆ Ｃｏｍｐｕｔｅｒ Ｖｉｓｉｏｎ、 ６０、 ２、 ｐｐ． ９１−１１０、 ２００４．

特許文献１ないし２に記載の方法は、偽造防止に有効な方法である。しかしながら、特別なインクを用いたり、紙に特別な繊維を漉き込んだりする必要があり、偽造防止のコストが高くなるという課題があった。また、真贋判定をしたい対象物が画像として捉えられている場合、すなわち、印影等の場合には、適用できないという課題があった。
また、画像を用いて照合を行う際には、登録画像（以下、モデル画像とも呼ぶ）と、照合画像で、画像の撮影条件が異なる場合がある。例えば、光学スキャナを用いて紙片の画像を採取する際には、紙片の搬送系の精度の原因で、登録画像と照合画像が互いに数画素平行移動したり、回転したりすることがある。また、撮像センサの精度の原因で、登録画像に比べ、照合画像は若干拡大したり、縮小したりすることがある。
特許文献３に記載の方法は、上述の課題または原因を考慮した方法である。すなわち、モデル画像から複数の特徴点を抽出し、モデル画像をさまざまなカメラパラメータで変形した画像を生成する（以下、変調画像と呼ぶ）。そして、変調画像から特徴点を抽出し、モデル画像と変調画像の特徴点を対応付け、照合画像とモデル画像および全変調画像との特徴点を照合する方法である。しかしながら、照合した特徴点は、独立した変調方式で変調した複数の画像から抽出した特徴点となる。また、特許文献４に記載の方法は、画像から抽出した複数の特徴点を利用した画像照合方法であり、特徴点同士の対応関係を、ＲＡＮＳＡＣ（ＲＡＮｄｏｍ Ｓａｍｐｌｅ Ｃｏｎｓｅｎｓｕｓ）法を用いて変換行列を求めることで算出する方法である。しかしながら、特許文献３と同様に、画像から抽出した各々の特徴点位置関係に制約はなく、画像から抽出した各々の特徴点の特徴量が類似している場合、正しく照合できないことがあった。しがたって、各々の特徴点の位置関係に制約はなく、画像から抽出した各々の特徴点の特徴量が類似している場合、正しく照合できないという課題があった。

さらに、特許文献５に記載の方法は、非特許文献１に記載のＳＩＦＴ特徴量を用いて、画像の回転に頑健な画像照合方法である。ＳＩＦＴ特徴は、画像の回転や縮小に頑健な特徴量であるが、特許文献３や特許文献４と同様に、画像から抽出した各々の特徴点位置関係に制約はないため、画像から抽出した各々の特徴点の特徴量が類似している場合、正しく照合できないことがあった。
特許文献６に開示されている方法は、登録画像と照合画像から複数の特徴点を抽出し、照合画像中の特徴点と対応する特徴点を登録画像中に求め、それらの座標に基づいて、両画像同士の回転や拡大・縮小を表現するアフィン変換関数を求める方法である。したがって、最初に選んだ特徴点の対を用いて求めたアフィン変換関数が、他の特徴点群を適切に表現できない場合、照合画像が登録画像と同じ物体を撮影したものであっても、偽物と判断される場合がある。
特許文献７に開示されている方法も、特許文献６と同様にアフィン変換関数を用いて、登録画像と照合画像の関係を記述する、画像照合方式である。しかしながら、特徴点の対応付けは、ランダムに選択した３点の候補点対を繰り返し選択することで行う。したがって、適切な対応点が選択される保証はなく、また、特に画像から抽出した各々の特徴点の特徴量が類似している場合、正しく照合できないことがあった。
以上を整理すると、本発明にて解決する主な課題は以下のとおりとなる。すなわち、本発明が解決しようとする第一の課題は、登録画像の特徴点と照合画像の特徴点を関係づける画像変換関数のパラメータを高精度に算出することである。また、本発明が解決しようとする第二の課題は、登録画像に対して照合画像が、部分的に回転していたり、拡大・縮小していたりしても、高精度に画像を照合することである。

本発明は、以上の点に鑑み、画像を高精度に照合する画像照合装置、画像照合方法、画像照合プログラムを提供することを目的とする。

本発明の第１の解決手段によると、
画像照合装置であって、
登録部と、
識別部と
を備え、
前記登録部は、
登録物体より取得した第１の画像から、複数の第１の特徴点の位置を求め、前記登録物体上に印刷しまたはデータベースに記憶し、
および、
前記識別部は、
照合物体より取得した第２の画像から求めた、複数の第２の特徴点の位置と、前記照合物体上に印刷されたまたは前記データベースに記憶された複数の前記第１の特徴点の位置とを照合し、
相互に対応する前記第１の特徴点と前記第２の特徴点の組を予め定められた複数組選択して、前記第１の特徴点と前記第２の特徴点とを相互に対応付ける画像変換関数を表す画像変換係数を算出し、前記画像変換係数の算出に用いる前記複数の組の少なくともひとつを他の組に入れ替えて若しくは削除して、または、他の組を追加して前記画像変換係数を繰り返し再算出し、
再算出した前記画像変換関数に従い算出した前記登録物体の特徴点の位置と前記照合物体の特徴点の位置の相違度または相関度に基づいて前記照合物体の真贋を判定し、
判定結果を表示部に表示させる、及び／又は、記憶部に記憶させる
画像照合装置が提供される。

本発明の第２の解決手段によると、
画像照合装置における画像照合方法であって、
前記画像照合装置は、
登録部と、
識別部と
を備え、
前記登録部は、
登録物体より取得した第１の画像から、複数の第１の特徴点の位置を求め、前記登録物体上に印刷しまたはデータベースに記憶し、
および、
前記識別部は、
照合物体より取得した第２の画像から求めた、複数の第２の特徴点の位置と、前記照合物体上に印刷されたまたは前記データベースに記憶された複数の前記第１の特徴点の位置とを照合し、
相互に対応する前記第１の特徴点と前記第２の特徴点の組を予め定められた複数組選択して、前記第１の特徴点と前記第２の特徴点とを相互に対応付ける画像変換関数を表す画像変換係数を算出し、前記画像変換係数の算出に用いる前記複数の組の少なくともひとつを他の組に入れ替えて若しくは削除して、または、他の組を追加して前記画像変換係数を繰り返し再算出し、
再算出した前記画像変換関数に従い算出した前記登録物体の特徴点の位置と前記照合物体の特徴点の位置の相違度または相関度に基づいて前記照合物体の真贋を判定し、
判定結果を表示部に表示させる、及び／又は、記憶部に記憶させる
画像照合方法が提供される。

本発明の第３の解決手段によると、
画像照合装置における画像照合プログラムであって、
前記画像照合装置の処理部が、登録物体より取得した第１の画像から、複数の第１の特徴点の位置を求め、前記登録物体上に印刷しまたはデータベースに記憶する手順と、
および、
前記処理部が、照合物体より取得した第２の画像から求めた、複数の第２の特徴点の位置と、前記照合物体上に印刷されたまたは前記データベースに記憶された複数の前記第１の特徴点の位置とを照合する手順と、
前記処理部が、相互に対応する前記第１の特徴点と前記第２の特徴点の組を予め定められた複数組選択して、前記第１の特徴点と前記第２の特徴点とを相互に対応付ける画像変換関数を表す画像変換係数を算出し、前記画像変換係数の算出に用いる前記複数の組の少なくともひとつを他の組に入れ替えて若しくは削除して、または、他の組を追加して前記画像変換係数を繰り返し再算出する手順と、
前記処理部が、再算出した前記画像変換関数に従い算出した前記登録物体の特徴点の位置と前記照合物体の特徴点の位置の相違度または相関度に基づいて前記照合物体の真贋を判定する手順と、
前記処理部が、判定結果を表示部に表示させる、及び／又は、記憶部に記憶させる手順と
をコンピュータに実行させるための画像照合プログラムが提供される。

本発明によると、画像を高精度に照合する画像照合装置、画像照合方法、画像照合プログラムを提供することができる。

本発明の各実施例における、透過光による紙紋画像の説明図の例である。 本発明の第一の実施例ならびに第二の実施例における画像照合装置の具体的な構成と、データの流れを示した図である。 本発明の各実施例における紙片の真贋判定のための画像照合システムの構成図の例である。 本発明の第一の実施例における、図２に示した画像照合装置および図３に示した画像照合システムを用いて、紙片の登録処理の流れを示した図である。 本発明の各実施例における、特徴点情報テーブルの説明図の例を示した図である。 本発明の各実施例における、特徴点情報、特徴点情報符号化手段、符号化特徴点情報、特徴点情報復号化手段、の関係を表した図である。 本発明の第一の実施例における、紙片の照合処理の流れ（１）を示した図である。 本発明の各実施例における、対応点テーブルの説明図の例である。 本発明の各実施例における、特徴点の対応関係を示した図である。 本発明の第二の実施例における、紙片の照合処理流れ（２）を示す図である。 本発明の第三の実施例における、反射光による印影画像の説明図の例である。 本発明の第三の実施例における、画像照合システムの構成図の例である。 本発明の第三の実施例における、印鑑登録手段の処理手順を表した図である。 本発明の第一の実施例における、紙片の照合処理の流れ（２）を示した図である。 特徴点の対応を表す説明図である。 対応点の入れ替えについての説明図（１）である。 対応点の入れ替えについての説明図（２）である。 印鑑ＤＢの説明図である。

Ａ．概要

上記目的を解決するため、本実施例の画像照合装置は、例えば、
登録したい物体の画像を取得する手段と、
取得した物体画像から当該物体の複数の特徴点を抽出する手段と、
抽出した複数の特徴点の物体上における位置および特徴量を符号化する手段と、
符号化した特徴点の位置および特徴量を記録する手段と、
記録した特徴点の位置と特徴量を読み取って復号化する手段と、
照合したい物体から特徴点を抽出し、登録した物体の特徴点と照合したい物体から抽出した複数の特徴点同士を対応付ける手段と、
対応付けた複数の特徴点群から照合画像上の残りの特徴点と登録画像上の残りの特徴点の画像変換関数を算出する手段と、
算出した画像変換関数のパラメータを更新する手段と、
更新した画像変換関数を用いて登録画像と照合画像の類似度を判定する手段
を備えることができる。

また、本実施例の他の画像照合装置は、例えば、
登録したい物体の画像を取得する手段と、
取得した物体画像から当該物体の複数の特徴点を抽出する手段と、
抽出した複数の特徴点の物体上における位置および特徴量を符号化する手段と、
符号化した特徴点の位置および特徴量を記録する手段と、
記録した特徴点の位置と特徴量を読み取って復号化する手段と、
照合したい物体から特徴点を抽出し、登録した物体の特徴点と照合したい物体から抽出した複数の特徴点同士を対応付ける手段と、
対応付けた複数の特徴点群から照合画像上の残りの特徴点と登録画像上の残りの特徴点の画像変換関数を算出する手段と、
画像変換関数で求めた特徴点の位置が大きく異なる特徴点群を集め、それらの特徴点群から上述と同様の手段で再帰的に画像変換関数のパラメータを算出する手段と、
算出した画像変換関数のパラメータを更新する手段と、
更新した画像変換関数を用いて登録画像と照合画像の類似度を判定する手段
を備えることができる。

Ｂ．紙片の真贋判定システム

１．システム
本発明の第一の実施例として、照合したい物体として紙片を考え、以下で説明する。本実施例を適用したシステムとしては、紙一枚一枚の真贋を判定する、紙片の真贋判定システムが考えられる。また、具体的な紙片としては、紙幣や株券などの有価証券、チケットやクーポン券、重要書類などが考えられる。なお、以下では主に「紙片」を例に挙げて説明するが、本発明および本実施例は、これに限らず、プラスチック等の化学系素材、木等の天然素材等適宜の素材を用いたシートやその他の物体に適用することができる。
本実施例では、透過光によって紙片の画像を得る場合について述べるが、反射光や非可視光領域の波長を用いて画像を得るようにしてもよい。また、磁気などの特別なセンサによって紙片の画像を得るようにしてもよい。

図１は、紙片の透過画像の説明図である。この画像は市販の透過型のスキャナを用いて、上質紙を６００ｄｐｉの解像度で撮像した画像である。目視しやすいように、濃淡を強調してある。
図２は、第一の実施例の画像照合処理装置の具体的な構成と、データの流れを示した図である。同図に示したように、画像照合システムは、大きく分けて紙片登録部２１と紙片識別部２２を備える。
紙片登録部２１において、透過画像取得部２０２では、紙片２０１の透過画像を取得し、特徴点抽出処理部２０３に送る。特徴点抽出処理部２０３では、入力された透過画像から特徴点を抽出し、特徴点情報符号化部２０４に送る。具体的な特徴点の抽出手順については、図４を用いて後述する。特徴点情報符号化部２０４は、入力された特徴点情報を符号化する。具体的な符号化手順については、図６を用いて後述する。特徴点情報印刷部２０５は、入力された符号化特徴量を、プリンタ３０２を用いて紙片２０１の上に印刷する。なお、紙片登録の際に、紙片２０１の全体の透過画像を取得するのではなく、紙片の一部の領域または複数の部分領域の透過画像を取得するようにしてもよい。この場合は、透過画像を取得する領域座標を画像照合システムに通知する手段が必要になるが、これは公知技術（例えば、特許文献８の［００４２］ないし［００４３］および図１１）を利用することができる。すなわち、取得する領域座標をバーコードに符号化し紙片２０１の上に印刷しておき、そのバーコードを画像照合システムが読み取って、読み取る領域座標を取得するようにすればよい。もしくは、事前に所定の領域座標を決め、画像照合システムに登録しておくようにしても良い。

紙片識別部２２において、紙片の照合時には、反射画像取得部２０７では、真贋判定の対象となる紙片２０６（以降、照合紙片と呼ぶ）の反射画像を取得し、特徴点情報復号化部２０８に送る。特徴点情報復号化部２０８は、入力された符号化特徴量を復号化し、対応点抽出部２１１に送る。具体的な、符号化特徴量の復号化手順については、図６を用いて後述する。また、紙片２０６の透過画像を透過画像取得部２０９で取得し、特徴点抽出処理部２１０に送る。特徴点抽出処理部２１０は、透過画像から特徴点を抽出し、対応点抽出部２１１に送る。紙片登録の際、紙片２０１の透過画像の一部の領域または複数の部分領域を用いた場合は、紙片照合の際も紙片２０６の透過画像の一部の領域または複数の部分領域を用いる。紙片登録部２１の処理の際と同様に、公知技術（例えば、特許文献８の［００４２］ないし［００４３］および図１１）を利用することができる。すなわち、バーコードに符号化して印刷された、取得すべき透過画像の領域座標を画像照合システムが読み取り、紙片２０６の透過画像の中から、その領域座標内の透過画像を抽出するようにすればよい。もしくは、事前に所定の領域座標を決め、画像照合システムに登録しておくようにしても良い。具体的な特徴点の抽出手段については、図４を用いて後述する。対応点抽出部２１１は、特徴点情報復号化部２０８で得られた登録紙片２０１の特徴点情報と、特徴抽出処理部２１０で得られた照合紙片２０６の特徴点情報の対応点を抽出し、対応点テーブル８０１（後述）として出力する。すなわち、登録用の紙片２０１の画像に含まれる複数の特徴点が、照合用の紙片２０６の画像のどの位置に対応するか、その対応関係を算出する。具体的な対応点の抽出方法については、図７を用いて後述する。次に、アフィンパラメータ推定部２１２は、入力された対応点テーブル８０１（後述）から、登録用の紙片２０１の画像を、回転、拡大・縮小することによって照合用の紙片２０６の画像ができたと仮定したとき、最も妥当性のある、回転角度や拡大率・縮小率（これらの数値を、以降ではアフィンパラメータと呼ぶ）を算出する。具体的なアフィンパラメータの算出方法については、図７を用いて後述する。次に、対応点ずれ量・分散算出部２１３では、入力されたアフィンパラメータから、対応点のずれ量・分散を算出する。続いて、対応点更新部２１４では、入力された対応点のずれ量・分散を用いて、アフィンパラメータの推定に用いる対応点の入れ替えを行い、アフィンパラメータを再計算する。評価値が収束するまで、対応点の入れ替えを続ける。つぎに、真贋判定部２１５では、登録紙片２０１の特徴点をアフィンパラメータで変換したときの特徴点座標と、照合紙片２０６から抽出したときの特徴点座標とを比較し、真贋判定を行う。特徴点座標同士の距離が小さければ真券と判断し、特徴点座標同士の距離が大きければ偽券と判断する。真贋判定結果は、表示装置に表示される、及び／又は、ＲＡＭに記憶される。ずれ量・分散の算出方法については、図９を用いて後述する。

図３は、紙片の真贋判定のための画像照合システムの構成図の例である。本システムは、少なくとも１台の、透過型スキャナ（３０１）と、真贋に用いる特徴量を符号化して印刷するプリンタ（３０２）と、真贋判定の計算等を行うＣＰＵ（３０３）と、ＣＰＵを操作するキーボード（３０４）と、符号化した特徴量を読み取る反射型スキャナ（３０９）と、真贋判定結果を操作者に表示する表示装置（３０６）と、真贋判定に用いるコンピュータソフトウェアを格納するＲＡＭ（３０７）と、ＣＰＵを操作するマウス（３０８）と、上述の各々のブロック間のデータを通信するためのバス（３０９）を備える。ＣＰＵ（３０３）は、図２に示した画像照合装置による処理を実行する。真贋判定に用いる特徴量を遠隔地のＤＢ（３１１）に格納し、それらをネットワーク（３１０）で繋ぐような構成にしてもよい。透過型スキャナ（３０１）、反射型スキャナ（３０９）には、市販のスキャナを用いることができる。

図５は、特徴点情報テーブルの説明図の例である。特徴点情報テーブル５０１は、特徴点の数５０２、特徴点ＩＤ５０３、特徴点のＸ座標５０４、特徴点のＹ座標５０５、特徴点の特徴量値５０６を含む。必要に応じて、特徴点の方向５０７を加えてもよい。特徴点の方向は、例えば、右方向を０度とし、半時計回りに回転角に応じて数を増やし、一周すると３６０度とする度数を用いることができる。
図８は、対応点テーブルの説明図の例である。対応点テーブル８０１は、対応点ＩＤ８０２、登録紙片上の特徴点ＩＤ８０３、登録紙片上の特徴点のＸ座標８０４、登録紙片上の特徴点のＹ座標８０５、照合紙片上の特徴点ＩＤ８０６、照合紙片上の特徴点のＸ座標８０７、照合紙片上の特徴点のＹ座標８０８、対応する特徴点同士の類似度８０９を含む。

２．処理

２−１．登録処理
図４は、図２に示した画像照合装置および図３に示した画像照合システムを用いて、紙片の登録処理の流れを示した図である。まず、ＣＰＵ（３０３）は、透過型スキャナ（３０１）により、登録したい紙片の透過画像を撮像し（４０２）、次に当該画像から、紙片の識別に用いる特徴量を抽出する。いま、画像の高さをＨ画素、画像の幅をＷ画素とし、ＣＰＵ（３０３）は、着目画素のＹ座標を表す変数Ｙを１からＨまで変えながら（４０３）、着目画素のＸ座標を表す変数Ｘを１からＷまで変え（４０４）以下の処理をくり返す。ＣＰＵ（３０３）は、座標（Ｘ、Ｙ）を中心とするＮ×Ｎ画素の画像を切り出し、部分画像ＳｕｂＩｍｇ（Ｘ、Ｙ）に格納する（４０５）。続いて、ＣＰＵ（３０３）は、ＳｕｂＩｍｇ（Ｘ、Ｙ）から画像特徴量を抽出する。画像特徴量としては、例えば、画像処理や映像処理で用いられる、濃淡値、色彩、マルチスペクトル情報や、テクスチャ情報、形状情報、空間的な配置情報や、あるいはそれらの組み合わせを用いることができる。以下では、画像に画像フィルタかけて特徴量を抽出する手順について説明する。ＣＰＵ（３０３）は、ＳｕｂＩｍｇ（Ｘ、Ｙ）に画像フィルをかけ、座標（Ｘ、Ｙ）における特徴量Ｆｔｒ（Ｘ、Ｙ）を抽出する（４０６）。画像フィルタには、公知の画像処理フィルタを用いることができる。例えば、田村秀行編著、“コンピュータ画像処理”、オーム社、に記載の勾配フィルタやラプラシアンフィルタ、ロバーツ・オペレータを用いてもよい。また、ＣＰＵ（３０３）は、非特許文献１に記載のＳＩＦＴを使うことで、特徴量に加え、特徴点の方向Ｄｉｒ（Ｘ、Ｙ）を抽出することもできる。続いて、ＣＰＵ（３０３）は、特徴量Ｆｔｒ（Ｘ、Ｙ）と予め定められた閾値とを比較し、特徴量Ｆｔｒ（Ｘ、Ｙ）のほうがその閾値より大きい場合には（４０７）、ＣＰＵ（３０３）は、特徴点情報テーブル５０１に、新しい特徴点ＩＤ５０３を採番して、座標（Ｘ、Ｙ）を追加し（４０８）、特徴量Ｆｔｒ（Ｘ、Ｙ）を追加し（４０９）、方向Ｄｉｒ（Ｘ、Ｙ）を追加する（４１０）。具体的には、図５を参照すると、ＣＰＵ（３０３）は、特徴量を追加するごとに、特徴点情報テーブル５０１の特徴点の数５０２の数を一つ増やし、特徴点情報テーブル５０１の行を一つ追加する。そして、ＣＰＵ（３０３）は、追加した行の特徴点ＩＤ５０３に新しく採番したＩＤを書き込み、Ｘ座標５０４には特徴点のＸ座標を、Ｙ座標５０５には特徴点のＹ座標を、特徴量値５０６には特徴量を、方向５０７には特徴点の方向を書き込む。処理４０７で用いる閾値は、本処理を適用するアプリケーションの要求精度に応じて適切に決めればよいし、また、登録後の紙片の劣化度合いに応じて適切に調整してもよい。透過画像の全ての点について上述の処理が終了したら、ＣＰＵ（３０３）は、特徴点情報テーブル５０１を走査しながら、特徴点情報を順に符号化する。符号化は、図６を用いて後述する。つぎに符号化した特徴点情報を登録用紙片に印刷する（４１２）。なお、ここでは符号化した特徴点情報（符号化特徴点情報）を紙片の上に印刷するとしたが、図３に示すような構成にして、遠隔地にあるＤＢ３１１に、ネットワーク３１０を経由して登録するようにしてもよい。

図６は、特徴点情報、特徴点情報符号化手段、符号化特徴点情報、特徴点情報復号化手段、の関係を表した図である。特徴点情報６０１は、図５に示した特徴点情報テーブル５０１の内容に対応して書かれている。すなわち、ＣＰＵ（３０３）は、特徴点情報テーブル５０１から特徴点の数５０２（図５の例では５）を取り出して、特徴点情報６０１に“５”と記載し、次に、特徴点情報テーブル５０１からＸ座標５０４を取り出して、特徴点情報６０１に“５６”と追記し、特徴点情報テーブル５０１からＹ座標５０５を取り出して、特徴点情報６０１に“１０”と追記し、特徴点情報テーブル５０１から、特徴量値５０６を取り出して、特徴点情報６０１に“２２０”と追記し、特徴点情報テーブル５０１から、特徴点の方向５０７を取り出して、特徴点情報６０１に追記する。
以上の処理を、ＣＰＵ（３０３）は、特徴点の数５０２の数だけ繰り返す。ここでは、特徴点情報６０１にデータを追記していく際、データの区切り目として、カンマ“、”を記載したが、他の手段でデータの区切り目が分かるようにしてもよい。次に、ＣＰＵ（３０３）は、特徴点情報符号化手段６０３を用いて、特徴点情報６０１を符号化し、符号化特徴点情報６０２に変換する。ここでは、ＪＩＳ Ｘ ０５１０に記載のＱＲコード（登録商標）を用いた例を示したが、他にも１次元のバーコードを用いたり、ＭＩＣＲ（Ｍａｇｎｅｔｉｃ Ｉｎｋ Ｃｈａｒａｃｔｅｒ Ｒｅｃｏｇｎｉｔｉｏｎ）の技術を用いて磁性体のインクで印刷するようにしたり、ＯＣＲ（Ｏｐｔｉｃａｌ Ｃｈａｒａｃｔｅｒ Ｒｅｃｏｇｎｉｔｉｏｎ）を用いて人が可読可能な文字で記載するようにしたりしてもよい。その場合は、反射画像取得部２０７は、必要に応じて、磁性体に反応するデバイス等に変更すればよい。符号化特徴点情報６０２を特徴点情報６０１に復号化するにも、上述のＪＩＳ Ｘ ０５１０に記載のＱＲコード（登録商標）や、１次元のバーコードを用いたり、ＭＩＣＲ、ＯＣＲなどの既存の手段を用いたりすることができる。

２−２．照合処理（１）
図７は、紙片の照合処理の流れ（１）を示した図である。まず、ＣＰＵ（３０３）は、反射画像取得７０２で、紙片上に印刷した符号化特徴点情報６０２を取得する。次に、ＣＰＵ（３０３）は、特徴点情報Ａ復号化７０３にて、符号化特徴点情報６０２から、復号化特徴点情報６０１を得る。これを特徴点群Ａとする。次に、ＣＰＵ（３０３）は、透過画像取得７０４で照合紙片の透過画像を取得し、取得した透過画像から特徴点群Ｂを抽出する（７０５）。次に、ＣＰＵ（３０３）は、特徴点群Ａと特徴点群Ｂに含まれる特徴点ａと特徴点ｂ同士の対応をとる。ＣＰＵ（３０３）は、特徴量点同士の対応にあたっては、各特徴点に付与されている特徴量値５０６を利用して、特徴量値５０６が近いもの同士を対応づける。また、ＣＰＵ（３０３）は、各特徴点について特徴点方向５０７が付与されている場合には、特徴点方向５０７を利用し、特徴点方向が近いもの同士を対応付けるようにしてもよい。また、ＣＰＵ（３０３）は、特徴量値５０６、特徴点方向５０７の両方を用いることにしてもよい。次に、ＣＰＵ（３０３）は、対応づけた特徴点同士の類似度が高い順に対応点を並び替える（７０７）。類似度の計算には、例えば、特徴点の特徴量値の差や特徴点の方向の差を用いることができる。並び替えた結果は、図８に示される対応点テーブル８０１の形で与えられる。次に、ＣＰＵ（３０３）は、画像変換パラメータの推定に用いる特徴点数Ｎを予め適宜定める（７０８）。なお、以下では、例としてアフィン変換を用いる場合を説明するので、それに応じて、Ｎ＝３とする。続いて、ＣＰＵ（３０３）は、対応点テーブル８０１から、上位Ｎ個の対応点を抜出し、集合Ｇに代入する（７０９）。次に、ＣＰＵ（３０３）は、変数Ｍに初期値Ｎ＋１を代入する（７１０）。なお、この例ではＭ＝４となる。続いてＣＰＵ（３０３）は、集合Ｇに含まれる対応点からアフィン変換パラメータを推定する（７１１）。アフィン変換とは、ある画像を回転させたり、拡大・縮小したり、平行移動したりして変形した画像を得ることである。ある画像上の座標（ｘ，ｙ）をもつ点を、３×３の行列を使って、変形後の画像上の座標（ｘ^’，ｙ^’）に変換する。すなわち、式１で与えられる。

ここで、ａ、ｂ、ｃ、ｄ、ｅ、ｆ、は変換パラメータであり、（ｘ，ｙ）と（ｘ^’，ｙ^’）の組が３つ与えられれば、６元連立方程式を解く要領で、計算によって求めることができる。集合Ｇから、登録画像上の特徴点１の座標（ｘ_１，ｙ_１）と照合画像上の特徴点１の座標（ｘ’_１，ｙ’_１）、登録画像上の特徴点２の座標（ｘ_２，ｙ_２）と照合画像上の特徴点２の座標（ｘ’_２，ｙ’_２）、登録画像上の特徴点３の座標（ｘ_３，ｙ_３）と照合画像上の特徴点３の座標（ｘ’_３，ｙ’_３）から、ａ、ｂ、ｃ、ｄ、ｅ、ｆを求めることができる。
なお、アフィン変換の代わりに、任意の幾何変換を用いるようにしても良い。例えば、射影変換を用いれば、回転、拡大・縮小、平行移動に加え、台形歪も扱うことができる。その際は、その幾何変換のパラメータを推定するのに必要な座標の数に応じて、特徴点数Ｎの値を適宜変更する。
次に、ＣＰＵ（３０３）は、アフィン変換パラメータの推定に用いなかった登録画像中の特徴点を、上述の手順で算出したアフィン変換パラメータを用いて変換し、照合画像中の座標を求め（７１２）、アフィン変換によって算出した、その座標に最も近い照合画像中の特徴点を探し、特徴点間の座標のずれを算出する（７１３）。続いて、ＣＰＵ（３０３）は、特徴点間の座標のずれの分散を算出し、変数Ｖに代入する（７１４）。アフィン変換によって求めた座標と、その座標に最も近い照合画像中の特徴点は、例えば、図９の点９２４と点９２６である。詳細は図９の説明で後述する。
次に、ＣＰＵ（３０３）は、Ｖ（ｉ）に、集合Ｇに含まれる対応点群から、ｉ番目の対応点をＭ番目の対応点に入れ替えて計算した分散を代入する（７１５）。

図１６は、対応点の入れ替えについての説明図（１）である。すなわち、ＣＰＵ（３０３）は、処理７１４の際、集合Ｇ（１６０１）には、対応点１（１６０２）、対応点２（１６０３）、対応点３（１６０４）の３つの対応点が含まれていたとすると、処理７１５では、集合Ｇ（１６０１）の対応点１（１６０２）を対応点Ｍ（１６０６）に置換した、集合Ｇ１（１６０５）、集合Ｇ（１６０１）の対応点２（１６０３）を対応点Ｍ（１６０６）に置換した、集合Ｇ２（１６０９）、集合Ｇ（１６０１）の対応点３（１６０４）を対応点Ｍ（１６０６）に置換した、集合Ｇ３（１６１３）の３つの集合を算出する。そして、ＣＰＵ（３０３）は、各々の集合に含まれる対応点から分散を算出し、Ｖ（ｉ）に代入する。そして、Ｖ（ｉ）がＶより小さい場合には（７１６）、ＶをＶ（ｉ）で入れ替え（７１７）、ＣＰＵ（３０３）は、集合Ｇのｉ番目の対応点とＭ番目の対応点を入れ替え（７１８）、Ｍの値を一つ増やし（７１９）、ステップ７１０の直後に戻る。
このようにして、ＣＰＵ（３０３）は、集合Ｇ、Ｇ１、Ｇ２、Ｇ３の各分散の中から一番小さい分散の値を求めることができる。
そして、ＣＰＵ（３０３）は、Ｖ（ｊ）とＶの差の絶対値を計算し、その値が予め定めた閾値Ｔｈ１以下の場合（７２０）、さらに、ＣＰＵ（３０３）にＶ（ｉ）を予め定めた閾値Ｔｈ２と比較する。Ｖ（ｉ）が閾値Ｔｈ２よりも小さい場合は、ＣＰＵ（３０３）は、照合紙片は真券と判断し（７２２）、照合処理を終了する（７２３）。そうでない場合は、ＣＰＵ（３０３）は、照合券は偽券と判断し、判定結果を表示装置（３０６）に表示する、及び／又は、ＲＡＭ（３０７）に記憶して（７２４）、照合処理を終了する（７２５）。

このように、本実施例によれば、対応点同士の類似度が高い順に並び替えた対応点から順番にアフィン変換パラメータを推定し、照合処理を行うことで、すべての対応点同士の組み合わせを計算することなく、効率的に照合処理を行うことができる。また、処理７１４にて、ずれの絶対値ではなく、ずれの分散を評価値に用いることで、パラメータの推定が安定してできるようになる。

２−３．照合処理（２）
図１４は、紙片の照合処理の流れ（２）を示した図である。
紙片の照合の流れは図１４のようにしてもよい。以下では、図７の流れと異なる点を中心に説明する。図１４の処理１４０１〜１４０９は、図７の処理７０１〜７０９にそれぞれ対応し、ＣＰＵ（３０３）は、同様の処理を実行する。続いてＣＰＵ（３０３）は、集合Ｇに含まれる対応点からアフィン変換パラメータを推定する（１４１０）。次に、ＣＰＵ（３０３）は、アフィン変換パラメータの推定に用いなかった照合画像中の特徴点を、上述の手順で算出したアフィン変換パラメータを用いて変換し、登録画像中の座標を求め（１４１１）、アフィン変換によって算出した、その座標にもっとも近い登録画像中の特徴点を探し、座標のずれを算出する（１４１２）。続いて、ＣＰＵ（３０３）は、座標のずれの分散を算出し、変数Ｖに代入する（１４１３）。
続いて、ＣＰＵ（３０３）は、処理１４０９で求めた集合Ｇに含まれる対応点群、すなわち、１、２、・・・、Ｎ番目の対応点の中から対応点ｉ（１≦ｉ≦Ｎ）ひとつだけを除外して計算した分散をＶ（ｉ）に代入する（１４１４）。

図１７は、対応点との入れ替えについての説明図（２）である。すなわち、例えば、処理１４１３の際は、図１７（Ａ）のように、集合Ｇ（１６１４）には、対応点１（１６１５）、対応点２（１６１６）、対応点３（１６１７）、・・・、対応点Ｎ−１（１６１８）、対応点Ｎ（１６１９）のＮ個の対応点が含まれていたとする。ＣＰＵ（３０３）は、処理１４１４では、集合Ｇ（１６１４）の対応点１（１６１５）を削除した、集合Ｇ１（１６２０）、集合Ｇ（１６１４）の対応点２（１６１６）を削除した、集合Ｇ２（１６２５）、集合Ｇ（１６１４）の対応点Ｎ（１６１９）を削除した、集合ＧＮ（１６３０）のＮ個の集合を算出する。そして、ＣＰＵ（３０３）は、各々の集合に含まれる対応点から各々の分散を算出し、Ｖ（ｉ）に代入する。
また、ＣＰＵ（３０３）は、処理１４０９で求めた集合Ｇに含まれる対応点群、すなわち、Ｎ＋１、Ｎ＋２、・・・、Ｐ番目の対応点の中から対応点ｊ（Ｎ＋１≦ｊ≦Ｐ）ひとつだけを追加して計算した分散をＶ（ｊ）に代入する（１４１５）。ただし、ｊ＝Ｎ＋１、Ｎ＋２、・・・、Ｐである。ここに、Ｐは、集合Ｇに含まれる対応点の個数（最大値）である。すなわち、処理１４１４の際は、集合Ｇ（１６１４）には、対応点１（１６３６）、対応点２（１６３７）、・・・、対応点Ｎ（１６３８）のＮ個の対応点が含まれていたとする。図１７（Ｂ）のように、ＣＰＵ（３０３）は、処理１４１５では、集合Ｇ（１６３５）に対応点Ｎ＋１（１６４３）を追加した集合ＧＮ＋１（１６３９）、集合Ｇ（１６３５）に対応点Ｎ＋２（１６４８）を追加した集合ＧＮ＋２（１６４４）、・・・、集合Ｇ（１６３５）に対応点Ｐを追加した集合ＧＰ（１６４９）、のＰ−Ｎ個の集合を算出する。
そして、ＣＰＵ（３０３）は、各々の集合に含まれる対応点から各々の分散を算出し、Ｖ（ｉ）に代入する。そして、ＣＰＵ（３０３）は、Ｖ（ｉ）とＶを比較し、ＶよりもＶ（ｉ）のほうが小さければ、ｉ番目の特徴点を除いた分散を計算してＶに代入し（１４１７）、集合Ｇからｉ番目の特徴点を除外し（１４１８）、Ｎの値をひとつ減らし（１４１９）、処理１４１０の直前に戻る（１４２０）。このようにして、ＣＰＵ（３０３）は、集合Ｇ、Ｇ１、Ｇ２、・・・、ＧＮの各手段の中から一番小さい分数の値を求めることができる。また、もし、Ｖ（ｊ）の値がＶよりも小さければ（１４２１）、ＣＰＵ（３０３）は、ｊ番目の特徴点を追加した分散を計算してＶに代入し（１４２２）、集合Ｇにｊ番目の特徴点を追加し（１４２３）、Ｎの値をひとつ増やし（１４２４）、処理１４１０の直前に戻る（１４２５）。このようにして、ＣＰＵ（３０３）は、集合Ｇ、ＧＮ＋１、ＧＮ＋２、・・・、ＧＰの各手段の中から一番小さい分数の値を求めることができる。そして、もし、特徴点ｊを追加したことによる分散の変化量、すなわち｜Ｖ（ｊ）−Ｖ｜が予め定めた閾値Ｔｈ１よりも小さく、かつ、特徴点ｉを減らしたことによる分散の変化量、すなわち｜Ｖ（ｉ）−Ｖ｜が閾値Ｔｈ１よりも小さければ（１４２６）、ＣＰＵ（３０３）は、Ｖを予め定めた閾値Ｔｈ２と比較する（１４２７）。閾値Ｔｈ２よりもＶが小さければ（１４２７）、ＣＰＵ（３０３）は、特徴点ｂを持つ紙片は真券と判定し（１２２８）、真贋判定処理を終了する（１４２９）。処理１４２６の判定結果が偽の場合は、ＣＰＵ（３０３）は、特徴点ｂをもつ紙片は偽券と判定し、判定結果を表示装置（３０６）に表示する、及び／又は、ＲＡＭ（３０７）に記憶して（１４３０）、処理を終了する（１４３１）。

すなわち、本実施例によると、分散Ｖを計算する特徴点の集合Ｇに、特徴点を追加したり、集合Ｇから特徴点を除外して、分散Ｖが最小になるように計算を繰り返し、分散Ｖの変動が閾値以下となった場合に、真券と判定し、閾値以上の場合は偽券と判定する。このような構成にすることにより、最初に、ある選択基準で選んだ特徴点の対を用いて求めたアフィン変換係数が、他の特徴点群を適切に表現できない場合、照合画像は真券を撮影したものであっても、偽券と判定される可能性があるという、公知技術の問題点を解決できる。また、別の公知技術にてランダムに特徴点の候補点対選択した場合に適切な対応点が選択されることがないことがおきる、という公知技術の課題を解決できる。

３．特徴点の対応関係
図９は、本発明の各実施例に用いる、特徴点の対応関係を示した図である。登録用紙片画像９０１ならびに照合紙片画像９０２には、特徴点抽出処理部２０３で抽出した、登録紙片画像上の特徴点１（９０３）、照合紙片画像上の特徴点１（９０４）、登録紙片画像上の特徴点１と照合紙片画像上の特徴点１の対応関係（９０５）、登録紙片画像上の特徴点２（９０６）、照合紙片画像上の特徴点２（９０７）、登録紙片画像上の特徴点２と照合紙片画像上の特徴点２の対応関係（９０８）、登録紙片画像上の特徴点３（９０９）、照合紙片画像上の特徴点３（９１０）、登録紙片画像上の特徴点３と照合紙片画像上の特徴点３の対応関係（９１１）が示されている。
また、登録紙片画像９１２ならびに照合紙片画像９１３には、登録紙片画像上の特徴点１（９１４）、照合紙片画像上の特徴点１（９１５）、登録紙片画像上の特徴点１と照合紙片画像上の特徴点１の対応関係（９１６）、登録紙片画像上の特徴点２（９１７）、照合紙片画像上の特徴点２（９１８）、登録紙片画像上の特徴点２と照合紙片画像上の特徴点２の対応関係（９１９）、登録紙片画像上の特徴点３（９２０）、照合紙片画像上の特徴点３（９２１）、登録紙片画像上の特徴点３と照合紙片画像上の特徴点３の対応関係（９２２）、登録紙片画像上の特徴点４（９２３）、照合紙片画像上の特徴点４（９２４）、登録紙片画像上の特徴点４と照合紙片画像上の特徴量４の対応関係（９２５）、照合紙片画像上の特徴点４に最も近い特徴点（９２６）、照合紙片画像上の特徴点４と照合紙片画像上の特徴点４に最も近い特徴点間の距離（９２７）、登録紙片画像上の特徴点５（９２８）、照合紙片画像上の特徴点５（９２９）、登録紙片画像上の特徴点５と照合紙片画像上の特徴点５の対応関係（９３０）、照合紙片画像上の特徴点５に最も近い特徴点（９３１）、照合紙片上の特徴点５と照合紙片上の特徴点５に最も近い特徴点間の距離（９３２）が示されている。
すなわち、変換行列９１６、９１９、９２２を用いてアフィン変換のパラメータを算出し、算出したアフィン変換のパラメータを用いて登録紙片画像上の特徴点４（９２３）の座標を変換した座標が、照合紙片画像上の特徴点４（９２４）である。同様に、登録紙片画像上の特徴点５（９２８）を変換した座標が、照合紙片画像上の特徴点５（９２９）である。図９には、説明の都合上、特徴点９２４、９２９が記されているが、この位置に必ずしも特徴点が抽出される必要はなく、座標を表している。
前述した、特徴点のずれ量とは、照合紙片画像上の特徴点の座標と、その照合紙片画像上の特徴点に最も近い特徴点の座標の差の平均値のことであり、例えば、差９２７や差９３２の平均値を指す。また、特徴点の分散とは、対応する特徴点同士の座標の差の分散のことである。具体的な算出手段を以下に示す。
いま、ｉ番目の特徴点の紙片画像上でのｘ座標をｐ_ｉ（ｘ）、ｙ座標をｐ_ｉ（ｙ）、ｉ番目の特徴点に対応する特徴点の紙片画像上でのｘ座標をｐ’_ｉ（ｘ）、ｙ座標をｐ’_ｉ（ｙ）、とすると、座標の差Δ_ｉは以下のようにして求めることができる。

そして、対応する特徴点の数がＬ個のとき、特徴点のずれ量は、以下のようにして算出することができる。

また、特徴点の分散は、以下のようにして算出することができる。

Ｃ．画像が局所的に歪んだ場合の紙片の真贋判定システム

本発明の第二の実施例として、第一の実施例において、登録紙片もしくは照合紙片の画像を採取する際、紙の搬送系や光学センサの特性により、画像が局所的に歪んだ場合の例を説明する。
画像照合処理装置の具体的な構成は、図２に示した第一の実施例と同じ構成を用いることができる。また、画像照合処理システムの物理的な構成についても、図３に示した第一の実施例と同じ構成を用いることができる。紙片を登録する流についても、図４に示した第一の実施例と同じ流を用いることができる。
図１０に、紙片の照合処理流れ（２）を示す。ただし、図１０は、第一の実施例の紙片照合の流れを説明した図７と大部分が同じなので、ここでは、異なる部分を中心に説明する。
図１０が図７と異なる点は、特に、ＣＰＵ（３０３）が、アフィン変換パラメータによって変換した特徴点の座標と、対応する特徴点座標のずれを算出した（１０１３）後、ＣＰＵ（３０３）は、ずれの大きな対応点を算出して集合Ｈに格納し（１０１４）、集合Ｈを新たな対応点として再帰的にアフィン変換パラメータを推定する（１０１５）ことである。なお、選択する対応点の個数は実施例１と同様である。
なお、図７の処理７０１〜７１３および処理７１４〜７２５が図１０の処理１００１〜１０１３および処理１０１６〜１０２７にそれぞれ対応し、ＣＰＵ（３０３）は、上述と同様の処理を実行する。

図１５は、特徴点の対応を表す説明図である。
以下では、ずれの大きな対応点を算出して集合Ｈに格納する様子を、図１５を用いて説明する。図１５の登録紙片画像１５０１ならびに照合紙片画像１５０２の上には、登録紙片画像上の特徴点１（１５０３）の座標をアフィン変換によって、照合紙片画像上に変換した特徴点１（１５０５）の座標、それらの特徴点の対の対応関係を示す１５０４、照合紙片画像上に変換した特徴点１（１５０５）に最も近い特徴点１’（１５０６）、特徴点１（１５０５）と特徴点１’（１５０６）の距離１５０７、が図示されている。同様に、登録紙片画像上の特徴点２（１５０８）の座標をアフィン変換によって、照合紙片画像上に変換した特徴点２（１５１０）の座標、それらの特徴点の対の対応関係を示す１５０９、照合紙片画像上に変換した特徴点２（１５１０）に最も近い特徴点２’（１５１１）、特徴点１（１５１０）と特徴点１’（１５１１）の距離１５１２、が図示されている。
また、登録紙片画像上の特徴点３（１５１３）の座標をアフィン変換によって、照合紙片画像上に変換した特徴点３（１５１５）の座標、それらの特徴点の対の対応関係を示す１５１４、照合紙片画像上に変換した特徴点３（１５１５）に最も近い特徴点３’（１５１６）、特徴３（１５１３）と特徴点３’（１５１６）の距離１５１７、が図示されている。さらに、登録紙片画像上の特徴点４（１５１８）の座標をアフィン変換によって、照合紙片画像上に変換した特徴点４（１５２０）の座標、それらの特徴点の対の対応関係を示す１５１９、照合紙片画像上に変換した特徴点４（１５２０）に最も近い特徴点４’（１５２１）、特徴点４（１５２０）と特徴点４’（１５２１）の距離１５２２、が図示されている。なお、変換した図面が煩雑になるので、アフィン変換の算出に用いた特徴点の対と対応関係を示す線分は省いてある。
対応した特徴点座標のずれを算出する処理（１０１３）とは、距離１５０７、１５１２、１５１７、１５２２を算出することである。そして、ずれ量の大きな対応点を算出して集合Ｈに格納する処理（１０１４）とは、距離１５０７、１５１２、１５１７、１５２２を、事前に決めた所定の閾値と比較し、閾値よりも大きいずれ量を持つ特徴点の対、例えば、特徴点１５１８と特徴点１５２１、特徴点１５１３と特徴点１５１６を集合Ｈに含めるということである。
このように第二の実施例では、局所的に特徴点同士の対応をとって照合処理を行うので、部分的に歪みを持つ画像同士の照合を行うことができる。

Ｄ．印鑑照合システム

本発明の第三の実施例として、印影の真贋を判定する、印鑑照合システムを以下で説明する。ここでは、一例として印影について説明するが、本発明および本実施例は、これに限らず適宜の画像に適用可能である。
図１１は、印影画像の説明図の例である。登録用印影画像１１０１について、特徴点抽出処理部２０３で抽出された、特徴点１１０３、特徴点１１０４、特徴点１１０５、特徴点１１０６、特徴点１１０７、特徴点１１０８、特徴点１１０９が示されている。また、照合用印影画像１１０２は、登録用印影画像１１０１を採取するのに使ったものと同じ印鑑であったが、異なる回転角度で採取された印鑑照合時の印影画像の例であり、特徴点抽出処理部２０３で抽出された、特徴点１１１０、特徴点１１１１、特徴点１１１２、特徴点１１１３、特徴点１１１４、特徴点１１１５、特徴点１１１６が示されている。特徴点１１０３が特徴点１１１０に、特徴点１１０４が特徴点１１１１に、特徴点１１０５が特徴点１１１２に、特徴点１１０６が特徴点１１１２に、特徴点１１０７が特徴点１１１３に、特徴点１１０８が特徴点１１１４に、特徴点１１０９が特徴点１１１４に、それぞれ対応する。
また、登録用印影像１１１７について、特徴抽出部２０３で抽出された、特徴点１１１９、特徴点１１２０、特徴点１１２１、特徴点１１２２、特徴点１１２３、特徴点１１２４、特徴点１１２５が示されている。また、照合用印影画像１１０８は、登録用印影画像１１１７を採取するのに使ったものと同じ印鑑であったが、異なる縮尺で採取された印鑑照合時の印影画像の例であり、特徴抽出処理２０３で抽出された、特徴点１１２６、特徴点１１２７、特徴点１１２８、特徴点１１２９、特徴点１１３０、特徴点１１３１、特徴点１１３２が示されている。特徴点１１１９が特徴点１１２６に、特徴点１１２０が特徴点１１２７に、特徴点１１２１が特徴点１１２８に、特徴点１１２２が特徴点１１２９に、特徴点１１２３が特徴点１１３０に、特徴点１１２４が特徴点１１３１に、特徴点１１２５が特徴点１１３２に、それぞれ対応する。

図１２は、第三の実施例の画像照合処理システムの具体的な構成と、データの流れを示した図である。ここでは、照合する画像として印影を用いた例を示す。印影の照合は、役所等の公的機関や銀行等で、申請書に押された印影が、届け出た印鑑の印影と同一であるかを判断する際に行われる。同図に示したように、画像照合システムは、大きく分けて印鑑登録部３１と印鑑照合部３２を備える。印鑑登録部３１において、印鑑登録者が持参した印鑑を紙等に押印した場合、反射画像取得部１２０４はその印影１２０２を取得し、特徴点抽出処理部１２０５へ送る。反射画像の取得には、市販のフラットベッドタイプのスキャナや、オープンタイプのスキャナ、デジタルスチルカメラ等を用いることができる。特徴点抽出処理部１２０５は、反射画像取得部１２０４で得られた反射画像特徴点を抽出する。特徴点の抽出には、第一の実施例の紙片登録手順に関する説明（図４参照）と同様の手順を用いることができる。抽出された特徴点の例を、図１１の特徴点１１０３、１１０４、・・・、１１０９等に示す。次に、特徴点情報符号化部１２０６は、特徴点抽出処理部１２０５を用いて抽出された特徴点を符号化する。具体的な符号化手順については、第一の実施例の特徴点の符号化手段の説明（図６参照）と同様の手順を用いることができる。また、印鑑所持者ＩＤ１２０１をＩＤ取得部１２０３で取得する。印鑑所持者ＩＤには、例えば、氏名や社会保険番号等の個人を特定できる番号を用いることができる。また、銀行に印鑑を登録する場合は口座番号等の印鑑所持者を特定できる番号を用いることもできる。続いて、ＤＢ登録部１２０７では、ＩＤ取得部１２０３で取得した印鑑所持者ＩＤ１２０１の情報と、特徴点情報符号化部１２０６で得られた符号化特徴点情報を紐づけて、印鑑ＤＢ１２０８へ送信し、印鑑ＤＢ１２０８に、その情報を記憶する。ここでは、特徴点の情報を印鑑ＤＢ１２０８に登録する形態を示したが、第一の実施例と同様に、紙片の上に印刷する形態にしてもよい。

図１８に印鑑ＤＢの説明図を示す。印鑑ＤＢ１２０８は、印鑑所持者ＩＤに対して、印影の符号化特徴点情報を記憶する。
印影の照合時には、印鑑照合部３２において、照合の対象となる印影１２１０（以降、照合印影と呼ぶ）の反射画像を反射画像取得部１２１４で取得し、特徴点抽出処理部１２１５に送る。特徴点抽出処理部１２１５では、照合印影の特徴点を抽出し、対応点抽出処理部１２１６に送る。また、印鑑所持者ＩＤ１２０９を、ＩＤ取得部１２１１で取得し、ＤＢ参照部１２１２に送る。ＤＢ参照部１２１２は、印鑑所持者ＩＤ１２０９をキーに印鑑ＤＢ１２０８を検索し、印鑑所持者ＩＤ１２０９に紐づけられた符号化特徴点情報を取得し、特徴点情報復号化部１２１３へ送信する。特徴点情報復号化部１２１３は、符号化特徴点情報を特徴点情報に復号化し、対応点抽出部１２１６へ送る。具体的な、符号化特徴点情報の復号化手順については、第一の実施例の説明（図６参照）を用いた手順を用いることができる。対応点抽出部１２１６は、印鑑登録の際に登録した印影の特徴点と、印鑑照合の際に得られた印影の特徴点の対応づけを行い、対応点テーブル８０１の形で出力する（図８参照）。対応点の抽出には、第一の実施例の説明（図７参照）での手順を用いることができる。次に、アフィンパラメータ推定部１２１７は、入力された対応点テーブル８０１から、登録用の印影１２０２の画像を、回転、拡大・縮小することによって照合用の印影１２１０の画像ができたと仮定したとき、最も妥当性のある、回転角度や拡大率・縮小率を算出する。具体的なアフィンパラメータの算出方法については、第一の実施例の説明（図７参照）での手順を用いることができる。次に、対応点ずれ量・分散算出部１２１８では、入力されたアフィンパラメータから、対応点のずれ量・分散を算出する。続いて、対応点更新部１２１９では、入力された対応点のずれ量・分散を用いて、アフィンパラメータの推定に用いる対応点の入れ替えを行い、アフィンパラメータを再計算する。このとき、対応点更新部１２１９は評価値が収束するまで、対応点の入れ替えを続ける。つぎに、真贋判定部１２２０では、登録印影１２０２の特徴点をアフィンパラメータで変換したときの特徴点座標と、照合印影１２１０から抽出したときの特徴点座標とを比較し、真贋判定を行う。このとき、特徴点座標同士の距離が小さければ真券と判断し、特徴点座標同士の距離が大きければ偽券と判断する。

Ｅ．実施例の効果

本発明または本実施例の第一の効果は、登録画像の特徴点と照合画像の特徴点を関係づける画像変換関数のパラメータを高精度に算出できることである。
また、本発明または本実施例の第二の効果は、登録画像に対して照合画像が、部分的に回転していたり、拡大・縮小していたりしても、高精度に画像を照合できることである。
また、本発明および本実施例の第三の効果は、照合したい物体の画像を撮像したときに発生する、画像の回転や歪に対して頑健な認証技術を提供できることである。

Ｆ．付記

以上では主に登録画像と照合画像の特徴点間のずれの分散を用いる場合について説明したが、これに限らず本発明または本実施例は、相関度、相違度等の適宜の指標を用いることができる。
また、本発明の画像照合方法又は画像照合装置・システムは、その各手順をコンピュータに実行させるための画像照合プログラム、画像照合プログラムを記録したコンピュータ読み取り可能な記録媒体、画像照合プログラムを含みコンピュータの内部メモリにロード可能なプログラム製品、そのプログラムを含むサーバ等のコンピュータ、等により提供されることができる。
なお、本発明は上記した実施例に限定されるものではなく、様々な変形例が含まれている。例えば、上記した実施例は本発明を分かりやすく説明するために詳細に説明したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。また、ある実施例の構成の一部を他の実施例の構成に置き換えることが可能であり、また、ある実施例の構成に他の実施例の構成を加えることも可能である。また、各実施例の構成の一部について、他の構成の追加・削除・置換をすることが可能である。
また、上記の各構成、機能、処理部、処理手段等は、それらの一部または全部を、例えば集積回路で設計する等によりハードウェアで実現してもよい。また、上記の各構成、機能等は、プロセッサがそれぞれの機能を実現するプログラムを解釈し、実行することによりソフトウェアで実現してもよい。各機能を実現するプログラム、テーブル、ファイル等の情報は、メモリや、ハードディスク、ＳＳＤ（Ｓｏｌｉｄ Ｓｔａｔｅ Ｄｒｉｖｅ）等の記録装置、または、ＩＣカード、ＳＤカード、ＤＶＤ等の記録媒体に置くことができる。
また、制御線や情報線は説明上必要と考えられるものを示しており、製品上必ずしも全ての制御線や情報線を示しているとは限らない。実際には殆ど全ての構成が相互に接続されていると考えてもよい。

１０１・・・透過光による紙紋画像、２０１・・・紙片、２０２・・・透過画像取得部、２０３・・・特徴点抽出処理部、２０４・・・特徴点情報符号化部、２０５・・・特徴点情報印刷部、２０６・・・紙片、２０７・・・反射画像取得部、２０８・・・特徴点情報複合化部、２０９・・・透過画像取得部、２１０・・・特徴点抽出処理部、２１１・・・対応点抽出部、２１２・・・アフィンパラメータ推定部、２１３・・・対応点ずれ量・分散算出部、２１４・・・対応点更新部、２１５・・・真贋判定部、３０１・・・透過型スキャナ、３０２・・・プリンタ、３０３・・・ＣＰＵ、３０４・・・キーボード、３０５・・・反射型スキャナ、３０６・・・表示装置、３０７・・・ＲＡＭ、３０８・・・マウス、３０９・・・バス、３１０・・・ネットワーク、３１１・・・ＤＢ、５０２・・・特徴点情報テーブル、５０３・・・特徴点ＩＤ、５０４・・・特徴点Ｘ座標、５０５・・・特徴点Ｙ座標、５０６・・・特徴値、５０７・・・特徴点方向、６０１・・・特徴点情報、６０２・・・符号化特徴点情報、６０３・・・特徴点情報符号化処理手段、６０４・・・特徴点情報複合化手段、８０１・・・対応点テーブル、８０２・・・対応点ＩＤ、８０３・・・登録紙片上の特徴点ＩＤ、８０４・・・登録紙片上の特徴点のＸ座標、８０５・・・登録紙片上の特徴点のＹ座標、８０６・・・照合紙片上の特徴点のＩＤ、８０７・・・照合紙片上の特徴点のＸ座標、８０８・・・照合紙片上の特徴点のＹ座標、８０９・・・対応点類似度、９０１・・・登録紙片画像、９０２・・・照合紙片画像、９０３・・・登録紙片画像上の特徴点１、９０４・・・照合紙片画像上の特徴点１、９０５・・・登録紙片画像上の特徴点１と照合紙片画像上の特徴点１の対応、９０６・・・登録紙片上の特徴点２、９０７・・・照合紙片上の特徴点２、９０８・・・登録紙片画像上の特徴点２と照合紙片画像上の特徴点２の対応、９０９・・・登録紙片画像上の特徴点３、９１０・・・照合紙片画像上の特徴点３、９１１・・・登録紙片画像上の特徴点３と照合紙片画像上の特徴点３の対応、９１２・・・登録紙片画像、９１３・・・照合紙片画像、９１４・・・登録紙片画像上の特徴点１、９１５・・・照合紙片画像上の特徴点１、９１６・・・登録紙片画像上の特徴点１と照合紙片画像上の特徴点１の対応、９１７・・・登録紙片画像上の特徴点２、９１８・・・照合紙片画像上の特徴点２、９１９・・・登録紙片画像上の特徴点２と照合紙片画像上の特徴点２の対応、９２０・・・登録紙片画像上の特徴点３、９２１・・・照合紙片画像上の特徴点３、９２２・・・登録紙片画像上の特徴点３と照合紙片画像上の特徴点３、９２３・・・登録紙片画像上の特徴点４、９２４・・・照合紙片画像上の特徴点４、９２５・・・登録紙片画像上の特徴点４と照合紙片画像上の特徴量４の対応、９２６・・・照合紙片画像上の特徴点４に最も近い特徴点、９２７・・・照合紙片画像上の特徴点４と照合紙片画像上の特徴点４に最も近い特徴点間の距離、９２８・・・登録紙片画像上の特徴点５、９２９・・・照合紙片画像上の特徴点５、９３０・・・登録紙片画像上の特徴点５と照合紙片画像上の特徴点５の対応、９３１・・・照合紙片画像上の特徴点５に最も近い特徴点、９３２・・・照合紙片画像上の特徴点５と照合紙片画像上の特徴点５に最も近い特徴点間の距離、１１０１・・・登録用印影画像、１１０２・・・照合用印影画像、１１０３・・・登録用印影画像中の特徴点１、１１０４・・・登録用印影画像中の特徴点２、１１０５・・・登録用印影画像中の特徴点３、１１０６・・・登録用印影画像中の特徴点４、１１０７・・・登録用印影画像中の特徴点５、１１０８・・・登録用印影画像中の特徴点６、１１０９・・・登録用印影画像中の特徴点７、１１１０・・・照合用印影画像中の特徴点１、１１１１・・・照合用印影画像中の特徴点２、１１１２・・・照合用印影画像中の特徴点３、１１１３・・・照合用印影画像中の特徴点４、１１１４・・・・照合用印影画像中の特徴点５、１１１５・・・照合用印影画像中の特徴点６、１１１６・・・照合用印影画像中の特徴点７、１１１７・・・登録用印影画像、１１１８・・・照合用印影画像、１１１９・・・登録用印影画像中の特徴点１、１１２０・・・登録用印影画像中の特徴点２、１１２１・・・登録用印影画像中の特徴点３、１１２２・・・登録用印影画像中の特徴点４、１１２３・・・登録用印影画像中の特徴点５、１１２４・・・登録用印影画像中の特徴点６、１１２５・・・登録用印影がオズ中の特徴点７、１１２６・・・照合用印影画像中の特徴点１、１１２７・・・照合用印影画像中のと、有頂天２、１１２８・・・照合用印影画像中の特徴点３、１１２９・・・照合用印影画像中の特徴点４、１１３０・・・照合用印影画像中の特徴点５、１１３１・・・照合用印影画像中の特徴点６、１１３２・・・照合用印影画像中の特徴点７、１２０１・・・印鑑所持者ＩＤ、１２０２・・・印影、１２０３・・・ＩＤ取得部、１２０４・・・反射画像取得部、１２０５・・・特徴点抽出処理部、１２０６・・・特徴点情報符号化部、１２０７・・・ＤＢ登録部、１２０８・・・印鑑ＤＢ、１２０９・・・印鑑所持者ＩＤ、１２１０・・・印影、１２１１・・・ＩＤ取得部、１２１２・・・ＤＢ参照部、１２１３・・・特徴点情報復号化部、１２１４・・・透過画像取得部、１２１５・・・特徴点抽出処理部、１２１６・・・対応点抽出部、１２１７・・・アフィンパラメータ推定部、１２１８・・・対応点ずれ量・分散算出部、１２１９・・・対応点更新部、１２２０・・・真贋判定部

Claims (15)

1. 画像照合装置であって、
   登録部と、
   識別部と
   を備え、
   前記登録部は、
   登録物体より取得した第１の画像から、複数の第１の特徴点の位置を求め、前記登録物体上に印刷しまたはデータベースに記憶し、
   および、
   前記識別部は、
   照合物体より取得した第２の画像から求めた、複数の第２の特徴点の位置と、前記照合物体上に印刷されたまたは前記データベースに記憶された複数の前記第１の特徴点の位置とを照合し、
   相互に対応する前記第１の特徴点と前記第２の特徴点の組を予め定められた複数組選択して、前記第１の特徴点と前記第２の特徴点とを相互に対応付ける画像変換関数を表す画像変換係数を算出し、前記画像変換係数の算出に用いる前記複数の組の少なくともひとつを他の組に入れ替えて若しくは削除して、または、他の組を追加して前記画像変換係数を繰り返し再算出し、
   再算出した前記画像変換関数に従い算出した前記登録物体の特徴点の位置と前記照合物体の特徴点の位置の相違度または相関度に基づいて前記照合物体の真贋を判定し、
   判定結果を表示部に表示させる、及び／又は、記憶部に記憶させる
   画像照合装置。
   
2. 請求項１に記載された画像照合装置において、
   前記登録部は、
   登録物体の透過または反射により前記第１の画像を取得し、
   取得した前記第１の画像から該登録物体中または該登録物体上の複数の第１の特徴点を抽出し、
   抽出した複数の前記第１の特徴点の位置および特徴量を含む第１の特徴点情報を求め、特徴量が予め定められた閾値より大きい前記第１の特徴点情報を特徴点情報テーブルに記憶し、
   前記特徴点情報テーブルを参照し、前記第１の特徴点情報を符号化して第１の符号化特徴点情報を求め、
   前記第１の符号化特徴点情報を前記登録物体上に印刷するまたはデータベースに記憶する
   ことを特徴とする画像照合装置。
   
3. 請求項１に記載された画像照合装置において、
   上記識別部は、
   前記第１の画像を取得した透過または反射と同様に、真贋判定をしたい照合物体の前記第２の画像を取得し、
   取得した前記第２の画像から該照合物体中または該照合物体上の複数の第２の特徴点を抽出し、
   抽出した複数の前記第２の特徴点の位置および特徴量を含む第２の特徴点情報を求め、
   前記照合物体上に印刷されたまたは前記データベースに記憶された前記第１の符号化特徴情報を読み取って復号化して第１の特徴点情報を求め、
   複数の前記第１の特徴点情報に含まれる複数の前記第１の情報点と、複数の前記第２の特徴点情報に含まれる複数の前記第２の情報点とを照合する
   ことを特徴とする画像照合装置。
   
4. 請求項１に記載された画像照合装置において、
   前記識別部は、
   前記登録物体の前記第１の特徴点と、該第１の特徴点に対応する前記照合物体の前記第２の特徴点を抽出し、類似度の大きい順に対応関係を対応点テーブルとして記憶し、
   前記対応点テーブルを参照して、前記画像変換関数に応じて定められる数の複数の対応データを選択し、選択した複数の対応データに基づいて、前記登録物体の前記第１の画像を前記照合物体の前記第２の画像に変換する前記画像変換関数の前記変換関数係数を算出し、
   前記第１の特徴点の位置を前記変換関数係数で変換した位置と、対応する前記第２の特徴点の位置との間のずれ量および該ずれ量の分散を算出し、
   前記ずれ量の分散を用いて、前記変換関数係数の推定に用いる前記対応データの少なくともひとつの入れ替え又は削除又は追加を行い、前記変換関数係数を再計算し、評価値が収束するまで、前記対応データの少なくともひとつの入れ替え又は削除又は追加を続け、
   前記第１の特徴点の位置を前記変換関数係数で変換した特徴点の位置と、前記第２の特徴点の位置との距離が予め定められた閾値より小さければ真と判断し、大きければ偽と判断する
   ことを特徴とする画像照合装置。
   
5. 請求項４に記載された画像照合装置において、
   前記識別部は、
   前記画像変換係数の推定に用いる前記第１の特徴点および前記第２の特徴点の前記対応データの組を、前記画像変換関数を求めるために必要な数であるＮ個抜出して集合Ｇとし、前記集合Ｇに含まれる複数の前記対応点から、ｉ（ｉ＝１〜Ｎ）番目の対応データをＮ＋１番目の対応データに入れ替えた複数の集合Ｇ（ｉ）について、前記第１の特徴点を画像変換した特徴点座標と、前記第１の特徴点に対応する前記第２の特徴点の特徴点座標との間の座標のずれの分散をそれぞれ算出し、
   算出した前記分散の一番小さい値が、予め定められた閾値よりも小さい場合は、前記照合物体は真と判断し、小さくない場合は偽と判断する
   ことを特徴とする画像照合装置。
   
6. 請求項４に記載された画像照合装置において、
   前記識別部は、
   前記画像変換係数の推定に用いる前記第１の特徴点および前記第２の特徴点の対応データの組を、前記画像変換関数を求めるために必要な数であるＮ個抜出して集合Ｇとし、前記集合Ｇに含まれる対応点の中から対応点ｉ（１≦ｉ≦Ｎ）ひとつを除外した複数の集合について、前記第１の特徴点を画像変換した特徴点座標と、前記第２の特徴点の特徴点座標との間の座標のずれを表す第１の分散Ｖ（ｉ）をそれぞれ計算し、
   前記集合Ｇに含まれる対応点の中から対応点ｊ（Ｎ＋１≦ｊ≦Ｐ、Ｐは、集合Ｇに含まれる対応データの個数）ひとつを追加した複数の集合について、前記第１の特徴点を画像変換した特徴点座標と、前記第１の特徴点に対応する前記第２の特徴点の特徴点座標の間の座標のずれを表す第２の分散Ｖ（ｊ）を計算し、
   算出した前記第１の分散Ｖ（ｉ）の一番小さい値および前記第２の分散Ｖ（ｊ）の一番小さい値が、それぞれ予め定められた閾値よりも小さければ真と判定し、小さくなければ偽と判定する
   ことを特徴とする画像照合装置。
   
7. 請求項４に記載された画像照合装置において、
   前記識別部は、
   前記ずれの比較的大きな対応点を、前記画像変換関数に応じて定められる数の複数選択して集合Ｈとし、前記集合Ｈを新たな対応点として再帰的に前記画像変換係数を推定することを特徴とする画像照合装置。
   
8. 請求項１に記載された画像照合装置において、
   前記登録部は、
   ＩＤ情報と、前記第１の特徴点情報を対応づけて前記データベースへ記憶または前記登録物体に印刷し、
   前記識別部は、
   前記データベースまたは前記印刷を参照し、前記ＩＤ情報に従い前記第１の符号化特徴点情報を取得する
   ことを特徴とする画像照合装置。
   
9. 請求項１に記載された画像照合装置において、
   前記登録部または前記識別部は、
   前記第１の画像を取得する際に、前記登録物体の全体の画像を取得すること、または、前記登録物体の一部の領域の画像を取得すること、または、前記登録物体の複数の部分領域の画像を取得することを特徴とする画像照合装置。
   
10. 請求項９に記載された画像照合装置において、
    前記登録部は、前記第１の画像を取得する領域をコードに符号化し、前記登録物体上に印刷しまたはデータベースに記憶し、
    前記識別部は、印刷されたコードまたは前記データベースに記憶されたコードから前記領域を求め、前記領域に従い前記照合画像から前記第２の画像を抽出することを特徴とする画像照合装置。
    
11. 請求項９に記載された画像照合装置において、
    前記登録部および前記識別部は、予め前記領域を定め、前記領域に従い、前記第１の画像および前記第２の画像をそれぞれ取得することを特徴とする画像照合装置。
    
12. 請求項１に記載された画像照合装置において、
    前記登録部または前記識別部は、
    前記第１の画像または前記第２の画像を、透過光若しくは反射光によって得ること、または、非可視光領域の波長を用いて得ること、または、磁気若しくは他のセンサによって得ることを特徴とする画像照合装置。
    
13. 請求項１に記載された画像照合装置において、
    前記登録部は、前記第１の特徴点情報を、他装置にある前記データベースにネットワークを経由して送信して登録することを特徴とする画像照合装置。
    
14. 画像照合装置における画像照合方法であって、
    前記画像照合装置は、
    登録部と、
    識別部と
    を備え、
    前記登録部は、
    登録物体より取得した第１の画像から、複数の第１の特徴点の位置を求め、前記登録物体上に印刷しまたはデータベースに記憶し、
    および、
    前記識別部は、
    照合物体より取得した第２の画像から求めた、複数の第２の特徴点の位置と、前記照合物体上に印刷されたまたは前記データベースに記憶された複数の前記第１の特徴点の位置とを照合し、
    相互に対応する前記第１の特徴点と前記第２の特徴点の組を予め定められた複数組選択して、前記第１の特徴点と前記第２の特徴点とを相互に対応付ける画像変換関数を表す画像変換係数を算出し、前記画像変換係数の算出に用いる前記複数の組の少なくともひとつを他の組に入れ替えて若しくは削除して、または、他の組を追加して前記画像変換係数を繰り返し再算出し、
    再算出した前記画像変換関数に従い算出した前記登録物体の特徴点の位置と前記照合物体の特徴点の位置の相違度または相関度に基づいて前記照合物体の真贋を判定し、
    判定結果を表示部に表示させる、及び／又は、記憶部に記憶させる
    画像照合方法。
    
15. 画像照合装置における画像照合プログラムであって、
    前記画像照合装置の処理部が、登録物体より取得した第１の画像から、複数の第１の特徴点の位置を求め、前記登録物体上に印刷しまたはデータベースに記憶する手順と、
    および、
    前記処理部が、照合物体より取得した第２の画像から求めた、複数の第２の特徴点の位置と、前記照合物体上に印刷されたまたは前記データベースに記憶された複数の前記第１の特徴点の位置とを照合する手順と、
    前記処理部が、相互に対応する前記第１の特徴点と前記第２の特徴点の組を予め定められた複数組選択して、前記第１の特徴点と前記第２の特徴点とを相互に対応付ける画像変換関数を表す画像変換係数を算出し、前記画像変換係数の算出に用いる前記複数の組の少なくともひとつを他の組に入れ替えて若しくは削除して、または、他の組を追加して前記画像変換係数を繰り返し再算出する手順と、
    前記処理部が、再算出した前記画像変換関数に従い算出した前記登録物体の特徴点の位置と前記照合物体の特徴点の位置の相違度または相関度に基づいて前記照合物体の真贋を判定する手順と、
    前記処理部が、判定結果を表示部に表示させる、及び／又は、記憶部に記憶させる手順と
    をコンピュータに実行させるための画像照合プログラム。

Images