(57)【特許請求の範囲】
【請求項１】
第１の期間に撮像された被写体の画像の画像データから第１の特徴量を抽出する第１の抽出手段と、
前記第１の特徴量の第１の平均値を、前記第１の期間よりも過去の第２の期間に撮像された前記被写体の画像の画像データから抽出された第２の特徴量の第２の平均値と、前記第１の期間よりも過去の期間であって、前記第２の期間とは異なる第３の期間に撮像された前記被写体の画像の画像データから抽出された第３の特徴量の第３の平均値とに基づいて推定する第１の推定手段と、
前記第１の平均値に対する前記第１の特徴量のばらつきの幅を示す第１の分布幅を、前記第２の平均値に対する前記第２の特徴量のばらつきの幅を示す第２の分布幅、および前記第３の平均値に対する前記第３の特徴量のばらつきの幅を示す第３の分布幅に基づいて推定する第２の推定手段と、
前記第１の特徴量の値が、前記第１の平均値および前記第１の分布幅に基づく所定の範囲内の値であるか否かを判定する判定手段と
を備える画像処理装置。
【請求項２】
前記第２の推定手段は、前記第２の分布幅および前記第３の分布幅のうちの大きい値を前記第１の分布幅とする
請求項１に記載の画像処理装置。
【請求項３】
前記第２の推定手段は、前記第１の分布幅をｄ３、前記第２の特徴量および前記第３の特徴量をｘ・ｙ、前記第２の平均値をｍ１、前記第３の平均値をｍ２、前記第２の特徴量の数をｎ１、前記第３の特徴量の数をｎ２とした場合に、次式
【数１】
0004462176.tif 000005
を計算することにより求められる値を前記第１の分布幅とする
請求項１に記載の画像処理装置。
【請求項４】
前記第２の期間に撮像された前記被写体の画像の画像データから前記第２の特徴量を抽出し、前記第３の期間に撮像された前記被写体の画像の画像データから前記第３の特徴量を抽出する第２の抽出手段をさらに備える
請求項１に記載の画像処理装置。
【請求項５】
特徴空間上における複数の前記第２の特徴量の分布と、複数の前記第３の特徴量の分布とを分離するための判別軸を算出する算出手段をさらに備え、
前記第１の抽出手段は、前記第１の期間に撮像された画像の画像データから抽出された特徴量と、前記判別軸との内積の値を前記第１の特徴量とし、
前記第２の抽出手段は、前記第２の期間に撮像された画像の画像データから抽出された特徴量と、前記判別軸との内積の値を前記第２の特徴量とし、前記第３の期間に撮像された画像の画像データから抽出された特徴量と、前記判別軸との内積の値を前記第３の特徴量とする
請求項４に記載の画像処理装置。
【請求項６】
第１の期間に撮像された被写体の画像の画像データから第１の特徴量を抽出し、
前記第１の特徴量の第１の平均値を、前記第１の期間よりも過去の第２の期間に撮像された前記被写体の画像の画像データから抽出された第２の特徴量の第２の平均値と、前記第１の期間よりも過去の期間であって、前記第２の期間とは異なる第３の期間に撮像された前記被写体の画像の画像データから抽出された第３の特徴量の第３の平均値とに基づいて推定し、
前記第１の平均値に対する前記第１の特徴量のばらつきの幅を示す第１の分布幅を、前記第２の平均値に対する前記第２の特徴量のばらつきの幅を示す第２の分布幅、および前記第３の平均値に対する前記第３の特徴量のばらつきの幅を示す第３の分布幅に基づいて推定し、
前記第１の特徴量の値が、前記第１の平均値および前記第１の分布幅に基づく所定の範囲内の値であるか否かを判定する
ステップを含む画像処理方法。
【請求項７】
第１の期間に撮像された被写体の画像の画像データから第１の特徴量を抽出する抽出ステップと、
前記第１の特徴量の第１の平均値を、前記第１の期間よりも過去の第２の期間に撮像された前記被写体の画像の画像データから抽出された第２の特徴量の第２の平均値と、前記第１の期間よりも過去の期間であって、前記第２の期間とは異なる第３の期間に撮像された前記被写体の画像の画像データから抽出された第３の特徴量の第３の平均値とに基づいて推定する第１の推定ステップと、
前記第１の平均値に対する前記第１の特徴量のばらつきの幅を示す第１の分布幅を、前記第２の平均値に対する前記第２の特徴量のばらつきの幅を示す第２の分布幅、および前記第３の平均値に対する前記第３の特徴量のばらつきの幅を示す第３の分布幅に基づいて推定する第２の推定ステップと、
前記第１の特徴量の値が、前記第１の平均値および前記第１の分布幅に基づく所定の範囲内の値であるか否かを判定する判定ステップと
をコンピュータに実行させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は画像処理装置および方法、並びにプログラムに関し、特に、個人の認証を行う場合に用いて好適な画像処理装置および方法、並びにプログラムに関する。
【背景技術】
【０００２】
被認証者の顔の画像（以下、顔画像とも称する）を利用して個人の認証を行う、顔画像認証が知られている。
【０００３】
従来、顔画像認証においては、図１に示すように、予め記録されている、過去に撮像された複数枚の被認証者の顔画像１１−１乃至顔画像１１−３と、顔画像認証時に新たに撮像された被認証者の顔画像１２とが比較されて、被認証者が本人であるか否かの認証が行われる。
【０００４】
過去の顔画像１１−１乃至顔画像１１−３の顔画像データのそれぞれ、および新たに撮像された顔画像１２の顔画像データには、所定の変換処理が施されて、それぞれの画像データから被認証者の顔の特徴の特徴量が抽出される。そして、顔画像１１−１乃至顔画像１１−３の顔画像データのそれぞれから抽出された特徴量のそれぞれからなる、特徴空間上の空間（以下、認証空間と称する）内に、顔画像１２の顔画像データから抽出された特徴量が位置するか否かが判定される。
【０００５】
顔画像１２の顔画像データから抽出された特徴量が、認証空間内に位置すると判定された場合、新たに撮像された顔画像１２の被認証者は本人であるとされて認証され、顔画像１２の顔画像データから抽出された特徴量が、認証空間内に位置しないと判定された場合、新たに撮像された顔画像１２の人物は、本人ではないと認証されて、その後の処理は拒絶される。
【０００６】
また、顔画像認証を行う場合に、過去に撮像された複数枚の顔画像を予め記録しておくことで、複数の顔画像の特徴量の平均を求め、被認証者の顔画像の平均的（中心的）な状態を推測（推定）したり、顔画像ごとに変化する特徴量のばらつきの幅を推測したりして、顔画像認証による認証の精度をより向上させることができる。
【０００７】
さらに、予め一般的な顔の各部位の経時変化の統計データと、過去に撮像された被認証者の顔画像（の画像データ）とを記憶しておき、統計データを用いて記憶している被認証者の顔画像の各部位を変形させて認証を行う認証装置もある（例えば、特許文献１参照）。この認証装置においては、統計データを用いて変形された過去の被認証者の顔画像と、認証時に新たに撮像された被認証者の顔画像とが比較されて認証が行われる。
【０００８】
【特許文献１】特開平０９−１３４４２８号公報
【発明の開示】
【発明が解決しようとする課題】
【０００９】
しかしながら、上述した技術では、時間の経過による被認証者の顔の変化、例えば、被認証者の髪形の変化や、加齢による顔の変化などは考慮されず、被認証者の顔が時間の経過とともにどのように変化するかを推測することはできなかった。
【００１０】
したがって、特徴空間上の認証空間の位置は、時間の経過とともに変化することはないので、過去に撮像された被認証者の顔の写真などを用いた、悪意の第３者によるなりすましを防止することが困難であった。また、被認証者の髪が伸びたり、加齢により被認証者の顔が変わったりすると、被認証者本人であっても本人であると認証されなくなってしまう可能性があった。
【００１１】
また、統計データを用いて過去の被認証者の顔画像を変形させる方法においては、統計データを用いて顔画像を変形させるので、被認証者の顔の経時変化にある程度は対応することができるが、処理量が多く認証に時間がかかってしまうため、実現は困難であった。
【００１２】
本発明は、このような状況に鑑みてなされたものであり、より簡単かつ確実に顔画像認証を行うことができるようにするものである。
【課題を解決するための手段】
【００１３】
本発明の一側面の画像処理装置は、第１の期間に撮像された被写体の画像の画像データから第１の特徴量を抽出する第１の抽出手段と、前記第１の特徴量の第１の平均値を、前記第１の期間よりも過去の第２の期間に撮像された前記被写体の画像の画像データから抽出された第２の特徴量の第２の平均値と、前記第１の期間よりも過去の期間であって、前記第２の期間とは異なる第３の期間に撮像された前記被写体の画像の画像データから抽出された第３の特徴量の第３の平均値とに基づいて推定する第１の推定手段と、前記第１の平均値に対する前記第１の特徴量のばらつきの幅を示す第１の分布幅を、前記第２の平均値に対する前記第２の特徴量のばらつきの幅を示す第２の分布幅、および前記第３の平均値に対する前記第３の特徴量のばらつきの幅を示す第３の分布幅に基づいて推定する第２の推定手段と、前記第１の特徴量の値が、前記第１の平均値および前記第１の分布幅に基づく所定の範囲内の値であるか否かを判定する判定手段とを備える。
【００１５】
前記第２の推定手段には、前記第２の分布幅および前記第３の分布幅のうちの大きい値を前記第１の分布幅とさせることができる。
【００１６】
前記第２の推定手段には、前記第１の分布幅をｄ３、前記第２の特徴量および前記第３の特徴量をｘ・ｙ、前記第２の平均値をｍ１、前記第３の平均値をｍ２、前記第２の特徴量の数をｎ１、前記第３の特徴量の数をｎ２とした場合に、次式
【数１】
0004462176.tif 000002
を計算することにより求められる値を前記第１の分布幅とさせることができる。
【００１７】
画像処理装置には、前記第２の期間に撮像された前記被写体の画像の画像データから前記第２の特徴量を抽出し、前記第３の期間に撮像された前記被写体の画像の画像データから前記第３の特徴量を抽出する第２の抽出手段をさらに設けることができる。
【００１８】
画像処理装置には、特徴空間上における複数の前記第２の特徴量の分布と、複数の前記第３の特徴量の分布とを分離するための判別軸を算出する算出手段をさらに設け、前記第１の抽出手段には、前記第１の期間に撮像された画像の画像データから抽出された特徴量と、前記判別軸との内積の値を前記第１の特徴量とさせ、前記第２の抽出手段には、前記第２の期間に撮像された画像の画像データから抽出された特徴量と、前記判別軸との内積の値を前記第２の特徴量とさせ、前記第３の期間に撮像された画像の画像データから抽出された特徴量と、前記判別軸との内積の値を前記第３の特徴量とさせることができる。
【００１９】
本発明の一側面の画像処理方法またはプログラムは、第１の期間に撮像された被写体の画像の画像データから第１の特徴量を抽出し、前記第１の特徴量の第１の平均値を、前記第１の期間よりも過去の第２の期間に撮像された前記被写体の画像の画像データから抽出された第２の特徴量の第２の平均値と、前記第１の期間よりも過去の期間であって、前記第２の期間とは異なる第３の期間に撮像された前記被写体の画像の画像データから抽出された第３の特徴量の第３の平均値とに基づいて推定し、前記第１の平均値に対する前記第１の特徴量のばらつきの幅を示す第１の分布幅を、前記第２の平均値に対する前記第２の特徴量のばらつきの幅を示す第２の分布幅、および前記第３の平均値に対する前記第３の特徴量のばらつきの幅を示す第３の分布幅に基づいて推定し、前記第１の特徴量の値が、前記第１の平均値および前記第１の分布幅に基づく所定の範囲内の値であるか否かを判定するステップを含む。
【００２０】
本発明の一側面においては、第１の期間に撮像された被写体の画像の画像データから第１の特徴量が抽出され、前記第１の特徴量の第１の平均値が、前記第１の期間よりも過去の第２の期間に撮像された前記被写体の画像の画像データから抽出された第２の特徴量の第２の平均値と、前記第１の期間よりも過去の期間であって、前記第２の期間とは異なる第３の期間に撮像された前記被写体の画像の画像データから抽出された第３の特徴量の第３の平均値とに基づいて推定され、前記第１の平均値に対する前記第１の特徴量のばらつきの幅を示す第１の分布幅が、前記第２の平均値に対する前記第２の特徴量のばらつきの幅を示す第２の分布幅、および前記第３の平均値に対する前記第３の特徴量のばらつきの幅を示す第３の分布幅に基づいて推定され、前記第１の特徴量の値が、前記第１の平均値および前記第１の分布幅に基づく所定の範囲内の値であるか否かが判定される。
【発明の効果】
【００２１】
本発明の一側面によれば、顔画像認証を行うことができる。特に、本発明の一側面によれば、より簡単かつ確実に顔画像認証を行うことができる。
【発明を実施するための最良の形態】
【００２２】
以下に本発明の実施の形態を説明するが、本発明の構成要件と、明細書又は図面に記載の実施の形態との対応関係を例示すると、次のようになる。この記載は、本発明をサポートする実施の形態が、明細書又は図面に記載されていることを確認するためのものである。従って、明細書又は図面中には記載されているが、本発明の構成要件に対応する実施の形態として、ここには記載されていない実施の形態があったとしても、そのことは、その実施の形態が、その構成要件に対応するものではないことを意味するものではない。逆に、実施の形態が構成要件に対応するものとしてここに記載されていたとしても、そのことは、その実施の形態が、その構成要件以外の構成要件には対応しないものであることを意味するものでもない。
【００２３】
本発明の一側面の画像処理装置は、第１の期間（例えば、図３の期間Ｔ３）に撮像された被写体の画像の画像データから第１の特徴量を抽出する第１の抽出手段（例えば、図２の抽出部８２および図９の変換部２１７）と、前記第１の特徴量の第１の平均値を、前記第１の期間よりも過去の第２の期間（例えば、図３の期間Ｔ１）に撮像された前記被写体の画像の画像データから抽出された第２の特徴量の第２の平均値と、前記第１の期間よりも過去の期間であって、前記第２の期間とは異なる第３の期間（例えば、図３の期間Ｔ２）に撮像された前記被写体の画像の画像データから抽出された第３の特徴量の第３の平均値とに基づいて推定する第１の推定手段（例えば、図９の分布中心推定部２１４）と、前記第１の平均値に対する前記第１の特徴量のばらつきの幅を示す第１の分布幅を、前記第２の平均値に対する前記第２の特徴量のばらつきの幅を示す第２の分布幅、および前記第３の平均値に対する前記第３の特徴量のばらつきの幅を示す第３の分布幅に基づいて推定する第２の推定手段（例えば、図９の分布範囲推定部２１６）と、前記第１の特徴量の値が、前記第１の平均値および前記第１の分布幅に基づく所定の範囲内の値であるか否かを判定する判定手段（例えば、図９の判定部２１８）とを備える。
【００２５】
画像処理装置には、前記第２の期間に撮像された前記被写体の画像の画像データから前記第２の特徴量を抽出し、前記第３の期間に撮像された前記被写体の画像の画像データから前記第３の特徴量を抽出する第２の抽出手段（例えば、図２の抽出部８４および図９の変換部２１３）をさらに設けることができる。
【００２６】
画像処理装置には、特徴空間上における複数の前記第２の特徴量の分布と、複数の前記第３の特徴量の分布とを分離するための判別軸を算出する算出手段（例えば、図９の判別分析部２１２）をさらに設け、前記第１の抽出手段（例えば、図１０のステップＳ１５の処理を実行する変換部２１７）には、前記第１の期間に撮像された画像の画像データから抽出された特徴量と、前記判別軸との内積の値を前記第１の特徴量とさせ、前記第２の抽出手段（例えば、図１４のステップＳ５３の処理を実行する変換部２１３）には、前記第２の期間に撮像された画像の画像データから抽出された特徴量と、前記判別軸との内積の値を前記第２の特徴量とさせ、前記第３の期間に撮像された画像の画像データから抽出された特徴量と、前記判別軸との内積の値を前記第３の特徴量とさせることができる。
【００２７】
本発明の一側面の画像処理方法またはプログラムは、第１の期間に撮像された被写体の画像の画像データから第１の特徴量を抽出し（例えば、図１０のステップＳ１２およびステップＳ１５）、前記第１の特徴量の第１の平均値を、前記第１の期間よりも過去の第２の期間に撮像された前記被写体の画像の画像データから抽出された第２の特徴量の第２の平均値と、前記第１の期間よりも過去の期間であって、前記第２の期間とは異なる第３の期間に撮像された前記被写体の画像の画像データから抽出された第３の特徴量の第３の平均値とに基づいて推定し（例えば、図１４のステップＳ５４）、前記第１の平均値に対する前記第１の特徴量のばらつきの幅を示す第１の分布幅を、前記第２の平均値に対する前記第２の特徴量のばらつきの幅を示す第２の分布幅、および前記第３の平均値に対する前記第３の特徴量のばらつきの幅を示す第３の分布幅に基づいて推定し（例えば、図１４のステップＳ５６）、前記第１の特徴量の値が、前記第１の平均値および前記第１の分布幅に基づく所定の範囲内の値であるか否かを判定する（例えば、図１０のステップＳ１６）ステップを含む。
【００２８】
以下、図面を参照して本発明を適用した実施の形態について説明する。
【００２９】
図２は、本発明を適用した認証装置の構成例を示すブロック図である。
【００３０】
認証装置７１は、人の顔を被写体として撮像し、撮像された被写体としての人が、所定の人物であるか否かの顔画像認証を行う。
【００３１】
認証装置７１は、撮像部８１、抽出部８２、顔画像メモリ８３、抽出部８４、推定部８５、認証部８６、および出力部８７から構成される。
【００３２】
撮像部８１は、例えば、カメラなどからなり、認証の対象となるユーザ（被認証者）の顔の画像（以下、顔画像とも称する）を撮像する。撮像部８１は、撮像により得られた顔画像データを抽出部８２に供給する。
【００３３】
抽出部８２は、撮像部８１から供給された顔画像データから、顔画像データに基づく顔画像の特徴の特徴量を抽出し、認証部８６に供給する。
【００３４】
顔画像メモリ８３は、例えば、半導体メモリなどの不揮発性メモリなどからなり、過去に撮像された、被認証者としてのユーザの顔画像の顔画像データを記録する。顔画像メモリ８３は、記録している顔画像データを抽出部８４に供給する。
【００３５】
抽出部８４は、顔画像メモリ８３から供給された顔画像データから、顔画像データに基づく顔画像の特徴の特徴量を抽出し、推定部８５に供給する。
【００３６】
推定部８５は、抽出部８４から供給された特徴量に基づいて、撮像部８１において新たに撮像された被認証者の顔画像の特徴量が存在すると推定される、特徴空間上の所定の空間の中心の位置（以下、推定中心位置と称する）、およびその空間の端から推定中心位置までの長さ（以下、推定幅と称する）を求める。推定部８５は、求めた推定中心位置および推定幅を認証部８６に供給する。
【００３７】
認証部８６は、推定部８５から供給された推定中心位置および推定幅により定まる特徴空間上の所定の空間内に、抽出部８２から供給された特徴量が位置するか否かを判定し、その判定結果を示す信号を出力部８７に供給する。
【００３８】
出力部８７は、認証部８６から供給された、判定結果を示す信号に基づいて、認証装置７１に接続されている図示せぬコンピュータ等の装置に、本人であると認証された旨の信号、または本人ではないと認証された旨の信号を出力する。
【００３９】
ところで、認証装置７１には、過去に撮像された被認証者の顔画像の顔画像データが予め記録されている。例えば、認証装置７１には、図３に示すように、連続する期間Ｔ１乃至期間Ｔ３のうち、過去の期間Ｔ１および期間Ｔ２に撮像された顔画像の顔画像データが記録されている。
【００４０】
なお、図において、横軸は時間を示し、時刻ｔ１は、現在時刻を示している。したがって、現在時刻ｔ１から見て、期間Ｔ１および期間Ｔ２は過去の期間であり、期間Ｔ３（但し、時刻ｔ１も期間Ｔ３に含まれるものとする）は未来の期間である。
【００４１】
図３では、期間Ｔ１乃至期間Ｔ３のそれぞれの期間の長さは、同じ長さとなっている。また、期間Ｔ１乃至期間Ｔ３のそれぞれは、例えば、一昨日、昨日、および今日、先々週、先週、および今週、または一昨年、昨年、および今年などの期間とされる。
【００４２】
認証装置７１は、記録している期間Ｔ１および期間Ｔ２に撮像された顔画像の顔画像データを用いて、時刻ｔ１において撮像された顔画像の被認証者の認証を行う。
【００４３】
一般に、過去の顔画像の顔画像データを用いて、被認証者の認証を行う場合、例えば、図４Ａに示すように、矢印Ａ１１により示される期間Ｔ１に撮像された複数の顔画像の特徴量の分布と、矢印Ａ１２により示される期間Ｔ２に撮像された複数の顔画像の特徴量の分布とを用いて、矢印Ａ１３により示される期間Ｔ３に撮像される顔画像の特徴量が分布すると推定される範囲を求める。
【００４４】
なお、図中、横軸は時間を示しており、矢印Ａ１１および矢印Ａ１２により示される特徴空間上の１つの点は、１つの顔画像データから抽出された特徴量を示している。
【００４５】
従来の技術では、期間Ｔ１に撮像された顔画像と、期間Ｔ２に撮像された顔画像とは特に区別されず、図４Ｂに示すように、特徴空間における、期間Ｔ３の顔画像の特徴量が分布すると推定される範囲（空間）１１１は、期間Ｔ１および期間Ｔ２に撮像された顔画像の特徴量の中心の位置に基づいて定められる。
【００４６】
すなわち、従来の技術においては、時刻ｔ１（図３）において撮像された顔画像の特徴量が、特徴空間における、期間Ｔ１および期間Ｔ２に撮像された顔画像の特徴量からなる空間１１１内に位置していれば、撮像された顔画像の被認証者は、本人であると認証されていた。
【００４７】
ところで、人間の顔は、人間が年をとるとともに変化するので、顔画像の特徴には、被認証者ごとに固有であり、かつ時間が経過してもその特徴量が一定である特徴と、被認証者ごとに固有ではあるが、時間の経過とともにその特徴量が変化する特徴とがある。
【００４８】
したがって、過去の顔画像の特徴量の中心に基づいて空間１１１を定める方法においては、顔画像データから抽出する特徴量が、時間の経過とともに変化してしまう場合、被認証者が本人であったとしても、抽出された特徴量が空間１１１内に位置せず、本人ではないと認証されてしまう恐れがあった。
【００４９】
そこで、認証装置７１においては、時間の経過とともに特徴量が変化する特徴に注目し、図４Ｃに示すように、期間Ｔ１に撮像された顔画像の特徴量が分布する範囲（空間）１２１と、期間Ｔ２に撮像された顔画像の特徴量が分布する範囲（空間）１２２とから、期間Ｔ３に撮像される顔画像の特徴量が分布する範囲（空間）１２３を推定する。
【００５０】
図４Ｃでは、空間１２１と空間１２２とは、交わっておらず、顔画像の特徴量が分布する範囲（空間）は、時間の経過とともに変化していることが分かる。そこで、認証装置７１においては、期間Ｔ１に撮像された顔画像の特徴量（以下、特徴量ｘ１とも称する）と、期間Ｔ２に撮像された顔画像の特徴量（以下、特徴量ｘ２とも称する）とを区別して、それらの特徴量に基づいて空間１２３の位置および大きさを推定する。
【００５１】
より具体的には、例えば、図５に示すように、期間Ｔ１および期間Ｔ２に撮像された複数の顔画像の顔画像データから特徴量を抽出した結果、それぞれ矢印Ａ３１および矢印Ａ３２に示す特徴量の分布が得られたとする。ここで、矢印Ａ３１および矢印Ａ３２により示される特徴空間上の１つの点は、１つの顔画像データから抽出された特徴量を示している。
【００５２】
認証装置７１は、過去の顔画像（期間Ｔ１および期間Ｔ２に撮像された顔画像）の顔画像データから、特徴量を抽出すると、矢印Ａ３３に示すように、例えば、線形判別法（フィッシャーの方法）を用いて、特徴量ｘ１の分布と、特徴量ｘ２の分布とを分離するための判別軸１４１（図５では、２次元のベクトル）を求める。
【００５３】
また、矢印Ａ３３に示される特徴空間においては、範囲（空間）１４２は、特徴量ｘ１が分布する範囲（空間）を示しており、範囲（空間）１４３は、特徴量ｘ２が分布する範囲（空間）を示している。
【００５４】
認証装置７１は、判別軸１４１を求めると、期間Ｔ３に撮像された顔画像の特徴量（以下、特徴量ｘ３とも称する）は、特徴量ｘ１および特徴量ｘ２と同程度の広がり（ばらつき）を持つものとして、図６に示す期間Ｔ３に撮像された顔画像の特徴量が分布すると推定される範囲（空間）１６１を定める。
【００５５】
すなわち、特徴量ｘ３の重心（中心）の位置から空間１６１の端までの距離は、特徴量ｘ１の重心（中心）の位置から空間１４２の端までの距離、および特徴量ｘ２の重心（中心）の位置から空間１４３の端までの距離と同程度の距離であるとして、認証装置７１は、特徴量ｘ３が分布すると推定される空間１６１を定める。
【００５６】
ここで、特徴量ｘ１乃至特徴量ｘ３のそれぞれは、ｋ次元（但し、１≦ｋ）（図６では２次元）のベクトルであるとすると、認証装置７１は、例えば、図７に示すように、特徴量ｘ１および特徴量ｘ２のそれぞれを、判別軸１４１上に正射影して１次元の特徴量に変換する（特徴量を縮退する）。例えば、認証装置７１は、特徴量ｘ１および特徴量ｘ２のそれぞれと、判別軸１４１との内積を求めることにより、ｋ次元の特徴量を１次元の特徴量に変換する。なお、図中、特徴空間上の１つの点は、１つの特徴量（特徴量ｘ１または特徴量ｘ２）を示している。
【００５７】
これにより、図８に示すように、ｋ次元の特徴空間が、１次元の特徴空間に変換される。空間１４２（図６）は、ｍ１を中心とする幅がｄ１である１次元の空間（範囲）に変換され、空間１４３（図６）は、ｍ２を中心とする幅がｄ２である１次元の空間（範囲）に変換される。
【００５８】
なお、図中、横軸は特徴量、すなわち、ｋ次元の特徴量と判別軸１４１との内積により求まる１次元の特徴量を示し、縦方向は、その特徴量を持つ顔画像の数を示している。また、曲線１８１は、期間Ｔ１に撮像された顔画像の１次元の特徴量の分布を示しており、曲線１８２は、期間Ｔ２に撮像された顔画像の１次元の特徴量の分布を示している。
【００５９】
認証装置７１は、曲線１８１の中心の位置ｍ１と、曲線１８２の中心の位置ｍ２とから期間Ｔ３に撮像される顔画像の１次元の特徴量の分布の中心の位置ｍ３（推定中心位置）を推定（予測）し、さらに、ｄ１およびｄ２に基づいて、期間Ｔ３に撮像される顔画像の特徴量の分布の推定中心位置ｍ３からの偏り（ばらつき）の幅ｄ３（推定幅）を推定する。そして、認証装置７１は、推定中心位置ｍ３を中心とする幅がｄ３である範囲（空間）を、期間Ｔ３に撮像される顔画像の１次元の特徴量が分布すると推定される範囲とする。
【００６０】
したがって、認証処理装置７１は、顔画像認証時に被認証者の顔画像を撮像すると、その顔画像の１次元の特徴量を求め、その特徴量の値が、ｍ３を中心とする幅がｄ３である範囲内の値であるか否かに基づいて、撮像された顔画像の被認証者が本人であるか否かの認証を行う。
【００６１】
図９は、図２の推定部８５および認証部８６のより詳細な構成例を示すブロック図である。
【００６２】
推定部８５は、分類部２１１、判別分析部２１２、変換部２１３、分布中心推定部２１４、分布偏り算出部２１５、および分布範囲推定部２１６から構成される。また、認証部８６は、変換部２１７および判定部２１８から構成される。
【００６３】
分類部２１１は、抽出部８４（図２）から供給されたｋ次元（但し、１≦ｋ）の特徴量を、顔画像が撮像された期間に基づくクラスに分類する。例えば、分類部２１１は、抽出部８４から供給された特徴量のうち、所定の期間内に撮像された顔画像の特徴量のそれぞれを同じクラスに分類する。
【００６４】
また、分類部２１１は、特徴量をクラスに分類すると、抽出部８４から供給された特徴量のそれぞれが属すクラスを示すクラス情報を生成し、生成したクラス情報および特徴量を、判別分析部２１２および変換部２１３に供給する。
【００６５】
判別分析部２１２は、分類部２１１から供給されたクラス情報および特徴量に基づいて判別軸を求め、求められた判別軸（を示す情報）を変換部２１３および認証部８６の変換部２１７に供給する。
【００６６】
変換部２１３は、判別分析部２１２から供給された判別軸を用いて、分類部２１１から供給されたｋ次元の特徴量を、１次元の特徴量に変換する。また、変換部２１３は、変換により得られた１次元の特徴量および分類部２１１から供給されたクラス情報を、分布中心推定部２１４および分布偏り算出部２１５に供給する。
【００６７】
分布中心推定部２１４は、変換部２１３から供給された特徴量およびクラス情報に基づいて、未来（例えば、期間Ｔ３）に撮像される顔画像の特徴量の分布の中心の位置である、推定中心位置を算出し、算出された推定中心位置（を示す情報）を、分布範囲推定部２１６および認証部８６の判定部２１８に供給する。
【００６８】
分布偏り算出部２１５は、変換部２１３から供給された特徴量およびクラス情報に基づいて、各クラスに属す特徴量の分布の偏り（ばらつき）の幅である分布幅を算出し、算出された分布幅（を示す情報）を分布範囲推定部２１６に供給する。
【００６９】
分布範囲推定部２１６は、分布中心推定部２１４から供給された推定中心位置、および分布偏り算出部２１５から供給された分布幅に基づいて、未来（例えば、期間Ｔ３）に撮像される顔画像の特徴量の分布の偏り（ばらつき）の幅である推定幅を算出し、算出された推定幅（を示す情報）を認証部８６の判定部２１８に供給する。
【００７０】
認証部８６の変換部２１７は、判別分析部２１２から供給された判別軸を用いて、抽出部８２（図２）から供給されたｋ次元の特徴量を、１次元の特徴量に変換する。また、変換部２１７は、変換により得られた１次元の特徴量を判定部２１８に供給する。
【００７１】
判定部２１８は、変換部２１７から供給された特徴量、分布中心推定部２１４から供給された推定中心位置、および分布範囲推定部２１６から供給された推定幅に基づいて、推定中心位置および推定幅により定まる特徴空間上の空間（範囲）（以下、推定分布範囲とも称する）内に、撮像部８１において撮像された顔画像の特徴量が位置するか否かを判定し、その判定結果を示す信号を出力部８７（図２）に供給する。
【００７２】
次に、図１０のフローチャートを参照して、認証装置７１による認証処理について説明する。
【００７３】
ステップＳ１１において、撮像部８１は、被写体としての被認証者の顔画像を撮像し、これにより得られた顔画像データを抽出部８２に供給する。
【００７４】
ステップＳ１２において、抽出部８２は、撮像部８１から供給された、撮像により得られた顔画像データからｋ次元の特徴量を抽出し、抽出された特徴量を変換部２１７に供給する。
【００７５】
例えば、抽出部８２は、図１１に示すように、PCM（Pulse Code Modulation）方式により、撮像部８１から供給された顔画像２５１の顔画像データをサンプリングし、これより得られたｋ個（ｎ×ｍ個）の画素の画素値のそれぞれを成分とするｋ次元のベクトルを特徴量とする。
【００７６】
図１１において、図中、右側の顔画像２５１上の１つの丸は、１つの画素を表わしており、矢印Ｐ０乃至矢印Ｐ２のそれぞれは、図中、一番上の左から１番目乃至３番目の画素のそれぞれを示している。
【００７７】
サンプリングにより得られた顔画像２５１上の被認証者の顔の部分の画素の数が、横方向にｎ個（図１１では６個）であり、縦方向にｍ個（図１１では６個）であるとすると、それぞれの画素の画素値をベクトルの成分とする、（ｎ×ｍ）次元のベクトルが得られる。例えば、矢印Ｐ０乃至矢印Ｐ２のそれぞれにより示される画素の画素値のそれぞれを、Ｐ０乃至Ｐ２とすると、Ｐ０乃至Ｐ２のそれぞれが、特徴量としてのベクトルの成分とされる。
【００７８】
また、例えば、図１２に示すように、PCM方式により、撮像部８１から供給された顔画像２５１の顔画像データをサンプリングして得られたｋ個（ｎ×ｎ個）の画素の画素値（顔画像データ）を、さらにDCT（Discrete Cosine Transformation）変換し、これにより得られたｋ個のDCT係数のそれぞれを成分とするｋ次元のベクトルを特徴量としてもよい。
【００７９】
図１２において、図中、中央の顔画像２５１上の１つの丸は、１つの画素を表わしている。サンプリングにより得られた顔画像２５１上の被認証者の顔の部分の画素の数が、横方向にｎ個（図１２では６個）であり、縦方向にｎ個（図１２では６個）であるとすると、サンプリングを行うことにより、（ｎ×ｎ）個の画素値が得られる。さらに、これらの画素の画素値をDCT変換すると、それぞれの画素に対応するDCT係数が得られるので、これらのDCT係数をベクトルの成分とする（ｎ×ｎ）次元のベクトルが得られる。
【００８０】
さらに、図１３に示すように、顔画像２５１上の顔の特徴点を検出し、検出された特徴点の座標の値を成分とするベクトルを特徴量としてもよい。
【００８１】
図１３において、顔画像２５１上の点Ｂ１乃至点Ｂ８のそれぞれは、検出された被認証者の顔の特徴点を示している。例えば、点Ｂ１および点Ｂ２は、被認証者の右目の右端および左端を示しており、点Ｂ３および点Ｂ４は、被認証者の左目の右端および左端を示している。また、点Ｂ５および点Ｂ６は、被認証者の鼻の上端（付け根）および下端（鼻の先）を示しており、点Ｂ７および点Ｂ８は、認証者の口の右端および左端を示している。
【００８２】
ここで、顔画像２５１上の所定の点を基準とする点Ｂｎ（但し、１≦ｎ≦８）の座標を（Ｘｎ，Ｙｎ）とすると、それぞれの点Ｂｎの座標の成分ＸｎおよびＹｎをベクトルの成分とする２ｎ次元（図１３では、１６次元）のベクトルが得られる。したがって、図１３の例においては、点Ｂ１乃至点Ｂ８のそれぞれの座標の成分からなるベクトル（Ｘ１，Ｙ１，Ｘ２，Ｙ２，・・・，Ｘ８，Ｙ８）が特徴量とされる。
【００８３】
図１０のフローチャートの説明に戻り、ステップＳ１２において、撮像により得られた顔画像データからｋ次元の特徴量が抽出されると、ステップＳ１３に進む。
【００８４】
ステップＳ１３において、抽出部８４は、顔画像メモリ８３に記録されている、過去の顔画像データからｋ次元の特徴量を抽出し、抽出された特徴量を推定部８５の分類部２１１に供給する。例えば、抽出部８４は、顔画像メモリ８３から顔画像データを取得して、抽出部８２がステップＳ１２において行う処理と同様の処理を行い、取得した顔画像データから特徴量を抽出する。
【００８５】
ステップＳ１４において、推定部８５は、分布範囲の推定処理を行う。なお、分布範囲の推定処理の詳細は後述するが、分布範囲の推定処理において、推定部８５は、クラスごとの特徴量の分布を分離するための判別軸を求め、認証部８６の変換部２１７に供給する。
【００８６】
また、分布範囲の推定処理において、推定部８５は、クラスごとの分布の中心の位置（分布中心位置）および分布幅を求め、求められた分布中心位置および分布幅に基づいて、未来（例えば、期間Ｔ３）に撮像される顔画像の特徴量の分布の中心である推定中心位置、およびその分布の偏りの幅である推定幅を算出し、算出された推定中心位置および推定幅を認証部８６の判定部２１８に供給する。
【００８７】
ステップＳ１５において、変換部２１７は、推定部８５（の判別分析部２１２）から供給された判別軸を用いて、抽出部８２から供給された、撮像により得られた顔画像のｋ次元の特徴量を１次元の特徴量に変換する。変換部２１７は、変換により得られた１次元の特徴量を判定部２１８に供給する。
【００８８】
例えば、撮像により得られた顔画像のｋ次元の特徴量（ベクトル）をｘとし、判別軸（ベクトル）をｙとすると、変換部２１７は、ｋ次元の特徴量ｘと、ｋ次元の判別軸ｙとの内積ｘ・ｙを求めることにより、ｋ次元の特徴量を１次元の特徴量に変換する。
【００８９】
ステップＳ１６において、判定部２１８は、変換部２１７から供給された１次元の特徴量ｘ・ｙの値が、推定部８５の分布中心推定部２１４から供給された推定中心位置、および推定部８５の分布範囲推定部２１６から供給された推定幅により定まる特徴空間上の推定分布範囲内の値であるか否かを判定する。そして、判定部２１８は、その判定の結果を示す信号を生成して出力部８７に供給する。
【００９０】
例えば、変換部２１７からの１次元の特徴量をｘ・ｙ、推定中心位置をｍ３、推定幅をｄ３とすると、判定部２１８は、１次元の特徴量ｘ・ｙが、式（１）に示される関係を満たすか否かを基に、特徴量ｘ・ｙが推定分布範囲内の値であるか否かを判定する。
【００９１】
ｍ３−ｄ３＜ｘ・ｙ≦ｍ３＋ｄ３                              ・・・（１）
【００９２】
例えば、判定部２１８は、１次元の特徴量が式（１）に示される関係を満たす場合、推定分布範囲内の値であると判定し、１次元の特徴量が式（１）に示される関係を満たさない場合、推定分布範囲内の値ではないと判定する。
【００９３】
ステップＳ１６において、特徴量が推定分布範囲内の値であると判定された場合、ステップＳ１７に進み、出力部８７は、判定部２１８から供給された判定の結果を示す信号に基づいて、被認証者が本人であると認証された旨の信号を生成する。出力部８７は、生成した本人であると認証された旨の信号を図示せぬコンピュータなどに出力し、認証処理は終了する。
【００９４】
一方、ステップＳ１６において、特徴量が推定分布範囲内の値でないと判定された場合、ステップＳ１８に進み、出力部８７は、判定部２１８から供給された判定の結果を示す信号に基づいて、被認証者が本人ではないと認証された旨の信号を生成する。出力部８７は、生成した本人ではないと認証された旨の信号を図示せぬコンピュータなどに出力し、認証処理は終了する。
【００９５】
このようにして、認証装置７１は、顔画像が撮像された期間に基づくクラスごとの特徴量の分布中心位置および分布幅を求める。そして、認証装置７１は、クラスごとの分布中心位置および分布幅を基に推定分布範囲を求め、撮像により得られた顔画像データから抽出された特徴量が、推定分布範囲内の値であるか否かに基づいて、被認証者が本人であるか否かの認証を行う。
【００９６】
このように、クラスごとの分布中心位置および分布幅を基に、時間の経過により変化する推定分布範囲を求めることで、所定の期間に撮像された顔画像の特徴量が分布する特徴空間上の範囲をより簡単かつ正確に推定することができる。したがって、新たに撮像された顔画像の顔画像データから抽出した特徴量が、推定分布範囲内の値であるか否かに基づいて被認証者の個人認証を行うことで、より簡単かつ確実に顔画像認証を行うことができる。これにより、過去に撮像された被認証者の顔の写真を用いたなりすましを防止することができ、また、被認証者の顔が時間の経過とともに変化する場合においても、被認証者が本人ではないと認証されるようなこともない。
【００９７】
次に、図１４のフローチャートを参照して、図１０のステップＳ１４の処理に対応する、分布範囲の推定処理について説明する。
【００９８】
ステップＳ５１において、分類部２１１は、抽出部８４から供給されたｋ次元の特徴量を、顔画像が撮像された期間に基づくクラスに分類する。分類部２１１は、特徴量をクラスに分類するとその分類の結果を示すクラス情報を生成し、生成したクラス情報および特徴量を、判別分析部２１２および変換部２１３に供給する。
【００９９】
例えば、顔画像メモリ８３に、図３に示した期間Ｔ１および期間Ｔ２に撮像された顔画像の顔画像データが記録されている場合、分類部２１１は、抽出部８４から供給された特徴量のうち、期間Ｔ１に撮像された顔画像の特徴量のそれぞれを、クラスＴ１に属す特徴量として分類し、期間Ｔ２に撮像された顔画像の特徴量のそれぞれを、クラスＴ２に属す特徴量として分類する。
【０１００】
ステップＳ５２において、判別分析部２１２は、分類部２１１から供給されたクラス情報および特徴量に基づいて判別軸を求め、求められた判別軸を変換部２１３および変換部２１７に供給する。
【０１０１】
例えば、特徴量ｘが、クラスＴ１またはクラスＴ２のいずれかのクラスに分類される場合、判別分析部２１２は、クラスＴ１に属す特徴量ｘを基に、クラスＴ１の変動を表わす変動行列ＳＣ１を求め、同様に、クラスＴ２に属す特徴量ｘを基に、クラスＴ２の変動を表わす変動行列ＳＣ２を求める。そして、判別分析部２１２は、求めた変動行列ＳＣ１および変動行列ＳＣ２を基に、クラス内変動行列ＳＷと、クラス間変動行列ＳＢとを求める。
【０１０２】
さらに、判別分析部２１２は、クラス内変動行列ＳＷの逆行列ＳＷ-1と、クラス間変動行列ＳＢとの積ＳＷ-1ＳＢの最大固有値に対応する固有ベクトルを求め、求められた固有ベクトルを判別軸とする。
【０１０３】
ステップＳ５３において、変換部２１３は、判別分析部２１２から供給されたｋ次元の判別軸を用いて、分類部２１１から供給された、過去の顔画像のｋ次元の特徴量を、１次元の特徴量に変換する。変換部２１３は、変換により得られた１次元の特徴量、および分類部２１１から供給されたクラス情報を、分布中心推定部２１４および分布偏り算出部２１５に供給する。
【０１０４】
例えば、過去の顔画像のｋ次元の特徴量（ベクトル）をｘとし、ｋ次元の判別軸（ベクトル）をｙとすると、変換部２１３はｋ次元の特徴量ｘと、判別軸ｙとの内積ｘ・ｙを求めることにより、ｋ次元の特徴量を１次元の特徴量に変換する。
【０１０５】
ステップＳ５４において、分布中心推定部２１４は、変換部２１３から供給された特徴量およびクラス情報に基づいて推定中心位置を算出し、算出された推定中心位置を分布範囲推定部２１６および判定部２１８に供給する。
【０１０６】
例えば、クラスＴｉ（但し、ｉ＝１またはｉ＝２）に属す特徴量ｘの集合をｘｉとすると、分布中心推定部２１４は、式（２）を計算することにより、クラスＴｉに属す特徴量ｘ（すなわち、集合ｘｉの元ｘ）が縮退された１次元の特徴量ｘ・ｙの分布の中心の位置ｍｉを求める。
【０１０７】
【数２】
0004462176.tif 000003
・・・（２）
【０１０８】
ここで、ｎｉは、クラスＴｉに属す特徴量ｘの数（集合ｘｉの元ｘの数）を表わしている。また、Σｘ・ｙは、クラスＴｉに属す特徴量ｘについてｘ・ｙを求め、その総和をとることを表わしている。例えば、分布中心推定部２１４が式（２）を計算することにより、図８に示したクラスＴ１の分布の中心の位置（以下、分布中心位置とも称する）ｍ１、およびクラスＴ２の分布中心位置ｍ２が求まる。
【０１０９】
また、期間Ｔ１乃至期間Ｔ３（図３）のそれぞれの長さは等しいので、期間Ｔ１の所定の基準となる時刻（例えば、期間Ｔ１の先頭の時刻）をｔとすると、期間Ｔｉ（但し、１≦ｉ≦３）の基準となる時刻（例えば、期間Ｔｉの先頭の時刻）は、ｉｔで表わすことができる。
【０１１０】
さらに、期間Ｔ１の分布中心位置がｍ１であり、期間Ｔ２の分布中心位置がｍ２であるので、期間Ｔ１から期間Ｔ２までの期間で分布中心位置は、ｍ２−ｍ１だけ変化することが分かる。そして、所定の期間Ｔ（ｉ−１）の基準となる時刻から、期間Ｔｉの基準となる時刻までの期間の長さは常にｔであるので、期間Ｔ（ｉ−１）から期間Ｔｉまでの期間で分布中心位置は、ｍ２−ｍ１だけ変化すると推定することができ、期間Ｔｉの分布中心位置ｍｉは、式（３）により表わすことができる。
【０１１１】
ｍｉ＝（（ｍ２−ｍ１）／ｔ）（Ｔｉ）＋２（ｍ１）−ｍ２      ・・・（３）
【０１１２】
この式（３）の変数Ｔｉに、期間Ｔ３＝３ｔを代入することにより、期間Ｔ３の推定中心位置（分布中心位置）ｍ３の値２ｍ１−ｍ１（＝（ｍ２−ｍ１）／ｔ）（３ｔ）＋２（ｍ１）−ｍ２）を求めることができる。分布中心推定部２１４は、式（３）を計算することにより、推定中心位置ｍ３を算出すると、算出された推定中心位置ｍ３を、分布範囲推定部２１６および判定部２１８に供給し、処理はステップＳ５５に進む。
【０１１３】
ステップＳ５５において、分布偏り算出部２１５は、変換部２１３から供給された特徴量およびクラス情報に基づいて、各クラスに属す特徴量の分布の分布幅を算出し、算出された分布幅を分布範囲推定部２１６に供給する。
【０１１４】
例えば、特徴量ｘが、クラスＴ１またはクラスＴ２のいずれかのクラスに分類された場合、クラスＴｉ（但し、ｉ＝１またはｉ＝２）に属す特徴量ｘの集合をｘｉとすると、分布偏り算出部２１５は、上述した式（２）を計算することにより、クラスＴｉに属す特徴量ｘ（すなわち、集合ｘｉの元ｘ）が縮退された１次元の特徴量ｘ・ｙの分布中心位置ｍｉを求める。
【０１１５】
そして、分布偏り算出部２１５は、各クラスＴｉに属す（特徴量ｘが縮退された）特徴量ｘ・ｙと、そのクラスＴｉの分布中心位置ｍｉとの差の絶対値のうち、最も大きい値をクラスＴｉの分布幅ｄｉとする。すなわち、分布偏り算出部２１５は、式（４）を計算することによって、クラスＴｉの分布幅ｄｉを求める。
【０１１６】
ｄｉ＝ｍａｘ（｜ｘ・ｙ−ｍｉ｜：ｘ∈ｘｉ）                  ・・・（４）
【０１１７】
ここで、ｍａｘ（Ｑ）は、括弧内の変数Ｑのうちの最大の値を示す。したがって、式（４）は、ｘを集合ｘｉの元として、｜ｘ・ｙ−ｍｉ｜のうち、最も大きい値をｄｉとすることを示している。
【０１１８】
これにより、例えば、図８に示したクラスＴ１の分布幅ｄ１、およびクラスＴ２の分布幅ｄ２が求まる。
【０１１９】
ステップＳ５６において、分布範囲推定部２１６は、分布偏り算出部２１５から供給された分布幅ｄｉに基づいて推定幅を算出する。そして、分布範囲推定部２１６は、算出された推定幅を認証部８６の判定部２１８に供給して、処理は図１０のステップＳ１５に進む。
【０１２０】
例えば、分布範囲推定部２１６は、分布偏り算出部２１５から分布幅ｄｉとして、クラスＴ１の分布幅ｄ１およびクラスＴ２の分布幅ｄ２が供給された場合、供給された分布幅ｄ１および分布幅ｄ２のうち、大きい分布幅（最大の分布幅）の値を、未来の期間Ｔ３の特徴量の分布の偏りの幅である推定幅ｄ３とする。すなわち、分布範囲推定部２１６は、式（５）を計算することにより、推定幅ｄ３を算出する。
【０１２１】
ｄ３＝ｍａｘ（ｄ１，ｄ２）                                  ・・・（５）
【０１２２】
ここで、ｍａｘ（ｄ１，ｄ２）は、ｄ１およびｄ２のうち、大きい方の値をｄ３とすることを示している。分布範囲推定部２１６が式（５）を計算することにより、例えば、図８に示した推定幅ｄ３が求まる。したがって、図８の例では、未来の期間Ｔ３に撮像される顔画像の１次元の特徴量ｘ・ｙは、推定中心位置ｍ３を中心とする幅がｄ３である空間（推定分布範囲）内に存在すると推定される。
【０１２３】
このようにして、推定部８５は、過去の顔画像のｋ次元の特徴量を１次元の特徴量に変換し、変換により得られた１次元の特徴量を基に、推定中心位置および推定幅を算出する。
【０１２４】
このように、ｋ次元の特徴量を１次元の特徴量に変換することで、変換により得られた１次元の特徴量を基に、より簡単に推定中心位置および推定幅を算出することができる。
【０１２５】
なお、推定幅の求め方の一例として、分布範囲推定部２１６が式（５）を計算し、各クラスＴｉの分布幅のうち、最も大きい分布幅を推定幅とする例について説明したが、各クラスＴｉの特徴量ｘ・ｙの標準偏差を用いて推定幅を求めるようにしてもよい。
【０１２６】
例えば、特徴量ｘが、クラスＴ１またはクラスＴ２のいずれかのクラスに分類された場合、分布範囲推定部２１６は、式（６）を計算することにより、推定幅ｄ３を算出する。
【０１２７】
ｄ３＝（ｅ）1/2                                            ・・・（６）
【０１２８】
なお、式（６）におけるｅは、式（７）により表わされる。
【０１２９】
【数３】
0004462176.tif 000004
・・・（７）
【０１３０】
ここで、右側のΣは、クラスＴｉに属す特徴量ｘについて、（ｘ・ｙ−ｍｉ）2の総和をとることを表わしており、左側のΣは、Σ（ｘ・ｙ−ｍｉ）2のｉを１から２まで変えて総和をとることを表わしている。また、ｎ１およびｎ２は、それぞれクラスＴ１およびクラスＴ２に属す特徴量ｘの数を表わしている。
【０１３１】
このように、特徴量の標準偏差を用いて推定幅を求めることにより、分布幅の最大値を推定幅とする場合と比べて、より正確な（より厳しい）推定幅を求めることができ、被認証者以外の人物が本人であると認証される可能性を低減することができる。
【０１３２】
また、分布中心振推定部２１４において、推定中心位置を算出し、分布偏り算出部２１５において、分布幅を算出すると説明したが、認証処理を行うごとに算出せずに、分布中心振推定部２１４および分布偏り算出部２１５が、算出された推定中心位置および分布幅を記憶するようにし、認証処理を行う場合に、記憶している推定中心位置および分布幅を必要に応じて用いるようにしてもよい。
【０１３３】
さらに、新たに撮像された顔画像により被認証者が認証されるごとに、その新たに撮像された顔画像の顔画像データを顔画像メモリ８３に記録しておき、次回以降に行う認証処理時に、過去の顔画像データとして用いるようにしてもよい。
【０１３４】
さらに、また、過去の期間として、顔画像メモリ８３に記録される過去の顔画像が撮像される期間は、連続する２つの期間に限定されない。すなわち、過去の顔画像として用いられる顔画像が撮像される期間は、複数の不連続の期間であってもよく、連続する期間であってもよい。また、それらの期間は２以上の期間であればよい。
【０１３５】
以上のように、クラスごとの分布中心位置および分布幅を基に推定分布範囲を求め、撮像された顔画像の顔画像データから抽出した特徴量が、推定分布範囲内の値であるか否かに基づいて被認証者の認証を行うようにしたので、より簡単かつ確実に顔画像認証を行うことができる。
【０１３６】
図１５は、上述した一連の処理をプログラムにより実行するパーソナルコンピュータの構成の例を示すブロック図である。パーソナルコンピュータ３０１のCPU（Central Processing Unit）３１１は、ROM（Read Only Memory）３１２、または記録部３１８に記録されているプログラムに従って各種の処理を実行する。RAM（Random Access Memory）３１３には、CPU３１１が実行するプログラムやデータなどが適宜記憶される。これらのCPU３１１、ROM３１２、およびRAM３１３は、バス３１４により相互に接続されている。
【０１３７】
CPU３１１にはまた、バス３１４を介して入出力インターフェース３１５が接続されている。入出力インターフェース３１５には、キーボード、マウス、マイクロホンなどよりなる入力部３１６、ディスプレイ、スピーカなどよりなる出力部３１７が接続されている。CPU３１１は、入力部３１６から入力される指令に対応して各種の処理を実行する。そして、CPU３１１は、処理の結果を出力部３１７に出力する。
【０１３８】
入出力インターフェース３１５に接続されている記録部３１８は、例えばハードディスクからなり、CPU３１１が実行するプログラムや各種のデータを記録する。通信部３１９は、インターネットやローカルエリアネットワークなどのネットワークを介して外部の装置と通信する。
【０１３９】
また、通信部３１９を介してプログラムを取得し、記録部３１８に記録してもよい。
【０１４０】
入出力インターフェース３１５に接続されているドライブ３２０は、磁気ディスク、光ディスク、光磁気ディスク、或いは半導体メモリなどのリムーバブルメディア３３１が装着されたとき、それらを駆動し、そこに記録されているプログラムやデータなどを取得する。取得されたプログラムやデータは、必要に応じて記録部３１８に転送され、記録される。
【０１４１】
上述した一連の処理は、ハードウエアにより実行させることもできるし、ソフトウエアにより実行させることもできる。一連の処理をソフトウエアにより実行させる場合には、そのソフトウエアを構成するプログラムが、専用のハードウエアに組み込まれているコンピュータ、または、各種のプログラムをインストールすることで、各種の機能を実行することが可能な、例えば汎用のパーソナルコンピュータなどに、プログラム記録媒体からインストールされる。
【０１４２】
コンピュータにインストールされ、コンピュータによって実行可能な状態とされるプログラムを格納するプログラム記録媒体は、図１５に示すように、磁気ディスク（フレキシブルディスクを含む）、光ディスク（CD-ROM(Compact Disc-Read Only Memory),DVD(Digital Versatile Disc)を含む）、光磁気ディスクを含む）、もしくは半導体メモリなどよりなるパッケージメディアであるリムーバブルメディア３３１、または、プログラムが一時的もしくは永続的に格納されるROM３１２や、記録部３１８を構成するハードディスクなどにより構成される。プログラム記録媒体へのプログラムの格納は、必要に応じてルータ、モデムなどのインターフェースである通信部３１９を介して、ローカルエリアネットワーク、インターネット、デジタル衛星放送といった、有線または無線の通信媒体を利用して行われる。
【０１４３】
なお、本明細書において、プログラム記録媒体に格納されるプログラムを記述するステップは、記載された順序に沿って時系列的に行われる処理はもちろん、必ずしも時系列的に処理されなくとも、並列的あるいは個別に実行される処理をも含むものである。
【０１４４】
なお、本発明の実施の形態は、上述した実施の形態に限定されるものではなく、本発明の要旨を逸脱しない範囲において種々の変更が可能である。
【図面の簡単な説明】
【０１４５】
【図１】従来の顔画像認識を説明するための図である。
【図２】認証装置の構成例を示すブロック図である。
【図３】顔画像が撮像される期間を説明するための図である。
【図４】顔画像の特徴量の分布を示す図である。
【図５】顔画像の特徴量の分布を示す図である。
【図６】顔画像の特徴量が分布すると推定される範囲を説明するための図である。
【図７】ｋ次元の特徴量の１次元の特徴量への変換を説明するための図である。
【図８】顔画像の特徴量が分布すると推定される範囲を説明するための図である。
【図９】推定部および認証部のより詳細な構成例を示すブロック図である。
【図１０】認証処理を説明するフローチャートである。
【図１１】顔画像データから抽出される特徴量の例を説明するための図である。
【図１２】顔画像データから抽出される特徴量の例を説明するための図である。
【図１３】顔画像データから抽出される特徴量の例を説明するための図である。
【図１４】分布範囲の推定処理を説明するフローチャートである。
【図１５】パーソナルコンピュータの構成例を示すブロック図である。
【符号の説明】
【０１４６】
７１  認証装置，  ８１  撮像部，  ８２  抽出部，  ８３  顔画像メモリ，  ８４  抽出部，  ８５  推定部，  ８６  認証部，  ２１２  判別分析部，  ２１３  変換部，  ２１４  分布中心推定部，  ２１５  分布偏り算出部，  ２１６  分布範囲推定部，  ２１７  変換部，  ２１８  判定部，  ３０１  パーソナルコンピュータ，  ３１１  CPU，  ３１２  ROM，  ３１３  RAM，  ３１８  記録部，  ３３１  リムーバブルメディア
【図１】
0004462176.tif 000006
【図２】
0004462176.tif 000007
【図３】
0004462176.tif 000008
【図４】
0004462176.tif 000009
【図５】
0004462176.tif 000010
【図６】
0004462176.tif 000011
【図７】
0004462176.tif 000012
【図８】
0004462176.tif 000013
【図９】
0004462176.tif 000014
【図１０】
0004462176.tif 000015
【図１１】
0004462176.tif 000016
【図１２】
0004462176.tif 000017
【図１３】
0004462176.tif 000018
【図１４】
0004462176.tif 000019
【図１５】
0004462176.tif 000020