(57)【要約】【課題】ランレングス符号化圧縮した画像データの離散的な画素値ヒストグラムから原画像の連続的な画素値ヒストグラムを近似復元する方法を提供する。近似復元した画素値ヒストグラムを用いて画像検索する方法および装置を提供する。【解決手段】圧縮画像データにおける画素値Ciに対応する画素値度数Liを算出することにより離散ヒストグラムデータを生成し、画素値Ci（i=1,2,…,M）に対して、当該画素値Ciに対応する前記離散ヒストグラムデータの前記画素値度数Liを、当該画素値Ciを中心とする標準偏差σの正規分布となるように当該画素値Ciの近傍の画素値xの画素値度数L（x）に分散配分することにより、前記原画像データの画素値の出現頻度を近似的に表す近似ヒストグラムデータを生成する。【選択図】        図３
2006139579.tif 000002
【特許請求の範囲】
【請求項１】
圧縮画像記憶手段に記憶された、原画像の画像データ（以下「原画像データ」という。）の各画素の画素値を量子化するとともにランレングス符号化することにより得られた圧縮画像の画像データ（以下「圧縮画像データ」という。）から、前記原画像の画素値ヒストグラムを近似的に復元するヒストグラム近似復元装置であって、
前記圧縮画像記憶手段に記憶された前記圧縮画像データの各画素値Ci（i=1,2,…,M，Mは圧縮画像データに含まれる全画素値の数）に対して、前記圧縮画像データの全体領域又は特定の部分領域における当該画素値に対応するランレングスの総和（以下「画素値度数」という。）Liを算出することにより離散ヒストグラムデータを生成する離散ヒストグラム生成手段と、
前記圧縮画像データの各画素値Ci（i=1,2,…,M）に対して、当該画素値Ciに対応する前記離散ヒストグラムデータの前記画素値度数Liを、当該画素値Ciを中心とする標準偏差σの正規分布となるように当該画素値Ciの近傍の画素値xの画素値度数L（x）に分散配分することにより、前記原画像データの画素値の出現頻度を近似的に表す近似ヒストグラムデータを生成する近似ヒストグラム生成手段と、
を有することを特徴とするヒストグラム近似復元装置。
【請求項２】
前記近似ヒストグラム生成手段は、
前記離散ヒストグラムデータ｛（Ci, Li）|i=1,2,…,M｝の各画素値度数Li（i=1,2,…,M）を、（数１）の正規分布関数G（x）で表される正規分布となるように画素値Ciの近傍の画素値xの画素値度数L（x）に分散配分することにより、前記近似ヒストグラムデータを生成すること
を特徴とする請求項１に記載のヒストグラム近似復元装置。
【数１】
2006139579.tif 000012
【請求項３】
前記近似ヒストグラム生成手段は、
前記圧縮画像データの各画素値Ci（i=1,2,…,M）に対して、当該画素値Ciに対応する前記離散ヒストグラムデータの前記画素値度数Liを、当該画素値Ciを中心とする標準偏差σ＝Ci／４の正規分布となるように当該画素値Ciの近傍の画素値xの画素値度数L（x）に分散配分することにより、前記原画像データの画素値の出現頻度を近似的に表す近似ヒストグラムデータを生成すること
を特徴とする請求項１又は２に記載のヒストグラム近似復元装置。
【請求項４】
前記圧縮画像データは、前記原画像データの各画素の画素値を量子化するとともに、各ランにおいて当該ランに対応する原画像データの画素値の分散が所定の閾値Γ以下となるようにランレングスを決定することによりランレングス符号化されたものであり、
前記近似ヒストグラム生成手段は、
前記圧縮画像データの各画素値Ci（i=1,2,…,M）に対して、当該画素値Ciに対応する前記離散ヒストグラムデータの前記画素値度数Liを、当該画素値Ciを中心とする標準偏差σ＝Γ１／２の正規分布となるように当該画素値Ciの近傍の画素値xの画素値度数L（x）に分散配分することにより、前記原画像データの画素値の出現頻度を近似的に表す近似ヒストグラムデータを生成すること
を特徴とする請求項１又は２に記載のヒストグラム近似復元装置。
【請求項５】
前記圧縮画像データは、前記原画像データの各画素の画素値を量子化するとともに、各ランにおいて当該ランに対応する原画像データの画素値の分散が所定の閾値Γ以下となるようにランレングスを決定することによりランレングス符号化されたものであり、
前記近似ヒストグラム生成手段は、  
前記標準偏差σが、Ci／４又は前記閾値Γの平方根のうちの何れか小さい値であること
を特徴とする請求項１又は２に記載のヒストグラム近似復元装置。
【請求項６】
原画像の画像データ（以下「原画像データ」という。）の各画素の画素値を量子化するとともにランレングス符号化することにより得られた圧縮画像の画像データ（以下「圧縮画像データ」という。）が複数記憶された画像データベース、
及び、基準画像の画素値ヒストグラムデータ（以下「基準ヒストグラムデータ」という。）が記憶された基準ヒストグラム記憶手段
を有し、前記画像データベースに記憶された前記圧縮画像データの中から前記基準画像に類似するものを検索する画像検索装置において、
前記画像データベースに記憶された前記圧縮画像データの各画素値Ci（i=1,2,…,M，Mは圧縮画像データに含まれる全画素値の数）に対して、前記圧縮画像データの全体領域又は特定の部分領域における当該画素値に対応するランレングスの総和（以下「画素値度数」という。）Liを算出することにより離散ヒストグラムデータを生成する離散ヒストグラム生成手段と、
前記圧縮画像データの各画素値Ci（i=1,2,…,M）に対して、当該画素値Ciに対応する前記離散ヒストグラムデータの前記画素値度数Liを、当該画素値Ciを中心とする標準偏差σの正規分布となるように当該画素値Ciの近傍の画素値xの画素値度数L（x）に分散配分することにより、前記原画像データの画素値の出現頻度を近似的に表す近似ヒストグラムデータを生成する近似ヒストグラム生成手段と、
前記基準ヒストグラム記憶手段に記憶されている前記基準ヒストグラムデータと前記近似ヒストグラムデータとの類似度を算出する類似度計算手段と、
前記各圧縮画像データの前記類似度に基づいて、前記基準画像に類似する一乃至複数の前記圧縮画像データを選出する画像選出手段と、
を備えていることを特徴とする画像検索装置。
【請求項７】
前記類似度計算手段は、前記近似ヒストグラムデータと前記基準ヒストグラムデータとの差分を前記類似度として算出すること
を特徴とする請求項６に記載の画像検索装置。
【請求項８】
原画像の画像データ（以下「原画像データ」という。）の各画素の画素値を量子化するとともにランレングス符号化することにより得られた圧縮画像の画像データ（以下「圧縮画像データ」という。）から、前記原画像の画素値ヒストグラムを近似的に復元するヒストグラム近似復元方法において、
圧縮画像記憶手段に記憶された前記圧縮画像データの各画素値Ci（i=1,2,…,M，Mは圧縮画像データに含まれる全画素値の数）に対して、前記圧縮画像データの全体領域又は特定の部分領域における当該画素値に対応するランレングスの総和（以下「画素値度数」という。）Liを算出することにより離散ヒストグラムデータを生成する第１ステップ、
及び、前記圧縮画像データの各画素値Ci（i=1,2,…,M）に対して、当該画素値Ciに対応する前記離散ヒストグラムデータの前記画素値度数Liを、当該画素値Ciを中心とする標準偏差σの正規分布となるように当該画素値Ciの近傍の画素値xの画素値度数L（x）に分散配分することにより、前記原画像データの画素値の出現頻度を近似的に表す近似ヒストグラムデータを生成する第２ステップ
を有することを特徴とするヒストグラム近似復元方法。
【請求項９】
前記第２ステップにおいて、
前記離散ヒストグラムデータ｛（Ci, Li）|i=1,2,…,M｝の各画素値度数Li（i=1,2,…,M）を、（数２）の正規分布関数G（x）で表される正規分布となるように画素値Ciの近傍の画素値xの画素値度数L（x）に分散配分することにより、前記近似ヒストグラムデータを生成すること
を特徴とする請求項８に記載のヒストグラム近似復元方法。
【数２】
2006139579.tif 000013
【請求項１０】
前記第２ステップにおいて、前記標準偏差σはCi／４であることを特徴とする請求項８又は９に記載のヒストグラム近似復元方法。
【請求項１１】
前記圧縮画像データは、原画像データの各画素の画素値を量子化するとともに、各ランにおいて当該ランに対応する原画像データの画素値の分散が所定の閾値Γ以下となるようにランレングスを決定することによりランレングス符号化されたものであり、
前記第２ステップにおいて、前記標準偏差σは前記閾値Γの平方根であること
を特徴とする請求項８又は９に記載のヒストグラム近似復元方法。
【請求項１２】
前記圧縮画像データは、原画像データの各画素の画素値を量子化するとともに、各ランにおいて当該ランに対応する原画像データの画素値の分散が所定の閾値Γ以下となるようにランレングスを決定することによりランレングス符号化されたものであり、
前記第２ステップにおいて、前記標準偏差σは、Ci／４又は前記閾値Γの平方根のうちの何れか小さい値であること
を特徴とする請求項８又は９に記載のヒストグラム近似復元方法。
【請求項１３】
原画像の画像データ（以下「原画像データ」という。）の各画素の画素値を量子化するとともにランレングス符号化することにより得られた圧縮画像の画像データ（以下「圧縮画像データ」という。）が複数記憶された画像データベース、
及び、基準画像の画素値ヒストグラムデータ（以下「基準ヒストグラムデータ」という。）が記憶された基準ヒストグラム記憶手段
を有するシステムにおいて、前記画像データベースに記憶された前記圧縮画像データの中から前記基準画像に類似するものを検索する画像検索方法であって、
前記画像データベースに記憶された前記圧縮画像データの各画素値Ci（i=1,2,…,M，Mは圧縮画像データに含まれる全画素値の数）に対して、前記圧縮画像データの全体領域又は特定の部分領域における当該画素値に対応するランレングスの総和（以下「画素値度数」という。）Liを算出することにより離散ヒストグラムデータを生成する離散ヒストグラム生成ステップ、
前記圧縮画像データの各画素値Ci（i=1,2,…,M）に対して、当該画素値Ciに対応する前記離散ヒストグラムデータの前記画素値度数Liを、当該画素値Ciを中心とする標準偏差σの正規分布となるように当該画素値Ciの近傍の画素値xの画素値度数L（x）に分散配分することにより、前記原画像データの画素値の出現頻度を近似的に表す近似ヒストグラムデータを生成する近似ヒストグラム生成ステップ、
及び、前記基準ヒストグラム記憶手段に記憶されている前記基準ヒストグラムデータと前記近似ヒストグラムデータとの類似度を算出する類似度計算ステップ
からなる画像比較過程を、前記画像データベースに記憶されているすべての圧縮画像データに対して実行する第１ステップと、
前記各圧縮画像データの前記類似度に基づいて前記基準画像に類似する一乃至複数の前記圧縮画像データを選出する第２ステップと、
を有することを特徴とする画像検索方法。
【請求項１４】
前記類似度計算ステップにおいて、
前記近似ヒストグラムデータと前記基準ヒストグラムデータとの差分を前記類似度として算出すること
を特徴とする請求項１３に記載の画像検索方法。
【請求項１５】
請求項８乃至請求項１２のヒストグラム近似復元方法をコンピュータに実行させることを特徴とするプログラム。
【請求項１６】
請求項１３又は請求項１４の画像検索方法をコンピュータに実行させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明はランレングス符号化された画像データから、原画像の画素値のヒストグラムを近似的に復元するヒストグラム近似復元技術、及び、近似復元された画素値ヒストグラムを用いて画像検索を行う画像検索技術に関する。
【背景技術】
【０００２】
多数の画像データの中から基準画像と同一の画像データを探し出す画像検索が、各種のアプリケーションで必要とされる。圧縮画像についても同様である。従来、圧縮画像について画像検索を行うためには、圧縮画像を元の画像データに戻して（復号して）、復号された画像データを基準画像のデータと直接比較する必要があったが、復号に時間を要するので、高速検索には不向きであった。そこで、近年は画像データを直接比較するのではなく、検索画像中の各画素値の出現頻度（画素値ヒストグラム）と基準画像の画素値ヒストグラムの類似度を求めて、画像検索を行う方法が各種提案されている（例えば、特許文献１）。
【特許文献１】特開２００２−１５０２８７号公報
【特許文献２】特開２００４−３１２６９３号公報
【特許文献３】特開２００３−３１２６９３号公報
【非特許文献１】鎌田清一郎，「ヒルベルト走査を利用した濃淡画像の情報圧縮に関する考察」，電子情報通信学会論文誌，Vol.J80-D-II，No.2，pp.426-433，１９９７年２月
【発明の開示】
【発明が解決しようとする課題】
【０００３】
しかしながら、ランレングス符号化された画像データは、画素値を量子化しているので、この画像データから画素値ヒストグラムを作成すると、図５に示すような離散的なヒストグラムになる。一方、基準画像のヒストグラムは図７に示すような連続的なヒストグラムになる。このような、外見が全く異なるヒストグラム同士を直接比較しても、一致点が少ないので、基準画像に類似する画像データを発見することは出来なかった。
【０００４】
そこで、本発明は、ランレングス符号化圧縮した画像データの離散的な画素値ヒストグラムから原画像の連続的な画素値ヒストグラムを近似復元する方法を提供することを目的とする。また、近似復元した画素値ヒストグラムを用いて画像検索する方法および装置を提供することを目的とする。
【課題を解決するための手段】
【０００５】
本発明に係るヒストグラム近似復元装置の第１の構成は、圧縮画像記憶手段に記憶された、原画像の画像データ（以下「原画像データ」という。）の各画素の画素値を量子化するとともにランレングス符号化することにより得られた圧縮画像の画像データ（以下「圧縮画像データ」という。）から、前記原画像の画素値ヒストグラムを近似的に復元するヒストグラム近似復元装置であって、前記圧縮画像記憶手段に記憶された前記圧縮画像データの各画素値Ci（i=1,2,…,M，Mは圧縮画像データに含まれる全画素値の数）に対して、前記圧縮画像データの全体領域又は特定の部分領域における当該画素値に対応するランレングスの総和（以下「画素値度数」という。）Liを算出することにより離散ヒストグラムデータを生成する離散ヒストグラム生成手段と、前記圧縮画像データの各画素値Ci（i=1,2,…,M）に対して、当該画素値Ciに対応する前記離散ヒストグラムデータの前記画素値度数Liを、当該画素値Ciを中心とする標準偏差σの正規分布となるように当該画素値Ciの近傍の画素値xの画素値度数L（x）に分散配分することにより、前記原画像データの画素値の出現頻度を近似的に表す近似ヒストグラムデータを生成する近似ヒストグラム生成手段と、を有することを特徴とする。
【０００６】
この構成によれば、離散ヒストグラム生成手段は、圧縮画像データの各画素値Ciに対し、上記画素値度数Liを算出することにより、圧縮画像データの復元を行うことなく圧縮画像データの全体領域又は特定部分領域における画素値のヒストグラム（離散ヒストグラムデータ）を算出することができる。すなわち、復号動作が伴わないので、圧縮画像データから高速に離散ヒストグラムデータを再生できる。更に、近似ヒストグラム生成手段は、画素値に対して離散的な画素値度数を持った離散ヒストグラムデータに対して、上記分散配分を行うことによって、画素値に対して連続的な画素値度数を持つ近似ヒストグラムデータに変換する。これにより原画像データに近いヒストグラムデータを復元することができる。
【０００７】
本発明に係るヒストグラム近似復元装置の第２の構成は、前記第１の構成において、前記近似ヒストグラム生成手段は、前記離散ヒストグラムデータ｛（Ci, Li）|i=1,2,…,M｝の各画素値度数Li（i=1,2,…,M）を、（数１）の正規分布関数G（x）で表される正規分布となるように画素値Ciの近傍の画素値xの画素値度数L（x）に分散配分することにより、前記近似ヒストグラムデータを生成することを特徴とする。
【０００８】
【数１】
2006139579.tif 000003
【０００９】
この構成によれば、上記正規分布関数G（x）は離散ヒストグラムデータ｛（Ci, Li）|i=1,2,…,M｝の各画素値Ci、画素値度数Liおよび標準偏差σの関数なので、標準偏差σを仮定するだけで、原画像データに近いヒストグラムデータを復元することができる。
【００１０】
本発明に係るヒストグラム近似復元装置の第３の構成は、前記第１又は２の構成において、前記近似ヒストグラム生成手段は、前記圧縮画像データの各画素値Ci（i=1,2,…,M）に対して、当該画素値Ciに対応する前記離散ヒストグラムデータの前記画素値度数Liを、当該画素値Ciを中心とする標準偏差σ＝Ci／４の正規分布となるように当該画素値Ciの近傍の画素値xの画素値度数L（x）に分散配分することにより、前記原画像データの画素値の出現頻度を近似的に表す近似ヒストグラムデータを生成することを特徴とする。
【００１１】
この構成によれば、標準偏差σの値をCi／４とするので、離散ヒストグラムデータ｛（Ci, Li）|i=1,2,…,M｝から一意的に上記正規分布を決定することができる。これにより、ヒストグラムデータの復元を更に高速化できる。
【００１２】
本発明に係るヒストグラム近似復元装置の第４の構成は、前記第１又は２の構成において、前記圧縮画像データは、前記原画像データの各画素の画素値を量子化するとともに、各ランにおいて当該ランに対応する原画像データの画素値の分散が所定の閾値Γ以下となるようにランレングスを決定することによりランレングス符号化されたものであり、前記近似ヒストグラム生成手段は、前記圧縮画像データの各画素値Ci（i=1,2,…,M）に対して、当該画素値Ciに対応する前記離散ヒストグラムデータの前記画素値度数Liを、当該画素値Ciを中心とする標準偏差σ＝Γ１／２の正規分布となるように当該画素値Ciの近傍の画素値xの画素値度数L（x）に分散配分することにより、前記原画像データの画素値の出現頻度を近似的に表す近似ヒストグラムデータを生成することを特徴とする。
【００１３】
この構成によれば、各ランに対応する原画像データの画素値の分散の閾値Γの平方根を標準偏差σの値とするので、標準偏差σの精度の良い推定値が得られ、ヒストグラムデータの近似復元の精度が向上する。
【００１４】
本発明に係るヒストグラム近似復元装置の第５の構成は、前記第１の構成において、前記圧縮画像データは、前記原画像データの各画素の画素値を量子化するとともに、各ランにおいて当該ランに対応する原画像データの画素値の分散が所定の閾値Γ以下となるようにランレングスを決定することによりランレングス符号化されたものであり、前記近似ヒストグラム生成手段は、前記標準偏差σが、Ci／４又は前記閾値Γの平方根のうちの何れか小さい値であることを特徴とする。
【００１５】
この構成によれば、標準偏差σの値をCi／４又は前記閾値Γの平方根のうちの何れか小さい値とするので、標準偏差σの値が過大に推定されるのを避けることができる。これにより、ヒストグラムデータの近似復元の精度が向上する。
【００１６】
本発明に係る画像検索装置の第１の構成は、原画像の画像データ（以下「原画像データ」という。）の各画素の画素値を量子化するとともにランレングス符号化することにより得られた圧縮画像の画像データ（以下「圧縮画像データ」という。）が複数記憶された画像データベース、及び、基準画像の画素値ヒストグラムデータ（以下「基準ヒストグラムデータ」という。）が記憶された基準ヒストグラム記憶手段を有し、前記画像データベースに記憶された前記圧縮画像データの中から前記基準画像に類似するものを検索する画像検索装置において、前記画像データベースに記憶された前記圧縮画像データの各画素値Ci（i=1,2,…,M，Mは圧縮画像データに含まれる全画素値の数）に対して、前記圧縮画像データの全体領域又は特定の部分領域における当該画素値に対応するランレングスの総和（以下「画素値度数」という。）Liを算出することにより離散ヒストグラムデータを生成する離散ヒストグラム生成手段と、前記圧縮画像データの各画素値Ci（i=1,2,…,M）に対して、当該画素値Ciに対応する前記離散ヒストグラムデータの前記画素値度数Liを、当該画素値Ciを中心とする標準偏差σの正規分布となるように当該画素値Ciの近傍の画素値xの画素値度数L（x）に分散配分することにより、前記原画像データの画素値の出現頻度を近似的に表す近似ヒストグラムデータを生成する近似ヒストグラム生成手段と、前記基準ヒストグラム記憶手段に記憶されている前記基準ヒストグラムデータと前記近似ヒストグラムデータとの類似度を算出する類似度計算手段と、前記各圧縮画像データの前記類似度に基づいて、前記基準画像に類似する一乃至複数の前記圧縮画像データを選出する画像選出手段と、を備えていることを特徴とする。
【００１７】
この構成によれば、離散ヒストグラム生成手段は、圧縮画像データの各画素値Ciに対し、上記画素値度数Liを算出することにより、圧縮画像データの復元を行うことなく圧縮画像データの全体領域又は特定部分領域における画素値のヒストグラム（離散ヒストグラムデータ）を算出することができる。すなわち、基準ヒストグラムデータ、圧縮画像データから高速に離散ヒストグラムデータを再生できる。更に、近似ヒストグラム生成手段は、画素値に対して離散的な画素値度数を持った離散ヒストグラムデータに対して、上記分散配分を行うことによって、画素値に対して連続的な画素値度数を持つ近似ヒストグラムデータに変換する。そして、前記近似ヒストグラムデータと基準ヒストグラムデータとの類似度が算出される。これにより、圧縮画像データの復元を行うことなく原画像データに近いヒストグラムデータを復元して、基準ヒストグラムデータとの類似度を算出できるので、圧縮画像データに対する画像検索を高速化できる。
【００１８】
本発明に係る画像検索装置の第２の構成は、前記第１の構成において、前記類似度計算手段は、前記近似ヒストグラムデータと前記基準ヒストグラムデータとの差分を前記類似度として算出することを特徴とする。
【００１９】
この構成によれば、類似度計算手段は、前記近似ヒストグラムデータと前記基準ヒストグラムデータとの差分を前記類似度として算出するので、類似度の算出が単純になり、画像検索に要する時間を更に短縮できる。
【００２０】
本発明に係るヒストグラム近似復元方法の第１の構成は、原画像の画像データ（以下「原画像データ」という。）の各画素の画素値を量子化するとともにランレングス符号化することにより得られた圧縮画像の画像データ（以下「圧縮画像データ」という。）から、前記原画像の画素値ヒストグラムを近似的に復元するヒストグラム近似復元方法において、圧縮画像記憶手段に記憶された前記圧縮画像データの各画素値Ci（i=1,2,…,M，Mは圧縮画像データに含まれる全画素値の数）に対して、前記圧縮画像データの全体領域又は特定の部分領域における当該画素値に対応するランレングスの総和（以下「画素値度数」という。）Liを算出することにより離散ヒストグラムデータを生成する第１ステップ、及び、前記圧縮画像データの各画素値Ci（i=1,2,…,M）に対して、当該画素値Ciに対応する前記離散ヒストグラムデータの前記画素値度数Liを、当該画素値Ciを中心とする標準偏差σの正規分布となるように当該画素値Ciの近傍の画素値xの画素値度数L（x）に分散配分することにより、前記原画像データの画素値の出現頻度を近似的に表す近似ヒストグラムデータを生成する第２ステップを有することを特徴とする。
【００２１】
本発明に係るヒストグラム近似復元方法の第２の構成は、前記第１の構成において、前記第２ステップにおいて、前記離散ヒストグラムデータ｛（Ci, Li）|i=1,2,…,M｝の各画素値度数Li（i=1,2,…,M）を、（数２）の正規分布関数G（x）で表される正規分布となるように画素値Ciの近傍の画素値xの画素値度数L（x）に分散配分することにより、前記近似ヒストグラムデータを生成することを特徴とする。
【００２２】
【数２】
2006139579.tif 000004
【００２３】
本発明に係るヒストグラム近似復元方法の第３の構成は、前記第１又は２の構成において、前記第２ステップで、前記標準偏差σはCi／４であることを特徴とする。
【００２４】
本発明に係るヒストグラム近似復元方法の第４の構成は、前記第１又は２の構成において、前記圧縮画像データは、原画像データの各画素の画素値を量子化するとともに、各ランにおいて当該ランに対応する原画像データの画素値の分散が所定の閾値Γ以下となるようにランレングスを決定することによりランレングス符号化されたものであり、前記第２ステップにおいて、前記標準偏差σは前記閾値Γの平方根であることを特徴とする。
【００２５】
本発明に係るヒストグラム近似復元方法の第５の構成は、前記第１又は２の構成において、前記圧縮画像データは、原画像データの各画素の画素値を量子化するとともに、各ランにおいて当該ランに対応する原画像データの画素値の分散が所定の閾値Γ以下となるようにランレングスを決定することによりランレングス符号化されたものであり、前記第２ステップにおいて、前記標準偏差σは、Ci／４又は前記閾値Γの平方根のうちの何れか小さい値であることを特徴とする。
【００２６】
本発明に係る画像検索方法の第１の構成は、原画像の画像データ（以下「原画像データ」という。）の各画素の画素値を量子化するとともにランレングス符号化することにより得られた圧縮画像の画像データ（以下「圧縮画像データ」という。）が複数記憶された画像データベース、及び、基準画像の画素値ヒストグラムデータ（以下「基準ヒストグラムデータ」という。）が記憶された基準ヒストグラム記憶手段を有するシステムにおいて、前記画像データベースに記憶された前記圧縮画像データの中から前記基準画像に類似するものを検索する画像検索方法であって、前記画像データベースに記憶された前記圧縮画像データの各画素値Ci（i=1,2,…,M，Mは圧縮画像データに含まれる全画素値の数）に対して、前記圧縮画像データの全体領域又は特定の部分領域における当該画素値に対応するランレングスの総和（以下「画素値度数」という。）Liを算出することにより離散ヒストグラムデータを生成する離散ヒストグラム生成ステップ、前記圧縮画像データの各画素値Ci（i=1,2,…,M）に対して、当該画素値Ciに対応する前記離散ヒストグラムデータの前記画素値度数Liを、当該画素値Ciを中心とする標準偏差σの正規分布となるように当該画素値Ciの近傍の画素値xの画素値度数L（x）に分散配分することにより、前記原画像データの画素値の出現頻度を近似的に表す近似ヒストグラムデータを生成する近似ヒストグラム生成ステップ、及び、前記基準ヒストグラム記憶手段に記憶されている前記基準ヒストグラムデータと前記近似ヒストグラムデータとの類似度を算出する類似度計算ステップからなる画像比較過程を、前記画像データベースに記憶されているすべての圧縮画像データに対して実行する第１ステップと、前記各圧縮画像データの前記類似度に基づいて前記基準画像に類似する一乃至複数の前記圧縮画像データを選出する第２ステップと、を有することを特徴とする。
【００２７】
本発明に係る画像検索方法の第２の構成は、前記第１の構成において、前記類似度計算ステップで、前記近似ヒストグラムデータと前記基準ヒストグラムデータとの差分を前記類似度として算出することを特徴とする。
【００２８】
本発明に係るプログラムは、上記第１乃至第５の構成のヒストグラム近似復元方法をコンピュータに実行されることを特徴とする。
【００２９】
本発明に係るプログラムは、上記第１又は第２の構成の画像検索方法をコンピュータに実行させることを特徴とする。
【発明の効果】
【００３０】
以上説明したように、本発明はランレングス符号化圧縮した画像データの離散的な画素値ヒストグラムの近傍に画素が正規分布するように画素を配分して、離散的な画素値ヒストグラムから原画像の連続的な画素値ヒストグラムを近似復元するので、原画像の画素値ヒストグラムを効率よく、かつ精度良く復元する方法を提供することができる。また、近似復元した画素値ヒストグラムと基準画像の画素値ヒストグラムを比較して類似度を求めるので画像データの検索を高速処理できる。このような効果により、本発明は画像検索の高速化と高精度化に資するものである。
【発明を実施するための最良の形態】
【００３１】
以下、本発明を実施するための最良の形態について、図面を参照しながら説明する。
【実施例１】
【００３２】
図１は、本発明の実施例１に係るヒストグラム近似復元装置及び画像検索装置の構成を表すブロック図である。
【００３３】
本実施例の画像検索装置１は、ヒストグラム近似復元装置２、フレーム・メモリ３、領域抽出手段４、領域圧縮画像記憶手段５、ヒストグラム作成手段６、基準ヒストグラム記憶手段７、画像データベース（圧縮画像記憶手段）８、領域抽出手段９、類似度計算手段１０、及び画像選出手段１１を備えている。
【００３４】
尚、本実施例の画像検索装置１は、ＣＣＤカメラ２０により撮像された微粒子の画像と類似する画像を画像データベース８に記憶されている微粒子画像と照合することにより、微粒子の種類を判別するシステム等に用いられる。
【００３５】
フレーム・メモリ３は、ＣＣＤカメラ２０により撮像され、ＡＤ変換器２１によりデジタル化された画像を記憶するメモリである。領域抽出手段４は、フレーム・メモリ３に記憶された画像中から部分画像（例えば、微粒子部分の画像）を抽出し、領域圧縮画像記憶手段５に保存する。
【００３６】
ヒストグラム作成手段６は、領域圧縮画像記憶手段５に記憶された部分画像のデータ（以下、「基準画像」という。）に含まれる各画素の画素値についてヒストグラム（以下、「基準ヒストグラムデータ」という。）を生成して、基準ヒストグラム記憶手段７に保存する。
【００３７】
画像データベース８には、複数の圧縮画像データが記憶されている。ここで、圧縮画像データは、原画像の画像データ（以下「原画像データ」という。）の各画素の画素値を量子化するとともにランレングス符号化することにより得られた圧縮画像の画像データである。尚、このランレングス符号化の詳細に関しては後述する。
【００３８】
領域抽出手段９は、画像データベース８から圧縮画像データを読み出して、その全部を出力し又は一部の領域を切り出して出力する。これは、例えば、圧縮画像データの一部領域に微粒子画像がある場合などに、その微粒子画像の部分のみを抽出するためのものである。
【００３９】
ヒストグラム近似復元装置２は、領域抽出手段９により画像データベース８から読み出された圧縮画像データ又はその部分データ（以下「領域圧縮画像データ」という。）から、原画像又はその部分画像の画素値ヒストグラムを近似的に復元する。ヒストグラム近似復元装置２は、領域圧縮画像記憶手段１２、離散ヒストグラム生成手段１３、及び近似ヒストグラム生成手段１４を備えている。
【００４０】
領域圧縮画像記憶手段１２は、領域抽出手段９が出力する領域圧縮画像データを一時的に記憶する。離散ヒストグラム生成手段１３は、領域圧縮画像記憶手段１２に記憶された領域圧縮画像データから離散ヒストグラムデータを生成する。離散ヒストグラムデータについては後述する。近似ヒストグラム生成手段１４は、領域圧縮画像記憶手段１２が出力する離散ヒストグラムデータに対して、画素値度数に分散をもたせることによって近似ヒストグラムデータを生成する。
【００４１】
類似度計算手段１０は、基準ヒストグラム記憶手段７に記憶されている基準ヒストグラムデータと、近似ヒストグラム生成手段１４が出力する近似ヒストグラムデータとの類似度を算出する。画像選出手段１１は、各圧縮画像データの前記類似度に基づいて、基準画像に類似する一乃至複数の圧縮画像データを選出し、それを出力装置２２に出力する。尚、出力装置２２には、磁気ディスクなどの外部記憶装置、ネットワーク上に出力する通信インタフェース等が使用される。
【００４２】
尚、本実施例の画像検索装置１は、ＬＳＩチップとしてハードウェア的に構成してもよいが、汎用のコンピュータにプログラムをロードすることによって各構成部分を機能モジュールとして構成してもよい。
【００４３】
次に、画像データベース８に記憶された圧縮画像データについて説明する。画像データベース８に記憶された圧縮画像データは、原画像データの各画素の画素値を量子化するとともに、各ランにおいて当該ランに対応する原画像データの画素値の分散が所定の閾値Γ以下となるようにランレングスを決定することによりランレングス符号化されたものである。この圧縮アルゴリズムについては、特許文献２，３又は非特許文献１に詳細に記載されている。ここではその圧縮アルゴリズムの一例を説明するに止める。
【００４４】
まず、２次元の原画像の画像データ（以下「原画像データ」という。）を走査して１次元の原画像データとする。この１次元化された原画像データを{xi ; i=1,2,…,N}と記す。ここで、Ｎは原画像データの画素数である。
【００４５】
次に、この１次元化された原画像データに対して、以下のような再帰的な２分割処理を実行する。
【００４６】
まず、区間ｌの１次元画素列を{xi ; i=1,2,…,N}とする。区間Ｌを区間Ｌ１と区間Ｌ２に分割する。図２にその分割例を示す。区間Ｌ，Ｌ１，Ｌ２の平均値は、それぞれ式（３），式（４），式（５）により表される。但し、区間Ｌ１の画素数をＮ１とする。
【００４７】
【数３】
2006139579.tif 000005
【００４８】
再帰的な２分割を行う上で、区間Ｌ，Ｌ１，Ｌ２の累積自乗誤差を、それぞれｅ，ｅ１，ｅ２とすると、これらは式（６），式（７），式（８）により表される。
【００４９】
【数４】
2006139579.tif 000006
【００５０】
二分割点は、次の評価式が最小となるＮ１（１＜Ｎ１＜Ｎ）とする。
【００５１】
【数５】
2006139579.tif 000007
【００５２】
この２分割操作を、Evalが所定の閾値Γ以下且つＮがＮｍａｘ以下となるまで再帰的に繰り返して行う。
【００５３】
以上のようなアルゴリズムにより、１次元化された原画像データは複数の区間に分割される。分割された各区間をランとし、その区間（ラン）における画素の平均値をそのランにおける画素値としてランレングス符号化を行う。このランレングス符号化によって、１次元化された原画像データは、各ランにおいて画素値が量子化され、各ランにおいて当該ランに対応する原画像データの画素値の分散が所定の閾値Γ以下となる。このようにして得られた圧縮画像データが、画像データベース８に複数保存されている。
【００５４】
図３は、ランレングス符号化により圧縮された画像データの例を示す模式図である。図３では、横軸が走査の順を表し、縦軸が画素値（輝度）を表す。図中に括弧書きした数字の組は、前の数字が画素値（ここでは輝度）、後の数字がランレングスを表す。一般に、圧縮画像データは下記の式（１０）で表される。但し、Ｃｉは画素値、ｌｉはランレングス、Ｎｌはランの数である。
【００５５】
【数６】
2006139579.tif 000008
【００５６】
次に、上述のように構成された本実施例に係る画像検索装置１について、その動作を説明する。
【００５７】
最初に、前処理として、撮像された画像データから基準ヒストグラムデータの作成が行われる。まず、ＣＣＤカメラ２０により、微粒子等の被写体が撮像される。ＣＣＤカメラ２０から出力される画像データは、ＡＤ変換器２１においてデジタル化され、フレーム・メモリ３に保存される。領域抽出手段４は、フレーム・メモリ３に記憶された画像から目的の領域を抽出して基準画像として領域圧縮画像記憶手段５に保存する。尚、領域抽出手段４による領域抽出方法は、種々の公知の方法を使用することができるが、本発明とは直接関係がないので説明は省略する。
【００５８】
次に、ヒストグラム作成手段６は、領域圧縮画像記憶手段５に記憶された基準画像から画素値のヒストグラムを作成し、基準ヒストグラムデータとして基準ヒストグラム記憶手段７に保存する。
【００５９】
以上のような前処理が終了した後に、画像データベース８に記憶された圧縮画像データから、基準画像に類似する画像の圧縮画像データを検索する画像検索処理を以下のように実行する。
【００６０】
図４は、画像検索処理の流れを表すフローチャートである。ます、最初に画像検索装置１は、画像データベース８から読み出す圧縮画像データの番号ｉを１に初期化する（Ｓ１）。
【００６１】
次に、領域抽出手段９は、画像データベース８からｉ番目の圧縮画像データを読み出す（Ｓ２）。そして、読み出した圧縮画像データから所定の領域（領域圧縮画像データ）を抽出する（Ｓ３）。ここで、抽出する領域は目的に応じて圧縮画像の一部分とされたり、圧縮画像全体とされたりする。抽出された領域圧縮画像データは、領域圧縮画像記憶手段１２に保存される。
【００６２】
次に、離散ヒストグラム生成手段１３は、領域圧縮画像記憶手段１２に記憶された領域圧縮画像データから離散ヒストグラムデータを生成する（Ｓ４）。
【００６３】
図５は、領域圧縮画像データから作成した基準ヒストグラムである。図５で、横軸は画素値（輝度）、縦軸は画素の出現頻度を表す。画素の出現頻度は、画素値が等しい区間のランレングスの総和である。例えば、画素値７０の画素の出現頻度とは、領域圧縮画像データ全体から、画素値７０を持つデータを全て拾い出して、それらのデータのランレングスを合計したものである。領域圧縮画像データにおいて画素値がＣｉであるランの集合をＬ（Ｃｉ）と記す。基準ヒストグラムは、一般的に、次式で表される。但し、Ｌｉは画素値度数（ランレングスの総和）、Ｍは画素値の個数である。
【００６４】
【数７】
2006139579.tif 000009
【００６５】
なお、画素の出現頻度は各画素値の画素値度数を一画面の画素の総数で除して規格化した値で表してもよい。
【００６６】
次に、近似ヒストグラム生成手段１４は、領域圧縮画像記憶手段１２が出力する離散ヒストグラムデータに対して、画素値度数に分散をもたせることによって近似ヒストグラムデータを生成する（Ｓ５）。
【００６７】
図５の画素値ヒストグラムの画素値Ｃｉは量子化された値であるが、原画像では画素値Ｃｉの近傍に画素が分布していたと考えられる。そこで、これらの画素が画素値Ｃｉを中心に正規分布に従って分布していると仮定し、各画素値Ｃｉの画素値度数Ｌｉを式（１３）の分布関数G(x)に従って、画素値Ｃｉの近傍に分散配分する。但し、ｘは画素値、σは標準偏差を表す。
【００６８】
【数８】
2006139579.tif 000010
【００６９】
離散的な画素値ヒストグラムの全ての画素値Ｃｉについて、式（１３）により画素値度数Ｌｉを画素値Ｃｉの近傍に配分すると、図６に示すように、画素値Ｃｉの近傍に分布する画素を示すヒストグラム３１〜３４が描ける。なお、標準偏差σは画素値Ｃｉの４分の１（Ｃｉ／４）を選ぶと良い結果が得られた。最後にこれらのヒストグラムを全て重畳すると、全ての画素値についての連続的な近似ヒストグラム３５が得られる。
【００７０】
次に、類似度計算手段１０は、近似ヒストグラム生成手段１４が出力する近似ヒストグラムデータと、基準ヒストグラム記憶手段７に記憶されている基準ヒストグラムデータとの類似度Ｓｉを算出する（Ｓ６）。類似度Ｓｉは、近似ヒストグラムデータと基準ヒストグラムデータとの差分により計算される。すなわち、近似ヒストグラムデータをＨｉ（ｘ）、基準ヒストグラムデータをＨ０（ｘ）とすると、類似度Ｓｉは式（１４）により表される。但し、ｘ，ｘｉは画素値を表す。
【００７１】
【数９】
2006139579.tif 000011
【００７２】
次に、画像検索装置１は、画像データベース８に記憶されているすべての圧縮画像データについて、類似度Ｓｉの計算が完了したか否かを判定し（Ｓ７）、未了であれば、ｉを１だけ増加させて（Ｓ８）、ステップＳ２の処理に戻る。
【００７３】
ステップＳ７で、すべての圧縮画像データについて類似度Ｓｉの計算が完了した場合、画像選出手段１１は、各圧縮画像データの前記類似度に基づいて、基準画像に類似する一乃至複数の圧縮画像データを選出し（Ｓ９）、それを出力装置２２に出力して（Ｓ１０）、画像検索処理を終了する。
【００７４】
尚、画像選出手段１１による圧縮画像データを選出の方法としては、類似度Ｓｉが最小の圧縮画像データを選出する方法や、類似度Ｓｉが所定の閾値以下の圧縮画像データをすべて選出する方法等が考えられる。何れの方法を採用するかは、用途に応じて決めればよい。
【００７５】
以上のように、本実施例に係る画像検索装置１では、圧縮画像データから直接近似復元した近似ヒストグラムデータと基準画像の基準ヒストグラムデータの類似度を求め、画像の検索を行うことが出来る。画像のヒストグラムは生の画像データに比べてデータ量が格段少ないので、画像データを直接検索する場合に比べて、１０倍以上の速さで高速に検索することができる。
【００７６】
尚、特許文献２及び非特許文献１に記載されたように、原画像データの各画素の画素値を量子化するとともに、各ランにおいて当該ランに対応する原画像データの画素値の分散が所定の閾値Γ以下となるようにランレングスを決定することによりランレングス符号化された画像データの画素値ヒストグラムを近似復元する場合には、前記閾値Γの平方根を前記標準偏差σとすることもできる。
【００７７】
あるいは、前記閾値Γの平方根又は画素値Ｃｉの４分の１（Ｃｉ／４）のいずれか小さい値を標準偏差σとすることもできる。
【図面の簡単な説明】
【００７８】
【図１】本発明の実施例１に係るヒストグラム近似復元装置及び画像検索装置の構成を表すブロック図である。
【図２】区間分割の例を示す図である。
【図３】ランレングス符号化圧縮された画像データの例を示す模式図である。
【図４】画像検索処理の流れを表すフローチャートである。
【図５】図１の画像データから作成した離散的な画素値ヒストグラムである。
【図６】近似復元した連続的な画素値ヒストグラムである。
【図７】基準画像の基準ヒストグラムの例である。
【符号の説明】
【００７９】
１  画像検索装置
２  ヒストグラム近似復元装置
３  フレーム・メモリ
４  領域抽出手段
５  領域圧縮画像記憶手段
６  ヒストグラム作成手段
７  基準ヒストグラム記憶手段
８  画像データベース（圧縮画像記憶手段）
９  領域抽出手段
１０  類似度計算手段
１１  画像選出手段
１２  領域圧縮画像記憶手段
１３  離散ヒストグラム生成手段
１４  近似ヒストグラム生成手段
２０  ＣＣＤカメラ
２１  ＡＤ変換器
２２  出力装置
３１〜３４  正規分布に従って配分された画素を示すヒストグラム  
３５  全ての画素値についての連続的なヒストグラム  
【図１】
2006139579.tif 000014
【図２】
2006139579.tif 000015
【図３】
2006139579.tif 000016
【図４】
2006139579.tif 000017
【図５】
2006139579.tif 000018
【図６】
2006139579.tif 000019
【図７】
2006139579.tif 000020