JP2004201290A - 画像処理方法および画像処理装置、並びに画像処理プログラムおよび記録媒体 - Google Patents

Abstract

【課題】 デジタル画像の画像ぼけを検出し補正を行うことで、効率的画像処理の演算および適応的な画質改善を行うことが可能な、画像処理方法および画像処理装置、並びに画像処理プログラムおよびそれを記録した記録媒体を提供する。
【解決手段】 画像処理装置は、入力デジタル画像ファイルのＤＣＴ係数の配列から複数のＤＣＴ係数を読み込むＤＣＴ係数読み込み部１と、上記複数のＤＣＴ係数から、全てのブロックの中で最小のぼけ量を持つブロックのぼけ量を数値として表す第１のぼけ指標を計算する第１のぼけ指標計算部２と、上記複数のＤＣＴ係数から、上記デジタル画像の全体としてのぼけ量を数値として表す第２のぼけ指標を計算する第２のぼけ指標計算部３と、上記入力デジタル画像ファイルに対し、上記第１のぼけ指標および上記第２のぼけ指標の値に基づいて選択的にフィルタを適用するフィルタ処理部４とを備えている。
【選択図】 図９

Description

本発明は、デジタル画像の処理を行う画像処理方法および画像処理装置、並びに画像処理プログラムおよびそれを記録した記録媒体に関するものであり、特に、デジタル画像の画像ぼけの量を予測および低減する画像処理方法および画像処理装置、並びに画像処理プログラムおよびそれを記録した記録媒体に関するものである。

画像のデジタル表示は、画像のアナログ表示、例えば、ネガフィルム、ポジスライドなどに対して多様な利点をもたらす。例えば、デジタル画像は、画像の劣化なしに電子的に格納し、送り、コピーすることができる。物体および事象のデジタル画像は、デジタルカメラなどの旧来のデバイスによって取り込んだり、ネガからスキャニングしたり、あるいはコンピュータで作成したりすることができる。デジタル画像は、また、カラーやコントラストなどの画像特性を変更したり、カメラのフラッシュに起因する赤目現象や画像ぼけなどのような画像の欠陥または人為的現象を低減させたりするために、デジタル的に処理および操作することができる。

画像ぼけは、画像のディテールおよび鮮鋭度の損失による画質の劣化を認識する上で最も重要な要素の１つである。画像ぼけは、例えば、被写体に焦点が合わせられていないカメラレンズ、前景および／または背景の絞りぼけ、カメラに対する被写体の相対的な運き、さらには大気中のもやの影響を含む多くの事の結果として発生しうるが、カラー写真などの旧来のアナログ画像における画像ぼけは、歴史的に目立って低減することができなかった。しかしながら、デジタル画像の画像ぼけは目立って低減させることができる。画像ぼけを低減させるための既存の処理技術は、典型的には、ぼけ予測部とぼけ低減部との両方を含んでいる。ぼけ予測部が所定の閾値を超える画像ぼけを検出すれば、画像がぼけ低減部によって処理される。画像ぼけを低減させるためのいくつかの既存の処理技術は、画像ぼけを定量化し、それに従って補正を調整している。画像ぼけを低減させるための他の技術は、画像ぼけの予測量が所定の閾値を下回るまで、画像を繰り返しぼけ低減部に通して循環させる反復プロセスを用いている。

画像ぼけを予測する１つの既存のシステムは、空間領域で行われている。このシステムでは、鋭いエッジを使用して点拡がり関数が生成され、この関数から画像ぼけの量が予測される。このシステムが頻繁に直面する問題は、実世界の画像の多くが、画像ぼけを正確に予測するのに十分な鋭いエッジを欠いているということである。

画像ぼけを予測する他のシステムは、周波数領域における画像のパワースペクトルから画像ぼけの量を予測している。このシステムは、ノンゼロ周波数におけるパワーが概ね周波数の二乗に逆比例して変化するとすれば、ほとんどの実世界の非ぼけ画像のパワースペクトル比較的不変であるという仮定を利用している。逆に、ぼけ画像は、周波数が増加するに従って、パワーがより急速に低下する特性を示す。この現象を、ぼけ画像および非ぼけ画像の両方のパワースペクトルを示す図１で説明する。図１から分かるように、ぼけ画像のパワーは、非常に低い周波数では非ぼけ画像のパワーとほぼ等しいが、およそ０．０５サイクル／サンプル以上の周波数では急速に低下して非ぼけ画像のパワーより低くなる。画像のパワースペクトルから画像ぼけを予測するシステムは、シーン統計量を使用するので、一般に、鋭いエッジ等の画像構成から得られる点拡がり関数に基づいて画像ぼけを予測するシステムよりも信頼性が高いと考えられる。しかしながら、画像のパワースペクトルから画像ぼけを予測するシステムは、計算上複雑であり、このことは、写真プリンタなどのような多くの応用例における実施を困難にしている。

画像ぼけを予測するためのさらに他のシステムは、チャン（Zhang）等の特許文献１に開示されている。このシステムは、デジタル画像を格納しているデジタル画像ファイルに含まれる全てのノンゼロＤＣＴ（離散コサイン変換）係数の重み付け出現ヒストグラムからデジタル画像の画像ぼけの量を予測している。

また、特許文献２には、ＤＣＴ係数を用いてぼけ領域を検出する方法が開示されている。
米国特許第６，２９８，１４５ Ｂ１号（２００１年１０月２日発行） 特開平１０−２８５５８７号公報（１９９７年３月３１日公開）

しかしながら、特許文献１のシステムは、特許文献１で述べられているように、「（ＤＣＴ）係数の値が、それらの表す画像のタイプと密接に関連する値になるまで、（ＤＣＴ）係数を直接使用しない」。代わりに、ＤＣＴ成分が前述の出現ヒストグラムを構築するのに使用される。このヒストグラムは、次に数学的に操作することが必要であり、このことは計算的な複雑さを付加する。

また、特許文献２の方法は、画像のぼけ領域を検出することは可能であるが、ぼけ補正を行う手段を有していないため、画像のぼけを補正することができないという問題点を有している。

そこで、望まれているのは、計算上効率的に画像ぼけを予測および低減するシステムである。すなわち、本発明の目的は、画像処理の演算および画質改善を効率的に行うことが可能な画像処理方法および画像処理装置、並びに画像処理プログラムおよびそれを記録した記録媒体を提供することにある。

本発明に係る画像処理方法は、上記課題を解決するために、デジタルファイル中にＤＣＴ係数の配列として格納されたデジタル画像であり、上記ＤＣＴ係数の配列が周波数成分の格子を表す複数のブロック内に配列され、上記各ブロックが画像の部分を表すデジタル画像に対し、デジタル画像のぼけを検出および修正する画像処理方法であって、（ａ）上記ＤＣＴ係数の配列から、各ブロックから少なくとも２つのＤＣＴ係数を読み込むように複数のＤＣＴ係数を読み込むステップと、（ｂ）上記複数のＤＣＴ係数から、全てのブロックの中で最小のぼけ量を持つブロックのぼけ量を数値として表す第１のぼけ指標を計算するステップと、（ｃ）上記複数のＤＣＴ係数から、上記デジタル画像の全体としてのぼけ量を数値として表す第２のぼけ指標を計算するステップと、（ｄ）上記デジタル画像に対し、上記第１のぼけ指標および上記第２のぼけ指標の値に基づいて選択的にフィルタを適用するステップとを含むことを特徴としている。

また、本発明に係る画像処理方法は、上記画像処理方法において、上記選択的にフィルタを適用するステップでは、複数の異なる種類のフィルタの中から１つのフィルタを選択してそのフィルタを上記デジタル画像に適用する方法であってもよい。

また、本発明に係る画像処理方法は、上記画像処理方法において、格子内の代表的な周波数成分を選択するステップをさらに含み、上記第１のぼけ指標は、画像内の全てのブロックにおける上記周波数成分でのＤＣＴ係数の中の最大値である方法であってもよい。

また、本発明に係る画像処理方法は、上記画像処理方法において、格子内の２つの代表的な周波数成分を選択し、これら２つの代表的な周波数成分におけるＤＣＴ係数の標準偏差を画像全体で計算するステップをさらに含み、上記第２のぼけ指標は、これら標準偏差の比である方法であってもよい。

また、本発明に係る画像処理方法は、上記画像処理方法において、（ｅ）周波数成分の格子を、周波数が増加する複数領域に分割するステップと、（ｆ）各ブロックの各領域におけるＤＣＴ係数をまとめるステップとをさらに含み、第１のぼけ指標を計算するステップおよび第２のぼけ指標を計算するステップでは、上記まとめられたＤＣＴ係数から、上記第１のぼけ指標および第２のぼけ指標を計算する方法であってもよい。

また、本発明に係る画像処理方法は、上記画像処理方法において、上記各領域におけるＤＣＴ係数をまとめるステップは、上記各領域内におけるＤＣＴ係数の平均値を得るステップを含む方法であってもよい。

また、本発明に係る画像処理方法は、上記画像処理方法において、第１のぼけ指標が、上記複数領域の中から選択された１つの領域内における最大の平均ＤＣＴ係数である方法であってもよい。

また、本発明に係る画像処理方法は、上記画像処理方法において、複数領域の中から選択された２つの領域の各々においてＤＣＴ係数の標準偏差を得るステップを含む方法であってもよい。

また、本発明に係る画像処理方法は、上記画像処理方法において第２のぼけ指標が、得られた標準偏差同士の比である方法であってもよい。

本発明に係る画像処理方法は、上記課題を解決するために、デジタルファイル中にＤＣＴ係数の配列として格納されたデジタル画像であり、上記ＤＣＴ係数の配列が水平周波数成分および垂直周波数成分の格子を表す複数のブロック内に配列され、上記各ブロックが画像の部分を表すデジタル画像に対し、デジタル画像のぼけを検出および修正する画像処理方法であって、（ａ）ＤＣＴ係数の２次元のヒストグラムを得るステップと、（ｂ）上記２次元のヒストグラムを１次元のヒストグラムにまとめるステップと、（ｃ）上記１次元のヒストグラムから少なくとも１つのぼけ指標を計算するステップと、（ｄ）上記デジタル画像に対し、上記少なくとも１つのぼけ指標に基づいて選択的にフィルタを適用するステップとを含むことを特徴としている。

また、本発明に係る画像処理方法は、上記画像処理方法において、２次元のヒストグラムを１次元のヒストグラムにまとめるステップでは、２次元ヒストグラムを周波数が概ね増加する複数領域に分割して、各ブロックの各領域内のＤＣＴ係数を平均化することによって１次元のヒストグラムを得る方法であってもよい。

また、本発明に係る画像処理方法は、上記画像処理方法において、上記ぼけ指標が、上記複数領域から選択された１つの領域内における最大の平均ＤＣＴ係数である方法であってもよい。

また、本発明に係る画像処理方法は、上記画像処理方法において、上記ぼけ指標が、上記複数領域から選択された２つの領域間における平均ＤＣＴ係数の比である方法であってもよい。

また、本発明に係る画像処理方法は、上記画像処理方法において、第１のぼけ指標および第２のぼけ指標を計算するステップをさらに含む方法であってもよい。

また、本発明に係る画像処理方法は、上記画像処理方法において、第１のぼけ指標および第２のぼけ指標を用いる方法であってもよい。

また、本発明に係る画像処理方法は、上記画像処理方法において、上記第１のぼけ指標が、上記複数領域から選択された１つの領域内における最大の平均ＤＣＴ係数であり、上記第２のぼけ指標が、上記複数領域から選択された２つの領域間における平均ＤＣＴ値の比である方法であってもよい。

また、本発明に係る画像処理方法は、上記画像処理方法において、上記フィルタが、異なる複数種類のフィルタの中から第１のぼけ指標および第２のぼけ指標に基づいて選択される方法であってもよい。

ところで、上記画像処理方法の各ステップは、ハードウェアで実現された画像処理装置に実行させてもよいし、プログラムの指示によってコンピュータに実行させてもよい。具体的には、本発明に係る画像処理プログラムは、上記画像処理方法の各ステップをコンピュータに実行させる画像処理プログラムであり、本発明に係るコンピュータ読み取り可能な記録媒体には、当該画像処理プログラムが記録されている。

本発明に係る画像処理装置は、上記課題を解決するために、デジタルファイル中にＤＣＴ係数の配列として格納されたデジタル画像であり、上記ＤＣＴ係数の配列が周波数成分の格子を表す複数のブロック内に配列され、上記各ブロックが画像の部分を表すデジタル画像に対し、デジタル画像のぼけを検出および修正する画像処理装置であって、上記ＤＣＴ係数の配列から、各ブロックから少なくとも２つのＤＣＴ係数を読み込むように複数のＤＣＴ係数を読み込むＤＣＴ係数読み込み部と、上記複数のＤＣＴ係数から、全てのブロックの中で最小のぼけ量を持つブロックのぼけ量を数値として表す第１のぼけ指標を計算する第１のぼけ指標計算部と、上記複数のＤＣＴ係数から、上記デジタル画像の全体としてのぼけ量を数値として表す第２のぼけ指標を計算する第２のぼけ指標計算部と、上記デジタル画像に対し、上記第１のぼけ指標および上記第２のぼけ指標の値に基づいて選択的にフィルタを適用するフィルタ処理部とを備えることを特徴としている。

本発明に係る画像処理装置は、上記課題を解決するために、デジタルファイル中にＤＣＴ係数の配列として格納されたデジタル画像であり、上記ＤＣＴ係数の配列が水平周波数成分および垂直周波数成分の格子を表す複数のブロック内に配列され、上記各ブロックが画像の部分を表すデジタル画像に対し、デジタル画像のぼけを検出および修正する画像処理装置であって、ＤＣＴ係数の２次元のヒストグラムを作成する２次元ヒストグラム作成部と、上記２次元のヒストグラムを１次元のヒストグラムに変換するヒストグラム変換部と、上記１次元のヒストグラムから少なくとも１つのぼけ指標を計算するぼけ指標計算部と、上記デジタル画像に対し、上記少なくとも１つのぼけ指標に基づいて選択的にフィルタを適用するフィルタ処理部とを備えることを特徴としている。

本発明に係る画像処理方法は、以上のように、（ａ）上記ＤＣＴ係数の配列で構成される、各ブロックから少なくとも２つのＤＣＴ係数を読み込むように複数のＤＣＴ係数を読み込むステップと、（ｂ）上記複数のＤＣＴ係数から、全てのブロックの中で最小のぼけ量を持つブロックのぼけ量を数値として表す第１のぼけ指標を計算するステップと、（ｃ）上記複数のＤＣＴ係数から、上記デジタル画像の全体としてのぼけ量を数値として表す第２のぼけ指標を計算するステップと、（ｄ）上記デジタル画像に対し、上記第１のぼけ指標および上記第２のぼけ指標の値に基づいて選択的にフィルタを適用するステップとを含んでいるので、ＤＣＴ係数を用いて符号化された画像データを復号化する際に、新規に複雑な処理を追加することなくぼけ検出およびぼけ補正を行うことができるので、処理の演算効率が向上し、且つ第１および第２のぼけ指標に基づき選択的にフィルタを適用することができるので、ぼけの程度に応じて適応的に画質改善を行うことができる。それゆえ、画像処理の演算を効率的に行い、且つ適応的に画質改善を行うことができるという効果を奏する。

また、本発明に係る画像処理方法は、以上のように、（ａ）ＤＣＴ係数の２次元のヒストグラムを得るステップと、（ｂ）上記２次元のヒストグラムを１次元のヒストグラムにまとめるステップと、（ｃ）上記１次元のヒストグラムから少なくとも１つのぼけ指標を計算するステップと、（ｄ）上記デジタル画像に対し、上記少なくとも１つのぼけ指標に基づいて選択的にフィルタを適用するステップとを含んでいるので、２次元ヒストグラムのままぼけ指標を計算するのではなく１次元にまとめて計算できるので、演算効率が向上し、且つ少なくとも１つのぼけ指標に基づいて選択的にフィルタを適用することができるので、適応的に画質改善を行うことができる。それゆえ、画像処理の演算を効率的に行い、且つ適応的に画質改善を行うことが可能な画像処理方法を提供できるという効果を奏する。

また、本発明に係る画像処理プログラムは、上記画像処理方法の各ステップをコンピュータに実行させるものであるので、上記画像処理方法と同様に、画像処理の演算を効率的に行い、且つ適応的に画質改善を行うことができるという効果を奏する。

また、本発明に係るコンピュータ読み取り可能な記録媒体は、上記画像処理プログラムを記録したものであるので、上記画像処理プログラムと同様に、画像処理の演算を効率的に行い、且つ適応的に画質改善を行うことができるという効果を奏する。

本発明に係る画像処理装置は、以上のように、上記ＤＣＴ係数の配列から、各ブロックから少なくとも２つのＤＣＴ係数を読み込むように複数のＤＣＴ係数を読み込むＤＣＴ係数読み込み部と、上記複数のＤＣＴ係数から、全てのブロックの中で最小のぼけ量を持つブロックのぼけ量を数値として表す第１のぼけ指標を計算する第１のぼけ指標計算部と、上記複数のＤＣＴ係数から、上記デジタル画像の全体としてのぼけ量を数値として表す第２のぼけ指標を計算する第２のぼけ指標計算部と、上記デジタル画像に対し、上記第１のぼけ指標および上記第２のぼけ指標の値に基づいて選択的にフィルタを適用するフィルタ処理部とを備えているので、ＤＣＴ係数を用いて符号化された画像データを復号化する際に、新規に複雑な処理を追加することなくぼけ検出およびぼけ補正を行うことができるので、処理の演算効率が向上し、且つ第１および第２のぼけ指標に基づき選択的にフィルタを適用することができるので、ぼけの程度に応じて適応的に画質改善を行うことができる。それゆえ、画像処理の演算を効率的に行い、且つ適応的に画質改善を行うことが可能な画像処理方法を提供できるという効果を奏する。

本発明に係る画像処理方法は、以上のように、ＤＣＴ係数の２次元のヒストグラムを作成する２次元ヒストグラム作成部と、上記２次元のヒストグラムを１次元のヒストグラムに変換するヒストグラム変換部と、上記１次元のヒストグラムから少なくとも１つのぼけ指標を計算するぼけ指標計算部と、上記デジタル画像に対し、上記少なくとも１つのぼけ指標に基づいて選択的にフィルタを適用するフィルタ処理部とを備えているので、２次元ヒストグラムのままぼけ指標を計算するのではなく、１次元にまとめて計算できるので、演算効率が向上し、且つ少なくとも１つのぼけ指標に基づいて選択的にフィルタを適用することができるので、適応的に画質改善を行うことができる。それゆえ、画像処理の演算を効率的に行い、且つ適応的に画質改善を行うことが可能な画像処理方法を提供できるという効果を奏する。

図２は、入力画像がぼけ予測モジュールに送られ、ぼけが存在するか否かが判定され、ぼけが存在すればぼけの量が予測される典型的なぼけ予測・低減システムの一形態を示す。図２において、実線の矢印はこのシステムにおいて必須の信号の流れを示し、破線の矢印は任意の信号の流れを示す。図２に示すように、本発明の画像処理装置の実施の一形態としてのぼけ予測・低減システム２０は、ぼけ予測部２１とぼけ低減部２２とを備えている。

ぼけ予測部（ぼけ予測モジュール）２１は、ぼけ低減プロセスを適用すべきか否かを判定し、適用すべきであれば、ぼけをどの程度低減させるべきかを決定する。ぼけ低減部２２は、ぼけ予測部２１の決定に従って画像ぼけを低減させる。ぼけ予測部２１およびぼけ低減部２２による処理は、望みに応じて繰り返してもよい。例えば、図２に矢印で示すように、入力画像を、ぼけ閾値に到達するまでぼけ予測・低減システム２０に繰り返し通して繰り返し処理を行ってもよい。

ほとんどのデジタル画像は、次式

で表される順方向の離散コサイン変換（ＤＣＴ）関数を用いて画像を符号化し、かつ、次式

で表される逆ＤＣＴ変換関数を用いて画像を復号する画像形式に格納される。

上記順方向離散コサイン変換は、空間領域におけるデジタル画像の画素情報を周波数領域におけるエネルギーまたはパワーの情報に変換する。周波数領域へのデジタル画像の変換によって、画像が顕著に圧縮可能となり、その結果として格納スペースが節約できる。例として、空間領域に格納された白黒のデジタル画像は、全ての画素が、グレイスケール上における画素の濃淡を表す任意のビット長の数で符号化されることを必要とする。便宜上、デジタル画像は、通常、デジタル画像ファイルの中に格納されるときに、８×８画素のブロックに組織化される。したがって、白黒画像に関する画素情報をデジタル画像ファイルに格納するのに８ビットグレイスケールが使用される場合、デジタル画像ファイルは、連続した８×８の配列に組織化され、かつ、対応する画素の濃淡を表す０〜２５５の間の値を各々がとる、一連の数を含む。これは、格納に多くの格納空スペースを必要とする。より高いビット長が使用されている場合、さらに大きいスペースが必要である。

空間領域に格納されたカラー画像は、画素の色情報を適切に保存するために、［Ｒ，Ｇ，Ｂ］または［Ｌ，Ｃｒ，Ｃｂ］などのような３つ（または何らかの数)の成分によって各画素を表すことができるが、この場合、さらに多くの格納スペースを必要とする。例えば、２４ビットのＹ，Ｃｒ，Ｃｂのデジタル画像ファイルがカラー画像を格納している場合、デジタルファイルは、３つの個別の８×８配列を含み、各８×８配列は選択された色座標系の１つの成分を表す。各配列の各数は、対応する画素の対応するカラー座標の濃淡を表す８ビットの数、０〜２５５になる。

順方向ＤＣＴ変換は、各配列を、空間領域の画素情報から、離散的な垂直周波数成分および水平周波数成分の周波数領域におけるエネルギー情報にそれぞれ変換する。空間領域に画像を格納する画像形式のように、周波数領域に画像を格納する画像形式は、典型的には、画像のパワー情報を８×８のブロックに組織化する。言い換えれば、画像のパワー情報は、図３に示すような８×８の格子内に図示している８つの離散的な垂直周波数成分および８つの離散的な水平周波数成分で格納することができる。図３に示すＣ００項１０は、画像のＤＣ（直流）周波数成分、すなわち、ゼロの水平周波数およびゼロの垂直周波数における画像パワーを表す。Ｃ０７項１２は、ゼロの垂直周波数および最高の水平周波数における画像のパワーを表し、Ｃ７０項１４は、ゼロの水平周波数および最高の垂直周波数における画像のパワーを表す。Ｃ７７項１６は、最高の垂直周波数および最高の水平周波数での画像のパワーを表す。

実際の画像中における大多数のエネルギーまたはパワーは低周波ブロックに含まれるので、周波数領域で符号化されたデジタル画像はより容易に圧縮できる。従って、デジタル画像ファイルが高い垂直周波数および水平周波数の画像情報を切り捨てられるものであれば、まだ、画像のディテールの大部分を保持しながら、表示のための空間領域に画像を復号できるであろう。１つの一般的に使用されているデジタル画像形式は、順方向ＤＣＴ変換を使用してＪＰＥＧ(Joint Photographic Experts Group)形式に画像を特徴付けているが、他のデジタル画像形式もまたこの変換を使用している。ＪＰＥＧ形式は、画質の顕著な劣化なしに画像を劇的に圧縮することを可能にするので、画像を特徴付けるのに広く使用されている。

上述した順方向ＤＣＴ変換は、パワースペクトル分析で使用される連続したフーリエ変換と際立って類似している。この事実から、本願発明者は、画像ぼけの量が画像のパワースペクトルから予測できるのと丁度同じように、ほとんどの圧縮されたデジタル画像ファイルに既に埋め込まれたＤＣＴ情報から画像ぼけを認識できるであろうと推論した。したがって、ＤＣＴ係数を含むデジタル画像の画像ぼけを予測するときに、パワースペクトル分析の実行にかかわる計算の複雑さが回避可能であろう。

この考えで、本願発明者は、デジタル画像の画像ぼけの量を予測するのに使用可能な、デジタル画像ファイル中に既に埋め込まれたさまざまな情報が本当にあると断定した。例として、図４および図５は、「ぼけ指標」と呼ばれる、そのような画像ぼけの尺度の１つを示す。図４は、複数のブロックを処理した結果として得られる非ぼけ画像に関する、各周波数におけるＤＣＴ値の標準偏差を示す表である。図５は、複数のブロックを処理した結果として得られるぼけ画像に関する、各周波数におけるＤＣＴ値の標準偏差を示す表である。これら２つの表の比較から分かるように、低周波でのぼけ画像と非ぼけ画像との標準偏差は、ほとんど同じである。しかしながら、周波数が垂直成分および水平成分の両方で増加するのに従って、ぼけ画像のための標準偏差が非ぼけ画像のものより劇的に低下する。したがって、デジタル画像ファイル中におけるＤＣＴ係数の標準偏差は、ぼけ指標として効果的に使用できる。画像ぼけを効果的に予測できることがわかった他の指標は、最大のＤＣＴ係数、最小のＤＣＴ係数、分散、および尖度またはモーメント(Kurtosis)を含んでいる。さらに、画像ぼけは、ＤＣＴ係数の標準偏差を２次元（水平および垂直）のＤＣＴ周波数成分の関数として表す表面の下の体積として図６に図示している画像ブロック内における係数の標準偏差の合計から予測できる。

ぼけ指標は、単一のＤＣＴ係数（生の、または、修正された）であってもよく、１セットの複数のＤＣＴ係数(生の、または、修正された)であってもよく、他のＤＣＴ係数に対する１つ以上のＤＣＴ係数の比較結果（生の、または、修正された）であってもよい。どの場合でも、ゼロ以外あるいはノンゼロのＤＣＴ係数の大きさの変化は、ぼけの予測および／または低減における変化をもたらしうる。

なお、ＪＰＥＧなどの典型的なデジタル画像形式では、デジタル画像ファイルは、画像の空間的な一部分に対する空間的なドメインの画素情報に各ブロックが対応する、ＤＣＴ係数の多くのブロックを含みうる。したがって、画像の異なる複数の空間部分からのＤＣＴ係数の１つ以上のセット(「生」の、または、統計的に修正された）から適切なぼけ指標を選択することで、上記技術は、画像の個々の領域での画像ぼけの検出を可能にする。多くの画像が画像の一部分の被写体部分に焦点を合わせる一方、画像の残りをぼかすよう意図されているので、これは強力なツールである。画像内における何らかの領域がぼけていないかを示すのに使用可能なぼけ指標は、所定の周波数におけるデジタル画像の最大のＤＣＴ係数である。ぼけの存在のために１以上のブロックのＤＣＴ係数を含む画像の個々の領域を調べる能力でもって、そのような部分的にぼけている画像を特定できるので、そのような画像をうかつに修正してしまうことはない。

図７は、ぼけ予測・低減技術の一例を示す。この例では、２次元(「２Ｄ」)のヒストグラムがデジタル画像ファイルから構築されているので、画像ファイル中の各ＤＣＴ係数はその対応する周波数成分と関連付けられている。そして、上記周波数成分は、２次元である。すなわち、上記周波数成分は、垂直周波数値および水平周波数値の両方を示す。

デジタル画像ファイル内の各画像ブロックは、図３で例示した２次元（８×８）ＤＣＴ周波数成分に対応する６４個の係数を含みうる。２次元ヒストグラムは、画像ファイルの中のＤＣＴ係数の各々を、それに関連付けられた特定の２次元周波数成分に従って組織化する。ＪＰＥＧなどのようなほとんどのデジタル画像符号化スキームは、順方向ＤＣＴ変換からＤＣＴ係数を計算した後にさらにＤＣＴ係数を修正している。ＪＰＥＧは、例えば、量子化テーブル（Ｑテーブル）に応じて各ブロック内のそれぞれのＤＣＴ係数を量子化した後、より効率的な圧縮を提供するためにハフマン符号化スキームで値を符号化する。従って、このプロセスでＪＰＥＧファイルが使用されている場合、２次元ヒストグラムが得られる前にハフマン符号化および量子化を取り除くべきである。

そこで、この２Ｄヒストグラムを、半径方向の周波数パワースペクトルと概ね同等な１次元（「１Ｄ」）のヒストグラムに変換してもよい。マッピング関数に従って２Ｄヒストグラムを１Ｄヒストグラムに変換してもよい。図８は、１つの適切なマッピング関数を示す。図５から分かるように、ブロック内でのそれらの位置に応じて、６４個のＤＣＴ周波数成分が複数のグループ（例えば０〜７のグループ）に分割される。これらのグループは、ＤＣ周波数成分に対して同心円状の領域へ大まかに組織化できる。この特定のマッピング関数内の複数領域中のＤＣＴ係数の各々が、所定の画像ブロックにおいて互いに近接した値を持つことが一般に観測されたので、この特定のマッピング関数内の領域が選択される。２次元ヒストグラムを１次元のヒストグラムに変換する１つの技術は、所定の画像ブロック内の各領域に対応する複数のＤＣＴ係数をまとめて平均化し、中央値を与えるものである。代わりに、単に、その領域に対する値を表すように各領域内の１つの値を選択してもよい。本質的には、これは、２次元ヒストグラムを無次元のヒストグラムに変換するか、さもなければ、データ量を低減させる。０〜７の値は、一般に増加する半径方向の周波数成分を表す。

ここで説明した複数ステップは、逆の順序で実行してもよい。すなわち、各画像ブロックはデジタル画像ファイルから読み込まれるので、図８に示すマッピング関数は、まず、各領域０〜７に対する中央値を得ることによってブロックを１次元に変換するのに使用することもできる。そこで、１次元のヒストグラムは、追加画像ブロックが読み込まれ、変換され、それらの中央値がそれらの関連領域に分類されるように、構成することもできる。

ＤＣＴ係数の異なるセットの選択は、望みに応じて部分的に重複していてもよいが、重複していないことが好ましい。これらのセットは、好ましくは、より高い周波数やより低い周波数などのような異なる複数の周波数成分を含む。異なる周波数値は、大きさの変化の比較を周波数の関数として容易に可能とする。

１次元のヒストグラムから、１つ以上のぼけ指標を得ることもできる。使用可能な第１のぼけ指標は、領域５内の最大のＤＣＴ係数（あるいは何らかの他の領域または係数)である。本願発明者は、第１の指標に対するこの領域内の１００の閾値が、デジタルカメラから得られる画像のぼけに対する適当な指標であると決定した。最大のＤＣＴ係数が閾値より小さいか否かを検査することによって、フォーカスには画像の何らかの部分が焦点のあった状態であるか否かを正確に判定することができる。すなわち、領域５内の最大のＤＣＴ係数が１００未満であるなら、典型的には画像全体がぼけているであろう。逆に、領域５内の最大のＤＣＴ係数が１００以上であれば、おそらく画像の少なくとも一部分が焦点のあった状態であろう。

使用可能な第２のぼけ指標は、第３の領域の標準偏差に対する第１の領域の標準偏差の比（または何らかの他の２つ以上の領域間の比較）である。ぼけ画像内におけるＤＣＴ係数の標準偏差は、周波数の増加に応じて非ぼけ画像よりも急速に低下する傾向がある。したがって、比較的高い比はぼけ画像を示すが、比較的低い比は焦点の合った画像を示す。本願発明者は、このぼけ指標のための５の閾値が、デジタルカメラからの画像のための画像ぼけの適切な指標であると決定した。標準偏差の指標は画像全体にわたって情報を平均化するので、このぼけ指標は、画像における焦点のあった孤立部分を検出する傾向がないが、代わりに全体的な画像ぼけを評価する傾向があるであろう。

前述した技術は、第１の指標として領域５内の最大のＤＣＴ係数を使用し、第２のぼけ指標として領域３のＤＣＴ係数の比率に対する領域１のＤＣＴ係数の標準偏差の比を使用していたが、他の技術は、他の領域からの他のぼけ指標を使用してもよい。同様に、最小の係数値などのような他の統計的な値もまた、ぼけ指標として使用できる。また、上述した技術は１次元のヒストグラムからぼけ指標を計算していたが、ヒストグラムなしにぼけ指標を導き出してもよい。例えば、選択されるぼけ指標は、領域５の最大の係数と、領域１の係数の平均に対する領域３の係数の平均の比とであってもよい。この場合、単に、データがデジタル画像ファイルから読み込まれるにしたがって、これらのぼけ指標を計算および更新することができる。

図７に示す２つのぼけ指標から異なるぼけ予測を計算してもよい。第１に、第１の指標が１００未満であり、第２のぼけ指標が５以上であれば、画像全体はかなりぼけているであろう。この条件が成立している場合、画像を強調するために最強の反ぼけフィルタを適用することができる。

第２に、第１の指標が１００未満であり、第２のぼけ指標が５未満であれば、画像全体がわずかにぼけているであろう。この場合、画像を強調するために中程度の反ぼけフィルタを適用することができる。

第３に、第１の指標が１００以上であり、第２のぼけ指標が３以上５未満であれば、画像中の物体に焦点が合っているが、画像はわずかに鮮鋭度を欠くであろう。そこで、この場合、画像を強調するために標準の鮮鋭度強調フィルタを適用することができる。

第４に、第１の指標が１００以上であり、第２のぼけ指標が５以上であれば、画像は鮮鋭度をかなり欠くであろう。この場合、強い鮮鋭度強調フィルタを適用することができる。

本実施形態では、第１のぼけ指標が閾値以上か否かの判定結果に基づいて、「反ぼけフィルタ」と「鮮鋭度強調フィルタ」とを切り替える。すなわち、第１のぼけ指標が閾値未満と判定された画像に対しては「反ぼけフィルタ」を、第１のぼけ指標が閾値以上と判定された画像に対しては「鮮鋭度強調フィルタ」を使用する。また、第２のぼけ指標によっても、それぞれのフィルタ特性を変更する。

この段落で説明したフィルタは、良く知られており、それらの設計は当業者にとって容易である。「反ぼけフィルタ」は、ぼけた画像を修復する目的で使用するフィルタである。「反ぼけフィルタ」と「鮮鋭度強調フィルタ」とのフィルタ特性の違いは、「反ぼけフィルタ」は、ぼけ補正をするフィルタ、「鮮鋭度強調フィルタ」は、画質改善の観点から鮮鋭度強調を行うフィルタである。そのため、両者の違いは、強調する周波数帯域にある。「反ぼけフィルタ」は、典型的に画像全体がぼけているであろうという予測に基づき、適応的な補正を行うために、あまり高周波領域を強調せず、人間がその画像に対してぼけていると感じ易い比較的中周波である領域に対し強調を行い（強調するバンドを限定している）である。これに対し、「鮮鋭度強調フィルタ」は、高周波領域に行くに従って強調度合いを高く設定するフィルタである。

なお、図７では、４つの説明したフィルタの中の１つを適用した後に終了する技術を示しているが、補正された画像で上記処理を繰り返し、さらなるぼけ補正が必要であるかどうかを確認することも可能である。

次に、本発明に係る画像処理装置の実施の一形態について、図９に基づいて説明する。

本形態の画像処理装置は、入力デジタル画像ファイルに対し、該入力デジタル画像ファイル中にＤＣＴ係数の配列として格納されたデジタル画像のぼけを検出する処理を行い、その検出結果に基づいてデジタル画像を修正する画像処理装置である。入力デジタル画像ファイルにおいて、上記ＤＣＴ係数の配列が水平周波数成分および垂直周波数成分の格子を表す複数のブロック内に配列され、上記各ブロックが画像の部分を表す。

図９に示すように、本形態の画像処理装置は、入力デジタル画像ファイルのＤＣＴ係数の配列から複数のＤＣＴ係数を読み込むＤＣＴ係数読み込み部１と、上記複数のＤＣＴ係数から第１のぼけ指標を計算する第１のぼけ指標計算部２と、上記複数のＤＣＴ係数から第２のぼけ指標を計算する第２のぼけ指標計算部３と、上記入力デジタル画像ファイルに対し、上記第１のぼけ指標および上記第２のぼけ指標の値に基づいて選択的にフィルタを適用するフィルタ処理部４とを備えている。

ＤＣＴ係数読み込み部１は、入力デジタル画像ファイルの各ブロックから少なくとも２つのＤＣＴ係数を読み込む。また、本形態では、第１のぼけ指標は、全てのブロックの中で最小のぼけ量を持つブロックのぼけ量を数値として表し、第２のぼけ指標は、上記デジタル画像の全体としてのぼけ量を数値として表す。

フィルタ処理部４は、反ぼけ効果（ぼけをなくす効果）および鮮鋭度強調効果の異なる４つのフィルタ処理部、第１の反ぼけフィルタ処理部６、第２の反ぼけフィルタ処理部７、第１の鮮鋭度強調フィルタ処理部８、および第２の鮮鋭度強調フィルタ処理部９を備えている。第１の反ぼけフィルタ処理部６は、標準より強い反ぼけフィルタ処理を行うものであり、第２の反ぼけフィルタ処理部７は、標準の反ぼけフィルタ処理を行うものである。第１の鮮鋭度強調フィルタ処理部８は、標準の鮮鋭度強調フィルタ処理を行うものであり、第２の鮮鋭度強調フィルタ処理部９は、標準より強い鮮鋭度強調フィルタ処理を行う
ものである。

そして、フィルタ処理部４は、フィルタ選択部５を備えている。このフィルタ選択部５が、上記第１のぼけ指標および上記第２のぼけ指標の値に基づいて、４つのフィルタ処理部、第１の反ぼけフィルタ処理部６、第２の反ぼけフィルタ処理部７、第１の鮮鋭度強調フィルタ処理部８、および第２の鮮鋭度強調フィルタ処理部９の何れかを選択し、選択したフィルタ処理部に対して入力デジタル画像ファイルを送り、フィルタ処理させるようになっている。そして、４つのフィルタ処理部の何れかでフィルタ処理されたデジタル画像ファイルが、出力デジタル画像ファイルとして出力される。

図９の画像処理装置を構成する各部、すなわち、ＤＣＴ係数読み込み部１、第１のぼけ指標計算部２、第２のぼけ指標計算部３、フィルタ選択部５、第１の反ぼけフィルタ処理部６、第２の反ぼけフィルタ処理部７、第１の鮮鋭度強調フィルタ処理部８、および第２の鮮鋭度強調フィルタ処理部９は、ＣＰＵが記憶装置に格納されたプログラムを実行し、図示しない入出力回路などの周辺回路を制御することによって実現される機能ブロックである。

次に、本発明に係る画像処理装置の他の実施の形態について、図１０に基づいて説明する。なお、説明の便宜上、図９に示す各部材と同一の機能を有する部材には、同一の符号を付記し、その説明を省略する。

図１０に示すように、本形態の画像処理装置は、ＤＣＴ係数読み込み部１に代えて、ＤＣＴ係数の２次元のヒストグラムを作成する２次元ヒストグラム作成部１１と、上記２次元のヒストグラムを１次元のヒストグラムに変換するヒストグラム変換部１２とを備え、第１のぼけ指標計算部２および第２のぼけ指標計算部３が、上記１次元のヒストグラムからぼけ指標（第１のぼけ指標または第２のぼけ指標）を計算するようになっている以外は、前記の図９に示す画像処理装置と同様である。

図１０の画像処理装置を構成する各部、すなわち、第１のぼけ指標計算部２、第２のぼけ指標計算部３、フィルタ選択部５、第１の反ぼけフィルタ処理部６、第２の反ぼけフィルタ処理部７、第１の鮮鋭度強調フィルタ処理部８、第２の鮮鋭度強調フィルタ処理部９、２次元ヒストグラム作成部１１、およびヒストグラム変換部１２は、ＣＰＵが記憶装置に格納されたプログラムを実行し、図示しない入出力回路などの周辺回路を制御することによって実現される機能ブロックである。

なお、上記実施形態では、画像処理装置を構成する各部が、「ＣＰＵなどの演算手段がＲＯＭやＲＡＭなどの記録媒体に格納されたプログラムコードを実行することで実現される機能ブロックである」場合を例にして説明したが、同様の処理を行うハードウェアで実現してもよい。また、処理の一部を行うハードウェアと、当該ハードウェアの制御や残余の処理を行うプログラムコードを実行する上記演算手段とを組み合わせても実現することもできる。さらに、上記各部材のうち、ハードウェアとして説明した部材であっても、処理の一部を行うハードウェアと、当該ハードウェアの制御や残余の処理を行うプログラムコードを実行する上記演算手段とを組み合わせても実現することもできる。なお、上記演算手段は、単体であってもよいし、装置内部のバスや種々の通信路を介して接続された複数の演算手段が共同してプログラムコードを実行してもよい。

上記演算手段によって直接実行可能なプログラムコード自体、または、後述する解凍などの処理によってプログラムコードを生成可能なデータとしてのプログラムは、当該プログラム（プログラムコードまたは上記データ）を記録媒体に格納し、当該記録媒体を配付したり、あるいは、上記プログラムを、有線または無線の通信路を介して伝送するための通信手段で送信したりして配付され、上記演算手段で実行される。

なお、通信路を介して伝送する場合、通信路を構成する各伝送媒体が、プログラムを示す信号列を伝搬し合うことによって、当該通信路を介して、上記プログラムが伝送される。また、信号列を伝送する際、送信装置が、プログラムを示す信号列により搬送波を変調することによって、上記信号列を搬送波に重畳してもよい。この場合、受信装置が搬送波を復調することによって信号列が復元される。一方、上記信号列を伝送する際、送信装置が、デジタルデータ列としての信号列をパケット分割して伝送してもよい。この場合、受信装置は、受信したパケット群を連結して、上記信号列を復元する。また、送信装置が、信号列を送信する際、時分割／周波数分割／符号分割などの方法で、信号列を他の信号列と多重化して伝送してもよい。この場合、受信装置は、多重化された信号列から、個々の信号列を抽出して復元する。いずれの場合であっても、通信路を介してプログラムを伝送できれば、同様の効果が得られる。

ここで、プログラムを配付する際の記録媒体は、取外し可能である方が好ましいが、プログラムを配付した後の記録媒体は、取外し可能か否かを問わない。また、上記記録媒体は、プログラムが記憶されていれば、書換え（書き込み）可能か否か、揮発性か否か、記録方法および形状を問わない。記録媒体の一例として、磁気テープやカセットテープなどのテープ、あるいは、フロッピー（登録商標）ディスクやハードディスクなどの磁気ディスク、または、ＣＤ−ＲＯＭや光磁気ディスク（ＭＯ）、ミニディスク（ＭＤ）やデジタルビデオディスク（ＤＶＤ）などのディスクが挙げられる。また、記録媒体は、ＩＣカードや光カードのようなカード、あるいは、マスクＲＯＭやＥＰＲＯＭ、ＥＥＰＲＯＭまたはフラッシュＲＯＭなどのような半導体メモリであってもよい。あるいは、ＣＰＵなどの演算手段内に形成されたメモリであってもよい。

なお、上記プログラムコードは、上記各処理の全手順を上記演算手段へ指示するコードであってもよいし、所定の手順で呼び出すことで、上記各処理の一部または全部を実行可能な基本プログラム（例えば、オペレーティングシステムやライブラリなど）が既に存在していれば、当該基本プログラムの呼び出しを上記演算手段へ指示するコードやポインタなどで、上記全手順の一部または全部を置き換えてもよい。

また、上記記録媒体にプログラムを格納する際の形式は、例えば、実メモリに配置した状態のように、演算手段がアクセスして実行可能な格納形式であってもよいし、実メモリに配置する前で、演算手段が常時アクセス可能なローカルな記録媒体（例えば、実メモリやハードディスクなど）にインストールした後の格納形式、あるいは、ネットワークや搬送可能な記録媒体などから上記ローカルな記録媒体にインストールする前の格納形式などであってもよい。また、プログラムは、コンパイル後のオブジェクトコードに限るものではなく、ソースコードや、インタプリトまたはコンパイルの途中で生成される中間コードとして格納されていてもよい。いずれの場合であっても、圧縮された情報の解凍、符号化された情報の復号、インタプリト、コンパイル、リンク、または、実メモリへの配置などの処理、あるいは、各処理の組み合わせによって、上記演算手段が実行可能な形式に変換可能であれば、プログラムを記録媒体に格納する際の形式に拘わらず、同様の効果を得ることができる。

以上の明細書で使用した用語および表現は、あくまでも本発明を説明するものであり、本発明を制限するものではなく、そのような用語および表現の使用において図示および説明した特徴の等価物を排除する意図はなく、本発明の範囲は請求項によって定義および限定されるものと解釈されるべきである。

本発明に係る画像処理方法および画像処理装置、並びに画像処理プログラムおよびそれを記録した記録媒体は、ＪＰＥＧ形式などのようなデジタル画像をＤＣＴ係数の配列として格納する形式のデジタル画像ファイルに対して画像のぼけを低減する画像処理を効率的に行うことを可能にする。したがって、本発明は、デジタルカメラ等のような画像入力装置、インクジェットプリンタや熱転写プリンタ等のような画像出力装置、画像処理プログラムをインストールしたコンピュータ・システム等のような画像処理システム等の用途において有用である。

ぼけ画像および非ぼけ画像のそれぞれのパワースペクトルを示すグラフである。 典型的なぼけ予測・低減システムを示すブロック図である。 増加する垂直周波数および水平周波数の格子に配置される典型的な２ＤのＤＣＴ係数ブロックを示す。 非ぼけ画像の１ブロック内におけるＤＣＴ係数の標準偏差を示す表である。 ぼけ画像の１ブロック内におけるＤＣＴ係数の標準偏差を示す表である。 画像内におけるＤＣＴ係数の標準偏差を垂直ＤＣＴ成分および水平ＤＣＴ成分の関数として表す表面を示す。 本発明に係る方法の一例を示す。 本発明に関して使用可能なマッピング関数の一例を示す。 本発明の実施の一形態に係る画像処理装置の構成を示すブロック図である。 本発明の他の実施の形態に係る画像処理装置の構成を示すブロック図である。

符号の説明

１ ＤＣＴ係数読み込み部
２ 第１のぼけ指標計算部（ぼけ指標計算部）
３ 第２のぼけ指標計算部（ぼけ指標計算部）
４ フィルタ処理部
５ フィルタ選択部
６ 第１の反ぼけフィルタ処理部
７ 第２の反ぼけフィルタ処理部
８ 第１の鮮鋭度強調フィルタ処理部
９ 第２の鮮鋭度強調フィルタ処理部
１１ ２次元ヒストグラム作成部
１２ ヒストグラム変換部
２０ ぼけ予測・低減システム（画像処理装置）
２１ ぼけ予測部
２２ ぼけ低減部

Claims (21)

1. デジタルファイル中にＤＣＴ係数の配列として格納されたデジタル画像であり、上記ＤＣＴ係数の配列が周波数成分の格子を表す複数のブロック内に配列され、上記各ブロックが画像の部分を表すデジタル画像に対し、デジタル画像のぼけを検出および修正する画像処理方法であって、
   （ａ）上記ＤＣＴ係数の配列で構成される各ブロックから少なくとも２つのＤＣＴ係数を読み込むように複数のＤＣＴ係数を読み込むステップと、
   （ｂ）上記複数のＤＣＴ係数から、全てのブロックの中で最小のぼけ量を持つブロックのぼけ量を数値として表す第１のぼけ指標を計算するステップと、
   （ｃ）上記複数のＤＣＴ係数から、上記デジタル画像の全体としてのぼけ量を数値として表す第２のぼけ指標を計算するステップと、
   （ｄ）上記デジタル画像に対し、上記第１のぼけ指標および上記第２のぼけ指標の値に基づいて選択的にフィルタを適用するステップとを含むことを特徴とする画像処理方法。
2. 上記選択的にフィルタを適用するステップでは、複数の異なる種類のフィルタの中から１つのフィルタを選択してそのフィルタを上記デジタル画像に適用することを特徴とする請求項１記載の画像処理方法。
3. 格子内の代表的な周波数成分を選択するステップをさらに含み、
   上記第１のぼけ指標は、画像内の全てのブロックにおける上記周波数成分でのＤＣＴ係数の中の最大値であることを特徴とする請求項１記載の画像処理方法。
4. 格子内の２つの代表的な周波数成分を選択し、これら２つの代表的な周波数成分におけるＤＣＴ係数の標準偏差を画像全体で計算するステップをさらに含み、
   上記第２のぼけ指標は、これら標準偏差の比であることを特徴とする請求項１記載の画像処理方法。
5. （ｅ）周波数成分の格子を、周波数が増加する複数領域に分割するステップと、
   （ｆ）各ブロックの各領域におけるＤＣＴ係数をまとめるステップとをさらに含み、
   第１のぼけ指標を計算するステップおよび第２のぼけ指標を計算するステップでは、上記まとめられたＤＣＴ係数から、上記第１のぼけ指標および第２のぼけ指標を計算することを特徴とする請求項１記載の画像処理方法。
6. 上記各領域におけるＤＣＴ係数をまとめるステップは、上記各領域内におけるＤＣＴ係数の平均値を得るステップを含むことを特徴とする請求項５記載の画像処理方法。
7. 第１のぼけ指標が、上記複数領域の中から選択された１つの領域内における最大の平均ＤＣＴ係数であることを特徴とする請求項６記載の画像処理方法。
8. 複数領域の中から選択された２つの領域の各々においてＤＣＴ係数の標準偏差を得るステップを含むことを特徴とする請求項６記載の画像処理方法。
9. 第２のぼけ指標が、得られた標準偏差同士の比であることを特徴とする請求項８記載の画像処理方法。
10. デジタルファイル中にＤＣＴ係数の配列として格納されたデジタル画像であり、上記ＤＣＴ係数の配列が水平周波数成分および垂直周波数成分の格子を表す複数のブロック内に配列され、上記各ブロックが画像の部分を表すデジタル画像に対し、デジタル画像のぼけを検出および修正する画像処理方法であって、
    （ａ）ＤＣＴ係数の２次元のヒストグラムを得るステップと、
    （ｂ）上記２次元のヒストグラムを１次元のヒストグラムにまとめるステップと、
    （ｃ）上記１次元のヒストグラムから少なくとも１つのぼけ指標を計算するステップと、
    （ｄ）上記デジタル画像に対し、上記少なくとも１つのぼけ指標に基づいて選択的にフィルタを適用するステップとを含むことを特徴とする画像処理方法。
11. ２次元のヒストグラムを１次元のヒストグラムにまとめるステップでは、２次元ヒストグラムを周波数が概ね増加する複数領域に分割して、各ブロックの各領域内のＤＣＴ係数を平均化することによって１次元のヒストグラムを得ることを特徴とする請求項１０記載の画像処理方法。
12. 上記ぼけ指標が、上記複数領域から選択された１つの領域内における最大の平均ＤＣＴ係数であることを特徴とする請求項１１記載の画像処理方法。
13. 上記ぼけ指標が、上記複数領域から選択された２つの領域間における平均ＤＣＴ係数の比であることを特徴とする請求項１１記載の画像処理方法。
14. 第１のぼけ指標および第２のぼけ指標を計算するステップをさらに含むことを特徴とする請求項１１記載の画像処理方法。
15. 第１のぼけ指標および第２のぼけ指標を用いることを特徴とする請求項１１記載の画像処理方法。
16. 上記第１のぼけ指標が、上記複数領域から選択された１つの領域内における最大の平均ＤＣＴ係数であり、
    上記第２のぼけ指標が、上記複数領域から選択された２つの領域間における平均ＤＣＴ値の比であることを特徴とする請求項１４または１５記載の画像処理方法。
17. 上記フィルタが、異なる複数種類のフィルタの中から第１のぼけ指標および第２のぼけ指標に基づいて選択されることを特徴とする請求項１６記載の画像処理方法。
18. 請求項１ないし１７のいずれか１項に記載の画像処理方法における各ステップをコンピュータに実行させるための画像処理プログラム。
19. 請求項１８記載の画像処理プログラムを記録したコンピュータ読み取り可能な記録媒体。
20. デジタルファイル中にＤＣＴ係数の配列として格納されたデジタル画像であり、上記ＤＣＴ係数の配列が周波数成分の格子を表す複数のブロック内に配列され、上記各ブロックが画像の部分を表すデジタル画像に対し、デジタル画像のぼけを検出および修正する画像処理装置であって、
    上記ＤＣＴ係数の配列から、各ブロックから少なくとも２つのＤＣＴ係数を読み込むように複数のＤＣＴ係数を読み込むＤＣＴ係数読み込み部と、
    上記複数のＤＣＴ係数から、全てのブロックの中で最小のぼけ量を持つブロックのぼけ量を数値として表す第１のぼけ指標を計算する第１のぼけ指標計算部と、
    上記複数のＤＣＴ係数から、上記デジタル画像の全体としてのぼけ量を数値として表す第２のぼけ指標を計算する第２のぼけ指標計算部と、
    上記デジタル画像に対し、上記第１のぼけ指標および上記第２のぼけ指標の値に基づいて選択的にフィルタを適用するフィルタ処理部とを備えることを特徴とする画像処理装置。
21. デジタルファイル中にＤＣＴ係数の配列として格納されたデジタル画像であり、上記ＤＣＴ係数の配列が水平周波数成分および垂直周波数成分の格子を表す複数のブロック内に配列され、上記各ブロックが画像の部分を表すデジタル画像に対し、デジタル画像のぼけを検出および修正する画像処理装置であって、
    ＤＣＴ係数の２次元のヒストグラムを作成する２次元ヒストグラム作成部と、
    上記２次元のヒストグラムを１次元のヒストグラムに変換するヒストグラム変換部と、
    上記１次元のヒストグラムから少なくとも１つのぼけ指標を計算するぼけ指標計算部と、
    上記デジタル画像に対し、上記少なくとも１つのぼけ指標に基づいて選択的にフィルタを適用するフィルタ処理部とを備えることを特徴とする画像処理装置。

Images