(57)【要約】  【課題】  ウェーブレット縮退を利用したノイズ処理方式であり、エッジ部分のボケを抑制しつつノイズ低減処理を行い、鮮明な画像が得られる画像処理装置を提供する。  【解決手段】  画像処理装置１は、画像に対してウェーブレット変換を行うウェーブレット変換部１０３と、ウェーブレット変換によって得られる信号成分のうち低域成分に基づいて、信号成分をクラスタリングするクラスタ解析処理部１０４と、クラスタ解析によって得られたクラスタごとに信号成分の統計的な特徴量を算出する特徴量算出部１０６と、特徴量に応じてコアリング閾値を決定するコアリング閾値決定部１０７と、コアリング閾値決定部１０７にて決定されたコアリング閾値を用いてコアリング処理を行うコアリング処理部１０８と、コアリング処理を行った信号成分に対して逆ウェーブレット変換を行う逆ウェーブレット変換部１０９とを備えた構成を有する。【選択図】  図１
2009098742.tif 000002
【特許請求の範囲】
【請求項１】
画像に対してウェーブレット変換を行うウェーブレット変換部と、
ウェーブレット変換によって得られる信号成分のうち低域成分に基づいて、前記信号成分をクラスタリングするクラスタ解析処理部と、
クラスタ解析によって得られたクラスタごとに前記信号成分の統計的な特徴量を算出する特徴量算出部と、
前記特徴量に応じてコアリング閾値を決定するコアリング閾値決定部と、
前記コアリング閾値決定部にて決定されたコアリング閾値を用いてコアリング処理を行うコアリング処理部と、
前記コアリング処理を行った信号成分に対して逆ウェーブレット変換を行う逆ウェーブレット変換部と、
を備える画像処理装置。
【請求項２】
前記特徴量算出部は、前記ウェーブレット変換によって得られた前記信号成分のうち高域成分に基づいて、前記特徴量を算出する請求項１に記載の画像処理装置。
【請求項３】
前記特徴量算出部は、前記信号成分の分散値を前記特徴量として算出し、
前記コアリング閾値決定部は、前記分散値が大きい場合に前記分散値が小さい場合より大きいコアリング閾値を決定する請求項１または２に記載の画像処理装置。
【請求項４】
前記コアリング閾値決定部は、ウェーブレット変換の単位となる信号成分が異なるクラスタに分類される場合に、前記単位に含まれる信号成分に対し、前記信号成分が同じクラスタに分類される場合より低いコアリング閾値を決定する請求項１〜３のいずれかに記載の画像処理装置。
【請求項５】
前記ウェーブレット変換部は、複数段階のウェーブレット変換を行い、
前記クラスタ解析処理部は、各段階において算出された信号成分をクラスタリングし、
前記コアリング閾値決定部は、各段階におけるコアリング閾値を決定する請求項１〜４のいずれかに記載の画像処理装置。
【請求項６】
請求項１〜５のいずれかに記載の画像処理装置を有する携帯電話。
【請求項７】
画像に対してウェーブレット変換を行うステップと、
ウェーブレット変換によって得られる信号成分のうち低域成分に基づいて、前記信号成分をクラスタリングするステップと、
クラスタ解析によって得られたクラスタごとに前記信号成分の統計的な特徴量を算出するステップと、
前記特徴量に応じてコアリング閾値を決定するステップと、
前記コアリング閾値を用いてコアリング処理を行うステップと、
前記コアリング処理を行った信号成分に対して逆ウェーブレット変換を行うステップと、
を備える画像処理方法。
【請求項８】
画像のノイズを低減する処理を行うためのプログラムであって、コンピュータに、
画像に対してウェーブレット変換を行うステップと、
ウェーブレット変換によって得られる信号成分のうち低域成分に基づいて、前記信号成分をクラスタリングするステップと、
クラスタ解析によって得られたクラスタごとに前記信号成分の統計的な特徴量を算出するステップと、
前記特徴量に応じてコアリング閾値を決定するステップと、
前記コアリング閾値を用いてコアリング処理を行うステップと、
前記コアリング処理を行った信号成分に対して逆ウェーブレット変換を行うステップと、
を備えるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ウェーブレット変換を用いてノイズを低減する画像処理装置に関する。
【背景技術】
【０００２】
従来、ノイズを低減する画像処理方法としては、ローパスフィルタと呼ばれる平滑化フィルタによってノイズ（高周波成分）を抑制する方法が知られている。ローパスフィルタは、周辺画素との加算平均をとることによって平滑化を行うフィルタである。
【０００３】
また、別の方法として、多重解像度変換の一つであるウェーブレット変換を利用するノイズ低減方法が知られている。この方法では、画像信号を、ウェーブレット変換によって低域周波数帯と高域周波数帯とのサブバンドに分離し、得られた周波数信号に対して、同様の処理を繰り返すことで、周波数帯域毎のサブバンドに分離された周波数帯域信号に分離する。ウェーブレット変換によって得られた高域成分のうち、値の小さいウェーブレット変換係数はノイズ成分であるという前提の下で、所定の閾値よりも小さい高域成分を０に置き換えるコアリング処理を行う。その後、コアリング処理された高域成分と低域成分にウェーブレット逆変換を行って信号を再構成し、ノイズが低減された画像を得る。このノイズ低減方法は、１９９４年Ｄｏｎｏｈｏらによって提案されたウェーブレット縮退という手法を用いたものとして非特許文献１に紹介されている。
【非特許文献１】中野宏毅、山本鎮男、吉田靖夫「ウェーブレットによる信号処理と画像処理」共立出版  P.101
【発明の開示】
【発明が解決しようとする課題】
【０００４】
上記した平滑化フィルタを用いるノイズ低減方法では、一の画像に対して一律にノイズ低減の強さが決定されることになる。従って、この方法は、画素同士の相関が高い場合には効果が大きいが、画素同士の相関が高くない場合には、エッジのボケなどが発生し、輪郭部分が不鮮明となる場合があった。
【０００５】
本発明は、上記課題を解決されるためになされたもので、ウェーブレット縮退を利用したノイズ処理方式であり、エッジ部分のボケを抑制しつつノイズ低減処理を行い、鮮明な画像が得られる画像処理装置を提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明の画像処理装置は、画像に対してウェーブレット変換を行うウェーブレット変換部と、ウェーブレット変換によって得られる信号成分のうち低域成分に基づいて、前記信号成分をクラスタリングするクラスタ解析処理部と、クラスタ解析によって得られたクラスタごとに前記信号成分の統計的な特徴量を算出する特徴量算出部と、前記特徴量に応じてコアリング閾値を決定するコアリング閾値決定部と、前記コアリング閾値決定部にて決定されたコアリング閾値を用いてコアリング処理を行うコアリング処理部と、前記コアリング処理を行った信号成分に対して逆ウェーブレット変換を行う逆ウェーブレット変換部とを備えた構成を有する。
【０００７】
このようにウェーブレット変換によって得られた信号成分のうちの低域成分を用いてクラスタ解析を行うことにより、細かな特徴の違いではなく、比較的大きな性質の違いにより信号成分を分類するので、急峻なエッジの有無に応じて信号成分をクラスタリングできる。そして、クラスタごとにコアリング閾値を決定することにより、輪郭のボケを抑制しつつ、適切にノイズを低減できる。
【０００８】
本発明の画像処理装置において、前記特徴量算出部は、前記ウェーブレット変換によって得られた前記信号成分のうち高域成分に基づいて、前記特徴量を算出する構成を有する。
【０００９】
この構成により、高域成分の統計的な特徴量を計算することにより、信号成分に含まれるノイズの大きさを算出できるので、ノイズの大きさに応じてコアリング閾値を決定することができる。
【００１０】
本発明の画像処理装置において、前記特徴量算出部は、前記信号成分の分散値を前記特徴量として算出し、前記コアリング閾値決定部は、前記分散値が大きい場合に前記分散値が小さい場合より大きいコアリング閾値を決定する構成を有する。
【００１１】
信号成分の分散値が大きければノイズ成分が大きく、分散値が小さければノイズ成分が小さいと推定されるので、信号成分の分散値の大きさに応じてコアリング閾値を決定することにより、適切にノイズを低減することができる。
【００１２】
本発明の画像処理装置において、前記コアリング閾値決定部は、ウェーブレット変換の単位となる矩形領域における信号成分が異なるクラスタに分類される場合に、前記単位に含まれる信号成分に対し、前記信号成分が同じクラスタに分類される場合より低いコアリング閾値を決定する構成を有する。
【００１３】
ウェーブレット変換の単位となる信号成分が複数のクラスタに分類される場合には、当該単位にエッジが含まれると推定される。本発明の構成により、ウェーブレット変換の単位となる信号が異なるクラスタに分類される場合には、コアリング閾値を低くすることにより、エッジを保存して画像の鮮鋭さを保つことができる。
【００１４】
本発明の画像処理装置において、前記ウェーブレット変換部は、複数段階のウェーブレット変換を行い、前記クラスタ解析処理部は、各段階において算出された信号成分をクラスタリングし、前記コアリング閾値決定部は、各段階におけるコアリング閾値を決定する構成を有する。
【００１５】
このようにウェーブレット変換の各段階において適切なコアリング閾値を決定することにより、画像の鮮鋭さを保ちつつ、ノイズを低減することができる。
【００１６】
本発明の携帯電話は、上記の画像処理装置を有する。
【００１７】
この構成により、携帯電話に搭載されたカメラにて撮像した画像に対し、輪郭のボケを抑制しつつ適切にノイズを低減する画像処理を行える。
【００１８】
本発明の画像処理方法は、画像に対してウェーブレット変換を行うステップと、ウェーブレット変換によって得られる信号成分のうち低域成分に基づいて、前記信号成分をクラスタリングするステップと、クラスタ解析によって得られたクラスタごとに前記信号成分の統計的な特徴量を算出するステップと、前記特徴量に応じてコアリング閾値を決定するステップと、前記コアリング閾値を用いてコアリング処理を行うステップと、前記コアリング処理を行った信号成分に対して逆ウェーブレット変換を行うステップとを備えた構成を有する。
【００１９】
この構成により、上記した画像処理装置と同様に、急峻なエッジの有無に応じて信号成分をクラスタリングし、クラスタごとにコアリング閾値を決定するので、輪郭のボケを抑制しつつ、適切にノイズを低減できる。なお、本発明の画像処理装置の各種の構成を本発明の画像処理方法に適用することも可能である。
【００２０】
本発明のプログラムは、画像のノイズを低減する処理を行うためのプログラムであって、コンピュータに、画像に対してウェーブレット変換を行うステップと、ウェーブレット変換によって得られる信号成分のうち低域成分に基づいて、前記信号成分をクラスタリングするステップと、クラスタ解析によって得られたクラスタごとに前記信号成分の統計的な特徴量を算出するステップと、前記特徴量に応じてコアリング閾値を決定するステップと、前記コアリング閾値を用いてコアリング処理を行うステップと、前記コアリング処理を行った信号成分に対して逆ウェーブレット変換を行うステップとを備えた構成を有する。
【００２１】
この構成により、上記した画像処理装置と同様に、急峻なエッジの有無に応じて信号成分をクラスタリングし、クラスタごとにコアリング閾値を決定するので、輪郭のボケを抑制しつつ、適切にノイズを低減できる。なお、本発明の画像処理装置の各種の構成を本発明のプログラムに適用することも可能である。
【発明の効果】
【００２２】
本発明によれば、ウェーブレット変換によって得られた信号成分のうちの低域成分を用いてクラスタ解析を行うことにより、細かな特徴の違いではなく、比較的大きな性質の違いにより信号成分を分類するので、急峻なエッジの有無に応じて信号成分をクラスタリングでき、クラスタごとにコアリング閾値を決定することにより、輪郭のボケを抑制しつつ、適切にノイズを低減できるというすぐれた効果を有する。
【発明を実施するための最良の形態】
【００２３】
以下、本発明の実施の形態に係る画像処理装置について、図面を参照しながら説明する。
図１は、実施の形態の画像処理装置１の構成を示すブロック図である。図１に示すように、画像処理装置１は、画像を撮影する撮影素子１０１と、撮影素子１０１から得られた画像を格納する画像メモリ１０２と、画像メモリ１０２から読み出された画像に対し、ウェーブレット変換を行うウェーブレット変換部１０３と、ウェーブレット変換部１０３にて生成された信号成分のうちの高域成分に対してコアリング処理を行うコアリング処理部１０８と、コアリング処理された高域成分とウェーブレット変換により帯域分割され生成された低域成分を基に信号を再構成するウェーブレット逆変換部１０９と、ウェーブレット逆変換部１０９により再構成された画像を表示させる画像表示部１１０とを備えている。
【００２４】
また、本実施の形態の画像処理装置１は、コアリング処理部１０８にてコアリングを行う際に用いる閾値を決定するための構成として、クラスタ解析処理部１０４と、ラベリング部１０５と、特徴量算出部１０６と、コアリング閾値決定部１０７とを備えている。
【００２５】
以下、実施の形態の画像処理装置１の各構成について順に説明する。
図２は、ウェーブレット変換部１０３の構成を示す図である。図２では、レベル１のウェーブレット変換の例を示す。なお、レベルとはサブバンド分解した回数である。ウェーブレット変換部１０３は、多重解像度変換の１つであるウェーブレット変換を利用し、観測信号をローパスフィルタ（以下、「ＬＰＦ」という）２０１とハイパスフィルタ（以下、「ＨＰＦ」という）２０２によって、低域成分、高域成分に分解（サブバンド分解）する。ウェーブレット変換部１０３は、サブバンド分解された低域成分、高域成分を二次元の画像信号に置き換え、さらにＬＰＦ２０３，２０５、ＨＰＦ２０４，２０６によってサブバンド分解することにより、水平、垂直、斜めの高域信号成分と低域信号成分を得る。
【００２６】
図３は、二次元の水平・垂直２画素［Ａ，Ｂ，Ｃ，Ｄ］を例として、各高域成分と、低域成分の具体的な計算方法を示す図である。図３に示すように、レベル１の低域成分は、ＬＬ＝Ａ＋Ｂ＋Ｃ＋Ｄによって計算することができる。各高域成分は、ＨＬ＝Ａ−Ｂ＋Ｃ−Ｄ、ＬＨ＝Ａ＋Ｂ−Ｃ−Ｄ、ＨＨ＝Ａ−Ｂ−Ｃ＋Ｄによって計算することができる。
【００２７】
ウェーブレット変換部１０３は、上記処理によって得られた低域信号成分を、上記と同様に順次サブバンド分解することにより、上位レベルのウェーブレット変換を行うことができる。図４は、ウェーブレット変換により、レベル２までサブバンド分解した例を示す図である。上位レベルの各成分も、下位レベルの低域成分ＬＬと同様にサブバンド分解することが可能である。
【００２８】
ウェーブレット変換の具体的手法として、帯域分離の精度は低いがフィルタ係数が「１」と「−１」との２つのみで構成されるＨａａｒの基底と呼ばれる方法が挙げられる。
【００２９】
［クラスタ解析処理部］
次に、クラスタ解析処理部１０４について説明する。クラスタ解析処理部１０４は、ウェーブレット変換後の信号の性質の類似度に基づいて、クラスタに分類する機能を有する。具体的には、クラスタ解析処理部１０４は、ウェーブレット変換部１０３にてサブバンド分解された低域成分に対して、クラスタ解析を行う。クラスタ解析の際、原画像を用いるのではなく低域成分を用いることにより、細かな特徴の違いではなく、比較的大きな性質の違いにより信号成分の分類を行うことができる。これにより、急峻なエッジを検出し、エッジの有無に応じて信号成分を分類することが可能となる。なお、クラスタ解析の手法は、一般的に知られているＫ−ｍｅａｎｓ法などを用いてもよいし、本発明に適した独自のアルゴリズムを適用してもよく、手法については特に限定しない。
【００３０】
図５は、クラスタ解析の例を示す図である。図５に示す例では、レベル１までのウェーブレット変換を行い、低域成分係数が二つのクラスタα、βに分類されている。また、分類された各クラスタに所属する係数の数を、クラスタのメンバー数と呼ぶことにする。図５に示す例では、クラスタαのメンバー数は２、クラスタβのメンバー数は２である。クラスタ解析処理部１０４は、各レベルの低域成分についてクラスタ解析を行う。例えばレベル３までのウェーブレット変換を考えた場合、レベル１、レベル２、レベル３のそれぞれで得られる低域成分に対し、クラスタ解析を行う。このように、ウェーブレット変換後の低域成分を解析することにより、帯域別に詳細な画像解析を行うことが可能となる。クラスタ解析処理部１０４は、各低域成分がどのクラスタに分類されたかを示すクラスタ解析の結果を、記録部１１１に保存する。
【００３１】
［ラベリング部］
ラベリング部１０５は、低域成分のクラスタ解析結果に基づいて高域成分にラベリングを行う。具体的には、ラベリング部１０５は、同じ信号成分について求められた低域成分のクラスタ解析結果を高域成分にラベリングする。
【００３２】
図６は、４×４画素の原画像に対してレベル１のウェーブレット変換を行った例を示す図である。この低域成分に対してクラスタ解析を行った結果、図示のように低域成分係数６０１がα、βの２つのクラスタに分類された。ラベリング部１０５は、高域成分係数６０２に対して、低域成分係数６０１と同じラベリングを行う。例えば、低域成分ＬＬ００は画素Ｓ００，Ｓ０１，Ｓ１０，Ｓ１１に基づいて生成された信号成分であるが、高域成分ＨＬ００，ＬＨ００，ＨＨ００も同じ画素に基づいて生成されている。従って、高域成分ＨＬ００，ＬＨ００，ＨＨ００に、低域成分ＬＬ００と同じクラスタαをラベリングする。
【００３３】
ラベリング部１０５は、低域成分係数とそれが所属するクラスタの関係を、上記した記録部１１１から読み出す。以下、全ての高域成分に対して、低域成分と同様のラベリングを行い、図示のようなラベリング結果を得る。ラベリング部１０５は、ラベリングの結果を記録部１１１に保存する。
【００３４】
［特徴量算出部］
本発明における特徴量算出部１０６について説明する。特徴量算出部１０６は、ラベリング部１０５にてラベリングされた信号成分毎に特徴量を算出する機能を有する。特徴量算出部１０６は、高域成分に基づいて特徴量を算出する。ここで、レベル１のウェーブレット変換により得られた垂直高域成分、水平高域成分、斜め高域成分に対し、図６に示すようなラベリングがなされた場合を例として、特徴量算出部１０６の処理について説明する。特徴量算出部１０６は、垂直高域成分、水平高域成分、斜め高域成分のそれぞれについて分散値を算出し、その平均値を求めてその特徴量とする。垂直高域成分ＨＬでクラスタαがラベリングされた信号成分、水平高域成分ＬＨでクラスタαがラベリングされた信号成分、斜め高域成分ＨＨでクラスタαがラベリングされた信号成分のそれぞれにおいて、分散値Ｖａｒ＿ＨＬ＿α，Ｖａｒ＿ＬＨ＿α，Ｖａｒ＿ＨＨ＿αを計算し、その平均をクラスタαの特徴量Ｖａｒ＿αとする。同様に、クラスタβについても計算を行う。すなわち、特徴量は、次式（１）（２）で表される。
【数１】
2009098742.tif 000003
【数２】
2009098742.tif 000004
【００３５】
特徴量算出部１０６は、クラスタとそれに対応する特徴量を記録部１１１に保存する。なお、本実施の形態では、特徴量の例として分散値を用いているが、特徴量は分散値に限定されず、その他の統計量を用いてもよいし、複数の統計量を組み合わせて用いてもよい。
【００３６】
［コアリング閾値決定部］
次に、コアリング閾値決定部１０７について説明する。低域成分のクラスタ解析によって似通った性質同士の成分が分類されているが、クラスタ内の高域成分は均一ではない。特徴量算出部１０６にて求めた分散値は、クラスタに属するデータのばらつき具合を示している。すなわち、分散値が大きければばらつきが大きくノイズ成分が大きいと推測され、逆に分散値が小さければばらつきが小さくノイズ成分が小さいと推測される。従って、コアリング閾値決定部１０７は、分散値が大きい場合にはノイズ除去の強弱を制御するコアリング閾値を大きい方向に作用させ、分散値が小さい場合は、コアリング閾値を小さい方向に作用させる。
【００３７】
図７は、コアリング閾値決定部１０７が分散値に応じてコアリング閾値を決定するための関数７０１を示す図である。図７に示すように特徴量である分散値が大きくなるにつれてコアリング閾値を大きくすることにより、ノイズ成分の割合を考慮したノイズ除去が可能となる。なお、コアリング閾値を決定する関数７０１は、分散値が大きくなるにつれてコアリング閾値が大きくなる関係を満たせばよく、一次関数に限定されない。以上に、説明したコアリング閾値の決定法を閾値決定法aという。
【００３８】
また、コアリング閾値決定部１０７は、対象のクラスタにエッジが含まれるか否かに応じて、コアリング閾値を決定する。Ｈａａｒ基底を例とした場合、例えばレベルｎまでのウェーブレット変換、逆変換を行うには、最小で２n×２n画素の画像を必要とする。図８および図９は、例えばレベル２でのウェーブレット変換を例としてコアリング閾値を決定する処理を説明する図である。
【００３９】
図８は、各レベル（レベル１、レベル２）の低域成分８０１，８０２が一つのクラスタに所属する例を示す図である。この場合、対象領域はほぼ性質の似通った信号成分によって構成され、エッジが存在していないことが推測される。ここで、対象領域とは、レベルｎまでのウェーブレット変換、逆変換を想定した場合、２n×２n画素の矩形領域を指している。
【００４０】
図９は、各レベル（レベル１、レベル２）において、低域成分９０１，９０２が複数のクラスタに所属する例を示す図である。この場合、対象領域にはエッジが存在していると推測でき、エッジのボケを抑制するために、コアリング閾値を小さくする必要があることが分かる。
【００４１】
コアリング閾値決定部１０７は、ウェーブレット変換、逆変換の対象領域において、その低域成分のクラスタ分類結果とそれを基に計算された高域成分の分散値を考慮して決定する。
【００４２】
以上の閾値決定までの処理過程を、Ｈａａｒ基底を用いてレベル２までのウェーブレット変換を行った場合を例として説明する。例えば、レベル２までのウェーブレット変換で生成される低域成分を、クラスタ解析処理部１０４によりクラスタ解析を行い、図８に示すようにいずれのレベルにおいても低域成分が一つのクラスタに分類された場合、再構成対象領域にはエッジが存在しないことになる。従って、各レベルで生成される高域成分の分散値から閾値決定法ａによって算出された閾値を適用する。一方、クラスタ解析により解析を行い、図９に示すように各レベルの低域成分が複数のクラスタに分類された場合、エッジが存在する可能性が高い。従って、エッジ成分のボケを抑制するためにコアリング閾値を小さく設定する必要がある。この際のコアリング閾値の決定方法は、例えば、対象領域において一番多く所属するクラスタの高域成分から算出した分散値を特徴量とし、微調整用として事前に設定した係数（ｔ）を掛け合わせることにより、そのコアリング閾値とする。
【００４３】
図１０は、コアリング閾値の決定方法の一例を示す図である。図１０は、対象領域にα、β、γのクラスタが存在し、各クラスタに分類された信号成分のメンバー数と、分散値との関係を示している。メンバーの最も多いクラスタαの分散値に基づくコアリング閾値（ＴＨａ）に微調整用係数（ｔ）を掛け合わせた結果をコアリング閾値（ＴＨｂ）として適用する。この微調整用係数は、対象領域にエッジが存在していると推測されるので、エッジのボケ抑制の効果を目的としたものである。この決定手法を閾値決定法ｂとする。
【００４４】
コアリング閾値決定部１０７は、各レベルでコアリング閾値を決定することになるが、複数のレベルでコアリング閾値を決定する際に、（１）複数のクラスタに分類されたレベルのみ閾値決定法ｂを適用するか、（２）対象領域に、複数のクラスタに分類されたレベルが一つでもあれば、すべてのレベルにおいて閾値決定法ｂを用いるかは、ユーザが事前に選択し設定しておくことができる。ここで、レベル３までのウェーブレット変換を行った場合を例として説明する。レベル１において複数のクラスタに分類され、レベル２およびレベル３では一つのクラスタに分類された場合、レベル１のみ閾値決定法ｂを使い、それ以外のレベルでは閾値決定法ａを使うか、それとも、いずれかのレベルで、複数のクラスタに所属した対象領域に存在したということで、全レベルで閾値決定法ｂを適用するかは、ユーザが選択することができる。
【００４５】
図１１は、コアリング閾値決定部１０７の処理を示すフローチャートである。コアリング閾値決定部１０７は、まず、レベル１〜３のいずれのレベルにおいても低域成分が同一のクラスタに分類されているか否かを判定する（Ｓ１０）。すべてのレベルで同一クラスタに分類されていると判定された場合（Ｓ１０でＹＥＳ）、コアリング閾値決定部１０７は、各レベル１〜３において、分散値のみによる閾値決定法ａによりコアリング閾値を算出する（Ｓ１２）。すべてのレベルで同一クラスタに分類されていると判定されなかった場合（Ｓ１０でＮＯ）、複数に分類されたレベルのみ閾値決定法ｂを適用するか否かを判定する（Ｓ１４）。ここでは、ユーザによって予め設定された内容によって判定を行う。複数に分類されたレベルのみ閾値決定法ｂを適用する場合には（Ｓ１４でＹＥＳ）、複数に分類されたレベルのみエッジを考慮した閾値決定法ｂを適用し、それ以外のレベルに対しては閾値決定法ａを適用して閾値を決定する（Ｓ１６）。複数に分類されたレベルのみ閾値決定法ｂを適用するという設定でない場合には（Ｓ１４でＮＯ）、全てのレベルにおいてエッジを考慮した閾値決定法ｂを適用して閾値を決定する（Ｓ１８）。
【００４６】
［コアリング処理部］
次に、コアリング処理部１０８について説明する。コアリング処理部１０８は、ウェーブレット変換部１０３にて生成される高域成分信号と、コアリング閾値決定部１０７により決定されたコアリング閾値とを比較し、コアリング閾値より低い高域成分信号を０に置き換えることによりノイズを除去する。
【００４７】
図１２は、コアリング閾値をＴｈとしたときのコアリング入出力特性の例を示す図である。コアリング処理部１０８は、コアリング閾値決定部１０７で決定したコアリング閾値を適用することにより、対象領域のエッジの有無、ノイズの分布具合を考慮したコアリング処理を行うことができる。
【００４８】
［ウェーブレット逆変換部］
図１３は、ウェーブレット逆変換部１０９の構成を示す図である。ウェーブレット逆変換部１０９は、コアリング処理部１０８にてコアリング処理された高域成分と、ウェーブレット変換部１０３にて得られた低域成分に対してウェーブレット逆変換を行うことによって信号を再構成する。図１３に示すように、サブバンド分解過程で使用した基底関数から生成される再構成用ＬＰＦ１３０１，１３０３，１３０５、ＨＰＦ１３０２，１３０４，１３０６を用いて、帯域分割により得られたＬＬ，ＬＨ，ＨＬ，ＨＨの成分から信号を再構成することができる。
【００４９】
図１４（ａ）は、ウェーブレット逆変換の計算方法の例を示す図である。図１４（ａ）は、Ｈａａｒ基底を用いたレベル１のウェーブレット変換を行うことにより、サブバンド分解して得られた結果（ＬＬ，ＬＨ，ＨＬ，ＨＨ）から二次元の水平・垂直２画素［Ａ，Ｂ，Ｃ，Ｄ］再構成の計算方法を示している。図１４（ａ）に示す通り、Ａ＝（ＬＬ＋ＨＬ＋ＬＨ＋ＨＨ）／４，Ｂ＝（ＬＬ−ＨＬ＋ＬＨ−ＨＨ）／４，Ｃ＝（ＬＬ＋ＨＬ−ＬＨ−ＨＨ）／４，Ｄ＝（ＬＬ−ＨＬ−ＬＨ＋ＨＨ）／４によってウェーブレット逆変換を行うことができ、図１４（ｂ）に示すように再構成画像が得られる。
【００５０】
また、上位レベルから再構成されたＬＬ成分と下位レベルの高域成分を上記と同様の計算方法によって順次ウェーブレット逆変換を行うことにより、ノイズ成分の抑制を行いつつ、エッジのボケが抑制された鮮鋭な画像を得ることができる。
【００５１】
以上、本発明の実施の形態の画像処理装置１について説明したが、本発明は上記した実施の形態に限定されない。
【００５２】
上記した実施の形態では、画像処理装置１について説明したが、上記した画像処理装置１の各構成を実現するためのモジュールを有するプログラムも本発明の範囲に含まれる。
【００５３】
上記した実施の形態では、ウェーブレット変換における基底をＨａａｒの基底を用いた例について説明を行ったが、Ｈａａｒ基底以外の基底に基づいたウェーブレット変換、逆変換を利用した画像処理装置にも適用することが可能である。また、サブバンドへの分離として、分離するレベルも任意である。
【産業上の利用可能性】
【００５４】
本発明によれば、輪郭のボケを抑制しつつ、適切にノイズを低減できるというすぐれた効果を有し、ウェーブレット変換を利用してノイズを低減する画像処理装置、あるいは当該画像処理装置を搭載したカメラや携帯電話等として有用である。
【図面の簡単な説明】
【００５５】
【図１】本発明の実施の形態の画像処理装置のブロック図
【図２】ウェーブレット変換部の構成を示す図
【図３】Ｈａａｒ基底によるウェーブレット変換計算方法図
【図４】ウェーブレット変換によりレベル２までサブバンド分解した例を示す図
【図５】レベル１の低域成分をクラスタ解析した結果を示す図
【図６】レベル１の高域成分をラベリングした結果を示す図
【図７】閾値決定法ａを示す図
【図８】各レベルの低域成分が同一クラスタに分類された例を示す図
【図９】低域成分が複数のクラスタに分類された例を示す図
【図１０】閾値決定法ｂを示す図
【図１１】コアリング閾値決定部の処理を示すフローチャート
【図１２】コアリング入出力特性の一例
【図１３】ウェーブレット逆変換部の構成を示す図
【図１４】（ａ）ウェーブレット逆変換の計算方法の例を示す図（ｂ）ウェーブレット逆変換による再構成画像の例を示す図
【符号の説明】
【００５６】
１  画像処理装置
１０１  撮像素子
１０２  画像メモリ
１０３  ウェーブレット変換部
１０４  クラスタ解析処理部
１０５  ラベリング部
１０６  特徴量算出部
１０７  コアリング閾値決定部
１０８  コアリング処理部
１０９  ウェーブレット逆変換部
１１０  画像表示部
１１１  記録部
【図１】
2009098742.tif 000005
【図２】
2009098742.tif 000006
【図３】
2009098742.tif 000007
【図４】
2009098742.tif 000008
【図５】
2009098742.tif 000009
【図６】
2009098742.tif 000010
【図７】
2009098742.tif 000011
【図８】
2009098742.tif 000012
【図９】
2009098742.tif 000013
【図１０】
2009098742.tif 000014
【図１１】
2009098742.tif 000015
【図１２】
2009098742.tif 000016
【図１３】
2009098742.tif 000017
【図１４】
2009098742.tif 000018