(57)【要約】【課題】写真画像データのノイズ変動幅を特定し、効果的にノイズ抑制可能な写真画像処理の方法、プログラム、及び装置を提供する。【解決手段】平坦画像領域から回帰直線を求めるノイズ変動幅特定ステップと、回帰直線から任意の注目画素に対するノイズ変動幅とノイズ抑制強度を算出するノイズ抑制強度算出ステップと、注目画素と周辺画素の濃度差分値とノイズ抑制強度から加重平均フィルタのフィルタ係数を算出する濃度対応フィルタ係数生成ステップと、濃度成分画像データを加重平均フィルタ処理する濃度対応ノイズ抑制処理ステップと、濃度成分画像データを移動平均フィルタ処理する移動平均フィルタ処理ステップと、移動平均画素のノイズ変動幅を求める移動平均画素ノイズ変動幅特定ステップと、移動平均画素のノイズ変動幅から加重平均画素を補正するグラデーション補正処理ステップとを含む写真画像処理方法。【選択図】図１０
2010108008.tif 000002
【特許請求の範囲】
【請求項１】
加重平均フィルタ処理により写真画像データに含まれる粒状のノイズを抑制する写真画像処理方法であって、
前記写真画像データを構成する濃度成分画像データから複数の平坦画像領域を探索し、各平坦画像領域の平均濃度値を説明変数とし最大濃度値及び最小濃度値の夫々を従属変数とする二本の回帰直線により前記写真画像データに含まれる粒状のノイズの変動幅を求めるノイズ変動幅特定ステップと、
前記ノイズ変動幅特定ステップで得られた各回帰直線に基づいて、任意の注目画素に対するノイズ変動幅を求め、当該ノイズ変動幅からノイズ抑制強度を算出するノイズ抑制強度算出ステップと、
当該注目画素と周辺画素の濃度差分値と、前記ノイズ抑制強度算出ステップで得られたノイズ抑制強度に基づいて、濃度成分画像データに対応したサイズの加重平均フィルタのフィルタ係数を動的に算出する濃度対応フィルタ係数生成ステップと、
前記濃度対応フィルタ係数生成ステップで得られた各注目画素に対応するフィルタ係数に基づいて、前記濃度成分画像データを加重平均フィルタ処理する濃度対応ノイズ抑制処理ステップと、
前記濃度成分画像データを前記濃度成分画像データに対応したサイズの移動平均フィルタでフィルタ処理する移動平均フィルタ処理ステップと、
前記ノイズ変動幅特定ステップで得られた各回帰直線に基づいて、前記移動平均フィルタ処理ステップで得られた各注目画素領域内の移動平均画素のノイズ変動幅を求める移動平均画素ノイズ変動幅特定ステップと、
前記移動平均画素ノイズ変動幅特定ステップで得られた移動平均画素のノイズ変動幅に基づいて、前記濃度対応ノイズ抑制処理ステップで得られた各加重平均画素を補正するグラデーション補正処理ステップと、
を含む写真画像処理方法。
【請求項２】
前記グラデーション補正処理ステップは、
各注目画素領域の濃度最大値と濃度最小値の濃度差分値を求める濃度差分値算出ステップと、
前記移動平均画素ノイズ変動幅特定ステップで得られた移動平均画素のノイズ変動幅と前記濃度差分値算出ステップで得られた濃度差分値の比率に基づいて、各注目領域の平坦度を算出する平坦度算出ステップと、
を含み、前記平坦度算出ステップで得られた平坦度を重み係数として、前記濃度対応ノイズ抑制処理ステップで得られた各加重平均画素と対応する移動平均画素とを重み付け加算してグラデーション再生濃度画像データを算出する請求項１記載の写真画像処理方法。
【請求項３】
前記グラデーション補正処理ステップは、
各注目画素領域の濃度値の標準偏差を求める濃度標準偏差算出ステップと、
前記移動平均画素ノイズ変動幅特定ステップで得られた移動平均画素のノイズ変動幅に基づいて、各注目画素領域内の移動平均画素の濃度値の仮想標準偏差を算出し、前記仮想標準偏差と前記濃度標準偏差算出ステップで得られた濃度値の標準偏差との比率に基づいて、各注目領域の平坦度を算出する平坦度算出ステップと、
を含み、前記平坦度算出ステップで得られた平坦度を重み係数として、前記濃度対応ノイズ抑制処理ステップで得られた各加重平均画素と対応する移動平均画素とを重み付け加算してグラデーション再生濃度画像データを算出する請求項１記載の写真画像処理方法。
【請求項４】
前記フィルタ係数は、前記ノイズ抑制強度を仮想分散値とし、前記注目画素と周辺画素の濃度差分値を変数とする正規分布関数に基づいて算出されるものである請求項１から３何れかに記載の写真画像処理方法。
【請求項５】
前記濃度対応ノイズ抑制処理ステップは、
前記ノイズ抑制強度を仮想分散値とし、前記濃度対応ノイズ抑制処理ステップで得られた各加重平均画素と対応する注目画素との濃度差分値を変数とする正規分布関数に基づいて重み係数を算出し、当該重み係数で各加重平均画素と対応する注目画素とを重み付け加算して補正濃度画像データを算出する補正処理ステップをさらに含む請求項１から４の何れかに記載の写真画像処理方法。
【請求項６】
前記ノイズ変動幅特定ステップは、
前記写真画像データから濃度成分画像データを抽出する濃度成分抽出ステップと、
前記濃度成分抽出ステップで抽出された前記濃度成分画像データを所定サイズの画素ブロックに分割して、画素ブロック毎に平均濃度値、最大濃度値、最小濃度値、濃度ばらつきを含むブロック特性値を算出するブロック特性値算出ステップと、
濃度範囲を所定間隔で複数の区間に区分し、前記平均濃度値に基づいて各画素ブロックを対応する区分に割り付ける画素ブロック割付ステップと、
各区分に割り付けられた画素ブロックから前記濃度ばらつきが最小となる画素ブロックを代表画素ブロックとして抽出する画素ブロック抽出ステップと、
を備え、前記画素ブロック抽出ステップで抽出された複数の代表画素ブロックの各平均濃度値を説明変数とし、各最大濃度値及び各最小濃度値の夫々を従属変数とする二本の回帰直線を算出するものである請求項１から５の何れかに記載の写真画像処理方法。
【請求項７】
加重平均フィルタ処理により写真画像データに含まれる粒状のノイズを抑制する写真画像処理装置であって、
前記写真画像データを構成する濃度成分画像データから複数の平坦画像領域を探索し、各平坦画像領域の平均濃度値を説明変数とし最大濃度値及び最小濃度値の夫々を従属変数とする二本の回帰直線により前記写真画像データに含まれる粒状のノイズの変動幅を求めるノイズ変動幅特定部と、
前記ノイズ変動幅特定部で得られた各回帰直線に基づいて、任意の注目画素に対するノイズ変動幅を求め、当該ノイズ変動幅からノイズ抑制強度を算出するノイズ抑制強度算出部と、
当該注目画素と周辺画素の濃度差分値と、前記ノイズ抑制強度算出部で得られたノイズ抑制強度に基づいて、濃度成分画像データに対応したサイズの加重平均フィルタのフィルタ係数を動的に算出する濃度対応フィルタ係数生成部と、
前記濃度対応フィルタ係数生成部で得られた各注目画素に対応するフィルタ係数に基づいて、前記濃度成分画像データを加重平均フィルタ処理する濃度対応ノイズ抑制処理部と、
前記濃度成分画像データを前記濃度成分画像データに対応したサイズの移動平均フィルタでフィルタ処理する移動平均フィルタ処理部と、
前記ノイズ変動幅特定部で得られた各回帰直線に基づいて、前記移動平均フィルタ処理部で得られた各注目画素領域内の移動平均画素のノイズ変動幅を求める移動平均画素ノイズ変動幅特定部と、
前記移動平均画素ノイズ変動幅特定部で得られた移動平均画素のノイズ変動幅に基づいて、前記濃度対応ノイズ抑制処理部で得られた各加重平均画素を補正するグラデーション補正処理部と、
を含む写真画像処理装置。
【請求項８】
前記グラデーション補正処理部は、
各注目画素領域の濃度最大値と濃度最小値の濃度差分値を求める濃度差分値算出部と、
前記移動平均画素ノイズ変動幅特定部で得られた移動平均画素のノイズ変動幅と前記濃度差分値算出部で得られた濃度差分値の比率に基づいて、各注目領域の平坦度を算出する平坦度算出部と、
を含み、前記平坦度算出部で得られた平坦度を重み係数として、前記濃度対応ノイズ抑制処理部で得られた各加重平均画素と対応する移動平均画素とを重み付け加算してグラデーション再生濃度画像データを算出する請求項７記載の写真画像処理装置。
【請求項９】
前記グラデーション補正処理部は、
各注目画素領域の濃度値の標準偏差を求める濃度標準偏差算出部と、
前記移動平均画素ノイズ変動幅特定部で得られた移動平均画素のノイズ変動幅に基づいて、各注目画素領域内の移動平均画素の濃度値の仮想標準偏差を算出し、前記仮想標準偏差と前記濃度標準偏差算出部で得られた濃度値の標準偏差との比率に基づいて、各注目領域の平坦度を算出する平坦度算出部と、
を含み、前記平坦度算出部で得られた平坦度を重み係数として、前記濃度対応ノイズ抑制処理部で得られた各加重平均画素と対応する移動平均画素とを重み付け加算してグラデーション再生濃度画像データを算出する請求項７記載の写真画像処理装置。
【請求項１０】
加重平均フィルタ処理により写真画像データに含まれる粒状のノイズを抑制する写真画像処理プログラムであって、
前記写真画像データを構成する濃度成分画像データから複数の平坦画像領域を探索し、各平坦画像領域の平均濃度値を説明変数とし最大濃度値及び最小濃度値の夫々を従属変数とする二本の回帰直線により前記写真画像データに含まれる粒状のノイズの変動幅を求めるノイズ変動幅特定ステップと、
前記ノイズ変動幅特定ステップで得られた各回帰直線に基づいて、任意の注目画素に対するノイズ変動幅を求め、当該ノイズ変動幅からノイズ抑制強度を算出するノイズ抑制強度算出ステップと、
当該注目画素と周辺画素の濃度差分値と、前記ノイズ抑制強度算出ステップで得られたノイズ抑制強度に基づいて、濃度成分画像データに対応したサイズの加重平均フィルタのフィルタ係数を動的に算出する濃度対応フィルタ係数生成ステップと、
前記濃度対応フィルタ係数生成ステップで得られた各注目画素に対応するフィルタ係数に基づいて、前記濃度成分画像データを加重平均フィルタ処理する濃度対応ノイズ抑制処理ステップと、
前記濃度成分画像データを前記濃度成分画像データに対応したサイズの移動平均フィルタでフィルタ処理する移動平均フィルタ処理ステップと、
前記ノイズ変動幅特定ステップで得られた各回帰直線に基づいて、前記移動平均フィルタ処理ステップで得られた各注目画素領域内の移動平均画素のノイズ変動幅を求める移動平均画素ノイズ変動幅特定ステップと、
前記移動平均画素ノイズ変動幅特定ステップで得られた移動平均画素のノイズ変動幅に基づいて、前記濃度対応ノイズ抑制処理ステップで得られた各加重平均画素を補正するグラデーション補正処理ステップと、
をコンピュータに実行させる写真画像処理プログラム。
【請求項１１】
前記グラデーション補正処理ステップは、
各注目画素領域の濃度最大値と濃度最小値の濃度差分値を求める濃度差分値算出ステップと、
前記移動平均画素ノイズ変動幅特定ステップで得られた移動平均画素のノイズ変動幅と前記濃度差分値算出ステップで得られた濃度差分値の比率に基づいて、各注目領域の平坦度を算出する平坦度算出ステップと、
を含み、前記平坦度算出ステップで得られた平坦度を重み係数として、前記濃度対応ノイズ抑制処理ステップで得られた各加重平均画素と対応する移動平均画素とを重み付け加算してグラデーション再生濃度画像データを算出する請求項１０記載の写真画像処理プログラム。
【請求項１２】
前記グラデーション補正処理ステップは、
各注目画素領域の濃度値の標準偏差を求める濃度標準偏差算出ステップと、
前記移動平均画素ノイズ変動幅特定ステップで得られた移動平均画素のノイズ変動幅に基づいて、各注目画素領域内の移動平均画素の濃度値の仮想標準偏差を算出し、前記仮想標準偏差と前記濃度標準偏差算出ステップで得られた濃度値の標準偏差との比率に基づいて、各注目領域の平坦度を算出する平坦度算出ステップと、
を含み、前記平坦度算出ステップで得られた平坦度を重み係数として、前記濃度対応ノイズ抑制処理ステップで得られた各加重平均画素と対応する移動平均画素とを重み付け加算してグラデーション再生濃度画像データを算出する請求項１０記載の写真画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、写真画像データに含まれる粒状のノイズを抑制する写真画像処理方法、写真画像処理プログラム、及び写真画像処理装置に関する。
【背景技術】
【０００２】
デジタルカメラで撮影された画像を観察すると、画像に細かな粒状のノイズが重畳されていることが観察される。特に画像の暗部でざらつくようなノイズが目立ち、暗所でストロボ発光無しで撮影された画像には同様のノイズが画像全体に現れる。
【０００３】
これらのノイズは、撮像素子のリセットノイズや暗電流ノイズ等、素子の特性等に起因して発生するもので、撮影環境が悪い、つまり被写体の照度が低いほどノイズ成分が勝り、さらにゲイン調整によって明るさを補正する場合にノイズが強調されるようになる。
【０００４】
同様に、光学カメラによって撮影されたフィルム画像にも色素雲の大きさに基づいた粒状ムラが存在し、フィルムスキャナを介してデジタル画像に変換した写真画像データには細かな粒状のノイズが重畳されている。これらの粒状のノイズは、フィルム粒子が大きい感度の高いフィルムに撮像されたとき、あるいは、写真画像を拡大処理したときに顕著に現れる。
【０００５】
これらの粒状のノイズを抑制するために、所定のフィルタサイズのメジアンフィルタやガウシアンフィルタによるフィルタ処理、つまり、ある画素についてその周囲の画素を加味して求めた所定の平均値に置換する処理を画像全体に施すような処理があるが、画像の輪郭部分にメリハリを与える鮮鋭化処理の後にこのようなぼかし処理を行なうと、輪郭までぼやけてしまうので鮮鋭化した意味が無くなり、上述のぼかし処理を行った後に鮮鋭化処理する場合には画像中の細かい構造（ディテール）が消失するといった不都合がある。
【０００６】
そこで、写真画像データを所定画素数でブロック分割して、各領域の平均値と分散値を求め、分散値が小さい平坦な画像領域では強い平滑化処理を施し、分散値が大きい輪郭を含む画像領域では弱い平滑化処理を施すことが考えられるが、この方法によれば、輪郭を含む画像領域がブロック内に存在すると、一律に弱い平滑化処理が実行されるため、輪郭を中心とする所定幅の帯状の領域でノイズ抑制効果が低減されることになり、あまり好ましい結果が得られなかった。
【０００７】
特許文献１には、フィルムスキャナを介してデジタル画像に変換した写真画像データに加重平均フィルタ処理を行なう画像処理方法として、写真画像データを構成するＲＧＢ画素成分毎に、注目画素を順次設定するステップと、注目画素を中心とした加重平均フィルタリング処理のための演算領域内に位置する各周辺画素と注目画素との間の画素値の差分値を算定するステップと、差分値に応じて加重平均フィルタリング処理に用いられる重み係数を決定するステップと、重み係数を用いて加重平均フィルタリング処理を実行して注目画素の補正画素値を求めるステップを備えた方法が提案されている。
【０００８】
この方法によれば、注目画素が平坦領域に位置する場合と輪郭領域に位置する場合とでは、注目画素とその周辺画素の画素値に所定の関係があることから、特定の差分値に対して強い影響力を与えるような重み係数を設定することで、平坦領域はもとより輪郭及びその付近においても十分に粒状ノイズを抑制できるとともに写真がもつ細かいディテールと質感が維持される粒状ノイズ抑制が実現可能となる。
【特許文献１】特開２００６−３０２０２３号公報
【発明の開示】
【発明が解決しようとする課題】
【０００９】
しかし、粒状ノイズを抑制する場合には、原画像データの特性、つまり、原画像データの画素数や、原画像に含まれている粒状ノイズのレベル等に基づいて適切なサイズのフィルタやフィルタ係数を設定する必要があるが、上述した従来の技術では、その程度を把握することができず、一律のフィルタ係数やフィルタサイズで粒状抑制処理がなされていたために、得られた画像の品質が異なるという問題があった。
【００１０】
例えば、粒状ノイズのレベルが相対的に低い原画像に対して、一律の加重平均フィルタ処理を施すと、平坦部がのっぺりした抑揚の無い画像に仕上がり、粒状ノイズのレベルが相対的に高い原画像に対して、一律の加重平均フィルタ処理を施すと、十分に粒状ノイズを抑制することができなかった。
【００１１】
デジタルカメラで撮影された写真画像データは、通常、ＪＰＥＧ画像に圧縮される際に、色差データが１／２または１／４に低減された後に直交変換されるため、輝度データと比べて色差データのデータ量が減少していることと相俟って、色ノイズが輝度ノイズより大きく現れやすく、上述した特許文献１に記載された色成分毎に加重平均フィルタ処理する方法では、輝度ノイズと色ノイズを適切に除去できないという問題もあった。
【００１２】
この場合、画像データの輝度成分と色成分を分離して、各成分に含まれるノイズを除去することも考えられるが、人間の目が細かな輝度の違いには敏感であるのに対して、色相と彩度の細かな違いには比較的気付かないという特性から、輝度成分に対するノイズ抑制処理の過不足が画質に直接的に影響することを考慮しなければならない困難さがあった。
【００１３】
本発明は、上述の従来欠点に鑑み、写真画像データからノイズの変動幅を特定し、効果的にノイズを抑制できる写真画像処理方法、写真画像処理プログラム、及び写真画像処理装置を提供する点にある。
【課題を解決するための手段】
【００１４】
上述の目的を達成するため、本発明による写真画像処理方法の第一の特徴構成は、特許請求の範囲の書類の請求項１に記載した通り、加重平均フィルタ処理により写真画像データに含まれる粒状のノイズを抑制する写真画像処理方法であって、前記写真画像データを構成する濃度成分画像データから複数の平坦画像領域を探索し、各平坦画像領域の平均濃度値を説明変数とし最大濃度値及び最小濃度値の夫々を従属変数とする二本の回帰直線により前記写真画像データに含まれる粒状のノイズの変動幅を求めるノイズ変動幅特定ステップと、前記ノイズ変動幅特定ステップで得られた各回帰直線に基づいて、任意の注目画素に対するノイズ変動幅を求め、当該ノイズ変動幅からノイズ抑制強度を算出するノイズ抑制強度算出ステップと、当該注目画素と周辺画素の濃度差分値と、前記ノイズ抑制強度算出ステップで得られたノイズ抑制強度に基づいて、濃度成分画像データに対応したサイズの加重平均フィルタのフィルタ係数を動的に算出する濃度対応フィルタ係数生成ステップと、前記濃度対応フィルタ係数生成ステップで得られた各注目画素に対応するフィルタ係数に基づいて、前記濃度成分画像データを加重平均フィルタ処理する濃度対応ノイズ抑制処理ステップと、前記濃度成分画像データを前記濃度成分画像データに対応したサイズの移動平均フィルタでフィルタ処理する移動平均フィルタ処理ステップと、前記ノイズ変動幅特定ステップで得られた各回帰直線に基づいて、前記移動平均フィルタ処理ステップで得られた各注目画素領域内の移動平均画素のノイズ変動幅を求める移動平均画素ノイズ変動幅特定ステップと、前記移動平均画素ノイズ変動幅特定ステップで得られた移動平均画素のノイズ変動幅に基づいて、前記濃度対応ノイズ抑制処理ステップで得られた各加重平均画素を補正するグラデーション補正処理ステップと、を含む点にある。
【００１５】
ノイズ抑制強度算出ステップでは、ノイズ変動幅特定ステップで得られた、平坦画像領域から算出されたノイズの変動幅に基づいて、任意の注目画素に対するノイズの変動幅が算出され、当該ノイズ変動幅からノイズ抑制強度が算出されるため、注目画素とその周辺画素の濃度差が大きな輪郭画像領域や絵柄が細かい画像領域等において、画像の内容を適切に表現する上で必要となる画素が、ノイズとして認識されないよう制限されたノイズの抑制強度を算出できるようになる。
【００１６】
濃度対応ノイズ抑制処理ステップでは、ノイズ抑制強度算出ステップ及び濃度対応フィルタ係数生成ステップを介して得られた加重平均フィルタを利用して、加重平均フィルタ処理を実行するため、ノイズの変動幅を超えた加重平均フィルタ処理を実行することにより、元の画像に含まれていた絵柄や輪郭までノイズとして抑制して消してしまう不具合を軽減できるようになる。
【００１７】
グラデーション補正処理ステップでは、移動平均フィルタ処理ステップと移動平均画素ノイズ変動幅特定ステップを介して、移動平均フィルタ処理によりノイズ量に依存しない一律的に平滑化された画像における移動平均画素のノイズ変動幅に基づいて、濃度対応ノイズ抑制処理ステップで得られた各加重平均画素を補正するため、移動平均画素のノイズ変動幅を超える程度にまで強度にノイズ抑制することを回避できるようになる。
【００１８】
同第二の特徴構成は、同請求項２に記載した通り、上述の第一特徴構成に加えて、前記グラデーション補正処理ステップは、各注目画素領域の濃度最大値と濃度最小値の濃度差分値を求める濃度差分値算出ステップと、前記移動平均画素ノイズ変動幅特定ステップで得られた移動平均画素のノイズ変動幅と前記濃度差分値算出ステップで得られた濃度差分値の比率に基づいて、各注目領域の平坦度を算出する平坦度算出ステップと、を含み、前記平坦度算出ステップで得られた平坦度を重み係数として、前記濃度対応ノイズ抑制処理ステップで得られた各加重平均画素と対応する移動平均画素とを重み付け加算してグラデーション再生濃度画像データを算出する点にある。
【００１９】
上述の構成によれば、平坦度算出ステップでは、移動平均画素のノイズ変動幅と濃度差分値算出ステップから得られる注目画素の濃度差分値は、同サイズの画像領域における移動平均画素及び注目画素の最大濃度値と最小濃度値の濃度差分値を示しており、当該濃度差分値を対比させることにより、移動平均画素のノイズ変動幅を基準として平坦度が算出される。
【００２０】
したがって、グラデーション補正処理ステップでは、平坦度算出ステップから得られた平坦度を重み係数として、濃度対応ノイズ抑制処理ステップで得られた各加重平均画素と対応する移動平均画素とを重み付け加算するため、移動平均フィルタ処理によりノイズ量に依存しない一律に平滑化された画像の濃度値の分布を基準として、平坦画像領域とそうではない画像領域を分離して、適切にノイズ抑制することができるようになる。
【００２１】
同第三の特徴構成は、同請求項３に記載した通り、上述の第一特徴構成に加えて、前記グラデーション補正処理ステップは、各注目画素領域の濃度値の標準偏差を求める濃度標準偏差算出ステップと、前記移動平均画素ノイズ変動幅特定ステップで得られた移動平均画素のノイズ変動幅に基づいて、各注目画素領域内の移動平均画素の濃度値の仮想標準偏差を算出し、前記仮想標準偏差と前記濃度標準偏差算出ステップで得られた濃度値の標準偏差との比率に基づいて、各注目領域の平坦度を算出する平坦度算出ステップと、を含み、前記平坦度算出ステップで得られた平坦度を重み係数として、前記濃度対応ノイズ抑制処理ステップで得られた各加重平均画素と対応する移動平均画素とを重み付け加算してグラデーション再生濃度画像データを算出する点にある。
【００２２】
上述の構成によれば、平坦度算出ステップでは、移動平均画素のノイズ変動幅に基づいて算出された移動平均画素の濃度値の仮想標準偏差と、濃度標準偏差算出ステップから得られる注目画素の濃度値の標準偏差は、同サイズの画像領域における移動平均画素及び注目画素の濃度値の標準偏差を示しており、当該濃度値の標準偏差を対比させることにより、移動平均画素のノイズ変動幅を基準として平坦度が算出される。
【００２３】
したがって、グラデーション補正処理ステップでは、平坦度算出ステップから得られた平坦度を重み係数として、濃度対応ノイズ抑制処理ステップで得られた各加重平均画素と対応する移動平均画素とを重み付け加算するため、移動平均フィルタ処理によりノイズ量に依存しない一律に平滑化された画像の濃度値の分布をさらに滑らかに仮想し、当該仮想した濃度値の分布を基準として、平坦画像領域とそうではない画像領域を滑らかに分離して、適切にノイズ抑制することができるようになる。
【００２４】
また、上述の第二の特徴構成の濃度差分値算出ステップでは、各注目画素領域の濃度最大値と濃度最小値を特定する場合に、各注目画素領域の全画素に対して、当該特定する処理中の濃度最大値及び濃度最小値と当該処理対象の画素の濃度値の大小を比較し、当該比較結果に基づいて濃度最大値及び濃度最小値を逐次更新する処理が必要である。さらに、濃度差分値算出ステップでは、当該逐次更新処理により算出された濃度最大値と濃度最小値を利用して、濃度差分値を算出できるようになる。しかし、当該第三の特徴構成の濃度標準偏差算出ステップでは、各注目画素領域の全画素に対して、処理対象の画素の濃度値を累積的に加算及び乗算する処理を実行して、濃度値の標準偏差を算出することができる。したがって、第二の特徴構成と比較して、処理速度を向上することができるようになる。
【００２５】
同第四の特徴構成は、同請求項４に記載した通り、上述の第一から第三の何れかの特徴構成に加えて、前記フィルタ係数は、前記ノイズ抑制強度を仮想分散値とし、前記注目画素と周辺画素の濃度差分値を変数とする正規分布関数に基づいて算出されるものである点にある。
【００２６】
上述の構成によれば、濃度対応フィルタ係数生成ステップにより、ノイズ抑制強度算出ステップで得られたノイズ抑制強度を標準偏差とする正規分布関数に基づいてフィルタ係数が算出されるため、滑らかな特性を持つフィルタ係数を算出することができるようになる。
【００２７】
同第五の特徴構成は、同請求項５に記載した通り、上述の第一から第四の何れかの特徴構成に加えて、前記濃度対応ノイズ抑制処理ステップは、前記ノイズ抑制強度を仮想分散値とし、前記濃度対応ノイズ抑制処理ステップで得られた各加重平均画素と対応する注目画素との濃度差分値を変数とする正規分布関数に基づいて重み係数を算出し、当該重み係数で各加重平均画素と対応する注目画素とを重み付け加算して補正濃度画像データを算出する補正処理ステップをさらに含む点にある。
【００２８】
濃度対応ノイズ抑制処理ステップで得られた加重平均フィルタ処理後の濃度成分画像データと、当該加重平均フィルタ処理前の濃度成分画像データの差分は、除去したノイズ量を示す。この場合、当該ノイズ量がノイズ変動幅特定ステップで得られた平坦画像領域から算出されたノイズ変動幅よりも大きい場合は、過度にノイズ抑制処理が行われたと考えられるため、ノイズ抑制処理に制限をかける方が望ましい。
【００２９】
上述の構成によれば、補正処理ステップは、濃度対応ノイズ抑制処理ステップで得られた加重平均フィルタ処理後の濃度成分画像データの画素である加重平均画素と、当該加重平均画素に対応する加重平均フィルタ処理前の濃度成分画像データの注目画素との濃度差分値を変数とし、ノイズ量はノイズ抑制強度算出ステップで算出したノイズ抑制強度を仮想分散値とした正規分布関数に従うと仮定して、正規分布の特性である滑らかな特性をもつ重み係数で各加重平均画素と対応する注目画素とを重み付け加算するため、過度にノイズを抑制しないように補正された強度でノイズを抑制することができるようになる。
【００３０】
同第六の特徴構成は、同請求項６に記載した通り、上述の第一から第五の何れかの特徴構成に加えて、前記ノイズ変動幅特定ステップは、前記写真画像データから濃度成分画像データを抽出する濃度成分抽出ステップと、前記濃度成分抽出ステップで抽出された前記濃度成分画像データを所定サイズの画素ブロックに分割して、画素ブロック毎に平均濃度値、最大濃度値、最小濃度値、濃度ばらつきを含むブロック特性値を算出するブロック特性値算出ステップと、濃度範囲を所定間隔で複数の区間に区分し、前記平均濃度値に基づいて各画素ブロックを対応する区分に割り付ける画素ブロック割付ステップと、各区分に割り付けられた画素ブロックから前記濃度ばらつきが最小となる画素ブロックを代表画素ブロックとして抽出する画素ブロック抽出ステップと、を備え、前記画素ブロック抽出ステップで抽出された複数の代表画素ブロックの各平均濃度値を説明変数とし、各最大濃度値及び各最小濃度値の夫々を従属変数とする二本の回帰直線を算出するものである点にある。
【００３１】
ブロック特性値算出ステップでは、濃度成分抽出ステップで抽出された濃度成分画像データが所定サイズの画素ブロックに分割され、画素ブロック毎に平均濃度値、最大濃度値、最小濃度値、濃度ばらつきを含むブロック特性値が算出される。画素ブロックが画像の平坦部位に属するときには、最大濃度値と最小濃度値にそれほどの差が無く、濃度ばらつきも小さくなるが、画素ブロックが画像の輪郭部位に属するときには、最大濃度値と最小濃度値に大きな差が生じ、濃度ばらつきも比較的大きくなる。つまり、濃度ばらつきが小さい画素ブロックでは、粒状ノイズによるばらつきの影響が顕著に表されているといえる。
【００３２】
画素ブロック割付ステップでは、所定間隔で複数の区間に区分された濃度範囲の何れかの区分に、各画素ブロックがその平均濃度値に基づいて割り付けられる。画素ブロック抽出ステップでは、区分された濃度範囲毎に、濃度ばらつきが最小となる画素ブロックを代表画素ブロックとして抽出する。すなわち、代表画素ブロックの濃度ばらつきは、区分された濃度範囲で最も平坦な領域を含む画素ブロックのノイズの変動幅を示していることになる。
【００３３】
さらに、ノイズ変動幅特定ステップでは、画素ブロック抽出ステップで抽出した区分された濃度範囲毎のノイズの変動幅を利用して、回帰計算を実施するため、区分された濃度範囲ではなく、本来の濃度域として連続したノイズ変動幅を算出できるようになる。実際に発明者らが検証したところ、本発明により得られたノイズ変動幅が、粒状ノイズの変動幅と略一致することが判明している。
【００３４】
同第七の特徴構成は、本発明による写真画像処理装置の特徴構成であり、同請求項７に記載した通り、加重平均フィルタ処理により写真画像データに含まれる粒状のノイズを抑制する写真画像処理装置であって、前記写真画像データを構成する濃度成分画像データから複数の平坦画像領域を探索し、各平坦画像領域の平均濃度値を説明変数とし最大濃度値及び最小濃度値の夫々を従属変数とする二本の回帰直線により前記写真画像データに含まれる粒状のノイズの変動幅を求めるノイズ変動幅特定部と、前記ノイズ変動幅特定部で得られた各回帰直線に基づいて、任意の注目画素に対するノイズ変動幅を求め、当該ノイズ変動幅からノイズ抑制強度を算出するノイズ抑制強度算出部と、当該注目画素と周辺画素の濃度差分値と、前記ノイズ抑制強度算出部で得られたノイズ抑制強度に基づいて、濃度成分画像データに対応したサイズの加重平均フィルタのフィルタ係数を動的に算出する濃度対応フィルタ係数生成部と、前記濃度対応フィルタ係数生成部で得られた各注目画素に対応するフィルタ係数に基づいて、前記濃度成分画像データを加重平均フィルタ処理する濃度対応ノイズ抑制処理部と、前記濃度成分画像データを前記濃度成分画像データに対応したサイズの移動平均フィルタでフィルタ処理する移動平均フィルタ処理部と、前記ノイズ変動幅特定部で得られた各回帰直線に基づいて、前記移動平均フィルタ処理部で得られた各注目画素領域内の移動平均画素のノイズ変動幅を求める移動平均画素ノイズ変動幅特定部と、前記移動平均画素ノイズ変動幅特定部で得られた移動平均画素のノイズ変動幅に基づいて、前記濃度対応ノイズ抑制処理部で得られた各加重平均画素を補正するグラデーション補正処理部と、を含む点にある。
【００３５】
同第八の特徴構成は、同請求項８に記載した通り、上述の第七の特徴構成に加えて、前記グラデーション補正処理部は、各注目画素領域の濃度最大値と濃度最小値の濃度差分値を求める濃度差分値算出部と、前記移動平均画素ノイズ変動幅特定部で得られた移動平均画素のノイズ変動幅と前記濃度差分値算出部で得られた濃度差分値の比率に基づいて、各注目領域の平坦度を算出する平坦度算出部と、を含み、前記平坦度算出部で得られた平坦度を重み係数として、前記濃度対応ノイズ抑制処理部で得られた各加重平均画素と対応する移動平均画素とを重み付け加算してグラデーション再生濃度画像データを算出する点にある。
【００３６】
同第九の特徴構成は、同請求項９に記載した通り、上述の第七の特徴構成に加えて、前記グラデーション補正処理部は、各注目画素領域の濃度値の標準偏差を求める濃度標準偏差算出部と、前記移動平均画素ノイズ変動幅特定部で得られた移動平均画素のノイズ変動幅に基づいて、各注目画素領域内の移動平均画素の濃度値の仮想標準偏差を算出し、前記仮想標準偏差と前記濃度標準偏差算出部で得られた濃度値の標準偏差との比率に基づいて、各注目領域の平坦度を算出する平坦度算出部と、を含み、前記平坦度算出部で得られた平坦度を重み係数として、前記濃度対応ノイズ抑制処理部で得られた各加重平均画素と対応する移動平均画素とを重み付け加算してグラデーション再生濃度画像データを算出する点にある。
【００３７】
同第十の特徴構成は、本発明による写真画像処理プログラムの特徴構成であり、同請求項１０に記載した通り、加重平均フィルタ処理により写真画像データに含まれる粒状のノイズを抑制する写真画像処理プログラムであって、前記写真画像データを構成する濃度成分画像データから複数の平坦画像領域を探索し、各平坦画像領域の平均濃度値を説明変数とし最大濃度値及び最小濃度値の夫々を従属変数とする二本の回帰直線により前記写真画像データに含まれる粒状のノイズの変動幅を求めるノイズ変動幅特定ステップと、前記ノイズ変動幅特定ステップで得られた各回帰直線に基づいて、任意の注目画素に対するノイズ変動幅を求め、当該ノイズ変動幅からノイズ抑制強度を算出するノイズ抑制強度算出ステップと、当該注目画素と周辺画素の濃度差分値と、前記ノイズ抑制強度算出ステップで得られたノイズ抑制強度に基づいて、濃度成分画像データに対応したサイズの加重平均フィルタのフィルタ係数を動的に算出する濃度対応フィルタ係数生成ステップと、前記濃度対応フィルタ係数生成ステップで得られた各注目画素に対応するフィルタ係数に基づいて、前記濃度成分画像データを加重平均フィルタ処理する濃度対応ノイズ抑制処理ステップと、前記濃度成分画像データを前記濃度成分画像データに対応したサイズの移動平均フィルタでフィルタ処理する移動平均フィルタ処理ステップと、前記ノイズ変動幅特定ステップで得られた各回帰直線に基づいて、前記移動平均フィルタ処理ステップで得られた各注目画素領域内の移動平均画素のノイズ変動幅を求める移動平均画素ノイズ変動幅特定ステップと、前記移動平均画素ノイズ変動幅特定ステップで得られた移動平均画素のノイズ変動幅に基づいて、前記濃度対応ノイズ抑制処理ステップで得られた各加重平均画素を補正するグラデーション補正処理ステップと、をコンピュータに実行させる点にある。
【００３８】
同第十一の特徴構成は、同請求項１１に記載した通り、上述の第十の特徴構成に加えて、前記グラデーション補正処理ステップは、各注目画素領域の濃度最大値と濃度最小値の濃度差分値を求める濃度差分値算出ステップと、前記移動平均画素ノイズ変動幅特定ステップで得られた移動平均画素のノイズ変動幅と前記濃度差分値算出ステップで得られた濃度差分値の比率に基づいて、各注目領域の平坦度を算出する平坦度算出ステップと、を含み、前記平坦度算出ステップで得られた平坦度を重み係数として、前記濃度対応ノイズ抑制処理ステップで得られた各加重平均画素と対応する移動平均画素とを重み付け加算してグラデーション再生濃度画像データを算出する点にある。
【００３９】
同第十二の特徴構成は、同請求項１２に記載した通り、上述の第十の特徴構成に加えて、前記グラデーション補正処理ステップは、各注目画素領域の濃度値の標準偏差を求める濃度標準偏差算出ステップと、前記移動平均画素ノイズ変動幅特定ステップで得られた移動平均画素のノイズ変動幅に基づいて、各注目画素領域内の移動平均画素の濃度値の仮想標準偏差を算出し、前記仮想標準偏差と前記濃度標準偏差算出ステップで得られた濃度値の標準偏差との比率に基づいて、各注目領域の平坦度を算出する平坦度算出ステップと、を含み、前記平坦度算出ステップで得られた平坦度を重み係数として、前記濃度対応ノイズ抑制処理ステップで得られた各加重平均画素と対応する移動平均画素とを重み付け加算してグラデーション再生濃度画像データを算出する点にある。
【発明の効果】
【００４０】
以上説明した通り、本発明によれば、写真画像データからノイズの変動幅を特定し、効果的にノイズを抑制できる写真画像処理方法、写真画像処理プログラム、及び写真画像処理装置を提供することができるようになった。
【発明を実施するための最良の形態】
【００４１】
以下、本発明による写真画像処理方法、写真画像処理プログラム、及び写真画像処理装置の実施の形態について、第一の実施形態として説明する。
【００４２】
図１に示すように、写真画像処理装置１は、印画紙Ｐに対して出力画像データに基づいた露光処理を行ない、露光された印画紙を現像処理して写真プリントを生成出力する写真プリンタ２と、写真画像に対するプリントオーダ情報を設定入力するとともに、各種の画像補正処理を行ない、原画像から編集処理した出力画像データを写真プリンタ２に出力する操作ステーション３を備えて構成されている。
【００４３】
操作ステーション３は、現像済みの写真フィルムＦから画像を読み込むフィルムスキャナ３１と、デジタルスチルカメラ等で撮影された画像データが格納されたメモリカード等の画像データ記憶メディアＭから画像データを読み取るメディアドライバ３２と、コントローラ３３としての汎用コンピュータ等を備えている。
【００４４】
図１及び図２に示すように、写真プリンタ２は、ロール状の印画紙Ｐを収容した二系統の印画紙マガジン２１と、印画紙マガジン２１から引き出された印画紙Ｐを所定のプリントサイズに切断するシートカッター２２と、切断後の印画紙Ｐの背面にコマ番号等のプリント情報を印字するバックプリント部２３と、プリントデータに基づいて印画紙Ｐを露光する露光部２４と、露光後の印画紙Ｐを現像、漂白、定着するための各処理液が充填された複数の処理槽２５ａ、２５ｂ、２５ｃを備えた現像処理部２５が印画紙Ｐの搬送経路に沿って配置され、現像処理後に乾燥処理された印画紙Ｐが排出される横送りコンベア２６と、横送りコンベア２６に集積された複数枚の印画紙（写真プリント）Ｐがオーダー単位で仕分けられるソータ２７を備えている。
【００４５】
露光部２４には、搬送機構２８によって副走査方向に搬送される印画紙Ｐに対して、搬送方向に直交する主走査方向に前記プリントデータに基づき変調されたＲＧＢ三色のレーザ光線束を出力して露光する露光ヘッド２４ａが収容されている。
【００４６】
搬送経路に沿って配置された露光部２４や現像処理部２５に、所定のプロセス速度で印画紙Ｐを搬送する複数のローラ対でなる搬送機構２８が配置され、露光部２４の前後には印画紙Ｐを複列に搬送可能なチャッカー式搬送機構２８ａが設けられている。
【００４７】
操作ステーション３に設けられたコントローラ３３には、汎用のオペレーティングシステムの管理下で動作し、写真処理装置１の各種の画像処理や入出力制御を実行するための複数のアプリケーションプログラムがインストールされ、オペレータとの操作インターフェースとしてモニタ３４、キーボード３５、マウス３６等が接続されている。当該アプリケーションプログラムに本発明による画像処理プログラムが含まれる。
【００４８】
コントローラ３３は、そのハードウェア及びソフトウェアが協働して写真処理プロセスを実行するブロックで、以下に、各機能ブロックに分けて説明する。
【００４９】
図３に示すように、コントローラ３３は、フィルムスキャナ３１やメディアドライバ３２によって読み取られた原画像としての写真画像データを受け取り、所定の前処理を行なってメモリ４１に転送する画像入力部４０と、モニタ３４の画面にプリントオーダ情報や画像編集情報を表示するとともに、それらに対して必要なデータ入力のための操作用アイコンを表示するグラフィック操作画面を生成し、或いは表示されたグラフィック操作画面に対するキーボード３５やマウス３６からの入力操作に基づいて各種の制御コマンドを生成するグラフィックユーザーインターフェース部４２と、画像入力部４０から転送される写真画像データ及び画像処理部４７による補正処理後の写真画像データやそのときの補正パラメータ、更には設定されたプリントオーダ情報等が所定領域に区画されて格納されるメモリ４１と、プリントオーダ情報を生成するオーダー処理部４３と、メモリ４１に格納された各写真画像データに対してコマ画像毎または所定枚数のコマ画像に濃度補正処理やコントラスト補正処理等を行なう画像処理部４７を備えている。
【００５０】
さらに、グラフィックユーザーインターフェース部４２からの表示コマンドに基づいてメモリ４１に展開された画像データや各種の入出力用グラフィックデータ等をモニタ３４に表示処理するビデオＲＡＭ等を備えた表示制御部４６と、各種の補正処理が終了した最終の補正画像を写真プリンタ２に出力するためのプリントデータを生成するプリントデータ生成部４４と、顧客のオーダーに応じて最終の補正画像をＣＤ−Ｒ等の記憶媒体に書き込むためのファイル形式に変換するフォーマッタ部４５等を備えている。
【００５１】
フィルムスキャナ３１は、フィルムＦに記録された画像を低解像度ではあるものの高速で読み取るプレスキャンモードと、低速ではあるものの高解像度で読み取る本スキャンモードの二モードで作動するように構成され、プレスキャンモードで読み込まれた低解像度の画像に対して各種の補正処理が行なわれ、その際に前記メモリ４１に記憶された補正パラメータに基づいて本スキャンモードで読み込まれた高解像度の画像に対する最終の補正処理が実行されてプリンタ２に出力される。
【００５２】
同様に、メディアドライバ３２から読み込まれた画像ファイルには高解像度の撮影画像とそのサムネイル画像が含まれ、サムネイル画像に対して後述の各種の補正処理が行なわれ、その際にメモリ４１に記憶された補正パラメータに基づいて高解像度の撮影画像に対する最終の補正処理が実行される。尚、画像ファイルにサムネイル画像が含まれないときには、画像入力部４０で高解像度の撮影画像からサムネイル画像が生成されてメモリ４１に転送される。
【００５３】
このように、頻繁に試行錯誤される各種の編集処理が低解像度の画像に対して実行されることによりコントローラ３３の演算負荷が低減されるように構成されている。
【００５４】
画像処理部４７には、メモリ４１に格納された原画像である写真画像データに対して撮影レンズに起因する歪を補正する歪補正部５０と、粒状ノイズを抑制する粒状ノイズ抑制処理部５１と、画像のエッジを強調し、ノイズを抑制する鮮鋭化処理部５２と、自然なカラーを再現できるようにカラーバランスを調整するカラー補正部５３と、写真プリントのサイズに適した画像サイズに変換する拡縮処理部５４等の複数の画像処理ブロックを備えている。
【００５５】
粒状ノイズ抑制処理部５１は、写真フィルムＦから本スキャンモードで読み込まれた画像データに含まれる粒状ノイズを抑制するフィルム画像用の粒状ノイズ抑制処理部５１と、本発明による画像処理装置として機能し、メディアドライバ３２から読み込まれた高解像度の画像データに含まれる粒状ノイズを抑制するデジタルカメラ撮影画像用の粒状ノイズ抑制処理部５１を備えている。
【００５６】
以下、メディアドライバ３２から入力されるデジタルカメラで撮影された写真画像データの処理について詳述するが、写真フィルムＦから本スキャンモードで読み込まれた写真画像データも同様の処理を行なうことが可能である。
【００５７】
メディアドライバ３２から読み込まれた複数のサムネイル画像データ及び高解像度の画像データがメモリ４１に格納されると、モニタ３４の画面にサムネイル画像データに基づく数コマの写真画像とグラフィック操作画面が表示される。尚、通常、メディアに格納されているデジタル画像はＪＰＥＧ方式により圧縮されているため、逆変換して伸長処理した画像データが格納される。
【００５８】
モニタ３４の画面に表示された各コマ画像に対して、グラフィック操作画面に対するオペレータの操作を介してオーダー処理部４３によりプリント条件、つまり、プリント枚数やプリントサイズ等が設定される。
【００５９】
さらに、各コマ画像に対して、オペレータの操作を介して画像処理部４７により歪補正、鮮鋭化処理、カラー補正が実行され、このとき設定された画像の補正処理条件がメモリ４１に格納される。
【００６０】
モニタ３４の画面に表示された各コマ画像に対するプリント条件の設定や補正処理の操作が終了すると、次画面にスクロールして、全コマ画像に対して同様の処理が実行される。
【００６１】
全ての操作処理が終了し、グラフィック操作画面を介してプリント出力操作が行なわれると、メモリ４１に格納されている高解像度の写真画像データに対して、歪補正、粒状ノイズ抑制処理、鮮鋭化処理、カラー補正処理、拡縮処理等の画像処理が順番に実行され、プリントデータ生成部４４によって画像処理後の画像データがプリントデータに変換生成され、写真プリンタ２に出力される。
【００６２】
高解像度の写真画像データに対する歪補正、鮮鋭化処理、カラー補正処の夫々は、サムネイル画像に対して設定された補正処理条件、つまりメモリ４１に格納された補正処理条件に基づいて自動的に補正処理が実行される。
【００６３】
以下、本発明によるデジタルカメラ撮影画像用の粒状ノイズ抑制処理部５１について詳述する。
【００６４】
図４に示すように、粒状ノイズ抑制処理部５１は、写真画像データを構成する濃度成分画像データから複数の平坦画像領域を探索し、各平坦画像領域の平均濃度値を説明変数とし最大濃度値及び最小濃度値の夫々を従属変数とする二本の回帰直線により写真画像データに含まれる粒状のノイズの変動幅を求めるノイズ変動幅特定部７０と、ノイズ変動幅特定部７０で得られた各回帰直線に基づいて、任意の注目画素に対するノイズ変動幅を求め、当該ノイズ変動幅からノイズ抑制強度を算出するノイズ抑制強度算出部７１と、当該注目画素と周辺画素の濃度差分値と、ノイズ抑制強度算出部７１で得られたノイズ抑制強度に基づいて、濃度成分画像データに対応したサイズの加重平均フィルタのフィルタ係数を動的に算出する濃度対応フィルタ係数生成部７２と、濃度対応フィルタ係数生成部７２で得られた各注目画素に対応するフィルタ係数に基づいて、濃度成分画像データを加重平均フィルタ処理する濃度対応ノイズ抑制処理部７３と、を備えている。
【００６５】
さらに、粒状ノイズ抑制処理部５１は、注目画素と周辺画素の濃度差分値と、ノイズ抑制強度算出部７１で得られたノイズ抑制強度に基づいて、色差成分画像データに対応したサイズの加重平均フィルタのフィルタ係数を動的に算出する色差対応フィルタ係数生成部７５と、色差対応フィルタ係数生成部７５で得られた各注目画素に対応するフィルタ係数に基づいて、色差成分画像データを加重平均フィルタ処理する色差対応ノイズ抑制処理部７６を備えている。
【００６６】
さらに、粒状ノイズ抑制処理部５１は、後述する濃度対応ノイズ抑制処理部７３から得られたノイズ抑制された輝度成分画像データと色差対応ノイズ抑制処理部７６から得られたノイズ抑制された色差成分画像データから得られる各画素のＹＣｂＣｒ成分の夫々をＲＧＢ成分に変換するＲＧＢ変換処理部７４を備えている。
【００６７】
ノイズ変動幅特定部７０は、写真画像データから濃度成分画像データを抽出する濃度成分抽出部７０１と、濃度成分抽出部７０１で抽出された濃度成分画像データを所定サイズの画素ブロックに分割して、画素ブロック毎に平均濃度値、最大濃度値、最小濃度値、濃度ばらつきを含むブロック特性値を算出するブロック特性値算出部７０２と、濃度範囲を所定間隔で複数の区間に区分し、前記平均濃度値に基づいて各画素ブロックを対応する区分に割り付ける画素ブロック割付部７０３と、各区分に割り付けられた画素ブロックから前記濃度ばらつきが最小となる画素ブロックを代表画素ブロックとして抽出する画素ブロック抽出部７０４と、画素ブロック抽出部７０４で抽出された複数の代表画素ブロックの各平均濃度値を説明変数とし、各最大濃度値及び各最小濃度値の夫々を従属変数とする二本の回帰直線を算出する回帰直線算出部７０５を備えている。
【００６８】
詳述すると、濃度成分抽出部７０１は、以下の変換式に基づいて、メモリ４１から読み出した高解像度の原画像である写真画像データ（例えば、８００万画素であれば、横３２６４画素×縦２４４８画素程度である）を構成する各画素のＲ（赤）、Ｇ（緑）、Ｂ（青）の夫々の色成分（以下、「ＲＧＢ」と記す。）を、輝度色差成分に変換してメモリ４１に格納する。
Ｙ    ＝    ０．３３３３３Ｒ＋０．３３３３４Ｇ＋０．３３３３３Ｂ
Ｃｂ  ＝  −０．１６８７４Ｒ−０．３３１２６Ｇ＋０．５００００Ｂ
Ｃｒ  ＝    ０．５００００Ｒ−０．４１８６９Ｇ−０．０８１３１Ｂ
【００６９】
ブロック特性値算出部７０２は、図５（ａ）に示すように、濃度成分抽出部７０１で抽出された輝度成分画像データを濃度成分画像データとして、所定サイズの画素ブロックに分割し、画素ブロック毎に平均濃度値、最大濃度値、最小濃度値を算出する。続いて、図５（ｂ）に示すように、平均濃度値を閾値として、閾値より高濃度の画素Ｈの平均濃度値である高濃度側平均濃度値と、当該閾値より低濃度の画素Ｌの平均濃度値である低濃度側平均濃度値を算出し、さらに、画素ブロックの濃度ばらつきを示すものとして、高濃度側平均濃度値と低濃度側平均濃度値との差分を算出する。
【００７０】
このように、ブロック特性値算出部７０２は、所定サイズの画素ブロック毎に、平均濃度値、最大濃度値、最小濃度値、濃度ばらつきを含むブロック特性値を算出し、処理対象画像データの一行目の左端の画素ブロックから行方向に配列された各画素ブロックを、順次注目画素ブロックとしてブロック特性値を算出する処理を繰り返し、右端の画素ブロックまで処理が終了すると二行目の左端の画素ブロックから行方向に同様の処理を行ない、最終行の右端の画素ブロックまで処理を繰り返す。
【００７１】
画素ブロック割付部７０３は、濃度範囲を所定間隔で複数の区間に区分し、ブロック特性値算出部７０２で算出した平均濃度値に基づいて各画素ブロックを対応する区分に割り付ける。例えば、図５（ｃ）は、画素ブロック割付部７０３により２５６階調の濃度範囲を１６個の区分に区分し、横軸が当該区分した濃度範囲を示し、縦軸が濃度ばらつきを示すよう、各画素ブロックを割り付けた結果を示す。
【００７２】
画素ブロック抽出部７０４は、画素ブロック割付部７０３により各区分に割り付けられた画素ブロックから、濃度ばらつきが最小となる画素ブロックを代表画素ブロックとして抽出する。例えば、図５（ｃ）に示すように、画像ブロックが高濃度の濃度範囲と低濃度の濃度範囲の２つに集中して割り付けられる画像データの場合であっても、高濃度の濃度範囲と低濃度の濃度範囲の夫々において、濃度ばらつきが最小の画素ブロックを抽出する。
【００７３】
つまり、所定間隔で区分された濃度範囲の各区分内で、濃度ばらつきが少ない、最も平坦な領域の画素ブロックが、ノイズ抑制処理の代表画素ブロックとして抽出される。
【００７４】
回帰直線算出部７０５は、例えば、図５（ｄ）に示すように、画素ブロック抽出部７０４で抽出された複数の代表画素ブロックの各平均濃度値を説明変数（ｘ軸）とし、各最大濃度値及び各最小濃度値の夫々を従属変数とする二本の回帰直線（傾きａ、切片ｂ、相関係数ｒ）を最小自乗法等で算出し、回帰直線で区画される範囲を粒状のノイズの変動幅として特定する。
【００７５】
本願発明者は、上述の二本の回帰直線を算出した場合に、算出対象の画像の内容に関係なく、大抵の画像において相関係数ｒは０．９９以内に収まることを知見している。また、本願発明者は、例えば、絵柄が複雑な画像、アンダーコントラストが少ない画像、十分にある濃度範囲の画素ブロックが収集できない画像等に対して上述の二本の回帰直線を算出した場合は、十分な相関係数ｒが得られず、算出された二本の回帰直線が示す平坦領域の画素ブロックの平均濃度値と最大濃度値及び最小濃度値の関係は、線形性の関係が十分に高いとはいえないことを知見している。そこで、相関係数ｒの値に応じて、それぞれ以下のように補正を行う。
【００７６】
相関係数ｒが０．９９以上の場合は、十分な相関係数ｒが得られているとして回帰直線の補正は行わない。相関係数ｒが第一の閾値（例えば、０．９７）より大きく、０．９９未満の場合は、回帰直線の傾きａを補正する。例えば、図６（ａ）に示すように、代表画素ブロックの最大濃度値の平均値及び代表画素ブロックの最小濃度値の平均値を示す座標を中心として、回帰直線を回転させて傾きａを１に近づけるように補正し、それに伴なって切片ｂも補正する。
【００７７】
詳述すると、図６（ａ）に示すように、当該傾きａの１への近づき度合は、相関係数ｒが第一の閾値の場合に傾きａが１となり、相関係数ｒが０．９９の場合に傾きａを補正しないように、第一の閾値と０．９９の差分値（例えば、０．０２）に対する相関係数ｒの値に応じて決定し、即ち、以下の式で決定する。
補正後ａ  ＝  補正前ａ−（補正前ａ−１）×（ｒ−第一の閾値）／（０．９９−第一の閾値）
【００７８】
当該平均値を示す座標を（ＡｖｅＸ，ＡｖｅＹ）とするとき、補正後の切片ｂは以下の式で求められる。
補正後ｂ  ＝
補正前ｂ×ｒｔｏ＋（ＡｖｅＹ−補正後のａ×ＡｖｅＸ）×（１−ｒｔｏ）
ここに、ｒｔｏ  ＝  （ｒ−第一の閾値）／（第一の閾値−第二の閾値）
【００７９】
相関係数が第二の閾値（例えば、０．９５）以上であり、第一の閾値未満の場合は、回帰直線を回転させて傾きａを１に設定した後、さらに、回帰直線回転後の切片ｂを０に近づけるように設定する。詳述すると、図６（ｂ）に示すように、当該切片ｂの０への近づき度合は、相関係数ｒが第二の閾値の場合に切片ｂが０となり、相関係数ｒが第一の閾値の場合に、切片ｂが移動しないように、第一の閾値と第二の閾値の差分値（例えば、０．０２）に対する相関係数ｒの値に応じて決定し、即ち、以下の式で決定する。
補正後ｂ  ＝  補正前ｂ×（ｒ−第二の閾値）／（第一の閾値−第二の閾値）
【００８０】
相関係数ｒが第二の閾値未満の場合は、算出された回帰直線は信頼できないものとして、傾きａが１、切片ｂが０の直線に変更する。
【００８１】
このようにして、原画像データを構成する所定間隔で区分された濃度範囲の各区分内で、濃度ばらつきが少ない、最も平坦な領域の画素ブロックから、各画素の濃度に対する粒状ノイズの変動幅が把握されるようになる。
【００８２】
ノイズ抑制強度算出部７１は、ノイズ変動幅特定部７０で得られた平坦画像領域から算出された二本の回帰直線（以下、代表画素ブロックの最大濃度値から算出された回帰直線を示す式を「上限式」、代表画素ブロックの最小濃度値から算出された回帰直線を示す式を「下限式」と記す。）に基づいて、濃度成分画像データの任意の注目画素に対するノイズ変動幅を算出する。即ち、上述の二本の回帰直線のｘ軸の値は、二本の回帰直線の説明変数である代表画素ブロックの平均濃度値であるが、当該ｘ軸の値を注目画素の濃度値であるとみなして、回帰直線が示す平坦画像領域のノイズの変動幅が、当該注目画素の濃度値に対するノイズの変動幅を示すものとして算出する。
【００８３】
尚、以下では、特筆しない限り、上限式の傾きをＹｍａｘＡ、切片をＹｍａｘＢ、及び、下限式の傾きＹｍｉｎＡ、切片をＹｍｉｎＢと示す。注目画素の座標（ｘ,ｙ）が示す濃度値をｄａｔｘ，ｙとすると、注目画素の濃度値に対するノイズの変動幅の上限値ＭａｘＹｘ,ｙと下限値ＭｉｎＹｘ,ｙは、以下の式で示すことができる。
【数１】
2010108008.tif 000003
【００８４】
したがって、注目画素の濃度値に対するノイズの変動幅Ｗｘ，ｙは、上限値ＭａｘＹｘ,ｙと下限値ＭｉｎＹｘ,ｙの差分値として、以下の式で示すことができる。
【数２】
2010108008.tif 000004
【００８５】
ここで、当該算出したノイズの変動幅Ｗｘ，ｙは、ノイズの大きさを示している為、ノイズ抑制の強度設定指標として利用することができる。
【００８６】
例えば、図７（ａ）に示すように、所定サイズの背景画像等の平坦画像領域における濃度値の分布と、図７（ｂ）に示すように、図７（ａ）に示す平坦画像領域と被写体の輪郭部を含む輪郭画像領域とを含む所定サイズの画像領域における濃度値の分布とを、正規化して対比すると、平坦画像領域と輪郭画像領域を含む画像領域では、平坦画像領域及び輪郭画像領域の夫々を正規化した場合の標準偏差σの値よりも、画像領域全般における標準偏差σの方が大きくなることがわかる。即ち、注目する画像領域の絵柄によっては、注目する画像領域に含まれる全画素の濃度値から濃度値のばらつき度合を示す標準偏差σを算出しても信頼性が高いとは言えない。
【００８７】
そこで、濃度値のばらつき度合を示す標準偏差σは、最大濃度値と最小濃度値の差に関連があることに着目し、さらに、当該算出したノイズ変動幅Ｗｘ，ｙ内の濃度値の分布が正規分布に従うものとみなすことにより、注目する画像領域の絵柄に依存しない標準偏差σ ｘ,ｙを算出することができる。また、当該標準偏差標準偏差σ ｘ,ｙは、濃度値のばらつき度合を示しているため、ノイズ抑制する場合のノイズ抑制強度として利用する。
【００８８】
具体的には、ノイズ抑制強度算出部７１は、ノイズの変動幅Ｗｘ，ｙをσ化係数ＳｉｇｍａＣｏｅｆで割った値を標準偏差σ ｘ,ｙとして算出する。尚、当該算出したノイズ抑制強度としての標準偏差σ ｘ,ｙに基づいて、過度にノイズ抑制することを避けるため、最大許容値Ｌｅｖｅｌと最小許容値１を設け、以下の式で示すように、標準偏差σ ｘ,ｙを算出する。
【数３】
2010108008.tif 000005
【００８９】
σ化係数ＳｉｇｍａＣｏｅｆは、例えば、ノイズ変動幅の半値の約７割に当たる値が標準偏差σを示すものと仮定した場合は、２．８５（＝１／（０．５×０．７））となる。また、本願発明者は、当該標準偏差σの値は、１０程度の値であり、最大でも１６程度に収まることを知見しており、最大許容値は、１６程度に設定することが望ましい。ただし、σ化係数ＳｉｇｍａＣｏｅｆ、最大許容値、最小許容値の値は、上述の数値に限るものではなく、適宜設計事項として調整して構わない。
【００９０】
したがって、ノイズ抑制強度算出部７１では、ノイズ変動幅特定部７０で得られた、平坦画像領域から算出されたノイズの変動幅に基づいて、任意の注目画素（ｘ,ｙ）に対するノイズの変動幅Ｗｘ，ｙが算出され、当該ノイズ変動幅Ｗｘ，ｙからノイズ抑制強度σ ｘ,ｙが算出されるため、注目画素（ｘ,ｙ）とその周辺画素の濃度差が大きな輪郭画像領域や絵柄が細かい画像領域等において、画像の内容を適切に表現する上で必要となる画素が、ノイズとして認識されないよう制限されたノイズの抑制強度σ ｘ,ｙを算出できるようになる。
【００９１】
濃度対応フィルタ係数生成部７２は、後述する濃度対応ノイズ抑制処理部７３において、濃度成分画像データを加重平均フィルタ処理する場合に利用する加重平均フィルタのフィルタ係数を動的に算出する。
【００９２】
詳述すると、以下の式で示すように、注目画素（ｘ,ｙ）を中心とした所定サイズのフィルタ領域内の各周辺画素（ｘ＋ｉ,ｙ＋ｊ）へのフィルタ係数Ｗｇｔｉ,ｊを、注目画素の濃度値ｄａｔ ｘ,ｙと各周辺画素の濃度値ｄａｔ ｘ＋ｉ,ｙ＋ｊの差分値ｚ ｉ,ｊと、ノイズ抑制強度算出部７１で得られたノイズ抑制強度、即ち、注目画素（ｘ,ｙ）を中心とした注目領域内の周辺画素の濃度値のばらつきを示す標準偏差σ ｘ,ｙを標準偏差とする正規分布関数に基づいて算出する。
【数４】
2010108008.tif 000006
【００９３】
尚、ｉ,ｊは、注目画素の座標（ｘ,ｙ）を（０,０）とみなした場合の周辺画素の座標を示す。
【００９４】
フィルタ係数Ｗｇｔ  ｉ,ｊは、標準偏差σｘ,ｙが大きい場合、即ち、フィルタ領域内における周辺画素の濃度値のばらつきが大きい場合は、注目画素に対して周辺画素の影響度を大きくする必要があるため、フィルタ係数Ｗｇｔ  ｉ,ｊの値は、大きくなる。
【００９５】
濃度フィルタ係数生成部７２は、処理対象となる濃度成分画像データの一行目の左端の画素から行方向に配列された各画素を、順次注目画素として、所定サイズ（例えば、５×５程度）のフィルタ領域内の各周辺画素へのフィルタ係数の算出を繰り返し、右端の画素まで当該算出処理が終了すると二行目の左端の画素から行方向に同様の処理を行ない、最終行の右端の画素まで処理を繰り返す。
【００９６】
したがって、濃度対応フィルタ係数生成部７２では、ノイズ抑制強度算出部７１で得られたノイズ抑制強度σ ｘ,ｙを標準偏差とする正規分布関数に基づいてフィルタ係数Ｗｇｔ  ｉ,ｊが算出されるため、滑らかな特性を持つフィルタ係数Ｗｇｔ  ｉ,ｊを算出することができる。
【００９７】
濃度対応ノイズ抑制処理部７３は、濃度対応フィルタ係数生成部７２で得られた各注目画素に対応するフィルタ係数で構成される加重平均フィルタにより、処理対象となる濃度成分画像データの一行目の左端の画素から行方向に配列された各画素を、順次注目画素として加重平均フィルタ処理を繰り返し、右端の画素まで処理が終了すると二行目の左端の画素から行方向に同様の処理を行ない、最終行の右端の画素まで処理を繰り返す。加重平均フィルタ処理された濃度成分画像データの濃度値、即ち、濃度成分画像データの加重平均画素の濃度値Ｙａｖｅ ｘ,ｙは、以下の式で示すことができる。
【数５】
2010108008.tif 000007
【００９８】
ここで、ｎＹは、注目画素と周辺画素の最大距離を示し、例えば、フィルタサイズが５×５の場合には、ｎＹ＝５÷２≒２となる。
【００９９】
したがって、濃度対応ノイズ抑制処理部７３では、ノイズ抑制強度算出部７１及び濃度対応フィルタ係数生成部７２を介して得られた加重平均フィルタを利用して、加重平均フィルタ処理を実行するため、ノイズの変動幅Ｗｘ，ｙを超えた加重平均フィルタ処理を実行することにより、元の画像に含まれていた絵柄や輪郭までノイズとして抑制して消してしまう不具合を軽減できる。
【０１００】
色差フィルタ係数生成部７５は、後述する色差対応ノイズ抑制処理部７６において、濃度成分抽出部７０１で抽出された色差成分画像データを加重平均フィルタ処理する場合に利用する、加重平均フィルタのフィルタ係数を動的に算出する。
【０１０１】
詳述すると、色差フィルタ係数生成部７５は、濃度対応フィルタ係数生成部７２と同様に、数４で示した式に従い、注目画素（ｘ,ｙ）を中心とした所定サイズのフィルタ領域内の各周辺画素（ｘ＋ｉ,ｙ＋ｊ）へのフィルタ係数Ｗｇｔｉ,ｊを、注目画素の濃度値ｄａｔ ｘ,ｙと各周辺画素の濃度値ｄａｔ ｘ＋ｉ,ｙ＋ｊの差分値ｚ ｉ,ｊと、ノイズ抑制強度算出部７１で得られたノイズ抑制強度、即ち、注目画素（ｘ,ｙ）を中心とした注目領域内の周辺画素の濃度値のばらつきを示す標準偏差σ ｘ,ｙに基づいて算出する。
【０１０２】
色差フィルタ係数生成部７５は、処理対象となる色差成分画像データの一行目の左端の画素から行方向に配列された各画素を、順次注目画素として、濃度対応フィルタ係数生成部７２とは異なる所定サイズ（例えば、９×９程度）のフィルタ領域内の各周辺画素へのフィルタ係数の算出を繰り返し、右端の画素まで当該算出処理が終了すると二行目の左端の画素から行方向に同様の処理を行ない、最終行の右端の画素まで処理を繰り返す。
【０１０３】
ここで、色差フィルタ係数生成部７５は、濃度対応フィルタ係数生成部７２とは異なる所定サイズ（例えば、９×９程度）のフィルタ領域内の各周辺画素へのフィルタ係数の算出を行っているが、その理由は、色ノイズは比較的大きな画素サイズで現れるという特性に対応するためである。
【０１０４】
詳述すると、通常、ＪＰＥＧ方式では、ＹＣｒＣｂ変換された画素を８×８のブロック毎にＤＣＴ変換し、得られた周波数成分を量子化し、ジグザグスキャンにより配列された量子化データをハフマン符合化するが、各ブロック内の色差成分画像データはＤＣＴ変換する前に輝度成分画像データの１／４に間引き処理される。人間の目が細かな輝度の違いには敏感であるのに対して、色相と彩度の細かな違いには比較的気付かないという特性に基づいて、色差成分画像データを削減するためである。
【０１０５】
従って、デジタルカメラで撮影された写真画像に現れるざらつきの原因である粒状のノイズを輝度ノイズと色ノイズに分けると、輝度ノイズが略一画素単位で現れるのに対して、色ノイズは比較的大きな画素サイズで現れ、その傾向がＤＣＴ変換時のブロックサイズと相関を持つという特性が見られる。
【０１０６】
そのため、色ノイズを抑制するためには、フィルタサイズをＪＰＥＧによる圧縮のためのブロックサイズより大きなフィルタサイズに設定することが好ましく、本実施形態では、９×９或は１１×１１程度のフィルタサイズに設定する。
【０１０７】
また、色差フィルタ係数生成部７５は、濃度対応フィルタ係数生成部７２と同様に、濃度成分の変化に基づいて決定されたノイズ抑制強度を算出しているが、その理由は、人間の視覚の特性上、色相成分の変化よりも敏感である濃度成分の変化に基づいて決定されたノイズ抑制強度に基づいて、後述する色差成分画像データのノイズ抑制処理を行い、加重平均フィルタ処理後に色変化が小さく、周囲と同化した画像が生成されることを回避するためである。
【０１０８】
色差対応ノイズ抑制処理部７６は、色差フィルタ係数生成部７５で得られた各注目画素に対応するフィルタ係数で構成される加重平均フィルタにより、処理対象となる色差成分画像データの一行目の左端の画素から行方向に配列された各画素を、順次注目画素として加重平均フィルタ処理を繰り返し、右端の画素まで処理が終了すると二行目の左端の画素から行方向に同様の処理を行ない、最終行の右端の画素まで処理を繰り返す。加重平均フィルタ処理された色差成分画像データの加重平均画素値Ｃｂａｖｅ ｘ,ｙ, Ｃｒａｖｅ ｘ,ｙは、以下の式で示すことができる。
【数６】
2010108008.tif 000008
【０１０９】
尚、Ｃｂｄａｔ ｘ,ｙ, Ｃｒｄａｔａｖｅ ｘ,ｙは、加重平均フィルタ処理される前の色差成分画像データの画素値を示し、ｎＣは、注目画素と周辺画素の最大距離を示す。例えば、フィルタサイズが９×９の場合には、ｎＣ＝９÷２≒４となる。
【０１１０】
ＲＧＢ変換処理部７４は、濃度対応ノイズ抑制処理部７３で得られた加重平均フィルタ処理後の濃度成分画像データと、色差対応ノイズ抑制処理部７６で得られた加重平均フィルタ処理後の色差成分画像データで構成される、ノイズ抑制後の新たな輝度色差成分を、次式に基づいてＲＧＢ色成分に変換して出力する。
Ｒ  ＝  Ｙ＋１．２９６０４８Ｃｂ−０．２２９２７６Ｃｒ
Ｇ  ＝  Ｙ−０．８２００８７Ｃｂ−０．９４３４３０Ｃｒ
Ｂ  ＝  Ｙ−０．４７５９６１Ｃｂ＋１．１７２７０６Ｃｒ
【０１１１】
さらに、図４に示すように、濃度対応ノイズ抑制処理部７３は、ノイズ抑制強度算出部７１で得られたノイズ抑制強度を仮想分散値とし、濃度対応ノイズ抑制処理部７３で得られた各加重平均画素と対応する注目画素との濃度差分値を変数とする正規分布関数に基づいて重み係数を算出し、当該重み係数で各加重平均画素と対応する注目画素とを重み付け加算して補正濃度画像データを算出する補正処理部７３１を備えていても構わない。
【０１１２】
詳述すると、濃度対応ノイズ抑制処理部７３で得られた加重平均フィルタ処理後の濃度成分画像データと、当該加重平均フィルタ処理前の濃度成分画像データの差分は、除去したノイズ量を示す。この場合、当該ノイズ量がノイズ変動幅特定部７０で得られた平坦画像領域から算出されたノイズ変動幅よりも大きい場合は、過度にノイズ抑制処理が行われたと考えられるため、ノイズ抑制処理に制限をかける方が望ましい。
【０１１３】
そこで、補正処理部７３１は、以下の式に示すように、濃度対応ノイズ抑制処理部７３で得られた加重平均フィルタ処理後の濃度成分画像データの画素である加重平均画素と、当該加重平均画素に対応する加重平均フィルタ処理前の濃度成分画像データの注目画素との濃度差分値を変数ｔｍｐとし、ノイズ量はノイズ抑制強度算出部７１で算出したノイズ抑制強度が示す標準偏差σ ｘ,ｙの正規分布に従うと仮定して、正規分布を利用することにより滑らかな特性をもつ重み係数ｗｇｔを算出する。
【数７】
2010108008.tif 000009
【０１１４】
さらに、補正処理部７３１は、以下の式に示すように、当該算出した重み係数ｗｇｔで、各加重平均画素と対応する注目画素とを重み付け加算することにより、過度にノイズ抑制処理を行わないように抑制された補正濃度画像データの画素の濃度値Ｙａｖｅ´ ｘ,ｙを算出し、当該補正濃度画像データをＲＧＢ変換処理部７４に出力する。
【数８】
2010108008.tif 000010
【０１１５】
したがって、補正処理部７３１では、濃度対応ノイズ抑制処理部７３で得られた加重平均フィルタ処理後の濃度成分画像データの画素である加重平均画素と、当該加重平均画素に対応する加重平均フィルタ処理前の濃度成分画像データの注目画素との濃度差分値を変数とし、ノイズ量はノイズ抑制強度算出部７１で算出したノイズ抑制強度を仮想分散値とした正規分布関数に従うと仮定して、正規分布の特性である滑らかな特性をもつ重み係数で各加重平均画素と対応する注目画素とを重み付け加算するため、過度にノイズを抑制しないように補正された強度でノイズを抑制することができる。
【０１１６】
この場合、ＲＧＢ変換処理部７４は、補正処理部７３１から補正濃度画像データが入力されると、濃度対応ノイズ抑制処理部７３で得られた加重平均フィルタ処理後の濃度成分画像データと同様に取り扱い、当該濃度成分画像データと、色差対応ノイズ抑制処理部７６で得られた加重平均フィルタ処理後の色差成分画像データで構成される、ノイズ抑制後の新たな輝度色差成分を、上述の通り、ＲＧＢ色成分に変換して出力する。
【０１１７】
以下に、上述した粒状ノイズ抑制処理部５１による各処理の手順を、図８に示すフローチャートに基づいて説明する。
【０１１８】
まず、写真画像データを構成する各画素のＲＧＢ成分データが、ノイズ変動幅特定部７０に入力されると、ノイズ変動幅特定ステップ（Ｓ１）が実行される。
【０１１９】
ノイズ変動幅特定ステップ（Ｓ１）は、濃度成分抽出ステップ（Ｓ１１）と、ブロック特性値算出ステップ（Ｓ１２）と、画素ブロック割付ステップ（Ｓ１３）と、画素ブロック抽出ステップ（Ｓ１４）と、回帰直線算出ステップ（Ｓ１５）で構成される。
【０１２０】
まず、濃度成分抽出ステップが実行され、入力された各画素のＲＧＢ成分の夫々がＹＣｂＣｒ成分に変換され、輝度成分画像データと色差成分画像データが出力される。出力された輝度成分画像データは、濃度成分画像データとして扱われる（Ｓ１１）。
【０１２１】
出力された濃度成分画像データがブロック特性値算出部７０２に入力されると、ブロック特性値算出ステップが実行され、所定サイズの画素ブロックに分割された濃度成分画像データの画素ブロック毎に、平均濃度値、最大濃度値、最小濃度値、濃度ばらつきを含むブロック特性値が算出され、当該ブロック特性値と画素ブロックに分割された濃度成分画像データが画素ブロック割付部７０３へ出力される（Ｓ１２）。
【０１２２】
濃度成分画像データと画素ブロック毎のブロック特性値が画素ブロック割付部７０３に入力されると、画素ブロック割付ステップが実行され、濃度成分画像データの濃度範囲が所定間隔で複数の区間に区分され、入力された画素ブロック特性値の平均濃度値に基づいて、濃度成分画像データの各画素ブロックが対応する区分に割り付けられる（Ｓ１３）。
【０１２３】
続いて、画素ブロック抽出部７０４の画素ブロック抽出ステップが実行され、画素ブロック割付ステップ（Ｓ１３）で各区分に割り付けられた画素ブロックから、ブロック特性値の濃度ばらつきが最小となる最も平坦な領域の画素ブロックが、代表画素ブロックとして抽出される（Ｓ１４）。
【０１２４】
続いて、ノイズ変動幅特定部７０の回帰直線算出ステップが実行され、画素ブロック割付ステップ（Ｓ１４）で抽出された複数の平坦画像領域を示す代表画素ブロックの各平均濃度値を説明変数とし、各最大濃度値及び各最小濃度値の夫々を従属変数とする二本の回帰直線が算出され、回帰直線で区画される範囲が粒状のノイズの変動幅として特定される。さらに、当該算出された二本の回帰直線の相関係数に応じて、二本の回帰直線の傾きａ及び切片ｂの補正が行われる（Ｓ１５）。
【０１２５】
続いて、ノイズ抑制強度算出部７１のノイズ抑制強度算出ステップが実行され、ノイズ変動幅特定ステップ（Ｓ１）で平坦画像領域から算出された各回帰直線に基づいて、濃度成分画像データの任意の注目画素（ｘ，ｙ）に対するノイズ変動幅Ｗｘ，ｙが算出され、当該ノイズ変動幅Ｗｘ，ｙからノイズ抑制強度σ ｘ,ｙが算出される（Ｓ２）。
【０１２６】
続いて、濃度対応フィルタ係数生成部７２の濃度対応フィルタ係数生成ステップが実行され、注目画素と周辺画素の濃度差分値ｚ ｉ,ｊと、ノイズ抑制強度算出ステップ（Ｓ２）で得られたノイズ抑制強度σ ｘ,ｙに基づいて、濃度成分画像データに対応したサイズの加重平均フィルタのフィルタ係数Ｗｇｔｉ,ｊが動的に算出される（Ｓ３）。
【０１２７】
続いて、濃度対応ノイズ抑制処理部７３の濃度対応ノイズ抑制処理ステップが実行され、濃度対応フィルタ係数生成ステップ（Ｓ３）で得られた各注目画素に対応するフィルタ係数Ｗｇｔｉ,ｊで構成される加重平均フィルタにより、処理対象となる濃度成分画像データの一行目の左端の画素から行方向に配列された各画素を、順次注目画素として加重平均フィルタ処理が繰り返され、右端の画素まで処理が終了すると二行目の左端の画素から行方向に同様の処理が行なわれ、最終行の右端の画素まで処理が繰り返され、加重平均フィルタ処理された濃度成分画像データの濃度値、即ち、輝度成分画像データの加重平均画素の濃度値Ｙａｖｅ ｘ,ｙが算出される（Ｓ４）。
【０１２８】
また、濃度成分抽出ステップ（Ｓ１１）で抽出された濃度成分画像データと色差成分画像データと、ノイズ抑制強度算出ステップ（Ｓ２）で算出されたノイズ抑制強度σ ｘ,ｙが色差対応フィルタ係数生成部７５に入力されると、色差対応フィルタ係数生成ステップが実行され、注目画素と周辺画素の濃度差分値ｚ ｉ,ｊと、ノイズ抑制強度算出ステップ（Ｓ２）で得られたノイズ抑制強度σ ｘ,ｙに基づいて、色差成分画像データに対応したサイズの加重平均フィルタのフィルタ係数Ｗｇｔｉ,ｊが動的に算出される（Ｓ５）。
【０１２９】
続いて、色差対応ノイズ抑制処理部７６の色差対応ノイズ抑制処理ステップが実行され、色差対応フィルタ係数生成ステップ（Ｓ５）で得られた各注目画素に対応するフィルタ係数Ｗｇｔｉ,ｊで構成される加重平均フィルタにより、処理対象となる色差成分画像データの一行目の左端の画素から行方向に配列された各画素を、順次注目画素として加重平均フィルタ処理が繰り返され、右端の画素まで処理が終了すると二行目の左端の画素から行方向に同様の処理が行なわれ、最終行の右端の画素まで処理が繰り返され、加重平均フィルタ処理された色差成分画像データの加重平均画素値Ｃｂａｖｅ ｘ,ｙ, Ｃｒａｖｅ ｘ,ｙが算出される（Ｓ６）。
【０１３０】
また、濃度対応ノイズ抑制処理部７３に補正処理部７３１が備えられている場合、濃度対応ノイズ抑制処理ステップ（Ｓ４）は、補正処理ステップ（Ｓ４１）をさらに備えて構成されている。
【０１３１】
詳述すると、濃度対応ノイズ抑制処理ステップ（Ｓ４）で加重平均フィルタ処理された濃度成分画像データの濃度値、即ち、輝度成分画像データの加重平均画素の濃度値Ｙａｖｅ ｘ,ｙと、濃度成分抽出ステップ（Ｓ１１）で抽出された濃度成分画像データが補正処理部７３１に入力されると補正処理ステップが実行され、ノイズ抑制強度算出ステップ（Ｓ２）で得られたノイズ抑制強度σ ｘ,ｙを標準偏差とし、濃度対応ノイズ抑制処理ステップ（Ｓ４）で得られた各加重平均画素と対応する注目画素との濃度差分値を変数ｔｍｐとする正規分布関数に基づいて重み係数ｗｇｔが算出される。
【０１３２】
また、補正処理ステップでは、当該算出された重み係数ｗｇｔで、各加重平均画素と対応する注目画素とが重み付け加算され、過度にノイズ抑制処理が行われないように抑制された補正濃度画像データＹａｖｅ´ ｘ,ｙが算出され、当該補正濃度画像データＹａｖｅ´ ｘ,ｙがＲＧＢ変換処理部７４に出力される（Ｓ４１）。
【０１３３】
続いて、濃度対応ノイズ抑制処理ステップ（Ｓ４）で加重平均フィルタ処理された濃度成分画像データ、あるいは、補正処理ステップ（Ｓ４１）で算出された補正濃度画像データと、色差対応ノイズ抑制処理ステップ（Ｓ６）で得られた加重平均フィルタ処理後の色差成分画像データがＲＧＢ変換処理部７４に入力されると、ＲＧＢ変換処理ステップが実行され、当該入力された、濃度成分画像データ、あるいは、補正濃度画像データと、色差成分画像データと、で構成されるノイズ抑制後の新たな輝度色差成分が、ＲＧＢ色成分に変換されて出力され、粒状ノイズ抑制処理部５１による粒状ノイズの抑制処理が終了される（Ｓ８）。
【０１３４】
尚、上述の各ステップにおいて、各ステップ間の入出力データをメモリ４１に保存するよう構成し、各ステップは、各ステップ間で直接データの入出力を行わずに、メモリ４１に保存されたデータにアクセスするよう構成しても構わない。
【０１３５】
上述した粒状ノイズ抑制処理部５１による各処理は、コントローラ３３に備えたハードディスクにインストールされた本発明の写真画像処理プログラムが実行されることにより実現される。
【０１３６】
つまり、加重平均フィルタ処理により写真画像データに含まれる粒状のノイズを抑制する写真画像処理プログラムであって、写真画像データを構成する濃度成分画像データから複数の平坦画像領域を探索し、各平坦画像領域の平均濃度値を説明変数とし最大濃度値及び最小濃度値の夫々を従属変数とする二本の回帰直線により前記写真画像データに含まれる粒状のノイズの変動幅を求めるノイズ変動幅特定ステップ（Ｓ１）と、ノイズ変動幅特定ステップ（Ｓ１）で得られた各回帰直線に基づいて、任意の注目画素に対するノイズ変動幅Ｗｘ，ｙを求め、当該ノイズ変動幅Ｗｘ，ｙからノイズ抑制強度σ ｘ,ｙを算出するノイズ抑制強度算出ステップ（Ｓ２）と、当該注目画素と周辺画素の濃度差分値ｚ ｉ,ｊと、ノイズ抑制強度算出ステップ（Ｓ２）で得られたノイズ抑制強度σ ｘ,ｙに基づいて、濃度成分画像データに対応したサイズの加重平均フィルタのフィルタ係数Ｗｇｔｉ,ｊを動的に算出する濃度対応フィルタ係数生成ステップ（Ｓ３）と、濃度対応フィルタ係数生成ステップ（Ｓ３）で得られた各注目画素に対応するフィルタ係数Ｗｇｔｉ,ｊに基づいて、濃度成分画像データを加重平均フィルタ処理する濃度対応ノイズ抑制処理ステップ（Ｓ４）と、をコンピュータに実行させる写真画像処理プログラムが記憶されたＣＤやＤＶＤ等の記憶媒体を介してインストールされている。
【０１３７】
さらに、当該注目画素と周辺画素の濃度差分値ｚ ｉ,ｊと、ノイズ抑制強度算出ステップ（Ｓ２）で得られたノイズ抑制強度σ ｘ,ｙに基づいて、色差成分画像データに対応したサイズの加重平均フィルタのフィルタ係数を動的に算出する色差対応フィルタ係数生成ステップ（Ｓ５）と、色差対応フィルタ係数生成ステップ（Ｓ５）で得られた各注目画素に対応するフィルタ係数Ｗｇｔｉ,ｊに基づいて、色差成分画像データを加重平均フィルタ処理する色差対応ノイズ抑制処理ステップ（Ｓ６）をコンピュータに実行させる写真画像処理プログラムが、当該プログラムが記憶されたＣＤやＤＶＤ等の記憶媒体を介してインストールされていても構わない。
【０１３８】
さらに、ノイズ抑制強度σ ｘ,ｙを仮想分散値とし、濃度対応ノイズ抑制処理ステップ（Ｓ４）で得られた各加重平均画素と対応する注目画素との濃度差分値を変数ｔｍｐとする正規分布関数に基づいて重み係数ｗｇｔを算出し、当該重み係数ｗｇｔで各加重平均画素と対応する注目画素とを重み付け加算して補正濃度画像データを算出する補正処理ステップ（Ｓ４１）をコンピュータに実行させる写真画像処理プログラムが、当該プログラムが記憶されたＣＤやＤＶＤ等の記憶媒体を介してインストールされていても構わない。
【０１３９】
尚、処理対象画像データの端部領域では、実際の画素が存在しないため、端部画素と同じ画素値のダミー画素を生成して処理されるが、実際にプリントされる有効画像サイズは原画像サイズよりやや小さいサイズとなるので、プリント画像に影響を与えるものではない。
【０１４０】
第一の実施形態では、画素ブロックのブロック特性値に含まれる濃度ばらつきとして、各画素ブロックの平均濃度値を閾値として、当該閾値より高濃度の画素の平均濃度値である高濃度側平均濃度値と、当該閾値より低濃度の画素の平均濃度値である低濃度側平均濃度値との差分演算により算出される値を採用する場合を説明したが、本発明に適用可能な濃度ばらつきは、このような値に制限されるものではなく、画素ブロックを構成する画素群を母集団とする各画素の濃度値のばらつきを統計的手法等で表す値であればよく、標準偏差、平均絶対偏差、範囲等を採用することも可能である。
【０１４１】
尚、上述の濃度成分抽出部７０１及びＲＧＢ変換処理部７４で利用する変換式は、輝度成分をＲＧＢ色成分の平均値として導出しているが、例えば、Ｙ成分の変換式は、以下の式で示すように、適宜設計事項として変更しても構わない。
Ｙ    ＝    ０．２９９００Ｒ＋０．５８７００Ｇ＋０．１１４００Ｂ
【０１４２】
ただし、濃度成分抽出部７０１においてＲＧＢ色成分から輝度色差成分へ変換された結果のデータをＲＧＢ変換処理部７４に入力し、ＲＧＢ色成分に変換した結果が元のＲＧＢ色成分と同じになるように、整合性がとれた変換式に変更する必要があることは言うまでもない。
【０１４３】
以下、第一の実施形態とは別の実施形態として、第二の実施形態について説明する。
【０１４４】
第一の実施形態の構成に加えて、粒状ノイズ抑制処理部５１は、図９に示すように、濃度成分画像データを濃度成分画像データに対応したサイズの移動平均フィルタでフィルタ処理する移動平均フィルタ処理部８０と、ノイズ変動幅特定部７０で得られた各回帰直線に基づいて、移動平均フィルタ処理部８０で得られた各注目画素領域内の移動平均画素のノイズ変動幅を求める移動平均画素ノイズ変動幅特定部８１と、移動平均画素ノイズ変動幅特定部８１で得られた移動平均画素のノイズ変動幅に基づいて、濃度対応ノイズ抑制処理部７３で得られた各加重平均画素を補正するグラデーション補正処理部８２と、を備えている。
【０１４５】
グラデーション補正処理部８２は、各注目画素領域の濃度最大値と濃度最小値の濃度差分値を求める濃度差分値算出部８２１と、移動平均画素ノイズ変動幅特定部８１で得られた移動平均画素のノイズ変動幅と濃度差分値算出部８２１で得られた濃度差分値の比率に基づいて、各注目領域の平坦度を算出する平坦度算出部８２２と、平坦度算出部８２２で得られた平坦度を重み係数として、濃度対応ノイズ抑制処理部７３で得られた各加重平均画素と対応する移動平均画素とを重み付け加算してグラデーション再生濃度画像データを算出する重み付け加算処理部８２３を備えている。
【０１４６】
詳述すると、移動平均フィルタ処理部８０は、フィルタ係数を全て「１」とした濃度成分画像データに対応したサイズの移動平均フィルタを利用してフィルタ処理を行い、以下の式に示すように、濃度成分の移動平均画素の濃度値Ｙ２ａｖｅ ｘ,ｙを算出する。
【数９】
2010108008.tif 000011
【０１４７】
当該移動平均フィルタ処理された画像は、画像の絵柄に関係なく画像を一律に平滑化するものであるため、上述の加重平均フィルタ処理された画像よりも、視覚的に滑らかな画像となる。
【０１４８】
移動平均画素ノイズ変動幅特定部８１は、ノイズ変動幅特定部７０で得られた平坦画像領域から算出された二本の回帰直線が示す上限式と下限式に基づいて、移動平均フィルタ処理部８０で得られた各注目画素領域内の移動平均画素のノイズ変動幅を算出する。
【０１４９】
即ち、ノイズ変動幅特定部７０で得られた二本の回帰直線の説明変数は、複数の平坦画像領域の代表画素ブロックの平均濃度値であり、当該代表画素ブロック内の移動平均画素の濃度値ともいえる。したがって、当該説明変数が示すｘ軸の値を移動平均画素の濃度値であるとみなし、回帰直線が示す平坦画像領域のノイズの変動幅が、当該移動平均画素の濃度値に対するノイズの変動幅を示すものと仮定して算出する。
【０１５０】
移動平均画素の濃度値に対するノイズの変動幅の上限値ＭａｘＹ２ｘ,ｙと下限値ＭｉｎＹ２ｘ,ｙは、以下の式で示すことができる。
【数１０】
2010108008.tif 000012
【０１５１】
したがって、移動平均画素の濃度値に対するノイズの変動幅ＤｉｆＰ ｘ,ｙは、上限値ＭａｘＹ２ｘ,ｙと下限値ＭｉｎＹ２ｘ,ｙの差分値として、以下の式で示すことができる。
【数１１】
2010108008.tif 000013
【０１５２】
濃度差分値算出部８２１は、以下の式に示すように、濃度成分画像データから、上述のフィルタと同サイズの各注目画素領域の濃度最大値ｍａｘ ｘ,ｙと濃度最小値ｍｉｎ ｘ,ｙの濃度差分値ＤｉｆＲ ｘ,ｙを求める。
【数１２】
2010108008.tif 000014
【０１５３】
平坦度算出部８２２は、移動平均画素ノイズ変動幅特定部８１で得られた移動平均画素のノイズ変動幅ＤｉｆＰ ｘ,ｙと濃度差分値算出部８２１で得られた濃度差分値ＤｉｆＲ ｘ,ｙの比率に基づいて、以下の式に示すように、注目画素（ｘ,ｙ）に対する各注目領域の平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙを算出する。
【数１３】
2010108008.tif 000015
【０１５４】
移動平均画素のノイズ変動幅ＤｉｆＰ ｘ,ｙと注目画素の濃度差分値ＤｉｆＲ ｘ,ｙは、同サイズの画像領域における移動平均画素及び注目画素の最大濃度値と最小濃度値の濃度差分値を示しており、当該濃度差分値を対比させることにより、移動平均画素のノイズ変動幅ＤｉｆＰ ｘ,ｙを基準として平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙを算出している。
【０１５５】
即ち、注目領域の濃度差が移動平均画素の濃度差に近い場合は、注目領域は平坦な領域であると判定され、平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙは１に近づく値を示す。注目領域の濃度差が移動平均画素の濃度差とかけ離れている場合には、注目領域は平坦な領域ではないと判定され、平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙは０に近づく値を示す。また、平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙを上限値１で制限して、移動平均フィルタ処理された画像を超える程度に平坦な領域であると判定することを抑制している。
【０１５６】
重み付け加算処理部８２３は、平坦度算出部８２２で得られた平坦度を重み係数として、以下の式に示すように、濃度対応ノイズ抑制処理部７３で得られた各加重平均画素の濃度値Ｙａｖｅ ｘ,ｙと対応する移動平均画素の濃度値Ｙ２ａｖｅ ｘ,ｙとを重み付け加算して、平坦化された画像データであるグラデーション再生濃度画像データの濃度値Ｙａｖｅ´´ ｘ,ｙを算出する。
【数１４】
2010108008.tif 000016
【０１５７】
即ち、移動平均フィルタ処理により、ノイズ量に依存しない一律に平滑化された画像の濃度値の分布を基準として、平坦画像領域とそうではない画像領域を分離して、適切にノイズ抑制することができる。
【０１５８】
例えば、頬において光が当たっている箇所と当たっていない箇所が滑らかに変化する箇所のように、グラデーションの変化が小さく、画素値の小さな値が滑らかに変化する平坦画像領域に対しては、加重平均フィルタ処理を行わず、そうではない頬の輪郭部と背景画像を含む画像領域のように、平坦ではない画像領域に対しては、加重平均フィルタ処理を行うことにより、適度なノイズ抑制強度でノイズを抑制することができる。
【０１５９】
尚、重み付け加算処理部８２３は、平坦度算出部８２２で得られた平坦度を重み係数として、以下の式に示すように、補正処理部７３１で得られた補正濃度画像データの濃度値Ｙａｖｅ´ ｘ,ｙと対応する移動平均画素の濃度値Ｙ２ａｖｅ ｘ,ｙとを重み付け加算して、平坦化された画像データであるグラデーション再生濃度画像データの濃度値Ｙａｖｅ´´ ｘ,ｙを算出するように構成しても構わない。
【数１５】
2010108008.tif 000017
【０１６０】
即ち、平坦ではない領域に対しては、補正処理部７３１で得られる補正濃度画像データを利用することにより、過度にノイズを抑制しないように補正された強度で、ノイズを抑制することができる。
【０１６１】
この場合、ＲＧＢ変換処理部７４は、グラデーション補正処理部８２からグラデーション再生濃度画像データが入力されると、濃度対応ノイズ抑制処理部７３で得られた加重平均フィルタ処理後の濃度成分画像データと同様に取り扱い、当該濃度成分画像データと、色差対応ノイズ抑制処理部７６で得られた加重平均フィルタ処理後の色差成分画像データで構成される、ノイズ抑制後の新たな輝度色差成分を、上述の通り、ＲＧＢ色成分に変換して出力する。
【０１６２】
以下に、第二の実施形態における上述した粒状ノイズ抑制処理部５１による各処理の手順を、図１０に示すフローチャートに基づいて説明する。
【０１６３】
まず、写真画像データを構成する各画素のＲＧＢ成分データが、ノイズ変動幅特定部７０に入力されると、第一の実施形態で説明した通り、ノイズ変動幅特定ステップ（Ｓ１）が実行される。
【０１６４】
続いて、ノイズ変動幅特定部７０の濃度成分抽出ステップ（Ｓ１１）で得られた濃度成分画像データが移動平均画素ノイズ変動幅特定部８１に入力されると、移動平均フィルタ処理ステップが実行され、濃度成分画像データは、濃度成分画像データに対応したサイズの移動平均フィルタでフィルタ処理される（Ｓ２１）。
【０１６５】
続いて、移動平均画素ノイズ変動幅特定部８１の移動平均画素ノイズ変動幅特定ステップが実行され、ノイズ変動幅特定ステップ（Ｓ１）で得られた各回帰直線に基づいて、移動平均フィルタ処理ステップ（Ｓ２１）で得られた各注目画素領域内の移動平均画素の濃度値に対するノイズの変動幅ＤｉｆＰ ｘ,ｙが算出される（Ｓ２２）。
【０１６６】
続いて、グラデーション補正処理部８２のグラデーション補正処理ステップ（Ｓ２３）が実行される。
【０１６７】
グラデーション補正処理ステップ（Ｓ２３）は、濃度差分値算出ステップ（Ｓ２３１）と、平坦度算出ステップ（Ｓ２３２）と、重み付け加算処理ステップ（Ｓ２３３）で構成される。
【０１６８】
まず、濃度差分値算出ステップが実行され、濃度成分抽出ステップ（Ｓ１１）で得られた濃度成分画像データから、移動平均フィルタ処理ステップ（Ｓ２１）で利用した移動平均フィルタと同サイズの各注目画素領域の濃度最大値ｍａｘ ｘ,ｙと濃度最小値ｍｉｎ ｘ,ｙの濃度差分値ＤｉｆＲ ｘ,ｙが算出される（Ｓ２３１）。
【０１６９】
続いて、平坦度算出ステップが実行され、移動平均画素ノイズ変動幅特定ステップ（Ｓ２２）で得られた移動平均画素のノイズ変動幅ＤｉｆＰ ｘ,ｙと濃度差分値算出ステップ（Ｓ２３１）で得られた濃度差分値ＤｉｆＲ ｘ,ｙの比率に基づいて、注目画素（ｘ,ｙ）に対する各注目領域の平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙが算出される。（Ｓ２３２）。
【０１７０】
続いて、重み付け加算処理ステップが実行され、平坦度算出ステップ（Ｓ２３２）で得られた平坦度を重み係数として、濃度対応ノイズ抑制処理ステップ（Ｓ４）で得られた各加重平均画素の濃度値Ｙａｖｅ ｘ,ｙ、あるいは、補正処理部ステップ（Ｓ４１）で得られた補正濃度画像データの濃度値Ｙａｖｅ´ ｘ,ｙと対応する移動平均画素の濃度値Ｙ２ａｖｅ ｘ,ｙとが重み付け加算され、平坦化された画像データの濃度値Ｙａｖｅ´´ ｘ,ｙが算出される（Ｓ２３３）。
【０１７１】
続いて、濃度対応ノイズ抑制処理ステップ（Ｓ４）で加重平均フィルタ処理された濃度成分画像データ、あるいは、補正処理ステップ（Ｓ４１）で算出された補正濃度画像データ、あるいは、グラデーション補正処理ステップ（Ｓ２３）で算出されたグラデーション再生濃度画像データと、色差対応ノイズ抑制処理ステップ（Ｓ６）で得られた加重平均フィルタ処理後の色差成分画像データがＲＧＢ変換処理部７４に入力されると、ＲＧＢ変換処理ステップが実行され、当該入力された、濃度成分画像データ、あるいは、補正濃度画像データと、色差成分画像データと、で構成されるノイズ抑制後の新たな輝度色差成分が、ＲＧＢ色成分に変換されて出力され、粒状ノイズ抑制処理部５１による粒状ノイズの抑制処理が終了される（Ｓ８）。
【０１７２】
尚、上述の各ステップにおいて、各ステップ間の入出力データをメモリ４１に保存するよう構成し、各ステップは、各ステップ間で直接データの入出力を行わずに、メモリ４１に保存されたデータにアクセスするよう構成しても構わない。
【０１７３】
上述した粒状ノイズ抑制処理部５１による各処理は、コントローラ３３に備えたハードディスクにインストールされた本発明の写真画像処理プログラムが実行されることにより実現される。
【０１７４】
つまり、加重平均フィルタ処理により写真画像データに含まれる粒状のノイズを抑制する写真画像処理プログラムであって、写真画像データを構成する濃度成分画像データから複数の平坦画像領域を探索し、各平坦画像領域の平均濃度値を説明変数とし最大濃度値及び最小濃度値の夫々を従属変数とする二本の回帰直線により写真画像データに含まれる粒状のノイズの変動幅を求めるノイズ変動幅特定ステップ（Ｓ１）と、ノイズ変動幅特定ステップ（Ｓ１）で得られた各回帰直線に基づいて、任意の注目画素に対するノイズ変動幅を求め、当該ノイズ変動幅からノイズ抑制強度を算出するノイズ抑制強度算出ステップ（Ｓ２）と、当該注目画素と周辺画素の濃度差分値と、ノイズ抑制強度算出ステップ（Ｓ２）で得られたノイズ抑制強度に基づいて、濃度成分画像データに対応したサイズの加重平均フィルタのフィルタ係数を動的に算出する濃度対応フィルタ係数生成ステップ（Ｓ３）と、濃度対応フィルタ係数生成ステップ（Ｓ３）で得られた各注目画素に対応するフィルタ係数に基づいて、濃度成分画像データを加重平均フィルタ処理する濃度対応ノイズ抑制処理ステップ（Ｓ４）と、濃度成分画像データを濃度成分画像データに対応したサイズの移動平均フィルタでフィルタ処理する移動平均フィルタ処理ステップ（Ｓ２１）と、ノイズ変動幅特定ステップ（Ｓ１）で得られた各回帰直線に基づいて、移動平均フィルタ処理ステップ（Ｓ２１）で得られた各注目画素領域内の移動平均画素のノイズ変動幅を求める移動平均画素ノイズ変動幅特定ステップ（Ｓ２２）と、移動平均画素ノイズ変動幅特定ステップ（Ｓ２２）で得られた移動平均画素のノイズ変動幅に基づいて、濃度対応ノイズ抑制処理ステップ（Ｓ４）で得られた各加重平均画素を補正するグラデーション補正処理ステップ（Ｓ２３）と、をコンピュータに実行させる写真画像処理プログラムが記憶されたＣＤやＤＶＤ等の記憶媒体を介してインストールされている。
【０１７５】
さらに、グラデーション補正処理ステップ（Ｓ２３）は、各注目画素領域の濃度最大値と濃度最小値の濃度差分値を求める濃度差分値算出ステップ（Ｓ２３１）と、移動平均画素ノイズ変動幅特定ステップ（Ｓ２２）で得られた移動平均画素のノイズ変動幅と濃度差分値算出ステップ（Ｓ２３１）で得られた濃度差分値の比率に基づいて、各注目領域の平坦度を算出する平坦度算出ステップ（Ｓ２３２）と、を含み、平坦度算出ステップ（Ｓ２３２）で得られた平坦度を重み係数として、濃度対応ノイズ抑制処理ステップ（Ｓ４）で得られた各加重平均画素と対応する移動平均画素とを重み付け加算してグラデーション再生濃度画像データを算出するコンピュータに実行させる写真画像処理プログラムが、当該プログラムが記憶されたＣＤやＤＶＤ等の記憶媒体を介してインストールされていても構わない。
【０１７６】
尚、処理対象画像データの端部領域では、実際の画素が存在しないため、端部画素と同じ画素値のダミー画素を生成して処理されるが、実際にプリントされる有効画像サイズは原画像サイズよりやや小さいサイズとなるので、プリント画像に影響を与えるものではない。
【０１７７】
尚、第二の実施形態の構成を変更し、平坦度算出部８２２は、以下の式に示すように、平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙを調整するための平坦率調整値ＦｌａｔＳｈｉｆｔと平坦率正規化係数ＦｌａｔＢａｌａｎｃｅの二つの調整用変数を設けて、平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙを算出するように構成しても構わない。
【数１６】
2010108008.tif 000018
【０１７８】
尚、第一または第二の実施形態により出力された結果画像よりも、ノイズ抑制強度を弱めた画像を得たい場合に備えて、ＲＧＢ変換処理部７４は、結果画像の輝度色差成分をＲＧＢ色成分に変換する前に、以下の式に基づいて、当該結果画像の各注目画素の画素値と対応する原画像の画素値とを設定可能な変数ｒｔｏで融合処理するように構成しても構わない。
【数１７】
2010108008.tif 000019
【０１７９】
尚、Ｙａｎｓ ｘ,ｙは、当該融合処理した結果の画像データの注目画素（ｘ,ｙ）の輝度成分を示し、Ｃｂａｎｓ ｘ,ｙ及びＣｒａｎｓ ｘ,ｙは、当該融合処理した結果の画像データの注目画素（ｘ,ｙ）の色差成分を示す。また、Ｙｒｅｓ ｘ,ｙは、当該融合処理前の本発明により得られるノイズ抑制された画像データの注目画素（ｘ,ｙ）の輝度成分を示し、Ｃｂｒｅｓ ｘ,ｙ及びＣｒｒｅｓ ｘ,ｙは、当該融合処理前の本発明により得られるノイズ抑制された画像データの注目画素（ｘ,ｙ）の色差成分を示す。
【０１８０】
以下、第一及び第二の実施形態とは別の実施形態として、第三の実施形態について説明する。
【０１８１】
図１１に示すように、第二の実施形態の構成のグラデーション補正部８２の構成に代えて、グラデーション補正部８２は、各注目画素領域の濃度値の標準偏差を求める濃度標準偏差算出部８２４と、移動平均画素ノイズ変動幅特定部８１で得られた移動平均画素のノイズ変動幅に基づいて、各注目画素領域内の移動平均画素の濃度値の仮想標準偏差を算出し、仮想標準偏差と濃度標準偏差算出部８２４で得られた濃度値の標準偏差との比率に基づいて、各注目領域の平坦度を算出する平坦度算出部８２２と、平坦度算出部８２２で得られた平坦度を重み係数として、濃度対応ノイズ抑制処理部７３で得られた各加重平均画素と対応する移動平均画素とを重み付け加算してグラデーション再生濃度画像データを算出する重み付け加算処理部８２３を備えている。
【０１８２】
濃度標準偏差算出部８２４は、以下の式に示すように、濃度成分画像データから、上述の濃度対応ノイズ抑制処理部７３で用いられるフィルタと同サイズの各注目画素領域の濃度値の標準偏差Ｓｔｄｘ,ｙを求める。
【数１８】
2010108008.tif 000020
【０１８３】
ここで、ｎＹは、注目画素と周辺画素の最大距離を示し、例えば、フィルタサイズが１１×１１の場合には、ｎＹ＝１１÷２≒５となる。また、ｓｕｍＹは、注目画素領域内の全画素の濃度値を累積加算した値を示し、ｓｕｍＹＹは、注目画素領域内の全画素の濃度値の自乗値を累積加算した値を示す。
【０１８４】
標準偏差Ｓｔｄｘ,ｙは、注目画素領域内の濃度値のばらつきを表しており、注目画素領域に何らかの特徴的な絵柄がある場合は大きくなり、注目画素領域が平坦な画像である場合は小さくなる。ただし、標準偏差Ｓｔｄｘ,ｙは、注目画素領域にノイズが多い場合にも大きくなる傾向があり、標準偏差Ｓｔｄｘ,ｙの大小が、ノイズに因るものなのか絵柄に因るものなのかを見極める必要がある。
【０１８５】
そこで、平坦度算出部８２２は、上述の第二の実施形態の構成に代えて、まず、上述のノイズ抑制強度算出部７１と同様に、移動平均画素のノイズ変動幅ＤｉｆＰ ｘ,ｙ内の濃度値の分布が正規分布に従うものとみなすことにより、注目画素領域の絵柄に依存しない移動平均画素の濃度値の仮想標準偏差σ´ｘ,ｙを算出する。
【０１８６】
具体的には、平坦度算出部８２２は、上述のノイズ抑制強度算出部７１で定義されたσ化係数ＳｉｇｍａＣｏｅｆと最大許容値Ｌｅｖｅｌと最小許容値１とを利用して、以下の式に示すように、仮想標準偏差σ´ｘ,ｙを算出する。
【数１９】
2010108008.tif 000021
【０１８７】
さらに、平坦度算出部８２２は、当該算出された仮想標準偏差σ´ｘ,ｙと濃度標準偏差算出部８２４で得られた濃度値の標準偏差Ｓｔｄｘ,ｙとの比率に基づいて、以下の式に示すように、注目画素（ｘ,ｙ）に対する各注目領域の平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙを算出する。
【数２０】
2010108008.tif 000022
【０１８８】
移動平均画素の仮想標準偏差σ´ｘ,ｙと注目画素の標準偏差Ｓｔｄｘ,ｙは、同サイズの画像領域における移動平均画素の濃度値及び注目画素の濃度値の標準偏差を示しており、当該標準偏差を対比させることにより、平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙを算出している。
【０１８９】
即ち、注目領域の濃度値の標準偏差Ｓｔｄｘ,ｙが移動平均画素の仮想標準偏差σ´ｘ,ｙに近い場合は、注目領域は平坦な領域であると判定され、平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙは１に近づく値を示す。注目領域の濃度差が移動平均画素の濃度差とかけ離れている場合には、注目領域は平坦な領域ではないと判定され、平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙは０に近づく値を示す。また、平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙを上限値１で制限して、移動平均フィルタ処理された画像を超える程度に平坦な領域であると判定することを抑制している。
【０１９０】
重み付け加算処理部８２３は、第二の実施形態と同様に構成され、平坦度算出部８２２で得られた平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙを重み係数として、数１４または数１５で示した式に従い、濃度対応ノイズ抑制処理部７３で得られた各加重平均画素の濃度値Ｙａｖｅ ｘ,ｙまたは補正処理部７３１で得られた補正濃度画像データの濃度値Ｙａｖｅ´ ｘ,ｙと対応する移動平均画素の濃度値Ｙ２ａｖｅ ｘ,ｙとを重み付け加算して、平坦化された画像データであるグラデーション再生濃度画像データの濃度値Ｙａｖｅ´´ ｘ,ｙを算出する。
【０１９１】
即ち、移動平均フィルタ処理によりノイズ量に依存しない一律に平滑化された画像の濃度値の分布をさらに滑らかに仮想し、当該仮想した濃度値の分布を基準として、平坦画像領域とそうではない画像領域を滑らかに分離して、適切にノイズ抑制することができるようになる。
【０１９２】
例えば、頬において光が当たっている箇所と当たっていない箇所が滑らかに変化する箇所のように、グラデーションの変化が小さく、画素値の小さな値が滑らかに変化する平坦画像領域に対しては、加重平均フィルタ処理を行わず、そうではない頬の輪郭部と背景画像を含む画像領域のように、平坦ではない画像領域に対しては、加重平均フィルタ処理を行うことにより、適度なノイズ抑制強度でノイズを抑制することができる。
【０１９３】
また、第二の実施形態の構成によりノイズ抑制処理した結果の画像を凝視した場合には、フィルタサイズの画素領域により画像が区分されたような模様が見られることがあるが、第三の実施形態の構成により、平坦画像領域とそうではない画像領域がさらに滑らかに分離されてノイズ抑制処理されることにより、上述の模様が見られることが軽減されるようになる。
【０１９４】
また、第二の実施形態におけるグラデーション補正部８２の濃度差分値算出部８２１では、数１２の式に示すように、各注目画素領域の濃度最大値ｍａｘ ｘ,ｙと濃度最小値ｍｉｎ ｘ,ｙを特定する場合に、各注目画素領域の全画素に対して、当該特定する処理中の濃度最大値及び濃度最小値と当該処理対象の画素の濃度値の大小を比較し、当該比較結果に基づいて濃度最大値及び濃度最小値を逐次更新する処理が必要である。さらに、濃度差分値算出部８２１では、当該逐次更新処理により算出された濃度最大値ｍａｘ ｘ,ｙと濃度最小値ｍｉｎ ｘ,ｙを利用して、濃度差分値ＤｉｆＲ ｘ,ｙを算出できる。
【０１９５】
しかし、第三の実施形態におけるグラデーション補正部８２の濃度標準偏差算出部８２４では、数１８の式に示すように、各注目画素領域の全画素に対して、処理対象の画素の濃度値を累積的に加算及び乗算する処理を実行して、濃度値の標準偏差Ｓｔｄｘ,ｙを算出することができる。したがって、第二の実施形態と比較して、処理速度を向上することができるようになる。
【０１９６】
以下に、第三の実施形態における上述した粒状ノイズ抑制処理部５１による各処理の手順を、図１２に示すフローチャートに基づいて説明する。
【０１９７】
まず、写真画像データを構成する各画素のＲＧＢ成分データが、ノイズ変動幅特定部７０に入力されると、第二の実施形態で説明した通り、ノイズ変動幅特定ステップ（Ｓ１）から移動平均画素ノイズ変動幅特定ステップ（Ｓ２２）の処理が実行される。
【０１９８】
続いて、グラデーション補正処理部８２のグラデーション補正処理ステップ（Ｓ２３）が実行される。
【０１９９】
グラデーション補正処理ステップ（Ｓ２３）は、濃度標準偏差算出ステップ（Ｓ２３４）と、平坦度算出ステップ（Ｓ２３２）と、重み付け加算処理ステップ（Ｓ２３３）で構成される。
【０２００】
まず、濃度標準偏差算出ステップが実行され、濃度成分抽出ステップ（Ｓ１１）で得られた濃度成分画像データから、濃度対応ノイズ抑制処理ステップ（Ｓ４）で用いられるフィルタと同サイズの各注目画素領域の濃度値の標準偏差Ｓｔｄｘ,ｙが算出される（Ｓ２３４）。
【０２０１】
続いて、平坦度算出ステップが実行され、移動平均画素ノイズ変動幅特定ステップ（Ｓ２２）で得られた移動平均画素のノイズ変動幅ＤｉｆＰ ｘ,ｙに基づいて、各注目画素領域内の移動平均画素の濃度値の仮想標準偏差σ´ｘ,ｙが算出され、仮想標準偏差σ´ｘ,ｙと濃度標準偏差算出ステップ（Ｓ２３４）で得られた濃度値の標準偏差Ｓｔｄｘ,ｙとの比率に基づいて、各注目領域の平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙが算出される。（Ｓ２３２）。
【０２０２】
続いて、重み付け加算処理ステップが実行され、平坦度算出ステップ（Ｓ２３２）で得られた平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙを重み係数として、濃度対応ノイズ抑制処理ステップ（Ｓ４）で得られた各加重平均画素の濃度値Ｙａｖｅ ｘ,ｙ、あるいは、補正処理部ステップ（Ｓ４１）で得られた補正濃度画像データの濃度値Ｙａｖｅ´ ｘ,ｙと対応する移動平均画素の濃度値Ｙ２ａｖｅ ｘ,ｙとが重み付け加算され、平坦化された画像データの濃度値Ｙａｖｅ´´ ｘ,ｙが算出される（Ｓ２３３）。
【０２０３】
続いて、第二の実施形態で説明した通り、濃度対応ノイズ抑制処理ステップ（Ｓ４）からＲＧＢ変換処理ステップ（Ｓ８）の処理が実行される。
【０２０４】
上述した第三の実施形態における粒状ノイズ抑制処理部５１による各処理は、コントローラ３３に備えたハードディスクにインストールされた本発明の写真画像処理プログラムが実行されることにより実現される。
【０２０５】
つまり、加重平均フィルタ処理により写真画像データに含まれる粒状のノイズを抑制する写真画像処理プログラムであって、写真画像データを構成する濃度成分画像データから複数の平坦画像領域を探索し、各平坦画像領域の平均濃度値を説明変数とし最大濃度値及び最小濃度値の夫々を従属変数とする二本の回帰直線により写真画像データに含まれる粒状のノイズの変動幅を求めるノイズ変動幅特定ステップ（Ｓ１）と、ノイズ変動幅特定ステップ（Ｓ１）で得られた各回帰直線に基づいて、任意の注目画素に対するノイズ変動幅を求め、当該ノイズ変動幅からノイズ抑制強度を算出するノイズ抑制強度算出ステップ（Ｓ２）と、当該注目画素と周辺画素の濃度差分値と、ノイズ抑制強度算出ステップ（Ｓ２）で得られたノイズ抑制強度に基づいて、濃度成分画像データに対応したサイズの加重平均フィルタのフィルタ係数を動的に算出する濃度対応フィルタ係数生成ステップ（Ｓ３）と、濃度対応フィルタ係数生成ステップ（Ｓ３）で得られた各注目画素に対応するフィルタ係数に基づいて、濃度成分画像データを加重平均フィルタ処理する濃度対応ノイズ抑制処理ステップ（Ｓ４）と、濃度成分画像データを濃度成分画像データに対応したサイズの移動平均フィルタでフィルタ処理する移動平均フィルタ処理ステップ（Ｓ２１）と、ノイズ変動幅特定ステップ（Ｓ１）で得られた各回帰直線に基づいて、移動平均フィルタ処理ステップ（Ｓ２１）で得られた各注目画素領域内の移動平均画素のノイズ変動幅を求める移動平均画素ノイズ変動幅特定ステップ（Ｓ２２）と、移動平均画素ノイズ変動幅特定ステップ（Ｓ２２）で得られた移動平均画素のノイズ変動幅に基づいて、濃度対応ノイズ抑制処理ステップ（Ｓ４）で得られた各加重平均画素を補正するグラデーション補正処理ステップ（Ｓ２３）と、をコンピュータに実行させる写真画像処理プログラムが記憶されたＣＤやＤＶＤ等の記憶媒体を介してインストールされている。
【０２０６】
さらに、グラデーション補正処理ステップ（Ｓ２３）は、各注目画素領域の濃度値の標準偏差を求める濃度標準偏差算出ステップ（Ｓ２３４）と、移動平均画素ノイズ変動幅特定ステップ（Ｓ２２）で得られた移動平均画素のノイズ変動幅ＤｉｆＰ ｘ,ｙに基づいて、各注目画素領域内の移動平均画素の濃度値の仮想標準偏差σ´ｘ,ｙを算出し、仮想標準偏差σ´ｘ,ｙと濃度標準偏差算出ステップ（Ｓ２３４）で得られた濃度値の標準偏差Ｓｔｄｘ,ｙとの比率に基づいて、各注目領域の平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙを算出する平坦度算出ステップ（Ｓ２３２）と、を含み、平坦度算出ステップ（Ｓ２３２）で得られた平坦度を重み係数として、濃度対応ノイズ抑制処理ステップ（Ｓ４）で得られた各加重平均画素と対応する移動平均画素とを重み付け加算してグラデーション再生濃度画像データを算出するコンピュータに実行させる写真画像処理プログラムが、当該プログラムが記憶されたＣＤやＤＶＤ等の記憶媒体を介してインストールされていても構わない。
【０２０７】
尚、第三の実施形態の構成を変更し、平坦度算出部８２２は、以下の式に示すように、平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙを調整するための平坦率調整値ＦｌａｔＳｈｉｆｔと平坦率正規化係数ＦｌａｔＢａｌａｎｃｅの二つの調整用変数を設けて、平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙを算出するように構成しても構わない。
【数２１】
2010108008.tif 000023
【０２０８】
尚、上述の第二または第三の実施形態における重み付け加算処理部８２３は、平坦度算出部８２２で得られた平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙを重み係数として、濃度成分画像データの濃度値ｄａｔ ｘ,ｙと対応する移動平均画素の濃度値Ｙ２ａｖｅ ｘ,ｙとを重み付け加算して、平坦化された画像データであるグラデーション再生濃度画像データの濃度値Ｙａｖｅ´´ ｘ,ｙを算出するように構成しても構わない。
【０２０９】
この場合、上述の第二または第三の実施形態におけるグラデーション補正処理ステップ（Ｓ２３）の重み付け加算処理ステップ（Ｓ２３３）において、平坦度算出ステップ（Ｓ２３２）で得られた平坦度ＦｌａｔＲａｔｉｏ ｘ,ｙを重み係数として、濃度成分画像データの濃度値ｄａｔ ｘ,ｙと対応する移動平均画素の濃度値Ｙ２ａｖｅ ｘ,ｙとが重み付け加算され、平坦化された画像データの濃度値Ｙａｖｅ´´ ｘ,ｙが算出されるよう構成する必要があることは言うまでもない。
【０２１０】
尚、上述した実施形態は、本発明の一例に過ぎず、本発明の作用効果を奏する範囲において各ブロックの具体的構成等を適宜変更設計できることは言うまでもない。
【図面の簡単な説明】
【０２１１】
【図１】写真画像処理装置の外観説明図
【図２】写真プリンタの説明図
【図３】写真画像処理装置の機能ブロック構成図
【図４】粒状ノイズ抑制処理部の機能ブロック構成図
【図５】ノイズ変動幅特定ステップの説明図であり、（ａ）は所定サイズの画素ブロックに分割された濃度成分画像データを示す説明図、（ｂ）は濃度成分画像データの各画素ブロックを示す説明図、（ｃ）は複数の区分に分割された濃度範囲に画素ブロックを割り付けた一例を示す説明図、（ｄ）はノイズ変動幅を特定する一例を示す説明図
【図６】平坦画像領域から算出された回帰直線の相関係数に応じた補正の説明図であり、（ａ）は、回帰直線の傾きを補正を示す説明図、（ｂ）は、回帰直線の傾きと切片の補正を示す説明図
【図７】所定の画像領域における濃度値の分布の説明図であり、（ａ）は、平坦画像領域のみを含む画像領域の正規化された濃度値の分布を示す説明図、（ｂ）は、平坦画像領域と輪郭画像領域をともに含む画像領域の正規化された濃度値の分布を示す説明図
【図８】粒状ノイズ抑制処理部の処理手順の一例を説明するためのフローチャート
【図９】粒状ノイズ抑制処理部の図４とは別の一例を示す機能ブロック構成図
【図１０】図９の構成に基づいた粒状ノイズ抑制処理部の処理手順の一例を説明するためのフローチャート
【図１１】粒状ノイズ抑制処理部の図４及び図９とは別の一例を示す機能ブロック構成図
【図１２】図１１の構成に基づいた粒状ノイズ抑制処理部の処理手順の一例を説明するためのフローチャート
【符号の説明】
【０２１２】
１：写真画像処理装置
５１：粒状ノイズ抑制処理部
７０：ノイズ変動幅特定部
７１：ノイズ抑制強度算出部
７２：濃度対応フィルタ係数生成部
７３：濃度対応ノイズ抑制処理部
７４：ＲＧＢ変換処理部
７５：色差対応フィルタ係数生成部
７６：色差対応ノイズ抑制処理部
８０：移動平均フィルタ処理部
８１：移動平均画素ノイズ変動幅特定部
８２：グラデーション補正処理部
７０１：濃度成分抽出部（ノイズ変動幅特定部）
７０２：ブロック特性値算出部（ノイズ変動幅特定部）
７０３：画素ブロック割付部（ノイズ変動幅特定部）
７０４：画素ブロック抽出部（ノイズ変動幅特定部）
７０５：回帰直線算出部（ノイズ変動幅特定部）
７３１：補正処理部（濃度対応ノイズ抑制処理部）
８２１：濃度差分値算出部（グラデーション補正処理部）
８２２：平坦度算出部（グラデーション補正処理部）
８２３：重み付け加算処理部（グラデーション補正処理部）
８２４：濃度標準偏差算出部（グラデーション補正処理部）
Ｓ１：ノイズ変動幅特定ステップ
Ｓ２：ノイズ抑制強度算出ステップ
Ｓ３：濃度対応フィルタ係数生成ステップ
Ｓ４：濃度対応ノイズ抑制処理ステップ
Ｓ５：色差対応フィルタ係数生成ステップ
Ｓ６：色差対応ノイズ抑制処理ステップ
Ｓ８：ＲＧＢ変換処理ステップ
Ｓ１１：濃度成分抽出ステップ（ノイズ変動幅特定ステップ）
Ｓ１２：ブロック特性値算出ステップ（ノイズ変動幅特定ステップ）
Ｓ１３：画素ブロック割付ステップ（ノイズ変動幅特定ステップ）
Ｓ１４：画素ブロック抽出ステップ（ノイズ変動幅特定ステップ）
Ｓ１５：回帰直線算出ステップ（ノイズ変動幅特定ステップ）
Ｓ２１：移動平均フィルタ処理ステップ
Ｓ２２：移動平均画素ノイズ変動幅特定ステップ
Ｓ２３：グラデーション補正処理ステップ
Ｓ４１：補正処理ステップ（濃度対応ノイズ抑制処理ステップ）
Ｓ２３１：濃度差分値算出ステップ（グラデーション補正処理ステップ）
Ｓ２３２：平坦度算出ステップ（グラデーション補正処理ステップ）
Ｓ２３３：重み付け加算処理ステップ（グラデーション補正処理ステップ）
Ｓ２３４：濃度標準偏差算出ステップ（グラデーション補正処理ステップ）
【図１】
2010108008.tif 000024
【図２】
2010108008.tif 000025
【図３】
2010108008.tif 000026
【図４】
2010108008.tif 000027
【図５】
2010108008.tif 000028
【図６】
2010108008.tif 000029
【図７】
2010108008.tif 000030
【図８】
2010108008.tif 000031
【図９】
2010108008.tif 000032
【図１０】
2010108008.tif 000033
【図１１】
2010108008.tif 000034
【図１２】
2010108008.tif 000035