(57)【要約】【課題】本発明は画像ノイズ除去方法及び画像ノイズ除去装置を提供する。【解決手段】画像ノイズ除去方法は、取得された処理待ち画像に対して初期ノイズ除去処理を行い、初期ノイズ除去画像を取得し；処理待ち画像及び初期ノイズ除去画像にそれぞれ対応する所定エネルギーパラメータの値に基づいて、処理待ち画像上の各単位領域の中心画素に対応する残差量を計算し；及び、該残差量を用いて各単位領域に対応する重み行列を計算し、そして、該重み行列を用いて処理待ち画像に対して非局所平均化計算を行い、処理待ち画像へのノイズ除去処理を実現することを含む。【選択図】図１
2015225665.tif 000002
【特許請求の範囲】
【請求項１】
画像ノイズ除去方法であって、
取得された処理待ち画像に対して初期ノイズ除去処理を行い、初期ノイズ除去画像を取得し；
前記処理待ち画像及び前記初期ノイズ除去画像にそれぞれ対応する所定エネルギーパラメータの値に基づいて、前記処理待ち画像上の各単位領域の中心画素に対応する残差量を計算し；及び
前記残差量を用いて前記各単位領域に対応する重み行列を計算し、前記重み行列に基づいて前記処理待ち画像に対して非局所平均化計算を行うことにより、前記処理待ち画像に対してノイズ除去処理を行うことを含む、方法。
【請求項２】
請求項1に記載の方法であって、
前記重み行列を計算することは、
前記処理待ち画像上の任意の単位領域を選択し、前記任意の単位領域の、前記処理待ち画像における関連領域を取得し；及び
前記関連領域中の各関連単位領域と、前記任意の単位領域との間の距離及び前記残差量を基づいて、前記各関連単位領域に対応する重み値を計算して前記重み行列を得ることを含む、方法。
【請求項３】
請求項2に記載の方法であって、
【数１】
2015225665.tif 000019
を用いて、前記関連領域中の任意の関連単位領域に対応する重み値を計算し、
そのうち、w(n,m)は、任意の関連単位領域(n,m)に対応する重み値を示し、d(n,m)は、前記任意の関連単位領域と、前記任意の単位領域との間の距離を示し、residuals(n,m)は、前記任意の単位領域の中心画素に対応する残差量を示し、hは、所定制御係数を示す、方法。
【請求項４】
請求項3に記載の方法であって、
【数２】
2015225665.tif 000020
を用いて、前記任意の関連単位領域(n,m)と、前記任意の単位領域(i,j)との距離を計算し、
そのうち、2r+1は、前記任意の関連単位領域の辺長を示し、t、k及びTは、中間値を示す、方法。
【請求項５】
請求項2〜4のうちの任意の１つに記載の方法であって、
前記非局所平均化計算は、
前記関連領域中の各関連単位領域と、前記重み行列中の重み値との間の対応関係に基づいて、前記関連領域中の全ての画素の加重和を、前記任意の単位領域の中心画素に対してノイズ除去処理を行った後の結果として計算することを含む、方法。
【請求項６】
画像ノイズ除去装置であって、
取得された処理待ち画像に対して初期ノイズ除去処理を行い、初期ノイズ除去画像を得るための初期ノイズ除去ユニット；
前記処理待ち画像及び前記初期ノイズ除去画像にそれぞれ対応する所定エネルギーパラメータの値に基づいて、前記処理待ち画像上の各単位領域の中心画素に対応する残差量を計算するための残差量計算ユニット；及び
前記残差量を用いて前記各単位領域に対応する重み行列を計算し、前記重み行列に基づいて前記処理待ち画像に対して非局所平均化計算を行うことにより、前記処理待ち画像に対してノイズ除去処理を行うための非局所平均値計算ユニットを含む、装置。
【請求項７】
請求項6に記載の装置であって、
前記非局所平均値計算ユニットは、
前記処理待ち画像上の任意の単位領域を選択し、前記任意の単位領域の、前記処理待ち画像における関連領域を確定するための関連領域確定ユニット；及び
前記関連領域中の各関連単位領域と、前記任意の単位領域との間の距離及び前記残差量に基づいて、前記各関連単位領域に対応する重み値を計算して前記重み行列を得るための重み行列取得ユニットを含む、装置。
【請求項８】
請求項7に記載の装置であって、
前記重み行列取得ユニットは、
【数３】
2015225665.tif 000021
を用いて、前記関連領域中の任意の関連単位領域に対応する重み値を計算し、
そのうち、w(n,m)は、任意の関連単位領域(n,m)に対応する重み値を示し、d(n,m)は、前記任意の関連単位領域と、前記任意の単位領域との間の距離を示し、residuals(n,m)は、前記任意の単位領域の中心画素に対応する残差量を示し、hは、所定制御係数を示す、装置。
【請求項９】
請求項8に記載の装置であって、
【数４】
2015225665.tif 000022
を用いて、前記任意の関連単位領域(n,m)と、前記任意の単位領域(i,j)との間の距離を計算し、
そのうち、2r+1は、前記任意の関連単位領域の辺長を示し、t、k及びTは、中間値を示す、装置。
【請求項１０】
請求項7〜9のうちの任意の１つに記載の装置であって、
前記非局所平均値計算ユニットは、
前記関連領域中の各関連単位領域と、前記重み行列中の重み値との間の対応関係に基づいて、前記関連領域中の全ての画素の加重和を、前記任意の単位領域の中心画素に対してノイズ除去処理を行った後の結果として計算するための加重和計算ユニットを含む、装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像処理技術分野に関し、具体的に、画像ノイズ除去方法及び画像ノイズ除去装置に関する。
【背景技術】
【０００２】
デジタルカメラ及びスキャナの発展及び普及に伴い、デジタル画像を非常に容易に得ることができるようになる。しかし、実際に用いられるセンサは、そのすべてが、ノイズ、例えば、熱、電気エネルギー又は他の物質によるノイズにより一定程度干渉され得る。ノイズが信号の測定にひどく影響し得るので、得られた画像データは、信号及びノイズの両方を含む場合がある。画像に関連する様々なアプリケーション、例えば、医用画像解析、画像分割、及び対象物検出などは、基本的には、有効なノイズ抑制方法を用いて、より信頼できる結果を生成することを要する。よって、画像処理及びコンピュータビジョンでは、画像フィルタリングが、既に幅広く使われている重要な技術の１つになっている。そのうち、グレースケール画像が一種の重要且つ常用の画像とされているため、そのノイズのフィルタリング方法は、非常に重要である。
【０００３】
今のところ、グレースケール画像のノイズをフィルタリングするためのいくつかのノイズ除去方法が存在する。例えば、しきい値に基づくウェーブレット解析法、非局所平均化フィルタリング法、Gaussianフィルタリング法、及びバイラテラルフィルタリング法がある。
【０００４】
しきい値に基づくウェーブレット解析法は、ウェーブレット変換を原始画像に対して行うことで、該原始画像をウェーブレット域に転換し、そして、しきい値法を用いてマルチバンドのウェーブレット係数をフィルタリングし、一般的には、分解された第一級中の対角線詳細係数、水平詳細係数、及び垂直詳細係数をフィルタリングする。そのうち、一種のしきい値法は、ハードしきい値函数によるものであり、即ち、しきい値よりも小さい全ての詳細係数をセロに設定し、残りの詳細係数を保留するとのことである。最後に、しきい値法による処理を行った後に、ウェーブレット逆変換により、全てのウェーブレット係数を画像域に戻す。ところが、しきい値に基づくウェーブレット解析法により、ノイズを抑えることができるが、その同時に画像の幾つかの詳細も抑えらてしまう。
【０００５】
非局所平均化フィルタリング法は、一種の非線形のエッジ保護フィルタリング法であり、それは、各出力画素を入力画素の加重和として計算する。１つの出力画素に貢献する一組の入力画素が入力画像における１つの大領域からものであるため、非局所と称される。これに対して、局所平均化フィルタリング法の１つのキー特徴は、小画像ブロックの間の距離に基づいて重みを決定することである。非局所平均化フィルタリング法により、画像の詳細を保留しながら高周波数のGaussianノイズを抑制することができる。ところが、このフィルタリング法は、ノイズが多く且つそのうちのGaussianノイズの割合が小さいリアルなシナリオに適しない。
【０００６】
Gaussianフィルタリング法は、一種の加重平均法である。各出力画素は、該画素の隣接画素の加重平均として設定され、原始画素の輝度値は、最大の重みを取り、隣接画素は、該原始画素との距離の増加に基づいて比較的小さい重みを取る。ところが、このような画像フィルタリングを行った後に、ノイズが減らされ得ると同時に、画像の詳細も減らされてしまう。
【０００７】
バイラテラルフィルタリング法は、一種のエッジ保護及びノイズ除去のフィルタリング法である。画像中の各画素の輝度値は、その周辺の画素の輝度値の加重平均により置換される。この重みは、Gaussian分布に基づくものである。キーとなるのは、この重みがユークリッド距離のみならず、輻射度の差にも頼りにする。それは、系統的に各画素を遍歴し、対応する隣接画素に重みを与えることで、明らかなエッジを保留することができる。ところが、画像がノイズの影響をひどく受けている場合、この方法により、エッジが偽物又は不自然に見えるようにさせてしまう。
【発明の概要】
【発明が解決しようとする課題】
【０００８】
上述に鑑みて、本発明は、少なくとも、従来のノイズ除去技術によりノイズ除去後の画像が自然に見えるようにさせることができないといった問題を解決するための画像ノイズ除去方法、画像ノイズ除去装置、及び電子装置を提供することを課題とする。
【課題を解決するための手段】
【０００９】
本発明の一側面によれば、画像ノイズ除去方法が提供され、それは、取得された処理待ち画像に対して初期ノイズ除去処理を行い、初期ノイズ除去画像を取得し；前記処理待ち画像及び前記初期ノイズ除去画像にそれぞれ対応する所定エネルギーパラメータの値に基づいて、前記処理待ち画像上の各単位領域の中心画素に対応する残差量を計算し；及び、前記残差量を用いて、前記各単位領域に対応する重み行列を計算し、そして、前記重み行列に基づいて、前記処理待ち画像に対して非局所平均化計算を行い、前記処理待ち画像に対するノイズ除去処理を実現することを含む。
【００１０】
本発明の他の側面によれば、画像ノイズ除去装置が提供され、それは、取得された処理待ち画像に対して初期ノイズ除去処理を行い、初期ノイズ除去画像を得るための初期ノイズ除去ユニット；前記処理待ち画像及び前記初期ノイズ除去画像にそれぞれ対応する所定エネルギーパラメータの値に基づいて、前記処理待ち画像上の各単位領域の中心像素に対応する残差量を計算するための残差量計算ユニット；及び、前記残差量を用いて、前記各単位領域に対応する重み行列を計算し、そして、前記重み行列に基づいて、前記処理待ち画像に対して非局所平均化計算を行い、前記処理待ち画像に対するノイズ除去処理を実現するための非局所平均値計算ユニットを含む。
【００１１】
本発明の他の側面によれば、電子装置が提供され、それは、上述のような画像ノイズ除去装置を含む。
【００１２】
本発明の他の側面によれば、コンピュータ読み取り可能な指令コードを含むプログラムが提供され、それは、実行される時に、コンピュータに、上述のような画像ノイズ除去方法を実行させることができる。
【００１３】
また、本発明の他の側面によれば、コンピュータ読み取り可能な記憶媒体が提供され、それは、上述のようなプログラムを記憶している。
【００１４】
上述のような画像ノイズ除去方法、画像ノイズ除去装置、及び電子装置によれば、各単位領域の中心画素の残差量を用いて各単位領域の重み行列を計算し、そして、この重み行列に基づいて処理待ち画像に対して非局所平均化計算を行うことで、少なくとも、次のような有益な効果の１つを実現することができ、即ち、画像の詳細な特徴をできるだけ保留できること、及び、ノイズ除去後の画像がより自然に見えるようになることである。
【図面の簡単な説明】
【００１５】
【図１】本発明の一実施例における画像ノイズ除去方法のフローチャートである。
【図２】本発明の他の実施例における画像ノイズ除去方法のフローチャートである。
【図３】本発明の一実施例における各画素の周辺領域を示す図である。
【図４】本発明の一実施例における原始画像に対しての2D-DCT変換を示す図である。
【図５】本発明の一実施例における2D-DCT域での２次元Wienerフィルタリングを示す図である。
【図６】本発明の一実施例における原始画像の初期ノイズ除去後の様子を示す図である。
【図７】本発明の一実施例における関連単位領域の重み行列の取得を示す図である。
【図８】本発明の一実施例における画像ノイズ除去装置のブロック図である。
【発明を実施するための形態】
【００１６】
以下、添付した図面を参照しながら本発明を実施するための形態を詳しく説明する。
【００１７】
図1は、本発明の一実施例における画像ノイズ除去方法のフローチャートである。
【００１８】
図1に示すように、本発明の実施例による画像ノイズ除去方法は、次のステップを含んでもよい。
【００１９】
ステップ102：取得された処理待ち画像（原始画像又は第１の画像とも言う）に対して初期ノイズ除去処理を行うことで、初期ノイズ除去画像（第２の画像とも言う）を取得し；
ステップ104：処理待ち画像及び初期ノイズ除去画像にそれぞれ対応する所定エネルギーパラメータの値に基づいて、処理待ち画像上の各単位領域の中心画素に対応する残差量を計算し；
ステップ106：残差量を用いて各単位領域に対応する重み行列を計算し、そして、重み行列に基づいて処理待ち画像に対して非局所平均化計算を行い、処理待ち画像に対するノイズ除去処理を実現する。
【００２０】
本実施例では、初期ノイズ除去画像及び原始画像を用いて残差量を得ることができ、該残差量は、ノイズ除去前後の画像情報の差であると理解されても良い。該残差量に基づいて単位領域の重み行列を計算し、そして、該重み行列に基づいて原始画像に対して非局所平均化計算を行うことで、現有技術において原始画像に対してノイズ除去を直接行うことにより生じた、画像詳細情報が無視される及びエッジに歪みが出るといった問題を避けることができる。
【００２１】
以下、図2を参照しながら本発明の他の実施例における画像ノイズ除去方法をより詳細に説明する。
【００２２】
図2に示すように、ステップ202では、選定された画素の周辺（neighboring）領域に対して２次元離散コサイン変換（2D-DCT）を行う。
【００２３】
画素x（i、j）について、図４に示すように、その周辺領域（図3に示すような周辺領域）を離散コサイン変換域に変換し、その後、後述の２次元Wienerフィルタリングを行う。
【００２４】
２次元離散コサイン変換及びその逆変換のプロセスは、式1.1.1乃至式1.1.3を参照することができる。
【数１】
2015225665.tif 000003
【００２５】
そのうち、x,u=0,1,...,M-1であり、y,v=0,1,...,N-1である。
【数２】
2015225665.tif 000004
【００２６】
そのうち、M及びNは、それぞれ、原始画像の高さ及び幅である。
【００２７】
ステップ204では、原始画像に対して２次元Wienerフィルタリング処理を行う。
【００２８】
本実施例では、２次元Wienerフィルタが、各画素の局所周辺領域に基づいて推定された統計データによる自己適応Wiener法を採用して、次の式により局所平均値mean(i,j)及び分散値（variance）var(i,j)を計算する。
【数３】
2015225665.tif 000005
【００２９】
そのうち、（k,m）は、N×Mの大小（サイズ）の周辺領域中の任意のアドレスである。
【００３０】
その後、式（1.2.3）に基づいて画像中の画素に対して２次元Wienerフィルタリングを行う。
【数４】
2015225665.tif 000006
【００３１】
そのうち、nvは、ノイズの分散である。
【００３２】
ステップ206では、式（1.3.1）を採用して画像の２次元離散コサイン変換域において二次元Wienerフィルタリングを行う。
【数５】
2015225665.tif 000007
【００３３】
そのうち、d(u,v)は、D(u,v)（２次元離散コサイン変換域）中の１つの画素であり、２次元離散コサイン変換域の２次元Wienerフィルタリング前後の対比は、図5に示されている。
【００３４】
２次元離散コサイン変換域中の各画素に対して２次元Wienerフィルタリングを行った後に、図６に示すように、初期ノイズ除去画像を得ることができる。
【００３５】
ステップ208では、原始画像及び初期ノイズ除去画像に基づいて、画素x（i,j）の残差量を計算する。
【００３６】
具体的には、残差量は、原始画像のグレースケール値と、初期去除画像のグレースケール値との差の絶対値である。
【００３７】
ステップ210では、改良された非局所平均化フィルタリング法を用いて残差量に対して処理を行い、各像素の重み行列を計算する。なお、本発明の実施例による、改進された非局所平均化フィルタリング法については、後述する。
【００３８】
各出力画素は、入力画素の加重和と定義され、次の式（2.2.1）に示すようである。
【数６】
2015225665.tif 000008
【００３９】
そのうち、x(n,m)は、画素（n,m）の輝度値であり、w(n,m)は、輝度値x(n,m)に対応する重みであり、w(n,m)の値は、画像領域の間の距離及び残差量により決定される。例えば、図7中の画素x(n,m)の重みは、残留量residuals(n,m)及び現在の画像ブロック（任意の単位領域）と、その参考画像ブロック（関連領域中の複数の関連単位領域）と間の距離に基づいて決定される。次の式は、重みw(n、m)の計算過程を示している。
【数７】
2015225665.tif 000009
【００４０】
ここでは、３つの係数、即ち、捜索半径R、関連単位領域の半径r、制御係数hを設定する。捜索半径Ｒは、加重和を計算するための画素の範囲、例えば、R×Rを意味する。関連単位領域（参考画像ブロックとも言う）の半径rは、図７に示すように、関連単位領域の半分のサイズを意味する。制御係数hは、ノイズ除去効果を制御するためのものであり、例えば、より多くのノイズを抑制する必要がある場合、制御係数hの値を上げてもよい。なお、この制御係数の範囲は、1以上且つ100以下である。
【００４１】
式（2.2.4）から分かるように、残差量が大きい場合、該画素のノイズが比較的大きいとのことを示し、この場合、その対応する重みが比較的小さい（それが採用される程度が小さい）。重みが小さくなった後に、ノイズの抑制が大きくなる。よって、同一の制御係数hの条件の下で、重みの異なる値が、異なる程度のノイズ抑制を表し、また、制御係数hの値を変更することで、ノイズ抑制の程度を変更することもできる。
【００４２】
式（2.2.2）及び式（2.2.1）による処理を行った後に、各画素（i,j）について、次のように示すR×Rの重み行列を得ることができる。
【数８】
2015225665.tif 000010
【００４３】
画素（i,j）の重みが最大値であるため、次のように、画素（i,j）の重み値の再設定及び重み行列の正規化を行う。
【数９】
2015225665.tif 000011
【００４４】
ステップ212では、関連領域中の各関連単位領域と、行列中の重み値との間の対応関係に基づいて、前記関連領域中の全ての画素の加重和を、任意の単位領域の中心画素に対してノイズ除去処理を行った後の結果として計算する。
【００４５】
このステップでは、R×Rの重み行列が原始画像の加重和を計算するために用いられる。原始画像の画素（i,j）は、OI（i,j）と定義される。
【００４６】
重み行列に基づいて原始画像（即ち、処理待ち画像）に対して非局所平均化計算を計算し、これにより、式（2.2.1）に示すように、画像の最後のフィルタリング結果NLM-OI(i,j)を得ることができる。
【数１０】
2015225665.tif 000012
【００４７】
図8は、本発明の実施例における画像ノイズ除去装置のブロック図である。
【００４８】
図8に示すように、本発明の実施例における画像ノイズ除去装置800は、次のようなものを含んでもよい。
【００４９】
初期ノイズ除去ユニット802：取得された処理待ち画像に対して初期ノイズ除去処理を行い、初期ノイズ除去画像を取得し；
残差量計算ユニット804：処理待ち画像及びノイズ除去画像にそれぞれ対応する所定エネルギーパラメータの値に基づいて、前記処理待ち画像上の各単位領域の中心画素に対応する残差量を計算し；
非局所平均値計算ユニット806：残差量を用いて、各単位領域に対応する重み行列を計算し、そして、重み行列に基づいて、処理待ち画像に対して非局所平均化計算を行い、処理待ち画像へのノイズ除去処理を実現する。
【００５０】
そのうち、非局所平均値計算ユニット806は、次のようなものを含む。
【００５１】
関連領域確定ユニット8062：処理待ち画像上の任意の１つの単位領域を選択し、そして、前記任意の１つの単位領域の、処理待ち画像における関連領域を確定する。
【００５２】
重み行列取得ユニット8064：関連領域中の各関連単位領域と、前記任意の１つの単位領域との距離値及び前記残差量に基づいて、各関連単位領域に対応する重み値を計算し、重み行列を得る。
【００５３】
重み行列取得ユニット8064は、次の式を用いて、前記関連領域中の任意の１つの関連単位領域に対応する重み値を計算する。
【数１１】
2015225665.tif 000013
【００５４】
そのうち、w(n,m)は、任意の１つの関連単位領域(n,m)に対応する重み値であり、d(n,m)は、前記任意の１つの関連単位領域と、前記任意の１つの単位領域と間の距離値であり、residuals(n,m)は、前記任意の１つの単位領域の中心画素に対する残差量であり、hは、所定制御係数である。
【００５５】
重み行列取得ユニット8064は、次の式を用いて前記任意の関連単位領域(n,m)と、前記任意の単位領域(i,j)との間の距離の値を計算する。
【数１２】
2015225665.tif 000014
【００５６】
そのうち、2r+1は、前記任意の１つの関連単位領域の辺長であり、t、k及びTは、中間値である。
【００５７】
非局所平均値計算ユニット806は、次のものを含んでも良い。
【００５８】
加重和計算ユニット8066：前記関連領域中の各関連単位領域と、前記重み行列中の重み値との対応関係に基づいて、前記関連領域中の全ての画素の加重和を、前記任意の１つの単位領域の中心画素に対してノイズ除去処理を行った後の結果として計算する。
【００５９】
画像ノイズ除去装置はさらに、変換ユニット808を含んでもよく、それは、前記初期ノイズ除去処理を行ってもらうために、前記処理待ち画像を２次元離散コサイン変換域に変換する。即ち、変換ユニット808は、前記初期ノイズ除去処理を行う前に、前記処理待ち画像を２次元離散コサイン変換域に変換する。そのうち、前記初期ノイズ除去ユニット802は、Wienerフィルタリング・アルゴリズムに基づいて前記初期ノイズ除去処理を行う。
【００６０】
本発明の画像ノイズ除去装置によれば、一種の改進された局所非平均化ノイズ除去方法を提供し、そして、初期ノイズ除去画像と原始画像との間の残差量に基づいて、原始画像に対してノイズ除去を行い、画像のより多くの詳細的な情報を保留することができ、これにより、ノイズ除去後の画像がより自然に見えるようになる。
【００６１】
また、本発明の実施例はさらに電子装置を提供し、この電子装置は、上述のような画像ノイズ除去装置を含む。本発明の実施例による上述の電子装置の具体的な実現方式の中では、上述の電子装置は、次のような装置のうちの任意の１つであってもよく、即ち、コンピュータ、タブレットＰＣ、ＰＤＡ、マルチメディアプレイヤー、携帯電話、及び電子ブックなどである。そのうち、この電子装置は、上述の画像ノイズ除去装置の各種の機能及び技術的効果を有するため、ここでは、その詳しい説明を省略する。
【００６２】
上述の本発明の実施例による画像ノイズ除去装置中の各構成ユニット、サブユニット、モジュールなどは、ソフトウェア、ファームウェア、ハードウェア又はその任意の組み合わせの方式で構成されてもよい。ソフトウェア又はファームウェアで実現する場合、記憶媒体又はネットワークから、専用ハードウェア構造を有するマシン（例えば、コンピュータ）に、該ソフトウェア又はファームウェアを構成するプログラムをインストールしてもよく、また、該マシンは、各種のプログラムがインストールされているときに、上述の各構成ユニット、サブユニットの各種の機能を実行することができる。
【００６３】
また、本発明はさらに、マシン（コンピュータ）可読な指令コードを含むプログラムに関する。
【００６４】
上述の指令コードは、マシンにより読み取りされて実行される時に、上述の本発明の実施例による画像イズ除去方法を実行することができる。それ相応に、このようなプログラムを記憶した、例えば磁気ディスク、光ディスク、磁気光ディスク、半導体記憶器などの各種の記憶媒体も、本発明に含まれている。
【００６５】
また、本発明の各実施例における方法は、明細書に述べた又は図面に示した時間順序で実行されるに限定されず、他の順序で、並列に又は独立して実行されてもよい。よって、本明細書に記載された方法の実行順序は、本発明の技術的範囲を限定しない。
【００６６】
また、本発明の上述の方法における各処理ステップは、各種のマシン可読な記憶媒体に記憶されているコンピュータ実行可能なプログラムで実現することもできる。
【００６７】
また、本発明の目的は、次のような方式で実現されてもよく、即ち、上述のコンピュータ実行可能なプログラムコードを記憶した記憶媒体が、直接又は間接的にシステム又は装置に提供され、そして、該システム又は装置中のコンピュータ又は中央処理ユニット（CPU）が上述のプログラムコードを読み出して実行するといった方式で実現されてもよい。
【００６８】
このとき、該システム又は装置は、プログラムを実行する機能を有すれば、本発明の実施形態は、プログラムに限定されず、また、該プログラムは、任意の形式、例えば、オブジェクト指向プログラム、インタープリタ実行可能なプログラム、又は、ＯＳ（operating  system）に提供されるスクリプト・プログラムなどであってもよい。
【００６９】
上述のようなマシン可読な記憶媒体は、各種の記憶器、記憶ユニット、半導体装置、ディスクユニット（例えば光、磁気及び磁気光ディスク）、及び、情報を記憶可能な他の媒体を含んでもよいが、これらに限定されない。
【００７０】
また、上述の各実施例を含む実施形態に関し、更に以下の付記を開示する。
【００７１】
（付記1）
画像ノイズ除去方法であって、
取得された処理待ち画像に対して初期ノイズ除去処理を行い、初期ノイズ除去画像を取得し；
前記処理待ち画像及び前記初期ノイズ除去画像にそれぞれ対応する所定エネルギーパラメータの値に基づいて、前記処理待ち画像上の各単位領域の中心画素に対応する残差量を計算し；及び
前記残差量を用いて、前記各単位領域に対応する重み行列を計算し、そして、前記重み行列に基づいて、前記処理待ち画像に対して非局所平均化計算を行い、前記処理待ち画像へのノイズ除去処理を実現することを含む、方法。
【００７２】
（付記2）
付記1に記載の画像ノイズ除去方法であって、
前記重み行列を計算することは、
前記処理待ち画像上の任意の単位領域を選択し、そして、前記任意の単位領域の、前記処理待ち画像における関連領域を確定し；及び
前記関連領域中の各関連単位領域と、前記任意の単位領域との間の距離及び前記残差量に基づいて、前記各関連単位領域に対応する重み値を計算して前記重み行列を得ることを含む、方法。
【００７３】
（付記3）
付記2に記載の画像ノイズ除去方法であって、
【数１３】
2015225665.tif 000015
を用いて、前記関連領域中の任意の関連単位領域の重み値を計算し、
そのうち、w(n,m)は、任意の関連単位領域(n、m)に対応する重み値であり、d(n,m)は、前記任意の関連単位領域と、前記任意の単位領域との間の距離であり、residuals(n,m)は、前記任意の単位領域の中心画素に対応する残差量であり、hは、所定制御係数である、方法。
【００７４】
（付記4）
付記3に記載の画像ノイズ除去方法であって、
【数１４】
2015225665.tif 000016
を用いて、前記任意の関連単位領域(n,m)と、前記任意の単位領域(i,j)との間の距離を計算し、
そのうち、2r+1は、前記任意の関連単位領域の辺長であり、t、k及びTは、中間値である、方法。
【００７５】
（付記5）
付記2〜4のうちの任意の１つに記載の画像ノイズ除去方法であって、
前記非局所平均化計算は、
前記関連領域中の各関連単位領域と、前記重み行列中の重み値との間の対応関係に基づいて、前記関連領域中の全ての画素の加重和を、前記任意の単位領域の中心画素に対してノイズ除去処理を行った後の結果として計算すことを含む、方法。
【００７６】
（付記6）
付記1に記載の画像ノイズ除去方法であって、
前記所定エネルギーパラメータは、グレースケール値を含む、方法。
【００７７】
（付記7）
付記1〜6のうちの任意の１つに記載の画像ノイズ除去方法であって、さらに、
前記処理待ち画像を２次元離散コサイン変換域に変換してから、前記初期ノイズ除去処理を行うことを含む、方法。
【００７８】
（付記8）
付記7に記載の画像ノイズ除去方法であって、
Wienerフィルタリング・アルゴリズムに基づいて、前記初期ノイズ除去処理を行う、方法。
【００７９】
（付記9）
画像ノイズ除去装置であって、
取得された処理待ち画像に対して初期ノイズ除去処理を行い、初期ノイズ除去画像を得るための初期ノイズ除去ユニット；
前記処理待ち画像及び前記初期ノイズ除去画像にそれぞれ対応する所定エネルギーパラメータの値に基づいて、前記処理待ち画像上の各単位領域の中心画素に対応する残差量を計算するための残差量計算ユニット；及び
前記残差量を用いて、前記各単位領域に対応する重み行列を計算し、そして、前記重み行列に基づいて、前記処理待ち画像に対して非局所平均化計算を行い、前記処理待ち画像へのノイズ除去処理を実現するための非局所平均値計算ユニットを含む、装置。
【００８０】
（付記10）
付記9に記載の画像ノイズ除去装置であって、
前記非局所平均値計算ユニットは、
前記処理待ち画像上の任意の単位領域を選択し、そして、前記任意の単位領域の、前記処理待ち画像における関連領域を得るための関連領域確定ユニット；及び
前記関連領域中の各関連単位領域と、前記任意の単位領域と間の距離及び前記残差量に基づいて、前記各関連単位領域に対応する重み値を計算して前記重み行列を得るための重み行列取得ユニットを含む、装置。
【００８１】
（付記11）
付記10に記載の画像ノイズ除去装置であって、
前記重み行列取得ユニットは、
【数１５】
2015225665.tif 000017
を用いて、前記関連領域中の任意の関連単位領域に対応する重み値を計算し、
そのうち、w(n,m)は、任意の関連単位領域(n,m)に対応する重み値であり、d(n,m)は、前記任意の関連単位領域と、前記任意の単位領域との間の距離をであり、residuals(n,m)は、前記任意の単位領域の中心画素に対応する残差量であり、hは、所定制御係数である、装置。
【００８２】
（付記12）
付記11に記載の画像ノイズ除去装置であって、
【数１６】
2015225665.tif 000018
を用いて、前記任意の関連単位領域(n,m)と、前記任意の単位領域(i,j)との間の距離を計算し、
そのうち、2r+1は、前記任意の関連単位領域の辺長であり、t、k及びTは、中間値である、装置。
【００８３】
（付記13）
付記10に記載の画像ノイズ除去装置であって、
前記非局所平均値計算ユニットは、
前記関連領域中の各関連単位領域と、前記重み行列中の重み値との間の対応関係に基づいて、前記関連領域中の全ての画素の加重和を、前記任意の単位領域の中心画素に対してノイズ除去処理を行った後の結果として計算する加重和計算ユニットを含む、装置。
【００８４】
（付記14）
付記9に記載の画像ノイズ除去装置であって、
前記所定エネルギーパラメータは、グレースケール値を含む、装置。
【００８５】
（付記15）
付記9〜14のうちの任意の１つに記載の画像ノイズ除去装置であって、さらに、
前記初期ノイズ除去処理を行う前に、前記処理待ち画像を２次元離散コサイン変換域に変換するための変換ユニットを含む、装置。
【００８６】
（付記16）
付記15に記載の画像ノイズ除去装置であって、
前記初期ノイズ除去ユニットは、Wienerフィルタリングアルゴリズムに基づいて、前記初期ノイズ除去処理を行う、装置。
【００８７】
（付記17）
電子装置であって、
前記電子装置は、付記9〜16のうちの任意の１つに記載の画像ノイズ除去装置を含む、で電子装置。
【００８８】
以上、本発明の好ましい実施形態を説明したが、本発明はこのような実施形態に限定されず、本発明の趣旨を離脱しない限り、本発明に対するあらゆる変更は本発明の技術的範囲に属する。
【図１】
2015225665.tif 000023
【図２】
2015225665.tif 000024
【図３】
2015225665.tif 000025
【図４】
2015225665.tif 000026
【図５】
2015225665.tif 000027
【図６】
2015225665.tif 000028
【図７】
2015225665.tif 000029
【図８】
2015225665.tif 000030