(57)【要約】【課題】効率的且つ高精度に空間・階調削減画像を得る。【解決手段】空間・階調削減装置１は、原画像をｊ個（ｊ≧１）のパラメータを用いて階調削減処理してｊ枚の階調削減画像を生成する階調削減部１０と、階調削減画像をｋ個（ｋ≧１）のパラメータを用いて空間縮小処理してｊ×ｋ枚（但し、ｊ×ｋ≧２）の空間・階調削減画像を生成する空間縮小部２０と、ｊ×ｋ枚の空間・階調削減画像を前記原画像と比較して、最も劣化の少ない空間・階調削減画像を選択する最適化部３０と、を備える。【選択図】図２
2014238657.tif 000002
【特許請求の範囲】
【請求項１】
原画像の階調削減処理及び空間縮小処理を行って空間・階調削減画像を生成する空間・階調削減装置であって、
原画像をｊ個（ｊ≧１）のパラメータを用いて階調削減処理してｊ枚の階調削減画像を生成する階調削減部と、
前記階調削減画像をｋ個（ｋ≧１）のパラメータを用いて空間縮小処理してｊ×ｋ枚（但し、ｊ×ｋ≧２）の空間・階調削減画像を生成する空間縮小部と、
前記ｊ×ｋ枚の空間・階調削減画像を前記原画像と比較して、最も劣化の少ない空間・階調削減画像を選択する最適化部と、
を備えることを特徴とする空間・階調削減装置。
【請求項２】
前記階調削減部は、
ｊ個の雑音閾値を用いて、原画像から孤立点を除去したｊ枚の雑音除去画像を生成するか、又はｊ個のグラデーション閾値を用いて、原画像のグラデーション領域のみからなるｊ枚のグラデーション画像を生成する画像生成部と、
前記雑音除去画像又は前記原画像のヒストグラムに対して、前記グラデーション画像の階調値の度数が高くなるように重み付けしたヒストグラムを訓練データとして用いるLloyd-Max法により、前記ｊ枚の階調削減画像を生成する階調決定部と、
を備えることを特徴とする請求項１に記載の空間・階調削減装置。
【請求項３】
前記雑音閾値又は前記グラデーション閾値は、原画像を周波数分解して複数の周波数帯域成分を生成し、該周波数帯域成分を解析することにより決定されることを特徴とする、請求項２に記載の空間・階調削減装置。
【請求項４】
前記雑音閾値又は前記グラデーション閾値は、空間縮小率又は階調削減率と対応付けた対応表に基づいて決定され、
前記対応表は、前記空間縮小率が高くなるほど、前記雑音閾値又は前記グラデーション閾値が小さくなるように対応付けているか、あるいは、前記階調削減率が高くなるほど、前記雑音閾値又は前記グラデーション閾値が大きくなるように対応付けていることを特徴とする、請求項２に記載の空間・階調削減装置。
【請求項５】
前記空間縮小部は、解像度の縮小劣化過程を模擬した関数、又はウェーブレット分解を用いて前記空間縮小画像を生成することを特徴とする、請求項１から４のいずれか一項に記載の空間・階調削減装置。
【請求項６】
前記最適化部は、前記原画像に対する空間・階調削減画像のパワースペクトル差分値及びヒストグラム差分値を算出し、両者の合計値が最も小さい空間・階調削減画像を選択することを特徴とする、請求項１から５のいずれか一項に記載の空間・階調削減装置。
【請求項７】
コンピュータを、請求項１から６のいずれか一項に記載の空間・階調削減装置として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、原画像の解像度（標本化周波数）の縮小（空間縮小）及び階調数の削減（階調削減）を行い、空間・階調削減画像を生成する空間・階調削減装置及びプログラムに関する。
【背景技術】
【０００２】
従来、原画像の空間縮小処理を行う技術や、原画像の階調削減処理を行う技術が知られている。空間縮小処理は、例えば画素の単純間引きやウェーブレット分解などにより行うことができる。階調削減処理は、例えば階調数の単純削減やLloyd-Max法などにより行うことができる。Lloyd-Max法とは、画像の階調値(輝度値)ごとの度数を示すヒストグラムを調べ、階調削減を行う際に、度数が高い階調値ほど多くの階調を割り当てるようにする手法である（非特許文献１及び２参照）。
【０００３】
その他、空間縮小した後でＲＧＢ空間からＹＩＱ空間へ変換後、量子化及び符号化を行うことを特徴とし、色空間をＲＧＢからＹＩＱに変換することで，人の目につきにくい色差成分等の情報量を削減する方法が知られている（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開平３−２６７８８０号公報
【非特許文献】
【０００５】
【非特許文献１】S.P.Lloyd, “Least squares quantization in PCM”, IEEE Trans. Information Theory, vol.IT-28, pp.129-136, March 1982
【非特許文献２】J.Max, “Quantizing for minimum distortion”, IEEE Trans. Information Theory, vol.IT-7, pp.7-12, March 1960
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかし、原画像の空間縮小処理及び階調削減処理を行うには、空間縮小処理と階調削減処理とを別々に適用する必要がある。このため、効率的な処理を行うことができず、また、空間縮小処理と階調削減処理とを一元的に最適化することはできなかった。
【０００７】
かかる事情に鑑みてなされた本発明の目的は、効率的且つ高精度に空間・階調削減画像を得ることができる空間・階調削減装置及びプログラムを提供することにある。
【課題を解決するための手段】
【０００８】
上記課題を解決するため、本発明に係る空間・階調削減装置は、原画像の階調削減処理及び空間縮小処理を行って空間・階調削減画像を生成する空間・階調削減装置であって、原画像をｊ個（ｊ≧１）のパラメータを用いて階調削減処理してｊ枚の階調削減画像を生成する階調削減部と、前記階調削減画像をｋ個（ｋ≧１）のパラメータを用いて空間縮小処理してｊ×ｋ枚（但し、ｊ×ｋ≧２）の空間・階調削減画像を生成する空間縮小部と、前記ｊ×ｋ枚の空間・階調削減画像を前記原画像と比較して、最も劣化の少ない空間・階調削減画像を選択する最適化部と、を備えることを特徴とする。
【０００９】
さらに、本発明に係る空間・階調削減装置において、前記階調削減部は、ｊ個の雑音閾値を用いて、原画像から孤立点を除去したｊ枚の雑音除去画像を生成するか、又はｊ個のグラデーション閾値を用いて、原画像のグラデーション領域のみからなるｊ枚のグラデーション画像を生成する画像生成部と、前記雑音除去画像又は前記原画像のヒストグラムに対して、前記グラデーション画像の階調値の度数が高くなるように重み付けしたヒストグラムを訓練データとして用いるLloyd-Max法により、前記ｊ枚の階調削減画像を生成する階調決定部と、を備えることを特徴とする。
【００１０】
さらに、本発明に係る空間・階調削減装置において、前記雑音閾値又は前記グラデーション閾値は、原画像を周波数分解して複数の周波数帯域成分を生成し、該周波数帯域成分を解析することにより決定されることを特徴とする。
【００１１】
さらに、本発明に係る空間・階調削減装置において、前記雑音閾値又は前記グラデーション閾値は、空間縮小率又は階調削減率と対応付けた対応表に基づいて決定され、前記対応表は、前記空間縮小率が高くなるほど、前記雑音閾値又は前記グラデーション閾値が小さくなるように対応付けているか、あるいは、前記階調削減率が高くなるほど、前記雑音閾値又は前記グラデーション閾値が大きくなるように対応付けていることを特徴とする。
【００１２】
さらに、本発明に係る空間・階調削減装置において、前記空間縮小部は、解像度の縮小劣化過程を模擬した関数、又はウェーブレット分解を用いて前記空間縮小画像を生成することを特徴とする。
【００１３】
さらに、本発明に係る空間・階調削減装置において、前記最適化部は、前記原画像に対する空間・階調削減画像のパワースペクトル差分値及びヒストグラム差分値を算出し、両者の合計値が最も小さい空間・階調削減画像を選択することを特徴とする。
【００１４】
また、上記課題を解決するため、本発明に係るプログラムは、コンピュータを、上記空間・階調削減装置として機能させることを特徴とする。
【発明の効果】
【００１５】
本発明によれば、効率的且つ高精度に空間・階調削減画像を得ることができるようになる。
【図面の簡単な説明】
【００１６】
【図１】本発明に係る空間・階調削減処理の概要を説明する図である。
【図２】本発明の第１の実施形態に係る空間・階調削減装置の構成例を示すブロック図である。
【図３】本発明の第１の実施形態に係る空間・階調削減装置における階調削減部の構成例を示すブロック図である。
【図４】本発明の第１の実施形態に係る空間・階調削減装置における周波数分解部の処理を説明する図である。
【図５】本発明の第１の実施形態に係る空間・階調削減装置における雑音除去画像生成部の構成例を示すブロック図である。
【図６】本発明の第１の実施形態に係る空間・階調削減装置における孤立点検出部の構成例を示すブロック図である。
【図７】本発明の第１の実施形態に係る空間・階調削減装置におけるグラデーション画像生成部の構成例を示すブロック図である。
【図８】本発明の第１の実施形態に係る空間・階調削減装置における階調削減初期情報生成部の第１の構成例を示すブロック図である。
【図９】本発明の第１の実施形態に係る空間・階調削減装置における階調削減初期情報生成部の第２の構成例を示すブロック図である。
【図１０】本発明の第１の実施形態に係る空間・階調削減装置における階調削減初期情報生成部の第３の構成例を示すブロック図である。
【図１１】本発明の第１の実施形態に係る空間・階調削減装置における空間縮小部の処理例を説明する図である。
【図１２】本発明の第１の実施形態に係る空間・階調削減装置における最適化部の処理例を説明する図である。
【図１３】本発明の第２の実施形態に係る空間・階調削減装置の構成例を示すブロック図である。
【図１４】本発明の第２の実施形態に係る空間・階調削減装置における階調削減部の構成例を示すブロック図である。
【発明を実施するための形態】
【００１７】
まず、本発明の概要について説明する。以下の説明において、原画像を階調削減処理し、且つ空間縮小処理することを「空間・階調削減処理」といい、空間・階調削減処理された画像のことを「空間・階調削減画像」という。
【００１８】
図１は、本発明に係る空間・階調削減処理の概要を説明する図である。図１に示す例では、「８Ｋ×４Ｋ解像度、１０ｂｉｔ階調」の原画像から「４Ｋ×２Ｋ解像度、８ｂｉｔ階調」の空間・階調削減画像を生成する際の処理概要を示している。本発明に係る空間・階調削減装置は、原画像をまず階調削減処理して「８Ｋ×４Ｋ解像度、８ｂｉｔ階調」の複数の階調削減画像を生成する。次に、この複数の階調削減画像を「４Ｋ×２Ｋ解像度、８ｂｉｔ階調」の画像に空間縮小し、複数の空間・階調削減画像を生成する。そして、原画像を基準として、階調削減処理及び空間縮小処理に用いたパラメータを一元的に最適化し、最適な空間・階調削減画像を選択する。
【００１９】
以下、本発明の実施形態について、詳細に説明する。
【００２０】
（第１の実施形態）
図２は、本発明の第１の実施形態に係る空間・階調削減装置の構成例を示すブロック図である。図２に示すように、空間・階調削減装置１は、階調削減部１０と、空間縮小部２０と、最適化部３０とを備える。
【００２１】
［階調削減部について］
まず、階調削減部１０について説明する。階調削減部１０は、原画像をｊ個（ｊ≧１）のパラメータを用いて階調削減処理してｊ枚の階調削減画像を生成し、空間縮小部２０に出力する。階調削減処理は、従来の階調数の単純削減やLloyd-Max法などにより行うことができる。ただし、階調数の単純削減では画像劣化が大きい。また、Lloyd-Max法では度数が高い階調値に多くの階調を割り当てて階調削減を行うため、原画像と階調削減画像との誤差を最小化するという観点からは最適な階調削減方法といえるが、視覚的かつ信号処理の観点からは最適な階調削減とはいえない。なぜなら、Lloyd-Max法では例えば雑音が支配的な階調値においても、その度数が高ければ多くの階調を割り当てて階調削減を行ってしまうからである。そこで、以下に説明するように雑音、グラデーション、又は双方を考慮して階調削減処理を行うのが好適である。
【００２２】
図３は、階調削減部１０の構成例を示すブロック図である。図３に示す例では、階調削減部１０は、周波数分解部１１と、雑音閾値決定部１２と、雑音除去画像生成部１３と、グラデーション閾値決定部１４と、グラデーション画像生成部１５と、階調削減初期情報生成部１６と、階調決定部１７とを備える。ここでは、雑音及びグラデーションの双方を考慮する場合について説明するが、雑音閾値決定部１２及び雑音除去画像生成部１３を備えないで原画像のグラデーションのみを考慮するようにしてもよく、また、グラデーション閾値決定部１４及びグラデーション画像生成部１５を備えないで原画像の雑音のみを考慮するようにしてもよい。
【００２３】
周波数分解部１１は、原画像を周波数分解（例えば、ウェーブレットパケット分解）して、複数の周波数帯域成分を生成し、雑音閾値決定部１２、雑音除去画像生成部１３、グラデーション閾値決定部１４、及びグラデーション画像生成部１５に出力する。周波数帯域成分を解析することにより、後述するように雑音閾値及びグラデーション閾値を決定する。以下では、説明の便宜上、周波数分解成分のように周波数帯域毎に分解された画像を構成する最小単位を「要素」と称し、原画像のように周波数分解されていない画像を構成する最小単位を「画素」と称して、両者を区別することとする。
【００２４】
図４は、周波数分解部１１による空間方向の周波数分解処理を説明する図である。図４では、原画像を空間方向に２階ウェーブレットパケット分解した場合の分解図を示しており、太線は１階ウェーブレットパケット分解の各周波数帯域を示している。図４に示すウェーブレットパケット分解の例では、低周波側、高周波側を空間方向に均等に周波数分解している。この分解図では空間周波数成分を、水平方向においては右側ほど高周波成分とし、垂直方向においては下側ほど高周波成分としている。よって、図４に示すように、左上は空間最低周波数帯域成分（周波数分解部１１により生成された複数の周波数分解成分のうち最も周波数が低い帯域成分）となり、右下は空間最高周波数帯域成分（周波数分解部１１により生成された複数の周波数分解成分のうち最も周波数が高い帯域成分）となる。なお、空間方向に加えて時間方向に周波数分解を行ってもよい。
【００２５】
図４では説明の便宜上、デシメーション有りでウェーブレットパケット分解を行った場合の分解図を示しているが、後述する孤立点検出部１３１にて原画像の画素単位で孤立点（雑音）を除去することができるようにするために、周波数分解部１１はデシメーション無し（すなわち、周波数帯域成分の画像サイズの縮小無し）で周波数分解を行い、生成される各周波数帯域成分の画像サイズを同一としてもよい。
【００２６】
雑音閾値決定部１２は、周波数分解部１１により生成された周波数分解成分の全帯域成分に対する空間最高周波数帯域成分のパワーの割合（以下、「空間最高周波数帯域成分の割合」という）を、各要素位置について算出する。そして、雑音閾値決定部１２は、空間最高周波数帯域成分の割合の平均値を基準として、複数の空間最高周波数帯域成分の割合を雑音閾値として決定し、雑音除去画像生成部１３に出力する。例えば、０％から、空間最高周波数帯域成分の割合の平均値までをａ等分した複数の値を雑音閾値として決定する。この場合、空間最高周波数帯域成分の割合の平均値が１０％、ａ＝５とすると、２％、４％、６％、８％、１０％が雑音閾値となる。雑音閾値決定部１２により決定された雑音閾値は、補助情報（パラメータ情報）として最適化部３０に出力される。
【００２７】
雑音除去画像生成部１３は、周波数分解部１１により生成された周波数帯域成分、及び雑音閾値決定部１２により決定された雑音閾値を用いて原画像の孤立点を検出し、該孤立点を雑音とみなす。そして、原画像から雑音が除去された雑音除去画像を生成し、階調削減初期情報生成部１６に出力する。
【００２８】
図５は、雑音除去画像生成部１３の構成例を示すブロック図である。図５に示すように、雑音除去画像生成部１３は、孤立点検出部１３１と、再構成部１３２とを備える。
【００２９】
孤立点検出部１３１は、周波数分解部１１により生成された周波数帯域成分を解析して原画像の孤立点を検出する。具体的には、孤立点検出部１３１は、空間最高周波数帯域成分の割合を各要素位置について算出する。そして、孤立点検出部１３１は、雑音閾値決定部１２により決定された複数の雑音閾値のうちの１つを選択し、算出した空間最高周波数帯域成分の割合が選択した雑音閾値Ｔｈ１を超える要素の値を１とし、雑音閾値Ｔｈ１以下の要素の値を０とする２値化画像Ｂを生成する。
【００３０】
孤立点検出部１３１は、２値化画像Ｂの要素値が１である要素について、該要素を中心とする所定の判定領域内の要素値の合計値と、所定の閾値Ｔｈ２とを比較する。２値化画像Ｂの判定領域内の要素値の合計値が閾値Ｔｈ２を超える場合には、当該要素は孤立点要素ではないと判定し、２値化画像Ｂの判定領域内の要素値の合計値が閾値Ｔｈ２以下である場合には、当該要素を孤立点要素と判定し、孤立点要素に対応する原画像の画素を孤立点と判定する。
【００３１】
図６は、孤立点検出部１３１による孤立点判定処理を説明する図である。図６に示す例では、判定領域は３×３要素である。閾値Ｔｈ２を１とすると、２値化画像Ｂの要素値が１となる要素素の周囲の８要素の要素値が０であるときのみ、２値化画像Ｂの要素値が１である要素を孤立要素とみなす。よって、閾値Ｔｈ２＝１の場合、図中の要素Ｐ１は孤立点要素であると判定され、要素Ｐ２と要素Ｐ３は孤立点要素と判定されない。
【００３２】
そして、孤立点検出部１３１は、各周波数帯域成分について、孤立点要素であると判定した要素位置の要素値を０とし、再構成部１３２に出力する。ただし、直流成分に雑音が含まれる場合には孤立点要素として検出されないおそれがあるため、直流成分を有する空間最低周波数帯域成分については階調値を変更しないようにするのが好適である。
【００３３】
再構成部１３２は、孤立点検出部１３１から入力される孤立点検出後の周波数帯域成分を用いて再構成処理を行い、雑音除去画像を生成し、階調削減初期情報生成部１６に出力する。例えば、周波数分解部１１においてｎ階ウェーブレットパケット分解を行って周波数帯域成分を生成していた場合には、孤立点検出後の周波数帯域成分をｎ階ウェーブレットパケット再構成する。
【００３４】
グラデーション閾値決定部１４は、周波数分解部１１により生成された周波数分解成分の全帯域成分に対する空間最低周波数帯域成分のパワーの割合（以下、「空間最低周波数帯域成分の割合」という）を、各要素位置について算出する。そして、グラデーション閾値決定部１４は、空間最低周波数帯域成分の割合の平均値を基準として、複数の空間最低周波数帯域成分の割合をグラデーション閾値として決定し、グラデーション画像生成部１５に出力する。例えば、空間最低周波数帯域成分の割合の平均値から１００％までをｂ等分した複数の値をグラデーション閾値として決定する。この場合、空間最低周波数帯域成分の割合平均値が９０％、ｂ＝５とすると、９０％、９２％、９４％、９６％、９８％がグラデーション閾値となる。グラデーション閾値決定部１４により決定されたグラデーション閾値は、補助情報として最適化部３０に出力される。
【００３５】
グラデーション画像生成部１５は、周波数分解部１１により生成された周波数帯域成分、及びグラデーション閾値決定部１４により決定されたグラデーション閾値を用いて原画像のグラデーション領域を検出し、該グラデーション領域のみからなるグラデーション画像を生成する。そして、グラデーション画像を階調削減初期情報生成部１６に出力する。
【００３６】
図７は、グラデーション画像生成部１５の構成例を示すブロック図である。図７に示すように、グラデーション画像生成部１５は、グラデーション領域検出部１５１と、再構成部１５２とを備える。
【００３７】
グラデーション領域検出部１５１は、周波数分解部１１により生成された周波数帯域成分を解析して原画像のグラデーション領域を検出する。具体的には、グラデーション領域検出部１５１は、空間最低周波数帯域成分の割合を各要素位置について算出する。そして、グラデーション領域検出部１５１は、グラデーション閾値決定部１４により決定された複数のグラデーション閾値のうちの１つを選択し、算出した空間最低周波数帯域成分の割合が選択したグラデーション閾値を超える領域をグラデーション領域と決定する。そして、グラデーション領域のみからなる画像を、再構成部１５２に出力する。
【００３８】
再構成部１５２は、グラデーション領域検出部１５１から入力されるグラデーション領域検出後の周波数帯域成分を用いて再構成処理を行い、グラデーション画像を生成し、階調削減初期情報生成部１６に出力する。例えば、周波数分解部１１においてｎ階ウェーブレットパケット分解を行って周波数帯域成分を生成していた場合には、グラデーション領域検出後の周波数帯域成分をｎ階ウェーブレット再構成する。
【００３９】
階調削減初期情報生成部１６は、雑音除去画像生成部１３により生成された雑音除去画像、及びグラデーション画像生成部１５により生成されたグラデーション画像に基づいて訓練データを生成するとともに、階調削減率（階調削減ビット数）に基づく階調変換テーブルを生成する。そして、訓練データ及び階調変換テーブルを階調決定部１７に出力する。
【００４０】
図８は、階調削減初期情報生成部１６の構成例を示すブロック図である。図８に示すように、階調削減初期情報生成部１６は、ヒストグラム生成部１６１と、重み付け部１６２と、階調変換テーブル生成部１６３とを備える。
【００４１】
ヒストグラム生成部１６１は、雑音除去画像生成部１３により生成された雑音除去画像について階調値ごとの度数を示すヒストグラムを生成し、重み付け部１６２に出力する。
【００４２】
重み付け部１６２は、グラデーション画像生成部１５により生成されたグラデーション画像の階調値を検出し、ヒストグラム生成部１６１により生成された雑音除去画像のヒストグラムに対して、グラデーション画像の階調値の度数が高くなるように重み付けする。そして、重み付けしたヒストグラムを訓練データとして階調決定部１７に出力する。例えば、グラデーション画像の階調値の度数を、所定の１を超える値（例えば、１．２）を乗じた値に変更する。なお、度数の高いグラデーション画像の階調値ほど重みを大きくするようにしてもよい。重み付けしたヒストグラムを訓練データとすることにより、階調決定部１７にて階調削減を行う際に、グラデーション領域により多くの階調を割り当てることができる。
【００４３】
階調変換テーブル生成部１６３は、階調削減率に応じて階調変換テーブルを生成し、階調決定部１７に出力する。階調変換テーブルは、原画像の色深度がｎビット（階調数２ｎ）、階調削減ビット数がｍビットの場合、ｎビットを（ｎ−ｍ）ビットに線形変換する階調変換テーブルである。
【００４４】
図９は、階調削減部１０がグラデーション閾値決定部１４及びグラデーション画像生成部１５を備えない場合の階調削減初期情報生成部１６の構成例を示すブロック図である。この場合、ヒストグラム生成部１６１は、雑音除去画像生成部１３により生成された雑音除去画像について階調値ごとの度数を示すヒストグラムを訓練データとして生成する。階調変換テーブル生成部１６３は、階調削減率に応じて階調変換テーブルを生成する。
【００４５】
図１０は、階調削減部１０が雑音閾値決定部１２及び雑音除去画像生成部１３を備えない場合の階調削減初期情報生成部１６の構成例を示すブロック図である。この場合、ヒストグラム生成部１６１は、原画像について階調値ごとの度数を示すヒストグラムを生成し、重み付け部１６２に出力する。重み付け部１６２は、グラデーション画像生成部１５により生成されたグラデーション画像の階調値を検出し、ヒストグラム生成部１６１により生成された原画像のヒストグラムに対して、グラデーション画像の度数が高くなるように重み付けし、重み付けしたヒストグラムを訓練データとする。階調変換テーブル生成部１６３は、階調削減率に応じて階調変換テーブルを生成する。
【００４６】
階調決定部１７は、階調削減初期情報生成部１６により生成された訓練データ及び階調変換テーブルを初期値として用いるLloyd-Max法により、原画像の階調を削減した階調削減画像を生成する。階調決定部１７は、訓練データの度数に応じて各階調の量子化ステップを決定して階調削減画像を生成した後、階調削減画像を逆階調変換して原画像と同じ階調数の画像を生成する。そして、原画像と逆階調変換した画像との差分値が所定の閾値以下となるまで、訓練データ、及び階調変換テーブルを更新し、階調削減処理を繰り返し行う。Lloyd-Max法によれば、度数の高い階調ほど量子化ステップを小さくするため、原画像と階調削減画像との誤差を小さくすることができる。
【００４７】
また、階調決定部１７は、生成した階調削減部を元の階調数に復元するための逆量子化テーブルを生成し、補助情報として最適化部３０に出力する。ここで、逆量子化テーブルとは、階調削減画像を生成した際の階調変換テーブルの入力値を出力値とし、階調変換テーブルの出力値を入力値としたテーブルである。
【００４８】
雑音除去画像生成部１３は、雑音閾値決定部１２により決定された雑音閾値ごとに雑音除去画像を生成する。また、グラデーション画像生成部１５は、グラデーション閾値決定部１４により決定されたグラデーション閾値ごとにグラデーション画像を生成する。よって、雑音閾値決定部１２により決定された雑音閾値がａ個、グラデーション閾値決定部１４により決定されたグラデーション閾値がｂ個の場合には、階調削減初期情報生成部１６によりａ×ｂ種類の訓練データが生成され、階調決定部１７によりｊ＝ａ×ｂ枚の階調削減画像（階調削減画像群）が生成される。
【００４９】
［空間縮小部について］
次に、空間縮小部２０について説明する。空間縮小部２０は、階調削減部１０により生成されたｊ枚の階調削減画像をそれぞれｋ個（ｋ≧１）のパラメータを用いて空間縮小処理してｊ×ｋ枚（但し、ｊ×ｋ≧２）の空間・階調削減画像を生成し、最適化部３０に出力する。
【００５０】
図１１は、空間縮小部２０の構成例を示すブロック図である。図１１に示す例では、空間縮小部２０は、平滑化部２１と、縮小部２２とを備える。
【００５１】
平滑化部２１は、階調削減画像群をそれぞれ平滑化処理して平滑化画像群を生成し、縮小部２２に出力する。平滑化は、解像度の縮小劣化過程（ぼやけ）を模擬した関数を用いて階調拡大画像を畳み込むことで処理する。解像度の縮小劣化過程を模擬した関数として、例えば、点拡がり関数（ＰＳＦ：Point spread function）又はガウシアン関数を用いることができる。平滑化部２１は、タップ数、分散値、ゲイン情報などを変更したｋ個の関数（パラメータ）を用いて平滑化処理を行う。また、平滑化部２１は、用いた関数のタップ数、分散値、ゲイン情報を、補助情報として最適化部３０に出力する。
【００５２】
縮小部２２は、平滑化部２１により生成された平滑化画像群をそれぞれ目標解像度になるように空間縮小処理し、空間・階調削減画像群を生成し、最適化部３０に出力する。空間縮小処理は、単純間引きにより行ってもよいし、平滑化画像を目標解像度になるようウェーブレット分解し、その空間低周波数成分を空間・階調削減画像としてもよい。ウェーブレット分解を行った場合には、そのフィルタの情報も補助情報として最適化部３０に出力される。
【００５３】
［最適化部について］
次に、最適化部３０について説明する。最適化部３０は、空間縮小部２０により生成されたｊ×ｋ枚（但し、ｊ×ｋ≧２）の空間・階調削減画像から、原画像と比較して最も劣化の少ない（最も差分の少ない）画像を選択して出力する。例えば、原画像とパワースペクトルの概形、ヒストグラムの概形、又はその両者の概形が最も近い画像を最も劣化の少ない画像とすることができる。また、最適化部３０は、階調削減部１０及び空間縮小部２０により生成された補助情報を入力し、選択した空間・階調削減画像に対応する補助情報を、最適補助情報として出力する。最適補助情報を出力することにより、該最適補助情報及び空間・階調削減画像を用いて、元の原画像を復元することができるようになる。
【００５４】
図１２は、最適化部３０の処理の一例を説明する図である。図１２に示す例では、最適化部３０は、原画像に対する各空間・階調削減画像のパワースペクトル差分値及びヒストグラム差分値を算出し、両者の合計値が最も小さい空間・階調削減画像を最も劣化の少ない画像として選択する。
【００５５】
具体的には、最適化部３０は、パワースペクトル差分値を算出するために、原画像及び空間・階調削減画像群をそれぞれｎ階ウェーブレット分解する。そして、原画像及び空間・階調削減画像の組み合わせごとに、各ｎ階帯域のＲＭＳ（Root Mean Square）パワー値の差分合計値のα倍の値を算出する。なお、図１２では空間・階調削減画像の水平方向の解像度が原画像の１／２倍の場合を示している。
【００５６】
また、最適化部３０は、ヒストグラム差分値を算出するために、原画像及び空間・階調削減画像群のヒストグラムを求め、各ヒストグラムのレベルをｍレベルに量子化する。そして、原画像及び空間・階調削減画像の組み合わせごとに、各ｍレベルの頻度の差合計値のβ倍の値を計算する。なお、図１２では空間・階調削減画像の階調数が原画像の１／４倍の場合を示している。ここで、α及びβは外部から設定可能であり、空間縮小処理の精度を重視する場合はαに大きな重みを与え、階調削減処理の精度を重視する場合はβに大きな重みを与えることができる。
【００５７】
最後に、最適化部３０は、パワースペクトル差分値及びヒストグラム差分値の和をパワースペクトル・ヒストグラム差分値として算出し、この値が最小となる空間・階調削減画像を選択する。
【００５８】
以上説明したように、空間・階調削減装置１は、階調削減部１０により原画像をｊ個（ｊ≧１）のパラメータを用いて階調削減処理してｊ枚の階調削減画像を生成し、空間縮小部２０によりｊ枚の階調削減画像をｋ個（ｋ≧１）のパラメータを用いて空間縮小処理してｊ×ｋ枚（但し、ｊ×ｋ≧２）の空間・階調削減画像を生成する。そして、最適化部３０によりｊ×ｋ枚の空間・階調削減画像を前記原画像と比較して、最も劣化の少ない空間・階調削減画像を選択する。このため、空間・階調削減装置１によれば、階調削減処理と空間縮小処理とを個別に行う場合よりも、効率的且つ高精度に空間・階調削減画像を得ることができる。
【００５９】
また、前記階調削減部１０は、ｊ個の雑音閾値を用いて、原画像から孤立点を除去したｊ枚の雑音除去画像を生成するか、又はｊ個のグラデーション閾値を用いて、原画像のグラデーション領域のみからなるｊ枚のグラデーション画像を生成する画像生成部（雑音除去画像生成部１３又はグラデーション画像生成部１５）と、雑音除去画像又は原画像のヒストグラムに対して、グラデーション画像の階調値の度数が高くなるように重み付けしたヒストグラムを訓練データとして用いるLloyd-Max法により、原画像の階調を削減したｊ枚の階調削減画像を生成する階調決定部１７と、を備えるのが好適である。雑音又はグラデーションを考慮して階調変換することにより、視覚的かつ信号処理的に優れ、画像ごとに最適化された階調削減を行うことができる。また、空間縮小前に階調削減を行うことで、不要な雑音成分が折り返し成分として空間削減画像に含まれるのを防止することができる。
【００６０】
なお、雑音及びグラデーションの双方を考慮する場合には、雑音除去画像生成部１３によりａ個の雑音閾値を用いて、原画像から孤立点を除去したａ枚の雑音除去画像を生成し、ｂ個のグラデーション閾値を用いて、原画像のグラデーション領域のみからなるｂ枚のグラデーション画像を生成し、階調決定部１７によりｊ＝ａ×ｂ枚の階調削減画像を生成する。
【００６１】
（第２の実施形態）
次に、本発明による第２の実施形態について説明する。図１３は、本発明の第２の実施形態に係る空間・階調削減装置の構成を示すブロック図である。図１３に示すように、空間・階調削減装置２は、階調削減閾値決定部４０と、階調削減部１００と、空間縮小部２０と、最適化部３０とを備える。第２の実施形態の空間・階調削減装置２は、第１の実施形態の空間・階調削減装置１と比較して、階調削減閾値決定部４０を更に設けて、第１の実施形態の空間・階調削減装置１と異なる方法で階調削減時の閾値を決定する。空間縮小部２０及び最適化部３０については第１の実施形態と同一であるため、説明を省略する。
【００６２】
階調削減閾値決定部４０は、空間縮小率及び／又は階調削減率に基づいて、雑音閾値及び／又はグラデーション閾値を決定し、階調削減部１００に出力する。ここでは、空間縮小率及び階調削減率に基づいて、雑音閾値及びグラデーション閾値を決定する場合について説明する。
【００６３】
階調削減閾値決定部４０は、空間縮小率及び階調削減率と、雑音閾値及びグラデーション閾値とを対応付けた対応表を予め有する。対応表の例を表１に示す。空間縮小率と、雑音閾値及びグラデーション閾値との関係については、空間縮小率が高くなるほど、雑音閾値及びグラデーション閾値を小さくする。空間縮小率が高い場合は、雑音成分の折り返しによる影響は少なく、またグラデーション領域検出を高い精度で行う必要はないためである。また、階調削減率と、雑音閾値及びグラデーション閾値との関係については、階調削減率が高くなるほど、雑音閾値及びグラデーション閾値を大きくするのが好適である。階調削減率が高い場合は、雑音に階調を割り当てないように雑音除去閾値を高く設定する必要があり、また疑似輪郭防止のためにグラデーション領域検出閾値を高く設定する必要があるからである。つまり、表１では、Ｎ１＜Ｎ２＜Ｎ３＜Ｎ４、Ｇ１＜Ｇ２＜Ｇ３＜Ｇ４とする。
【００６４】
【表１】
2014238657.tif 000003
【００６５】
表１の対応表を用いて雑音閾値及びグラデーション閾値を決定する場合、空間縮小処理に対する重みα、及び階調削減処理に対する重みβを用いるのが好適である。例えば、空間縮小率が１／２、階調削減率が１/８の場合に表１を用いて雑音閾値Ｎ及びグラデーション閾値Ｇを求めると、Ｎ＝（αＮ４＋βＮ３）／２、Ｇ＝（αＧ４＋βＧ３）／２となる。空間縮小処理の精度を重視する場合はαに大きな重みを与え、階調削減処理の精度を重視する場合はβに大きな重みを与える。なお、雑音閾値及びグラデーション閾値の組を複数組選択するようにしてもよい。
【００６６】
図１４は、階調削減部１００の構成例を示すブロック図である。図１４に示す例では、階調削減部１０は、周波数分解部１１と、雑音除去画像生成部１３と、グラデーション画像生成部１５と、階調削減初期情報生成部１６と、階調決定部１７とを備える。各構成部の詳細な処理内容は第１の実施形態で説明した通りであるため、ここでは概略のみ述べる。
【００６７】
雑音除去画像生成部１３は、周波数分解部１１により生成された周波数帯域成分、及び階調削減閾値決定部４０により決定された雑音閾値を用いて原画像の孤立点を検出し、該孤立点を雑音とみなす。そして、原画像から雑音が除去された雑音除去画像を生成し、階調削減初期情報生成部１６に出力する。
【００６８】
グラデーション画像生成部１５は、周波数分解部１１により生成された周波数帯域成分、及び階調削減閾値決定部４０により決定されたグラデーション閾値を用いて原画像のグラデーション領域を検出し、該グラデーション領域のみからなるグラデーション画像を生成する。そして、グラデーション画像を階調削減初期情報生成部１６に出力する。
【００６９】
階調削減初期情報生成部１６は、雑音除去画像生成部１３により生成された雑音除去画像、及びグラデーション画像生成部１５により生成されたグラデーション画像に基づいて訓練データを生成するとともに、階調削減率に基づく階調変換テーブルを生成する。そして、訓練データ及び階調変換テーブルを階調決定部１７に出力する。
【００７０】
階調決定部１７は、階調削減初期情報生成部１６により生成された訓練データ及び階調変換テーブルを初期値として用いるLloyd-Max法により、原画像の階調を削減した階調削減画像を生成する。また、階調決定部１７は逆量子化テーブルを生成し、補助情報として最適化部３０に出力する。
【００７１】
以上説明したように、空間・階調削減装置２は、階調削減閾値決定部４０を更に備え、雑音閾値又は前記グラデーション閾値を、該閾値と、空間縮小率又は階調削減率とを対応付けた対応表に基づいて決定する。このため、空間・階調削減装置２によれば、第１の実施形態の空間・階調削減装置１と同様の効果に加え、雑音閾値又は前記グラデーション閾値を迅速に決定して処理時間を短縮することができるようになる。
【００７２】
なお、上述した空間・階調削減装置１，２として機能させるためにコンピュータを用いることができ、そのようなコンピュータは、空間・階調削減装置１，２の各機能を実現する処理内容を記述したプログラムを当該コンピュータの記憶部に格納しておき、当該コンピュータのＣＰＵによってこのプログラムを読み出して実行させることで実現することができる。なお、このプログラムは、コンピュータ読取り可能な記録媒体に記録することができる。
【００７３】
上述の実施形態は、代表的な例として説明したが、本発明の趣旨及び範囲内で、多くの変更及び置換ができることは当業者に明らかである。したがって、本発明は上述の実施形態によって制限するものと解するべきではなく、特許請求の範囲から逸脱することなく、種々の変形や変更が可能である。例えば、実施形態に記載の複数の構成ブロックを１つに組み合わせたり、あるいは１つの構成ブロックを分割したりすることが可能である。
【産業上の利用可能性】
【００７４】
このように、本発明は原画像の階調削減処理及び空間縮小処理を行って空間・階調削減画像を生成する任意の用途に有用である。
【符号の説明】
【００７５】
１，２  空間・階調削減装置
１０    階調削減部
１１    周波数分解部
１２    雑音閾値決定部
１３    雑音除去画像生成部
１４    グラデーション閾値決定部
１５    グラデーション画像生成部
１６    階調削減初期情報生成部
１７    階調決定部
２０    空間縮小部
２１    平滑化部
２２    縮小部
３０    最適化部
４０    階調削減閾値決定部
１００  階調削減部
１３１  孤立点検出部
１３２  再構成部
１５１  グラデーション領域検出部
１５２  再構成部
１６１  ヒストグラム生成部
１６２  重み付け部
１６３  階調変換テーブル生成部
【図１】
2014238657.tif 000004
【図２】
2014238657.tif 000005
【図３】
2014238657.tif 000006
【図４】
2014238657.tif 000007
【図５】
2014238657.tif 000008
【図６】
2014238657.tif 000009
【図７】
2014238657.tif 000010
【図８】
2014238657.tif 000011
【図９】
2014238657.tif 000012
【図１０】
2014238657.tif 000013
【図１１】
2014238657.tif 000014
【図１２】
2014238657.tif 000015
【図１３】
2014238657.tif 000016
【図１４】
2014238657.tif 000017