ノイズ低減装置、ノイズ低減方法、ノイズ低減プログラム、記録媒体
【課題】エッジなど人間が画像を知覚する上で重要な要素を保存しつつノイズ低減を行い、カラー画像における非加法性のノイズを適切に低減する、ノイズ低減理装置等を提供すること。
【解決手段】カラー画像の注目画素の画素値に、注目画素を含む近傍領域の画素の画素値を用いてノイズ低減処理を施すノイズ低減装置100であって、所定の色空間において注目画素との類似度が閾値以上の類似領域に含まれる画素を近傍領域から抽出し、類似領域の画素の代表値ベクトルを算出する代表値ベクトル算出手段103と、注目画素の注目画素ベクトルと前記代表値ベクトルとの差分ベクトルを、色空間の特定方向に射影する差分射影手段105と、代表値ベクトルに、特定方向に射影された前記差分ベクトルを加えたベクトルの要素を、注目画素の画素値に置き換える画素値修正手段106と、を有することを特徴とする。
【解決手段】カラー画像の注目画素の画素値に、注目画素を含む近傍領域の画素の画素値を用いてノイズ低減処理を施すノイズ低減装置100であって、所定の色空間において注目画素との類似度が閾値以上の類似領域に含まれる画素を近傍領域から抽出し、類似領域の画素の代表値ベクトルを算出する代表値ベクトル算出手段103と、注目画素の注目画素ベクトルと前記代表値ベクトルとの差分ベクトルを、色空間の特定方向に射影する差分射影手段105と、代表値ベクトルに、特定方向に射影された前記差分ベクトルを加えたベクトルの要素を、注目画素の画素値に置き換える画素値修正手段106と、を有することを特徴とする。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、カラー画像のノイズを低減するノイズ低減理装置、ノイズ低減方法、プログラム及び記憶媒体に関する。
【背景技術】
【0002】
スキャナやデジタルカメラ等の撮像機器により撮影した画像には、撮像素子及び回路の特性上、暗電流ノイズ、熱雑音、ショットノイズなどのノイズが含まれる。高画質の画像を得るためには、これらのノイズを低減する処理を施さなければならない。しかし、単純にローパスフィルタを用いてノイズを低減すると、エッジなど人間が画像を知覚する上で重要な要素も同時に失われて画質が劣化する。そこで、画像の領域の特性に応じて、適応的にノイズ低減を行う技術が必要となる。
【0003】
そのようなノイズ低減技術の一つとして、εフィルタが提案されている(例えば、非特許文献1参照。)。εフィルタでは、注目画素との信号差が一定閾値以下となる周辺画素をフィルタ処理に用いることで、エッジなどの信号差の大きい成分を保存しつつノイズを低減することを可能としている。
【0004】
また、バイラテラルフィルタを用いたノイズ低減技術が提案されている(例えば、非特許文献2参照。)。バイラテラルフィルタは、注目画素との信号差及び空間距離によってフィルタの重み係数を作成しフィルタ処理を行う方法であり、εフィルタ同様、エッジの保存とノイズ低減を両立させることができる。なお、バイラテラルフィルタに準ずる原理によるノイズ低減は他にも提案されている(例えば、特許文献1,2参照。)。
しかしながら、εフィルタやバイラテラルフィルタには、背景部分のような振幅の小さい信号成分が失われやすいという問題がある。そこで、この問題を解決する技術が提案されている(例えば、特許文献3参照。)。特許文献3には、εフィルタ出力値と注目画素値を注目画素近傍の空間周波数に応じた割合で加重加算する技術が開示されている。
【0005】
しかしながら、特許文献3記載の方法でも、ショットノイズなど非加法性のノイズに対応できず、また、カラー画像に対するノイズ低減性能が悪い、という問題がある。カラー画像における非加法性ノイズの低減技術が提案されている(例えば、非特許文献3参照。)が、この方法は画像のセグメンテーション、統計量の算出に多大な処理コストがかかるという問題がある。
【発明の概要】
【発明が解決しようとする課題】
【0006】
本発明は上記従来技術の問題を鑑みたものであり、エッジなど人間が画像を知覚する上で重要な要素を保存しつつノイズ低減を行い、さらに振幅の小さい信号を保存し、カラー画像における非加法性のノイズを高精度に少ない処理コストで低減する、ノイズ低減理装置、ノイズ低減方法、プログラム及び記憶媒体を提供することを目的とする。
【課題を解決するための手段】
【0007】
カラー画像の注目画素の画素値に、注目画素を含む近傍領域の画素の画素値を用いてノイズ低減処理を施すノイズ低減装置であって、所定の色空間において注目画素との類似度が閾値以上の類似領域に含まれる画素を前記近傍領域から抽出し、前記類似領域の画素の代表値ベクトルを算出する代表値ベクトル算出手段と、注目画素の注目画素ベクトルと前記代表値ベクトルとの差分ベクトルを、前記色空間の特定方向に射影する差分射影手段と、前記代表値ベクトルに、前記特定方向に射影された前記差分ベクトルを加えたベクトルの要素を、注目画素の画素値に置き換える画素値修正手段と、を有することを特徴とする。
【発明の効果】
【0008】
エッジなど人間が画像を知覚する上で重要な要素を保存しつつノイズ低減を行い、カラー画像における非加法性のノイズを適切に低減する、ノイズ低減理装置、ノイズ低減方法、プログラム及び記憶媒体を提供することができる。
【図面の簡単な説明】
【0009】
【図1】デジタルカメラの概略構成図の一例である。
【図2】ノイズ低減装置の機能ブロック図の一例である。
【図3】ノイズ低減処理の手順を示すフローチャート図の一例である。
【図4】平均、分散の算出を模式的に説明する図の一例である。
【図5】ステップS103の詳細フローの一例を示すフローチャート図である。
【図6】信号値(画素値)に対するノイズの分散をプロットした図の一例である。
【図7】ノイズの分散を推定するための関数の一例を示す図である。
【図8】差分ベクトルsの変換を模式的に説明する図の一例である。
【図9】ステップS105の詳細フローの一例を示すフローチャート図である。
【図10】特定方向eと類似領域の平均値のRGB比率一定方向の関係を示す図である。
【図11】ノイズ低減装置の機能ブロック図の一例である(実施例2)。
【図12】ノイズ低減処理の手順を示すフローチャート図の一例である(実施例2)。
【図13】ノイズ低減装置の機能ブロック図の一例である(実施例3)。
【図14】ノイズ低減処理の手順を示すフローチャート図の一例である(実施例3)。
【図15】ステップS201の詳細フローの一例を示すフローチャート図である。
【図16】ノイズ低減装置の機能ブロック図の一例である(実施例4)。
【図17】ノイズ低減処理の手順を示すフローチャート図の一例である(実施例4)。
【図18】フィルタサイズが5×5の場合のフィルタの一例を示す図である。
【図19】ノイズ低減装置の機能ブロック図の一例である(実施例5)。
【図20】ノイズ低減処理の手順を示すフローチャート図の一例を示す(実施例5)。
【図21】差分ベクトルの二次元方向への射影を模式的に示す図の一例である。
【図22】差分ベクトルの三次元方向への射影を模式的に示す図の一例である。
【発明を実施するための形態】
【0010】
以下、本発明を実施するための形態について図面を参照しながら実施例を挙げて説明する。
【実施例1】
【0011】
本実施例では、一つの画素がRGBの3つの成分を持つカラー画像から、信号値に大きさが依存する非加法性のノイズを低減するノイズ低減装置100について説明する。以下、信号値を画素値という場合がある。本実施例のノイズ低減装置100は、デジタルカメラ、スキャナ、カメラ付き携帯などカラー画像を撮影する機器であれば好適に適用できる。また、コンピュータが記憶媒体から読み出したり、ネットワークを介して取得したカラー画像に本実施例のノイズ低減処理を施すこともできる。
【0012】
図1は、デジタルカメラの概略構成図の一例を示す。デジタルカメラの構成そのものは公知のものを利用できる。デジタルカメラは、例えば、像を撮影する撮像部11と、撮像部11により得られた画像信号に処理を施す信号処理部18と、デジカメ全体を制御する主制御部17、画像データを記憶するフレームメモリ19、及び、その他の部品と接続されるI/F(インタフェース)21、を有する。I/F21 には画像を表示する表示部22と、画像データを保存したり画像データを読み出すための外部メモリ23、及び、記憶媒体24を装着するメモリカード装着部24とが接続されている。ノイズ低減処理は、例えば信号処理部18と主制御部17に実装されるが、実装例は図に限定されるものではない。
【0013】
撮像部11は、レンズ12、絞り13、電子シャッター14、光電変換素子15及び前処理部16を有する。光電変換素子15は、例えば、CMOS(Complementary Metal Oxide Semiconductor)やCCD(Charge Coupled Device)が使用される。撮像部11は不図示のカラーフィルタ(原色系又は補色系のいずれでもよい)を有し、例えばセル状に配列されたR、G、B各色毎に1つの光電変換素子15が配置されている。
【0014】
前処理部16は、プリアンプやAGC(自動利得制御)等のアナログ信号処理やA/D変換器を備えており、光電変換素子15により出力されたアナログ映像信号に対して増幅・クランプ等の前処理を施した後、アナログ映像信号をデジタル映像信号に変換する。
【0015】
信号処理部18は、DSP等を有し、撮像部11から得られたデジタル信号に対して色分解、ホワイトバランス調整、γ補正等の種々の画像処理を施す。信号処理部18は処理した画像データをフレームメモリ19に格納したり、格納した画像データをフレームメモリ19から読み出し、ノイズ低減などの画像処理を施すことができる。フレームメモリ19には、VRAM、SRAM、DRAM等の半導体メモリが使用される。
【0016】
フレームメモリ19から読み出された画像データは、信号処理部18において画像圧縮等の信号処理が施された後、I/F21を介して外部メモリ23又はメモリカード装着部24に装着された記憶媒体25に記憶される。外部メモリ23は、フラッシュメモリ等の不揮発メモリであり、記憶媒体25はUSBメモリ、SDメモリカード、光磁気ディスクなどの脱着可能な可搬型の不揮発メモリである。IF21はこれらを切り替えて利用できる。また、不図示の通信ユニットにより画像データをサーバに送信したり、サーバから画像データを受信して外部メモリ23に記憶してもよい。この場合、通信ユニットは、携帯電話網や無線/有線LAN網に接続して、画像データを送受信する。
【0017】
フレームメモリ19、外部メモリ23又は記憶媒体25から読み出された画像データは、表示部22に表示される。表示部22は、液晶ディスプレイや有機ELディスプレイ等、デジタルカメラの筐体に設置されている。表示部22はタッチパネルを一体に有する。このタッチパネルと、筐体に設けられたキーボードにより、ユーザの操作を入力することができる。
【0018】
主制御部17、又は、主制御部17と信号処理部18は、マイコンやLSIなどにより構成される。したがって、例えば、主制御部、信号処理部及びフレームメモリ19が、ノイズ低減装置100を構成する主要部の一例である。マイコンは、主にCPU、RAM、EEPROM、ASIC、等がバスで接続されたコンピュータを実体とし、EEPROMに記憶されたノイズ低減プログラム20をCPUが実行することで、以下で説明するノイズ低減処理が画像データに施される。また、ASIC等のハードウェアによりノイズ低減処理を実装してもよい。ノイズ低減プログラム20は、予め、EEPROMに記憶された状態でデジタルカメラと共に出荷される。また、記憶媒体25にノイズ低減プログラム20を記憶させて配布し、I/F21を介してノイズ低減プログラム20をEEPROMに読み込むことができる。また、ネットワークを介してノイズ低減プログラム20をEEPROMにダウンロードしてもよい。
【0019】
〔ノイズ低減処理〕
図2はノイズ低減装置100の機能ブロック図の一例を、図3はノイズ低減処理の手順を示すフローチャート図の一例をそれぞれ示す。
【0020】
ノイズ低減装置100は画像データ(以下、処理対象画像31という)を入力とし、ノイズ低減画像32を出力する。ノイズ低減装置100は、画像取得部101、注目画素選択部102、平均・分散算出部103、ノイズ分散推定部104、差分変換部105、加算部106、及び、画素選択済み判定部107を有する。各機能ブロックは、CPUがノイズ低減プログラム20を実行することで実現される。以下、各機能ブロックについて、図3のフローチャート図に沿って説明する。
【0021】
<S101>ステップS101では、画像取得部101がフレームメモリ19から処理対象画像31を取得する。処理対象画像31は、デジタルカメラが直接撮影した以外の画像データでもよい。
【0022】
<S102>ステップS102では、注目画素選択部102が、処理対象画像31の複数の画素の中から注目画素を選択する。たとえばラスタ走査により順次、注目画素を選択すればよく、全ての画素を選択する方法であれば走査方法はどのようなものでもよい。
【0023】
<S103>ステップS103では、平均・分散算出部103が、注目画素の近傍領域における画素の画素値の平均、分散を算出する。ここではεフィルタやバイラテラルフィルタなどの方法を取ることができる。一例として次の方法を示す。
【0024】
図4は、平均、分散の算出を模式的に説明する図の一例である。近傍画素は実際にはカラーである。まず、平均・分散算出部103は、注目画素を中心とした特定ウィンドウ内の近傍画素を選択する。特定ウィンドウの大きさは予め定められていてもよいし、ユーザが設定したノイズ低減処理の強度に応じて可変となっていてもよい。図では(9×9)画素の矩形が特定ウィンドウである。
【0025】
そして、平均・分散算出部103が、注目画素と各近傍画素の類似度を算出する。類似度として、例えばRGB空間でのユークリッド距離の逆数を用いることができる。類似度が特定閾値以上の近傍画素を、平均、分散の演算に用いることとする。こうすることで、注目画素と色合いの近似した近傍画素を抽出し、エッジを保存しながらノイズを低減できる。なお、類似度は画素単位で算出し、平均値はR、G、B毎に算出する。
【0026】
図4では、注目画素を中心とする球内にある画素が、注目画素との類似度が特定閾値以上の近傍画素となる。以下、特定半径の球内を類似領域という。
【0027】
平均・分散算出部103は、次の公式を用いて分散を算出する。
【0028】
【数1】
ここでVar(X)は確率変数Xの分散、Cov(X, Y)は確率変数XとYの共分散、E(X)は確率変数Xの期待値を示す。すなわち、画素値を二乗したデータ、異なる色成分の画素値を乗算したデータ、元の画素値、を組み合わせることで、類似領域の画素の分散・共分散(フィルタ処理)を算出することができる。なお、処理コスト削減の面から共分散は算出せず0と置いてもよい。
【0029】
図5は、ステップS103の詳細フローの一例を示す。xは、近傍画素の画素値であり、RGBの各成分を要素として持つ3次元ベクトルである。平均・分散算出部103は、変数μ、Σ、nに初期値であるゼロを設定する(S1301)。
【0030】
ついで、平均・分散算出部103は、特定ウィンドウの近傍画素の画素値を取得し(S1302)、注目画素との類似度を算出する(S1303)。
【0031】
平均・分散算出部103は、類似度が特定閾値以上か否かを判定する(S1304)。類似度が特定閾値以上の場合(S1304のYes)、その近傍画素が平均、分散の算出対象なので、平均・分散算出部103は、変数μにxを、変数ΣにxxTを、変数nにn+1を、それぞれ加える(S1305)。なお、「T」は転置を表す。
【0032】
ついで、平均・分散算出部103は、特定ウィンドウ内の全ての近傍画素を選択したか否かを判定する(S1306)。特定ウィンドウ内の全ての近傍画素を選択していない場合(S1305のNo)、ステップS1032に戻り上記の処理を繰り返す。
【0033】
特定ウィンドウ内の全ての近傍画素を選択した場合(S1305のYes)、平均・分散算出部103は、「μ/n」、「Σ/n−μμT」をそれぞれ算出する(S1307)。二乗和から平均を引く一般的な分散算出と同様の結果が、少ない計算量で得られる。
【0034】
図5の処理の結果、μが類似領域の平均値ベクトルとして、Σが分散共分散行列として得られる。Σの対角成分は類似領域の分散である。なお、類似領域の平均でなく、中央値や最頻値など種々の代表値を採用することができる。平均等の代表値に、輝度やホワイトバランスを修正する定数を加えたり乗じてもよい。
【0035】
<S103の変形例>
ステップS103の平均、分散の算出方法には種々のバリエーションを取りうる。例えば、実施例1では類似領域の画素に対し単純に平均・分散を求めたが、類似領域の画素の画素値に重み付けをして平均・分散を求めてもよい。バイラテラルフィルタのように、注目画素との類似度と相対位置をフィルタの重み係数に反映させることも有効である。例えば、平均・分散算出部103は、類似度が大きいほど大きな重み係数を画素値に与えて平均・分散を算出する。
【0036】
また、上記実施例ではRGB空間でのユークリッド距離を類似度の算出に用いたが、類似度算出にLnノルム(nは自然数)を用いてもよい。また、色空間には、L*a*b*表色系、XYZ表色系等、適宜選択できる。この場合、平均・分散算出部103は、特定ウィンドウ内の画素の画素値(RGB)を選択した表色系の画素値に変換する。変換公式は既知である。
【0037】
平均・分散算出部103は、このRGB空間とは別の色空間でのノルムを用いることもできる。たとえば、L*a*b*などの等色空間のユークリッドノルムを用いると、人間の知覚に近い色の類似度を得ることができる。
【0038】
<S104>図3に戻りステップS104では、ノイズ分散推定部104が、ノイズの分散を推定する。
図6は信号値(画素値)に対するノイズの分散をプロットした図の一例である。ショットノイズでは信号値とノイズの分散に比例関係が見られる。このように、一般にノイズの分散は画素値に依存する。
【0039】
したがって、画素値からノイズの分散を推定する関数を導出しておくことが可能になる。なお、ノイズの分散はRGBの各成分ごとに求めるので、ノイズの分散を推定する関数もR、G、Bの各成分毎に導出しておく。
【0040】
注目画素のノイズを低減するためには、注目画素に含まれる平均的なノイズを知りたいので、ノイズ分散推定部104は、類似領域のノイズの分散を求める。図7は、類似領域のノイズの分散を推定するための関数の一例を示す。ノイズ分散推定部104は、平均値ベクトルμの各成分に基づき図7の関数を利用して、類似領域のノイズの分散σ2を求める。
【0041】
これを各色成分に適用して各成分のノイズの分散を求める。求めたノイズの分散を対角成分に置いた行列をノイズの分散共分散行列Σnoizeとする。多くの場合はノイズの分散共分散行列の非対角成分(共分散)を0と置くことができる。ノイズの共分散が無視できない場合は、分散推定間数と同様に共分散の推定関数を同定し、それを用いて共分散を推定し分散共分散行列を算出する。
【0042】
<S105>図3に戻りステップS105では、差分変換部105が注目画素gの注目画素ベクトルgから類似領域の平均値ベクトルμを引き、得られた差分(以下、差分ベクトルsという)の方向と大きさを変換する。
【0043】
図8は、差分ベクトルsの変換を模式的に説明する図の一例である。図8(a)に示すように、差分変換部105は、まず、特定方向eに差分ベクトルsを射影する。特定方向eとして、類似領域の平均値ベクトルμの成分(R、G、B)と同じRGB比率となる一定方向(以下、RGB一定比率方向という)を取ることができる。この方向を特定方向eとすることで、注目画素の画素値が類似領域色合いとずれることを防ぐことができる。また、射影することで、差分ベクトルsが、類似領域の色合いの方向に変換され、かつ、長さが低減されるので、類似領域の画素と異なる色合いをもたらすノイズやその大きさを低減する効果が得られる。
【0044】
そして、更に好ましくは、差分変換部105は、類似領域のノイズの分散共分散行列Σnoizeを用いて射影した差分ベクトルsに明示的にノイズを低減する処理を施す。図8(b)に示すように、差分変換部105は、射影したベクトルの大きさに(テクスチャ量)/(ノイズ量+テクスチャ量)の値を乗じる。これにより、ノイズ量及びテクスチャ量に応じてノイズ低減の強度を適応的に制御することができる。
【0045】
なお、ノイズ量はノイズの分散共分散行列Σnoizeのうち特定方向eの分散、テクスチャ量は類似領域の分散共分散行列からノイズの分散共分散行列Σnoizeを引いて得られる分散共分散行列の特定方向eの分散である。一般に分散共分散行列がΣであるとき、特定方向e(単位ベクトル)の分散σe2は次式で求めることができる。
【0046】
【数2】
テクスチャ量は類似領域の分散共分散行列からノイズの分散共分散行列Σnoizeを引いて得られるものであるから、実質的には射影した差分ベクトルsの大きさに(テクスチャ量)/(ノイズ量+テクスチャ量)={1−(ノイズの分散のe方向の大きさ)/(類似領域の分散のe方向の大きさ)}を乗算すればよい。
【0047】
図9は、ステップS105の詳細フローの一例を示す。図9の手順は一例に過ぎず、同様の結果を得られる限り各ステップの内容・順序を入れ替えることができる。注目画素ベクトルgは注目画素の画素値(RGBのベクトル)、Σtotal はステップS1305のΣと同じもので類似領域の分散共分散行列である。
【0048】
まず、差分変換部105は、差分ベクトルsを算出する(S1501)。
s ← g−μ
ついで、差分変換部105は、差分ベクトルsを特定方向eに射影する(S1502)。
s ← eeTs
ついで、差分変換部105は、ノイズの分散共分散行列Σnoizeの特定方向eの大きさと、類似領域の分散共分散行列Σtotalの特定方向eの大きさをそれぞれ求め、両者の比率を1から引くことで、重みwを算出する(S1503)。すなわち、非加法性のノイズを低減できる。
w ← 1−(eTΣnoisee)/(eTΣtotale)
差分変換部105は、射影した差分ベクトルsの大きさを変換する(S1503)。
s ← ws
以上の処理で、色合いを損なうことなく、ノイズ量及びテクスチャ量に応じてノイズ低減の強度を適応的に制御することができる。
【0049】
<S105の変形例>
差分ベクトルsをRGB比率一定の方向に射影するのでなく、射影方向に輝度方向を用いることも有用である。この場合はわずかに色が変わる可能性があるが、RGB空間において輝度方向は一定なので、特定方向eが類似領域毎に(すなわち注目画素に)依存せず一定となり、処理の高速化が図れる。輝度とRGBには、「輝度 = 0.299R + 0.587G + 0.114B」の関係があるので、例えばRGB空間の(0.3 0.6 0.1)と原点を通る方向が特定方向eである。
【0050】
また、射影したベクトルに乗ずる係数も、ノイズ量が多くなれば小さくなり、テクスチャ量が多くなれば大きくなるような特性を持たせることができれば、よりノイズを効率的に低減できる。たとえば、任意の定数αとβを用いて、「α×(ノイズ量)/(ノイズ量+テクスチャ量)+β」とする。α、βを1以上とすれば、重みwをより小さくでき、ノイズを大きく低減することになる。α、βを1より小さくすれば、重みwをより大きくでき、テクスチャやエッジを保存しやすくできる。αとβをパラメータとして、設計者が実現したい画質に合わせて設定することができる。
【0051】
<S106>図3に戻り、ステップS106では、加算部106が変換した差分ベクトルsを類似領域の平均値ベクトルμに加算する。差分ベクトルsを射影して平均値ベクトルμに加算することで、フィルタ処理で削られる高周波を付加となり、振幅の小さい新合成分を保存することが可能となる。なお、実施例2にて説明するように、はじめに図9と同一内容の線形変換の行列を算出し、それと平均値ベクトルμをかけ合わせてもよい。加算部106は、加算して得られたベクトルの要素RGBにより、注目画素の画素値を置き換える。
【0052】
<S107>ステップS107では、画素選択済み判定部107が処理対象画像31のすべての画素を選択したか否かを判定する。まだ選択されていない画素がある場合は、ステップS102に進み、全ての画素を選択した場合は、ノイズ低減装置100がノイズ低減画像32を出力する。ノイズ低減画像32は、例えばフレームメモリ19にいったん記憶され、その後、外部メモリ23又は記憶媒体24に記憶される。
【0053】
〔RGB以外のカラー画像のノイズ低減〕
RGB以外のカラー画像にも、本実施例のノイズ低減方法を好適に適用することができる。RGB以外のカラー画像の場合は、ノイズ低減装置100の前段に、色空間をRGBに変換する処理部を、ノイズ低減装置100の後段に、RGBから元の色空間に戻す処理部を追加すればよい。これらの変換公式は既知である。
【0054】
または、ノイズ低減装置100に、色変換を行わずに処理対象画像31を入力し、色空間の変換が必要となる処理ステップ(類似度算出S103、特定方向への射影S105)でのみ部分的に色変換を行ってもよい。
【0055】
〔ベイヤ配列上での実行〕
また、カラーで撮影する場合、光電変換素子15毎に赤、青、緑の輝度を表すベイヤ(Bayer)配列型のフィルタを有することが多い。このようなベイヤ配列型のフィルタを介して撮影された処理対象画像31に対しても、本発明を適用することができる。
【0056】
例えば、ノイズ低減装置100の前処理部16に、ベイヤ配列のデータをRGBデータに補間する(同時化する)回路を設け、ノイズ低減装置100の後段にRGBデータをベイヤ配列上に間引く処理部を設けることができる。また、例えば、ベイヤ配列型のデータを直接ノイズ低減装置100に入力し、ひとつの画素に対してRGBが必要となる箇所(類似度算出S103、特定方向への射影S105)でのみ一時的にRGB補間(同時化)を行うこともできる。
【0057】
以上説明したように、本実施例のノイズ低減装置100は、ノイズ低減処理による画素値変化を色空間の特定方向eに制約しているため、色ずれを起こさず色ノイズを効果的に低減することができる。また、ノイズ量及びテクスチャ量に応じてノイズ低減の強度を可変にするので、信号値に大きさが依存する非加法性ノイズを、元の画像のテクスチャの量を勘案して適切に低減することができる。
【実施例2】
【0058】
実施例1では、差分ベクトルsを特定方向eに射影して類似領域の平均値ベクトルμに加算するという処理を行った。ここで、特定方向eに、類似領域の平均のRGB比率一定の方向を用いる場合、特定方向eは類似領域の画素の平均値ベクトルμの方向と同じである。
【0059】
図10(a)は、特定方向eとRGB比率一定方向の関係を示す図である。すなわち、平均値ベクトルμの方向と特定方向eが一致する。したがって、差分ベクトルsをRGB一定比率方向に射影して類似領域の画素の平均値ベクトルμに加算するという処理は、平均値ベクトルμをスカラー倍することに等しい。
【0060】
本実施例では、このことを利用し、図3のステップS105及びS106の処理を、適切な乗算数ηを求める処理と、類似領域の画素の平均値ベクトルμにスカラーである乗算数ηを乗算する処理で代替したノイズ低減装置100について説明する。
【0061】
図11はノイズ低減装置100の機能ブロック図の一例を、図12はノイズ低減処理の手順を示すフローチャート図の一例をそれぞれ示す。なお、図11において図2と同一部には同一の符号を付しその説明は省略し、図12において図3と同一部には同一の符号を付しその説明は省略する。
【0062】
<S105−2>ステップS150−2において、乗算数算出部1051は、乗算数ηを算出する。乗算数ηの算出には次式を用いることができる。
【0063】
【数3】
ここで平均値ベクトルμはRGB一定比率方向を、gは注目画素ベクトルを、 Σnoiseはノイズの分散共分散行列Σnoiseを、Σtotalは類似領域の分散共分散行列Σtotalを、それぞれ表す。これは、上述の例でe=μ/|μ|とおいて式を変形することで導くことができるため、同一の処理結果となる。本実施例の処理が処理コストが低いという利点がある。
【0064】
<S106−2>ステップS160−2において、乗算部1061は平均値ベクトルμをスカラーη倍する。
図10(b)は乗算後の類似領域の平均(RGBのベクトル)を示す。乗算前の類似領域のRGB一定比率方向を乗算数η倍することでμ←η×μが得られ、乗算後のμが実施例1の「μ+s」に等しいことがわかる。
【0065】
本実施例のノイズ低減装置100は、より小さい処理コストで実施例1と同等の効果を得ることができる。
【0066】
〔実施例1と実施例2の等価性〕
実施例2において、特定方向eにRGB比率一定方向を用いる場合は、類似領域の平均値ベクトルμを乗算数η倍することで、差分ベクトルsを特定方向eに射影して類似領域の平均値ベクトルμに加算するという処理と等しいと説明した。両者が等価であることを説明する。
【0067】
実施例1の処理は次式で表すことができる。
【0068】
【数4】
式(6)は、差分ベクトル「g−μ」を線形変換してμに加算する処理を意味する。なお、vについては式(5)に示したとおりである。
【0069】
式(6)の右辺は次のように変形することができる。
【0070】
【数5】
式(7)は、μに係数を乗じる処理になっている。また、式(7)の中カッコ{}の中の変数は、式(4)の乗算数ηと同じである。したがって、特定方向eをRGB比率一定方向とする場合、差分ベクトルsを特定方向eに射影して類似領域の平均値ベクトルμに加算するという処理は、類似領域の平均μをη倍する処理と等しいことがわかる。したがって、実施例1と実施例2は等価であることがわかる。
【実施例3】
【0071】
本実施例では、類似領域の画素の最大分散方向に差分ベクトルsを射影しノイズを低減するノイズ低減装置100について説明する。
【0072】
図13は本実施例におけるノイズ低減装置100の機能ブロック図の一例である。図13において図2と同一部には同一の符号を付しその説明は省略する。図13の機能ブロック図では、最大分散方向算出部201を有する点で図2と異なる。また、差分変換部202は実施例1の差分変換部105と異なる動作をする。
【0073】
図14はノイズ低減処理の手順を示すフローチャート図の一例を示す。図14において図3と同一ステップには同一の符号を付しその説明は省略する。
【0074】
<S201>
ステップS104までの処理により、類似領域の平均値ベクトルμ、Σtotal 、及び、Σnoiseが算出されている。ステップS201では、最大分散方向算出部201が、類似領域の画素の本来の画素値の色空間における最大分散方向を算出する。具体的には次の手順で行う。
【0075】
図15は、ステップS201の詳細フローの一例を示す。
【0076】
まず、最大分散方向算出部201は、類似領域の画素の分散共分散行列Σtotalから、ノイズの分散共分散行列Σnoiseを減算する(S2011)。
Σtotal −Σnoise
最大分散方向算出部201は、得られた行列の固有値と固有ベクトルを求める(S2012)。分散共分散行列の固有ベクトルは、与えられたデータの分散がどちらの方向を向いているかを示す。また、固有値は、その方向に対してどの程度の大きさで分散が広がっているか、つまりどの程度ばらついているかを示す。したがって、分散共分散行列の固有値が最大になる固有ベクトルを求めることで、最大分散方向が得られる(S2013)。
【0077】
<S202>
ステップS202では、差分変換部202は、最大分散方向を特定方向eとして実施例1と同様の処理を行う。すなわち、差分変換部202は、注目画素ベクトルgから類似領域の平均値ベクトルμを引き、得られた差分ベクトルsを、最大分散方向に射影する。そして、類似領域の分散共分散行列Σtotalを用いて、射影した差分ベクトルsにノイズを低減する処理を施す。
【0078】
最大分散方向に差分ベクトルsを射影することで、RGB比率一定方向や輝度方向と異なり、元の画像の性質をより強く反映することができる。ノイズ低減による微細信号の劣化を低減することができる。また、エッジを保存しやすくできる。
【実施例4】
【0079】
本実施例では、ノイズ量、テクスチャ量の算出に画素値の分散共分散行列からではなく、空間周波数を用いるノイズ低減装置100について説明する。
【0080】
図16は、本実施例における機能ブロック図の一例を示す。図16において図2と同一部には同一の符号を付しその説明は省略する。図16では、平均・分散出部103の代わりに平均・周波数算出部301を、ノイズ分散推定部104の代わりにノイズ量算出部302を、を有する。また、差分変換部303は、差分変換部105と同等の作用を提供する。
【0081】
図17は、ノイズ低減処理の手順を示すフローチャート図の一例を示す。図17において図3と同一ステップには同一の符号を付しその説明は省略する。
【0082】
<S301>
ステップS301では、平均・周波数算出部301が類似領域の画素の平均値ベクトルμと特定周波数の強度ベクトルを算出する。平均の算出方法は実施例1と同じなので省略する。特定周波数の強度ベクトルを算出するため、平均・周波数算出部301は、類似領域の画素の特定周波数の強度を各色(RGB)成分のフィルタ応答により求める。
【0083】
図18は、フィルタサイズが5×5の場合のフィルタの一例を示す。なお、フィルタ演算には平均値演算に用いた画素のみを使用する。用いるフィルタを調整することで、処理対象画像31から所望の周波数成分を抽出できる。例えば、エッジを抽出するフィルタや平滑化するフィルタが知られている。図18は、高周波成分を抽出するフィルタの例である。
【0084】
処理対象画像31に含まれるノイズは空間周波数に依存する性質を持つ場合があることが知られている。実施例1のように画素値の平均からノイズを推定すると、処理対象画像31の全ての周波数成分からノイズを推定することになり、周波数成分によっては適切なノイズ低減でないことが生じることがある。
【0085】
これに対し、この空間周波数による強度調節の方法は、低減したいノイズの周波数に合わせてノイズ低減強度を調節する側面と、残したいテクスチャの周波数に合わせてノイズ低減強度を調節する側面、の二つ側面を有する。すなわち、ある周波数帯におけるノイズとテクスチャの量をノイズ低減強度に反映させることができる。分散によるノイズの推定が全周波数帯を対象とするのに対して、特定の周波数帯を対象とすることが可能となる。したがって、特定の周波数帯のノイズとテクスチャの量をノイズ低減強度に反映させればいいことが分かっている場合に効果的なノイズ低減処理を実現できる。
【0086】
<S302>
ステップS302では、ノイズ強度算出部302がノイズ強度を算出する。ノイズ強度を算出するため、類似領域の画素の平均からノイズの特定周波数の強度を返す関数を予め同定しておく。この関数は、図7と同様なので省略する。これを各色成分に適用してノイズの特定周波数の強度ベクトルを求める。
【0087】
<S303>
ステップS303では、差分変換部303が差分ベクトルsを変換する。以降の手順は実施例1と同様である。差分変換部303は特定方向eに差分ベクトルsを射影してから、その大きさを変える。特定方向eとしては、RGB比率一定方向、輝度方向、最大分散方向等とすることができる。
【0088】
差分変換部303は、ステップS302で算出された、ノイズの特定周波数の強度ベクトルを特定方向eに射影したベクトルの大きさをノイズ量とする。また、特定方向eに射影された差分ベクトルsからノイズの特定周波数の強度ベクトルを減算して得られたベクトルの大きさをテクスチャ量とする。このノイズ量とテクスチャ量を用いて、例えば、ステップS1503のようにして差分ベクトルsの射影ベクトルの大きさを変更する。なお、特定方向eをRGB比率一定方向とした場合は、実施例2のように、類似領域の平均値ベクトルμをη倍してもよい。
【0089】
本実施例のノイズ低減装置100は、フィルタにより特定周波数の成分を抽出しノイズの強度を算出するので、所望の周波数成分のノイズを低減しやすくできる。
【実施例5】
【0090】
これまでの実施例では、差分ベクトルsの射影方向を特定方向eに変換したが、本実施例では、差分ベクトルsを、特定方向eと特定方向eに垂直な方向のそれぞれのノイズ量に応じて変換するノイズ低減装置100について説明する。いわば、差分ベクトルsをRGB空間における二次元部分空間にて変換することになる。
【0091】
図19は、本実施例における機能ブロック図の一例を示す。図19において図2と同一部には同一の符号を付しその説明は省略する。図19では、差分変換部401を有する。
【0092】
図20は、ノイズ低減処理の手順を示すフローチャート図の一例を示す。図20において図3と同一ステップには同一の符号を付しその説明は省略する。
【0093】
<S401>
処理ステップS401では、差分変換部401が、差分ベクトルsの、特定方向eの分散と、特定方向eに垂直な方向e⊥の分散、それぞれついて大きさを変換する。
図21は、差分ベクトルsの二次元方向への射影を模式的に示す図の一例である。図21(a)に示すように、差分ベクトルsと特定方向eを含む平面が一意に定まる。また、ノイズ量はRGB成分毎に求められている。特定方向eの分散の求め方は実施例1と同じである。特定方向eに垂直な方向e⊥のノイズ量は、ノイズの分散共分散行列Σnoizeのうち特定方向eに垂直な方向e⊥の分散である。特定方向eに垂直な方向e⊥のテクスチャ量は、類似領域の分散共分散行列Σtotalからノイズの分散共分散行列Σnoizeを引いて得られる分散共分散行列の特定方向eに垂直な方向e⊥の分散である。
【0094】
図21(b)は、それぞれの方向において(テクスチャ量)/(ノイズ量+テクスチャ量)で大きさを変換した差分ベクトルsの一例を示す。特定方向eと、垂直な方向e⊥の2つの方向に対し、差分ベクトルsを各方向のノイズ量に応じて変換することで、特定方向eに近づける方向で変換できる。したがって、特定方向eと垂直な方向のテクスチャ量がノイズ量に比して大きければ、射影される差分ベクトルsに注目画素の情報を多く残すことができる。逆に、特定方向eと垂直な方向のテクスチャ量がノイズ量に比して小さければ、射影される差分ベクトルsを類似領域の平均値ベクトルμに近くすることでノイズを低減できる。なお、特定方向eはRGB比率一定方向、輝度方向、又は、最大分散方向のいずれでもよい。
【0095】
また、例えば特定方向eを輝度方向に取った場合、特定方向eに垂直な方向e⊥の差分ベクトルsの大きさを特定方向eの大きさよりも高い比率で縮小することができる。輝度方向と垂直な方向は人間の眼の感度が低いので、特定方向eに垂直方向の差分ベクトルsをその方向のノイズ量より大きく圧縮しても、人間の視覚が感じる変化は小さい。したがって、人間の視覚に重要な輝度成分を保存したまま、効率的にカラーノイズを低減することができる。
【0096】
本実施例のノイズ低減装置100によれば、差分ベクトルsをRGB空間における二次元部分空間にて変換することで、ノイズ量に応じて注目画素の情報を適切に圧縮しながら、カラー画像のノイズを低減できる。
【実施例6】
【0097】
実施例1では差分ベクトルsを特定方向eの一次元方向に射影し、実施例5では差分ベクトルsを二次元部分空間のベクトルに変換することでノイズを低減した。本実施例では、差分ベクトルsを三次元部分空間のベクトルに変換することでノイズを低減するノイズ低減装置100について説明する。すなわち、本実施例のノイズ低減装置100は、三次元の色空間(RGB色空間)におけるノイズの低減を一般化したものである。
【0098】
なお、本実施例の機能ブロック図は、図2又は図19と同様であり、特に差分変換部105又は差分変換部401を三次元空間の処理に拡張したものである。また、ノイズ低減の手順は図3又は図20と等しい。
【0099】
実施例1,5を一般化すると、差分ベクトルsを線形変換し、類似領域の平均値ベクトルに加算する処理となる。具体的には、類似領域の平均値ベクトルμ、注目画素値ベクトルg、線形変換の行列Wを用いて次の式で出力画素値fを算出する処理に相当する。
【0100】
【数6】
行列Wはたとえば次のように設定することができる。
【0101】
【数7】
Σnoiseはノイズの分散共分散行列、Σtextureはテクスチャの分散共分散行列である。テクスチャの分散共分散行列Σtextureは、類似領域の分散共分散行列Σtotal からノイズの分散共分散行列Σnoiseを引いたものである。したがって、式(9)から、行列Wは、三次元空間のそれぞれの方向で(テクスチャ量)/(ノイズ量+テクスチャ量)の演算を可能とする行列であることがわかる。
【0102】
式(8)は、差分ベクトル(g−μ)を三次元空間で線形変換し、平均値ベクトルμに加算する処理であるので、色空間における所望の3方向、それぞれについてノイズ量とテクスチャ量のバランスを勘案して、差分ベクトルsを三次元空間のベクトルに変換することができる。
【0103】
図22は、差分ベクトルsの三次元方向への射影を模式的に示す図の一例である。方向e1, e2, e3はWの固有ベクトルの方向である。注目画素gが含むノイズ量を求めておけば、図示するように、方向e1, e2, e3それぞれについて、(テクスチャ量)/(ノイズ量+テクスチャ量)に対応した変換が可能になる。変換後の各方向の分散から出力画素値のベクトルfが求められる。
【0104】
ノイズ量に対してテクスチャ量が多い方向では注目画素値に近い値となり、出力画素値に注目画素の情報を多く残すことができる。その逆に、ノイズ量に対してテクスチャ量が少ない方向では類似領域の平均に近い値となるので適切にノイズを低減できる。
【0105】
また、特定方向e0(単位ベクトル)の分散を相対的に多く残すことで、人間の知覚上有用な成分を選択的に多く残すことも可能である。そのためには線形変換の行列Wを次のように設定する。
【0106】
【数8】
ここでIは3×3の単位行列である。αは設計パラメータであり、特定方向e0の分散の強調の程度を指定する。αが大きいほど、特定方向e0の分散を強調しやすくできる。特定方向e0としてはこれまで説明したようにRGB比率一定方向、輝度方向、最大分散方向等を取ることができる。
【0107】
また、より設計の自由度を高く持たせて、3×3のパラメータ行列A, Bを用いて、次のように設定してもよい。
【0108】
【数9】
式(11)は、実施例1の<S105の変形例>で説明した、「α×(ノイズ量)/(ノイズ量+テクスチャ量)+β」と同等である。すなわち、AとBを調整することで、射影したベクトルに乗ずる係数も、ノイズ量が多くなれば小さくなり、テクスチャ量が多くなれば大きくなるような特性を持たせることができる。
【0109】
本実施例のノイズ低減装置100は、RGB空間のノイズの低減を、三次元部分空間に一般化して、各方向のノイズ量の大きさに応じて、注目画素の情報を多く残したり、類似領域の平均に近い値としたりすることができる。
【0110】
以上説明したように、本実施形態のノイズ低減装置100は、注目画素から類似度が閾値内の画素をフィルタ処理することで、エッジなど人間が画像を知覚する上で重要な要素を保存しつつノイズ低減を行うことができ、差分ベクトルを例えば平均値ベクトルμに加算することで振幅の小さい信号を保存でき、カラー画像における非加法性のノイズを推定してノイズ量とテクスチャ量に応じてフィルタ処理の出力を変換することで、信号値に大きさが依存する非加法性ノイズを、元の画像のテクスチャの量を勘案して適切に低減することができる。また、ノイズ低減処理による画素値変化を色空間の特定方向に制約しているため、色ずれを起こさず色ノイズを効果的に低減することができる。
【符号の説明】
【0111】
11 撮像部
20 ノイズ低減プログラム
24 記憶媒体
31 処理対象画像
32 ノイズ低減画像
100 ノイズ低減装置
101 画像取得部
102 注目画素選択部
103 平均・分散推定部
104 ノイズ分散推定部
105 差分変換部
106 加算部
107 画素選択済み判定部
【先行技術文献】
【特許文献】
【0112】
【特許文献1】特開2008−205737号公報
【特許文献2】特開2007−288439号公報
【特許文献3】特開2002−259965号公報
【非特許文献】
【0113】
【非特許文献1】原島 博, 小田島 薫, 鹿喰 善明, 宮川 洋, "ε-分離非線形ディジタルフィルタとその応用," 電子情報通信学会論文誌 A , Vol.J65-A, No.4, pp.297-304, 1982.
【非特許文献2】CC. Tomasi and R. Manduchi, "Bilateral Filtering for Gray and Color Images," Proc. Sixth Int'l Conf. Computer Vision, pp. 839-846, 1998.
【非特許文献3】Ce Liu, et al., "Automatic Estimation and Removal of Noise from a Single Image," IEEE Trans. Pattern Analysis and Machine Intelligence, vol. 30, no. 2, pp. 299-314, 2008.
【技術分野】
【0001】
本発明は、カラー画像のノイズを低減するノイズ低減理装置、ノイズ低減方法、プログラム及び記憶媒体に関する。
【背景技術】
【0002】
スキャナやデジタルカメラ等の撮像機器により撮影した画像には、撮像素子及び回路の特性上、暗電流ノイズ、熱雑音、ショットノイズなどのノイズが含まれる。高画質の画像を得るためには、これらのノイズを低減する処理を施さなければならない。しかし、単純にローパスフィルタを用いてノイズを低減すると、エッジなど人間が画像を知覚する上で重要な要素も同時に失われて画質が劣化する。そこで、画像の領域の特性に応じて、適応的にノイズ低減を行う技術が必要となる。
【0003】
そのようなノイズ低減技術の一つとして、εフィルタが提案されている(例えば、非特許文献1参照。)。εフィルタでは、注目画素との信号差が一定閾値以下となる周辺画素をフィルタ処理に用いることで、エッジなどの信号差の大きい成分を保存しつつノイズを低減することを可能としている。
【0004】
また、バイラテラルフィルタを用いたノイズ低減技術が提案されている(例えば、非特許文献2参照。)。バイラテラルフィルタは、注目画素との信号差及び空間距離によってフィルタの重み係数を作成しフィルタ処理を行う方法であり、εフィルタ同様、エッジの保存とノイズ低減を両立させることができる。なお、バイラテラルフィルタに準ずる原理によるノイズ低減は他にも提案されている(例えば、特許文献1,2参照。)。
しかしながら、εフィルタやバイラテラルフィルタには、背景部分のような振幅の小さい信号成分が失われやすいという問題がある。そこで、この問題を解決する技術が提案されている(例えば、特許文献3参照。)。特許文献3には、εフィルタ出力値と注目画素値を注目画素近傍の空間周波数に応じた割合で加重加算する技術が開示されている。
【0005】
しかしながら、特許文献3記載の方法でも、ショットノイズなど非加法性のノイズに対応できず、また、カラー画像に対するノイズ低減性能が悪い、という問題がある。カラー画像における非加法性ノイズの低減技術が提案されている(例えば、非特許文献3参照。)が、この方法は画像のセグメンテーション、統計量の算出に多大な処理コストがかかるという問題がある。
【発明の概要】
【発明が解決しようとする課題】
【0006】
本発明は上記従来技術の問題を鑑みたものであり、エッジなど人間が画像を知覚する上で重要な要素を保存しつつノイズ低減を行い、さらに振幅の小さい信号を保存し、カラー画像における非加法性のノイズを高精度に少ない処理コストで低減する、ノイズ低減理装置、ノイズ低減方法、プログラム及び記憶媒体を提供することを目的とする。
【課題を解決するための手段】
【0007】
カラー画像の注目画素の画素値に、注目画素を含む近傍領域の画素の画素値を用いてノイズ低減処理を施すノイズ低減装置であって、所定の色空間において注目画素との類似度が閾値以上の類似領域に含まれる画素を前記近傍領域から抽出し、前記類似領域の画素の代表値ベクトルを算出する代表値ベクトル算出手段と、注目画素の注目画素ベクトルと前記代表値ベクトルとの差分ベクトルを、前記色空間の特定方向に射影する差分射影手段と、前記代表値ベクトルに、前記特定方向に射影された前記差分ベクトルを加えたベクトルの要素を、注目画素の画素値に置き換える画素値修正手段と、を有することを特徴とする。
【発明の効果】
【0008】
エッジなど人間が画像を知覚する上で重要な要素を保存しつつノイズ低減を行い、カラー画像における非加法性のノイズを適切に低減する、ノイズ低減理装置、ノイズ低減方法、プログラム及び記憶媒体を提供することができる。
【図面の簡単な説明】
【0009】
【図1】デジタルカメラの概略構成図の一例である。
【図2】ノイズ低減装置の機能ブロック図の一例である。
【図3】ノイズ低減処理の手順を示すフローチャート図の一例である。
【図4】平均、分散の算出を模式的に説明する図の一例である。
【図5】ステップS103の詳細フローの一例を示すフローチャート図である。
【図6】信号値(画素値)に対するノイズの分散をプロットした図の一例である。
【図7】ノイズの分散を推定するための関数の一例を示す図である。
【図8】差分ベクトルsの変換を模式的に説明する図の一例である。
【図9】ステップS105の詳細フローの一例を示すフローチャート図である。
【図10】特定方向eと類似領域の平均値のRGB比率一定方向の関係を示す図である。
【図11】ノイズ低減装置の機能ブロック図の一例である(実施例2)。
【図12】ノイズ低減処理の手順を示すフローチャート図の一例である(実施例2)。
【図13】ノイズ低減装置の機能ブロック図の一例である(実施例3)。
【図14】ノイズ低減処理の手順を示すフローチャート図の一例である(実施例3)。
【図15】ステップS201の詳細フローの一例を示すフローチャート図である。
【図16】ノイズ低減装置の機能ブロック図の一例である(実施例4)。
【図17】ノイズ低減処理の手順を示すフローチャート図の一例である(実施例4)。
【図18】フィルタサイズが5×5の場合のフィルタの一例を示す図である。
【図19】ノイズ低減装置の機能ブロック図の一例である(実施例5)。
【図20】ノイズ低減処理の手順を示すフローチャート図の一例を示す(実施例5)。
【図21】差分ベクトルの二次元方向への射影を模式的に示す図の一例である。
【図22】差分ベクトルの三次元方向への射影を模式的に示す図の一例である。
【発明を実施するための形態】
【0010】
以下、本発明を実施するための形態について図面を参照しながら実施例を挙げて説明する。
【実施例1】
【0011】
本実施例では、一つの画素がRGBの3つの成分を持つカラー画像から、信号値に大きさが依存する非加法性のノイズを低減するノイズ低減装置100について説明する。以下、信号値を画素値という場合がある。本実施例のノイズ低減装置100は、デジタルカメラ、スキャナ、カメラ付き携帯などカラー画像を撮影する機器であれば好適に適用できる。また、コンピュータが記憶媒体から読み出したり、ネットワークを介して取得したカラー画像に本実施例のノイズ低減処理を施すこともできる。
【0012】
図1は、デジタルカメラの概略構成図の一例を示す。デジタルカメラの構成そのものは公知のものを利用できる。デジタルカメラは、例えば、像を撮影する撮像部11と、撮像部11により得られた画像信号に処理を施す信号処理部18と、デジカメ全体を制御する主制御部17、画像データを記憶するフレームメモリ19、及び、その他の部品と接続されるI/F(インタフェース)21、を有する。I/F21 には画像を表示する表示部22と、画像データを保存したり画像データを読み出すための外部メモリ23、及び、記憶媒体24を装着するメモリカード装着部24とが接続されている。ノイズ低減処理は、例えば信号処理部18と主制御部17に実装されるが、実装例は図に限定されるものではない。
【0013】
撮像部11は、レンズ12、絞り13、電子シャッター14、光電変換素子15及び前処理部16を有する。光電変換素子15は、例えば、CMOS(Complementary Metal Oxide Semiconductor)やCCD(Charge Coupled Device)が使用される。撮像部11は不図示のカラーフィルタ(原色系又は補色系のいずれでもよい)を有し、例えばセル状に配列されたR、G、B各色毎に1つの光電変換素子15が配置されている。
【0014】
前処理部16は、プリアンプやAGC(自動利得制御)等のアナログ信号処理やA/D変換器を備えており、光電変換素子15により出力されたアナログ映像信号に対して増幅・クランプ等の前処理を施した後、アナログ映像信号をデジタル映像信号に変換する。
【0015】
信号処理部18は、DSP等を有し、撮像部11から得られたデジタル信号に対して色分解、ホワイトバランス調整、γ補正等の種々の画像処理を施す。信号処理部18は処理した画像データをフレームメモリ19に格納したり、格納した画像データをフレームメモリ19から読み出し、ノイズ低減などの画像処理を施すことができる。フレームメモリ19には、VRAM、SRAM、DRAM等の半導体メモリが使用される。
【0016】
フレームメモリ19から読み出された画像データは、信号処理部18において画像圧縮等の信号処理が施された後、I/F21を介して外部メモリ23又はメモリカード装着部24に装着された記憶媒体25に記憶される。外部メモリ23は、フラッシュメモリ等の不揮発メモリであり、記憶媒体25はUSBメモリ、SDメモリカード、光磁気ディスクなどの脱着可能な可搬型の不揮発メモリである。IF21はこれらを切り替えて利用できる。また、不図示の通信ユニットにより画像データをサーバに送信したり、サーバから画像データを受信して外部メモリ23に記憶してもよい。この場合、通信ユニットは、携帯電話網や無線/有線LAN網に接続して、画像データを送受信する。
【0017】
フレームメモリ19、外部メモリ23又は記憶媒体25から読み出された画像データは、表示部22に表示される。表示部22は、液晶ディスプレイや有機ELディスプレイ等、デジタルカメラの筐体に設置されている。表示部22はタッチパネルを一体に有する。このタッチパネルと、筐体に設けられたキーボードにより、ユーザの操作を入力することができる。
【0018】
主制御部17、又は、主制御部17と信号処理部18は、マイコンやLSIなどにより構成される。したがって、例えば、主制御部、信号処理部及びフレームメモリ19が、ノイズ低減装置100を構成する主要部の一例である。マイコンは、主にCPU、RAM、EEPROM、ASIC、等がバスで接続されたコンピュータを実体とし、EEPROMに記憶されたノイズ低減プログラム20をCPUが実行することで、以下で説明するノイズ低減処理が画像データに施される。また、ASIC等のハードウェアによりノイズ低減処理を実装してもよい。ノイズ低減プログラム20は、予め、EEPROMに記憶された状態でデジタルカメラと共に出荷される。また、記憶媒体25にノイズ低減プログラム20を記憶させて配布し、I/F21を介してノイズ低減プログラム20をEEPROMに読み込むことができる。また、ネットワークを介してノイズ低減プログラム20をEEPROMにダウンロードしてもよい。
【0019】
〔ノイズ低減処理〕
図2はノイズ低減装置100の機能ブロック図の一例を、図3はノイズ低減処理の手順を示すフローチャート図の一例をそれぞれ示す。
【0020】
ノイズ低減装置100は画像データ(以下、処理対象画像31という)を入力とし、ノイズ低減画像32を出力する。ノイズ低減装置100は、画像取得部101、注目画素選択部102、平均・分散算出部103、ノイズ分散推定部104、差分変換部105、加算部106、及び、画素選択済み判定部107を有する。各機能ブロックは、CPUがノイズ低減プログラム20を実行することで実現される。以下、各機能ブロックについて、図3のフローチャート図に沿って説明する。
【0021】
<S101>ステップS101では、画像取得部101がフレームメモリ19から処理対象画像31を取得する。処理対象画像31は、デジタルカメラが直接撮影した以外の画像データでもよい。
【0022】
<S102>ステップS102では、注目画素選択部102が、処理対象画像31の複数の画素の中から注目画素を選択する。たとえばラスタ走査により順次、注目画素を選択すればよく、全ての画素を選択する方法であれば走査方法はどのようなものでもよい。
【0023】
<S103>ステップS103では、平均・分散算出部103が、注目画素の近傍領域における画素の画素値の平均、分散を算出する。ここではεフィルタやバイラテラルフィルタなどの方法を取ることができる。一例として次の方法を示す。
【0024】
図4は、平均、分散の算出を模式的に説明する図の一例である。近傍画素は実際にはカラーである。まず、平均・分散算出部103は、注目画素を中心とした特定ウィンドウ内の近傍画素を選択する。特定ウィンドウの大きさは予め定められていてもよいし、ユーザが設定したノイズ低減処理の強度に応じて可変となっていてもよい。図では(9×9)画素の矩形が特定ウィンドウである。
【0025】
そして、平均・分散算出部103が、注目画素と各近傍画素の類似度を算出する。類似度として、例えばRGB空間でのユークリッド距離の逆数を用いることができる。類似度が特定閾値以上の近傍画素を、平均、分散の演算に用いることとする。こうすることで、注目画素と色合いの近似した近傍画素を抽出し、エッジを保存しながらノイズを低減できる。なお、類似度は画素単位で算出し、平均値はR、G、B毎に算出する。
【0026】
図4では、注目画素を中心とする球内にある画素が、注目画素との類似度が特定閾値以上の近傍画素となる。以下、特定半径の球内を類似領域という。
【0027】
平均・分散算出部103は、次の公式を用いて分散を算出する。
【0028】
【数1】
ここでVar(X)は確率変数Xの分散、Cov(X, Y)は確率変数XとYの共分散、E(X)は確率変数Xの期待値を示す。すなわち、画素値を二乗したデータ、異なる色成分の画素値を乗算したデータ、元の画素値、を組み合わせることで、類似領域の画素の分散・共分散(フィルタ処理)を算出することができる。なお、処理コスト削減の面から共分散は算出せず0と置いてもよい。
【0029】
図5は、ステップS103の詳細フローの一例を示す。xは、近傍画素の画素値であり、RGBの各成分を要素として持つ3次元ベクトルである。平均・分散算出部103は、変数μ、Σ、nに初期値であるゼロを設定する(S1301)。
【0030】
ついで、平均・分散算出部103は、特定ウィンドウの近傍画素の画素値を取得し(S1302)、注目画素との類似度を算出する(S1303)。
【0031】
平均・分散算出部103は、類似度が特定閾値以上か否かを判定する(S1304)。類似度が特定閾値以上の場合(S1304のYes)、その近傍画素が平均、分散の算出対象なので、平均・分散算出部103は、変数μにxを、変数ΣにxxTを、変数nにn+1を、それぞれ加える(S1305)。なお、「T」は転置を表す。
【0032】
ついで、平均・分散算出部103は、特定ウィンドウ内の全ての近傍画素を選択したか否かを判定する(S1306)。特定ウィンドウ内の全ての近傍画素を選択していない場合(S1305のNo)、ステップS1032に戻り上記の処理を繰り返す。
【0033】
特定ウィンドウ内の全ての近傍画素を選択した場合(S1305のYes)、平均・分散算出部103は、「μ/n」、「Σ/n−μμT」をそれぞれ算出する(S1307)。二乗和から平均を引く一般的な分散算出と同様の結果が、少ない計算量で得られる。
【0034】
図5の処理の結果、μが類似領域の平均値ベクトルとして、Σが分散共分散行列として得られる。Σの対角成分は類似領域の分散である。なお、類似領域の平均でなく、中央値や最頻値など種々の代表値を採用することができる。平均等の代表値に、輝度やホワイトバランスを修正する定数を加えたり乗じてもよい。
【0035】
<S103の変形例>
ステップS103の平均、分散の算出方法には種々のバリエーションを取りうる。例えば、実施例1では類似領域の画素に対し単純に平均・分散を求めたが、類似領域の画素の画素値に重み付けをして平均・分散を求めてもよい。バイラテラルフィルタのように、注目画素との類似度と相対位置をフィルタの重み係数に反映させることも有効である。例えば、平均・分散算出部103は、類似度が大きいほど大きな重み係数を画素値に与えて平均・分散を算出する。
【0036】
また、上記実施例ではRGB空間でのユークリッド距離を類似度の算出に用いたが、類似度算出にLnノルム(nは自然数)を用いてもよい。また、色空間には、L*a*b*表色系、XYZ表色系等、適宜選択できる。この場合、平均・分散算出部103は、特定ウィンドウ内の画素の画素値(RGB)を選択した表色系の画素値に変換する。変換公式は既知である。
【0037】
平均・分散算出部103は、このRGB空間とは別の色空間でのノルムを用いることもできる。たとえば、L*a*b*などの等色空間のユークリッドノルムを用いると、人間の知覚に近い色の類似度を得ることができる。
【0038】
<S104>図3に戻りステップS104では、ノイズ分散推定部104が、ノイズの分散を推定する。
図6は信号値(画素値)に対するノイズの分散をプロットした図の一例である。ショットノイズでは信号値とノイズの分散に比例関係が見られる。このように、一般にノイズの分散は画素値に依存する。
【0039】
したがって、画素値からノイズの分散を推定する関数を導出しておくことが可能になる。なお、ノイズの分散はRGBの各成分ごとに求めるので、ノイズの分散を推定する関数もR、G、Bの各成分毎に導出しておく。
【0040】
注目画素のノイズを低減するためには、注目画素に含まれる平均的なノイズを知りたいので、ノイズ分散推定部104は、類似領域のノイズの分散を求める。図7は、類似領域のノイズの分散を推定するための関数の一例を示す。ノイズ分散推定部104は、平均値ベクトルμの各成分に基づき図7の関数を利用して、類似領域のノイズの分散σ2を求める。
【0041】
これを各色成分に適用して各成分のノイズの分散を求める。求めたノイズの分散を対角成分に置いた行列をノイズの分散共分散行列Σnoizeとする。多くの場合はノイズの分散共分散行列の非対角成分(共分散)を0と置くことができる。ノイズの共分散が無視できない場合は、分散推定間数と同様に共分散の推定関数を同定し、それを用いて共分散を推定し分散共分散行列を算出する。
【0042】
<S105>図3に戻りステップS105では、差分変換部105が注目画素gの注目画素ベクトルgから類似領域の平均値ベクトルμを引き、得られた差分(以下、差分ベクトルsという)の方向と大きさを変換する。
【0043】
図8は、差分ベクトルsの変換を模式的に説明する図の一例である。図8(a)に示すように、差分変換部105は、まず、特定方向eに差分ベクトルsを射影する。特定方向eとして、類似領域の平均値ベクトルμの成分(R、G、B)と同じRGB比率となる一定方向(以下、RGB一定比率方向という)を取ることができる。この方向を特定方向eとすることで、注目画素の画素値が類似領域色合いとずれることを防ぐことができる。また、射影することで、差分ベクトルsが、類似領域の色合いの方向に変換され、かつ、長さが低減されるので、類似領域の画素と異なる色合いをもたらすノイズやその大きさを低減する効果が得られる。
【0044】
そして、更に好ましくは、差分変換部105は、類似領域のノイズの分散共分散行列Σnoizeを用いて射影した差分ベクトルsに明示的にノイズを低減する処理を施す。図8(b)に示すように、差分変換部105は、射影したベクトルの大きさに(テクスチャ量)/(ノイズ量+テクスチャ量)の値を乗じる。これにより、ノイズ量及びテクスチャ量に応じてノイズ低減の強度を適応的に制御することができる。
【0045】
なお、ノイズ量はノイズの分散共分散行列Σnoizeのうち特定方向eの分散、テクスチャ量は類似領域の分散共分散行列からノイズの分散共分散行列Σnoizeを引いて得られる分散共分散行列の特定方向eの分散である。一般に分散共分散行列がΣであるとき、特定方向e(単位ベクトル)の分散σe2は次式で求めることができる。
【0046】
【数2】
テクスチャ量は類似領域の分散共分散行列からノイズの分散共分散行列Σnoizeを引いて得られるものであるから、実質的には射影した差分ベクトルsの大きさに(テクスチャ量)/(ノイズ量+テクスチャ量)={1−(ノイズの分散のe方向の大きさ)/(類似領域の分散のe方向の大きさ)}を乗算すればよい。
【0047】
図9は、ステップS105の詳細フローの一例を示す。図9の手順は一例に過ぎず、同様の結果を得られる限り各ステップの内容・順序を入れ替えることができる。注目画素ベクトルgは注目画素の画素値(RGBのベクトル)、Σtotal はステップS1305のΣと同じもので類似領域の分散共分散行列である。
【0048】
まず、差分変換部105は、差分ベクトルsを算出する(S1501)。
s ← g−μ
ついで、差分変換部105は、差分ベクトルsを特定方向eに射影する(S1502)。
s ← eeTs
ついで、差分変換部105は、ノイズの分散共分散行列Σnoizeの特定方向eの大きさと、類似領域の分散共分散行列Σtotalの特定方向eの大きさをそれぞれ求め、両者の比率を1から引くことで、重みwを算出する(S1503)。すなわち、非加法性のノイズを低減できる。
w ← 1−(eTΣnoisee)/(eTΣtotale)
差分変換部105は、射影した差分ベクトルsの大きさを変換する(S1503)。
s ← ws
以上の処理で、色合いを損なうことなく、ノイズ量及びテクスチャ量に応じてノイズ低減の強度を適応的に制御することができる。
【0049】
<S105の変形例>
差分ベクトルsをRGB比率一定の方向に射影するのでなく、射影方向に輝度方向を用いることも有用である。この場合はわずかに色が変わる可能性があるが、RGB空間において輝度方向は一定なので、特定方向eが類似領域毎に(すなわち注目画素に)依存せず一定となり、処理の高速化が図れる。輝度とRGBには、「輝度 = 0.299R + 0.587G + 0.114B」の関係があるので、例えばRGB空間の(0.3 0.6 0.1)と原点を通る方向が特定方向eである。
【0050】
また、射影したベクトルに乗ずる係数も、ノイズ量が多くなれば小さくなり、テクスチャ量が多くなれば大きくなるような特性を持たせることができれば、よりノイズを効率的に低減できる。たとえば、任意の定数αとβを用いて、「α×(ノイズ量)/(ノイズ量+テクスチャ量)+β」とする。α、βを1以上とすれば、重みwをより小さくでき、ノイズを大きく低減することになる。α、βを1より小さくすれば、重みwをより大きくでき、テクスチャやエッジを保存しやすくできる。αとβをパラメータとして、設計者が実現したい画質に合わせて設定することができる。
【0051】
<S106>図3に戻り、ステップS106では、加算部106が変換した差分ベクトルsを類似領域の平均値ベクトルμに加算する。差分ベクトルsを射影して平均値ベクトルμに加算することで、フィルタ処理で削られる高周波を付加となり、振幅の小さい新合成分を保存することが可能となる。なお、実施例2にて説明するように、はじめに図9と同一内容の線形変換の行列を算出し、それと平均値ベクトルμをかけ合わせてもよい。加算部106は、加算して得られたベクトルの要素RGBにより、注目画素の画素値を置き換える。
【0052】
<S107>ステップS107では、画素選択済み判定部107が処理対象画像31のすべての画素を選択したか否かを判定する。まだ選択されていない画素がある場合は、ステップS102に進み、全ての画素を選択した場合は、ノイズ低減装置100がノイズ低減画像32を出力する。ノイズ低減画像32は、例えばフレームメモリ19にいったん記憶され、その後、外部メモリ23又は記憶媒体24に記憶される。
【0053】
〔RGB以外のカラー画像のノイズ低減〕
RGB以外のカラー画像にも、本実施例のノイズ低減方法を好適に適用することができる。RGB以外のカラー画像の場合は、ノイズ低減装置100の前段に、色空間をRGBに変換する処理部を、ノイズ低減装置100の後段に、RGBから元の色空間に戻す処理部を追加すればよい。これらの変換公式は既知である。
【0054】
または、ノイズ低減装置100に、色変換を行わずに処理対象画像31を入力し、色空間の変換が必要となる処理ステップ(類似度算出S103、特定方向への射影S105)でのみ部分的に色変換を行ってもよい。
【0055】
〔ベイヤ配列上での実行〕
また、カラーで撮影する場合、光電変換素子15毎に赤、青、緑の輝度を表すベイヤ(Bayer)配列型のフィルタを有することが多い。このようなベイヤ配列型のフィルタを介して撮影された処理対象画像31に対しても、本発明を適用することができる。
【0056】
例えば、ノイズ低減装置100の前処理部16に、ベイヤ配列のデータをRGBデータに補間する(同時化する)回路を設け、ノイズ低減装置100の後段にRGBデータをベイヤ配列上に間引く処理部を設けることができる。また、例えば、ベイヤ配列型のデータを直接ノイズ低減装置100に入力し、ひとつの画素に対してRGBが必要となる箇所(類似度算出S103、特定方向への射影S105)でのみ一時的にRGB補間(同時化)を行うこともできる。
【0057】
以上説明したように、本実施例のノイズ低減装置100は、ノイズ低減処理による画素値変化を色空間の特定方向eに制約しているため、色ずれを起こさず色ノイズを効果的に低減することができる。また、ノイズ量及びテクスチャ量に応じてノイズ低減の強度を可変にするので、信号値に大きさが依存する非加法性ノイズを、元の画像のテクスチャの量を勘案して適切に低減することができる。
【実施例2】
【0058】
実施例1では、差分ベクトルsを特定方向eに射影して類似領域の平均値ベクトルμに加算するという処理を行った。ここで、特定方向eに、類似領域の平均のRGB比率一定の方向を用いる場合、特定方向eは類似領域の画素の平均値ベクトルμの方向と同じである。
【0059】
図10(a)は、特定方向eとRGB比率一定方向の関係を示す図である。すなわち、平均値ベクトルμの方向と特定方向eが一致する。したがって、差分ベクトルsをRGB一定比率方向に射影して類似領域の画素の平均値ベクトルμに加算するという処理は、平均値ベクトルμをスカラー倍することに等しい。
【0060】
本実施例では、このことを利用し、図3のステップS105及びS106の処理を、適切な乗算数ηを求める処理と、類似領域の画素の平均値ベクトルμにスカラーである乗算数ηを乗算する処理で代替したノイズ低減装置100について説明する。
【0061】
図11はノイズ低減装置100の機能ブロック図の一例を、図12はノイズ低減処理の手順を示すフローチャート図の一例をそれぞれ示す。なお、図11において図2と同一部には同一の符号を付しその説明は省略し、図12において図3と同一部には同一の符号を付しその説明は省略する。
【0062】
<S105−2>ステップS150−2において、乗算数算出部1051は、乗算数ηを算出する。乗算数ηの算出には次式を用いることができる。
【0063】
【数3】
ここで平均値ベクトルμはRGB一定比率方向を、gは注目画素ベクトルを、 Σnoiseはノイズの分散共分散行列Σnoiseを、Σtotalは類似領域の分散共分散行列Σtotalを、それぞれ表す。これは、上述の例でe=μ/|μ|とおいて式を変形することで導くことができるため、同一の処理結果となる。本実施例の処理が処理コストが低いという利点がある。
【0064】
<S106−2>ステップS160−2において、乗算部1061は平均値ベクトルμをスカラーη倍する。
図10(b)は乗算後の類似領域の平均(RGBのベクトル)を示す。乗算前の類似領域のRGB一定比率方向を乗算数η倍することでμ←η×μが得られ、乗算後のμが実施例1の「μ+s」に等しいことがわかる。
【0065】
本実施例のノイズ低減装置100は、より小さい処理コストで実施例1と同等の効果を得ることができる。
【0066】
〔実施例1と実施例2の等価性〕
実施例2において、特定方向eにRGB比率一定方向を用いる場合は、類似領域の平均値ベクトルμを乗算数η倍することで、差分ベクトルsを特定方向eに射影して類似領域の平均値ベクトルμに加算するという処理と等しいと説明した。両者が等価であることを説明する。
【0067】
実施例1の処理は次式で表すことができる。
【0068】
【数4】
式(6)は、差分ベクトル「g−μ」を線形変換してμに加算する処理を意味する。なお、vについては式(5)に示したとおりである。
【0069】
式(6)の右辺は次のように変形することができる。
【0070】
【数5】
式(7)は、μに係数を乗じる処理になっている。また、式(7)の中カッコ{}の中の変数は、式(4)の乗算数ηと同じである。したがって、特定方向eをRGB比率一定方向とする場合、差分ベクトルsを特定方向eに射影して類似領域の平均値ベクトルμに加算するという処理は、類似領域の平均μをη倍する処理と等しいことがわかる。したがって、実施例1と実施例2は等価であることがわかる。
【実施例3】
【0071】
本実施例では、類似領域の画素の最大分散方向に差分ベクトルsを射影しノイズを低減するノイズ低減装置100について説明する。
【0072】
図13は本実施例におけるノイズ低減装置100の機能ブロック図の一例である。図13において図2と同一部には同一の符号を付しその説明は省略する。図13の機能ブロック図では、最大分散方向算出部201を有する点で図2と異なる。また、差分変換部202は実施例1の差分変換部105と異なる動作をする。
【0073】
図14はノイズ低減処理の手順を示すフローチャート図の一例を示す。図14において図3と同一ステップには同一の符号を付しその説明は省略する。
【0074】
<S201>
ステップS104までの処理により、類似領域の平均値ベクトルμ、Σtotal 、及び、Σnoiseが算出されている。ステップS201では、最大分散方向算出部201が、類似領域の画素の本来の画素値の色空間における最大分散方向を算出する。具体的には次の手順で行う。
【0075】
図15は、ステップS201の詳細フローの一例を示す。
【0076】
まず、最大分散方向算出部201は、類似領域の画素の分散共分散行列Σtotalから、ノイズの分散共分散行列Σnoiseを減算する(S2011)。
Σtotal −Σnoise
最大分散方向算出部201は、得られた行列の固有値と固有ベクトルを求める(S2012)。分散共分散行列の固有ベクトルは、与えられたデータの分散がどちらの方向を向いているかを示す。また、固有値は、その方向に対してどの程度の大きさで分散が広がっているか、つまりどの程度ばらついているかを示す。したがって、分散共分散行列の固有値が最大になる固有ベクトルを求めることで、最大分散方向が得られる(S2013)。
【0077】
<S202>
ステップS202では、差分変換部202は、最大分散方向を特定方向eとして実施例1と同様の処理を行う。すなわち、差分変換部202は、注目画素ベクトルgから類似領域の平均値ベクトルμを引き、得られた差分ベクトルsを、最大分散方向に射影する。そして、類似領域の分散共分散行列Σtotalを用いて、射影した差分ベクトルsにノイズを低減する処理を施す。
【0078】
最大分散方向に差分ベクトルsを射影することで、RGB比率一定方向や輝度方向と異なり、元の画像の性質をより強く反映することができる。ノイズ低減による微細信号の劣化を低減することができる。また、エッジを保存しやすくできる。
【実施例4】
【0079】
本実施例では、ノイズ量、テクスチャ量の算出に画素値の分散共分散行列からではなく、空間周波数を用いるノイズ低減装置100について説明する。
【0080】
図16は、本実施例における機能ブロック図の一例を示す。図16において図2と同一部には同一の符号を付しその説明は省略する。図16では、平均・分散出部103の代わりに平均・周波数算出部301を、ノイズ分散推定部104の代わりにノイズ量算出部302を、を有する。また、差分変換部303は、差分変換部105と同等の作用を提供する。
【0081】
図17は、ノイズ低減処理の手順を示すフローチャート図の一例を示す。図17において図3と同一ステップには同一の符号を付しその説明は省略する。
【0082】
<S301>
ステップS301では、平均・周波数算出部301が類似領域の画素の平均値ベクトルμと特定周波数の強度ベクトルを算出する。平均の算出方法は実施例1と同じなので省略する。特定周波数の強度ベクトルを算出するため、平均・周波数算出部301は、類似領域の画素の特定周波数の強度を各色(RGB)成分のフィルタ応答により求める。
【0083】
図18は、フィルタサイズが5×5の場合のフィルタの一例を示す。なお、フィルタ演算には平均値演算に用いた画素のみを使用する。用いるフィルタを調整することで、処理対象画像31から所望の周波数成分を抽出できる。例えば、エッジを抽出するフィルタや平滑化するフィルタが知られている。図18は、高周波成分を抽出するフィルタの例である。
【0084】
処理対象画像31に含まれるノイズは空間周波数に依存する性質を持つ場合があることが知られている。実施例1のように画素値の平均からノイズを推定すると、処理対象画像31の全ての周波数成分からノイズを推定することになり、周波数成分によっては適切なノイズ低減でないことが生じることがある。
【0085】
これに対し、この空間周波数による強度調節の方法は、低減したいノイズの周波数に合わせてノイズ低減強度を調節する側面と、残したいテクスチャの周波数に合わせてノイズ低減強度を調節する側面、の二つ側面を有する。すなわち、ある周波数帯におけるノイズとテクスチャの量をノイズ低減強度に反映させることができる。分散によるノイズの推定が全周波数帯を対象とするのに対して、特定の周波数帯を対象とすることが可能となる。したがって、特定の周波数帯のノイズとテクスチャの量をノイズ低減強度に反映させればいいことが分かっている場合に効果的なノイズ低減処理を実現できる。
【0086】
<S302>
ステップS302では、ノイズ強度算出部302がノイズ強度を算出する。ノイズ強度を算出するため、類似領域の画素の平均からノイズの特定周波数の強度を返す関数を予め同定しておく。この関数は、図7と同様なので省略する。これを各色成分に適用してノイズの特定周波数の強度ベクトルを求める。
【0087】
<S303>
ステップS303では、差分変換部303が差分ベクトルsを変換する。以降の手順は実施例1と同様である。差分変換部303は特定方向eに差分ベクトルsを射影してから、その大きさを変える。特定方向eとしては、RGB比率一定方向、輝度方向、最大分散方向等とすることができる。
【0088】
差分変換部303は、ステップS302で算出された、ノイズの特定周波数の強度ベクトルを特定方向eに射影したベクトルの大きさをノイズ量とする。また、特定方向eに射影された差分ベクトルsからノイズの特定周波数の強度ベクトルを減算して得られたベクトルの大きさをテクスチャ量とする。このノイズ量とテクスチャ量を用いて、例えば、ステップS1503のようにして差分ベクトルsの射影ベクトルの大きさを変更する。なお、特定方向eをRGB比率一定方向とした場合は、実施例2のように、類似領域の平均値ベクトルμをη倍してもよい。
【0089】
本実施例のノイズ低減装置100は、フィルタにより特定周波数の成分を抽出しノイズの強度を算出するので、所望の周波数成分のノイズを低減しやすくできる。
【実施例5】
【0090】
これまでの実施例では、差分ベクトルsの射影方向を特定方向eに変換したが、本実施例では、差分ベクトルsを、特定方向eと特定方向eに垂直な方向のそれぞれのノイズ量に応じて変換するノイズ低減装置100について説明する。いわば、差分ベクトルsをRGB空間における二次元部分空間にて変換することになる。
【0091】
図19は、本実施例における機能ブロック図の一例を示す。図19において図2と同一部には同一の符号を付しその説明は省略する。図19では、差分変換部401を有する。
【0092】
図20は、ノイズ低減処理の手順を示すフローチャート図の一例を示す。図20において図3と同一ステップには同一の符号を付しその説明は省略する。
【0093】
<S401>
処理ステップS401では、差分変換部401が、差分ベクトルsの、特定方向eの分散と、特定方向eに垂直な方向e⊥の分散、それぞれついて大きさを変換する。
図21は、差分ベクトルsの二次元方向への射影を模式的に示す図の一例である。図21(a)に示すように、差分ベクトルsと特定方向eを含む平面が一意に定まる。また、ノイズ量はRGB成分毎に求められている。特定方向eの分散の求め方は実施例1と同じである。特定方向eに垂直な方向e⊥のノイズ量は、ノイズの分散共分散行列Σnoizeのうち特定方向eに垂直な方向e⊥の分散である。特定方向eに垂直な方向e⊥のテクスチャ量は、類似領域の分散共分散行列Σtotalからノイズの分散共分散行列Σnoizeを引いて得られる分散共分散行列の特定方向eに垂直な方向e⊥の分散である。
【0094】
図21(b)は、それぞれの方向において(テクスチャ量)/(ノイズ量+テクスチャ量)で大きさを変換した差分ベクトルsの一例を示す。特定方向eと、垂直な方向e⊥の2つの方向に対し、差分ベクトルsを各方向のノイズ量に応じて変換することで、特定方向eに近づける方向で変換できる。したがって、特定方向eと垂直な方向のテクスチャ量がノイズ量に比して大きければ、射影される差分ベクトルsに注目画素の情報を多く残すことができる。逆に、特定方向eと垂直な方向のテクスチャ量がノイズ量に比して小さければ、射影される差分ベクトルsを類似領域の平均値ベクトルμに近くすることでノイズを低減できる。なお、特定方向eはRGB比率一定方向、輝度方向、又は、最大分散方向のいずれでもよい。
【0095】
また、例えば特定方向eを輝度方向に取った場合、特定方向eに垂直な方向e⊥の差分ベクトルsの大きさを特定方向eの大きさよりも高い比率で縮小することができる。輝度方向と垂直な方向は人間の眼の感度が低いので、特定方向eに垂直方向の差分ベクトルsをその方向のノイズ量より大きく圧縮しても、人間の視覚が感じる変化は小さい。したがって、人間の視覚に重要な輝度成分を保存したまま、効率的にカラーノイズを低減することができる。
【0096】
本実施例のノイズ低減装置100によれば、差分ベクトルsをRGB空間における二次元部分空間にて変換することで、ノイズ量に応じて注目画素の情報を適切に圧縮しながら、カラー画像のノイズを低減できる。
【実施例6】
【0097】
実施例1では差分ベクトルsを特定方向eの一次元方向に射影し、実施例5では差分ベクトルsを二次元部分空間のベクトルに変換することでノイズを低減した。本実施例では、差分ベクトルsを三次元部分空間のベクトルに変換することでノイズを低減するノイズ低減装置100について説明する。すなわち、本実施例のノイズ低減装置100は、三次元の色空間(RGB色空間)におけるノイズの低減を一般化したものである。
【0098】
なお、本実施例の機能ブロック図は、図2又は図19と同様であり、特に差分変換部105又は差分変換部401を三次元空間の処理に拡張したものである。また、ノイズ低減の手順は図3又は図20と等しい。
【0099】
実施例1,5を一般化すると、差分ベクトルsを線形変換し、類似領域の平均値ベクトルに加算する処理となる。具体的には、類似領域の平均値ベクトルμ、注目画素値ベクトルg、線形変換の行列Wを用いて次の式で出力画素値fを算出する処理に相当する。
【0100】
【数6】
行列Wはたとえば次のように設定することができる。
【0101】
【数7】
Σnoiseはノイズの分散共分散行列、Σtextureはテクスチャの分散共分散行列である。テクスチャの分散共分散行列Σtextureは、類似領域の分散共分散行列Σtotal からノイズの分散共分散行列Σnoiseを引いたものである。したがって、式(9)から、行列Wは、三次元空間のそれぞれの方向で(テクスチャ量)/(ノイズ量+テクスチャ量)の演算を可能とする行列であることがわかる。
【0102】
式(8)は、差分ベクトル(g−μ)を三次元空間で線形変換し、平均値ベクトルμに加算する処理であるので、色空間における所望の3方向、それぞれについてノイズ量とテクスチャ量のバランスを勘案して、差分ベクトルsを三次元空間のベクトルに変換することができる。
【0103】
図22は、差分ベクトルsの三次元方向への射影を模式的に示す図の一例である。方向e1, e2, e3はWの固有ベクトルの方向である。注目画素gが含むノイズ量を求めておけば、図示するように、方向e1, e2, e3それぞれについて、(テクスチャ量)/(ノイズ量+テクスチャ量)に対応した変換が可能になる。変換後の各方向の分散から出力画素値のベクトルfが求められる。
【0104】
ノイズ量に対してテクスチャ量が多い方向では注目画素値に近い値となり、出力画素値に注目画素の情報を多く残すことができる。その逆に、ノイズ量に対してテクスチャ量が少ない方向では類似領域の平均に近い値となるので適切にノイズを低減できる。
【0105】
また、特定方向e0(単位ベクトル)の分散を相対的に多く残すことで、人間の知覚上有用な成分を選択的に多く残すことも可能である。そのためには線形変換の行列Wを次のように設定する。
【0106】
【数8】
ここでIは3×3の単位行列である。αは設計パラメータであり、特定方向e0の分散の強調の程度を指定する。αが大きいほど、特定方向e0の分散を強調しやすくできる。特定方向e0としてはこれまで説明したようにRGB比率一定方向、輝度方向、最大分散方向等を取ることができる。
【0107】
また、より設計の自由度を高く持たせて、3×3のパラメータ行列A, Bを用いて、次のように設定してもよい。
【0108】
【数9】
式(11)は、実施例1の<S105の変形例>で説明した、「α×(ノイズ量)/(ノイズ量+テクスチャ量)+β」と同等である。すなわち、AとBを調整することで、射影したベクトルに乗ずる係数も、ノイズ量が多くなれば小さくなり、テクスチャ量が多くなれば大きくなるような特性を持たせることができる。
【0109】
本実施例のノイズ低減装置100は、RGB空間のノイズの低減を、三次元部分空間に一般化して、各方向のノイズ量の大きさに応じて、注目画素の情報を多く残したり、類似領域の平均に近い値としたりすることができる。
【0110】
以上説明したように、本実施形態のノイズ低減装置100は、注目画素から類似度が閾値内の画素をフィルタ処理することで、エッジなど人間が画像を知覚する上で重要な要素を保存しつつノイズ低減を行うことができ、差分ベクトルを例えば平均値ベクトルμに加算することで振幅の小さい信号を保存でき、カラー画像における非加法性のノイズを推定してノイズ量とテクスチャ量に応じてフィルタ処理の出力を変換することで、信号値に大きさが依存する非加法性ノイズを、元の画像のテクスチャの量を勘案して適切に低減することができる。また、ノイズ低減処理による画素値変化を色空間の特定方向に制約しているため、色ずれを起こさず色ノイズを効果的に低減することができる。
【符号の説明】
【0111】
11 撮像部
20 ノイズ低減プログラム
24 記憶媒体
31 処理対象画像
32 ノイズ低減画像
100 ノイズ低減装置
101 画像取得部
102 注目画素選択部
103 平均・分散推定部
104 ノイズ分散推定部
105 差分変換部
106 加算部
107 画素選択済み判定部
【先行技術文献】
【特許文献】
【0112】
【特許文献1】特開2008−205737号公報
【特許文献2】特開2007−288439号公報
【特許文献3】特開2002−259965号公報
【非特許文献】
【0113】
【非特許文献1】原島 博, 小田島 薫, 鹿喰 善明, 宮川 洋, "ε-分離非線形ディジタルフィルタとその応用," 電子情報通信学会論文誌 A , Vol.J65-A, No.4, pp.297-304, 1982.
【非特許文献2】CC. Tomasi and R. Manduchi, "Bilateral Filtering for Gray and Color Images," Proc. Sixth Int'l Conf. Computer Vision, pp. 839-846, 1998.
【非特許文献3】Ce Liu, et al., "Automatic Estimation and Removal of Noise from a Single Image," IEEE Trans. Pattern Analysis and Machine Intelligence, vol. 30, no. 2, pp. 299-314, 2008.
【特許請求の範囲】
【請求項1】
カラー画像の注目画素の画素値に、注目画素を含む近傍領域の画素の画素値を用いてノイズ低減処理を施すノイズ低減装置であって、
所定の色空間において注目画素との類似度が閾値以上である類似領域に含まれる画素を前記近傍領域から抽出し、前記類似領域の画素の代表値ベクトルを算出する代表値ベクトル算出手段と、
注目画素の注目画素ベクトルと前記代表値ベクトルとの差分ベクトルを、前記色空間の特定方向に射影する差分射影手段と、
前記代表値ベクトルに、前記特定方向に射影された前記差分ベクトルを加えたベクトルの要素を、注目画素の画素値に置き換える画素値修正手段と、
を有することを特徴とするノイズ低減装置。
【請求項2】
カラー画像の注目画素の画素値に、注目画素を含む近傍領域の画素の画素値を用いてノイズ低減処理を施すノイズ低減装置であって、
所定の色空間において注目画素との類似度が閾値以上の類似領域に含まれる画素を前記近傍領域から抽出し、前記類似領域の画素の代表値ベクトルを算出する代表値ベクトル算出手段と、
前記代表値ベクトルに含まれるノイズ量及び注目画素の画素値に基づき、乗算数を算出する乗算数算出手段と、
前記代表値ベクトルに前記乗算数を乗算する乗算手段と、
前記乗算数が乗算された前記代表値ベクトルの要素を、注目画素の画素値に置き換える画素値修正手段と、
を有することを特徴とするノイズ低減装置。
【請求項3】
前記差分射影手段は、前記代表値ベクトルに含まれるノイズ量及びテクスチャ量に応じた変換係数を用いて、前記差分ベクトルに線形変換を施す、
ことを特徴とする請求項1記載のノイズ低減装置。
【請求項4】
前記差分射影手段は、前記特定方向に射影された前記差分ベクトルの大きさを、
前記代表値ベクトルに含まれるノイズ量及びテクスチャ量に応じて変更する、
ことを特徴とする請求項3記載のノイズ低減装置。
【請求項5】
前記差分射影手段は、前記代表値ベクトルに含まれるノイズ量及びテクスチャ量に応じて、前記差分ベクトルの長さを、前記特定方向、及び、前記特定方向に対し垂直方向、のそれぞれの方向おいて変更する、
ことを特徴とする請求項3記載のノイズ低減装置。
【請求項6】
前記特定方向は、前記代表値ベクトルの要素とRGB比率が一定のRGB比率一定定方向である、
ことを特徴とする請求項4又は5記載のノイズ低減装置。
【請求項7】
前記特定方向は、輝度方向である、
ことを特徴とする請求項4又は5記載のノイズ低減装置。
【請求項8】
前記類似領域の画素の画素値の最大分散方向を求める最大分散方向算出手段を有し、
前記特定方向を該最大分散方向とする、
ことを特徴とする請求項4又は5記載のノイズ低減装置。
【請求項9】
前記代表値ベクトルは、前記類似領域の画素の画素値の平均情報を有し、
画素値に含まれるノイズ量の分散と画素値の関係を登録したマップに基づき、前記平均情報に対応したノイズ量を前記代表値ベクトルに含まれるノイズ量と推定するノイズ量推定手段、
を有することを特徴とする請求項3記載のノイズ低減装置。
【請求項10】
前記ノイズ量推定手段は、前記代表値ベクトルの分散から前記類似領域のノイズ量の分散を減じてテクスチャ量を算出する、
ことを特徴とする請求項9記載のノイズ低減装置。
【請求項11】
前記ノイズ量推定手段は、前記類似領域の画素に周波数変換を施し、周波数変換された画素の前記平均情報に基づき前記代表値ベクトルに含まれるノイズ量を推定すると共に、前記テクスチャ量を算出する、
ことを特徴とする請求項10記載のノイズ低減装置。
【請求項12】
カラー画像の注目画素の画素値に、注目画素を含む近傍領域の画素の画素値を用いてノイズ低減処理を施すノイズ低減方法であって、
代表値ベクトル算出手段が、所定の色空間において注目画素との類似度が閾値以上の類似領域に含まれる画素を前記近傍領域から抽出し、前記類似領域の画素の代表値ベクトルを算出するステップと、
差分射影手段が、注目画素の注目画素ベクトルと前記代表値ベクトルとの差分ベクトルを、前記色空間の特定方向に射影するステップと、
画素値修正手段が、前記代表値ベクトルに、前記特定方向に射影された前記差分ベクトルを加えたベクトルの要素を、注目画素の画素値に置き換えるステップと、
を有することを特徴とするノイズ低減方法。
【請求項13】
カラー画像の注目画素の画素値に、注目画素を含む近傍領域の画素の画素値を用いてノイズ低減処理を施すノイズ低減方法であって、
代表値ベクトル算出手段が、所定の色空間において注目画素との類似度が閾値以上の類似領域に含まれる画素を前記近傍領域から抽出し、前記類似領域の画素の代表値ベクトルを算出するステップと、
乗算数算出手段が、前記代表値ベクトルに含まれるノイズ量及び注目画素の画素値に基づき乗算数を算出するステップと、
乗算手段が、前記代表値ベクトルに前記乗算数を乗算するステップと、
画素値修正手段が、前記乗算数が乗算された前記代表値ベクトルの要素を、注目画素の画素値に置き換えるステップと、
を有することを特徴とするノイズ低減方法。
【請求項14】
カラー画像の注目画素の画素値に、注目画素を含む近傍領域の画素の画素値を用いてノイズ低減処理を施すノイズ低減プログラムであって、
コンピュータに、
前記カラー画像を入力するステップと、
所定の色空間において注目画素との類似度が閾値以上の類似領域に含まれる画素を前記近傍領域から抽出し、前記類似領域の画素の代表値ベクトルを算出するステップと、
注目画素の注目画素ベクトルと前記代表値ベクトルとの差分ベクトルを、前記色空間の特定方向に射影するステップと、
前記代表値ベクトルに、前記特定方向に射影された前記差分ベクトルを加えたベクトルの要素を、注目画素の画素値に置き換えるステップと、
画素値が置き換えられた前記カラー画像を記憶装置に記憶するステップと、
を実行させるノイズ低減プログラム。
【請求項15】
カラー画像の注目画素の画素値に、注目画素を含む近傍領域の画素の画素値を用いてノイズ低減処理を施すノイズ低減プログラムであって、
コンピュータに、
前記カラー画像を入力するステップと、
所定の色空間において注目画素との類似度が閾値以上の類似領域に含まれる画素を前記近傍領域から抽出し、前記類似領域の画素の代表値ベクトルを算出するステップと、
前記代表値ベクトルに含まれるノイズ量及び注目画素の画素値に基づき乗算数を算出するステップと、
前記代表値ベクトルに前記乗算数を乗算するステップと、
前記乗算数が乗算された前記代表値ベクトルの要素を、注目画素の画素値に置き換えるステップと、
画素値が置き換えられた前記カラー画像を記憶装置に記憶するステップと、
を実行させるノイズ低減プログラム。
【請求項16】
請求項14又は15いずれか記載のノイズ低減プログラムを記憶した記憶媒体。
【請求項1】
カラー画像の注目画素の画素値に、注目画素を含む近傍領域の画素の画素値を用いてノイズ低減処理を施すノイズ低減装置であって、
所定の色空間において注目画素との類似度が閾値以上である類似領域に含まれる画素を前記近傍領域から抽出し、前記類似領域の画素の代表値ベクトルを算出する代表値ベクトル算出手段と、
注目画素の注目画素ベクトルと前記代表値ベクトルとの差分ベクトルを、前記色空間の特定方向に射影する差分射影手段と、
前記代表値ベクトルに、前記特定方向に射影された前記差分ベクトルを加えたベクトルの要素を、注目画素の画素値に置き換える画素値修正手段と、
を有することを特徴とするノイズ低減装置。
【請求項2】
カラー画像の注目画素の画素値に、注目画素を含む近傍領域の画素の画素値を用いてノイズ低減処理を施すノイズ低減装置であって、
所定の色空間において注目画素との類似度が閾値以上の類似領域に含まれる画素を前記近傍領域から抽出し、前記類似領域の画素の代表値ベクトルを算出する代表値ベクトル算出手段と、
前記代表値ベクトルに含まれるノイズ量及び注目画素の画素値に基づき、乗算数を算出する乗算数算出手段と、
前記代表値ベクトルに前記乗算数を乗算する乗算手段と、
前記乗算数が乗算された前記代表値ベクトルの要素を、注目画素の画素値に置き換える画素値修正手段と、
を有することを特徴とするノイズ低減装置。
【請求項3】
前記差分射影手段は、前記代表値ベクトルに含まれるノイズ量及びテクスチャ量に応じた変換係数を用いて、前記差分ベクトルに線形変換を施す、
ことを特徴とする請求項1記載のノイズ低減装置。
【請求項4】
前記差分射影手段は、前記特定方向に射影された前記差分ベクトルの大きさを、
前記代表値ベクトルに含まれるノイズ量及びテクスチャ量に応じて変更する、
ことを特徴とする請求項3記載のノイズ低減装置。
【請求項5】
前記差分射影手段は、前記代表値ベクトルに含まれるノイズ量及びテクスチャ量に応じて、前記差分ベクトルの長さを、前記特定方向、及び、前記特定方向に対し垂直方向、のそれぞれの方向おいて変更する、
ことを特徴とする請求項3記載のノイズ低減装置。
【請求項6】
前記特定方向は、前記代表値ベクトルの要素とRGB比率が一定のRGB比率一定定方向である、
ことを特徴とする請求項4又は5記載のノイズ低減装置。
【請求項7】
前記特定方向は、輝度方向である、
ことを特徴とする請求項4又は5記載のノイズ低減装置。
【請求項8】
前記類似領域の画素の画素値の最大分散方向を求める最大分散方向算出手段を有し、
前記特定方向を該最大分散方向とする、
ことを特徴とする請求項4又は5記載のノイズ低減装置。
【請求項9】
前記代表値ベクトルは、前記類似領域の画素の画素値の平均情報を有し、
画素値に含まれるノイズ量の分散と画素値の関係を登録したマップに基づき、前記平均情報に対応したノイズ量を前記代表値ベクトルに含まれるノイズ量と推定するノイズ量推定手段、
を有することを特徴とする請求項3記載のノイズ低減装置。
【請求項10】
前記ノイズ量推定手段は、前記代表値ベクトルの分散から前記類似領域のノイズ量の分散を減じてテクスチャ量を算出する、
ことを特徴とする請求項9記載のノイズ低減装置。
【請求項11】
前記ノイズ量推定手段は、前記類似領域の画素に周波数変換を施し、周波数変換された画素の前記平均情報に基づき前記代表値ベクトルに含まれるノイズ量を推定すると共に、前記テクスチャ量を算出する、
ことを特徴とする請求項10記載のノイズ低減装置。
【請求項12】
カラー画像の注目画素の画素値に、注目画素を含む近傍領域の画素の画素値を用いてノイズ低減処理を施すノイズ低減方法であって、
代表値ベクトル算出手段が、所定の色空間において注目画素との類似度が閾値以上の類似領域に含まれる画素を前記近傍領域から抽出し、前記類似領域の画素の代表値ベクトルを算出するステップと、
差分射影手段が、注目画素の注目画素ベクトルと前記代表値ベクトルとの差分ベクトルを、前記色空間の特定方向に射影するステップと、
画素値修正手段が、前記代表値ベクトルに、前記特定方向に射影された前記差分ベクトルを加えたベクトルの要素を、注目画素の画素値に置き換えるステップと、
を有することを特徴とするノイズ低減方法。
【請求項13】
カラー画像の注目画素の画素値に、注目画素を含む近傍領域の画素の画素値を用いてノイズ低減処理を施すノイズ低減方法であって、
代表値ベクトル算出手段が、所定の色空間において注目画素との類似度が閾値以上の類似領域に含まれる画素を前記近傍領域から抽出し、前記類似領域の画素の代表値ベクトルを算出するステップと、
乗算数算出手段が、前記代表値ベクトルに含まれるノイズ量及び注目画素の画素値に基づき乗算数を算出するステップと、
乗算手段が、前記代表値ベクトルに前記乗算数を乗算するステップと、
画素値修正手段が、前記乗算数が乗算された前記代表値ベクトルの要素を、注目画素の画素値に置き換えるステップと、
を有することを特徴とするノイズ低減方法。
【請求項14】
カラー画像の注目画素の画素値に、注目画素を含む近傍領域の画素の画素値を用いてノイズ低減処理を施すノイズ低減プログラムであって、
コンピュータに、
前記カラー画像を入力するステップと、
所定の色空間において注目画素との類似度が閾値以上の類似領域に含まれる画素を前記近傍領域から抽出し、前記類似領域の画素の代表値ベクトルを算出するステップと、
注目画素の注目画素ベクトルと前記代表値ベクトルとの差分ベクトルを、前記色空間の特定方向に射影するステップと、
前記代表値ベクトルに、前記特定方向に射影された前記差分ベクトルを加えたベクトルの要素を、注目画素の画素値に置き換えるステップと、
画素値が置き換えられた前記カラー画像を記憶装置に記憶するステップと、
を実行させるノイズ低減プログラム。
【請求項15】
カラー画像の注目画素の画素値に、注目画素を含む近傍領域の画素の画素値を用いてノイズ低減処理を施すノイズ低減プログラムであって、
コンピュータに、
前記カラー画像を入力するステップと、
所定の色空間において注目画素との類似度が閾値以上の類似領域に含まれる画素を前記近傍領域から抽出し、前記類似領域の画素の代表値ベクトルを算出するステップと、
前記代表値ベクトルに含まれるノイズ量及び注目画素の画素値に基づき乗算数を算出するステップと、
前記代表値ベクトルに前記乗算数を乗算するステップと、
前記乗算数が乗算された前記代表値ベクトルの要素を、注目画素の画素値に置き換えるステップと、
画素値が置き換えられた前記カラー画像を記憶装置に記憶するステップと、
を実行させるノイズ低減プログラム。
【請求項16】
請求項14又は15いずれか記載のノイズ低減プログラムを記憶した記憶媒体。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【公開番号】特開2010−218110(P2010−218110A)
【公開日】平成22年9月30日(2010.9.30)
【国際特許分類】
【出願番号】特願2009−62876(P2009−62876)
【出願日】平成21年3月16日(2009.3.16)
【公序良俗違反の表示】
(特許庁注:以下のものは登録商標)
1.EEPROM
【出願人】(000006747)株式会社リコー (37,907)
【Fターム(参考)】
【公開日】平成22年9月30日(2010.9.30)
【国際特許分類】
【出願日】平成21年3月16日(2009.3.16)
【公序良俗違反の表示】
(特許庁注:以下のものは登録商標)
1.EEPROM
【出願人】(000006747)株式会社リコー (37,907)
【Fターム(参考)】
[ Back to top ]