JP2017224337A - 勾配ヒストグラムに基づいて画像記述子を変換する方法および関連する画像処理装置 - Google Patents

勾配ヒストグラムに基づいて画像記述子を変換する方法および関連する画像処理装置 Download PDF

Info

Publication number
JP2017224337A
JP2017224337A JP2017154781A JP2017154781A JP2017224337A JP 2017224337 A JP2017224337 A JP 2017224337A JP 2017154781 A JP2017154781 A JP 2017154781A JP 2017154781 A JP2017154781 A JP 2017154781A JP 2017224337 A JP2017224337 A JP 2017224337A
Authority
JP
Japan
Prior art keywords
histogram
gradient
descriptor
elements
transformed
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2017154781A
Other languages
English (en)
Other versions
JP6429134B2 (ja
Inventor
パスカラキス、スタブロス
Paschalakis Stavros
ボーバー、ミロスラウ
Bober Miroslaw
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sisvel Technology SRL
Original Assignee
Sisvel Technology SRL
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sisvel Technology SRL filed Critical Sisvel Technology SRL
Publication of JP2017224337A publication Critical patent/JP2017224337A/ja
Application granted granted Critical
Publication of JP6429134B2 publication Critical patent/JP6429134B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T1/00General purpose image data processing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/40Analysis of texture
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/40Analysis of texture
    • G06T7/41Analysis of texture based on statistical description of texture
    • G06T7/44Analysis of texture based on statistical description of texture using image operators, e.g. filters, edge density metrics or local histograms
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/42Global feature extraction by analysis of the whole pattern, e.g. using frequency domain transformations or autocorrelation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/46Descriptors for shape, contour or point-related descriptors, e.g. scale invariant feature transform [SIFT] or bags of words [BoW]; Salient regional features
    • G06V10/462Salient features, e.g. scale invariant feature transforms [SIFT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/50Extraction of image or video features by performing operations within image blocks; by using histograms, e.g. histogram of oriented gradients [HoG]; by summing image-intensity values; Projection analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/74Image or video pattern matching; Proximity measures in feature spaces
    • G06V10/75Organisation of the matching processes, e.g. simultaneous or sequential comparisons of image or video features; Coarse-fine approaches, e.g. multi-scale approaches; using context analysis; Selection of dictionaries
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20004Adaptive image processing
    • G06T2207/20012Locally adaptive
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20021Dividing image into blocks, subimages or windows

Abstract

【課題】複数のヒストグラムビンhiを含む勾配ヒストグラムに基づく画像記述子を、複数の値vjからなるセットを含む変換された勾配ヒストグラム記述子に変換する方法を提供する。【解決手段】方法は、複数の値vjからなるセットを算出するべく少なくとも1つの基準を用いる段階を含む。少なくとも1つの基準は、ヒストグラムビンhiの配置基準を含む。【選択図】図8

Description

この発明は、ロバストで、判別可能で、スケーラブルで、かつ、簡潔な画像記述子の算出に関係する。
さらに特には、本発明は、画像の関心点周りのローカル領域内で部分領域の勾配のヒストグラムを計算することにより算出される画像記述子に言及する。
画像記述子を、少しの例を挙げれば、対象認識、内容に基づく画像取得、および、画像登録を含む、多くのコンピューター画像アプリケーションに広く適用できることを見出した。この種類の画像記述子の最も広く知られている例の1つは、Scale Invariant Feature Transform(SIFT)記述子である。
手短に言えば、SIFTのやり方では、以下のようにローカル画像記述子が形成される。最初に、縮尺および配向に対して不変である安定した画像キーポイントを識別して、かつ、位置付けるべく、複数の画像の縮尺および位置に亘る検索が実行される。そして、各キーポイントに対して、1または複数の優勢な配向がローカル画像勾配に基づいて決定される。各キーポイントの割り当てられた配向、縮尺および位置に関して、次のローカル記述子の算出を実行できるようにする。これによって、これらの変換に対する不変性を達成する。
そして、以下のようにキーポイント周りのローカル画像記述子が形成される。最初に、キーポイント周りの領域において、勾配の大きさおよび配向の情報が画像サンプルポイントで算出される。そして、これらのサンプルが、n×nの部分領域の全体の内容を概括する配向ヒストグラムへと蓄積される。
実例としてのみ、キーポイント記述子の例が図1aおよび1bに示されている。図1aは、ローカル領域Rの4×4の部分領域SRへの細分を示し、図1bは、各配向ヒストグラムhに対して、そのヒストグラムエントリの大きさに対応する各矢印の長さと共に、配向の360°範囲の8つのビンhへの細分を示している。
従って、図1aで示されたようなローカル画像記述子は、4×4×8=128の要素を有している。SIFTのやり方は、International Journal of Computer Vision, 60, 2 (2004), 91‐110ページ内での、David G. Loweによる「Distinctive image features from scale−invariant keypoints」でさらに詳しく示されている。
安定した画像キーポイントの検出に対する異なるメカニズムと、キーポイント周りのローカル領域の細分に対する異なる手法と、部分領域勾配ヒストグラムの算出に対する異なる手法とを用いた、SIFTのやり方の多くの代案およびバリエーションが存在する。
例えば、図2aおよび図2bは夫々、SIFTのやり方で使用されているデカルト空間細分の代案として、IEEE Transactions of Pattern Analysis and Machine Intelligence 27(10):1615‐1630内でのK. MikolajczykおよびC. Schmidによる「A performance evaluation of local descriptors」で説明されているGradient Location Orientation Histogram (GLOH)、IEEE Transactions of Pattern Analysis and Machine Intelligence 27(10):1615‐1630内でのChandrasekhar他、MikolajczykおよびC. Schmidによる「A performance evaluation of local descriptors」で説明されているUncompressed Histogram of Gradients (UHoG)、および、International Journal on Computer Vision, Vol. 94, No. 5, May 2011内でのChandrasekhar他による「Compressed Histogram of Gradients: A Low−Bitrate Descriptor」で説明されているUncompressed Histogram of Gradients (UHoG)のような、他の技術のlog−polar空間細分特性を示している。
他の例として、図3aおよび図3bは、勾配配向の360°範囲をSIFTのやり方で使用されるビンへと細分することの代案として、勾配のXおよびY成分の2次元空間をUHoGの特性であるビンへと細分することに基づいた、勾配ヒストグラムの算出へのアプローチを示している。
本発明がロバストで、判別可能で、スケーラブルかつ簡潔な画像記述子の算出を実行することに基づく画像記述子を生成する技術の例として唯一、上述の先行技術が本明細書で考慮されている。
そのような画像記述子が、これより前に議論したように、多くのコンピューター画像アプリケーションにおいて広く適用され得ることを見出したけれども、これらのストレージおよび伝送のコストは、これらのサイズによってバイトで規定されるように、特定アプリケーション領域で一般に高いものと考えられる。何故かと言うと、画像におけるキーポイントに対するローカル画像記述子のサイズが相対的に低いものであってもよいけれども、画像記述子全体が数百のそのようなキーポイントおよび関連したローカル記述子を備えるからである。これは、画像記述子全体が実際の画像から抽出されるJPEGの圧縮版と同程度のサイズを有し得ることを意味している。
記述子サイズのこのレベルが問題であると考えられる1つのそのようなアプリケーション領域は、移動端末を使用した画像検索である。異なる構成がこのアプリケーション領域で実現可能であるけれども、1つの典型的な構成は、携帯電話等の移動端末クライアントによる関心対象の画像のキャプチャ、そのクライアントによる画像記述子の自動抽出、その画像記述子を処理してその関心対象に関する識別または追加情報等の適切な応答を提供するであろうサーバへと無線通信ネットワークを介してその画像記述子を伝送すること、および、そのクライアントへの前述の応答の返送を必要とする。従って、無線ネットワーク上でクライアントからサーバへと送信される情報量の最小化が望ましいことは明らかである。そのようなアプリケーションの恩恵を得るべく、そのような画像記述子の圧縮における、かなりの量の開発が行われてきた。
キーポイント記述子に基づいた勾配のヒストグラムの圧縮への最も簡潔なアプローチは、ヒストグラムビンの値のスカラー量子化によるものであり、それは、夫々のビンの値の表現において個々に用いられるビット数を減らすことを意味している。実際には、記述子の弁別力を大幅に損なわずにとても高い圧縮率を達成することは困難であるので、このアプローチは一般に用いられていない。例えば、ビンごとに8ビットでSIFT記述子ヒストグラムビンを符号化することが一般に用いられているが、バイトでのサイズが無線ネットワーク上で伝送するには大きすぎると一般に考えられている画像記述子をもたらす結果となっている。他方で、ビンごとに例えば1または2だけのわずかなビットへのスカラー量子化は、画像記述子の弁別力を危険にさらすことが判明してきた。
それゆえ、より複雑な圧縮スキームが提案されてきた。そのようなスキームのレビューは、Proceedings of International Conference on Pattern Recognition (ICPR), Istanbul, Turkey, August 2010内でのV.Chandrasekhar他による「Survey of SIFT compression schemes」で示されている。
手短に言えば、代表的な有限数のベクトル重心の1つに変換することによってビンの値が併せて量子化されるベクトル量子化を主要テーマとするスキームは、特に有名であって、ツリー構造化および生成物ベクトル量子化等の様々な形態で調査されてきた。そのようなアプローチの欠点は、それらが相対的に高い計算の複雑性を必要とし、複数の重心のストレージ用に数百キロバイトから数メガバイトあるいはそれ以上の極めて多くの必要メモリを必要とすることであり、複数の重心の数は数千から数百万に及び、そして複数の重心の決定は計算上複雑なトレーニングフェーズを必要とする。
全ての可能な入力ベクトルを含む空間内でタイプの一様な格子を形成し、与えられた入力ベクトルのいずれに対しても、最も近いタイプの指標によって符号化することにより、ビンの値が再度併せて量子化されるタイプ・コーディングを主要テーマとするスキームも、徹底的に調査されてきた。そのようなアプローチの必要メモリは、ベクトル量子化アプローチに比べて減らされるが、高い圧縮率での認識性能の観点で、合成圧縮記述子が量子化されたベクトル記述子にあまり匹敵しないことも判明してきた。全体的にみれば、タイプ・コーディングに関連する計算コストは簡単なスカラー量子化に比べて著しく高い。
他の圧縮スキームは、例えば128−dimensional SIFTキーポイント記述子のようなキーポイント記述子上で、PCA等として知られている次元低減方法を利用しており、合成次元のスカラー量子化がこれに続く。そのようなアプローチにある重要な問題は、それらが高い計算の複雑性とオーバートレーニングのハイリスクとを必要とすることである。
要するに、記述子に基づく勾配のヒストグラムの圧縮、および、ロバストで、判別可能で、スケーラブルかつ簡潔な画像記述子の生成に対する既存のアプローチは、特定の欠点を示している。
記述子要素のスカラー量子化等の簡単なアプローチは、とても低い計算の複雑性と必要メモリとの恩恵を有しているが、高い圧縮率における記述子の弁別力を危険にさらすことが判明してきた。
より複雑なアプローチは、高い圧縮率でより良い性能を達成することが示されてきたが、異なる欠点を被る。ベクトル量子化アプローチは、計算の複雑性および必要メモリを大幅に増大させてきた。タイプ・コーディングアプローチは、増大した複雑性を必要とし、ベクトル量子化アプローチの必要メモリによる負担を受けないながらも、そのようなアプローチと比較して性能が低いことも判明してきた。さらに、ベクトル量子化およびタイプ・コーディングのいずれのアプローチも、圧縮された領域における次元低減にはあまり適していない。PCA等として知られている次元低減技術に基づくアプローチも使用されてきたが、高い計算の複雑性およびオーバートレーニングのハイリスクを被りもしている。
それゆえ、本発明の目的は、画像記述子のサイズを低減できるようにする勾配ヒストグラムに基づいた画像記述子を変換するための方法を示すことである。
本発明のさらなる目的は、特に移動端末アプリケーションにおける無線ネットワークを介した格納および伝送のコストを低減できるようにする勾配ヒストグラムに基づいた画像記述子を変換するための方法を示すことである。
本発明のさらなる目的は、画像キーポイント周りの部分領域に関する複数の勾配ヒストグラムによって構成された画像記述子のサイズを低減できるようにする勾配ヒストグラムに基づいた、画像記述子を変換するための方法を示すことである。
本発明のこれらの、および、他の目的は、本説明の不可欠な部分である添付の特許請求の範囲でクレームされているように、勾配ヒストグラムに基づいた画像記述子を変換するための方法および関連する画像処理装置を通して達成される。
合成において、本発明は、勾配の前述のヒストグラムの変換に基づいた勾配のヒストグラムを用いる画像記述子から、ロバストで、判別可能で、スケーラブルかつ簡潔な画像記述子を算出することに関係しており、前述の変換は、その中に含まれた顕著でロバストな情報を、分布の形状、および、それらのビンの値の間における関係の形態で捉えている。より具体的には、本発明は、ヒストグラムビン間の相関に対応する値を算出することで、勾配のヒストグラムを変換している。
本発明による方法は、画像キーポイント周りの部分領域に関する勾配のヒストグラムを処理する異なる手法に基づいた、ロバストで、判別可能でスケーラブルかつ簡潔な画像記述子の計算にも関係している。
本発明のさらなる特徴は、添付の特許請求の範囲において提示されており、本説明の不可欠な部分として意図されている。
上記の目的は、特に添付の図面を参照することで、勾配ヒストグラムに基づいた画像記述子を変換するための方法の以下の詳細な説明から、より明らかになるであろう。
先行技術のキーポイント記述子の例を示している。 先行技術のキーポイント記述子の例を示している。 GLOHの先行技術において用いられている空間細分スキームを示している。 UHoGの先行技術において用いられている空間細分スキームを示している。 UHoGの技術の勾配ヒストグラムの算出へのアプローチを示している。 UHoGの技術の勾配ヒストグラムの算出へのアプローチを示している。 本発明によるキーポイント記述子の例を示している。 本発明によるキーポイント記述子の例を示している。 画像キーポイント周りの部分領域に対応する勾配のヒストグラムを処理する方法を示している。 画像キーポイント周りの部分領域に対応する勾配のヒストグラムを処理する方法を示している。 画像キーポイント周りの部分領域に対応する勾配のヒストグラムを処理する方法を示している。 本発明の実施形態のフロー図を表している。 本発明の実施形態のフロー図を表している。 本発明の実施形態のフロー図を表している。 本発明の実施形態のフロー図を表している。 本発明の実施形態のフロー図を表している。 本発明の実施形態のフロー図を表している。 図2bの空間細分スキームを使用することによる、本発明による画像記述子の処理を示している。 図3aの勾配ヒストグラム細分スキームを使用することによる、本発明による画像記述子の処理を示している。 本発明による方法を実行するのに適した画像処理装置を示している。
例として、いかなる制限もせずに、画像記述子、特に、ローカル領域Rを4×4の部分領域SRへと細分する例を示している図4a、および、部分領域に対する8つのビンの勾配ヒストグラムの例を示している図4bを備える、図4aおよび図4bで示されるようなSIFT画像記述子から、ロバストで、判別可能で、スケーラブルかつ簡潔な画像記述子を算出するための、本発明の複数の態様が示される。
hは図4bの勾配ヒストグラムを表すものとし、ビンhにおいて、iはビンの指標であり、以下の関係式(1)で表される。
i=0...n−1・・・関係式(1)
本発明による勾配ヒストグラム変換は、勾配ヒストグラムビン間の相関に対応する値vを計算する。
複数のビンは、その間で相関に対応する値が算出されるものであり、図4bの勾配ヒストグラムの場合に、それらの角度分離へと変換する、それらの位置または近接に従って選択されてもよい。
勾配ヒストグラムビンの間の相関に対応する値の例として、以下の数式(1)で示されるように、隣接するビンの差は、ヒストグラムの全体的な形状および極値の位置を捉えている。
=h−h
=h−h
=h−h
=h−h
=h−h
=h−h
=h−h
=h−h・・・数式(1)
以下の数式(2)で示されるような90°の、または、勾配ヒストグラム特性およびビン細分性が許容するように、より概して90°と同程度の角度差を備える勾配に対応するビン間の差が、互いにほぼ直角の勾配間の相関についての情報を捉えている。
=h−h
=h−h
=h−h
=h−h
=h−h
=h−h
=h−h
=h−h・・・数式(2)
以下の数式(3)で示されるような180°の、または、勾配ヒストグラム特性およびビン細分性が許容するように、より概して180°と同程度の角度差を備える勾配に対応するビン間の差が、水平配向内および垂直配向内等の配向内での勾配間における相関についての重要なハイレベルの情報を捉えている。
=h−h
=h−h
=h−h
=h−h・・・数式(3)
以下の数式(4)で示される関数は、2よりも多い勾配ビンを結合し、どちらか一方の方向における垂直勾配に対する、どちらか一方の方向における水平勾配等の、特定の配向間の勾配の相関についての情報を捉えている。
=(h+h)−(h+h
=(h+h)−(h+h)・・・数式(4)
以下の数式(5)で示される関数は、全体のヒストグラムの形状の代表となる判別可能な値を作り出すべく、勾配ヒストグラムの全てのビンを結合している。
=(h+h+h+h)−(h+h+h+h
=(h+h+h+h)−(h+h+h+h)・・・数式(5)
その間で相関に対応する値が算出されるビンは、異なる抽出条件に従って選択されてもよい。例えば、その間で相関に対応する値が算出されるビンは、合成値の弁別力を増大するべく選択されてもよい。そのような関数は、例えば、以下の数式(6)で示されるような関数の形をとってもよい。
=(h+h+h+h)−(h+h+h+h
=(h+h+h+h)−(h+h+h+h
=(h+h)−(h+h
=(h+h)−(h+h)・・・数式(6)
他の例として、その間で相関に対応する値が算出されるビンは、ランダム選択処理を通して選択されてもよい。そのような関数の例は、以下の数式(7)で示される。
=2h−h−h
=h+h+h−3h・・・数式(7)
ヒストグラム変換の他の例として、以下の数式(8)で示される関数のセットは可逆変換を構成しており、それによって、vが全ヒストグラムカウントを代表する値となっている。
=h−h
=h−h
=h−h
=h−h
=h−h
=h−h
=h−h
=h+h+h+h+h+h+h+h・・・数式(8)
数式(1)〜(8)の関数は、ビン間の和または差、または、ビンの和の間の差として算出される勾配ヒストグラムビン間の相関を示しているが、これは限定的なことではなく、ビン間の比、ビンの和の間の比、ビンの生成物間の比、または、ビンの生成物間の差等のように、異なる線形または非線形演算処理が使用されてもよい。
一般に、勾配ヒストグラムは、勾配ヒストグラムビン間の相関に対応する値vを計算することによって変換されてもよく、前述の勾配ヒストグラムビンの選択は、それらの近接度や合成値の弁別力等の1または複数の基準、または、その基準の組合せ次第である。
さらに、それらの近接に従ってビンが選択される場合、勾配ヒストグラムは、例えば、数式(1)のような隣接するビン間のみの相関や、数式(3)のような180°に分離されたビン間のみの相関等に基づいて、単一角度分離を示す勾配ヒストグラムビン間の相関に対応する値vを計算することにより、変換されてもよい。代替的に、勾配ヒストグラムは、例えば、数式(1)で示されるようないくつかの相関や、数式(2)で示されるようないくつかの相関等に基づいて、複数の角度分離を示す勾配ヒストグラムビン間の相関に対応する値vを計算することにより、変換されてもよい。全体的にみれば、以下の関係式(2)によって互いに関係付けられ得る勾配ヒストグラムの次元nと比較して、合成変換記述子の次元kと異なるタイプの情報の混合に対して特殊なタイプの情報の量の均衡を保っているので、変換選択処理はかなり重要である。
k=n
OR
k<n
OR
k>n・・・関係式(2)
下文に、勾配ヒストグラムhの変換から算出される、以下の関係式(3)を満たす値vは、k個の要素の変換された勾配ヒストグラム記述子vとして一括して称されるであろう。
j=0...k−1・・・関係式(3)
これより前に見られるように、単一の画像記述子Hは、画像キーポイントKP周りのローカル領域Rの部分領域SRにそれぞれ対応する、勾配hの複数のヒストグラムを構成してもよい。図4aは、ローカル領域Rを4×4の部分領域SRへと細分する例を示しており、勾配hのヒストグラムによってそれぞれ描かれている。そのような画像記述子は、各部分領域SRの勾配ヒストグラムを、上述したように変換された勾配ヒストグラム記述子vpへと個々に変換することで、新しい画像記述子Vへと変換されてもよい。
上述の勾配ヒストグラム変換処理は、簡潔な画像記述子を生成する点で多くの利点を有する。最初の利点は、変換の適切な選択を通して、上述の処理により、例えば8つの勾配ビンをこれより前で示された数式(3)および(4)を組み合わせた6つの要素と置き換えることによって、最初の勾配ヒストグラムより少ない要素を使用する勾配ヒストグラムの、顕著な形状特性および判別可能な情報を符号化できることであり、その結果、画像コレスポンデンスを確立、または、検証する点で依然としてとても成功である、より小さな画像記述子をもたらし、かつ、高い認識性能を達成する。
他の利点は、最初の勾配ヒストグラムと異なり、変換された勾配ヒストグラム記述子が、簡単で粗いスカラー量子化による圧縮に適している点であり、それによって、各記述子要素が、1または2ビットと同じくらい低い、たった数ビットに、個々に量子化され、実際に、特定の勾配ビン間の相関、または、特定のビンの他の特定のビンに対する優劣性の粗い測定のみを提供する。しかしながら、前述の粗い測定は、画像コレスポンデンスを確立、または、検証する点で成功であり、かつ、同一のスカラー量子化条件下、および、ベクトル量子化およびタイプ・コーディング等のより複雑な圧縮スキーム下の記述子にも基づいた最初の勾配ヒストグラムより高い認識性能を達成している、ロバストで、判別可能で簡潔な画像記述子をやはり提供する。
例えば、スカラー量子化を使用すると、値vは、例えば以下の関係式(4)等でqレベルに個々に量子化されるけれども、これは限定的なことではなく、異なる数の量子化レベルは、合成記述子において望ましい簡潔さを達成するべく用いられてもよい。
q=2
OR
q=3
OR
q=4・・・関係式(4)
この量子化は、各値vの値を閾値のセットと比較することによって実行されてもよい。量子化閾値を決定するには多数の選択肢がある。量子化閾値は、例えば、各vのダイナミックレンジの均一な分割によって決定されてもよい。代替的に、量子化閾値は、例えば一様分布のように、qの量子化レベル間の値の特定の長期分布を達成するべくセットされてもよい。さらに、量子化閾値は、全ての部分領域SRにおいて算出される全ての値vに対して同一であってもよく、全ての部分領域SRに亘って、同一の指標jに対して同一であってもよく、全てのjに亘って、各部分領域SRに対して同一であってもよく、各部分領域SRおよび各指標jに対して異なっていてもよい。
従って、上述の勾配ヒストグラム変換処理は、同様の、または、改良された性能特性を達成しながら、全体的な計算の複雑性および必要メモリにおける関連する低減と共に、これより前に見られた複雑な圧縮スキーム等に対する必要性を排除している。
結果として続く第3の利点は、変換された勾配ヒストグラム記述子の次元が、圧縮前だけでなく、圧縮後にも非常にスケーラブルである点と、必要ならば、その個々に量子化された要素の1または複数を簡単に除去することで、アプリケーションまたは伝送チャネル特性によって、変換された勾配ヒストグラム記述子の次元が簡単にさらに減少されてもよい点であり、それは、ベクトル量子化またはタイプ・コーディング等のより複雑な圧縮技術で簡単に達成されるものではない。
それと同時に、上述のような各部分領域の勾配ヒストグラムの個々の変換および圧縮は、非常に有利である。なぜならば、もし必要とされるならば、1または複数の個々に変換され量子化された部分領域ヒストグラムを簡単に除去することによって、全体の変換された画像記述子のスケーラビリティを助けることになるからである。
しかしながら、上述の変換された勾配ヒストグラム記述子が、これより前に議論されたより複雑な圧縮スキーム、または、勾配に基づく記述子の最初のヒストグラムに適用可能な如何なる圧縮スキームと矛盾していないこと、および、もし望まれるならば、そのような圧縮スキームと併せて用いられてもなおよいこと、も強調されるべきである。
これより前に見られるように、画像記述子Hに基づく勾配の単一のヒストグラムは、画像キーポイントKP周りのローカル領域Rの部分領域SRにそれぞれ対応する勾配hの複数のヒストグラムを構成してもよく、そのような画像記述子Hは、上述のような各部分領域SRの勾配ヒストグラムの個々の変換によって変換されてもよい。より具体的には、各部分領域勾配ヒストグラムは、他の部分領域勾配ヒストグラムと同一の方式で変換されてもよいが、代替的に、少なくとも1つの部分領域勾配ヒストグラムが、他の部分領域勾配ヒストグラムと異なって変換されてもよい。
例えば、いくつかの部分領域ヒストグラムを他の部分領域ヒストグラムと異なって変換することの文脈において、勾配のいずれのヒストグラムの変換関数も、勾配のその空間的に隣接するヒストグラムの少なくともいくつかに於いて使用される変換関数と少なくとも部分的に異ならなくてはいけないという要件が、画像記述子の変換処理に課されてもよい。これは、図5で示されている。
より具体的には、図5におけるAおよびBは、互いに完全に、または、部分的に異なる変換関数のセットを参照している。例えば、数式(1)の関数をAと、数式(2)の関数をBとみなすことで、関数の完全に異なるセットがもたらされる一方で、数式(3)の関数および数式(1)の偶数要素をAと、数式(3)の関数および数式(1)の奇数要素をBとみなすことで、関数の部分的に異なるセットがもたらされる。図5においては、A(または代わりにB)によって変換された各部分領域ヒストグラムに対し、B(または代わりにA)によってその4つの水平および垂直に隣接するヒストグラムが変換されるように、変換関数AおよびBが部分領域ヒストグラムに適用される。変換関数の複数のセットのこの利用が少しばかり実施の複雑性を増大させながらも、全ての部分領域ヒストグラムにおいて同一の変換関数を利用する上で大きな利点を有している。SIFT記述子を用いれば、勾配に基づく画像記述子の他のヒストグラムと同様に、特に水平または垂直に接続された隣接する勾配ヒストグラムが、変換された勾配ヒストグラム記述子においても維持される重要な相関性を示していることは、当業者に知られている。その後のいずれかの粗い量子化と組合せて、そのような隣接するヒストグラム内で同一の変換を適用することにより、隣接するヒストグラムが同一の要素を持つ可能性が高まるという結果がもたらされる。この問題は、上述のような異なる変換関数を利用することで緩和され、一様化と、結果として続く、記述子の弁別力を増大させる。
いくつかの部分領域ヒストグラムを他の部分領域ヒストグラムと異なって変換することの文脈における他の例として、いくつかの部分領域の変換された勾配ヒストグラム記述子は、キーポイントの中央または周りの領域に対する部分領域の位置等の特定の決定要素に従って、他の部分領域の変換された勾配ヒストグラム記述子と異なる要素の数を構成してもよい。
図6はこれの1つの例を示しており、12個の境界部分領域SRが上述のように要素kへと変換される一方で、画像キーポイントKPを直接囲んでいる中央の4つの部分領域SRは上述のように要素kへと変換され、以下の関係式(5)を満たす。
>k・・・関係式(5)
それゆえ、これは、全体的な変換された画像記述子に対する中央部分領域SRcの情報内容を増大させており、前述の部分領域SRcは、当業者によって、画像コレスポンデンスを確立または検証する点でより重要であると一般に考えられている。
いくつかの部分領域ヒストグラムを他の部分領域ヒストグラムと異なって変換することの文脈における他の例として、いくつかの部分領域の変換された勾配ヒストグラム記述子は、キーポイント周りの領域の中央に対する部分領域の位置等の特定決定要素に従って、他の部分領域の変換された勾配ヒストグラム記述子と異なって圧縮されてもよい。
図7は、これの1つの例を示しており、12個の境界部分領域SRが上述のように変換され、qレベルへと量子化される一方で、画像キーポイントKPを直接囲む中央の4つの部分領域SRは上述のように変換され、qレベルへと量子化され、以下の関係式(6)を満たす。
>q・・・関係式(6)
従って、これは全体的な変換された画像記述子の中で中央部分領域SRの表現精度を向上させ、前述の部分領域SRは、当業者によって、画像コレスポンデンスを確立または検証する点でより重要であると一般に考えられている。
さらに、特異的要素の変換関数またはそれらの元となる統計等の特定の決定要素に従って、異なる圧縮特性が、変換された勾配ヒストグラム記述子の特異的要素だけに適用されてもよい。
全体的にみれば、異なる部分領域に対する適切な変換関数の選択、および、異なる部分領域および/または部分領域要素に対する適切な圧縮パラメータの選択等の要因がかなり重要である。なぜならば、その複雑性およびストレージ/伝送コストに対して、それらが記述子の情報内容および弁別力を制御するからである。
本発明の実施形態は、図8から図13を用いて今説明される。
[第1の実施形態]
本発明の第1の実施形態は、図8において示されている。例えば図4aおよび図4bに従って、勾配ヒストグラムhを構成する勾配ヒストグラムに基づいた記述子Hが、それらの近接度に従う選択された勾配ヒストグラムビン間の相関に対応する値を算出することでその次元を低減するべく、各勾配ヒストグラムhを変換することで処理される。
より具体的には、図8において、以下の関係式(7)を満たす各ヒストグラムhが、順に処理される。
p=0...N−1 (N=16)・・・関係式(7)
ステップS100において、pは0にセットされる。
そして、ステップS110において、n個(n=8)のビンのヒストグラムhが、隣接するビン間の相関を捉える単一の角度分離、すなわち数式(9)の関数を有するビン間の相関を捉えるべく選択された関数のセットによって、以下の関係式(8)と共に、k個(k=4)の要素の記述子vに変換される。
k<n・・・関係式(8)
=h−h
=h−h
=h−h
=h−h・・・数式(9)
代替的に、数式(9)の関数は、90°で分離されたビン間の相関を捉える関数、または、異なる要素の数kを備える記述子vをもたらす関数等の、他の単一の角度分離を有するビン間の相関を捉える、異なる関数によって置換されてもよいが、なお以下の関係式(9)を満たす。
k<n・・・関係式(9)
そして、ステップS120においてpの値が考察され、もしHにおける最後のヒストグラムの指標と等しいならば処理は終わり、そうでないならばステップS130で1増やされ、処理はステップS110に戻る。
従って、合成記述子Vは、大幅にHより少ない要素を使用するHの顕著な形状特性および判別可能な情報を符号化し、画像コレスポンデンスを確立または検証する点、および、高い認識性能を達成する点で依然としてとても成功である、より小さな記述子という結果をもたらす。
例えば、非負の値、および/または、特定のダイナミックレンジ、例えば8ビットに変換するべく、任意で、Vの個々の要素の値が適切になるよう調整され、かつ、変えられてもよい。
2つの記述子が異なる画像における同一のキーポイントから抽出されたかどうかを決定する、および、画像間の画像コレスポンデンスを確立するための、合成記述子のその次の処理は、本発明の範囲を超えており、最初の勾配ヒストグラム記述子に関しては、例えばInternational Journal of Computer Vision, 60, 2 (2004), 91‐110ページ内のDavid G. Loweによる「Distinctive image features from scale−invariant keypoints」で説明されているものと類似した方式で、または、任意の他の適切な方式で進んでもよい。
[第2の実施形態]
本発明の第2の実施形態は、図9において示されている。例えば図4aおよび図4bによる、勾配ヒストグラムhを構成している勾配ヒストグラムに基づく記述子Hは、それらの近接に関する複数の基準に従って選択された勾配ヒストグラムビン間の相関に対応する値を算出することによって次元を低減するべく、各勾配ヒストグラムhを変換することで処理される。
より具体的には、図9において、以下の関係式(10)を満たす各ヒストグラムhが、順に処理される。
p=0...N−1 (N=16)・・・関係式(10)
ステップS200において、pは0にセットされる。
そして、ステップS210において、n個のビン(n=8)のヒストグラムhが、複数の角度分離を有するビン間の相関を捉えるべく選択された関数のセット、すなわち、隣接するビン間の相関、および、180°で分離されたビン間の相関を捉える、数式(10)の関数によって、以下の関係式(11)を満たす、k個の要素(k=6)の記述子vに変換される。
k<n・・・関係式(11)
=h−h
=h−h
=h−h
=h−h
=h−h
=h−h・・・数式(10)
代替的に、数式(10)の関数は、90°で分離されたビン間の相関、および、180°で分離されたビン間の相関を捉える関数、または、異なる要素の数kを備える記述子vという結果をもたらす関数等の、他の複数の角度分離を有するビン間の相関を捉える異なる関数によって置換されてもよいが、なお以下の関係式(12)を満たす。
k<n・・・関係式(12)
そして、ステップS220においてpの値が考察され、もしHにおける最後のヒストグラムの指標と等しければ処理は終了し、そうでないならばステップS230で1増やされ、処理はステップS210に戻る。
従って、合成記述子Vは、Hより少ない要素を使用するHの、異なるタイプの顕著な形状特性および判別可能な情報を符号化し、画像コレスポンデンスを確立または検証し、かつ、高い認識性能を達成する点で依然としてとても成功である、より小さな記述子という結果をもたらす。
そして、合成記述子のその次の処理は、最初の実施形態に類似した方式で進んでもよい。
[第3の実施形態]
本発明の第3の実施形態は、図10において示されている。例えば図4aおよび図4bによる、勾配ヒストグラムhを構成している勾配ヒストグラムに基づく記述子Hは、それらの近接に関する複数の基準に従って選択された勾配ヒストグラムビン間の相関に対応する値を算出することによって、各勾配ヒストグラムhを変換することで処理される。
より具体的には、図10において、以下の関係式(13)を満たす各ヒストグラムhが順に処理される。
p=0...N−1 (N=16)・・・関係式(13)
ステップS300において、pは0にセットされる。
そして、ステップS310において、n個のビン(n=8)のヒストグラムhが、複数の角度分離を有するビン間の相関を捉えるべく選択された関数のセット、すなわち、数式(11)の関数によって、k個の要素(k=8)の記述子vに変換される。
=h−h
=h−h
=h−h
=h−h
=h−h
=h−h
=(h+h)−(h+h
=(h+h+h+h)−(h+h+h+h)・・・数式(11)
代替的に、数式(11)の関数は、他の複数の角度分離を有するビン間の相関を捉える異なる関数、または、異なる要素の数kを備える記述子vという結果をもたらす関数によって置換されてもよく、以下の関係式(14)を満たす。
k=n
OR
k<n
OR
k>n・・・関係式(14)
そして、ステップS320において、k個の要素の記述子vの各要素がqレベルへと個々に量子化され、量子化された記述子
Figure 2017224337
を与える(以下、同様に上に波線「〜」がある「v」および「V」はそれぞれ「v〜」および「V〜」と表す)。この実施形態において、我々はq=3をセットするけれども、これは限定的ではなく、合成記述子における望ましい簡潔さを達成するべく、例えばq=2、q=4等の、量子化レベルの異なる数が用いられてもよい。この量子化は、各要素vpjの値を閾値のセットと比較することによって実行されてもよい。量子化閾値を決定する場合に多数の選択肢がある。量子化閾値は、例えばvpjのダイナミックレンジの均一な分割によって決定されてもよい。代替的に、例えば一様分布といった、qの量子化レベル間の値の特定の長期分布を達成するべく、量子化閾値がセットされてもよい。さらに、量子化閾値は、全ての要素vpjに対して同一であってもよく、全てのpに亘る同一の指標jに対して同一であってもよく、全てのjに亘る同一の指標pに対して同一であってもよく、各vpjに対して異なっていてもよい。
そして、ステップS330において、pの値は考察され、もしHにおける最後のヒストグラムの指標と等しければ処理は終了し、そうでないならばステップS340において1増やされ、処理はステップS310に戻る。
明らかであるが、代替的な実施は、例えば最初に全体の記述子Vを算出し、そして量子化された記述子V〜を生成するべく量子化を進めることによって、本発明の範囲から逸脱せずに、図10と比較して適切に順序または特定の動作を変化させてもよい。
従って、合成記述子V〜は、Hの異なるタイプの顕著な形状特性および判別可能な情報を符号化する。粗いスカラー量子化は、特定の勾配ビン間の相関、または、他の特定のビンに対する特定のビンの優劣性の粗い測定のみを提供する、ほんのわずかなレベルへと各記述子要素を個々に量子化する。そして、画像コレスポンデンスを確立または検証する点で成功であり、同一のスカラー量子化条件下で、かつ、より複雑な圧縮スキーム下でも、最初の勾配ヒストグラムに基づく記述子より高い認識性能を達成する、ロバストで、判別可能で簡潔な画像記述子という結果をもたらす。
そして、合成記述子のその次の処理は、前の実施形態と類似した方式で進んでもよい。
[第4の実施形態]
本発明の第4の実施形態は、図11において示されている。例えば図4aおよび図4bに従って、勾配ヒストグラムhを構成している勾配ヒストグラムに基づく記述子Hは、隣接するhである勾配ヒストグラムのいくつかの処理において用いられる変換関数のセットとの違いを示す変換関数のセットで、各勾配ヒストグラムhを変換することによって処理される。
より具体的には、図11において、以下の関係式(15)を満たす各ヒストグラムhが順に処理される。
p=0...N−1 (N=16)・・・関係式(15)
ステップS400において、pは0にセットされる。
そして、ステップS410において、hの処理に適切な関数を選択すべく、pの値が考察される。この選択は、例えば図5において示されたラインに沿って進む。AおよびBと示された2つの変換のセットが用いられ、p=0,2,5,7,8,10,13,15の場合にAが選択され、および、p=1,3,4,6,9,11,12,14の場合にBが選択されるB。しかしながら、これは限定的ではなく、勾配ヒストグラムの処理において用いられる変換関数のセットが、隣接する前述のヒストグラムである勾配ヒストグラムのいくつかの処理において用いられる変換関数のセットとの違いを示すべく、異なる空間配置と共に変換関数の異なる数のセットが使用されてもよい。
そして、ステップS420において、適切に選択された変換関数のセットを使用して、n個のビン(n=8)のヒストグラムhがk個の要素(k=8)の記述子vに変換される。変換関数AおよびBのセットは、夫々以下に示される数式(12)および(13)により規定されてもよい。各関数のセットは複数の角度分離を有するビン間の相関を捉えるべく選択され、AとBとの間で共通する関数はない。
=h−h
=h−h
=h−h
=h−h
=h−h
=h−h
=(h+h)−(h+h
=(h+h+h+h)−(h+h+h+h)・・・数式(12)
=h−h
=h−h
=h−h
=h−h
=h−h
=h−h
=(h+h)−(h+h
=(h+h+h+h)−(h+h+h+h)・・・数式(13)
代替的に、数式(12)および(13)において示されるような関数のどちらか、または、両方のセットが、他の複数の角度分離を有するビン間の相関を捉える異なる関数、または、異なる要素の数kを備える記述子という結果をもたらす関数によって置換されてもよく、以下の関係式(16)を満たす。
k=n
OR
k<n
OR
k>n・・・関係式(16)
さらに、関数AおよびBのセットは、共通する関数を1つも含まなくてもよく、また、いくつかの共通する関数を含んでもよい。さらに、相関に対応する値がその間で算出されるビンは、関数AおよびBのセットのどちらか、または、両方に対して、異なる抽出条件に従い、選択されてもよい。そのような基準は、これより前に見られるように、合成値の弁別力を向上すべく、相関に対応する値がその間で算出されるであろうビンを選択すること、または、相関に対応する値がランダム選択処理を通してその間で算出されるであろうビンを選択することを含んでもよい。さらにまた、前述の基準のいずれかの混合が、関数AおよびBのセットのどちらか、または、両方の選択において用いられてもよい。
そして、ステップS430において、k個の要素の記述子vの各要素がqレベルへと個々に量子化され、第3の実施形態に関して類似した方式で量子化された記述子v〜を与える。
そして、ステップS440において、pの値が考察され、もしHにおける最後のヒストグラムの指標と等しければ処理は終了し、そうでないならばステップS450において1増やし、処理はステップS410に戻る。
明らかであるが、代替的な実施は、例えば最初に全体の記述子Vを算出し、そして量子化された記述子V〜を生成するべく量子化を進めることによって、本発明の範囲から逸脱せずに、図11と比較して適切に順序または特定の動作を変化させてもよい。
従って、合成記述子V〜はHの判別可能な情報を符号化する。この実施形態のように複数のセットの変換関数を利用することは、全ての部分領域ヒストグラムにおいて同一変換関数を利用することに対して大きな利点を有する。勾配のヒストグラムに基づく画像記述子を用いれば、隣接する勾配ヒストグラムは変換された勾配ヒストグラム記述子において維持され、粗いスカラー量子化との組合せで、隣接する変換されたヒストグラムが同一の要素を持つ可能性を高めるという結果をもたらす、有意な相関を示すことは、当業者に知られている。この問題は、一様化を向上し、かつ、その結果として記述子の弁別力を向上するこの実施形態により、異なる変換関数を利用することで緩和される。
そして、合成記述子のその次の処理は、前の実施形態と類似した方式で進めてもよい。
[第5の実施形態]
本発明の第5の実施形態は、図12において示されている。例えば図4aおよび図4bに従って、勾配ヒストグラムhを構成している勾配ヒストグラムに基づく記述子Hは、勾配ヒストグラムビン間の相関に対応する値を算出することで各勾配ヒストグラムhを変換することによって、かつ、いくつかの部分領域の変換された勾配ヒストグラム記述子が他の部分領域の変換された勾配ヒストグラム記述子と異なる要素の数を含むべく、処理される。
より具体的には、図12において、以下の関係式(17)を満たす各ヒストグラムhが順に処理される。
p=0...N−1 (N=16)・・・関係式(17)
ステップS500において、pは0にセットされる。
そして、ステップS510において、hの処理のための適切な次元を選択するべく、pの値が考察される。この選択は、例えば図6において示されたラインに沿って進む。p=5,6,9,10を満たす中央部分領域ヒストグラムは、合成記述子が各(k=3)の要素kを有するべく変換され、p=0,1,2,3,4,7,8,11,12,13,14,15を満たす境界部分領域ヒストグラムは、合成記述子が各(k=2)の要素kを有するべく変換される。しかしながら、これは限定的ではなく、異なる数の可能な次元が、異なる空間配置と共に、使用されてもよい。
そして、ステップS520において、n個のビン(n=8)のヒストグラムhが、以下に示される数式(14)の関数に従って、勾配ヒストグラムビン間の相関に対応する値を算出することによって、記述子vに変換される。
=h−h (中央および境界部分領域の両方に対して用いられる)
=h−h (中央および境界部分領域の両方に対して用いられる)
=(h+h)−(h+h) (中央境界部分領域に対してのみ用いられる)・・・数式(14)
明らかであるが、この実施形態は、いずれかの前の実施形態と組み合わされてもよく、数式(14)の変換関数は、これより前に見られるような選択条件のいずれかの混合に従って選択された変換関数によって置換されてもよく、および/または、異なるセットの変換関数が異なる部分領域ヒストグラムに対して用いられてもよい。
そして、ステップS530において、第3および第4の実施形態に対するものと類似した方式で量子化された記述子v〜を与えて、k個の要素の記述子vの各要素がqレベルへと個々に量子化される。
そして、ステップS540において、pの値が考察され、もしHにおける最後のヒストグラムの指標と等しければ処理は終了し、そうでないならばステップS550において1増やし、処理はステップS510に戻る。
明らかであるが、代替的な実施は、例えば最初に全体の記述子Vを算出し、そして量子化された記述子V〜を生成するべく量子化を進めることによって、本発明の範囲から逸脱せずに、図12と比較して適切に順序または特定の動作を変化させてもよい。
従って、合成記述子V〜は、境界部分領域ヒストグラムに対する中央部分領域ヒストグラム等の画像コレスポンデンスを確立または検証する点でより重要であり得る、Hのこれら部分をより顕著にする方式で、これら部分に対して増大した次元の表現を可能にすることによって、Hの判別可能な情報を符号化する。
[第6の実施形態]
本発明の第6の実施形態は、図13において示されている。例えば図4aおよび図4bに従って、勾配ヒストグラムhを構成している勾配ヒストグラムに基づく記述子Hは、勾配ヒストグラムビン間の相関に対応する値を算出することで各勾配ヒストグラムhを変換することによって、かつ、いくつかの変換された勾配ヒストグラム記述子要素が他の変換された勾配ヒストグラム記述子要素と異なるレベルの数へと量子化されるべく、処理される。
より具体的には、図13において、以下の関係式(18)を満たすヒストグラムhが順に処理される。
p=0...N−1 (N=16)・・・関係式(18)
ステップS600において、pは0にセットされる。
そして、ステップS610において、n個のビン(n=8)のヒストグラムhが、以下の数式(15)で示されるような、ビン間の相関を捉えるべく選択された関数のセットによって、k個の要素(k=8)の記述子vに変換される。
=h−h
=h−h
=h−h
=h−h
=h−h
=h−h
=(h+h)−(h+h
=(h+h+h+h)−(h+h+h+h)・・・数式(15)
明らかであるが、この実施形態は、いずれの前の実施形態と組み合わされてもよく、数式(15)の変換関数は、これより前に見られるような選択条件のいずれかの混合に従って選択された変換関数によって置換されてもよく、および/または、変換関数の異なるセットが異なる部分領域ヒストグラムに対して用いられてもよく、および/または、異なる変換された記述子の次元が異なる部分領域ヒストグラムに対して用いられてもよい。
そして、ステップS620において、各要素vpjの量子化に対する量子化レベルの適切な数を選択するべく、pの値が考察される。この選択は、例えば図7において示されたラインに沿って進む。p=5,6,9,10を満たす中央部分領域記述子はqレベル(q=4)へと量子化され、p=0,1,2,3,4,7,8,11,12,13,14,15を満たす境界部分領域記述子はqレベル(q=2)へと量子化される。しかしながら、これは限定的ではなく、異なる数の可能な量子化レベルが、異なる空間配置と共に、使用されてもよい。さらに、これは図13において示されていないけれども、各要素vpjに対する量子化レベルの数は、pの値、すなわち要素の部分領域位置に代えて、または、それに加えて、jの値に従い、すなわち、記述子要素の特殊なタイプに従い、決定されてもよい。
そして、ステップS630において、第3、第4および第5の実施形態に関して類似した方式で量子化された記述子v〜を与え、k個の要素の記述子vの要素が量子化レベルの適切な数へと個々に量子化される。
そして、ステップS640において、pの値が考察され、もしHにおける最後のヒストグラムの指標と等しければ処理は終了し、そうでないならばステップS650において1増やし、処理はステップS610に戻る。
明らかであるが、代替的な実施は例えば最初に全体の記述子Vを算出し、そして量子化された記述子V〜を生成するべく量子化に進むことによって、本発明の範囲から逸脱せずに、図13と比較して適切に順序または特定の動作を変化させてもよい。
従って、合成記述子V〜は、画像コレスポンデンスを確立または検証する点でより重要であろう、Hのこれら部分またはVのこれらの要素に対して、より高度な表現精度を与える方式で、Hの判別可能な情報を符号化する。
図4aおよび図4bにおいて示されているようなSIFT画像記述子から、ロバストで、判別可能で、スケーラブルかつ簡潔な画像記述子を算出すべく、本発明の態様および実施形態が詳細に示されているけれども、本発明は勾配のヒストグラムに基づく他の画像記述子に適用可能である。例えば、図14は、隣接する勾配ヒストグラムの処理において用いられる変換関数のセットとの違いを示す変換関数のセットと共に勾配ヒストグラムの夫々を変換することによって、本発明の第4の実施形態と類似した方式で、前述の変換関数のセットに対応する図14におけるシンボルA、BおよびCと共に、log−polar空間細分を使用する勾配ヒストグラムに基づいた記述子の処理を示している。他の例として、勾配のXおよびY成分の2次元空間の、図15において示されるようなビンへの細分に基づいた勾配のヒストグラムに対して、本発明の前の実施形態と類似した方式における適切な1つの変換関数のセットは以下の数式(16)である。
=h−h
=h−h
=h−h
=h−h
=h−h
=h−h
=h−h
=h−h
=h−(h+h+h+h+h+h+h+h)・・・数式(16)
さらに、本発明は、適切に処理された勾配のヒストグラムに基づく画像記述子に適用可能でもある。そのような適切な処理は、例えば、本発明に従う処理より前に、部分領域勾配ヒストグラムの組合せを必要としてもよい。
説明的な目的のためだけに図4aを参照し、p=0...15を備える部分領域ヒストグラムhを構成する勾配ヒストグラム記述子Hに対して、p=0、1、4、5を備えるヒストグラムが、それらのビンの値を平均化することによって単一の部分領域ヒストグラムへと組み合わされてもよく、類似した組合せがp=2、3、6、7およびp=8、9、12、13およびp=10、11、14、15を備えるヒストグラムに対して実行されてもよい。これは、本発明の前の任意の実施形態に従ってその後に処理され得る減らされた次元を備える、勾配ヒストグラム記述子という結果をもたらす。勾配ヒストグラムはまた、和等の、それらのビンの値の代替的な関数によって組み合わされてもよい。
代替的に、または加えて、そのような適切な処理は、例えば、部分領域勾配ヒストグラム内でのビンの合併を必要としてもよい。説明的な目的のためだけに図4bを参照し、勾配ヒストグラムhに対して、隣接するビンが、平均化、加算、メディアンまたは任意の適切な関数によって、単一のビンへと合併されてもよい。これは、本発明の前の任意の実施形態に従ってその後に処理され得る減らされた次元を備える、勾配ヒストグラム記述子という結果をもたらす。
例としてのみ、図16は、本発明に従う方法を実行する概念的な処理装置を示す。より具体的には、処理装置1100は、画像またはビデオデータ等の画像データを構成する入力、勾配のヒストグラムに基づいて予め算出された記述子、本発明の方法に従って予め算出された簡潔な記述子、プログラミング命令、またはユーザ入力を、ユーザ入力装置、メディアリーダ、または、送信されたシグナルのレシーバの形をとってもよい入力装置1000から受け取る。処理装置1100は、他の処理ブロック、揮発性メモリ1120、不揮発性メモリ1130、任意で、勾配のヒストグラムに基づく記述子を生成すべく構成された記述子抽出ブロック1140、本発明に従う方法を実行すべく構成された簡潔な記述子抽出ブロック1150、および、任意で、例えば画像コレスポンデンスを確立または検証するべく前述の簡潔な記述子を処理するように構成された簡潔な記述子処理ブロック1160の動作を制御する、中央演算処理装置1110のメイン処理ブロックを構成する。処理装置1100は、画像ディスプレイユニット、メディアライタ、または、画像またはビデオデータ等の注釈付画像データ、確立または検証された画像コレスポンデンス等の処理情報、または、本発明の方法に従って算出された簡潔な記述子を含み得る出力を提供する、シグナルの伝送器の形をとり得る出力装置1900に接続される。図16において示されている処理ブロックおよび構成が単に概念的であって、本発明に従う方法を実施する全ての装置に厳密に対応しなくてもよいことが理解されるべきである。
本明細書で例として説明された、勾配ヒストグラムに基づく画像記述子を変換する方法および関連する画像処理装置は、発明的思想の新規な真意から逸脱することなく多くの可能なバリエーションが施されてもよい。本発明の実際的な実施において、説明された詳細が異なる形状を有してもよく、または、他の技術的に等価な要素と置換されてもよいことも明らかである。
それゆえ、本発明が、勾配ヒストグラムに基づく画像記述子を変換する方法および関連する画像処理装置に制限されないが、以下の特許請求の範囲において明確に特定されるような発明的な思想から逸脱することなく、多くの変更、改良または等価な部分および要素の置換が施されてもよいことが容易に理解され得る。

Claims (19)

  1. 複数のヒストグラムビン(h)を含む勾配ヒストグラム(h)に基づく画像記述子を、複数の値(v)からなるセットを構成する変換された勾配ヒストグラム記述子(v)に変換する方法であって、
    前記複数の値(v)からなるセットを算出するべく少なくとも1つの基準を用いる段階を含み、前記少なくとも1つの基準は、前記複数のヒストグラムビン(h)の配置基準を含む方法。
  2. 前記配置基準は、単一の角度分離を示す複数の勾配ヒストグラムビン(h)間の複数の値(v)を算出することを規定する請求項1に記載の方法。
  3. 前記配置基準は、複数の角度分離を示す複数の勾配ヒストグラムビン(h)間の複数の値(v)を算出することを規定する請求項1に記載の方法。
  4. 前記複数の値(v)の少なくとも1つは、隣接する複数のヒストグラムビン(h)間で算出される請求項2または3に記載の方法。
  5. 前記複数の値(v)は複数のレベル(q)へと個々に量子化され、前記複数のレベルは、複数の値(v)の夫々に対して同一であるか、または、少なくとも2つの値(v)に対して異なる請求項1から4のいずれか一項に記載の方法。
  6. 前記複数の値(v)は、複数のヒストグラムビン(h)間の和または差、複数のヒストグラムビン(h)の和の間の差、複数のヒストグラムビン(h)間の線形または非線形演算処理、複数のヒストグラムビン(h)間の比、複数のヒストグラムビン(h)の和の間の比、複数のヒストグラムビン(h)の生成物間の比、複数のヒストグラムビン(h)の生成物間の差のうちの、1または複数の動作に従って算出される請求項1から5のいずれか一項に記載の方法。
  7. 前記勾配ヒストグラム(h)は複数の要素の第1の数(n)の次元を有し、前記複数の値(v)からなるセットは複数の要素の第2の数(k)の次元を有し、複数の要素の前記第1の数(n)は複数の要素の前記第2の数(k)と等しい請求項1から6のいずれか一項に記載の方法。
  8. 前記勾配ヒストグラム(h)は複数の要素の第1の数(n)の次元を有し、前記複数の値(v)からなるセットは複数の要素の第2の数(k)の次元を有し、複数の要素の前記第1の数(n)は複数の要素の前記第2の数(k)より大きい請求項1から6のいずれか一項に記載の方法。
  9. 前記勾配ヒストグラム(h)は複数の要素の第1の数(n)の次元を有し、前記複数の値(v)からなるセットは複数の要素の第2の数(k)の次元を有し、複数の要素の前記第1の数(n)は複数の要素の前記第2の数(k)より小さい請求項1から6のいずれか一項に記載の方法。
  10. 複数の勾配ヒストグラム(h)に基づく画像記述子(H)を、変換された画像記述子(V)に変換する方法であって、
    前記複数の勾配ヒストグラム(h)の各々は、画像キーポイント(KP)周りのローカル領域(R)の部分領域(SR)に関係し、前記複数の勾配ヒストグラム(h)の夫々は、請求項1から9のいずれか一項に従って、変換された勾配ヒストグラム記述子(V)へと変換される複数のヒストグラムビン(h)を含む、方法。
  11. 複数の部分領域勾配ヒストグラムのうちの少なくとも1つの部分領域勾配ヒストグラムは、他の複数の部分領域勾配ヒストグラムと異なって変換される請求項10に記載の方法。
  12. 複数の部分領域勾配ヒストグラムは、単一の部分領域ヒストグラムにおいて、特に、それらの複数のビンの値を合計する、または、平均化することによって、組み合わされる請求項11に記載の方法。
  13. 変換された前記少なくとも1つの部分領域勾配ヒストグラムを得るべく使用される動作は、複数の隣接する部分領域の少なくとも1つの部分領域の、複数の変換された部分領域勾配ヒストグラムを得るべく使用される動作と、少なくとも部分的に異なる請求項11に記載の方法。
  14. 変換された前記少なくとも1つの部分領域勾配ヒストグラムは、複数の隣接する部分領域の少なくとも1つの部分領域の、複数の前記変換された部分領域勾配ヒストグラムとは異なる数の複数の要素を含む請求項11に記載の方法。
  15. 前記画像キーポイント(KP)を直接囲む複数の部分領域(SR)に対する複数の変換された部分領域勾配ヒストグラムは、複数の残っている部分領域(SR)の複数の変換された勾配ヒストグラムの複数の要素(k)より多くの複数の要素(k)を含む請求項14に記載の方法。
  16. 前記画像キーポイント(KP)を直接囲む複数の部分領域(SR)の前記複数の変換された勾配ヒストグラムの複数の要素は、複数の残っている部分領域(SR)の複数の変換された勾配ヒストグラムが量子化されるレベル(q)の数よりも多い数のレベル(q)で量子化される請求項15に記載の方法。
  17. 前記量子化は、前記複数の部分領域勾配ヒストグラムの複数の特異的要素に対してのみ適用される請求項16に記載の方法。
  18. 前記変換された画像記述子(V)は、量子化されて変換された画像記述子(V〜)を生成するべく、量子化処理が施される請求項10に記載の方法。
  19. 請求項1から18のいずれか一項に記載の方法を実行する手段を備える画像処理装置。
JP2017154781A 2012-07-09 2017-08-09 勾配ヒストグラムに基づいて画像記述子を変換する方法および関連する画像処理装置 Active JP6429134B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
ITTO2012A000602 2012-07-09
IT000602A ITTO20120602A1 (it) 2012-07-09 2012-07-09 Method for transforming an image descriptor based on a gradient histogram and relative image processing apparatus.

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP2015520897A Division JP6292454B2 (ja) 2012-07-09 2013-07-01 勾配ヒストグラムに基づいて画像記述子を変換する方法および関連する画像処理装置

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2018196276A Division JP6714669B2 (ja) 2012-07-09 2018-10-18 勾配ヒストグラムに基づいて画像記述子を変換する方法および関連する画像処理装置

Publications (2)

Publication Number Publication Date
JP2017224337A true JP2017224337A (ja) 2017-12-21
JP6429134B2 JP6429134B2 (ja) 2018-11-28

Family

ID=46833121

Family Applications (3)

Application Number Title Priority Date Filing Date
JP2015520897A Active JP6292454B2 (ja) 2012-07-09 2013-07-01 勾配ヒストグラムに基づいて画像記述子を変換する方法および関連する画像処理装置
JP2017154781A Active JP6429134B2 (ja) 2012-07-09 2017-08-09 勾配ヒストグラムに基づいて画像記述子を変換する方法および関連する画像処理装置
JP2018196276A Expired - Fee Related JP6714669B2 (ja) 2012-07-09 2018-10-18 勾配ヒストグラムに基づいて画像記述子を変換する方法および関連する画像処理装置

Family Applications Before (1)

Application Number Title Priority Date Filing Date
JP2015520897A Active JP6292454B2 (ja) 2012-07-09 2013-07-01 勾配ヒストグラムに基づいて画像記述子を変換する方法および関連する画像処理装置

Family Applications After (1)

Application Number Title Priority Date Filing Date
JP2018196276A Expired - Fee Related JP6714669B2 (ja) 2012-07-09 2018-10-18 勾配ヒストグラムに基づいて画像記述子を変換する方法および関連する画像処理装置

Country Status (13)

Country Link
US (2) US9454710B2 (ja)
EP (2) EP3702965B1 (ja)
JP (3) JP6292454B2 (ja)
KR (1) KR102193177B1 (ja)
CN (2) CN108520265B (ja)
AR (1) AR091674A1 (ja)
BR (1) BR112014031268A2 (ja)
ES (1) ES2811532T3 (ja)
HK (1) HK1259139A1 (ja)
IT (1) ITTO20120602A1 (ja)
RU (1) RU2661795C2 (ja)
TW (1) TWI532012B (ja)
WO (1) WO2014009198A1 (ja)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101710967B1 (ko) * 2013-01-16 2017-02-28 후아웨이 테크놀러지 컴퍼니 리미티드 이미지의 히스토그램 맵을 컨텍스트 기반 인코딩하기 위한 방법
ITTO20130629A1 (it) 2013-07-24 2015-01-25 Sisvel Technology Srl Method for encoding an image descriptor based on a gradient histogram and relative image processing apparatus
EP3671472B1 (en) 2017-09-25 2022-05-04 Huawei Cloud Computing Technologies Co., Ltd. Data access method and device
JP7087695B2 (ja) * 2018-06-07 2022-06-21 株式会社リコー 学習装置および学習方法
KR20230017549A (ko) * 2021-07-28 2023-02-06 주식회사 넥스트칩 특징점에 대한 기술자를 생성하기 위한 전자 장치 및 그 동작 방법
WO2023043885A1 (en) * 2021-09-15 2023-03-23 Beijing Dajia Internet Information Technology Co., Ltd., Sign prediction for block-based video coding
CN114331925B (zh) * 2022-03-15 2022-07-19 北京锐影医疗技术有限公司 一种多尺度相对梯度直方图均衡化方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000187731A (ja) * 1998-12-21 2000-07-04 Ricoh Co Ltd 画像特徴抽出方法およびその方法の各工程をコンピュータに実行させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体
US20090285459A1 (en) * 2008-05-15 2009-11-19 Gaurav Aggarwal Fingerprint representation using gradient histograms
US20110052050A1 (en) * 2009-09-02 2011-03-03 Kabushiki Kaisha Toyota Chuo Kenkyusho Image processing apparatus, image processing method, and computer readable storage medium

Family Cites Families (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61239376A (ja) 1985-04-16 1986-10-24 Fujitsu Ltd 画像間の角度差検出装置
GB2352076B (en) * 1999-07-15 2003-12-17 Mitsubishi Electric Inf Tech Method and apparatus for representing and searching for an object in an image
JP3550681B2 (ja) * 1999-12-10 2004-08-04 日本電気株式会社 画像検索装置及び方法、並びに類似画像検索プログラムを格納した記憶媒体
US6865295B2 (en) 2001-05-11 2005-03-08 Koninklijke Philips Electronics N.V. Palette-based histogram matching with recursive histogram vector generation
JP2003150965A (ja) * 2001-11-14 2003-05-23 Sanyo Electric Co Ltd 形状認識方法及びその装置
KR100959053B1 (ko) * 2003-01-13 2010-05-20 한국전자통신연구원 복수의 이미지 프레임을 갖는 비디오 시퀀스 검색을 위한비선형 양자화 및 유사도 매칭 방법
US7840081B2 (en) * 2004-09-23 2010-11-23 Mitsubishi Denki Kabushiki Kaisha Methods of representing and analysing images
US7715623B2 (en) * 2005-11-14 2010-05-11 Siemens Medical Solutions Usa, Inc. Diffusion distance for histogram comparison
CN100530239C (zh) * 2007-01-25 2009-08-19 复旦大学 基于特征匹配与跟踪的视频稳定方法
JP4629740B2 (ja) * 2008-01-16 2011-02-09 富士フイルム株式会社 撮像装置および方法並びにプログラム
KR101247891B1 (ko) 2008-04-28 2013-03-26 고리츠다이가쿠호징 오사카후리츠다이가쿠 물체 인식용 화상 데이터베이스의 작성 방법, 처리 장치 및 처리용 프로그램
US8532438B2 (en) 2008-05-09 2013-09-10 Empire Technology Development Llc Matching images with shape descriptors
JP4547639B2 (ja) * 2008-08-26 2010-09-22 ソニー株式会社 画像処理装置および方法、並びにプログラム
US8625921B1 (en) * 2008-09-26 2014-01-07 Google Inc. Method for image processing using local statistics convolution
WO2010055399A1 (en) * 2008-11-12 2010-05-20 Nokia Corporation Method and apparatus for representing and identifying feature descriptors utilizing a compressed histogram of gradients
CN101739397A (zh) * 2008-11-17 2010-06-16 新疆亚奥数码科技有限公司 基于mpeg-7的图像检索系统
CN101493891B (zh) * 2009-02-27 2011-08-31 天津大学 基于sift的具有镜面翻转不变性的特征提取和描述方法
US8773430B2 (en) 2009-12-09 2014-07-08 Thomson Licensing Method for distinguishing a 3D image from a 2D image and for identifying the presence of a 3D image format by feature correspondence determination
US8582889B2 (en) 2010-01-08 2013-11-12 Qualcomm Incorporated Scale space normalization technique for improved feature detection in uniform and non-uniform illumination changes
CN101794395B (zh) 2010-03-11 2012-04-25 合肥金诺数码科技股份有限公司 一种基于Sift算法的图像匹配定位方法
US9530073B2 (en) * 2010-04-20 2016-12-27 Qualcomm Incorporated Efficient descriptor extraction over multiple levels of an image scale space
US8625902B2 (en) * 2010-07-30 2014-01-07 Qualcomm Incorporated Object recognition using incremental feature extraction
US8428397B1 (en) * 2010-08-26 2013-04-23 Adobe Systems Incorporated Systems and methods for large scale, high-dimensional searches
US8538164B2 (en) * 2010-10-25 2013-09-17 Microsoft Corporation Image patch descriptors
JP5289412B2 (ja) * 2010-11-05 2013-09-11 株式会社デンソーアイティーラボラトリ 局所特徴量算出装置及び方法、並びに対応点探索装置及び方法
KR101127793B1 (ko) * 2011-01-04 2012-03-23 위드로봇 주식회사 이미지 인식 시스템 및 그 제공방법
CN102306383B (zh) * 2011-07-01 2013-05-29 北京交通大学 适用于宽基线图像密集匹配的描述符的构建方法
JP2013097583A (ja) * 2011-11-01 2013-05-20 Nec Corp 特徴量生成装置、方法及びプログラム
WO2013076364A1 (en) * 2011-11-21 2013-05-30 Nokia Corporation Method for image processing and an apparatus
EP2801055B1 (en) * 2012-01-02 2016-04-20 Telecom Italia S.p.A. Method and system for image analysis

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000187731A (ja) * 1998-12-21 2000-07-04 Ricoh Co Ltd 画像特徴抽出方法およびその方法の各工程をコンピュータに実行させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体
US20090285459A1 (en) * 2008-05-15 2009-11-19 Gaurav Aggarwal Fingerprint representation using gradient histograms
JP2011521355A (ja) * 2008-05-15 2011-07-21 インターナショナル・ビジネス・マシーンズ・コーポレーション 勾配ヒストグラムを用いた指紋表現
US20110052050A1 (en) * 2009-09-02 2011-03-03 Kabushiki Kaisha Toyota Chuo Kenkyusho Image processing apparatus, image processing method, and computer readable storage medium
JP2011053953A (ja) * 2009-09-02 2011-03-17 Toyota Central R&D Labs Inc 画像処理装置及びプログラム

Also Published As

Publication number Publication date
KR102193177B1 (ko) 2020-12-21
CN104428793B (zh) 2018-05-08
CN108520265B (zh) 2021-09-07
ES2811532T3 (es) 2021-03-12
BR112014031268A2 (pt) 2017-06-27
JP2019040612A (ja) 2019-03-14
WO2014009198A1 (en) 2014-01-16
AR091674A1 (es) 2015-02-18
RU2661795C2 (ru) 2018-07-19
TWI532012B (zh) 2016-05-01
KR20150031333A (ko) 2015-03-23
EP3702965B1 (en) 2021-08-04
US20170300775A1 (en) 2017-10-19
JP2015522194A (ja) 2015-08-03
HK1259139A1 (zh) 2019-11-29
EP2870568A1 (en) 2015-05-13
JP6714669B2 (ja) 2020-06-24
US10713523B2 (en) 2020-07-14
EP3702965A1 (en) 2020-09-02
US9454710B2 (en) 2016-09-27
EP2870568B1 (en) 2020-06-24
CN108520265A (zh) 2018-09-11
ITTO20120602A1 (it) 2014-01-10
JP6292454B2 (ja) 2018-03-14
CN104428793A (zh) 2015-03-18
JP6429134B2 (ja) 2018-11-28
US20150332116A1 (en) 2015-11-19
TW201407542A (zh) 2014-02-16
RU2015104066A (ru) 2016-08-27

Similar Documents

Publication Publication Date Title
JP6429134B2 (ja) 勾配ヒストグラムに基づいて画像記述子を変換する方法および関連する画像処理装置
CN103026368B (zh) 使用增量特征提取的对象辨识
Tsai et al. Location coding for mobile image retrieval
JP6044547B2 (ja) 局所特徴量抽出装置、局所特徴量抽出方法、及びプログラム
US8687892B2 (en) Generating a binary descriptor representing an image patch
US9747308B2 (en) Method and apparatus for searching an image, and computer-readable recording medium for executing the method
JP6035173B2 (ja) 画像検索システム及び画像検索方法
US9779320B2 (en) Image processing apparatus and method for encoding an image descriptor based on a gradient histogram
JP6042778B2 (ja) 画像に基づくバイナリ局所特徴ベクトルを用いた検索装置、システム、プログラム及び方法
CN110442749B (zh) 视频帧处理方法及装置
KR20140104789A (ko) 정지영상의 식별자 생성 방법 및 이를 수행하는 영상 처리 장치
KR101081459B1 (ko) 고차원 이진 데이터 검색 장치 및 방법
JP2012150590A (ja) 情報処理装置及び情報処理プログラム

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20180817

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20180828

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20180926

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20181010

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20181018

R150 Certificate of patent or registration of utility model

Ref document number: 6429134

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: R3D02

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250