(57)【特許請求の範囲】
【請求項１】
原画像を表す原画像データに対して画像処理を施す画像処理装置において、
前記原画像データを多重解像度変換して、該原画像データから複数の周波数帯域毎の画像を表す画像データを得、該複数の周波数帯域のうち所定周波数帯域における前記原画像のｘ方向高周波ｙ方向低周波画像データ、ｘ方向低周波ｙ方向高周波画像データおよび前記所定周波数帯域の次に低周波の周波数帯域におけるｘ方向高周波ｙ方向高周波画像データのそれぞれの画素値を、該ｘ方向高周波ｙ方向低周波画像データ、該ｘ方向低周波ｙ方向高周波画像データおよび前記所定周波数帯域の次に低周波の周波数帯域におけるｘ方向高周波ｙ方向高周波画像データのそれぞれに対応する第１から第３の閾値と比較し、前記ｘ方向高周波ｙ方向低周波画像データの画素値が前記第１の閾値未満となり、前記ｘ方向低周波ｙ方向高周波画像データの画素値が前記第２の閾値未満となり、かつ前記所定周波数帯域の次に低周波の周波数帯域におけるｘ方向高周波ｙ方向高周波画像データの画素値が前記第３の閾値未満となる画素を抽出し、前記複数の画像データのうち前記所定周波数帯域におけるｘ方向高周波ｙ方向高周波画像データの前記抽出された画素に対応する画素の画素値の分散値を算出し、該分散値を第４の閾値と比較することにより前記原画像の粒状が多いか少ないかを推定する粒状推定手段と、
前記粒状が少ない場合に鮮鋭度を重視した補間演算を、前記粒状が多い場合に滑らかさを重視した補間演算を前記原画像データに対して施して処理済み画像データを得る画像処理手段とを備えたことを特徴とする画像処理装置。
【請求項２】
原画像を表す原画像データに対して画像処理を施す画像処理装置において、
前記原画像データを多重解像度変換して、該原画像データから複数の周波数帯域毎の画像を表す画像データを得、該複数の周波数帯域のうち所定周波数帯域における前記原画像のｘ方向高周波ｙ方向低周波画像データ、ｘ方向低周波ｙ方向高周波画像データおよび前記所定周波数帯域の次に低周波の周波数帯域におけるｘ方向高周波ｙ方向高周波画像データのそれぞれの画素値を、該ｘ方向高周波ｙ方向低周波画像データ、該ｘ方向低周波ｙ方向高周波画像データおよび前記所定周波数帯域の次に低周波の周波数帯域におけるｘ方向高周波ｙ方向高周波画像データのそれぞれに対応する第１から第３の閾値と比較し、前記ｘ方向高周波ｙ方向低周波画像データの画素値が前記第１の閾値未満となり、前記ｘ方向低周波ｙ方向高周波画像データの画素値が前記第２の閾値未満となり、かつ前記所定周波数帯域の次に低周波の周波数帯域におけるｘ方向高周波ｙ方向高周波画像データの画素値が前記第３の閾値未満となる画素を抽出し、前記複数の画像データのうち前記所定周波数帯域におけるｘ方向高周波ｙ方向高周波画像データの前記抽出された画素に対応する画素の画素値の分散値を算出し、該分散値を第４の閾値と比較することにより前記原画像の粒状が多いか少ないかを推定する粒状推定手段と、
前記粒状が多い場合に粒状を低減する処理を前記原画像データに対して施して処理済み画像データを得る画像処理手段とを備えたことを特徴とする画像処理装置。
【請求項３】
原画像を表す原画像データに対して画像処理を施す画像処理装置において、
前記原画像データを多重解像度変換して、該原画像データから複数の周波数帯域毎の画像を表す画像データを得、前記複数の周波数帯域のうち所定周波数帯域の次に低周波の周波数帯域における前記原画像のｘ方向高周波ｙ方向低周波画像データおよびｘ方向低周波ｙ方向高周波画像データのそれぞれの画素値を、該ｘ方向高周波ｙ方向低周波画像データおよび該ｘ方向低周波ｙ方向高周波画像データのそれぞれに対応する第１および第２の閾値と比較し、前記所定周波数帯域の次に低周波の周波数帯域における前記ｘ方向高周波ｙ方向低周波画像データの画素値が前記第１の閾値を越え、かつ前記ｘ方向低周波ｙ方向高周波画像データの画素値が前記第２の閾値を越える画素を抽出し、前記所定周波数帯域の次に低周波の周波数帯域における前記ｘ方向高周波ｙ方向低周波画像データに対する前記所定周波数帯域における前記原画像のｘ方向高周波ｙ方向低周波画像データの前記抽出された画素の画素値の比および前記所定周波数帯域の次に低周波の周波数帯域における前記ｘ方向低周波ｙ方向高周波画像データに対する前記所定周波数帯域における前記原画像のｘ方向低周波ｙ方向高周波画像データの前記抽出された画素の画素値の比のいずれか大きい方の比を算出し、該比に応じて前記原画像の鮮鋭度が高いか低いかを推定する鮮鋭度推定手段と、
前記鮮鋭度が低い場合に鮮鋭度を重視した補間演算を、前記鮮鋭度が高い場合に滑らかさを重視した補間演算を前記原画像データに対して施して処理済み画像データを得る画像処理手段とを備えたことを特徴とする画像処理装置。
【請求項４】
原画像を表す原画像データに対して画像処理を施す画像処理装置において、
前記原画像データを多重解像度変換して、該原画像データから複数の周波数帯域毎の画像を表す画像データを得、前記複数の周波数帯域のうち所定周波数帯域の次に低周波の周波数帯域における前記原画像のｘ方向高周波ｙ方向低周波画像データおよびｘ方向低周波ｙ方向高周波画像データのそれぞれの画素値を、該ｘ方向高周波ｙ方向低周波画像データおよび該ｘ方向低周波ｙ方向高周波画像データのそれぞれに対応する第１および第２の閾値と比較し、前記所定周波数帯域の次に低周波の周波数帯域における前記ｘ方向高周波ｙ方向低周波画像データの画素値が前記第１の閾値を越え、かつ前記ｘ方向低周波ｙ方向高周波画像データの画素値が前記第２の閾値を越える画素を抽出し、前記所定周波数帯域の次に低周波の周波数帯域における前記ｘ方向高周波ｙ方向低周波画像データに対する前記所定周波数帯域における前記原画像のｘ方向高周波ｙ方向低周波画像データの前記抽出された画素の画素値の比および前記所定周波数帯域の次に低周波の周波数帯域における前記ｘ方向低周波ｙ方向高周波画像データに対する前記所定周波数帯域における前記原画像のｘ方向低周波ｙ方向高周波画像データの前記抽出された画素の画素値の比のいずれか大きい方の比を算出し、該比に応じて前記原画像の鮮鋭度を推定する鮮鋭度推定手段と、
前記鮮鋭度が低いほど鮮鋭度を強調する処理を前記原画像データに対して施して処理済み画像データを得る画像処理手段とを備えたことを特徴とする画像処理装置。
【請求項５】
複数の原画像を表す原画像データを合成して合成画像データを得る画像合成装置において、
請求項１から４のいずれか１項記載の画像処理装置により前記各原画像データに対して画像処理を施して処理済み画像データを得る手段と、
該各処理済み画像データを合成して前記合成画像データを得る合成手段とを備えたことを特徴とする画像合成装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、画像データに対して画像処理を施す画像処理装置および複数の画像を合成する画像合成装置に関するものである。
【０００２】
【従来の技術】
画像を表す画像データを得、この画像データに適切な画像処理を施した後、画像を再生することが種々の分野で行われている。このような画像を再生する出力装置は、システムとして接続された入力装置の特性と出力装置の特性とに応じて、画像データを補正するなどの画像処理を施して画像を再生している。ここで、出力装置に入力装置が１台のみしか接続されていない場合は、出力装置において行われる画像処理はその入力装置の特性に対応したもののみを行えばよく、出力装置に複数の入力装置が接続されている場合は、各入力装置の特性に応じた画像処理を画像データに対して施せばよい。
【０００３】
しかしながら、複数の入力装置と出力装置とを接続した場合、入力装置は種々変更される場合があり、また、出力装置に入力される画像データはいずれの入力装置において得られたものであるかが特定できないため、入力装置の特性と画像処理の内容とを対応させて最適な画像処理を行うことは困難である。このため、例えば特許第2660170号においては、入力装置から入力された画像データについては全て汎用的な階調となるように画像データを変換し、出力装置においては自らの出力特性のみを補正して適正な再生画像を得るようにした画像処理装置が提案されている。
【０００４】
また、入力装置において得られる画像は、その装置に応じて粒状特性（ノイズ成分の状態）や鮮鋭度も種々異なるものとなっている。このため、得られた画像データに対して入力装置の粒状特性および／または鮮鋭度特性に応じてスムージング処理、鮮鋭度強調処理などの画像処理を施して、より高画質の画像を再生することが行われている。また、画像データを補間することにより画像を拡大、縮小して再生する場合もあり、この場合も入力装置の粒状特性、鮮鋭度特性に応じて、鮮鋭度を重視した補間演算（例えばCubicスプライン補間演算）、滑らかさを重視した補間演算（例えばＢスプライン補間演算）を施している。
【０００５】
一方、上述した画像データを処理するための方法として、画像を複数の周波数帯域毎の多重解像度画像に変換し、各周波数帯域の画像に対して所定の処理を行って、再度これを逆多重解像度変換することにより、最終的な処理済み画像を得るための多重解像度変換なる方法が提案されている。この場合の所定の処理としては、具体的には、ノイズ除去のための高周波の分離、さらにはノイズの多い周波数帯域のデータを削減することによる圧縮処理などが挙げられる。この多重解像度変換の方法としてはウェーブレット変換、ラプラシアンピラミッドあるいはフーリエ変換等の方法が知られている。とくに、ウェーブレット変換は、信号の周波数解析方法の１つであるが、同じく周波数解析方法として広く用いられているフーリエ変換に比べ、信号の局所的な変化情報を検出しやすいという点で優れていることから、近年あらゆる信号処理の分野で脚光を浴びている（OLIVIER RIOUL and MARTIN VETTERLI;Wavelets and Signal Processing,IEEE SP MAGAZINE,P.14-38,OCTOBER 1991、Stephane Mallat;Zero-Crossings of a   Wavelet Transform,IEEE TRANSACTIONS ON INFORMATION THEORY,VOL.37,NO.4,P.1019-1033,JULY 1991 、特開平6-350989号、同6-350990号、同7-23228号、同7-23229号、同7-79350号、特願平8-14510号など）。
【０００６】
【発明が解決しようとする課題】
ところで、近年のコンピュータネットワークの進歩により、画像の出力装置には種々の入力装置が接続され、さらにネットワーク上を様々な処理が施された画像データが転送されている。このため、出力装置に入力される画像データはいずれの入力装置において得られたものであるか、さらにその画像データの粒状あるいは鮮鋭度、階調について画像処理が施されているのか否かが全く分からない状態でデータがやり取りされている。このような場合、上記特許第2660170号に記載された画像処理装置を用いれば、入力装置の種類に拘わらず汎用的な階調となるように画像データに対して画像処理を施すことができる。
【０００７】
しかしながら、画像データにより表される画像の粒状（ノイズ成分）あるいは鮮鋭度については、画像データがいずれの入力装置により得られたものであるかが不明であると、その入力装置の特性に最適なスムージング処理、鮮鋭度強調処理、補間演算などの画像処理を行うことができない。この場合、適当に定めた１つの手法により全ての画像データに対して同一の粒状、鮮鋭度に関する画像処理を施すことが考えられるが、画像によってはその画像処理が最適ではなく、例えば粒状が多いにも拘わらず鮮鋭度が強調されてしまうおそれがある。また、鮮鋭度強調の程度が不足したぼけた画像となったり、元々鮮鋭度が高いにも拘わらずさらに鮮鋭度が強調されてオーバーシュート、アンダーシュートなどのアーチファクトが発生した画像となるおそれがある。とくに、複数の画像を合成して再生する場合にこの手法を用いると、合成に用いる画像データがそれぞれ異なる粒状特性、鮮鋭度特性を有するものであっても同一の画像処理が施されてしまうため、合成された画像は、部分毎に粒状や鮮鋭度のバランスが異なって見た目に違和感のあるものとなってしまう。また、一旦画像を再生し、オペレータがその画像を観察しながら人為的に画像データを修正することも考えられるが、非常に手間がかかり、オペレータの負担が大きい。
【０００８】
本発明は上記事情に鑑みなされたものであり、画像の粒状および／または鮮鋭度に拘わらずその画像に最適な画像処理を施すことができる画像処理装置および画像合成装置を提供することを目的とするものである。
【０００９】
【課題を解決するための手段】
本発明による第１の画像処理方法は、原画像を表す原画像データに対して画像処理を施す画像処理装置において、
前記原画像データから前記原画像の粒状が多いか少ないかを推定する推定手段と、
前記粒状が少ない場合に鮮鋭度を重視した補間演算を、前記粒状が多い場合に滑らかさを重視した補間演算を前記原画像データに対して施して処理済み画像データを得る画像処理手段とを備えたことを特徴とするものである。
【００１０】
本発明による第２の画像処理装置は、原画像を表す原画像データに対して画像処理を施す画像処理装置において、
前記原画像データから前記原画像の粒状が多いか少ないかを推定する推定手段と、
前記粒状が多い場合に粒状を低減する処理を前記原画像データに対して施して処理済み画像データを得る画像処理手段とを備えたことを特徴とするものである。
【００１１】
本発明による第３の画像処理装置は、原画像を表す原画像データに対して画像処理を施す画像処理装置において、
前記原画像データを多重解像度変換して、該原画像データから複数の周波数帯域毎の画像を表す画像データを得、該複数の周波数帯域のうち所定周波数帯域における前記原画像のｘ方向高周波ｙ方向低周波画像データ、ｘ方向低周波ｙ方向高周波画像データおよび前記所定周波数帯域の次に低周波の周波数帯域におけるｘ方向高周波ｙ方向高周波画像データのそれぞれの画素値を、該ｘ方向高周波ｙ方向低周波画像データ、該ｘ方向低周波ｙ方向高周波画像データおよび前記所定周波数帯域の次に低周波の周波数帯域におけるｘ方向高周波ｙ方向高周波画像データのそれぞれに対応する第１から第３の閾値と比較し、前記ｘ方向高周波ｙ方向低周波画像データの画素値が前記第１の閾値未満となり、前記ｘ方向低周波ｙ方向高周波画像データの画素値が前記第２の閾値未満となり、かつ前記所定周波数帯域の次に低周波の周波数帯域におけるｘ方向高周波ｙ方向高周波画像データの画素値が前記第３の閾値未満となる画素を抽出し、前記複数の画像データのうち前記所定周波数帯域におけるｘ方向高周波ｙ方向高周波画像データの前記抽出された画素に対応する画素の画素値の分散値を算出し、該分散値を第４の閾値と比較することにより前記原画像の粒状が多いか少ないかを推定する粒状推定手段と、
前記粒状が少ない場合に鮮鋭度を重視した補間演算を、前記粒状が多い場合に滑らかさを重視した補間演算を前記原画像データに対して施して処理済み画像データを得る画像処理手段とを備えたことを特徴とするものである。
【００１２】
本発明による第４の画像処理装置は、原画像を表す原画像データに対して画像処理を施す画像処理装置において、
前記原画像データを多重解像度変換して、該原画像データから複数の周波数帯域毎の画像を表す画像データを得、該複数の周波数帯域のうち所定周波数帯域における前記原画像のｘ方向高周波ｙ方向低周波画像データ、ｘ方向低周波ｙ方向高周波画像データおよび前記所定周波数帯域の次に低周波の周波数帯域におけるｘ方向高周波ｙ方向高周波画像データのそれぞれの画素値を、該ｘ方向高周波ｙ方向低周波画像データ、該ｘ方向低周波ｙ方向高周波画像データおよび前記所定周波数帯域の次に低周波の周波数帯域におけるｘ方向高周波ｙ方向高周波画像データのそれぞれに対応する第１から第３の閾値と比較し、前記ｘ方向高周波ｙ方向低周波画像データの画素値が前記第１の閾値未満となり、前記ｘ方向低周波ｙ方向高周波画像データの画素値が前記第２の閾値未満となり、かつ前記所定周波数帯域の次に低周波の周波数帯域におけるｘ方向高周波ｙ方向高周波画像データの画素値が前記第３の閾値未満となる画素を抽出し、前記複数の画像データのうち前記所定周波数帯域におけるｘ方向高周波ｙ方向高周波画像データの前記抽出された画素に対応する画素の画素値の分散値を算出し、該分散値を第４の閾値と比較することにより前記原画像の粒状が多いか少ないかを推定する粒状推定手段と、
前記粒状が多い場合に粒状を低減する処理を前記原画像データに対して施して処理済み画像データを得る画像処理手段とを備えたことを特徴とするものである。
【００１３】
本発明による第５の画像処理装置は、原画像を表す原画像データに対して画像処理を施す画像処理装置において、
前記原画像データを多重解像度変換して、該原画像データから複数の周波数帯域毎の画像を表す画像データを得、前記複数の周波数帯域のうち所定周波数帯域の次に低周波の周波数帯域における前記原画像のｘ方向高周波ｙ方向低周波画像データおよびｘ方向低周波ｙ方向高周波画像データのそれぞれの画素値を、該ｘ方向高周波ｙ方向低周波画像データおよび該ｘ方向低周波ｙ方向高周波画像データのそれぞれに対応する第１および第２の閾値と比較し、前記所定周波数帯域の次に低周波の周波数帯域における前記ｘ方向高周波ｙ方向低周波画像データの画素値が前記第１の閾値を越え、かつ前記ｘ方向低周波ｙ方向高周波画像データの画素値が前記第２の閾値を越える画素を抽出し、前記所定周波数帯域の次に低周波の周波数帯域における前記ｘ方向高周波ｙ方向低周波画像データに対する前記所定周波数帯域における前記原画像のｘ方向高周波ｙ方向低周波画像データの前記抽出された画素の画素値の比および前記所定周波数帯域の次に低周波の周波数帯域における前記ｘ方向低周波ｙ方向高周波画像データに対する前記所定周波数帯域における前記原画像のｘ方向低周波ｙ方向高周波画像データの前記抽出された画素の画素値の比のいずれか大きい方の比を算出し、該比に応じて前記原画像の鮮鋭度が高いか低いかを推定する鮮鋭度推定手段と、
前記鮮鋭度が低い場合に鮮鋭度を重視した補間演算を、前記鮮鋭度が高い場合に滑らかさを重視した補間演算を前記原画像データに対して施して処理済み画像データを得る画像処理手段とを備えたことを特徴とするものである。
【００１４】
本発明による第６の画像処理装置は、原画像を表す原画像データに対して画像処理を施す画像処理装置において、
前記原画像データを多重解像度変換して、該原画像データから複数の周波数帯域毎の画像を表す画像データを得、前記複数の周波数帯域のうち所定周波数帯域の次に低周波の周波数帯域における前記原画像のｘ方向高周波ｙ方向低周波画像データおよびｘ方向低周波ｙ方向高周波画像データのそれぞれの画素値を、該ｘ方向高周波ｙ方向低周波画像データおよび該ｘ方向低周波ｙ方向高周波画像データのそれぞれに対応する第１および第２の閾値と比較し、前記所定周波数帯域の次に低周波の周波数帯域における前記ｘ方向高周波ｙ方向低周波画像データの画素値が前記第１の閾値を越え、かつ前記ｘ方向低周波ｙ方向高周波画像データの画素値が前記第２の閾値を越える画素を抽出し、前記所定周波数帯域の次に低周波の周波数帯域における前記ｘ方向高周波ｙ方向低周波画像データに対する前記所定周波数帯域における前記原画像のｘ方向高周波ｙ方向低周波画像データの前記抽出された画素の画素値の比および前記所定周波数帯域の次に低周波の周波数帯域における前記ｘ方向低周波ｙ方向高周波画像データに対する前記所定周波数帯域における前記原画像のｘ方向低周波ｙ方向高周波画像データの前記抽出された画素の画素値の比のいずれか大きい方の比を算出し、該比に応じて前記原画像の鮮鋭度を推定する鮮鋭度推定手段と、
前記鮮鋭度が低いほど鮮鋭度を強調する処理を前記原画像データに対して施して処理済み画像データを得る画像処理手段とを備えたことを特徴とするものである。
【００１５】
ここで、「鮮鋭度を重視した補間演算」としては、例えばCubicスプライン補間演算を、「滑らかさ重視の補間演算」としては、例えばＢスプライン補間演算を用いることができる。
【００１６】
本発明による画像合成装置は、複数の原画像を表す原画像データを合成して合成画像データを得る画像合成装置において、
本発明による第１から第６の画像処理装置により前記各原画像データに対して画像処理を施して処理済み画像データを得る手段と、
該各処理済み画像データを合成して前記合成画像データを得る合成手段とを備えたことを特徴とするものである。
【００１７】
なお、本発明による画像処理装置および画像合成装置において行われる処理を画像処理方法および画像合成方法として提供してもよい。
【００１８】
【発明の効果】
本発明による画像処理装置は、原画像データに基づいて原画像の粒状および／または鮮鋭度を推定し、推定された粒状および／または鮮鋭度に応じて、原画像データに対して異なる画像処理を施すようにしたものである。このため、粒状および／または鮮鋭度の推定結果を用いて、原画像の粒状量、鮮鋭度に応じたスムージング処理、鮮鋭度強調処理、補間演算などを原画像データに対して施すことができる。例えば、粒状が多いと推定された場合は、原画像データに対して粒状を低減するためのスムージング処理や、補間を行う場合は滑らかさを重視したＢスプライン補間演算を施し、粒状が少ないと推定された場合は、原画像データに対してスムージング処理を施さないか、補間を行う場合は鮮鋭度を重視したCubicスプライン補間演算を施すことができる。一方、原画像における鮮鋭度が高いと推定された部分にはそれほど鮮鋭度を強調せず、鮮鋭度が低いと推定された部分には鮮鋭度をより強調するように鮮鋭度強調処理を施すことができる。また、画像全体として鮮鋭度が高いときに原画像を拡大、縮小する場合は滑らかさを重視したＢスプライン補間演算を行い、鮮鋭度が低いときには鮮鋭度を重視したCubicスプライン補間演算を行うことができる。そしてこれにより、その原画像がいずれの入力装置から入力されたものであるか、すでに画像処理が施されものであるかが不明な場合であっても、その原画像に最適な粒状および／または鮮鋭度に関する画像処理を施すことができる。
【００１９】
また、本発明の画像合成方法および装置は、複数の原画像データに対して本発明の画像処理方法および装置により画像処理を施して処理済み画像データを得、この処理済み画像データを合成して合成画像データを得るようにしたものである。したがって、原画像の粒状および／または鮮鋭度に応じて最適な画像処理が施された処理済み画像データに基づいて合成画像データを得ることができるため、合成画像データを再生することにより、全ての原画像データに対して同一の画像処理を施した場合に生じる、部分毎の粒状バランス、鮮鋭度バランスの崩れをなくすことができ、これにより見た目に違和感のない画像を得ることができる。
【００２０】
【発明の実施の形態】
以下図面を参照して本発明の実施形態について説明する。
【００２１】
図１は本発明の実施形態による画像処理装置の構成を示す概略ブロック図である。図１に示すように本実施形態による画像処理装置１は、画像入力手段２から入力された原画像データＳにより表される画像の粒状量を推定する粒状推定手段３と、画像入力手段２から入力された原画像データＳにより表される画像の鮮鋭度を推定する鮮鋭度推定手段７と、粒状推定手段３および鮮鋭度推定手段７において推定された粒状量および／または鮮鋭度に基づいて、画像処理手段５における画像処理方法の選択および／または画像処理のパラメータの設定を行う画像処理方法選択および／またはパラメータ設定手段４と、画像処理方法選択および／またはパラメータ設定手段４において設定された画像処理方法および／またはパラメータにしたがって、原画像データＳに対して画像処理を施して処理済み画像データＳ′を得る画像処理手段５と、処理済み画像データＳ′をプリンタ、ＣＲＴあるいはネットワーク上に出力する画像出力手段６とを備えてなる。
【００２２】
図２は粒状推定手段３の構成を示す概略ブロック図である。図２に示すように粒状推定手段３は、原画像データＳに対してウェーブレット変換を施して、原画像データＳにより表される原画像を多重解像度変換し、複数の周波数帯域毎の画像を表す画像データを得るウェーブレット変換手段１０と、ウェーブレット変換手段１０において得られた画像データＨＬ０，ＬＨ０，ＨＨ１（詳細は後述する）をそれぞれ所定の閾値と比較する比較手段１１と、比較手段１１による比較結果に基づいて、画像データＨＨ０の各画素の分散値Ｂを算出する分散値算出手段１２と、算出された分散値Ｂに基づいて原画像データＳの粒状量を推定する推定手段１３とを備える。
【００２３】
ウェーブレット変換手段１０においては、以下のようにして原画像データＳに対してウェーブレット変換が施される。図３は原画像データＳをウェーブレット変換する状態を模式的に示す図である。なお、本実施形態においては、図３に示すように紙面に対する横方向をｘ方向、縦方向をｙ方向として説明する。本実施形態のウェーブレット変換においては、下記に示すHaar-Wavelet基底を使用するとともに、ウェーブレット変換において通常行われている画素の間引き（サブサンプリング）を実施せず、ウェーブレット変換により得られる多重解像度画像の解像度を原画像と同一とする。
【００２４】
Haar-Wavelet：h(x)=-1,1
g(x)=1/2,1/2
ここで、h(x)はｘ方向について高周波成分を抽出するウェーブレット基底であり、g(x)はｘ方向について低周波成分を抽出するウェーブレット基底である。
【００２５】
まず、原画像データＳに対してｘ方向にh(x)を畳み込み演算し、画像データＨ０を得る。
【００２６】
Ｈ０(x,y)＝h(x)＊Ｓ(x,y)（(x,y)は画素の位置）
なお、通常のウェーブレット変換においては、ここで画像データＨ０のｘ方向について画素数を１／２に間引くが、本実施形態においては間引かないものとする。
【００２７】
次に画像データＨ０に対してｙ方向にh(y)を畳み込み演算してｘ方向ｙ方向ともに最高周波数帯域の画像を表す画像データＨＨ０を得る。
【００２８】
ＨＨ０(x,y)＝h(y)＊Ｈ０(x,y)
以下、同様にして下記の演算により、画像データＨ０にg(y)を畳み込み演算した画像データＨＬ０、原画像データＳにg(x)を畳み込み演算した画像データＬ０にh(y)を畳み込み演算した画像データＬＨ０、および画像データＬ０にg(y)を畳み込み演算した画像データＬＬ１を求める。
【００２９】
ＨＬ０(x,y)＝g(y)＊Ｈ０(x,y)
Ｌ０(x,y)＝g(x)＊Ｓ(x,y)
ＬＨ０(x,y)＝h(y)＊Ｌ０(x,y)
ＬＬ１(x,y)＝g(y)＊Ｌ０(x,y)
さらに、次の階層のHaar-Wavelet基底であるh(x)′=-1,0,1を使用し、画像データＬＬ１に、h(x)′およびh(y)′を畳み込み演算して画像データＨＨ１を求める。ここで、画像データＬＬ１を求める際に画素を間引いていれば、ウェーブレット基底として上述したh(x)=-1,1を用いればよいが、本実施形態においては画素を間引いていないため、間引きに相当する数の値０を挿入したものである。
【００３０】
ＨＨ１(x,y)＝h(y)′＊｛h(x)′＊ＬＬ１(x,y)｝
なお、この後ウェーブレット変換を繰り返すことにより、さらに低周波数帯域の画像データを求めることができるが、粒状推定手段３においては画像データＨＨ０，ＨＬ０，ＬＨ０，ＨＨ１のみを用いるため、ここで演算を終了する。なお、本実施形態においては、ウェーブレット変換を行う際に画素を間引いていないため、画像データＨＨ０，ＨＬ０，ＬＨ０，ＬＬ１，ＨＨ１により表される画像は図３に示すように原画像と同一の解像度を有するものとなる。また、画像データＨＨ０、画像データＨＬ０、画像データＬＨ０および画像データＬＬ１は、それぞれ原画像の最高周波数帯域におけるｘ方向ｙ方向ともに高周波画像、ｘ方向高周波ｙ方向低周波画像、ｘ方向低周波ｙ方向高周波画像およびｘ方向ｙ方向ともに低周波画像を表すものとなる。
【００３１】
比較手段１１においては、画像データＨＬ０，ＬＨ０，ＨＨ１の画素値が閾値Ｔｈ１，Ｔｈ２，Ｔｈ３と比較され、ＨＬ０(x,y)＜Ｔｈ１かつＬＨ０(x,y)＜Ｔｈ２かつＨＨ１(x,y)＜Ｔｈ３の３条件を満たす画素(x,y)が抽出される。
【００３２】
分散値算出手段１２は、比較手段１１において抽出された画素(x,y)に対応する画像データＨＨ０の画素のみを用いて下記の式により画素値の分散値Ｂを算出する。なお、本実施形態においては、図３に示すように各画像データＨＨ０，ＨＬ０，ＬＨ０，ＬＬ０，ＨＨ１の解像度は同一であるため、画素(x,y)の対応付けを正確かつ簡易に行うことができる。
【００３３】
Ｂ＝Σ｛ＨＨ０(x,y)｝2／Ｍ（Ｍ：抽出された画素数）
ここで、原画像に含まれるエッジなどの境界線は、画素値として高い値を有するが、最高周波数帯域画像を表す画像データＨＨ０および低周波数帯域画像を表す画像データＨＬ０，ＬＨ０，ＨＨ１のいずれにおいても高い値を有する。一方、原画像に含まれる粒状（ノイズ成分）も境界線と同様に画素値として高い値を有するため、画像データＨＨ０においても同様に高い値を有する。しかしながら、画像データＨＬ０，ＬＨ０，ＨＨ１においては原画像において高周波成分をなす粒状は除去されてしまうため、粒状に対応する画素の画素値はそれほど大きなものとはならない。したがって、画像データＨＬ０，ＬＨ０，ＨＨ１の画素値と閾値Ｔｈ１，Ｔｈ２，Ｔｈ３との比較の結果により、閾値Ｔｈ１，Ｔｈ２，Ｔｈ３よりも値が小さい画素に対応する最高周波数帯域画像の画素を、粒状が含まれる可能性がある画素と見なして、この画素について画素値の分散値Ｂを算出するものである。ここで、平坦な画像部分（すなわち信号値が略一定）に粒状が含まれている場合は分散値Ｂの値は大きくなり、粒状がそれほど含まれていない場合は分散値Ｂの値は小さくなる。
【００３４】
したがって、推定手段１３は分散値算出手段１２において算出された分散値Ｂを閾値ａと比較し、Ｂ≧ａのとき原画像には粒状が多いと推定し、Ｂ＜ａのとき原画像には粒状は少ないと推定する。そして、この推定結果を画像処理方法選択および／またはパラメータ設定手段４に出力する。
【００３５】
図４は鮮鋭度推定手段７の構成を示す概略ブロック図である。図４に示すように鮮鋭度推定手段７は、粒状推定手段３のウェーブレット変換手段１０と同様に、原画像データＳに対してウェーブレット変換を施して、原画像データＳにより表される原画像を多重解像度変換し、複数の周波数帯域毎の画像を表す画像データを得るウェーブレット変換手段３０と、ウェーブレット変換手段３０において得られた画像データＨＬ１，ＬＨ１（詳細は後述する）をそれぞれ所定の閾値と比較する比較手段３１と、比較手段３１による比較結果に基づいて、画像データＨＬ１，ＬＨ１と画像データＨＬ０，ＬＨ０とにおける画素の画素値の比Ｃを算出する比算出手段３２と、算出された比Ｃに基づいて原画像データＳの鮮鋭度を推定する推定手段３３とを備える。
【００３６】
ウェーブレット変換手段３０においては、粒状推定手段３のウェーブレット変換手段１０と同様にして原画像データＳに対してウェーブレット変換が施されるが、画像データＬＬ１に対してさらに下記のようにしてウェーブレット変換を施すものである。
【００３７】
すなわち、上述したHaar-Wavelet基底h(x)=-1,1およびg(x)＝1/2,1/2の次の階層のHaar-Wavelet基底であるh(x)′=-1,0,1およびg(x)′＝1/2,0,1/2を使用し、画像データＬＬ１にh(x)′およびh(y)′を畳み込み演算した画像データＨＨ１、画像データＬＬ１にh(x)′およびg(y)′を畳み込み演算した画像データＨＬ１、画像データＬＬ１にg(x)′およびh(y)′を畳み込み演算した画像データＬＨ１を下記の演算により求める。ここで、画像データＬＬ１を求める際に画素を間引いていれば、ウェーブレット基底として上述したh(x)=-1,1およびg(x)=1/2,1/2を用いればよいが、本実施形態においては画素を間引いていないため、間引きに相当する数の値０を挿入したものである。
【００３８】
ＨＨ１(x,y)＝h(y)′＊｛h(x)′＊ＬＬ１(x,y)｝
ＨＬ１(x,y)＝g(y)′＊｛h(x)′＊ＬＬ１(x,y)｝
ＬＨ１(x,y)＝h(y)′＊｛g(x)′＊ＬＬ１(x,y)｝
ここで、画像データＨＨ１、画像データＨＬ１および画像データＬＨ１は、それぞれ原画像の最高周波数帯域の次の周波数帯域におけるｘ方向ｙ方向ともに高周波画像、ｘ方向高周波ｙ方向低周波画像およびｘ方向低周波ｙ方向高周波画像を表すものとなる。
【００３９】
なお、この後ウェーブレット変換を繰り返すことにより、さらに低周波数帯域の画像データを求めることができるが、鮮鋭度推定手段７においては画像データＨＬ０，ＬＨ０，ＨＬ１，ＬＨ１のみを用いるため、ここで演算を終了する。また、本実施形態においては、ウェーブレット変換を行う際に画素を間引いていないため、画像データＨＨ０，ＨＬ０，ＬＨ０，ＬＬ０，ＨＨ１，ＨＬ１，ＬＨ１，ＬＬ２により表される画像は図３に示すように原画像と同一の解像度を有するものとなる。ここで、ＬＬ２は画像データＬＬ１にg(x)′およびg(y)′を畳み込み演算した画像データである。
【００４０】
比較手段３１においては、画像データＨＬ１，ＬＨ１の画素値が閾値Ｔｈ４，Ｔｈ５と比較され、ＨＬ１(x,y)＞Ｔｈ４かつＬＨ１(x,y)＞Ｔｈ５の２条件を満たす画素(x,y)が抽出される。
【００４１】
比算出手段３２は、比較手段３１において抽出された画素(x,y)の画素値と、この画素(x,y)に対応する画像データＨＬ０，ＬＨ０の画素値との比のうちいずれか大きい方の比Ｃを下記の式（１）により算出する。なお、本実施形態においては、図３に示すように各画像データＨＬ０，ＬＨ０，ＨＬ１，ＬＨ１の解像度は同一であるため、画素(x,y)の対応付けを正確かつ簡易に行うことができる。
【００４２】
Ｃ(x,y)＝ｍａｘ（ＨＬ０／ＨＬ１，ＬＨ０／ＬＨ１）    （１）
ここで、原画像に含まれるエッジなど、原画像の鮮鋭度推定の基準となる部分は画素値として比較的高い値を有し、低周波数帯域画像を表す画像データＨＬ１，ＬＨ１においても高い値を有する。逆に低周波数帯域画像を表す画像データＨＬ１，ＬＨ１において画素値の値が低ければそこにエッジなどの鮮鋭度推定の基準となるものが存在しないと考えることができる。したがって、低周波数帯域画像における画素値と閾値Ｔｈ４，Ｔｈ５との比較の結果により、閾値Ｔｈ４，Ｔｈ５よりも値が大きい画素を鮮鋭度推定の基準となる画素と見なし、この画素の画素値とこの画素に対応する最高周波数帯域画像を表す画像データＨＬ０，ＬＨ０の画素の画素値との比を求める。
【００４３】
ここで、最高周波数帯域画像は、原画像における最高周波成分を担持することから比較的狭い範囲の信号の変化量を表し、低周波数帯域画像は最高周波数帯域画像よりも低周波成分を担持することから広い範囲の信号の変化量を表すものとなる。例えば、図５（ａ）に示すように、非常に急峻なエッジにおいては、狭い範囲Ａにおいても広い範囲Ｂにおいても、その信号の変化量ΔＡ，ΔＢは同一となる。これに対し、図５（ｂ）に示すように、なだらかに信号が変化している部分においては、周波数成分としては急峻なエッジよりも低い周波数帯域のものとなるため、狭い範囲Ａにおける信号の変化量ΔＡよりも広い範囲Ｂにおける信号の変化量ΔＢの方が大きくなる。したがって、図５（ａ）に示す急峻なエッジにおいては、エッジ部分に対応する画素の画素値は最高周波数帯域画像と低周波数帯域画像とにおいて略同一の値となり、図５（ｂ）に示すなだらかに信号が変化している部分においては、信号変化位置に対応する画素の画素値は、最高周波数帯域画像よりも低周波数帯域画像の方が大きな値となる。
【００４４】
このため、上記式（１）により求められる低周波数帯域画像の画素の画素値に対する最高周波数帯域画像の画素の画素値の比Ｃは、急峻なエッジすなわち鮮鋭度が高い部分ほどその値は１に近くなり、信号の変化が滑らかな部分すなわち鮮鋭度が低い部分ほどその値は１よりも小さくなる。
【００４５】
したがって、推定手段３３は比算出手段３２において算出された比Ｃを１と比較し、Ｃ＝１のとき（または１に近いとき）その部分は鮮鋭度が高いと推定し、Ｃ＜＜１のときその部分は鮮鋭度が低いと推定する。そして、この推定結果を画像処理方法選択および／またはパラメータ設定手段４に出力する。
【００４６】
画像処理方法選択および／またはパラメータ設定手段４は、粒状推定手段３および／または鮮鋭度推定手段７による原画像の粒状量および／または鮮鋭度の推定結果に基づいて、画像処理方法の選択および／またはパラメータの設定を行う。具体的には、粒状量については、粒状が多い場合は、粒状を低減させるスムージング処理を選択し、あるいは原画像を拡大、縮小する場合には、滑らかさを重視したＢスプライン補間演算を選択するとともにこれのパラメータを決定する。また、粒状が少ない場合は、スムージング処理は実行しないか、あるいは原画像を拡大、縮小する場合は、鮮鋭度重視のCubicスプライン補間演算を選択する。
【００４７】
ここで、Cubicスプライン補間演算およびＢスプライン補間演算について説明する。本実施形態において使用される原画像データＳは、等間隔の周期でサンプリングされた一方向に配列されたサンプリング点（画素）Ｘk-2，Ｘk-1 ，Ｘk ，Ｘk+1 ，Ｘk-2…に対応したデジタル画像データ（Ｓk-2，Ｓk-1 ，Ｓk，Ｓk+1 Ｓk+2 …）であるとする。
【００４８】
Cubicスプライン補間演算は、オリジナルのサンプリング点（画素）Ｘk 〜Ｘk+1 間に設けられた補間点Ｘp の補間データＹ′を表す３次のCubicスプライン補間演算式（２）における補間データＹk-1 ，Ｙk ，Ｙk+1 ，Ｙk+2 にそれぞれ対応する補間係数ｃk-1 ，ｃk ，ｃk+1 ，ｃk+2 を、下記にそれぞれ示す演算により求めるものである。
【００４９】
Ｙ′＝ｃk-1 Ｙk-1 ＋ｃk Ｙk ＋ｃk+1 Ｙk+1 ＋ｃk+2 Ｙk+2   （２）
ｃk-1 ＝（−ｔ3 ＋２ｔ2 −ｔ）／２
ｃk ＝（３ｔ3 −５ｔ2 ＋２）／２
ｃk+1 ＝（−３ｔ3 ＋４ｔ2 ＋ｔ）／２
ｃk+2 ＝（ｔ3 −ｔ2）／２
（但し、ｔ（０≦ｔ≦１）は格子間隔を１とし、画素Ｘk を基準としたときの補間点Ｘp の画素Ｘk+1 方向への位置を示す。）
Ｂスプライン補間演算は、オリジナルのサンプリング点Ｘk 〜Ｘk+1 間に設けられた補間点Ｘp の補間データＹ′を表す３次のＢスプライン補間演算式（３）における補間データＹk-1、Ｙk 、Ｙk+1 、Ｙk+2 にそれぞれ対応する補間係数ｂk-1 ，ｂk ，ｂk+1 ，ｂk+2 を、下記にそれぞれ示す演算により求めるものである。
【００５０】
Ｙ′＝ｂk-1 Ｙk-1 ＋ｂk Ｙk ＋ｂk+1 Ｙk+1 ＋ｂk+2 Ｙk+2   （３）
ｂk-1 ＝（−ｔ3 ＋３ｔ2 −３ｔ＋１）／６
ｂk ＝（３ｔ3 −６ｔ2 ＋４）／６
ｂk+1 ＝（−３ｔ3 ＋３ｔ2 ＋３ｔ＋１）／６
ｂk+2 ＝ｔ3 ／６
（但し、ｔ（０≦ｔ≦１）は格子間隔を１とし、画素Ｘk を基準としたときの補間点Ｘp の画素Ｘk+1 方向への位置を示す。）そして、画像処理選択および／またはパラメータ設定手段４においては、ｃk-1 ，ｃk ，ｃk+1 ，ｃk+2 あるいはｂk-1 ，ｂk ，ｂk+1 ，ｂk+2 がパラメータとして設定される。
【００５１】
一方、鮮鋭度については、鮮鋭度推定手段７による原画像の鮮鋭度の推定結果、すなわち鮮鋭度の高さに基づいて、下記の式（４）に示すアンシャープマスキング処理の強調係数Sk(x,y)の値をパラメータとして設定する。
【００５２】
Ｓ′(x,y)＝Ｓ(x,y)＋Sk(x,y)×｛Ｓ(x,y)−Ｓus(x,y)｝      （４）
Ｓ′(x,y)：処理済み画像データ
Ｓus(x,y)：画像データＳ(x,y)を５×５のマスクにてぼかした画像データ
ここで、図６に示すように、Ｃ(x,y)＝１のときは０、Ｃ(x,y)＝０．５のときは１、Ｃ(x,y)＝０のときは３の値を採るような関数ｆ（Ｃ）を設定し、Sk(x,y)＝ｆ（Ｃ）としてシャープネスの強調係数Sk(x,y)を求める。なお、上記閾値Ｔｈ１，Ｔｈ２との比較において抽出されなかった画素(x,y)は、鮮鋭度推定の基準となり得ない平坦な領域（信号値が一定）であり、シャープネス強調は不要であるため、その部分についてはSk(x,y)＝０とすればよい。
【００５３】
なお、この際、上記閾値Ｔｈ４，Ｔｈ５との比較において抽出された画素の領域（以下抽出領域とする）と抽出されなかった画素の領域（以下非抽出領域とする）とにおいて画像が不連続となることを防止するために、非抽出領域であっても、その領域からあらかじめ設定した距離δ内に抽出領域が存在する場合は、最も近い位置にある抽出領域における強調係数Sk(x1,y1)および非抽出領域と抽出領域との距離Δに基づいて、下記の式（５）に示すように強調係数Sk(x,y)を設定することが好ましい。
【００５４】
Sk(x,y)＝Sk(x1,y1)×｛（δ−Δ）／δ｝                （５）
また、画像処理方法選択および／またはパラメータ設定手段４は、原画像を拡大、縮小する場合においては、鮮鋭度が高いと推定された場合には、滑らかさを重視したＢスプライン補間演算を選択するとともにこれのパラメータを決定する。一方、鮮鋭度が低いと推定された場合には、鮮鋭度重視のCubicスプライン補間演算を選択する。
【００５５】
画像処理手段５は、画像処理選択および／またはパラメータ設定手段４において選択された画像処理および／またはパラメータに基づいて、画像処理および／または補間演算を原画像データＳに対して施して処理済み画像データＳ′を得る。
【００５６】
次いで、本実施形態の動作について説明する。画像入力手段２から入力された原画像データＳは、粒状推定手段３および／または鮮鋭度推定手段７に入力される。粒状推定手段３においては、ウェーブレット変換により原画像データＳが複数の周波数帯域毎の画像を表す画像データに分解され、上述したようにして最高周波数帯域画像を表す画像データＨＨ０の分散値Ｂが算出され、原画像の粒状量が推定される。一方、鮮鋭度推定手段７においては、粒状推定手段３と同様に、ウェーブレット変換により原画像データＳが複数の周波数帯域毎の画像を表す画像データに分解され、上述した式（１）により比Ｃが算出され、原画像の鮮鋭度が推定される。
【００５７】
画像処理方法選択および／またはパラメータ設定手段４は、粒状推定手段３において推定された粒状量に基づいて画像処理手段５において原画像データＳに施す画像処理を選択する、あるいは原画像を拡大、縮小する場合は、補間演算の種類を選択するとともに補間演算のパラメータを設定する。また、鮮鋭度推定手段７において推定された鮮鋭度に基づいて、画像処理手段５において原画像データＳに施すアンシャープマスキング処理の強調係数Sk(x,y)を設定する、あるいは原画像を拡大、縮小する場合は、補間演算の種類を選択するとともに補間演算のパラメータを設定する。
【００５８】
画像処理手段５は、画像処理方法選択および／またはパラメータ設定手段４において選択された画像処理および／または補間演算を原画像データＳに対して施して処理済み画像データＳ′を得る。処理済み画像データＳ′は画像出力手段６よりプリンタ、ＣＲＴに出力されて可視像として再生されるか、またはネットワーク上に転送される。
【００５９】
このように、本実施形態においては、原画像を表す原画像データＳから原画像の粒状量および／または鮮鋭度を推定しているため、この推定結果を用いて画像処理手段５においてその粒状量および／または鮮鋭度に応じたスムージング処理、アンシャープマスキング処理などの画像処理および／または補間演算などを原画像データＳに対して施すことができる。これにより、その原画像がいずれの入力装置から入力されたものであるか、すでに画像処理が施されものであるかが不明な場合であっても、その原画像に最適な粒状および／または鮮鋭度に関する画像処理を施すことができる。
【００６０】
次いで、本発明の他の実施形態について説明する。
【００６１】
図７は図１に示す画像処理装置１を適用した画像合成装置の構成を示すブロック図である。図７に示すように本発明の他の実施形態は、複数の原画像データＳを上記画像処理装置１に入力し、原画像データＳのそれぞれについて上記実施形態と同様に最適な画像処理を施し、編集指示入力手段２１から入力された指示に基づいて、画像合成手段２０において処理が施された原画像データＳ′を合成して合成画像データＧを得、この合成画像データＧを画像出力手段６からプリンタなどに出力するものである。
【００６２】
このように、複数の原画像を合成する場合に、本実施形態による画像処理装置により各原画像データＳの粒状量および／または鮮鋭度を推定し、推定された粒状量および／または鮮鋭度に基づいて各原画像データＳに対して上記実施形態と同様に画像処理を施すことにより、それぞれ粒状および／または鮮鋭度に関する最適な画像処理が施された処理済み画像データＳ′を得ることができる。ここで、画像処理装置１に入力される各原画像データＳにより表される原画像の粒状が多い場合にはその粒状を低減するようなスムージング処理などの画像処理が施され、あるいは原画像を拡大、縮小する場合には上記Ｂスプライン補間演算が施される。また、粒状が少ない場合にはスムージング処理が施されないか、原画像を拡大、縮小する場合には上記Cubicスプライン補間演算が施される。一方、画像処理装置１に入力される各原画像データＳにより表される原画像の鮮鋭度が高いと推定された部分にはその鮮鋭度をそれほど強調しないような処理が施され、鮮鋭度が低いと推定された部分には鮮鋭度が強調される処理が施される。また、画像全体として鮮鋭度が高いときに原画像を拡大、縮小する場合は滑らかさを重視したＢスプライン補間演算を行い、鮮鋭度が低いときには鮮鋭度を重視したCubicスプライン補間演算を行うことができる。
【００６３】
したがって、処理済み画像データＳ′により表される各画像は、いずれの画像においても粒状および／または鮮鋭度が最適に補正されたものとなる。このため、画像合成手段２０において処理済み画像データＳ′を合成して得られる合成画像データＧを再生することにより、全ての原画像データＳに対して同一の画像処理を施した場合に生じる、部分毎の粒状バランス、鮮鋭度バランスの崩れをなくすことができ、これにより見た目に違和感のない画像を得ることができる。
【００６４】
なお、上記画像合成装置の実施形態においては、編集指示入力手段２１からの指示を画像合成手段２０に入力しているが、画像処理装置１にも入力し、指示された領域についてのみ画像処理を施すようにしてもよい。
【００６５】
また、上記各実施形態においては、原画像データＳの全画素についてウェーブレット変換を施しているが、主要な被写体が含まれる可能性の高い原画像の中央部分のみ、あるいは任意に選択された領域についてのみウェーブレット変換を施して、演算時間を短縮するようにしてもよい。
【００６６】
また、上記各実施形態においては、ウェーブレット変換を施す際に画素を間引いていないが、画素を間引いて低周波数帯域ほど縮小された画像となるようにウェーブレット変換を施すようにしてもよい。
【００６７】
さらに、上記各実施形態においては、ウェーブレット変換を用いて原画像を多重解像度空間に変換しているが、ラプラシアンピラミッドあるいはフーリエ変換を用いてもよい。
【００６８】
さらにまた、上記各実施形態においては、比較手段１１において、最高周波数帯域画像を表す画像データＨＨ０よりも１段階低周波数帯域の画像データＨＬ０，ＬＨ０，ＬＬ０と閾値Ｔｈ１，Ｔｈ２，Ｔｈ３とを比較しているが、ウェーブレット変換手段１０においてさらに低周波数帯域の画像を表す画像データを求め、この画像データと閾値Ｔｈ１，Ｔｈ２，Ｔｈ３とを比較するようにしてもよい。
【００６９】
また、上記各実施形態においては、比較手段３１において、最高周波数帯域画像を表す画像データＨＬ０，ＬＨ０よりも１段階低周波数帯域の画像データＨＬ１，ＬＨ１と閾値Ｔｈ４，Ｔｈ５とを比較しているが、ウェーブレット変換手段３０においてさらに低周波数帯域の画像を表す画像データを求め、この画像データと閾値Ｔｈ４，Ｔｈ５とを比較するようにしてもよい。
【００７０】
さらに、上記実施形態においては、粒状推定手段３および鮮鋭度推定手段７の各々にウェーブレット変換手段１０，３０を設けているが、ウェーブレット変換手段を１つのみ設け、単一のウェーブレット変換手段の出力を粒状推定手段３および鮮鋭度推定手段７の双方にて用いるようにしてもよい。
【００７１】
さらにまた、上記実施形態の粒状推定手段３および鮮鋭度推定手段７において行われる粒状の推定および鮮鋭度の推定は、上記ウェーブレット変換を用いた手法に限定されるものではなく、種々の手法を用いることができるものである。
【図面の簡単な説明】
【図１】  本発明の実施形態による画像処理装置の構成を示す概略ブロック図
【図２】  本発明の実施形態による粒状推定手段の構成を示す概略ブロック図
【図３】  ウェーブレット変換の状態を示す図
【図４】  本発明の実施形態による鮮鋭度推定手段の構成を示す概略ブロック図
【図５】  高周波数帯域画像と低周波数帯域画像とにおける信号値の大小を説明するための図
【図６】  強調係数Skを設定するための関数を示すグラフ
【図７】  本発明の他の実施形態による画像処理装置の構成を示す概略ブロック図
【符号の説明】
１    画像処理装置
２    画像入力手段
３    粒状推定手段
４    画像処理方法選択および／またはパラメータ設定手段
５    画像処理手段
６    画像出力手段
７    鮮鋭度推定手段
１０    ウェーブレット変換手段
１１    比較手段
１２    分散値算出手段
１３    推定手段
２０    画像合成手段
２１    編集指示入力手段
３０    ウェーブレット変換手段
３１    比較手段
３２    比算出手段
３３    推定手段
【図１】
0003759310.tif 000002
【図２】
0003759310.tif 000003
【図３】
0003759310.tif 000004
【図４】
0003759310.tif 000005
【図５】
0003759310.tif 000006
【図６】
0003759310.tif 000007
【図７】
0003759310.tif 000008