"入力画像から当該入力画像を特定するための局所特徴点、並びに当該局所特徴点の特徴量である局所特徴量の組を複数を抽出する画像処理装置であって、与えられた画像の局所特徴点と当該局所特徴点の局所特徴量とを抽出する局所特徴量抽出手段と、前記入力画像に対して、予め設定した角度での回転処理を行なう回転処理手段と、該回転処理手段による回転処理の前と後のそれぞれの画像に対して前記局所特徴量抽出手段を適用することで得られた局所特徴点及び特徴量の安定性を評価し、安定性の高い順に、予め設定された個数までの局所特徴点及び局所特徴量の組を、前記入力画像に対する有効な局所特徴点並びに局所特徴量として決定する決定手段とを有することを特徴とする画像処理装置。前記入力画像に対し、予め設定された縮小率による縮小処理を再帰的に実行することで複数の段階の縮小画像を生成する縮小画像生成手段とを備え、前記回転手段は、各縮小画像について回転処理を行ない、前記決定手段は、各縮小画像毎に、有効な局所特徴点並びに局所特徴量を決定することを特徴とする請求項1に記載の画像処理装置。前記決定手段で決定する有効な局所特徴点並びに局所特徴量の個数は、各縮小画像のサイズに依存することを特徴とする請求項1に記載の画像処理装置。コンピュータに読み込ませ実行させることで、前記コンピュータを、請求項1の画像処理装置が有する各手段として機能させるためのプログラム。請求項4に記載のプログラムを格納したコンピュータが読み込み可能な記憶媒体。入力画像から当該入力画像を特定するための局所特徴点、並びに当該局所特徴点の特徴量である局所特徴量の組を複数を抽出する画像処理装置の制御方法であって、局所特徴量抽出手段が、与えられた画像の局所特徴点と当該局所特徴点の局所特徴量とを抽出する局所特徴量抽出工程と、回転処理手段が、前記入力画像に対して、予め設定した角度での回転処理を行なう回転処理工程と、決定手段が、該回転処理工程による回転処理の前と後のそれぞれの画像に対して前記局所特徴量抽出工程を適用することで得られた局所特徴点及び特徴量の安定性を評価し、安定性の高い順に、予め設定された個数までの局所特徴点及び局所特徴量の組を、前記入力画像に対する有効な局所特徴点並びに局所特徴量として決定する決定工程とを有することを特徴とする画像処理装置の制御方法。携帯電話、スマートフォン、デジタルカメラ等、撮影機能を持つ携帯デバイスが増えているが、それぞれが持つ演算能力や搭載メモリは異なり、局所特徴量の算出処理やその特徴量記述サイズはこれらの機器のリソースに応じたものが好ましい。例えば、演算能力が高くなく、搭載メモリもそれほど多くないデジタルカメラでは局所特徴量抽出処理も簡素なもので且つ特徴量記述サイズが少ない事が好ましい。他方、PCやサーバーにおいては、携帯端末ほど局所特徴量の算出処理やその特徴量記述サイズの制約は厳しくないし、演算能力も高いと言える。この様に、各端末、PC/サーバーにおいて求められる局所特徴量の算出処理やその特徴量記述サイズが異なる事があり、この様な場合においても、支障をきたさないことが望まれる。換言すれば、局所特徴量の算出処理やそのトータルの局所特徴量記述サイズ異なったとしても、類似比較処理における精度が低下しない事が好ましいと言える。そのためには、サーバー上などリッチな環境ではトータルサイズの大きな局所特徴量を、プアな端末機器ではトータルサイズの小さな局所特徴量の記述を行い、これらの局所特徴量が包含関係部分一致関係にある事が望まれる。しかしながら、特許文献1では画像変動に対してロバストな安定した局所特徴点を絞り込む考えに留まっており、局所特徴点数を一回の処理で、スケーラブルに記述する事は出来ていない。本発明は係る問題点に鑑みなされたものであり、信頼性の高い局所特徴点及び局所特徴量を得る技術を提供しようとするものである。この課題を解決するため、例えば本発明の画像処理装置は以下の構成を備える。すなわち、入力画像から当該入力画像を特定するための局所特徴点、並びに当該局所特徴点の特徴量である局所特徴量の組を複数を抽出する画像処理装置であって、与えられた画像の局所特徴点と当該局所特徴点の局所特徴量とを抽出する局所特徴量抽出手段と、前記入力画像に対して、予め設定した角度での回転処理を行なう回転処理手段と、該回転処理手段による回転処理の前と後のそれぞれの画像に対して前記局所特徴量抽出手段を適用することで得られた局所特徴点及び特徴量の安定性を評価し、安定性の高い順に、予め設定された個数までの局所特徴点及び局所特徴量の組を、前記入力画像に対する有効な局所特徴点並びに局所特徴量として決定する決定手段とを有する。本発明の態様によれば、本発明では、1枚の画像に対するトータルの局所特徴量サイズの支配的なファクターが局所特徴点の個数であると考え、包含関係という課題の解決を行う。即ち、局所特徴点の再現性を評価し、再現性の高いものから記述し局所特徴量の個数を制御することで、画像に対するトータルの局所特徴量のサイズ調整を行う。というのも、小さい記述サイズの局所特徴量の対象となる局所特徴点が、類似する画像の大きい記述サイズの局所特徴量の対象となる局所特徴点の一部であれば、局所特徴量記述サイズの大小にかかわらず安定した検索が可能となる。そこで、第一の方法は、解析対象の画像を回転或いは縮小する等の画像処理を行う事による、局所特徴点の位置および局所特徴量の特徴間距離を求め、これを安定性の指標とする。局所特徴点の位置の変動が大きいと明らかに局所特徴量は変化してしまうため、まず、局所特徴点の位置が定められた閾値範囲内に収まる特徴点に絞る。そして、回転或いは縮小する等の画像処理の前後での局所特徴量の変化、即ち局所特徴量間距離を安定性の指標とする。そして、距離の小さな局所特徴点ほど安定性の高いものと見なし、これを距離の昇順でソートし、局所特徴量サイズが指定されたサイズに達するまで、安定度の高い局所特徴量を優先に記述してゆく。第二の方法は、スケールの大きいより縮小度の大きな画像から抽出した局所特徴点を安定性の高い局所特徴量とするものである。即ち、局所特徴点を抽出したスケールを安定性の指標とする。スケールの大きい画像はローパスフィルタの掛かった画像で有り、低域周波数成分からなる画像である。この画像では、マクロなエッジコーナーが残り、このエッジコーナーでは安定した局所特徴点の抽出が可能となる。しかも、マクロなエッジコーナーであるため、スケールが仮に異なったとしても、局所特徴量の変化も少ない。従って、局所特徴量の記述は、低周波から高周波へソート或いは帯域ごとに管理して記述する。第三の方法は、指定された局所特徴量の記述サイズを各スケール画像に対して配分し、各スケール画像に対し配分されたサイズ数の安定した局所特徴点の局所特徴量を記述する方法である。例えば、指定された局所特徴量の記述サイズをスケール画像の面積の比で分割し、各スケールでの予め決定した局所特徴量とすればよい。各スケール画像における、記述すべき局所特徴点およびその局所特徴量は、第一の方法に安定性の指標に基づくとよい。以上の局所特徴記述に関して述べてきたが、比較処理に関しては、先行技術の様に記述されている局所特徴量同士を記述サイズを意識せずにマッチングする事が可能である。なぜならば、安定度の観点でソートされた局所特徴量が記述してあるため、理論的には、トータルの記述サイズの大きな局所特徴量は、同記述サイズの小さい局所特徴量を包含する関係にあるからである。本発明によれば、信頼性の高い局所特徴点及び局所特徴量を得ることができる。以下、添付図面に従い本発明に係る実施形態を詳細に説明する。<第1の実施形態>本第1の実施形態では、入力画像から局所特徴点を検出し、その特徴点の回転拡縮不変の局所特徴量を抽出する局所特徴抽出を用いて説明を行う。先ず、第1の実施形態の画像処理装置の構成を図13に示す。この画像処理装置は、CPUCentralProcessingUnit1301、入力装置1302、出力装置1303、記憶装置1304、RAMRandomAccessMemory1305、ROMReadOnlyMemory1306、BUS1307からなる。これらの働きについて説明する。まずCPU1301は、各種データ処理のための論理演算、判断等を行い、またBUS1307で接続されている各構成要素の制御を行う。次に入力装置1302は、アルファベットキー、ひらがなキー、カタカナキー、句点等の文字記号入力キー、カーソル移動を指示するカーソル移動キー等のような、各種の機能キーを備えたキーボードが接続されている。また、画像の局所特徴量を算出する画像を入力するカメラ等の撮像素子、もしくは撮像素子によって記録した画像を保持し入力するためのHDD等のストレージも挙げられる。さらに、GUIGraphicalUserInterface上で画面の制御可能位置を指し示し、機能の選択指示などを行うポインティング機器、例えばマウスやスティックポインタ等も挙げられる。出力装置1303は、液晶パネル等の各種表示装置である。記憶装置1304には、入出力データや処理プログラム等、各種情報が格納される。これらのデータ及びプログラムを格納する記憶媒体としては、ハードディスク、CDROM、DVDROM、フラッシュメモリ等を用いることができる。RAM1305は、各構成要素からの各種データを一時的に記憶する際に用いられる。ROM1306は、本実施形態で実行される処理プログラム等の制御プログラムを記憶する。これら全ての構成要素はBUS1307で接続されている。図1は、第1の実施形態における画像比較装置の構成例を示すブロック図である。図1に示すように、画像比較装置100は、局所特徴点/局所特徴量抽出部102、局所特徴点座標/局所特徴量比較部104から構成される。そして、局所特徴点座標/局所特徴量比較部104から比較結果105を出力する。クエリ画像および比較先画像データベースに登録しようとする画像から、局所特徴点/局所特徴量抽出部102が局所特徴点/局所特徴量を抽出し、局所特徴点座標/局所特徴量比較部104がその比較結果105を出力する。通常、検索処理装置では、比較対象画像の局所特徴点座標/局所特徴量を予め求めておき、登録画像と対応付けて画像特徴データベースに登録しておく。そして、検索時には、クエリ画像から局所特徴点座標/局所特徴量を抽出し、それと画像特徴データベース内の登録画像の局所特徴点座標/局所特徴量とを比し、類似度の降順にソートして検索結果を求める。本第1の実施形態では比較処理までの説明に留め、検索処理の説明は割愛する。図2は、局所特徴点/局所特徴量抽出部102の内部の構成例を示すブロック図である。処理対象画像201として、クエリ画像又は比較先画像が入力画像となる。縮小画像生成部203が画像入力部202に入力された処理対象画像201から縮小画像を生成し、局所特徴点抽出部204が縮小画像から局所特徴点を抽出する。そして、局所特徴量算出部205がその局所特徴点の局所特徴量を算出し、局所特徴量検定部206が局所特徴量を検定する。そして検定結果207は、処理対象画像がクエリ画像101の場合には、クエリ画像の局所特徴点座標/局所特徴量となり、処理対象画像が比較先画像103の場合には、比較対象画像の局所特徴点座標/局所特徴量となる。[画像の局所特徴点座標/局所特徴量抽出処理]まず、画像を登録する際に行う各部の動作について説明する。図3は、局所特徴点/局所特徴量抽出部102における、画像の登録処理の制御手順を表すフローチャートである。S301で、図2に示す画像入力部202が比較先画像103を読み込み、S302で、当該比較先画像103から輝度成分を抽出する。そして、抽出した輝度成分に基づき輝度成分画像を生成し、生成した輝度成分画像を縮小画像生成部203に渡す。次に、S303では、縮小画像生成部203が画像入力部202から渡された当該輝度成分画像を倍率縮小率pに従って順次縮小することを繰り返し、オリジナルのサイズの画像から段階的に縮小した、オリジナルの画像を含めてn枚の縮小画像を生成し、当該縮小画像を局所特徴点抽出部204に渡す。ここで、倍率p及び縮小画像の枚数nは予め決められているものとする。図4は、縮小画像を生成する縮小画像生成部203の処理の一例を示す図である。図4に示す例は、倍率pが2の-1/4乗、縮小画像の枚数nが9の場合である。もちろん、倍率pが必ずしも2の-1/4乗で無くとも良い。図4において、符号401は画像入力部202から渡された輝度成分画像である。符号402は当該輝度成分画像から倍率pに従って再帰的に4回の縮小処理を行って得られた縮小画像である。そして、符号403は当該輝度成分画像から倍率pに従って8回縮小された縮小画像である。この例では、縮小画像402は画像入力部202から渡された輝度成分画像が1/2に縮小された画像となり、縮小画像403は画像入力部202から渡された輝度成分画像が1/4に縮小された画像となる。尚、画像を縮小する方法は何れの方法でも良く、第1の実施形態では、線形補間による縮小方法により縮小画像を生成するものとする。次に、S304では、局所特徴点抽出部204が、縮小画像生成部203から渡されたn枚の縮小画像の各々に画像の回転があってもロバストrobustに抽出されるような局所的な特徴点局所特徴点を抽出する。抽出した局所特徴点は局所特徴量算出部205に渡される。この局所特徴点の抽出方法として、第1の実施形態ではHarris作用素を用いる非特許文献1:CHarrisandMJStephensAcombinedcornerandedgedetectorInAlveyVisionConferencepages1471521988参照。具体的には、Harris作用素を作用させて得られた出力画像H上の画素について、当該画素及び当該画素の8近傍にある画素合計9画素の画素値を調べる。そして、当該画素が局所極大になる当該9画素の中で当該画素の画素値が最大になる点を局所特徴点として抽出する。ここで、当該画素が局所極大になったときでも、当該画素の値がしきい値以下の場合には局所特徴点として抽出しないようにする。尚、局所特徴点を抽出可能な方法であれば、上述のHarris作用素による特徴点抽出方法に限らず、どのような特徴点抽出方法でも局所特徴点抽出部204に適用可能である。次に、S305で、局所特徴量算出部205が、局所特徴点抽出部204から渡された局所特徴点の各々に画像の回転があっても不変となるように定義された特徴量局所特徴量を算出する。ここで抽出された局所特徴量は局所特徴点座標と関連付けられた上で局所特徴量局所特徴点座標207として出力する。この局所特徴量の算出方法として、第1の実施形態ではLocalJet及びそれらの導関数の組み合わせを用いるJJKoenderinkandAJvanDoornRepresentationoflocalgeometryinthevisualsystemRiologicalCyberneticsvol55pp3673751987参照。具体的には、以下の式1により局所特徴量Vを算出する。ただし、式1の右辺で用いている記号は、以下に示す式2から式7で定義される。ここで、式2右辺のGx,yはガウス関数、Ix,yは画像の座標x,yにおける画素値であり、は畳み込み演算を表す記号である。また、式3は式2で定義された変数Lのxに関する偏導関数、式4は当該変数Lのyに関する偏導関数である。式5は式3で定義された変数Lxのyに関する偏導関数、式6は式3で定義された変数Lxのxに関する偏導関数、式7は式4で定義されたLyのyに関する偏導関数である。尚、局所特徴量を算出可能な方法であれば、上述したような特徴量算出方法に限らず、どのような特徴量算出方法でも適用可能である。さて、第1の実施形態における検定はS306にて行うが、オリジナルの画像と、縮小画像に対し画像を45°回転の画像処理を行った画像の局所特徴量を比較し、局所特徴量の安定性を評価する。というのも、原理的に回転不変の局所特徴量ではあるが、実際に回転処理や画像縮小処理を行った場合に、ある程度の局所特徴量が起きてしまう事は事実であり、この影響は無視できないからである。本第1の実施形態では、画像の回転で行ったが、もちろん、拡大や縮小した画像を用いた検定を行っても良い。この場合、スケール画像の生成とは異なる率で拡大/縮小すると良い。S306の処理の詳細フローを示すのが、図5のフローチャートである。S306で、局所特徴量検定部206が、局所特徴量算出部205から渡された各々の局所特徴量の変動を求め変動の少ないものからソートする。S303では、縮小画像生成部203が画像入力部202から渡された当該輝度成分画像を倍率pに従って順次縮小し、縮小画像をn枚生成したが、S501にてこの縮小画を読み込み、S502にて画像を45度回転変換する。そして、S503にて、当該縮小画像を局所特徴点抽出部204に渡す。ただし、倍率p及び縮小画像の枚数nは前述した回転していない処理で用いた値と同じ値にする。次に、S504では、局所特徴点抽出部204が、縮小画像生成部203から渡されたn枚の縮小画像の各々に画像の回転があってもロバストに抽出されるような局所的な特徴点局所特徴点を抽出する。抽出した局所特徴点は局所特徴量算出部205に渡される。この局所特徴点の抽出方法として、回転していない処理と同様に、Harris作用素を用いる。具体的には、Harris作用素を作用させて得られた出力画像H上の画素について、当該画素及び当該画素の8近傍にある画素合計9画素の画素値を調べる。そして、当該画素が局所極大になる当該9画素の中で当該画素の画素値が最大になる点を局所特徴点として抽出する。ここで、当該画素が局所極大になったときでも、当該画素の値がしきい値以下の場合には局所特徴点として抽出しないようにする。尚、局所特徴点を抽出可能な方法であれば、上述のHarris作用素による特徴点抽出方法に限らず、どのような特徴点抽出方法でも局所特徴点抽出部204に適用可能である。S504にて、回転前後で、局所特徴点座標の変化が許容範囲に無いものを破棄し、局所特徴点を絞り込む処理を行う。この概念を表したものが図6である。左側のスケール画像が回転前のものである。右側のスケール画像は特徴抽出のために回転したスケール画像をわかり易い様に回転を戻したものである。点N2は回転前後で座標か許容範囲以内、その他の点N3、N4に関しては回転前後で対応する局所特徴点の位置が許容範囲を超えた変動をしている事をあらわす。次に、S505で、LocalJet及びそれらの導関数の組み合わせを用い、式1により局所特徴量を算出する。ここでは、回転変換した45°縮小画像から、S305で局所特徴量算出部205が算出した同じ方法により局所特徴量45°を算出する。ここで、これ以降の説明を簡略化するため、局所特徴量算出部205から渡された局所特徴量を局所特徴量0°と記載する。そしてS506では、局所特徴量0°と局所特徴量45°とのマハラノビス距離を比較し、着目する局所特徴点の局所特徴量0°とそのマハラノビス距離を対にして一時的に記憶する。尚、マハラノビス距離を計算するための分散共分散行列は、予め多くの局所特徴量から統計により求めたものを使用する。特徴間距離としてマハラノビス距離を用いる理由は、特徴空間における分布を考慮する事により、局所特徴量の安定性の評価の精度を高めるためである。もちろん精度を犠牲にしてでも高速化を望むのであれば、ユークリッド距離を用いても構わない。そして、S504にて絞り込んだ局所特徴点全てに対して、45度回転前後での局所特徴量変動をあらわすマハラノビス距離を求め、一次記憶する。そして、S507にて全ての局所特徴点座標/局所特徴量0°をマハラノビス距離の昇順でソートする。これは即ち、45°の回転があったとしても、局所特徴点の位置やその局所特徴量の変質に対する耐性の高い順番で局所特徴量を記憶している事に成る。図3のフローに戻り説明を続ける。S307にて、1画像当たりの局所特徴量記述局所特徴点座標局所特徴量サイズに制限を設ける場合に、上記ソートしてある局所特徴点と局所特徴量0°の対の累計サイズが特徴量記述サイズ個数の制限に達するまで記憶する。上記処理により、画像データベース不図示への登録対象の画像から局所特徴点/局所特徴量が抽出でき、これらが画像検索の際のインデックステーブルに登録されることになる。図12は、或る画像を画像データベースに登録する際に抽出された、その画像を特定するための局所特徴点/局所特徴量のインデックステーブルの一例を示している。図示の如く、画像データベースに登録する画像データへのポインタ或いは、登録する画像が画像ファイルである場合には、その画像ファイルのパス付きファイル名でも良いを格納するエリアと、スケール毎縮小画像毎のソートされた局所特徴点/局所特徴量のデータ列を格納するエリアで構成されている。実施形態ではn=9の例を示しているので、スケールの種類は図示の如く、SC1乃至SC9の9種類である。そして、fi1乃至fiMiはスケールを示し1乃至9の値、Mは予め設定された値で示されるソート済みの局所特徴点/局所特徴量で構成される。このように、第1の実施形態では、縮小画像を回転変換し局所特徴量の変動の程度を予め検定し、限られた局所特徴量のサイズ、より変動の少ない局所特徴量から順番に記憶する事が出来る。尚、第1の実施形態では、検定時の回転角度を45°として、45°回転検定だけを実施する場合を説明したが、45°回転検定だけでなく、例えば0°~90°までt°単位に回転検定を実施しても良い。ただし、tは予め定められた0<t<90を満たす値である。[画像間の類似度演算の方法]次に、図7のフローチャートを用いて類似指標算出方法として、局所特徴点座標/局所特徴量比較部104を説明する。ここで、クエリ画像101の局所特徴量をVq、局所特徴量に関連付けされている局所特徴点をQ,座標をQx’,y’と表わす。また、比較先画像103に存在する局所特徴量をVs、局所特徴量に関連付けされている局所特徴点をS、座標をSx,yとする。まず、S701で、最終投票数を表す変数VoteMaxを0に初期化する。次に、S702で、VqとVsとの特徴量間距離を全ての組み合わせについて計算し、最短距離対応点リストを作成する。即ち、閾値Tvを用い、計算した特徴量間の距離が閾値以下となり、かつ、最短距離となるようなVqとVsとの組み合わせ対応点を抽出し、最短距離対応点リストに登録する。これ以降、最短距離対応点リストに登録されたk番目の対応点について、当該対応点の局所特徴量をそれぞれVqkとVskと記載する。更にVqkとVskに対応付けられている局所特徴点をそれぞれQk、Sk、座標をQkx’k,y’k、Skxk,ykなどと添え字を合わせて記載する。またS702で作成された最短距離対応点リストに登録された対応点の組数をm組とする。次に、S703で、類似度算出処理の反復カウント数を表す変数Countを0に初期化する。次に、S704で、反復カウント数Countが予め定められた最大反復処理回数Rnを超えていないを判定する。ここで、超えている場合はS618へ進み、最終投票数VoteMaxを類似指標として出力し、この処理を終了する。また、S704で、反復カウント数Countが最大反復処理回数Rnを超えていない場合はS705へ進み、投票数を表す変数Voteを0に初期化する。次に、S706で、当該最短距離対応点リストから対応点の組の座標をランダムに2組抽出する。ここで、これらの座標をQ1x’1,y’1、S1x1,y1及びQ2x’2,y’2、S2x2,y2と記載する。次に、S707で、抽出したQ1x’1,y’1、S1x1,y1及びQ2x’2,y’2、S2x2,y2が式8に示す変換を満たしていると仮定し、式8中の変数a~fを求める。ただし、図6に示すS607では、変数a~dで構成される行列をMで示し、変数e~fで構成される行列をTで示している。ここで、第1の実施形態では、簡略化のため、相似変換だけを考える。このとき、上記式8は以下の式9のように書き換えられる。このとき、変数a、b、e、fはx’1、y’1、x1、y1、x’2、y’2、x2、y2を使って式10から式13で表される。次に、S708で、上述のS705で当該最短距離対応点リストからランダムに抽出された2組の点以外の点を選択するために、対応点選択変数kを3に初期化する。そして、S709で、対応点選択変数kが当該最短距離対応点リストに登録されている対応点の組数mを超えていないかを判定する。ここで、超えている場合はS715へ処理を移すが、これについては後述する。S709における判定で対応点選択変数kが当該最短距離対応点リストに登録されている対応点の組数mを超えていない場合はS710へ処理を移す。このS710では、上述のS605で当該最短距離対応点リストからランダムに抽出した2組の点S1x1,y1及びS2x2,y2以外の点を当該最短距離対応点リストから抽出する。第1の実施形態では、抽出された点をSkxk,ykと記載する。次に、S611で、Skxk,ykが式9を使って移される座標Sk’x’k,y’kを求める。その後、S712では、座標Sk’x’k,y’kと座標Qkx’k,y’kとの幾何学的距離をユークリッド距離で計算し、当該ユークリッド距離がしきい値Td以下であるか否かを判定する。当該ユークリッド距離がしきい値Td以下の場合はS713へ進み、投票数Voteをインクリメントし、S714へ処理を移す。また、当該ユークリッド距離がしきい値Tdより大きい場合は、何もせずにS714へ処理を移す。このS714では、対応点選択変数kをインクリメントし、S709に戻り、対応点選択変数kが当該最短距離対応点リストに登録されている対応点の組数mを超えるまで、上述の処理を繰り返す。次に、S709で、対応点選択変数kが当該最短距離対応点リストに登録されている対応点の組数mを超えた場合の処理であるS615を説明する。S715では、投票数Voteの値と最終投票数VoteMaxの値とを比較し、投票数Voteの値が最終投票数VoteMaxの値よりも大きい場合にはS616へ処理を移す。このS716では、最終投票数VoteMaxの値を投票数Voteの値で置き換えた後、S717で反復カウント数Countをインクリメントし、上述のS704に処理を戻す。また、S715で、投票数Voteの値が最終投票数VoteMaxの値以下の場合にはS717へ処理を移し、反復カウント数Countをインクリメントし、上述のS704に処理を戻す。尚、第1の実施形態における類似度の算出方法の説明では、相似変換だけを考えて説明したが、アフィン変換などその他の幾何学変換についても、S707でそれぞれに応じた変換行列を求めることにより、対応可能である。例えば、アフィン変換の場合には、まずS706で、ランダムに選択する対応点の組の座標数を3とする。次に、S707で、式9ではなく式8を使うこととし、S706で選択した3組の対応点合計6点を使って変数a~fを求めれば良い。また、第1の実施形態における類似度の算出方法の説明では、S719で類似指標として最終投票数VoteMaxを出力する方法を説明したが、本発明はこれに限らず、他の類似指標を計算するようにしても良い。即ち、S703以降の処理を行わずに、S702で作成された最短距離対応点リストに登録された対応点の組数mを類似指標として出力するようにする。これにより、一定の検索精度を保ちつつ、検索速度を向上させることが可能になる。以上のように、第1の実施形態の画像検索装置では、クエリ画像および登録画像から算出された局所特徴量の変動を検定し、変動の大きい局所特徴量を検索時に使わないように構成した。これにより、コンパクトな局所特徴量の記述で検索精度の低下を抑制することが可能になる。尚、上記実施形態では、画像自身を回転させた後に局所特徴量を算出し、算出した当該局所特徴量を使って検定した。しかし、局所特徴量算出部205で局所特徴量算出に使ったフィルタを45°回転させたフィルタ45°回転フィルタを予め準備しておき、当該45°回転フィルタを使って検定を実施しても良い。ここで、フィルタを45°回転させるとは、所定のフィルタ係数に対して幾何学的に回転変換を施すことを意味する。また、類似度演算処理において、クエリ画像および登録画像の局所特徴量の記述サイズ、即ち記述する有効な局所特徴点の個数に関しては、異なる大きさでも良い。なぜならば、安定度の観点でソートされた局所特徴量が記述してあるため、理論的には、記述サイズの大きな特徴量は、記述サイズの小さい特徴量を包含する関係にあるからである。従って、本第1の実施形態では、記述する有効な局所特徴点の個数に違いが有っても、その全てを距離演算に用いて良い。他方、高速に処理する場合には、クエリ画像と登録画像の記述特徴点のうち、局所特徴点の数の少ない方の画像の特徴点数Nを、特徴点の多い方の局所特徴点の距離演算の参照範囲個数とする方法もある。というのも、安定性の良いものからソートして記述しているので、対応する局所特徴点を比較する可能性が高く、参照個数を減らす事で高速化が図れるからである。<第2の実施形態>次に、図面を参照しながら本発明に係る第2の実施形態を説明する。第2の実施形態では、安定した局所特徴点局所特徴量の多くが低周波成分から得られる特徴があり、高周波成分は画像ノイズの影響が大きいという経験則を用いる。安定した局所特徴点を取りやすい低周波成分の画像、即ち縮小度合いの大きな画像から優先して採用する。というのも、縮小するぼかすほど、局所特徴点の数が減り、マクロな骨太な局所特徴点が残るからである。図4に、局所特徴量を抽出する際の段階的な縮小画像であるスケール画像を示す。スケール番号が大きいほど縮小処理を行っている。そこで、縮小度合いの大きい画像から得られる局所特徴点ほど安定性が高いという経験則を用い、スケール番号の大きさを安定性の高さ指標として処理を行う。即ち、第1の実施形態と同様、Harris作用素を作用させて得られた出力画像H上の画素について、当該画素及び当該画素の8近傍にある画素合計9画素の画素値を調べる。そして、当該画素が局所極大になる当該9画素の中で当該画素の画素値が最大になる点を局所特徴点として抽出する。ここで、当該画素が局所極大になったときでも、当該画素の値がしきい値以下の場合には局所特徴点として抽出しないようにする。局所特徴量の抽出に関する処理フローは第1の実施形態で用いた図3で同じであるが、S306およびS307の処理内容が異なる。本第2の実施形態のS306とS307を合わせた処理内容を図8のフローを元に処理の説明を行う。S801にて、第1の実施形態の時と同様、座標位置変化の許容範囲の局所特徴点に絞り込む図9参照。S802にいて、スケール番号に相当するカウンタの初期値を8にセットする。図4でスケールを9個用いる例なので、局所特徴点の得られるスケール番号8から処理が開始するためである。S803では、初期化のため処理結果をクリアする。S804では、スケール番号1では局所特徴点がない為に、N<2で処理を終了させる。他方、N<2を満たさなければS805にて、スケール番号Nの局所特徴点数Mを求める。カウンタKを0にリセットし、S807にてK<Mを満たすかどうか判断し、満たさない場合、即ち、スケール番号Nの局所特徴点数M分記述した場合にはS808にてNをディクリメントし、再びS804に戻る。S807にて、K<Mを満たす場合にはS809にて、次に記述した場合に許される記述サイズが超えるか判断し、超える場合には処理を終了する。他方、超えない場合には、S810にてスケールNのk1番目の局所特徴量処理結果に追記する。尚、スケールNにおける局所特徴点の並び順番K1に関しては、なるべく画像に対して均等な分布になるようにランダムにすれば良い。そして、S811にてKをインクリメントし、S807に戻る。上記処理により、S306の処理結果は、スケール番号の大きな局所特徴点から順番に、局所特徴点座標/局所特徴量が記述する事になる。この際に、画像特徴データベース107に登録する1画像のトータルの局所特徴量記述サイズに制限を設ける場合には、スケール番号の大きい局所特徴点から局所特徴点座標/局所特徴量の記述を行う。そして、累計サイズがトータルの局所特徴量記述サイズの制限に達するまで記憶し、画像の登録処理が終了する。尚、検索処理に関しては、第1の実施形態と同様に行う。以上のように、本第2の実施形態によれば、画像の登録処理時に安定した局所特徴点を取りやすい低周波成分の画像、即ち縮小度合いの大きな画像から優先して登録するように構成したので、サイズ制限のある局所特徴量記述において、検索精度の低下を抑制することが可能になる。<第3の実施形態>局所特徴量の抽出に関する処理フローは第1の実施形態で用いた図3で同じであるが、S306およびS307の処理内容が異なる。本第3の実施形態のS306とS307を合わせた処理内容を図11のフローを元に処理の説明を行う。複数のスケール帯域において、採用する局所特徴点数をスケールに依存して配分し、そのスケール内での特徴の安定度の高いものから記述を行うものである。本第3の実施形態では、面積比に応じて総局所特徴点数Sを割り振るが、もちろん固定数でも良い。総数S個の局所特徴量の記述をするにあたって、各スケールに対して総記述数をS、縮小比率をrとした場合、図10の様に、スケールが小さいほど少ない局所特徴点数Sを決定すれば良い。S1101にて、第1の実施形態の時と同様、座標位置変化の許容範囲の局所特徴点に絞り込む。S1102にいて、スケール番号に相当するカウンタの初期値を8にセットする。図4で示すように実施形態ではスケールを9個用いる例を説明しており、局所特徴点の得られるスケール番号8から処理が開始するためである。S1103では、図10に従って、各スケールに対する局所特徴記述を行う目標数M[8]を決定する。配列となっているのは、配列要素がスケールに対応しているからである。但し、M1、及び、M8の値は、局所特徴点が原理的に無いので、0である。S1104にて初期化のため処理結果をクリアする。S1105では、スケール番号1では局所特徴点がない為に、N<2で処理を終了させる。他方、N<2を満たさなければS1106にて、スケールNの局所特徴数Aと、記述目標数数MNの小さいほうをLIMへセットする。S1107にてカウンタKを0にリセットし、S1108にてK<LIMを満たすかどうか判断し、満たさない場合、即ち、スケール番号Nの局所特徴点数M記述した場合にはS1109にてNをディクリメントし、再びS1105に戻る。S1108にて、K<LIMを満たす場合にはS1110にて、記述した場合に許される記述サイズを超えるか判断し、超える場合には処理を終了する。他方、超えない場合には、S1111にてスケールNのK1番目に安定度が高い局所特徴量を追記する。この局所特徴量の安定性に関しては第1の実施形態1で述べた画像回転に対しても変動の少ない局所特徴点座標/局所特徴量を各スケール毎に管理し、そのスケール内でマハラノビスの距離の小さいものからソートして置けばよい。そして、S1112にてKをインクリメントし、S1108に戻る。この様に、第3の実施形態では、各スケール帯域に対し適切な局所特徴量サイズを分配し、各スケールにおいて回転させたフィルタを使用して局所特徴量の変動の程度を予め検定し、変動の少ないものから各スケールに割り振られたサイズの局所特徴量を記述する。これにより、安定した局所特徴点局所特徴量の多くが、低域成分から得られる特徴で有るという経験則と、画像の変化に対する変動が少ないという両面から選ばれた局所特徴量を所定のサイズで記述する事が可能となる。もちろん、各スケール内の特徴点から無作為に選んだ場合でも、全スケールの局所特徴量から無作為にS個を選んだ場合よりも、安定性の高い局所特徴量の記述が可能である。尚、検索処理に関しては、第1の実施形態と同様に行う。その他の実施例また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェアプログラムを、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータまたはCPUやMPU等がプログラムを読み出して実行する処理である。もちろん、本発明はデジタルカメラなどの撮像デバイスにおいても、撮像した画像データに対し、上述した実施形態の機能を実現するソフトウェアプログラムを、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータまたはCPUやMPU等がプログラムを読み出して実行する事が可能である。"