特許 6249173 画像処理方法、及び画像処理装置

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6249173

(24)【登録日】2017年12月1日

(45)【発行日】2017年12月20日

(54)【発明の名称】画像処理方法、及び画像処理装置

(51)【国際特許分類】

   G06T 5/00 20060101AFI20171211BHJP

【ＦＩ】

   G06T5/00 705

【請求項の数】10

【全頁数】16

(21)【出願番号】特願2014-512625(P2014-512625)

(86)(22)【出願日】2013年4月24日

(86)【国際出願番号】JP2013061969

(87)【国際公開番号】WO2013161839

(87)【国際公開日】20131031

【審査請求日】2016年3月3日

(31)【優先権主張番号】特願2012-100910(P2012-100910)

(32)【優先日】2012年4月26日

(33)【優先権主張国】JP

(73)【特許権者】

【識別番号】000004237

【氏名又は名称】日本電気株式会社

(74)【代理人】

【識別番号】100079005

【弁理士】

【氏名又は名称】宇高 克己

(72)【発明者】

【氏名】先崎 健太

(72)【発明者】

【氏名】塚田 正人

(72)【発明者】

【氏名】河口 裕明

【審査官】 板垣 有紀

(56)【参考文献】

【文献】 特開２００７−０１８３７９（ＪＰ，Ａ）

【文献】 特開２００２−１８３７２７（ＪＰ，Ａ）

【文献】 特開２０１１−０４１１８３（ＪＰ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０６Ｔ ５／００

Ｇ０６Ｔ ５／２０

Ｈ０４Ｎ １／４０９

Ｈ０４Ｎ ５／２１

Ｈ０４Ｎ ５／２３２

(57)【特許請求の範囲】

【請求項1】

注目画素を含み、範囲が順次狭くなる多階層の領域毎に画素の画素統計値とエッジ情報とを求め、
補正処理を行う当該階層の領域の画素統計値と前記当該階層の領域より広い階層の領域の画素統計値との差分情報を、前記当該階層の領域より広い階層の領域におけるエッジ情報を用いて補正し、
補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値を用いて、前記当該階層の領域の画素統計値を補正し、
前記当該階層の領域より広い階層の領域におけるエッジ情報を用いて、補正前の前記当該階層の領域の画素統計値と、補正後の前記当該階層の領域の画素統計値の差分情報とを再補正し、
最大範囲から最小範囲の領域になるまで、順次各階層において、前記当該階層の領域の画素統計値の補正と再補正を繰り返すことで、前記注目画素を補正する
画像処理方法。

【請求項2】

前記再補正は、
前記注目画素がエッジに含まれる場合には補正量を小さくし、
前記注目画素が平坦領域に含まれる場合には補正量を大きくし、
前記注目画素がテクスチャ領域に含まれる場合にはエッジ情報に応じた補正量により、補正を行う
請求項１に記載の画像処理方法。

【請求項3】

当該階層の領域の画素統計値と、当該階層の領域より広い階層の領域の画素統計値と、当該階層の領域より広い階層の領域におけるエッジ情報とを用いて、前記差分情報を算出する
請求項１又は請求項２に記載の画像処理方法。

【請求項4】

当該階層の領域より広い階層の領域におけるエッジ情報が所定の閾値を超える場合、当該階層の画素統計値の補正を行わない
請求項３に記載の画像処理方法。

【請求項5】

補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値とを用いて、当該階層の領域の画素統計値を補正する際、補正の強度を階層毎に変化させる
請求項１から請求項４のいずれかに記載の画像処理方法。

【請求項6】

補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値とを用いて、当該階層の領域の画素統計値を補正する際、当該階層の領域のノイズによる画素値の変動量に応じて補正の強度を変化させる
請求項５に記載の画像処理方法。

【請求項7】

補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値とを用いて、当該階層の領域の画素統計値を補正する際、当該階層の領域より広い領域のエッジ量に応じて補正の強度を変化させる
請求項５に記載の画像処理方法。

【請求項8】

補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値とを用いて、当該階層の領域の画素統計値を補正する際、前記差分情報がゼロに近づくほどゼロに近い値を出力し、前記差分情報が大きくなるほど補正を弱める補正を行う
請求項１から請求項７のいずれかに記載の画像処理方法。

【請求項9】

注目画素を含み、範囲が順次狭くなる多階層の領域毎に、その領域の画素の画素統計値を算出する画素統計値算出手段と、
注目画素を含み、範囲が順次狭くなる多階層の領域毎に、エッジ情報を算出するエッジ情報算出手段と、
補正処理を行う当該階層の領域の画素統計値と前記当該階層の領域より広い階層の領域の画素統計値との差分情報を、前記当該階層の領域より広い階層の領域におけるエッジ情報を用いて補正し、補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値を用いて、前記当該階層の領域の画素統計値を補正し、前記当該階層の領域より広い階層の領域におけるエッジ情報を用いて、補正前の前記当該階層の領域の画素統計値と、補正後の前記当該階層の領域の画素統計値の差分情報とを再補正し、最大範囲から最小範囲の領域になるまで、順次各階層において、前記当該階層の領域の画素統計値の補正と再補正を繰り返すことで、前記注目画素を補正する補正手段と
を有する画像処理装置。

【請求項10】

前記補正手段は、
前記注目画素がエッジに含まれる場合には補正量を小さくし、
前記注目画素が平坦領域に含まれる場合には補正量を大きくし、
前記注目画素がテクスチャ領域に含まれる場合にはエッジ情報に応じた補正量により、補正を行う
請求項９に記載の画像処理装置。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、画像処理方法及び画像処理装置に関する。

【背景技術】

【0002】

画像処理技術の中でも、画像に含まれるランダムなノイズを低減する技術は、撮像した画像をより鮮明に再現するためには、欠かせない技術である。ランダムノイズを低減する技術としては、例えば、特許文献１に開示された技術がある。

【0003】

特許文献１に記載された技術は、入力画像処理回路から出力されたカラーデジタル信号の主走査方向の任意の注目画素ｉについて、所定の数式に基づき移動平均画素数ｎを算出する複数の演算回路と、注目画素ｉ及び前後ｎ画素の参照画素ｊを選択出力する複数のビット選択回路と、注目画素ｉの出力レベルと参照画素ｊのそれぞれの出力レベルとの差分の絶対値を算出する複数の差分回路と、注目画素ｉを出力するとともに、複数の差分回路から出力される値と閾値記憶メモリ内の所定の閾値とを比較し、該比較結果に基づいて参照画素ｊを出力する複数の判定回路と、複数の判定回路から出力された出力信号の移動平均化処理を行う複数の演算回路とから構成される。

【0004】

すなわち、注目画素ｉの出力レベルと参照画素ｊの出力レベルとの差分の絶対値が閾値以下である場合にのみ、該参照画素ｊが移動平均化処理に加えられるような構成としているため、上記差分の絶対値が閾値以上に急峻に変化している部分については移動平均化処理から除外され、これにより、ノイズ成分を効果的に除去することができる。

【先行技術文献】

【特許文献】

【0005】

【特許文献1】特開２００２−５７９００号公報

【発明の概要】

【発明が解決しようとする課題】

【0006】

しかしながら、特許文献１の技術は、平滑化のフィルタサイズ以上の周期性を有する低周波ノイズを除去することができなかった。

【0007】

特に、画像中のエッジ成分やテクスチャ成分を保存しつつ、画像中のノイズを効果的に除去することができなかった。

【0008】

そこで、本発明は上記課題に鑑みて発明されたものであって、その目的は画像中のエッジ成分やテクスチャ成分を保存しつつ、画像中のノイズを効果的に除去できる画像処理方法及び画像処理装置を提供することにある。

【課題を解決するための手段】

【0009】

本発明は、注目画素を含み、範囲が順次狭くなる多階層の領域毎に画素の画素統計値とエッジ情報とを求め、当該階層の領域の画素統計値と当該階層の領域より広い階層の領域の画素統計値との差分情報を、前記エッジ情報を用いて補正し、補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値を用いて、当該階層の領域の画素統計値を補正し、前記エッジ情報を用いて、補正前の前記当該階層の領域の画素統計値と、補正後の前記当該階層の領域の画素統計値の差分情報とを再補正し、最大範囲から最小範囲の領域になるまで、順次各階層において、前記当該階層の領域の画素統計値の補正と再補正を繰り返すことで、前記注目画素を補正する画像処理方法である。

【0010】

本発明は、注目画素を含み、範囲が順次狭くなる多階層の領域毎に、その領域の画素の画素統計値を算出する画素統計値算出手段と、注目画素を含み、範囲が順次狭くなる多階層の領域毎に、エッジ情報を算出するエッジ情報算出手段と、補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値を用いて、当該階層の領域の画素統計値を補正し、前記エッジ情報を用いて、補正前の前記当該階層の領域の画素統計値と、補正後の前記当該階層の領域の画素統計値の差分情報とを再補正し、最大範囲から最小範囲の領域になるまで、順次各階層において、前記当該階層の領域の画素統計値の補正と再補正を繰り返すことで、前記注目画素を補正する補正手段とを有する画像処理装置である。

【発明の効果】

【0011】

本発明によれば、ノイズを効果的に除去することができる。

【図面の簡単な説明】

【0012】

【図1】図１は本発明の実施の形態である画像処理方法を説明するための図である。

【図2】図２はFunc関数の一例を示した図である。

【図3】図３は広領域におけるパラメータaの設定例を示す図である。

【図4】図４は本実施の形態における再補正関数の形状の一例を説明する為の図である。

【図5】図５は本実施の形態における再補正関数の形状の一例を説明する為の図である。

【図6】図６は本実施の形態における再補正関数の形状の一例を説明する為の図である。

【図7】図７は本実施の形態の画像処理装置のブロック図である。

【図8】図８は本実施の形態の画像処理装置のブロック図である。

【発明を実施するための形態】

【0013】

本発明の実施形態を説明する。

【0014】

図１は、本発明の実施の形態である画像処理方法を説明するための図である。尚、図１では、三階層の多重解像度処理を行う場合の処理の流れを示しているが、本発明は二階層にも対応できるし、四階層以上でも容易に拡張が可能である。

【0015】

図１は、本発明の実施の形態である画像処理方法を説明するための図である。

【0016】

本発明の実施の形態は、画素位置(x,y)（注目画素）を中心とする広領域における空間の画素統計値である空間平均値L3(x,y)、およびその領域でのエッジ情報、もしくはエッジ量E3(x,y)を用いて、中領域における空間平均値（L2(x,y)）を補正する。さらに、L2(x,y)と、補正された中領域における空間平均値（L2’(x,y)）の差分値R2(x,y)から、E3(x,y)に基づいて信号を抽出し、抽出した信号をL2’(x,y)に合成することで、L2(x,y)を再補正する。そして、再補正された中領域における空間平均値（L2”(x,y)）、中領域でのエッジ量（E2(x,y)）を用いて、狭領域の空間平均値（L1(x,y)）を補正する。この補正を順番に処理していくことで、入力画素値P_in(x,y)を補正し、出力画素値P_out(x,y)を得ることを特徴とする。

【0017】

ここで、画素統計値とは、対象領域の空間平均値であり、空間平均値は領域に存在する画素の相加平均値や、相乗平均値、加重平均値などである。以下の説明では、画素統計値を空間平均値とし、空間平均値を領域に存在する画素の相加平均値である場合を説明する。また、エッジ量、もしくはエッジ情報とは、注目画素（入力画素）を中心とした上下左右の領域間の画素の統計量（平均値、メディアンなど）の差分値で定義されるものとする。以下の説明では、画素統計値と同様に、エッジ量を空間平均値の差分である場合を説明する。

【0018】

以下、具体的な処理の説明を行うが、図１における各階層での処理は、補正量を決定するためのパラメータが異なるのみで、処理の流れは同一である。そこで，一例として，広領域における空間平均値L3(x,y)、およびその領域でのエッジ量E3(x,y)を用いて、中領域における空間平均値L2(x,y)を補正する処理の詳細を説明する。

【0019】

まず、画素位置(x,y)における広領域（範囲：-k3〜k3）の空間平均値L3(x,y)および中領域（範囲：-k2〜k2）の空間平均値L2(x,y)を式（１）および式（２）のように計算する。ここで、広領域、中領域をそれぞれｋ3、ｋ2にて縦横同じ画素数で範囲を指定しているが、必ずしも縦横同じ画素数でなくてもよい。

次に、広領域におけるエッジ量E3(x,y)を算出する。エッジ量の算出では，まず、縦方向のエッジ量EV3(x,y)と横方向のエッジ量EH3(x,y)とを、式（３）および式（４）のように算出し、これらを式（５）のように加算することで、広領域におけるエッジ量E3(x,y)を算出する。ここで、広領域をｋ3にて縦横同じ画素数で範囲を指定しているが、必ずしも縦横同じ画素数でなくてもよい。

続いて、算出されたエッジ量E3(x,y)から算出される合成重みα3(x,y)を用いて、式（６）のように広領域における空間平均値L3(x,y)を補正し、補正された広領域の空間平均値L3’(x,y)を算出する。合成重みα3(x,y)は、あらかじめ設定されている閾値hi3およびlo3を用いて、式（７）のように算出される。

最後に、算出された空間平均値L3’を用い、式（８）のように中領域の空間平均値L2(x,y)を補正する。

補正関数Funcの一例として、図２に示したものを用いる。例えば、画素位置(x,y)における中領域の空間平均値L2(x,y)の補正は、diffinを(L2(x,y)-L3”(x,y))とし、図２の補正関数によって得られる補正量diffoutをL2(x,y)に加算することで、実行される。図２の補正関数内のパラメータであるa,b,limitは、処理する解像度および補正する色成分毎に決定される。

【0020】

また、各階層において式（５）で算出されるエッジ量E3(x,y)を、ノイズ成分を抑制するFunc関数（補正関数）に反映させ、各階層のFunc関数（補正関数）を変更することで、各階層のノイズ成分を適応的に抑制する構成もできる。

【0021】

いま、エッジ量E3(x，y)に応じて式（９）に示すように値が変化する係数β3(x，y)を定義する。E3(x, y)の閾値であるhi3およびlo3は予め設定する値である。

式（９）で定義される係数β3(x，y)は、0〜1.0までの実数となる。係数β3(x，y)を利用して、以下の式（１０）のように、Func関数におけるパラメータaを設定する。エッジ量E3(x,y)とパラメータaの関係を図３に示す。

ここで、a_lo3は、エッジ量E3(x,y)が閾値lo3より小さい場合にパラメータaとして使用する値であり、a_hi3はエッジ量が閾値hi3より大きくなる場合のパラメータaに使用する値である。エッジ量E3(x,y)が閾値lo3からhi3までパラメータaは、a_hi3からa_lo3までの間の値となる。ここで、a_hi3は０以上の実数であり、a_lo3は，a_lo3>＝a_hi3の実数である。

【0022】

こうして得られたFunc関数を用い、中領域の空間平均値L2(x,y)を式（１１）のように補正する。

このようにして、式（１）から式（８）までを用いて、あるいは、式（１）から式（５）、および式（９）から式（１１）までを用いて、補正された中領域における空間平均値L2’(x,y)を得る。

【0023】

続いて、式（１２）に示すように、中領域における空間平均値L2(x,y)と、補正された中領域における空間平均値（L2’(x,y)）の差分値R2(x,y)を求める。すなわち、R2(x,y)はL2(x,y)に対する補正量を表す。さらに、式（１３）に示すように、R2(x,y)から、関数Fresidに基づいて信号を抽出し、抽出した信号をL2’(x,y)に合成して、再補正された中領域における空間平均値L2”(x,y)を得る。

式（１３）における関数F_residの一例を以下の式（１４）に示す。

α1<=α2であり、基本的には0.0<=α1<=α2<=1.0である。つまり、式（１４）は、注目画素がエッジに含まれるような場合（E3(x,y)>hi3）は、差分値R2(x,y)から高い割合で値を抽出することにより、最終的な補正量を小さくしてエッジの喪失を防ぐ。また、注目画素が平坦領域に含まれるような場合（E3(x,y)<lo3）は、差分値R2(x,y)から低い割合で値を抽出することにより、最終的な補正量を大きく保って平坦領域でのノイズ除去性能を維持する。さらに、注目画素がテクスチャ領域に含まれる場合は、差分値R2(x,y)から対象領域のエッジ量に比例して値を抽出することで、解像度感を保ちつつノイズを除去する。式（１４）において、関数F_residの入力r=1の場合の、エッジ量eと関数F_residの出力の関係を図４に示す。

【0024】

なお、式（１３）におけるα１及びα２は、各階層で変更してもよいし、そうでなくてもよい。

【0025】

また、関数F_residは式（１４）に限定されるものでなく、例えば、以下の式（１５）や式（１６）のような構成も可能である。

式（１５）、式（１６）において、関数F_residの入力r=1の場合の、エッジ量eと関数F_residの出力の関係を図５および図６に示す。

【0026】

更に、式（１３）を下記式（１７）に置き換えることによっても、同様な効果を得ることができる。

ここで、関数F_residは、式（１８）又は式（１９）となる。式（１８）におけるsign関数は入力の符号を出力する関数であり、式（１８）や式（１９）におけるλ’は式（２０）となる。

ここで、式（２０）のλはあらかじめ設定する値である。また、閾値hi、lo、パラメータα１、α２は式（１４）、（１５）、（１６）に示されたものと同様に、各階層で変更してもよいし、そうでなくてもよい。

【0027】

λ’についても式（２０）に限定されるものではなく、式（１４）に対する式（１５）および式（１６）と同様に代替も可能である。

【0028】

このように、第１の実施形態では、エッジ情報に基づいて、ノイズ抑制における再補正量の制御を行っている。実施の形態で述べた画像処理方法を用いることにより、ノイズの目立ちやすい平坦領域におけるノイズ除去性能を維持しつつ、エッジ領域・テクスチャ領域の解像度感を改善することができる。

【0029】

次に、本発明の実施の形態の画像処理装置を説明する。

【0030】

図７は本発明の実施の形態の画像処理装置のブロック図である。

【0031】

第１の実施形態は、式（１）から式（８）までを用いて補正値を得る方法と、式（１）から式（５）、および式（９）から式（１１）までを用いて補正値を得る方法がある。

【0032】

図７は式（１）から式（８）までを用いて補正値を得る場合の画像処理装置である。

【0033】

図７に示す本実施の形態の画像処理装置は、領域画素値抽出部１と、空間平均値算出部２と、補正部３と、再補正部４と、出力画像制御部５と、エッジ情報算出部６とを備える。

【0034】

領域画素値抽出部１は、画素位置(x,y)（注目画素）を中心とする広領域にある画素の画素値と、中領域にある画素の画素値と、狭領域にある画素の画素値と、入力画素値P_in(x,y) （注目画素）の画素値とを、それぞれのタイミングで抽出して空間平均値算出部２とエッジ情報算出部６に出力する。

【0035】

空間平均値算出部２は、領域画素値抽出部１からの各領域の画素値を受信し、その領域の空間平均値を算出する。そして、算出された広領域、中領域、狭領域の空間平均値は、補正部３と再補正部４に出力される。

【0036】

エッジ情報算出部６は、領域画素値抽出部１からの広領域に存在する画素の画素値に基づいて、広領域におけるエッジ量E3(x,y)を算出する。エッジ量の算出は、縦方向のエッジ量EV3(x,y)と横方向のエッジ量EH3(x,y)とを、式（３）および式（４）のように算出し、これらを式（５）のように加算することで、広領域におけるエッジ量E3(x,y)を算出する。同様に、中領域のエッジ量E2(x,y)、狭領域のエッジ量E1(x,y)を算出する。なお、横方向と縦方向のエッジ量について言及してあるが、斜め方向のエッジ量を算出し使用してもよい。

【0037】

補正部３は、エッジ情報算出部６により算出されたエッジ量E3 (x,y)から得られる合成重みα3 (x,y)を用いて、式（６）のように広領域における空間平均値L3(x,y)を補正し、補正された広領域の空間平均値L3”(x,y)を算出する。尚、合成重みα3 (x,y)は、あらかじめ設定されている閾値hi3およびlo3を用いて、式（７）のように算出される。

【0038】

更に、算出された空間平均値L3”を用い、式（８）のように中領域の空間平均値L2(x,y)を補正し、補正された空間平均値L2’(x,y)を得る。同様な補正を、空間平均値L1(x,y)、入力画素値P_in(x,y)にも行う。

【0039】

再補正部４は、補正前の空間平均値L2(x,y)と補正された空間平均値L2’(x,y)の差分値を計算し、エッジ情報算出部６により算出されたエッジ量E3 (x,y)に基づいて、式（１７）のようにL2’(x,y)を再補正し、再補正された空間平均値L2”(x,y)を得る。同様な再補正を、補正された空間平均値L1’(x,y)、補正された入力画素値P_in’(x,y)にも行う。

【0040】

出力画像制御部５は、順次補正された空間平均値が入力される毎に、次階層の領域の画素の画素値を抽出するように領域画素値抽出部１に指示を出す。また、再補正された空間平均値が入力される毎に、補正部３にその値をフィードバックする。そして、１画素のP_out(x,y)が入力されると、出力画素値としてP_out(x,y)を出力する。

【0041】

次に、式（１）から式（５）、および式（９）から式（１１）までを用いて補正値を得る場合の本発明の実施の形態の画像処理装置を図８に示す。

【0042】

図８に示す本発明の実施の形態の画像処理装置は、領域画素値抽出部１と、空間平均値算出部２と、補正部３と、再補正部４と、出力画像制御部５と、エッジ情報算出部６と、補正関数決定部７とを備える。

【0043】

領域画素値抽出部１、空間平均値算出部２、再補正部４、出力画像制御部５、エッジ情報算出部６は、図７に示す第１の実施形態の画像処理装置と同じ動作を行う。

【0044】

補正関数決定部７は、エッジ情報算出部６により算出されたエッジ量に基づいて、式（９）、（１０）に示すようにFunc関数（補正関数）におけるパラメータaを算出し、各階層のFunc関数（補正関数）を決定する。

【0045】

補正部３は、補正関数決定部７で決定された各階層のFunc関数（補正関数）を用いて各階層の空間平均値を補正する。

【0046】

以上の如く、本発明の実施の形態によれば、エッジ情報に基づいて、ノイズ抑制における再補正量の制御を行うことで、ノイズの目立ちやすい平坦領域におけるノイズ除去性能を維持しつつ、エッジ領域・テクスチャ領域の解像度感を改善することができる。

【0047】

尚、上述した説明からも明らかなように、各部をハードウェアで構成することも可能であるが、コンピュータプログラムにより実現することも可能である。この場合、プログラムメモリに格納されているプログラムで動作するプロセッサによって、上述した各実施の形態と同様の機能、動作を実現させる。また、上述した実施の形態の一部の機能のみをコンピュータプログラムにより実現することも可能である。

【0048】

また、上記の実施形態の一部又は全部は、以下の付記のようにも記載されうるが、以下には限られない。

【0049】

（付記１）
注目画素を含み、範囲が順次狭くなる多階層の領域毎に画素の画素統計値とエッジ情報とを求め、
当該階層の領域の画素統計値と当該階層の領域より広い階層の領域の画素統計値との差分情報を、前記エッジ情報を用いて補正し、
補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値を用いて、当該階層の領域の画素統計値を補正し、
前記エッジ情報を用いて、補正前の前記当該階層の領域の画素統計値と、補正後の前記当該階層の領域の画素統計値の差分情報とを再補正し、
最大範囲から最小範囲の領域になるまで、順次各階層において、前記当該階層の領域の画素統計値の補正と再補正を繰り返すことで、前記注目画素を補正する
画像処理方法。

【0050】

（付記２）
前記再補正は、
前記注目画素がエッジに含まれる場合には補正量を小さくし、
前記注目画素が平坦領域に含まれる場合には補正量を大きくし、
前記注目画素がテクスチャ領域に含まれる場合にはエッジ情報に応じた補正量により、
補正を行う
付記１に記載の画像処理方法。

【0051】

（付記３）
当該階層の領域の画素統計値と、当該階層の領域より広い階層の領域の画素統計値と、当該階層の領域より広い階層の領域におけるエッジ情報とを用いて、前記差分情報を算出する
付記１又は付記２に記載の画像処理方法。

【0052】

（付記４）
当該階層の領域より広い階層の領域におけるエッジ情報が所定の閾値を超える場合、当該階層の画素統計値の補正を行わない
付記３に記載の画像処理方法。

【0053】

（付記５）
補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値とを用いて、当該階層の領域の画素統計値を補正する際、補正の強度を階層毎に変化させる
付記１から付記４のいずれかに記載の画像処理方法。

【0054】

（付記６）
補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値とを用いて、当該階層の領域の画素統計値を補正する際、当該階層の領域のノイズによる画素値の変動量に応じて補正の強度を変化させる
付記５に記載の画像処理方法。

【0055】

（付記７）
補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値とを用いて、当該階層の領域の画素統計値を補正する際、当該階層の領域より広い領域のエッジ量に応じて補正の強度を変化させる
付記５に記載の画像処理方法。

【0056】

（付記８）
補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値とを用いて、当該階層の領域の画素統計値を補正する際、前記差分情報がゼロに近づくほどゼロに近い値を出力し、前記差分情報が大きくなるほど補正を弱める補正を行う
付記１から付記７のいずれかに記載の画像処理方法。

【0057】

（付記９）
前記画素統計値は、画素の空間平均値を用いる
付記１から付記８のいずれかに記載の画像処理方法。

【0058】

（付記１０）
前記空間平均値は、画素の相加平均値、相乗平均値、加重平均値のいずれかである
付記９に記載の画像処理方法。

【0059】

（付記１１）
注目画素を含み、範囲が順次狭くなる多階層の領域毎に、その領域の画素の画素統計値を算出する画素統計値算出手段と、
注目画素を含み、範囲が順次狭くなる多階層の領域毎に、エッジ情報を算出するエッジ情報算出手段と、
補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値を用いて、当該階層の領域の画素統計値を補正し、前記エッジ情報を用いて、補正前の前記当該階層の領域の画素統計値と、補正後の前記当該階層の領域の画素統計値の差分情報とを再補正し、最大範囲から最小範囲の領域になるまで、順次各階層において、前記当該階層の領域の画素統計値の補正と再補正を繰り返すことで、前記注目画素を補正する補正手段と
を有する画像処理装置。

【0060】

（付記１２）
前記補正手段は、
前記注目画素がエッジに含まれる場合には補正量を小さくし、
前記注目画素が平坦領域に含まれる場合には補正量を大きくし、
前記注目画素がテクスチャ領域に含まれる場合にはエッジ情報に応じた補正量により、
補正を行う
付記１１に記載の画像処理装置。

【0061】

（付記１３）
前記補正手段は、当該階層の領域の画素統計値と、当該階層の領域より広い階層の領域の画素統計値と、当該階層の領域より広い階層の領域におけるエッジ情報とを用いて、前記差分情報を算出する
付記１１又は付記１２に記載の画像処理装置。

【0062】

（付記１４）
前記補正手段は、当該階層の領域より広い階層の領域におけるエッジ情報が所定の閾値を超える場合、当該階層の画素統計値の補正を行わない
付記１３に記載の画像処理装置。

【0063】

（付記１５）
前記補正手段は、補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値とを用いて、当該階層の領域の画素統計値を補正する際、補正の強度を階層毎に変化させる
付記１１から付記１４のいずれかに記載の画像処理装置。

【0064】

（付記１６）
前記補正手段は、補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値とを用いて、当該階層の領域の画素統計値を補正する際、当該階層の領域のノイズによる画素値の変動量に応じて補正の強度を変化させる
付記１５に記載の画像処理装置。

【0065】

（付記１７）
前記補正手段は、補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値とを用いて、当該階層の領域の画素統計値を補正する際、当該階層の領域より広い領域のエッジ量に応じて補正の強度を変化させる
付記１５に記載の画像処理装置。

【0066】

（付記１８）
前記補正手段は、補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値とを用いて、当該階層の領域の画素統計値を補正する際、前記差分情報がゼロに近づくほどゼロに近い値を出力し、前記差分情報が大きくなるほど補正を弱める補正を行う
付記１１から付記１７のいずれかに記載の画像処理装置。

【0067】

（付記１９）
前記画素統計値は、画素の空間平均値を用いる
付記１１から付記１８のいずれかに記載の画像処理装置。

【0068】

（付記２０）
前記空間平均値は、画素の相加平均値、相乗平均値、加重平均値のいずれかである
付記１９に記載の画像処理装置。

【0069】

（付記２１）
注目画素を含み、範囲が順次狭くなる多階層の領域毎に、その領域の画素の画素統計値を算出する画素統計値算出処理と、
注目画素を含み、範囲が順次狭くなる多階層の領域毎に、エッジ情報を算出するエッジ情報算出処理と、
補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値を用いて、当該階層の領域の画素統計値を補正し、前記エッジ情報を用いて、補正前の前記当該階層の領域の画素統計値と、補正後の前記当該階層の領域の画素統計値の差分情報とを再補正し、最大範囲から最小範囲の領域になるまで、順次各階層において、前記当該階層の領域の画素統計値の補正と再補正を繰り返すことで、前記注目画素を補正する補正処理と
を、コンピュータに実行させるプログラム。

【0070】

（付記２２）
前記補正処理は、
前記注目画素がエッジに含まれる場合には補正量を小さくし、
前記注目画素が平坦領域に含まれる場合には補正量を大きくし、
前記注目画素がテクスチャ領域に含まれる場合にはエッジ情報に応じた補正量により、
補正を行う
付記２１に記載のプログラム。

【0071】

（付記２３）
前記補正処理は、当該階層の領域の画素統計値と、当該階層の領域より広い階層の領域の画素統計値と、当該階層の領域より広い階層の領域におけるエッジ情報とを用いて、前記差分情報を算出する
付記２１又は付記２２に記載のプログラム。

【0072】

（付記２４）
前記補正処理は、当該階層の領域より広い階層の領域におけるエッジ情報が所定の閾値を超える場合、当該階層の画素統計値の補正を行わない
付記２３に記載のプログラム。

【0073】

（付記２５）
前記補正処理は、補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値とを用いて、当該階層の領域の画素統計値を補正する際、補正の強度を階層毎に変化させる
付記１１から付記１４のいずれかに記載のプログラム。

【0074】

（付記２６）
前記補正処理は、補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値とを用いて、当該階層の領域の画素統計値を補正する際、当該階層の領域のノイズによる画素値の変動量に応じて補正の強度を変化させる
付記２５に記載のプログラム。

【0075】

（付記２７）
前記補正処理は、補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値とを用いて、当該階層の領域の画素統計値を補正する際、当該階層の領域より広い領域のエッジ量に応じて補正の強度を変化させる
付記２５に記載のプログラム。

【0076】

（付記２８）
前記補正処理は、補正後の前記差分情報と、前記当該階層の領域より広い領域の画素統計値とを用いて、当該階層の領域の画素統計値を補正する際、前記差分情報がゼロに近づくほどゼロに近い値を出力し、前記差分情報が大きくなるほど補正を弱める補正を行う
付記２１から付記２７のいずれかに記載のプログラム。

【0077】

（付記２９）
前記画素統計値は、画素の空間平均値を用いる
付記２１から付記２８のいずれかに記載のプログラム。

【0078】

（付記３０）
前記空間平均値は、画素の相加平均値、相乗平均値、加重平均値のいずれかである
付記２９に記載のプログラム。

【0079】

以上好ましい実施の形態をあげて本発明を説明したが、本発明は必ずしも上記実施の形態に限定されるものではなく、その技術的思想の範囲内において様々に変形し実施することが出来る。

【0080】

本出願は、２０１２年４月２６日に出願された日本出願特願２０１２−１００９１０号を基礎とする優先権を主張し、その開示の全てをここに取り込む。

【符号の説明】

【0081】

１ 領域画素値抽出部
２ 空間平均値算出部
３ 補正部
４ 再補正部
５ 出力画像制御部
６ エッジ情報算出部
７ 補正関数決定部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】