JP3846851B2 - 画像のマッチング処理方法及びその装置 - Google Patents

Description

【０００１】
【発明の属する技術分野】
本発明は、テンプレート画像を用いて、入力画像中から対象物を検出する画像のマッチング処理方法及びその装置に関するものである。
【０００２】
【従来の技術】
従来より、テンプレート画像を予め登録し、入力画像とテンプレート画像との間で、パターンマッチングを行って、テンプレート画像に類似する画像の、入力画像における位置を検出する技術は、周知である。
【０００３】
しかし、テンプレート画像に類似する画像の背景如何では、誤認識が発生しやすいため、この点を工夫した技術が、特開平５−２８２７３号公報として、開示されている。
【０００４】
このものでは、テンプレート画像とテンプレート画像に対応する画像との、類似値を、次式で定義する。
【数１】

【０００５】
つまり、テンプレート画像のエッジ法線ベクトルと、入力画像のエッジ法線ベクトルとのなす角Θの、内積（ｃｏｓΘ）を類似値の成分としている。
【０００６】
【発明が解決しようとする課題】
しかしながら、後に詳しく説明するように、対象物の像の周囲で、背景部分の輝度がばらつくと、内積の正負が逆転して、類似値が実体に合わなくなり、誤認識が発生しやすく、良好な認識結果が得にくいという問題点がある。
【０００７】
また、入力画像及びテンプレート画像のエッジ法線方向ベクトルに対して、類似値式が、非線形であり、テンプレート画像に係る処理と、入力画像のそれとを、同時に実行しなければならない。
【０００８】
そして、入力画像上にテンプレート画像を走査させ、各走査点毎に、入力画像と参照画像の相関計算を行わざるを得ず、計算量が膨大で、事実上、実時間処理は不可能である。
【０００９】
そこで本発明は、正確かつ先鋭な認識結果と処理の高速性とが得られる画像処理方法及びその装置を提供することを目的とする。
【００１０】
【課題を解決するための手段】
本発明においては、テンプレート画像と入力画像とのマッチングを、類似値のマップを用いて評価する。また、テンプレート画像と入力画像とのそれぞれについて、評価ベクトルを生成する。しかも、これらの評価ベクトルは、該当画像のエッジ法線方向ベクトルを、偶数倍角変換した成分を含む。
【００１１】
この構成により、認識結果が正確かつ先鋭で、処理が高速な画像処理方法が実現できる。
【００１２】
【発明の実施の形態】
請求項１記載の画像処理方法では、テンプレート画像と入力画像とのマッチングを、類似値のマップを用いて評価するものであり、テンプレート画像と入力画像とのそれぞれについて、評価ベクトルを生成し、しかも、これらの評価ベクトルは、該当画像のエッジ法線方向ベクトルを、偶数倍角変換した成分を含む。
【００１３】
これにより、背景部分の輝度のばらつきにより、テンプレート画像のエッジ法線ベクトルと、入力画像のエッジ法線ベクトルとのなす角Θの、内積（ｃｏｓΘ）の正負が反転する場合でも、類似値に影響が無く、正当にマッチングを評価できる。
【００１４】
請求項２記載の画像処理方法では、テンプレート画像及び入力画像の双方について、該当画像を入力し、該当画像のエッジ法線方向ベクトルを求めるステップと、エッジ法線方向ベクトルから評価ベクトルを生成するステップと、評価ベクトルを直交変換するステップとを備え、テンプレート画像及び入力画像のそれぞれについて得た、直交変換後の各評価ベクトルについて、対応スペクトルデータを積和計算するステップと、積和計算結果を、逆直交変換して類似値のマップを生成するステップとを有し、かつ、類似値の式、直交変換及び逆直交変換は、いずれも線形性をもつものである。
【００１５】
これにより、テンプレート画像に係るフーリエ変換値と、入力画像に係るフーリエ変換値は、同時に求める必要はない。つまり、テンプレート画像に係るフーリエ変換値を、入力画像のそれに先行して、求めておくことができ、処理負担を軽くして、一層高速化できる。
【００１６】
請求項３記載の画像処理方法では、直交変換後の各評価ベクトルを、圧縮して処理量を削減するステップを含む。
【００１７】
これにより、処理対象を、有効な成分（例えば、低周波成分）のみに削減して、さらに高速な処理を実現できる。
【００１８】
請求項４記載の画像処理方法では、テンプレート画像については、入力画像を入力するに先立ち、直交変換後の評価ベクトルを圧縮するまでのステップを行っておき、その結果を記録手段に保存しておく。
【００１９】
これにより、テンプレート画像側の処理は、記録手段から読み出すだけで良く、一層高速化できる。
【００２０】
請求項５記載の画像処理方法では、評価ベクトルは、ベクトルの長さについて正規化されている。
【００２１】
撮影条件により、入力画像のエッジの強度は変化し、ベクトルの長さが変化するが、このようにすることにより、長さの変化を受けないようにして、パターン抽出の安定性を高めることができる。
【００２２】
請求項６記載の画像処理方法では、テンプレート画像に係る評価ベクトルは、エッジ法線方向ベクトルの個数で正規化されている。
【００２３】
このように、ｎで割って正規化することにより、テンプレート画像における、エッジ個数の多寡によらず、同じ尺度で類似度合いを評価できる。
【００２４】
請求項７記載の画像処理方法では、積和計算を行う前に、直交変換の複素共役性を利用してデータ量を削減し、積和計算後にデータ量を復元する。
【００２５】
これにより、データ量を大幅に削減して、処理を高速化し、記憶容量を節約できる。
【００２６】
請求項８記載の画像処理方法では、テンプレート画像を各種サイズに拡縮し、各サイズについて評価ベクトルの加算処理を行う。
【００２７】
これにより、各種サイズ毎にマッチングを繰り返す必要がなくなり、高速化できる。
【００２８】
請求項９記載の画像処理方法では、テンプレート画像について、評価ベクトルの加算処理を、圧縮して処理量を削減するステップの後に行う。
【００２９】
これにより、加算処理対象を削減でき、一層高速化を図れる。
【００３０】
請求項１０記載の画像処理方法では、テンプレート画像は、顔を模式化したものである。
【００３１】
これにより、顔の全体的な位置だけでなく、目鼻口等の主要な顔部品の位置も把握できる。
【００３２】
請求項１１記載の画像処理方法では、類似値のピークを急峻にするピークパターンを用意し、このピークパターンのデータを直交変換したものを、積和計算に反映させる。
【００３３】
これにより、ピークパターンを反映して、より先鋭且つ安定に、テンプレートと類似する部分を入力画像中から検出できる。
【００３４】
請求項１２記載の画像処理方法では、テンプレート画像に依存するマスクパターンを作成し、このマスクパターンのデータを直交変換したものを、積和計算に反映させる。
【００３５】
これにより、テンプレート画像の形状以外の属性を加味して、より精密な検出ができる。
【００３６】
請求項１３記載の画像処理方法では、マスクパターンは、テンプレート画像における画像内部の画素数平均を示すものである。
【００３７】
これにより、簡単なマスクパターンで、テンプレート画像の属性を反映できる。
【００３８】
請求項１４記載の画像処理方法では、テンプレート画像について、元のテンプレート画像の評価ベクトルの正負符号を処理して、元のテンプレート画像に対し左右対称な画像の評価ベクトルを生成するステップを有し、生成した評価ベクトルを積和計算に反映させる。
【００３９】
これにより、テンプレート画像の記録量を節約でき、しかも、左右反転したテンプレート画像の評価ベクトルを、直接的な演算なしに、生成でき、高速化を図れる。
【００４０】
請求項１５記載の画像処理方法では、抽出した顔画像に基づいて、点双列相関係数のマップを生成し、顔部品の位置を求める。
【００４１】
これにより、顔部品の位置を、より正確に特定できる。
【００４２】
請求項１６記載の画像処理方法では、抽出した顔画像に基づき、マスクパターンを用いてｙ方向の投影値分布を求め、この分布から２つの極大点を求め、この極大点間を口範囲として出力する。
【００４３】
これにより、口範囲を、より正確に特定できる。
【００４４】
請求項１７記載の画像処理方法では、抽出した顔画像に基づいて、入力画像を、顔画像のみと、顔画像を除いた部分とに分離し、顔画像のみに電子透かしを埋込み、電子透かしを埋め込んだ顔画像と、顔画像を除いた部分とを、合成して出力する。
【００４５】
これにより、改ざんの対象になりやすい顔部分に重点的に透かしデータを埋め込める。
【００４６】
請求項１８記載の画像処理方法では、抽出した顔画像に基づいて、入力画像を、顔画像のみと、顔画像を除いた部分とに分離し、顔画像のみに編集を施し、編集した顔画像と、顔画像を除いた部分とを、合成して出力する。
【００４７】
これにより、顔以外の部分に影響を与えずに、顔画像のみ補正できる。
【００４８】
請求項３２記載の画像処理方法では、抽出した顔画像に基づいて、入力画像から顔画像を切り出し、切り出した顔画像から顔内部画像を抽出し、抽出した顔画像に基づいて顔画像の補正に役立つ特徴を求め、求めた特徴に基づいて補正関数を決定し、少なくとも切り出した顔画像について、決定した補正関数に基づく画像補正を施す。
【００４９】
これにより、顔画像の補正は、顔内部画像のみの特徴に依存し、顔画像でない画像の悪影響を受けないので、顔画像の見やすさを確実に向上できる。
【００５０】
（実施の形態１）
次に、本発明の画像処理装置の基本形について、図１から図８を参照しながら、説明する。
【００５１】
図１に示すように、実施の形態１に係る画像処理装置は、テンプレート画像処理部１００と、入力画像処理部２００との２つの系統の処理部を持ち、テンプレート画像と入力画像とのマッチングを類似値Ｌのマップを用いて評価する。また、この画像処理装置では、テンプレート画像処理部１００と入力画像処理部２００とにおいて、それぞれ、線形性を持つ直交変換を行い、積算した後に、逆直交変換して、類似値Ｌを求める。
【００５２】
ここで、以下すべての実施の形態において、このような直交変換として、ＦＦＴ（高速離散フーリエ変換）を使用することとするが、他に、Ｈａｒｔｌｅｙ変換や数論的変換などを用いることもでき、以下の説明において、「フーリエ変換」とあるのを、これらの変換に置き換えて適用することができる。
【００５３】
また、テンプレート画像処理部１００、入力画像処理部２００のいずれにおいても、エッジ法線方向ベクトルの内積を利用し、エッジ法線方向ベクトルの方向が近いほど、高い相関を出すようにしている。しかも、この内積は、偶数倍角表現を用いて評価される。以下簡単のため、偶数倍角の例として、２倍角の場合のみを説明するが、４倍角、６倍角などの、他の偶数倍角においても、本発明と同様の効果がある。
【００５４】
次に、テンプレート画像処理部１００について、説明する。エッジ抽出手段１は、テンプレート画像に対して、ｘ方向、ｙ方向のそれぞれについて、微分処理（エッジ抽出）を施し、テンプレート画像のエッジ法線方向ベクトルを出力する。
【００５５】
本例では、ｘ方向について、
【数２】

なるＳｏｂｅｌフィルタを用い、
ｙ方向について、
【数３】

なるＳｏｂｅｌフィルタを用いている。
【００５６】
これらのフィルタにより、次式で定義される、テンプレート画像のエッジ法線方向ベクトルが求められる。
【数４】

【００５７】
本例では、横断歩道付近を撮影した入力画像から、特定の姿勢をして横断歩道を歩いている人物の像を抽出するものとする。
【００５８】
このとき、人物のテンプレート画像は、例えば、図２（ａ）のような画像になる。また、図２（ａ）のテンプレート画像について、（数２）のフィルタ処理を施すと、図２（ｂ）のような結果（ｘ成分）が得られ、図２（ａ）のテンプレート画像について、（数３）のフィルタ処理を施すと、図２（ｃ）のような結果（ｙ成分）が得られる。
【００５９】
評価ベクトル生成手段２は、エッジ抽出手段１からテンプレート画像のエッジ法線方向ベクトルを、入力し、次に述べる処理を行って、テンプレート画像の評価ベクトルを直交変換手段３へ出力する。
【００６０】
まず、評価ベクトル生成手段２は、次式を用いて、テンプレート画像のエッジ法線方向ベクトルを、長さについて正規化する。
【数５】

【００６１】
ここで一般に、撮影条件により、入力画像のエッジの強度は、変化してしまう。しかし、入力画像とテンプレート画像との、それぞれのエッジ間の角度差（あるいは、この角度差により単調変化する従属関数値）は、撮影条件の影響を受けにくい。
【００６２】
そこで、本発明では、後述するように、入力画像処理部２００において、入力画像のエッジ法線方向ベクトルを、長さ１に正規化している。これにあわせて、テンプレート画像処理部１００においても、テンプレート画像のエッジ法線方向ベクトルも、長さ１に正規化している。
【００６３】
これにより、パターン抽出の安定性を高めることができる。なお通常、正規化する長さは、通常「１」が良いと考えられるが、他の定数を用いることもできる。
【００６４】
また、周知のように、三角関数については、次の倍角公式が成立する。
【数６】

【００６５】
そして、評価ベクトル生成手段２は、次式で定義される、テンプレート画像の評価ベクトルを求める。
【数７】

【００６６】
（数７）について説明する。まず、定数ａより小さなベクトルを、ゼロベクトルにしているのは、ノイズなどを除去するためである。
【００６７】
次に、この評価ベクトルのｘ，ｙ各成分を、ｎで割って正規化している点について、説明する。
【００６８】
一般に、テンプレート画像の形状は、任意であり、そのエッジの形状は様々である。例えば、図８（ａ）に示すように、エッジの個数が少ないこともあるし、図８（ｂ）のように、（図８（ａ）に比べれば）エッジの個数が多いこともある。そこで、本形態では、ｎで割って正規化することにより、テンプレート画像における、エッジ個数の多寡によらず、同じ尺度で類似度合いを評価できるように、工夫している。
【００６９】
但し、ｎで割る正規化処理は、必ず実施しなければならないものではなく、一種類のテンプレート画像のみしか用いない場合や、エッジの個数が同じテンプレート画像しか用いない場合では、ｎで割る正規化処理は、省略できる。
【００７０】
次に、（数７）のｘ，ｙ各成分が、（数５）のｘｙ各成分の、倍角に係る余弦・正弦の従属関数になっている点について、説明する。
【００７１】
ここで、従来の技術のように、テンプレートの評価ベクトルＴと、入力画像の評価ベクトルＩとのなす角を、Θとし、その内積、つまり、ｃｏｓΘを類似尺度として用いると、次のような問題がある。
【００７２】
例えば、テンプレート画像は図７（ａ）のとおりであり、入力画像は図７（ｂ）に示すようであるものとする。このうち、図７（ｂ）の背景部分について、対象物の像の左側部分が、対象物の像よりも明るく、対象物の像の右側部分が、対象物の像よりも暗くなっている。
【００７３】
画像のみで見れば、図７（ａ）のテンプレート画像の中心が、図７（ｂ）の入力画像の中心に一致するとき、対象物の像が完全に一致するので、このとき、類似値は、最大にならなければならない。そして、エッジ法線方向ベクトルは、対象部の像から外側に向くものを正とすると、図７（ｂ）の明るい背景部分でも、暗い背景部分でも、対象物の像から見て、同じ向き（外向き／内向き）を向かなければならない。
【００７４】
しかしながら、このとき、図７（ｂ）の背景部分の輝度が、対象物の左右でばらついていると、図７（ｂ）に矢印で示しているように、向きが、反対（明るい背景部分では対象物から外側向き、暗い背景部分では対象物の内側向き）になってしまう。
【００７５】
このような場合、本来、最大の類似値となるべき場合において、必ずしも類似値が高い値にならず、誤認識を招きやすい。
【００７６】
以上の点を、図６を用いて、さらに詳しく説明する。テンプレート画像の評価ベクトルＴと、入力画像の評価ベクトルＩとのなす角を、Θとし、その内積、つまり、ｃｏｓΘを類似値として用いる場合、以上説明したように、対象物の像の周囲に存在する、背景部分の輝度のばらつきにより、入力画像の評価ベクトルＩの方向には、図６のＩと、その正反対のＩ’の、２通りの可能性がある。
【００７７】
このため、類似尺度である、内積は、ｃｏｓΘと、ｃｏｓΘ’との、２通りがあり得ることになる。
【００７８】
しかも、Θ＋Θ’＝πであり、ｃｏｓΘ’＝ｃｏｓ（π−Θ）＝−ｃｏｓΘである。
【００７９】
つまり、ｃｏｓΘを類似尺度として用いると、本来、類似値を増やすように作用しなければならない場合において、逆に類似値を削減することになる場合があり、類似値を減らすように作用しなければならない場合において、逆に類似値を増やすことになる場合がある。
【００８０】
つまり、従来技術の類似値では、テンプレート画像と入力画像とのマッチングを、正しく評価できない。その結果、従来技術では、誤認識が発生しやすいし、認識ができても認識結果が先鋭でないという問題点がある。
【００８１】
そこで本発明では、Θの倍角の余弦（ｃｏｓ（２Θ））を、類似値の式に使用している。こうすると、ｃｏｓΘ’＝−ｃｏｓΘとなっていても、（数６）の倍角公式から、ｃｏｓ（２Θ’）＝ｃｏｓ（２Θ）となる。つまり、類似値を増やすように作用しなければならない場合には、背景部分に影響されずに、類似値は高くなる。したがって、背景部分の輝度のばらつきがあっても、正当にマッチングを評価することができる。以上の点は、２倍角だけでなく、４倍角や６倍角などでも同様に成立する。
【００８２】
しかも、偶数倍角の余弦を用いても、類似値が角Θに対して単調減少になる性質は変わらない。
【００８３】
従って、本発明によれば、偶数倍角評価により、背景の輝度条件に関わらず、安定してパターンを抽出できる。
【００８４】
より具体的には、本発明では、次式により、類似値を定義する。
【数８】

【００８５】
ここで、（数８）では、加算、積算のみからなるので、類似値は、入力画像及びテンプレート画像のそれぞれの評価ベクトルについて、線形である。したがって、（数８）をフーリエ変換すると、フーリエ変換の離散相関定理（参考文献：高速フーリエ変換 宮川洋訳 科学技術出版社）により、
【数９】

となる。
【００８６】
また、（数９）を、逆フーリエ変換すれば、（数８）の類似値が得られる。そして、（数９）を注意深く観察すると、次の２点が明らかとなろう。
【００８７】
まず、（１）直交変換した後の変換値においては、テンプレート画像に係るフーリエ変換値と、入力画像に係るフーリエ変換値は、単純に積和すればよい。
【００８８】
また、（２）テンプレート画像に係るフーリエ変換値と、入力画像に係るフーリエ変換値は、同時に求める必要はなく、テンプレート画像に係るフーリエ変換値を、入力画像のそれに先行して、求めておいてもかまわない。
【００８９】
そこで、本形態では、テンプレート画像処理部１００に、記録手段５を設け、入力画像の入力に先立ち、圧縮手段４の出力を記録しておくことにしている。これにより、入力画像が入力画像処理部２００に入力された後は、テンプレート画像処理部１００は、何らテンプレート画像の処理を行う必要がなく、入力画像処理部２００及び積算手段１０より後段の処理に、処理能力を集中させることができ、一層処理を高速化できる。
【００９０】
次に、評価ベクトル生成手段２よりも後段の構成を説明する。図１に示すように、テンプレート画像処理部１００において、評価ベクトル生成手段２から出力される、テンプレート画像の評価ベクトルは、直交変換手段３により、フーリエ変換され、圧縮手段４に出力される。
【００９１】
圧縮手段４は、フーリエ変換後の評価ベクトルを、削減して、記録手段５に格納する。図３に示すように、変換後の評価ベクトルは、ｘ，ｙ両方向について、高低さまざまな周波数成分を含んでいる。本発明者らの実験によれば、すべての周波数成分について処理を行わなくとも、低周波数成分（例えば、ｘｙ両方向について、低周波側半分ずつ等）について、処理を行えば、十分な精度が得られることがわかっている。なお、図３において、斜線を付していない領域（−ａ≦ｘ≦ａ，−ｂ≦ｙ≦ｂ）が、元の領域であり、斜線を付した領域（−ａ／２≦ｘ≦ａ／２，−ｂ／２≦ｙ≦ｂ／２）が削減後の領域である。即ち、処理量は、１／４倍となる。
【００９２】
このようにすれば、処理対象を削減して、さらに高速な処理を実現できる。
【００９３】
なお、圧縮手段４及び記録手段５は、データ量が小さいときや、高速性が要求されないときは、省略することも可能である。
【００９４】
次に、入力画像処理部２００について説明する。入力画像処理部２００は、テンプレート画像処理部１００とほぼ同じような処理を行う。即ち、エッジ抽出手段６は、（数２）、（数３）により、次式で定義される入力画像のエッジ法線方向ベクトルを出力する。
【数１０】

【００９５】
また、評価ベクトル生成手段７は、エッジ抽出手段６から入力画像のエッジ法線方向ベクトルを入力し、次の２つの式で定義される、入力画像の評価ベクトルを出力する。
【数１１】

【数１２】

【００９６】
テンプレート画像処理部１００と異なる点は、ｎで割る正規化処理を行っていない点だけである。即ち、偶数倍角による評価、長さ１への正規化処理、ノイズ除去処理は、テンプレート画像処理部１００と同様に、行う。
【００９７】
次に、評価ベクトル生成手段７よりも後段の構成を説明する。図１に示すように、入力画像処理部２００において、評価ベクトル生成手段７から出力される、入力画像の評価ベクトルは、直交変換手段８により、フーリエ変換され、圧縮手段９に出力される。
【００９８】
圧縮手段９は、フーリエ変換後の評価ベクトルを、削減して、積算手段１０へ出力する。ここで、圧縮手段９は、圧縮手段４と同一の周波数帯（本例では、ｘｙ両方向について、低周波側半分ずつ等）に、処理対象を削減する。
【００９９】
もちろん、圧縮手段９は、データ量が小さいときや、高速性が要求されないときは、省略することも可能であるが、テンプレート画像処理部１００において、圧縮手段４を省略するときは、圧縮手段９も同様に省略する。
【０１００】
次に、積算手段１０以降を説明する。さて、テンプレート画像処理部１００及び入力画像処理部２００の処理が完了すると、積算手段１０は、記録手段５と圧縮手段９から、テンプレート画像及び入力画像の、それぞれの評価ベクトルのフーリエ変換値を入力する。
【０１０１】
そこで、記録手段５は、（数９）による積和演算を行い、結果（類似値Ｌのフーリエ変換値）を逆直交変換手段１１へ出力する。
【０１０２】
逆直交変換手段１１は、類似値Ｌのフーリエ変換値を、逆フーリエ変換し、類似値ＬのマップＬ（ｘ，ｙ）をマップ処理手段１２へ出力する。マップ処理手段１２は、このマップＬ（ｘ，ｙ）から、値の高い点（ピーク）を抽出し、その位置と値を出力する。なお、マップ処理手段１２以降は、必要に応じて自由に構成することができる。
【０１０３】
次に、図４，図５を用いて、図２のテンプレート画像による処理例を説明する。まず、入力画像が図４（ａ）のようであったとすると、エッジ抽出手段６は、図４（ｂ）のようなｘ方向エッジ成分と、図４（ｃ）のようなｙ方向エッジ成分を抽出する。
【０１０４】
そして、上述した処理の結果、図５（ａ）のような類似値のマップＬ（ｘ，ｙ）が得られる。ここで、「最大値」と記載した矢印の先端部分が、このマップのピークであり、図５（ｂ）の入力画像と比較すれば明らかなように、正しい点を１点のみで先鋭に認識できていることが分かる。
【０１０５】
従来技術では、入力画像のサイズをＡ（＝２^γ）とし、テンプレート画像のサイズをＢとすると、テンプレート画像を順次入力画像上を走査し、各位置での（数１）のρ_voを求めるために、
積の回数：２ＡＢ
の計算回数が必要になる。ここで、計算回数は、計算コストが高い、積の回数で評価する。
【０１０６】
一方、本形態では、直交変換手段３、８による２回のＦＦＴ、積算手段１０の積和計算、及び、逆直交変換手段１１による１回の逆ＦＦＴが必要で、
積の回数：３｛（２γ−４）Ａ＋４｝＋２Ａ
の計算回数ですむ。
【０１０７】
これらの計算回数を比較すると、例えば、Ａ＝２５６×２５６＝２¹⁶とし、Ｂ＝６０×６０とした場合、本形態による積の計算回数は、従来技術の積の計算回数の約１／１００になり、非常に高速な処理が可能となる。
【０１０８】
従来技術に係る（数１）のように、非線形な式では、フーリエ変換の離散相関定理を利用することができない。
【０１０９】
したがって、従来技術では、本形態の図１に示しているように、テンプレート画像に係る処理を、入力画像のそれに先行させることができない。つまり、従来技術では、テンプレート画像と入力画像の両方の処理を同時に行わざるを得ず、この点においても、本形態は、従来技術よりも処理が高速となる。
【０１１０】
（実施の形態２）
本形態では、図９に示すように、図１の要素に、共役圧縮手段１３及び共役復元手段１４を追加している。共役圧縮手段１３は、記録手段５から読み出される、テンプレート画像の、評価ベクトルのフーリエ変換値を、フーリエ変換の複素共役性を利用して、さらに半分に削減する。
【０１１１】
この点について、説明する。実数のフーリエ変換により得られるスペクトルについては、次式が成り立つ。
【数１３】

【０１１２】
即ち、ある点のスペクトル値は、ｕｖ座標系において、その点に対称な位置のスペクトル値の複素共役に等しい。この性質を利用すれば、図１０（ａ）に示すように、処理量を、ｕｖ座標系の半分に削減することができる。
【０１１３】
さらに、ここでは、圧縮手段４及び圧縮手段９による圧縮が行われるため、図１０（ｂ）に示すように、元のデータに対して、データ量を１／８に削減することができ、処理の一層の高速化と記憶容量の節約を図ることができる。
【０１１４】
このため、本形態では、共役圧縮手段１３が記録手段５から読み出したフーリエ変換値を半分にして、積算手段１０へ出力し、共役復元手段１４が積算手段１０からの出力を、倍にする処理を行って、逆直交変換手段１１へ出力することとしている。
【０１１５】
以上が、本発明の基本である。以下さらに、以上の基本を応用して、有用な技術を開示する。
【０１１６】
（実施の形態３）
本形態では、第１例と、第２例の２つの例を説明する。いずれの例においても、拡縮されたテンプレート画像を用いて、顔抽出を行う技術を開示する。なお、この技術によれば、サイズが異なる、相似形のテンプレートによる処理を、効率良く、かつ、高速に、実施できる。なお、この技術の用途は、顔抽出だけに限られるものではない。
【０１１７】
「第１例」
さて、図１１（ａ）は、本発明の実施の形態３（第１例）によるテンプレート画像処理部のブロック図である。図示しているように、図１又は図９のテンプレート画像処理部１００に若干の変更を加えてある。
【０１１８】
まず、テンプレート画像が、テンプレート画像処理部１０１に入力されると、拡縮手段１５で拡縮されてから、エッジ抽出手段１へ拡縮後のテンプレート画像がエッジ抽出手段１へ出力される。
【０１１９】
また、評価ベクトル生成手段２がテンプレート画像の評価ベクトルを出力すると、これを加算手段１６が加算処理し、加算処理後のデータが、直交変換手段３へ出力される。
【０１２０】
ここで、加算手段１６は、複数のサイズ範囲毎に、次式により加算処理を行う。
【数１４】

【０１２１】
実施の形態１，２では、１つのテンプレート画像だけで処理を行っていたが、本形態では、複数テンプレートで、複数大きさのテンプレートを処理する。
【０１２２】
その上、ある範囲毎に、テンプレートの大きさを分け、テンプレート処理結果を重ねる。
【０１２３】
図１２は、入力画像の例を示し、この入力画像には、人物の顔の像が含まれている。このような場合、図１３（ａ）のような、真っ直ぐな顔のテンプレート画像を用意する。また、図１３（ｄ）のような、ある角度傾けた顔のテンプレート画像を用意する。傾けた顔について、どのような角度の範囲で用意するかは、適宜選択できる。
【０１２４】
そして、本形態では、図１３（ａ）、（ｄ）のテンプレート画像のエッジを抽出すると、それぞれ図１３（ｂ）、（ｅ）のようになる。また、図１３（ａ）又は（ｄ）の画像を入力したとき、図１３（ｃ）、（ｆ）のような、エッジ画像に対応する評価ベクトルを生成するものである。
【０１２５】
そして、実施の形態１，２と同様に、テンプレート画像に関するデータを、記録手段５に記録しておく。
【０１２６】
入力画像が入力されたならば、実施の形態１，２と同様、入力画像処理部２００及び積算手段１０以降の処理を行う。
【０１２７】
なお、各テンプレート画像のすべてのサイズ範囲の重ねテンプレートデータに対して、それぞれの類似値マップを求める。
【０１２８】
これにより、例えば、図１４に示すような、検出結果が得られる。図１４から、明らかなように、顔の全体的な位置だけでなく、その大きさや、目鼻口等の主要な顔部品の位置も把握することができる。
【０１２９】
このとき、マップ処理手段１２の後段において、図１４のように、その処理結果を、ディスプレイに表示したり、プリンタに出力したりするとよい。
【０１３０】
「第１例」によると、次の効果がある。
（１）あるサイズ範囲のテンプレートデータを加算し、重ねたテンプレートで処理を行っても、テンプレートに類似する部分は、高い類似値を示すことが多いため、すべてのサイズの類似値を求めた場合、テンプレートの個数をＮとし、すべてのサイズの個数をＭとした場合、積和計算部と逆直交変換部の処理をＮＭ回繰り返さなければいけないのに対し、重ねテンプレートでは、重ね範囲の幅をＨとすると処理回数は、ＮＭ／Ｈとなり、効率化できる。従って、顔抽出などを行う際に、より高速な処理が可能となる。
【０１３１】
（２）また、単に、顔の位置データを出力するだけでなく、図１４のように、顔候補位置にテンプレートを重ねることで、目、鼻、口の大まかな候補領域を抽出できる。
【０１３２】
後の形態で述べるように、これらの領域の画像を更に細かく処理することで、より正確な目鼻口等の顔部品の位置を抽出することも可能である。
【０１３３】
「第２例」
本例では、図１１（ｂ）に示すように、テンプレート画像処理部１０２を構成する。即ち、「第１例」に対して、加算手段１６の位置を、圧縮手段４と記録手段５との間へ移動している。こうすると、加算手段１６は、次式で加算処理をするように構成する。
【数１５】

【０１３４】
フーリエ変換前に線形であった、類似値式は、変換後も線形である。したがって、「第１例」から「第２例」のように、加算手段１６の位置を変更できる。
【０１３５】
しかもこうすると、加算手段１６の処理対象が、圧縮手段４によって圧縮されているため、「第１例」よりも、加算手段１６の処理対象を削減でき、一層高速化を図ることができるため、好適である。
【０１３６】
（実施の形態４）
図１５および図１６を用いて、実施の形態４を説明する。本形態では、実施の形態１，２において説明した類似値マップの極大点の値を、より強いピーク値にする効率的な技術を開示する。
【０１３７】
さて一般に、類似値マップでは、テンプレート画像と合致する部分に、ピークが、現れる。本形態では、極大点を含めた周辺のピークパターンｐを類似値マップにフィルターとしてかけ、類似値マップ中のピークパターンに類似する部分の値を増幅する。
【０１３８】
即ち、図１５に示すように、本形態では、実施の形態１を示す、図１の構成に加え、ピークパターン処理部３００を追加している。
【０１３９】
また、図１６は、このピークパターンのためのマスクを例示する。図１６に指名しているように、このピークパターンでは、平均値が、０になるように、正規化してある。
【０１４０】
そして、ピークパターン処理部３００において、直交変換手段１７は、このピークパターンをフーリエ変換し、圧縮手段１８はフーリエ変換値を圧縮し、圧縮後のデータが記録手段１９に記録される。
【０１４１】
また、マスクを使用しているため、類似値式は、（数８）そのものではなく、マスクを反映した次式による。
【数１６】

【０１４２】
したがって、積算手段１０は、記録手段５、記録手段１９及び圧縮手段９からデータを読み込んで、積和演算を行い、ピークパターンで修正された、類似値のフーリエ変換値を出力する。
【０１４３】
なお、類似値マップＬにピークパターンフィルタを掛けるには、次式によればよいが、このようにすると、大量の積和計算が必要となり、効率的でない。
【数１７】

【０１４４】
一方、本形態のように、（数１７）のような、大量の計算を行わず、（数１６）により、簡便且つ正確に処理を行う。
【０１４５】
したがって、本形態によれば、効率的に、類似値マップのピーク点を増幅できる。また、ピークパターンを反映して、より先鋭且つ安定に、テンプレートと類似する部分を入力画像中から検出することができる。
【０１４６】
（実施の形態５）
本形態では、入力画像のエッジ部分とテンプレート画像のエッジ部分での類似値に加え、テンプレート画像の領域内の画素平均も、類似判定に加えている。
【０１４７】
構成は、図１７に例示している。即ち、図１５に似ており、実施の形態１に係る図１に対し、マスクパターン処理部４００を追加している。
【０１４８】
但し、図１５と異なり、マスクパターン処理部４００は、ピークパターンを入力するのではなく、テンプレート画像を入力し、この画像に依存するマスクパターンを生成するマスクパターン生成手段２０を設けている点が異なる。
【０１４９】
また、図１５と同様、マスクパターン生成手段２０の出力は、直交変換手段２１でフーリエ変換され、圧縮手段２２で圧縮されてから、記録手段２３に記録される。
【０１５０】
また、マスクを使用しているため、類似値式は、（数８）そのものではなく、マスクを反映した次式による。
【数１８】

【０１５１】
このまま積算すると、大量の積和計算が必要となり、効率的でないのは、実施の形態４で述べたのと、同様である。
【０１５２】
これをフーリエ変換すると、
【数１９】

となり、非常に簡潔に演算を行える。
【０１５３】
したがって、積算手段１０は、（数１９）による積和演算を行う。
【０１５４】
さて次に、テンプレート画像とマスクパターンとの関係について、図１８を用いて説明する。ここでは、テンプレート画像の領域内の画素平均を、類似判定に加えるため、マスクパターン生成手段２０は、図１８（ａ）のようなテンプレート画像に対し、図１８（ｂ）のようなマスクパターンを生成する。
【０１５５】
より詳しくは、図１８（ａ）のテンプレート画像の場合、画素平均を求めたい内部（丸の内側）の各点に１／Ｎの値を設定し、その他を０とする。ここで、Ｎは、内部の点数であり、マスクパターンの全ての点の値を加えると、１になる。
【０１５６】
本形態によれば、画像内部の画素値の平均値も類似値に加えることができ、より正確に、かつ効率よく、入力画像中から対象物を抽出することが可能になる。
【０１５７】
なお、各画素の２乗値の平均を求めるためには、入力画像処理部で、入力画像の各画素の２乗をしたデータを作成し、同様の処理を行うことで可能となる。
【０１５８】
従って、領域内の平均だけでなく、分散値なども、効率よく求めることができる。
【０１５９】
（実施の形態６）
本形態は、テンプレート画像に対し、左右対称な画像について、要領良く処理できる技術を開示する。
【０１６０】
その構成は、図１９に例示している。即ち、実施の形態１に係る図１の構成に加え、記録手段５と積算手段１０との間に、対称ベクトル生成手段２４を追加している。また、本形態では、実施の形態１と同様の、類似値式に（数８）を用いる。
【０１６１】
さて、次に左右対称となるテンプレート画像の取り扱いについて、説明する。例えば、図２０（ａ）のテンプレート画像が元の画像だとすると、それを左右反転させたテンプレート画像は、図２０（ｂ）のようになる。
【０１６２】
また、これらのテンプレート画像についての、エッジ法線方向ベクトルの関係は、次式に示すとおりである。
【数２０】

【０１６３】
そして、左右反転させたテンプレート画像の評価ベクトルは、次式のようになる。
【数２１】

【０１６４】
また、フーリエ変換については、
【数２２】

なる関係があるから、
（数２１）に（数２２）を適用して、
【数２３】

【０１６５】
つまり、左右反転させたテンプレート画像の評価ベクトルは、元のテンプレート画像の評価ベクトルから、正負反転するなど、簡単な処理で生成できる。
【０１６６】
したがって、図１９において、対称ベクトル生成手段２４が記録手段５の元のテンプレート画像の評価ベクトルについて、（数２３）を適用するだけで、左右反転させたテンプレート画像の評価ベクトルを得ることができる。
【０１６７】
この点、単純に考えて、図２０（ａ）の画像から、図２０（ｂ）の画像を生成し、図２０（ｂ）の画像から、評価ベクトルを再計算するなど、煩雑な処理を行う必要はないのである。
【０１６８】
これにより、左右反転したテンプレート画像の評価ベクトルを、直接的な演算なしに、生成することができ、高速化を図ることができる。また、テンプレート画像について、その左右反転させたテンプレート画像をわざわざ記録しておく必要が無くなるので、記録容量の節約を図ることができる。
【０１６９】
（実施の形態７）
本形態の形態では、実施の形態３で述べた顔抽出の処理に、目眉の抽出処理を追加する。
【０１７０】
実施の形態３で述べた処理により、図２２（ａ）の入力画像から、図２２（ｂ）のように、大まかに目眉候補領域を抽出できる。
【０１７１】
この目眉候補領域の画像を対象に、図２３（ａ）〜図２３（ｄ）に示す、点双列相関係数フィルタを各点に掛け、点双列相関値のマップを作成し、マップ中の相関値の高くなる点を、それぞれ、目の中心位置３００２、眉の中心位置３００３とする。
【０１７２】
点双列相関係数ηは、次式で定義される（参考文献：多変量解析ハンドブック 現代数学社 １７頁）。
【０１７３】
図２３において、図２３（ａ）は全領域の位置関係を示し、図２３（ｂ）は、全領域マスク、図２３（ｃ）は領域１マスク、図２３（ｄ）は領域２マスクを示す。
【０１７４】
この点双列相関係数によるフィルタ形状を、図２３（ａ）に示すようにすると、図２２（ｃ）に示すように、眉中心３００２や目中心３００３を抽出できることが期待される。
【０１７５】
次に、点双列相関係数によるフィルタ処理について、説明する。
【０１７６】
まず、（数２４）のうち、主要な成分は、次式のとおりである。
【数２５】

【０１７７】
それぞれの成分をフーリエ変換すると、
【数２６】

【数２７】

【数２８】

【０１７８】
以上の処理を実現するには、例えば、図２１のように構成する。まず、全領域、領域１、領域２の各マスクは、直交変換手段５１〜５３でフーリエ変換される。
【０１７９】
また、目眉候補領域抽出手段５４には、入力画像と、実施の形態３で述べた顔抽出の結果がマップ処理手段１２を介して入力される。目眉候補領域抽出手段５４は、これらの入力に基づいて、図２２（ａ）から図２２（ｂ）のような、目眉候補領域のみを取り出す。この目眉候補領域のデータは、直交変換手段５５によりそのままフーリエ変換されると共に、２乗手段５６によって２乗されたのち、直交変換手段５７によりフーリエ変換される。
【０１８０】
そして、ηマップ作成手段６５の前段にある、逆直交変換手段６２には、乗算手段５８を介して（数２７）に示すデータが入力される。又同様に、逆直交変換手段６３には、乗算手段６０を介して（数２６）に示すデータが入力され逆直交変換手段６４には、乗算手段６１を介して（数２８）に示すデータが入力される。そして、これらの逆直交変換手段６２〜６４は、入力したデータを、逆フーリエ変換し、ηマップ作成手段６５に出力する。
【０１８１】
次に、ηマップ作成手段６５は、逆直交変換手段６２〜６４からの入力に対して、（数２４）の演算を行い、点双列相関係数のマップη（ｘ，ｙ）を出力する。
【０１８２】
そして、目眉中心抽出手段６６は、ηマップ作成手段６５が出力するマップη（ｘ，ｙ）の中から、値の高い２点を抽出し、目眉それぞれの中心として出力する。
【０１８３】
以上の構成により、図２３（ａ）のフィルタのサイズが１５×１５（画素）とし、入力画像の画素数をＮとすると、大まかに、１５×１５×Ｎ＝２２５Ｎ回の掛け算処理が必要になる。
【０１８４】
一方、本形態によると、積計算回数は、大まかに、Ｎ＋｛（２γ−４）Ｎ＋４｝×５＋１２Ｎ＋Ｎ＋Ｎ＝５Ｎ（２γ−１）＋２０回である。但し、簡単にするため、√計算を掛け算１回の処理と同等とみなす。
【０１８５】
つまり、γが２２以下で、（本形態の計算量）＜（従来技術の計算量）となる。
【０１８６】
ここで、Ｎ＝２²²以上の巨大な画像は、通常用いないので、本形態の処理の方が、従来の処理よりも演算回数が少なく、高速であることが分かる。
【０１８７】
このように、本形態によれば、より効率よく、点双列相関係数フィルタ処理を行える。
【０１８８】
なお、図２３に示したマスクを種々変更することにより、目尻、口端、鼻孔、黒目領域など、他の顔部品の位置を求めることができる。
【０１８９】
（実施の形態８）
本形態では、実施の形態３による、顔抽出の機能を拡張し、顔画像から顔器官である口範囲を取り出す技術を開示する。
【０１９０】
さて、実施の形態３において、図１４を用いて説明したように、図２５（ａ）の入力画像から、図２５（ｂ）に示すような、口候補領域を取り出すことができる。
【０１９１】
取り出した、口候補領域を、Ｙ軸に投影（Ｘ軸に沿って、画素値の総和ｈをプロット）すると、おおむね図２５（ｃ）に示すような、グラフになる。
【０１９２】
ここで、この総和ｈを定義すると、次式のようになる。
【数２９】

【０１９３】
この総和を、直交変換を用いて効率良く求めるには、図２５（ｄ）のような、マスクを用意する。そして、このマスクを含めて、（数２９）を書き換えると、次式のとおりである。
【数３０】

【０１９４】
これをフーリエ変換すると、
【数３１】

【０１９５】
即ち、この投影値を求めるには、入力画像とマスクパターンをフーリエ変換して、（数３１）の計算を行い、その結果を逆フーリエ変換すればよい。そうすると、投影値ｈのマップｈ（ｘ，ｙ）を求めることができる。
【０１９６】
以上の結果から、例えば、図２４のように構成すると良い。図２４に示すように、入力画像とマスクパターンを、それぞれ直交変換手段２５、２６でフーリエ変換する。
【０１９７】
そして、変換値を、積算手段２７で乗算し、逆直交変換手段２８で逆フーリエ変換する。すると、投影値のマップｈ（ｘ，ｙ）が得られる。
【０１９８】
また、得られたマップｈ（ｘ，ｙ）は、投影データ取出手段２９に展開され、極大点抽出手段３０は、このマップから２つの極大点を求めて、これらの極大点間を口範囲として出力する。
【０１９９】
勿論、実施の形態１等と同様に、マスクパターンを前もって、フーリエ変換しておいても良く、入力画像の処理に処理能力を集中でき、処理は高速である。
【０２００】
本形態によれば、顔器官である口の位置（より詳しくは、上下両方の唇の位置）を正確に求めることができる。
【０２０１】
（実施の形態９）
本形態では、その第１、第２例により、実施の形態３による、顔抽出の機能を拡張し、顔画像に所定の処理を施す技術を開示する。
【０２０２】
「第１例」
本例は、抽出した顔画像に電子透かしを埋め込む。即ち、図２６のように構成する。図１１のマップ処理手段１２の後段に、結果利用部６００を設ける。
【０２０３】
結果利用部６００は、入力画像そのものと、マップ処理手段１２による顔位置を参照して、入力画像を、顔画像と、顔画像でない部分とに、分離する、顔画像切出手段３１を有する。顔画像部分は、顔画像切出手段３１から電子透かし埋込手段３２へ出力され、所定の電子透かしが顔画像部分に埋め込まれる。ここで、この電子透かし自体は、周知のものであって良い。
【０２０４】
一方、顔画像を除いた部分は、顔画像切出手段３１から画像合成手段３３へそのまま出力される。
【０２０５】
画像合成手段３３は、電子透かしを埋め込んだ顔画像と、顔画像を除いた部分とを合成し、電子透かし埋め込み済みの、画像を出力する。
【０２０６】
したがって、例えば、図２７（ａ）の入力画像について、電子透かしを顔部分に埋め込んだ、図２７（ｂ）の出力画像が得られる。
【０２０７】
なお、電子透かしとしては、被写体のモデルの氏名や撮影日などが、好適である。
【０２０８】
本例により、改ざんの対象になりやすい顔部分に重点的に透かしデータを簡単に埋め込むことができる。
【０２０９】
「第２例」
本例は、抽出した顔画像のみに、特定の編集を加える。即ち、図２８のように構成する。図１１のマップ処理手段１２の後段に、結果利用部６０１を設ける。
【０２１０】
結果利用部６０１は、入力画像そのものと、マップ処理手段１２による顔位置を参照して、入力画像を、顔画像と、顔画像でない部分とに、分離する、顔画像切出手段３１を有する。顔画像部分は、顔画像切出手段３１から画像補正手段３４へ出力され、所定の編集が顔画像部分にのみ適用される。ここで、この編集自体は、周知のものであって良い。
【０２１１】
一方、顔画像を除いた部分は、顔画像切出手段３１から画像合成手段３３へそのまま出力される。
【０２１２】
画像合成手段３３は、編集済みの顔画像と、顔画像を除いた部分とを合成し、変種済みの、画像を出力する。
【０２１３】
したがって、例えば、図２９（ａ）の入力画像について、図２９（ｂ）のように、顔領域を抽出し、さらに、図２９（ｃ）のように、顔部分だけを編集した、出力画像が得られる。
【０２１４】
なお、画像補正手段３４による編集としては、逆光で撮影された人物像があるとき、顔が暗くなりすぎて見づらい様な場合、顔部分の画像に対して、顔色を白っぽくするなどが考えられるが、任意に選択して差し支えない。
【０２１５】
本例により、簡単に、顔部分のみを選択的に、つまり、顔以外の部分に影響を与えずに、画像補正ができる。
【０２１６】
「第３例」
本例では、抽出した顔画像を、より見やすくすべく補正する。即ち、図３０のように構成する。まず、図１１のマップ処理手段１２の後段に、結果利用部６０２を設ける。
【０２１７】
結果利用部６０２は、入力画像そのものを、顔画像と、顔画像でない部分とに、分離する、顔画像切出手段３１を有する。切り出された顔画像は、顔内部領域抽出手段３５と画像補正手段３４に出力される。
【０２１８】
顔内部領域抽出手段３５は、実施の形態３で述べた要領により、テンプレートを用いて、切り出された顔画像から、目鼻口等の主要な顔部品の位置を求め、顔内部画像（全て顔の内部に位置する画像）を抽出する。本例では、顔内部画像は、顔中心（例えば、鼻の中央）を中心とする、一定サイズの矩形領域内の画像としたが、この矩形領域のサイズを変更してもよいし、顔中心からやや外れた位置を中心とするようにしても差し支えない。
【０２１９】
抽出された顔内部画像は、画像特徴抽出手段３６へ出力される。画像特徴抽出手段３６は、顔画像の補正に役立つ特徴を求める。本例では、画像特徴抽出手段３６は、顔内部画像の輝度分布を求め、輝度ヒストグラムとして、補正関数決定手段３７へ出力する。勿論、画像特徴抽出手段３６が補正関数決定手段３７へ輝度ヒストグラムにおける、輝度最小値と輝度最大値とを出力するようにしても良い。
【０２２０】
補正関数決定手段３７は、画像特徴抽出手段３６から入力した特徴を参照して、顔内部画像の明暗がより明確になるような補正関数を決定し、この補正関数を画像補正手段３４へ出力する。
【０２２１】
本例では、上述のように、画像特徴抽出手段３６は、輝度ヒストグラムを出力するので、補正関数決定手段３７は、このヒストグラムにあらわれる最小輝度が、ダイナミックレンジにおける最小の輝度となり、かつ、最大輝度が、ダイナミックレンジにおける、最大の輝度となるように、補正関数を求め、画像補正手段３４へ出力する。
【０２２２】
因みに、図３０の輝度ヒストグラムの横軸は輝度値であり、輝度表現できるダイナミックレンジの一部しか活用されていないことが理解されよう。そこで、ダイナミックレンジをフル活用するように、補正関数決定手段３７が補正関数を決定することによって、画像補正手段３４は、より見やすく、具体的には、見た目の階調豊かに、顔画像を、補正できるのである。この補正の結果、顔画像以外の部分は、白くとんだり、暗くなったりすることが考えられるが、顔画像は確実に見やすくなる。
【０２２３】
ここで、補正関数は、顔内部画像のみに基づいて、決定されるため、顔でない部分によって、補正がゆがめられることはなく、補正後において、顔の見やすさが向上することを、保証できる。
【０２２４】
なお、図３０では、切り出した顔画像のみに補正を施すように表示しているが、入力画像全部に、同様の補正を施すこともできる。
【０２２５】
さらに、画像特徴抽出手段３６が求める特徴としては、上述したように、輝度など画像の明るさ／暗さを代表する指標を用いることができるし、その他、彩度平均や色相平均を用いても良い。
【０２２６】
彩度平均を用いるときは、補正関数決定手段３７は、画像補正手段３４に、例えば彩度増幅係数を出力するようにし、色相平均を用いるときは、補正関数決定手段３７は、画像補正手段３４に、例えば色相回転角を出力するようにするとよい。
【０２２７】
ここで、彩度平均を特徴として用いると、色が薄い顔画像を、より色鮮やかに補正できる。また、色相平均を用いると、赤みが少ない顔画像を、より赤みがからせることができる。いずれにしても、顔画像を、より自然な感じにすることができる。加えて、輝度、彩度平均又は色相平均から、二種以上組み合わせたものを、特徴としても良い。
【０２２８】
【発明の効果】
本発明によれば、撮影条件の違いや背景の違いの影響が少なく安定に、かつ、高速に、パターンを抽出できる。
【０２２９】
また、複数のテンプレートのパターンの抽出も高速にできる。
【０２３０】
類似値マップのピーク点を増幅でき、より安定に、テンプレートと類似する部分を入力画像中から検出できる。
【０２３１】
画像内部の属性を類似値に加味し、より正確・効率よく、処理できる。
【０２３２】
反転したテンプレート画像のデータを、間接的に生成でき、テンプレート画像のデータ量を減らせる。
【０２３３】
処理対象を合理的に削減できる。
【０２３４】
点双列相関係数フィルタ処理で、目眉の中心位置を抽出できる。
【０２３５】
顔器官である口の位置を、正確に求めることができる。
【０２３６】
改ざんの対象になりやすい顔部分に、重点的に透かしデータを埋め込むことができる。
【０２３７】
顔部分のみに選択的に、画像補正ができるようになる。
【０２３８】
顔内部画像の特徴にのみ基づいて、顔画像の見やすさを確実に向上できる。
【図面の簡単な説明】
【図１】本発明の実施の形態１における画像処理装置のブロック図
【図２】（ａ）同テンプレート画像の例示図
（ｂ）同テンプレート画像のエッジ抽出画像（ｘ成分）の例示図
（ｃ）同テンプレート画像のエッジ抽出画像（ｙ成分）の例示図
【図３】同評価ベクトルの圧縮処理説明図
【図４】（ａ）同入力画像の例示図
（ｂ）同入力画像のエッジ抽出画像（ｘ成分）の例示図
（ｃ）同入力画像のエッジ抽出画像（ｙ成分）の例示図
【図５】（ａ）同類似値マップの例示図
（ｂ）同入力画像の例示図
【図６】同内積の正負反転説明図
【図７】（ａ）同テンプレート画像の例示図
（ｂ）同入力画像の例示図
【図８】（ａ）同テンプレート画像の例示図
（ｂ）同テンプレート画像の例示図
【図９】本発明の実施の形態２における画像処理装置のブロック図
【図１０】（ａ）同共役性を示すグラフ
（ｂ）同共役性を示すグラフ
【図１１】（ａ）本発明の実施の形態３（第１例）におけるテンプレート画像処理部のブロック図
（ｂ）同（第２例）によるテンプレート画像処理部のブロック図
【図１２】同入力画像の例示図
【図１３】（ａ）同テンプレート画像の例示図
（ｂ）同エッジ抽出画像の例示図
（ｃ）同拡縮テンプレート画像の例示図
（ｄ）同テンプレート画像の例示図
（ｅ）同エッジ抽出画像の例示図
（ｆ）同拡縮テンプレート画像の例示図
【図１４】同顔抽出結果の例示図
【図１５】本発明の実施の形態４における画像処理装置のブロック図
【図１６】同ピークパターンの例示図
【図１７】本発明の実施の形態５における画像処理装置のブロック図
【図１８】（ａ）同テンプレート画像の例示図
（ｂ）同マスクパターンの例示図
【図１９】本発明の実施の形態６における画像処理装置のブロック図
【図２０】（ａ）同元のテンプレート画像の例示図
（ｂ）同左右反転後のテンプレート画像の例示図
【図２１】本発明の実施の形態６における画像処理装置の一部ブロック図
【図２２】（ａ）同入力画像の例示図
（ｂ）同目眉候補領域の例示図
（ｃ）同認識結果の例示図
【図２３】（ａ）本発明の実施の形態７におけるフィルタ形状の例示図
（ｂ）同全領域マスクの説明図
（ｃ）同領域１マスクの説明図
（ｄ）同領域２マスクの説明図
【図２４】本発明の実施の形態８における画像処理装置の一部ブロック図
【図２５】（ａ）同入力画像の例示図
（ｂ）同口候補領域の説明図
（ｃ）同投影値のグラフ
（ｄ）同マスクパターンの説明図
（ｅ）同投影値マップ画像の例示図
（ｆ）同投影値のグラフ
【図２６】本発明の実施の形態９（第１例）における画像処理装置の一部ブロック図
【図２７】（ａ）同入力画像の例示図
（ｂ）同出力画像の例示図
【図２８】本発明の実施の形態９（第２例）における画像処理装置の一部ブロック図
【図２９】（ａ）同入力画像の例示図
（ｂ）同出力画像の例示図
【図３０】本発明の実施の形態９（第３例）における画像処理装置の一部ブロック図
【符号の説明】
１ エッジ抽出手段
２ 評価ベクトル生成手段
３ 直交変換手段
４ 圧縮手段
５ 記録手段
６ エッジ抽出手段
７ 評価ベクトル生成手段
８ 直交変換手段
９ 圧縮手段
１０ 積算手段
１１ 逆直交変換手段
１２ マップ処理手段
１３ 共役圧縮手段
１４ 共役復元手段
１５ 拡縮手段
１６ 加算手段
１７ 直交変換手段
１８ 圧縮手段
１９ 記録手段
２０ マスクパターン生成手段
２１ 直交変換手段
２２ 圧縮手段
２３ 記録手段
２４ 対称ベクトル生成手段
２５、２６ 直交変換手段
２７ 積算手段
２８ 逆直交変換手段
２９ 投影データ取出手段
３０ 極大点抽出手段
３１ 顔画像切出手段
３２ 電子透かし埋込手段
３３ 画像合成手段
６０１ 結果利用部
３４ 画像補正手段
５１ 直交変換手段
５２ 直交変換手段
５３ 直交変換手段
５４ 目眉候補領域抽出手段
５５ 直交変換手段
５６ ２乗手段
５７ 直交変換手段
５８ 乗算手段
５９ 減算手段
６０ 乗算手段
６１ 乗算手段
６２、６３、６４ 逆直交変換手段
６５ ηマップ作成手段
６６ 目眉中心抽出手段
１００ テンプレート画像処理部
１０１、１０２ テンプレート画像処理部
２００ 入力画像処理部
３００ ピークパターン処理部
４００ マスクパターン処理部
５００ 投影画像処理部
６００ 結果利用部

Claims (12)

1. テンプレート画像を入力しエッジ法線方向ベクトルを求めるステップと、
   絶対値が所定値を超える前記エッジ法線方向ベクトルの個数ｎを算出し、前記個数ｎにより正規化された評価ベクトルを求める生成ステップと、
   前記評価ベクトルを直交変換するステップと、
   前記テンプレート画像について求められた直交変換後の評価ベクトルと、入力画像について求められた直交変換後の評価ベクトルとについて、対応スペクトルデータを積和計算する積算ステップと、
   積和計算結果を逆直交変換して類似値のマップを生成する逆直交変換ステップとを含む画像のマッチング処理方法。
2. 積和計算を行う前に、直交変換の複素共役性を利用してデータ量を削減し、積和計算後にデータ量を復元するステップをさらに含むことを特徴とする請求項１記載の画像のマッチング処理方法。
3. テンプレート画像を各種サイズへ拡縮し、各サイズについて評価ベクトルの加算処理を行うことを特徴とする請求項１から２のいずれか記載の画像のマッチング処理方法。
4. テンプレート画像について、評価ベクトルの加算処理を、圧縮して処理量を削減する後に行うことを特徴とする請求項３記載の画像のマッチング処理方法。
5. テンプレート画像は、顔を模式化したものであることを特徴とする請求項１から４のいずれか記載の画像のマッチング処理方法。
6. 類似値のピークを急峻にするピークパターンを用意し、このピークパターンのデータを直交変換したものを、積和計算に反映させることを特徴とする請求項１から５のいずれか記載の画像のマッチング処理方法。
7. テンプレート画像に依存するマスクパターンを作成し、このマスクパターンのデータを直交変換したものを、積和計算に反映させると共に、前記マスクパターンは、テンプレート画像における画像内部の画素数平均を示すものであることを特徴とする請求項１から６のいずれか記載の画像のマッチング処理方法。
8. テンプレート画像について、元のテンプレート画像の評価ベクトルの正負符号を処理して、元のテンプレート画像に対し左右対称な画像の評価ベクトルを生成するステップを有し、生成した評価ベクトルを積和計算に反映させることを特徴とする請求項１から７のいずれか記載の画像のマッチング処理方法。
9. 抽出した顔画像に基づいて、点双列相関係数のマップを生成し、顔部品の位置を求めることを特徴とする請求項５記載の画像のマッチング処理方法。
10. 抽出した顔画像に基づいて、入力画像を、顔画像のみと、顔画像を除いた部分とに分離し、顔画像のみに電子透かしを埋込み、電子透かしを埋め込んだ顔画像と、顔画像を除いた部分とを、合成して出力することを特徴とする請求項５記載の画像のマッチング処理方法。
11. 抽出した顔画像に基づいて、入力画像を、顔画像のみと、顔画像を除いた部分とに分離し、顔画像のみに編集を施し、編集した顔画像と、顔画像を除いた部分とを、合成して出力することを特徴とする請求項５記載の画像のマッチング処理方法。
12. 入力画像を入力し、入力画像のエッジ法線方向ベクトルを求め、エッジ法線方向ベクトルから評価ベクトルを生成し、評価ベクトルを直交変換する入力画像処理部と、
    テンプレート画像について求められた直交変換後の評価ベクトルと、前記入力画像処理装置により入力画像について求められた直交変換後の評価ベクトルとについて、対応スペクトルデータを積和計算する積算手段と、
    積和計算結果を、逆直交変換して類似値のマップを生成する逆直交変換手段とを有する画像処理装置であって、
    前記テンプレート画像について求められた前記直交変換後の評価ベクトルは、
    テンプレート画像を入力しエッジ法線方向ベクトルを求め、
    絶対値が所定値を超える前記エッジ法線方向ベクトルの個数ｎを算出し、前記個数ｎにより正規化された評価ベクトルを求め、
    前記評価ベクトルを直交変換することにより生成されることを特徴とする画像処理装置。

Images