JP2014165874A - 画像処理装置、画像処理方法、及びプログラム - Google Patents

Abstract

【課題】着目画素を基準にした類似度に応じて重みを算出して重みつきの平均値を計算すると、その平均値は着目画素自身に含まれるノイズの影響を受けので、ノイズのない真の画素値からずれた値に収束し、残留ノイズとして現れる。
【解決手段】着目画素に対する複数の参照画素の画素値の重み付き平均値と、着目画素の画素値と、着目画素の真の画素値と、の間に成り立っている相関関係に基づいて着目画素に対する新たな値を導出する。
【選択図】図４

Description

本発明は、画像処理装置、画像処理方法、およびプログラムに関する。より詳細には、本発明は、画像データのノイズを低減する技術に関する。

より暗いシーンあるいはより速く運動する被写体を撮像するために、カメラの感度を従来よりも高める事が望まれている。しかしながらカメラの感度を高めるとノイズが増えるため、ノイズをより強く減衰させるノイズ低減処理が必要となる。単純な平均化などによってノイズ低減処理を施すと解像感が低下するため、エッジやテクスチャを保存しつつノイズを低減する技術が公開されている（例えば特許文献１、特許文献２）。これらのノイズ低減技術は、ノイズを低減しようとしている着目画素と周辺の参照画素との類似度を計算し、その類似度に応じて重みを決定して重み付き平均値を計算する。そして重み付け平均値を用いて着目画素の値を決定するものである。

特開２０１０−２４５９７７号公報 特表２００７−５３６６６２号公報

しかしながら、このように着目画素を基準にした類似度に応じて重みを算出して重みつきの平均値を計算すると、その平均値は着目画素自身の画素値に含まれるノイズの影響を受ける。このため、参照画素の数を増やしたとしても着目画素の真の画素値には収束せず、ノイズのない真の画素値からずれた値に収束し、これが残留ノイズとして現れるという課題がある。

本発明に係る画像処理装置は、着目画素に対する複数の参照画素の画素値の重み付き平均値と、前記着目画素の画素値と、前記着目画素の真の画素値と、の間に成り立っている相関関係に基づいて前記着目画素に対する新たな値を導出する導出手段を備えることを特徴とする。

本発明によれば、残留ノイズを低減する事ができ、より良好なノイズ低減後の画像データを得る事ができる。

実施形態に係る画像処理装置の構成の一例を示すブロック図である。 着目画素と参照画素の例を示す図である。 non local means法における着目領域と参照領域との例を示す図である。 残留ノイズが生じる原因を示す概念図である。 残留ノイズを低減させた例を示す図である。 第１実施形態の画像処理部の構成の一例を示すブロック図である。 第１実施形態における処理の流れの一例を示すフローチャートである。 第２実施形態の画像処理部の構成の一例を示すブロック図である。 第２実施形態における処理の流れの一例を示すフローチャート図である。 第３実施形態の画像処理部の構成の一例を示すブロック図である。 第３実施形態における処理の流れの一例を示すフローチャートである。 第４実施形態の画像処理部の構成の一例を示すブロック図である。 第４実施形態における処理の流れの一例を示すフローチャート図である。

以下に図面を参照して、発明の実施の形態を詳しく説明する。なお、以下の実施の形態はあくまで例示であり、本発明の範囲を限定する趣旨のものではない。

＜第１実施形態＞
本実施形態における画像処理装置の構成について、図１を参照して説明する。

図１において、画像処理装置はＣＰＵ１０１、ＲＡＭ１０２、ＨＤＤ１０３、汎用インターフェース（Ｉ／Ｆ）１０４、モニタ１０８、メインバス１０９を備える。汎用Ｉ／Ｆはカメラなどの撮像装置１０５や、マウス、キーボードなどの入力装置１０６、メモリーカードなどの外部メモリ１０７をメインバス１０９に接続する。

以下では、ＣＰＵ１０１がＨＤＤ１０３に格納された各種ソフトウェア（コンピュータプログラム）を動作させることで実現する各種処理について述べる。

まず、ＣＰＵ１０１はＨＤＤ１０３に格納されている画像処理アプリケーションを起動し、ＲＡＭ１０２に展開するとともに、モニタ１０８にユーザインターフェース（ＵＩ）を表示する。続いて、ＨＤＤ１０３や外部メモリ１０７に格納されている各種データ、撮像装置１０５で撮像された画像データ、入力装置１０６からの指示データなどがＲＡＭ１０２に転送される。さらに、画像処理アプリケーション内の処理に従って、ＲＡＭ１０２に格納されているデータに対してＣＰＵ１０１からの指令に基づく各種演算が行れる。演算結果はモニタ１０８に表示されたり、ＨＤＤ１０３、外部メモリ１０７に格納されたりする。

上記の構成に於いて、ＣＰＵ１０１からの指令に基づき、画像処理アプリケーションにノイズ低減の対象となる画像データとノイズ低減パラメータとを入力し、ノイズ低減処理をかける。その上で、そのノイズ低減処理した画像データに対してさらに残留ノイズ低減処理を行う処理の詳細について説明する。まず、ノイズ低減技術について説明する。その後、ノイズ低減処理したデータにおける残留ノイズについて説明する。次に、その残留ノイズのノイズ低減処理について説明する。

＜ノイズ低減処理の説明＞
画像処理におけるノイズ低減処理は、一般に、ノイズを低減しようと着目している着目画素近傍に複数の参照画素を設定し、それぞれの参照画素の画素値と着目画素の画素値とを用いて平均値を計算する。そして、計算した平均値をノイズ低減後の着目画素の画素値とするものである。最も単純なノイズ低減処理は平滑化フィルタを掛ける処理である。例えば図２に示すように、着目画素を中心とした５×５サイズの領域にある全ての画素を参照画素群とし、それら参照画素群の画素値の平均値を計算し、着目画素の画素値とする。ただしこのような平滑化フィルタを掛ける処理では、エッジや細かなテクスチャがぼけて失われる。

エッジやテクスチャを保存するようなノイズ低減処理として最も簡単な処理は、平均値の計算に着目画素の画素値と参照画素の画素値との差分を考慮する処理である。差分がある閾値Ｓよりも小さい時にその参照画素を平均値の計算に使用し、逆に大きいときにはその参照画素は平均値の計算に使用しないと判定する。例えば、着目画素の画素値が５０であり、閾値Ｓを１０とするとき、画素値が４０から６０の範囲である参照画素は平均の計算に使用し、この範囲に無いものは平均の計算には使用しないと判定する。この判定により、エッジやテクスチャをある程度残すことができる。

上記の閾値を用いて平均値の算出に用いる参照画素を判定するノイズ低減処理よりもさらに高度な処理としては、例えばバイラテラルフィルタやNon-local means法が挙げられる。これらのエッジ保存型のノイズ低減処理は、一般に参照画素群の画素値の重み付き平均値を用いている。ここで、ノイズ低減処理を行おうとしている入力画像のi番目の画素を着目画素とし、また参照画素群におけるj番目の画素の画素値をIjとして、さらに着目画素に対応するノイズ低減処理後の画像のi番目の画素の画素値をIi,newとする。また、参照画素群に含まれる参照画素数をNpとする。このとき、ノイズ低減処理後の画像のi番目の画素値Ii,newは以下の式で書き表される。

ここで、Σは和をとる事を意味し、複数存在する参照画素を識別するインデックスjについて和を取ることを意味する。またwijは、i番目の着目画素におけるj番目の参照画素に対する重みである。なお、複数ある参照画素の中には着目画素が含まれる場合、上の式でj=iとなってもよい。

ノイズ低減処理方法の違いは参照画素に対する重みwijの計算方法の違いで説明される。一般に、重みwijは着目画素と参照画素との類似度に応じて決定される。例えば、先述した最も簡単な例では、重みwijは以下のようになる。

この処理においては類似度は画素値の差の絶対値を用いる。画素値の差の絶対値が小さい参照画素は、着目画素と参照画素は類似しているとみなし、着目画素の画素値を求めるための平均値算出に用いる。画素値の差の絶対値が大きい参照画素は、着目画素と参照画素は類似していないとみなし、着目画素の画素値を求めるための平均値算出に用いない。

また、バイラテラルフィルタでは、ｊ番目の参照画素に対する重みwijは以下の様に定義されている。

ここで、σIとσｒはノイズ低減の強さを決めるパラメータである。また、rijはi番目の着目画素とj番目の参照画素との画像上での距離である。画素値の差分が小さく距離的にも近い場合に着目画素に対する参照画素の類似度が高いとみなし、重みが大きくなるように定義した式となっている。

Non-local means法では、類似度の評価にブロックマッチングが取り入れられている。すなわち図３に示すように、着目画素周辺の領域を着目領域と設定し、また、参照画素周辺の領域を参照領域と設定する。そして着目領域に含まれる各画素と参照領域に含まれる各画素とで対応する画素毎に画素値の差分の二乗を計算し、さらにそれらの和を取って類似度Sijを計算している。式で表現すると以下のようになる。

ここでNbは参照領域画素数であり、Iikは、着目領域を構成しているk番目の画素の画素値であり、Ijkはj番目の参照画素に対応する参照領域を構成しているk番目の画素の画素値である。そしてｊ番目の参照画素に対する重みは以下の通り計算される。この類似度Sijは、着目領域と参照領域を構成する画素の画素値がそれぞれ完全に一致した場合は０となり、画素値の差が異なるほど大きな値となる。そして、重みwijは類似度Sijを使用して以下のように計算される。

ここでhはノイズ低減の強さを決めるパラメータである。なお、non-local meansという名前は、参照画素を画像の局所領域ではなく画像全体に設定する事からこう名付けられているが、参照画素を画像の局所領域に限定してもノイズ低減効果は得られる。

以上が代表的なエッジ保存型ノイズ低減技術の説明である。

＜残留ノイズの説明＞
次に、このようなエッジ保存型のノイズ低減技術である重み付き平均処理によって生じる残留ノイズについて、図４を用いて説明する。残留ノイズとは、重み付き平均処理を用いたノイズ低減処理を行う場合に、参照画素の数を増やしても低減されない（すなわち、残留する）特殊なノイズのことである。

図４は、従来技術のように着目画素と参照画素との類似度に応じた重み付き平均値を計算した場合に得られる期待値と真の画素値との間の差を示す概念図である。真の画素値とは、ノイズが全く含まれない理想値のことである。換言すれば、真の画素値とは、ノイズが完全に除去された値のことである。この図は簡単のために簡略化しているが、従来技術の期待値が真の画素値とならない事を示すには十分である。図４の横軸は画素値を示す。画素値分布は、各画素値に対する頻度の分布であり、重み関数は各画素値に対応する重みを示している。

まず、着目画素の真の画素値をIi,trueとする。ここで参照画素群の画素の真の画素値も全てIi,trueと同じ値だと仮定する。つまり、ノイズが発生していない場合には全ての画素値が同じである画像データにおいて、着目画素と参照画素群とについて考える。このとき、参照画素のノイズを含んだ画素値（すなわち、ノイズが発生した実際の画像データにおける参照画素の画素値）の分布は、真の画素値Ii,trueを中心に広がっている（図４の「画素値分布」参照）。ノイズを含んだ参照画素の画素値が、真の画素値Ii,trueを中心に広がって分布している理由は、各参照画素に対して、それぞればらついたノイズ量が含まれているからである。

ここで、ある着目画素の画素値Iiが、ノイズによって真の画素値Ii,trueよりも大きな値であったとする。上述したように、参照画素に対する重みは着目画素と参照画素との類似度に応じて決まる。着目画素と参照画素との類似度は、参照画素の画素値が着目画素の画素値Ｉｉから離れるに従って低くなるのが一般的である。従って重みも参照画素の画素値が着目画素の画素値Ｉiから離れるに従って減少する（図４の「重み分布」参照）。すなわち、参照画素に対する重みは、着目画素の画素値Ii近傍では高く、着目画素の画素値Iiから遠い部分では低くなる。この時計算される重み付き平均の値は、式（１）に示すように、参照画素の画素値と参照画素に対する重みとに基づく。図４からわかるように、参照画素の画素値が真の画素値Ii,true近傍かそれよりも小さい場合（図４のＡ）は、画素値分布の値が大きいためそのような値を持つ画素数はそれなりに多いものの対応する重みが小さい。また、参照画素の画素値が着目画素の画素値Ii近傍かそれよりも大きい場合（図４のＢ）には、重みは大きいもののその値を持つ画素数が少ない。結局、重み付き平均値に主に寄与するのは、図４において画素値分布における頻度と重みとがともに大きな値をもつ画素値である。したがって、重み付き平均値である着目画素の画素値Ii,newは、真の画素値Ii,trueと着目画素の画素値Iiとの間に分布することになる（図４の「従来技術結果の頻度分布」参照）。ここで、従来技術結果の確率分布における期待値をIi,new,expとする。

重み付き平均値のばらつき（分布の広がり）は参照画素の数を増やすほど減少していくが、その収束先は真の画素値Ii,trueではなく、期待値Ii,new,expとなる。この期待値Ii,new,expは、ノイズが含まれる着目画素の画素値Iiに依存し、真の画素値Ii,trueと期待値Ii,new,expとの差は参照画素の数を増やしても縮まる事はない。すなわち、重みを着目画素の画素値Iiに基づいて決める場合には、期待値Ii,new,expは、真の画素値Ii,trueとはならないのである。この真の画素値Ii,trueと、期待値Ii,new,expとの差が重み付き平均のばらつきよりも大きい場合に、この差が残留ノイズとして認識される。図４から判る通り、この残留ノイズは元々のノイズ（すなわち、着目画素の画素値Iiと真の画素値Ii,trueとの差）より小さい。したがって元々のノイズがそれほど大きくない場合には、残留ノイズはそれよりもさらに小さいために目立ちにくい。しかし元々のノイズ量が大きくなると、それにつれて残留ノイズも大きくなるため目立つようになる。特に、空や壁面など、均一な色が大面積で広がっている領域においてはこの残留ノイズが目立つため、この残留ノイズを低減する必要がある。以上が残留ノイズについての説明である。

なお、図４で説明した残留ノイズは、着目画素の画素値に基づいて重みを決定する場合に生じるものである。仮に着目画素の画素値に基づいて重みを決定しないような処理、例えば上述した平均フィルタを用いるような手法では、重みは定数となる。この場合、期待値Ii,new,expは真の画素値Ii,trueと一致する。したがって、図４の従来技術結果分布は、図４においては真の画素値Ii,trueを中心とする分布となる。換言すれば、図４で説明するような残留ノイズは着目画素の画素値に基づいて重みを決定する場合には避けることができない問題となる。

本実施形態ではこの残留ノイズを低減する処理を説明とする。

＜残留ノイズ低減処理の原理＞
次に、本実施形態の残留ノイズ低減処理の原理について説明を行う。先述したとおり、残留ノイズとはノイズ低減処理後の着目画素の期待値Ii,new,expと真の画素値Ii,trueとのずれである。そこで、以下では、ノイズ低減処理後の着目画素の期待値Ii,new,expと真の画素値Ii,trueとのずれを縮めることによって残留ノイズを低減するための原理について説明する。

まず、ノイズ低減処理後の着目画素の画素値Ii,newの期待値Ii,new,expは以下のように計算する。

ここで、式中の山括弧は統計的な期待値を計算する事を意味する。この式（６）に重み付き平均の式（１）を代入すると、

となる。重みwijは少なくとも着目画素の画素値Iiと参照画素の画素値Ijとに依存するとする。重みの具体的な計算方法は、個々のノイズ低減処理技術によるものである。あるノイズ低減処理技術において、重みの計算に１つ以上のパラメータξmethodを使用するならば、その重みはξmethodにも依存する。以下、ξmethodをノイズ低減パラメータと呼ぶ。ノイズ低減パラメータξmethodは、例えばバイラテラルフィルタではσIとσrが、non-local meansであればhがそれぞれ対応する。

次に、ノイズ低減処理後の着目画素の期待値Ii,new,expを統計的に計算する方法について説明する。なお、本実施形態において期待値を求めようとしている着目画素の画素値Iiについては、画像データから取得することができる。従って、着目画素の画素値Iiが定まっている条件の元で期待値Ii,new,expを計算する。ノイズ低減処理後の着目画素の画素値Ii,newの期待値Ii,new,expの計算は、次のように行われる。重みの計算に必要な画素は、図３が示す通り、着目画素と参照画素の他に、参照画素群には含まれないがいずれかの参照画素に対応する参照領域には含まれる画素である。ここで参照画素群には含まれないが重みの計算に必要な画素の画素値をIothersとする。重みの計算に関わる着目画素および参照画素以外のあらゆる画素の画素値Iothersと参照画素の画素値Ijとの確率分布関数をノイズ低減処理後の着目画素の画素値Ii,newに乗じて、IjとIothersとについて積分を取ることで計算される。ここで画素値Iothersはベクトルとしておく。参照画素の画素値Ijと、着目画素および参照画素以外のあらゆる画素の画素値Iothersと、の画素値分布は、それぞれの画素の真の画素値とノイズ量とによって決まる。これらの画素値分布を決める１つ以上のパラメータをまとめてσtrueであらわす。以下、σtrueを画素値分布モデルパラメータと呼ぶ。以上の点を踏まえ式（７）を変形すると、

となる。上記の通り、重みwijは、少なくとも着目画素の画素値Iiと参照画素の画素値Ijとノイズ低減パラメータξmethodとに依存するので、式（８）ではwij内にこれらの値を明示している。なお、式（８）でwij内にIothersを含めている理由は、たとえばnon local meansの場合には参照領域を構成する、参照画素以外の画素も重みの計算に用いられることを考慮しているからである。すなわち、重みの計算方法によっては、着目画素の画素値Iiと参照画素の画素値Ij以外の画素値Iothersに依存する場合もあるので、式（８）のwijには画素値Iothersも含めている。P(Ij,Iothers,σtrue)は、上述した参照画素の画素値Ijと着目画素および参照画素以外のあらゆる画素の画素値Iothersとの確率分布関数を示している。また、P(Ij,Iothers,σtrue)は、これらの画素値の分布に関わる画素値分布モデルパラメータσtrueをパラメータとして含めている。確率分布関数P(Ij,Iothers,σtrue)は、撮像センサーのノイズの統計的な性質と、真の画像の画素値などから具体的な形にモデル化できる。すなわち、画素値分布モデルパラメータσtrueに応じて具体的な形にモデル化できる。例えば、ノイズ分布として標準偏差σを持つガウス分布を仮定し、参照画素の画素値Ijと着目画素および参照画素以外のあらゆる画素の画素値Iothersとの真の画素値を、それぞれIj,trueと Iothers,trueと仮定する。すると、確率分布関数P(Ij,Iothers,σtrue)は以下のような式で表せる。

ここでNothersはIothersの要素数を表している。参照画素の真の画素値Ij,trueと 他の画素の真の画素値Iothers,trueとを仮定している理由は、画素値の分布は、真の画素値を中心にノイズの標準偏差σに応じてばらつくことを示すためである。式（９）においては、参照画素の真の画素値Ij,true,他の画素の真の画素値Iothers,true,およびノイズの標準偏差σが、式（８）で示す画素値分布モデルパラメータσtrueに対応している。式（９）はノイズ分布に標準偏差σを持つガウス分布を仮定した例を示したものであるが、他の分布に従う場合には、式（８）で示す画素値分布モデルパラメータσtrueは、その分布に従った関数や値に対応させることができる。

次に、ノイズ低減後の着目画素の画素値の期待値の計算に説明を戻す。

式（８）の重みWijに、個々の重み付け平均に応じた重み計算方法を代入する事で、個々の重み付け平均ごとの期待値Ii,new,expを計算できる。式（８）の右辺には様々な変数が表れている。しかし、積分によって消える変数（Ij, Iothers）を考慮する。すると、ノイズ低減処理後の着目画素の画素値の期待値Ii,new,expは、着目画素の画素値Ii、ノイズ低減パラメータξmethodそして画素値分布モデルパラメータσtrueに依存する事がわかる。すなわち、ノイズ低減処理後の着目画素の画素値の期待値Ii,new,expは、着目画素の画素値Ii、ノイズ低減パラメータξmethodと、画素値分布モデルパラメータσtrueとに応じて求められる。
これを形式的に表現すると、以下の様な関数hで表す事ができる。

上記式を考慮する前に、まずノイズ低減処理の根本的な考え方について言及する。ノイズ低減処理は参照画素のいくつかにその参照画素の真の画素値が着目画素の真の画素値と同じものが含まれていないと成立しない。特に、重み付け平均を行う場合には、着目画素の画素値と同じ画素値を有する参照画素に対して重みを重くするという手法を採用していることからも上記の考えは裏付けされる。したがって、着目画素の真の画素値Ii,trueと同じ値が、画素値分布モデルパラメータσtrueに含まれていると考えても差し支えない。その事を明示的に示すために、画素値分布モデルパラメータσtrueから真の画素値Ii,trueを分離して式（１０）を書きなおすと、以下の様な関数fであらわす事ができる。

式（１１）は、図４に示すグラフと合わせるとその意味が把握しやすい。式（１１）の着目画素の真の画素値Ii,trueとノイズ低減処理後の着目画素の画素の期待値Ii,new,expとは、図４のIi,trueとIi,new,expとにそれぞれ対応する。式（１１）の画素値分布モデルパラメータσtrueの一部が、図４の画素値分布に対応する。式（１１）のノイズ低減パラメータξmethodの一部が、図４の重み関数に対応する。

この式（１１）を、真の画素値Ii,trueについて解く。これにより、真の画素値Ii,trueを着目画素の画素値Iiと、ノイズ低減処理後の画素値の期待値Ii,new,expと、ノイズ低減パラメータξmethodと、画素値分布モデルパラメータσtrueとの関数gとして求める事ができる。

重みの計算方法、つまりノイズ低減処理技術の種類によっては、これらの計算を解析的に行って関数ｇを得る事は困難である事が多い。しかし、様々な近似を用いて大まかに評価する事は可能であるし、また、コンピュータを使って数値的に関数ｇを評価する事はいつでも可能である。

式（１２）は入力画像データがσtrueによって特定されるある画素値分布モデルに従っている場合、次のことを意味する。すなわち、着目画素の真の画素値Ii,trueが、ξmethodによって特定されるノイズ低減処理によるノイズ低減処理後の画素値の期待値Ii,new,expと着目画素の画素値Iiとから決まる事を意味する。ノイズ低減後の着目画素の画素値Ii,newは、図４の従来技術結果の頻度分布に示すように、ノイズ低減処理後の画素値の期待値Ii,new,expの周辺にいくらかばらついている。期待値は、試行を行ったとき，その結果として得られる数値の平均値のことであり、実際の処理結果は期待値の周辺にばらつくことは周知である。

式（１２）の意味を言い変えると、着目画素の画素値Iiとノイズ低減後の画素値Ii,newと着目画素の真の画素値Ii,trueにはノイズ低減処理方法に応じた相関関係が成りたっているということである。ノイズ低減処理方法に応じた相関関係とは、画素値の分布を特定するパラメータσtrueとノイズ低減処理を特定するパラメータξmethodとに応じた相関関係のことである。

本実施形態ではこの相関関係に基づいて残留ノイズの低減処理を行う。式（１２）によって真の画素値Ii,trueを推定にするには、着目画素の画素値Iiと、ノイズ低減パラメータξmethodと、画素値分布モデルパラメータσtrueの他に、ノイズ低減処理後の画素値の期待値Ii,new,expが必要である。上述したように、ノイズ低減後の着目画素の画素値Ii,newは、ノイズ低減処理後の画素値の期待値Ii,new,expの周辺にいくらかばらついている。したがって、厳密にはノイズ低減処理によって得られたIi,newが期待値Ii,new,expと一致しない場合があるが、ノイズ低減後の着目画素の画素値Ii,newは期待値Ii,new,expのよい推定値である。よって本実施形態では、実際に得られる重み付き平均値（すなわち、ノイズ低減後の画素値）Ii,newを期待値Ii,new,expだとみなして、重み付き平均値Ii,newを式（１２）の期待値Ii,new,expに代入する。これによって着目画素の画素値Ii、ノイズ低減処理後の画素値Ii,new、パラメータξmethodおよびσtrueから真の画素値Ii,trueの推定値を算出する事ができる。このようにして得た推定値には上述の残留ノイズは低減されている。この推定値に残る主なノイズは、ノイズ低減処理によって得られる重み付き平均値Ii,newがその期待値Ii,new,expに対して持つばらつきに起因するものとなる。

なお一般的には式（１２）の通りであるが、着目画素周辺領域や参照画素周辺領域が着目画素と同一の画素値を持ち、さらにノイズ低減処理がノイズ分散を考慮しているような場合には、関数gが画素値分布モデルパラメータσtrueに依存しない場合がある。つまり以下に示すより簡単な式になる場合がある。

さらに、関数gがIiとIi,new,expに対して線形であるか、あるいは近似的に線形である場合、式（１３）は以下のように表される場合もある。

ここで、γはノイズ低減パラメータξmethodから計算されるパラメータであり、１よりも大きな値を持つ。典型的には1~1.5程度の値である。つまり、式（１４）においては、真の画素値Ii,trueはノイズ低減処理前の画素値Iiとノイズ低減処理後の画素値Ii,newのIi,new側（重み付き平均値側）の外分点であると解釈できる。これは図４のIi,true、Ii,new,exp、Iiの大小関係からも理解できる。

前述の通り、関数ｇの評価には真の画素値の分布に対するモデル（画素値分布モデルパラメータσtrue）が必要である。例えば空や壁のようにどの画素も一定の値を持つのか、それともエッジが存在するのか、といったモデルが考えられる。厳密にいえば、ノイズ画像からは真の画素値の分布は知り得ないため、どのような画像に対しても誤差無く厳密に残留ノイズを消す事は不可能である。しかしながら、残留ノイズの量は元々のノイズよりも小さく、また残留ノイズが特に目立つのは主に画素値がほぼ一定であるような領域である。そのため、画素値一定という仮定をおいて評価した関数ｇを利用して残留ノイズ低減処理を行っても画質改善効果が得られる。より精度が必要であれば、画像の領域毎にその領域が一定の画素値を持つ部分であるか、あるいはエッジを有するかなどを推定し、その結果に基づいて適した関数ｇを選択する事が可能である。また、上述したように、関数gの評価は、あらかじめ様々な近似を用いて大まかに評価しておいてもよいし、また、コンピュータを使って評価していてもよい。あるいは、関数gに対応する４次元ルックアップテーブルを用意しておいて、かかるテーブルに基づいて処理を行っても良い。関数gは、ノイズ低減処理方法や画素値分布などに応じて変化するものである。よって、それぞれのノイズ低減処理方法（たとえばnon local mensやバイラテラルフィルタなど）や、画素値分布（エッジかベタかなど）に応じて適切な関数gを適用することができる。

また、上記の例では、ノイズ低減処理を単一の画像データに対して行う例を説明したが、本実施形態は、複数の画像データを用いたノイズ低減処理に適用することも可能である。たとえば、同一の被写体に対して連続して撮像した画像データを用いるノイズ低減処理にも適用することができる。

また、本実施形態は、辞書（あらかじめ用意された複数の画像）を利用するタイプのノイズ低減処理にも適用可能である。辞書を利用するタイプのノイズ低減処理では、参照画素と参照領域を事前に用意した辞書の中から選び、着目領域と比較する。このような場合、式（８）にある確率分布関数P(Ij,Iothers,σtrue)のうち、次の要素が辞書によって決まる。すなわち、参照画素Ijと、着目画素および参照画素以外のあらゆる画素の画素値Iohtersおよび画素値分布モデルパラメータσtrueの要素のうちの参照領域に対応する要素と、が辞書によって決まる。従って、それらの値の確率分布を考える必要がない。この事を式によって明記すると、以下のようになる。

ここで、dicは辞書の要素数を示す。また、It,othersは、Iothersの要素のうち着目領域に対応したものを指す。たとえば、non local meansの場合には、着目領域を構成する画素の画素値のうち、着目画素Iiを除いた画素に対応する。同様に、σt,trueもσtureの要素のうち、着目領域に関連したものを指す。上述したように、参照領域に対応する画素については辞書の中から選択されるので、式（１５）の確率分布関数P(It,others, σt,ture)には含まれていない。辞書を利用するタイプのノイズ低減処理においても、辞書と着目画素の真の画素値の間に相関がなければノイズ低減処理としてそもそも成り立たないことを考慮すると、式（１５）も式（１１）さらには式（１２）の形に帰着できる。ただし画素値分布モデルパラメータσtrueは、着目領域に関連したσt,trueである。従って、辞書を利用するタイプのノイズ低減処理においても同様に本実施形態で説明した残留ノイズ低減処理が適用可能である。

以上が本実施形態における残留ノイズ低減処理の原理であるが、最後に図５を用いて残留ノイズ低減処理の効果の実例を示す。図５(a)はノイズのない画像の５×５の領域の画素値を示したものである。この例では全ての画素値が100である。すなわち、図５(a)は、真の画素値が全て１００である画像データであるといえる。図５(b)は、図５(a)の画像データに標準偏差が約24のガウス分布に従うノイズがのった画像の、ある５×５の領域の画素値を示したものである。図５(b)の画像データは、ノイズ低減の対象となる画像データの一部である。この領域だけでノイズの標準偏差を計算すると25.6であった。図５(c)はノイズをのせた図５(b)の画像データにnon-local meansフィルタを適用して得た画像の、同じ５×５の領域の画素値を示したものである。non-local meansフィルタを適用した結果、ノイズの標準偏差は3.57にまで減り、ノイズが低減されている事がわかる。しかしながら、図５(b)で画素値が真の画素値である100よりも小さい画素は、non-local meansフィルタを適用した図５(c)でも100より小さい傾向にある。逆に、図５(b)で100よりも大きな画素値をもつ画素では、non-local meansフィルタを適用した図５(c)でも100より大きな画素値を持つ傾向がある事が確認できる。これが残留ノイズである。図５(d)は、本実施形態で説明した処理によって残留ノイズを低減させた結果である。具体的には式（１４）に従って処理を行い、パラメータγを1.167とした例である。例えば最も左上の画素でみると、Ii,new,expには図５(c)に含まれている102を、Iiには図５(b)に含まれている110を、それぞれ代入する。式１４は1.167×102 - 0.167×110 = 100.664となるから、残留ノイズ低減後の画素値はこれを四捨五入した101となる。このような処理の結果、図５(d)におけるノイズの標準偏差は0.73となり、図５(c)よりもさらにノイズが低減し、真の画素値である図５(a)に近づいている事が確認できる。

＜残留ノイズ低減処理＞
次に、具体的な残留ノイズ低減処理を行う画像処理装置およびその方法について、図６のブロック図と図７のフローチャートを用いて説明する。ここでは、前述のnon-local means法のように、着目画素と参照画素との類似度について、着目領域における各画素と参照領域における各画素との差分に基づいて類似度を算出する方法を例とする。

図６は、本実施形態にかかる画像処理装置の画像処理部６００の一例を示すブロック図であり、入力画像に対して画素毎に処理を行う。この処理を全ての画素に対して行うことで、残留ノイズの少ない画像を得ることができる。画像処理部６００は、重み付け平均値導出部６０１と、残留ノイズ低減部６０２とを有する。

重み付き平均値導出部６０１は、全ての参照領域の画素情報と、着目領域の画素情報と、ノイズ低減パラメータとを入力し、重み付き平均値を出力する。残留ノイズ低減部６０２は、重み付き平均値導出部で出力された重み付き平均値と、ノイズ低減パラメータと、着目画素の画素値と、必要に応じてノイズ標準偏差とを入力し、残留ノイズ低減後画素値を出力する。

画像処理全体としては、まず重み付き平均値導出部６０１においてノイズ低減処理を行い、その後、残留ノイズ低減部６０２において残留ノイズを低減する処理を行う。

図７のステップＳ７０１において、まず、重み付き平均値導出部６０１は、着目領域の画素情報と全ての参照領域の画素情報とを取得する。ここで、画素情報とは、例えば画素値や画素位置などの情報を含む。次に、ステップＳ７０２において、重み付き平均値導出部６０１は、ノイズ低減パラメータを取得する。なお、ステップＳ７０１とＳ７０２は逆の順序でも並行して行われても良い。ステップＳ７０３において、重み付き平均値導出部６０１は、ステップＳ７０１で得た着目領域の画素情報及び全ての参照領域の画素情報と、ステップＳ７０２で得たノイズ低減パラメータとに基づいて重み付き平均値を導出する。

ここで、着目領域と参照領域とは、それぞれ着目画素と参照画素について類似度を計算するために必要な領域を指す。着目領域および参照領域の具体的な形状は、類似度と重みをどのように計算するかに依存する。例えばnon-local meansフィルタの場合はブロックマッチングを行うため、そのブロックの形状が着目領域と参照領域の形状である。バイラテラルフィルタの場合は、着目画素と参照画素の画素値の差と画素間の距離とによって類似度を計算するため、着目領域には着目画素そのもの、および参照領域には参照画素そのものがそれぞれ対応する。

ノイズ低減パラメータには例えばノイズの標準偏差あるいは分散などが含まれうる。一般にノイズの標準偏差あるいは分散などは感度などの撮像パラメータに依存するが、撮像パラメータを考慮した値を取得する事はよく知られておりここでは詳しい説明を省略する。この重み付き平均値を導出するまでの処理については、＜ノイズ低減処理の説明＞で説明したような公知技術を用いて行うことができる。重み付き平均値導出部６０１は、導出した重み付き平均値を残留ノイズ低減部６０２に出力する。

次に、ステップＳ７０４において、残留ノイズ低減部６０２は、重み付き平均値導出部６０１で導出された重み付き平均値、着目画素の画素値、ノイズ低減パラメータ、さらには必要に応じてノイズの標準偏差を取得する。すなわち、上述したノイズ低減後の画素値Ii,newと、着目画素の画素値Iiと、ノイズ低減パラメータξmethodと、必要に応じてノイズの特性を表す画素値分布モデルパラメータσtrueとを取得する。

次にステップＳ７０５において、残留ノイズ低減部６０２は、先述した通り関数ｇに従って残留ノイズを低減した出力値（画素値）を導出する。上述したように、関数ｇは、ノイズ低減パラメータξmethodや画素値分布モデルパラメータσtrueに応じて変化するものであり、個々の具体的な処理については説明を省略する。いずれの処理によっても、着目画素の真の画素値Ii,trueはノイズ低減処理前の画素値Iiとノイズ低減処理後の画素値Ii,newのIi,new側の外分点となり、残留ノイズが低減される。

ステップＳ７０６において、画像処理部６００は、全ての画素を処理したかどうかを判定する。すなわち、全ての画素を着目画素としてステップＳ７０１からステップＳ７０５の処理を行ったかを判定する。画像処理部６００が処理していないと判定した場合、ステップＳ７０１に戻り、処理していると判定した場合、処理を終了する。

以上のとおり、本実施形態によれば、着目画素と参照画素との類似度に応じて重みを計算し、重み付き平均値を導出してノイズを低減する場合に生じる残留ノイズを低減する事ができる。したがって、より良好なノイズ低減後の画像データを得る事ができる。

＜第２実施形態＞
第２実施形態における画像処理装置の構成は第１実施形態の図１と同様であるため説明を省略し、以下、第１実施形態と異なる部分について説明する。

第１実施形態は、近似などを用いて関数ｇの形が解析的に評価できる場合の実施形態である。一方、第２実施形態は、関数ｇの形が解析的に評価できない場合について説明する。関数ｇの形が解析的に評価できない場合は、あらかじめ数値計算をして補正テーブルを作成しておく。そして、パラメータに応じて補正テーブルを参照して残留ノイズの低減を行う。

以下、第２実施形態の処理に関して、図８のブロック図と図９のフローチャートを用いて説明する。図８も図６と同様、入力画像に対して画素毎に処理を行う部分のブロック図である。図８に示す画像処理部８００は、重み付き平均値導出部８０１と、補正テーブル記憶部８０２と、残留ノイズ低減部８０３とを有する。図６と異なる点は、補正テーブル記憶部８０２を有している点である。残留ノイズ低減部８０３は補正テーブル記憶部に格納されている補正テーブルを参照して残留ノイズ低減処理を行う。

補正テーブルは、重み付き平均値、着目画素値、ノイズ低減パラメータ、および画素値分布モデルパラメータより残留ノイズ低減後の画素値を得るために用いられる。この補正テーブルは関数ｇを数値的に評価して作成したものである。あるいは補正テーブルは、シミュレーションによって関数ｇを求めて作成したものでもよい。シミュレーションによって関数ｇを求める場合は、例えば以下のようにすればよい。まずＣＧなどによって作られた、ノイズが非常に少ない既知の画像を示す正解画像データを生成する。次にその正解画像データに対してノイズを付加してノイズ付加画像データを生成し、さらに所定の公知技術によるノイズ低減処理を行いノイズ低減画像データを生成する。そして最後に、正解画像データとノイズ付加画像データとノイズ低減画像データとの間の相関を調べ、関数ｇの形状を決定する。正解画像データ、ノイズ付加画像データ、ノイズ低減画像データがそれぞれIi,true、Ii、Ii,new,expに対応する。

次に、具体的な残留ノイズ低減処理方法について、図８のブロック図と図９のフローチャートを用いて説明する。なお、上述したような手法によって作成された補正テーブルが、補正テーブル記憶部に格納されているものとする。ステップＳ９０１からステップＳ９０３までの処理については第１実施形態で説明したものと同様であるため、説明を省略する。ステップＳ９０４において、残留ノイズ低減部８０３は、重み付き平均値導出部８０１で導出された重み付き平均値、着目画素の画素値、ノイズ低減パラメータ、さらには必要に応じてノイズの標準偏差を取得する。そしてステップＳ９０５において、残留ノイズ低減部８０３はステップＳ９０４で得た画素値やパラメータを基に補正テーブル記憶部８０２から残留ノイズを補正するための値を取得する。ステップＳ９０５では、残留ノイズ低減部８０３は、例えば着目画素の真の画素値Ii,trueやノイズ低減後の画素の期待値Ii,new,expなどの値を取得する。次に、ステップＳ９０６において、残留ノイズ低減部８０３は、ステップＳ９０４及びステップＳ９０５で得られた情報に基づいて、先述した通り関数ｇに従って残留ノイズ低減後の画素値を導出する。

ステップＳ９０７において、画像処理部８００は、全ての画素を処理したかどうかを判定する。画像処理部８００は、全ての画素を処理していないと判定した場合、ステップＳ９０１に戻り、処理していると判定した場合、処理を終了する。

以上により、本実施形態によれば、関数ｇの形が解析的に評価できない場合であっても残留ノイズの少ない良好なノイズ低減後の画像データを取得する事ができる。

＜第３実施形態＞
第３実施形態では、画素値がほぼ一定である領域に対してのみ残留ノイズ低減を行う処理について説明する。第３実施形態における画像処理装置の構成は第１実施形態の図１と同様であるため説明を省略する。

以下、第３実施形態の処理に関して、図１０のブロック図と図１１のフローチャートを用いて説明する。図１０も図６と同様、入力画像に対して画素毎に処理を行う部分のブロック図である。図１０に示す画像処理部１０００は、重み付き平均値導出部１００１と、残留ノイズ低減部１００２とを有する。図６と異なる点は、残留ノイズ低減処理を行う残留ノイズ低減部１００２に、画素値一定領域判定情報を入力する点である。画素値一定領域判定情報は、着目画素周辺の領域において、各画素の画素値がほぼ一定であるかどうかを判定した結果を示す情報である。画素値がほぼ一定である領域かどうかの判定は、公知技術を用いればよく、例えば、着目画素周辺の画素値の分散を計算し、それがノイズの分散の２倍程度より小さければほぼ画素値が一定である領域だと判定する事ができる。

以下、図１１を用いて具体的な残留ノイズ低減処理方法について述べる。ステップＳ１１０１からステップＳ１１０３までは第１実施形態と同様の処理であるので、説明を省略する。

ステップＳ１１０４において、残留ノイズ低減部１００２は着目画素周辺の画素値がほぼ一定かどうかを判定した結果を示す情報を取得する。ステップＳ１１０５において、残留ノイズ低減部１００２は、ステップＳ１１０４で取得した情報に基づいて残留ノイズ低減処理を行うか否かを判定する。ステップＳ１１０５において、着目画素周辺の画素値がほぼ一定であると判定した場合、残留ノイズ低減部１００２は、処理をステップＳ１１０６に進める。一方、ステップＳ１１０５において、着目画素周辺の画素値がほぼ一定でないと判定した場合、残留ノイズ低減部１００２は、処理をステップＳ１１０８に進める。

画素値がほぼ一定であると判定した場合、ステップＳ１１０６において、残留ノイズ低減部１００２は、ステップＳ１１０３で導出された重み付き平均値、着目画素の画素値、ノイズ低減パラメータ、さらには必要に応じてノイズの標準偏差を取得する。ステップＳ１１０７において、残留ノイズ低減部１００２は、関数ｇに従って残留ノイズを低減した画素値を導出して、出力する。

一方、画素値が一定領域ではないと判定した場合には、ステップＳ１１０８において、残留ノイズ低減部１００２は、重み付き平均値を残留ノイズ低減後の画素値として出力する。ステップＳ１１０９において、画像処理部１０００は、全ての画素を処理したかどうかを判定する、画像処理部１０００は、全ての画素を処理していないと判定した場合、ステップＳ１１０１に戻り、処理したと判定した場合、処理を終了する。

本実施形態によれば、画素値がほぼ一定である領域を示す情報を用いることにより、残留ノイズが認識されやすい領域に対してのみ残留ノイズ低減処理を行うことができる。

＜第４実施形態＞
第３実施形態では、画素値がほぼ一定である領域に対してのみ残留ノイズ低減を行う処理について説明した。第４実施形態では、着目画素周辺のテクスチャを詳しく識別した情報が得られた場合に、それに応じて補正処理を切り替える例を示す。具体的には着目画素が(1)エッジ部、(2)画素値一定領域、(3)その他の領域と３つのカテゴリに分類される場合に、エッジ部か否かで補正に用いるテーブルを変更する例を示す。

第４実施形態における画像処理装置の構成は第１実施形態の図１と同様であるため説明を省略する。

以下、第４実施形態の処理に関して、図１２のブロック図と図１３のフローチャートを用いて説明する。図１２も図８と同様、入力画像に対して画素毎に処理を行う部分のブロック図である。図１２の画像処理部１２００は、重み付き平均値導出部１２０１と、補正テーブル記憶部１２０２と、残留ノイズ低減部１２０３とを有する。図１２が図８と異なる点は、残留ノイズの低減処理を行う残留ノイズ低減部１２０３が、着目画素が画素値一定の領域であるのか、エッジを含んだ領域であるのか、あるいはいずれでもないのかを判定した結果を示すテクスチャ判定情報が入力される事である。エッジである場合のテクスチャ判定情報にはエッジの両サイドの画素値あるいはそれらの差が含まれる。これは、画素値分布モデルパラメータσtrueの１つである。なお、補正テーブル記憶部１２０２に格納されている補正テーブルは、画素値一定領域とエッジ部とのそれぞれに対してシミュレーションによって関数ｇを評価して作成したものである。すなわち、補正テーブル記憶部１２０２には、画素値一定領域用の補正テーブルと、エッジ部用の補正テーブルとが格納されている。

以下、図１２を参照して具体的な残留ノイズ低減処理方法について述べる。ステップＳ１３０１からステップＳ１３０３までの処理については第１実施形態と同様であるので説明を省略する。

ステップＳ１３０４において、残留ノイズ低減部１２０３は着目画素周辺のテクスチャ判定情報を取得する。そしてステップＳ１３０５において、残留ノイズ低減部１２０３はテクスチャ判定情報に基づいて残留ノイズ補正処理を行うか否かを判定する。例えば、残留ノイズ低減部１２０３は、テクスチャ判定情報が、着目画素がエッジ部であるかまたは画素値一定領域であることを示す場合、残留ノイズ低減処理を行うと判定する。これらの判定は、たとえば公知の各種フィルタを入力画像データに対して適用することで行うことができる。一方、テクスチャ判定情報がその他の領域を示す場合、残留ノイズ低減処理を行わないと判定する。

残留ノイズ低減処理を行うと判定した場合、ステップＳ１３０６において、残留ノイズ低減部１２０３は、ステップＳ１３０３で導出された重み付き平均値、着目画素の画素値、ノイズ低減パラメータと、必要に応じてノイズ標準偏差を取得する。ステップＳ１３０７において、残留ノイズ低減部１２０３は、ステップＳ１３０４で得たテクスチャ判定情報を基にエッジ部かどうかを判定する。

エッジ部であると判定した場合、残留ノイズ低減部１２０３はステップＳ１３０８に処理を進める。ステップＳ１３０８において、残留ノイズ低減部１２０３は、エッジ部に対応する補正テーブルを参照して補正テーブル記憶部１２０２から残留ノイズを補正するための値を取得する。

エッジ部でないと判定した場合、残留ノイズ低減部１２０３はステップＳ１３０９に処理を進める。ステップＳ１３０９において、残留ノイズ低減部１２０３は、画素値一定領域に対応する補正テーブルを参照して、補正テーブル記憶部１２０２から残留ノイズを補正するための値を取得する。

ステップＳ１３１０において、残留ノイズ低減部１２０３は、ステップＳ１３０８またはステップＳ１３０９で得た値を用いて残留ノイズ低減後の画素値を導出して出力する。

ステップＳ１３０５にて、残留ノイズ補正処理を行わないと判定した場合には、残留ノイズ低減部１２０３は、ステップＳ１３１１において重み付き平均値を低減した画素値として出力する。

ステップＳ１３１１において、画像処理部１２００は全ての画素を処理したかどうかを判定する。全ての画素を処理していないと判定した場合、画像処理部１２００は、ステップＳ１３０１に戻り、処理していると判定した場合、処理を終了する。

本実施形態によれば、着目画素のテクスチャを示す情報を用いることにより、より好適な残留ノイズ処理を着目画素毎に選択して実行することが可能となる。

なお、第４実施形態ではその他の領域と判定される着目画素に対しては重み付き平均値を出力する例を説明したが、画素値一定領域だけでなくその他の領域と判定された着目画素に対しても画素値一定領域と共通の補正処理を行ってもよい。

＜その他の実施形態＞
本発明は、上述した実施形態の機能を実現するソフトウェアのプログラムコードを記録した記憶媒体を、システム或いは装置に供給することによっても実現できる。この場合、そのシステム或いは装置のコンピュータ（又はＣＰＵやＭＰＵ）がコンピュータが読み取り可能に記憶媒体に格納されたプログラムコードを読み出し実行することにより、上述した実施形態の機能を実現する。

Claims (14)

1. 着目画素に対する複数の参照画素の画素値の重み付き平均値と、前記着目画素の画素値と、前記着目画素の真の画素値と、の間に成り立っている相関関係に基づいて前記着目画素に対する新たな値を導出する導出手段を備えることを特徴とする画像処理装置。
2. 前記重み付き平均値は、着目画素の画素値に応じた重みによる重み付き平均によって得られることを特徴とする請求項１に記載の画像処理装置。
3. 前記導出手段は、さらに前記着目画素が有するノイズの特性に基づいて前記新たな画素値を導出することを特徴とする請求項１または２に記載の画像処理装置。
4. 前記ノイズの特性は前記ノイズの標準偏差であることを特徴とする請求項３に記載の画像処理装置。
5. 前記相関関係は、前記重み付き平均値を算出するためのノイズ低減パラメータと、画素値分布モデルパラメータとに応じて特定されることを特徴とする請求項１から４のいずれか一項に記載の画像処理装置。
6. 前記導出手段は、前記相関関係に基づいて予め作成されたテーブルを参照する事によって前記新しい値を導出する事を特徴とする請求項１から５のいずれか一項に記載の画像処理装置。
7. 前記導出手段は、前記着目画素の画素値と前記重み付き平均値との前記重み付き平均値側の外分点を前記新たな値として導出することを特徴とする請求項１から６のいずれか一項に記載の画像処理装置。
8. 画素値が一定の領域に前記着目画素が属するかを判定する判定手段と、
   前記判定手段によって画素値が一定の領域に前記着目画素が属すると判定された場合、前記導出手段によって導出された前記新しい値を前記着目画素の値として出力し、
   前記判定手段によって画素値が一定の領域に前記着目画素が属すると判定されない場合、前記重み付き平均値の値を前記着目画素の値として出力する出力手段と
   をさらに備えることを特徴とする請求項１から７のいずれか一項に記載の画像処理装置。
9. 着目画素の周辺のテクスチャを識別するテクスチャ判定情報を取得する取得手段をさらに備え、
   前記導出手段は、前記取得手段で取得したテクスチャ判定情報に基づいて前記新しい値を導出するかを決定することを特徴とする請求項１から７のいずれか一項に記載の画像処理装置。
10. 前記取得手段で取得したテクスチャ判定情報に基づいて、前記着目画素の値を、前記新しい値として出力するか、前記重み付き平均値の値として出力するかを判定する判定手段をさらに備えることを特徴とする請求項９に記載の画像処理装置。
11. 画像データにおいて画素ごとにノイズを低減する画像処理装置であって、
    前記画像データにおける着目画素と該着目画素近傍の参照画素との類似度に基づいて前記参照画素に対する重みを決定する決定手段と、
    前記重みに基づいて、前記参照画素の画素値を重み付き平均により前記着目画素に対応する重み付き平均値を算出する算出手段と、
    前記着目画素の画素値と前記重み付き平均値との相関関係に基づいて、前記着目画素の出力値を導出する導出手段とを有することを特徴とする画像処理装置。
12. 着目画素に対する複数の参照画素の画素値の重み付き平均値と、前記着目画素の画素値と、前記着目画素の真の画素値と、の間に成り立っている相関関係に基づいて前記着目画素に対する新たな値を導出する導出ステップを備えることを特徴とする画像処理方法。
13. 画像データにおいて画素ごとにノイズを低減する画像処理方法であって、
    前記画像データにおける着目画素と該着目画素近傍の参照画素との類似度に基づいて前記参照画素に対する重みを決定する決定ステップと、
    前記重みに基づいて、前記参照画素の画素値を重み付き平均により前記着目画素に対応する重み付き平均値を算出する算出ステップと、
    前記着目画素の画素値と前記重み付き平均値との相関関係に基づいて、前記着目画素の出力値を導出する導出ステップとを有することを特徴とする画像処理方法。
14. コンピュータを、請求項１から１１のいずれか一項に記載の画像処理装置として機能させるためのプログラム。

Images