(57)【要約】【課題】デジタルカメラ等の画像撮影装置において撮影画像を自動的に補正可能にする。【解決手段】デジタルカメラ１は、撮影画像を自動的に評価（露出条件評価、コントラスト評価、ボケ、焦点ボケの評価）する画像評価処理を実行する。また、撮影モードが指定されると、撮影画像を記録する記録媒体の残り容量が予め設定された容量以下である場合、記録媒体に記録された画像のうち、画像評価処理により算出された評価値が所定値以下である撮影画像の画質を補正する補正処理を実行する。そして、画質補正が施された画像の評価処理を再度行い、この画像評価処理により算出された評価値が所定値以下の画像を、記録媒体から削除する候補として決定し、決定された削除候補の画像の中から、ユーザによって指定された画像を削除する。【選択図】図３
2006050494.tif 000002
【特許請求の範囲】
【請求項１】
被写体を撮影して撮影画像を取得する撮像手段と、
前記取得された撮影画像を記録媒体に記録する記録手段と、
前記取得された撮影画像の画質を評価する第１の評価手段と、
前記記録媒体に記録された撮影画像のうち、前記第１の評価手段により算出された評価結果が所定の条件を満たす撮影画像の画質を補正する補正手段と、
前記補正手段により画質が補正された撮影画像の画質を評価する第２の評価手段と、
前記第２の評価手段による評価結果が所定の条件を満たす撮影画像を前記記録媒体から削除する候補として決定する決定手段と、
を備えることを特徴とする画像撮影装置。
【請求項２】
前記記録媒体の残り容量が予め設定された容量以下であるか否かを判定する判定手段を備え、
前記決定手段は、前記判定手段により、前記記録媒体の残り容量が予め設定された容量以下であると判定された場合に、前記第２の評価手段による評価結果が所定の条件を満たす撮影画像を前記記録媒体から削除する候補として決定することを特徴とする請求項１に記載の画像撮影装置。
【請求項３】
前記第１の評価手段及び第２の評価手段は、撮影画像の輝度又は濃度のヒストグラム分布に基づいて、当該撮影画像の画質を評価することを特徴とする請求項１又は２に記載の画像撮影装置。
【請求項４】
前記第１の評価手段及び第２の評価手段は、撮影画像の周波数特性に基づいて、当該撮影画像の画質を評価することを特徴とする請求項１又は２に記載の画像撮影装置。
【請求項５】
前記第１の評価手段及び第２の評価手段は、撮影画像のブレ又は焦点ボケの度合いに基づいて当該撮影画像の画質を評価することを特徴とする請求項１又は２に記載の画像撮影装置。
【請求項６】
前記補正手段は、撮影画像に対してガンマ補正を施すことを特徴とする請求項１〜５の何れか一項に記載の画像撮影装置。
【請求項７】
前記補正手段は、撮影画像の輝度又は濃度のヒストグラム分布を所定の変換式により変換することにより、当該撮影画像を補正することを特徴とする請求項１〜５の何れか一項に記載の画像撮影装置。
【請求項８】
前記補正手段は、撮影画像の輝度又は濃度のヒストグラム分布の平坦化を行うことにより、当該撮影画像を補正することを特徴とする請求項１〜５の何れか一項に記載の画像撮影装置。
【請求項９】
前記補正手段は、画像を平滑化するための空間フィルタを用いて撮影画像を補正することを特徴とする請求項１〜５の何れか一項に記載の画像撮影装置。
【請求項１０】
前記補正手段は、所定の微分フィルタを用いて撮影画像を鮮鋭化する処理を行うことを特徴とする請求項１〜５の何れか一項に記載の画像撮影装置。
【請求項１１】
前記補正手段は、撮影画像を周波数変換し、周波数変換された撮影画像に対して、鮮鋭化処理又はぼかし処理を施すことを特徴とする請求項１〜５の何れか一項に記載の画像撮影装置。
【請求項１２】
前記補正手段は、点像分布関数又は線像分布関数を用いて、撮影画像の直線ブレ又は焦点ボケを補正することを特徴とする請求項１〜５の何れか一項に記載の画像撮影装置。
【請求項１３】
画像劣化が生じる方式で見本画像を撮影する劣化画像撮像手段と、
前記劣化画像撮像手段により得られた画像の劣化成分を算出する算出手段と、
画像劣化の度合い別に、前記算出手段により算出された画像の劣化成分の情報を記憶する記憶手段と、を備え、
前記撮像手段により取得された撮影画像の評価結果が所定の条件を満たす場合、前記補正手段は、当該撮影画像に対応する劣化成分を前記記憶手段から読み出し、その読み出された劣化成分に基づいて当該撮影画像を補正することを特徴とする請求項１〜１２の何れか一項に記載の画像撮影装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像撮影装置に関する。
【背景技術】
【０００２】
静止画像や動画像を撮像し、メモリーカード等の記録媒体に記録するデジタルカメラでは、フィルム代や現像代などを気にすることなく撮影画像を記録することができる。最近では、記録媒体の価格が下がり、記録容量が大幅に増え、多数の画像が記録可能になっている（例えば、特許文献１参照。）。記録媒体の容量が一杯になった場合は、記録媒体に記録された画像データをＰＣ（Personal Computer）等に転送してから記録媒体を空にするか、記録済みの画像の中から、ユーザが画像の良し悪しを判断して不要な画像を消去するなどして、記録媒体の空き領域が確保されている。
【特許文献１】特許第２７８７７８１号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
しかしながら、従来、デジタルカメラの記録媒体に記録された画像の中から不要な画像を削除する場合、ユーザが記録済み画像の画質等を目視で評価し、その評価結果に基づいて不要な画像を決めて削除していたため、実際は画質が良い画像であっても削除される可能性があった。
【０００４】
本発明の課題は、デジタルカメラ等の画像撮影装置において撮影画像を自動的に補正することである。
【課題を解決するための手段】
【０００５】
上記課題を解決するため、請求項１に記載の画像撮影装置は、被写体を撮影して撮影画像を取得する撮像手段と、前記取得された撮影画像を記録媒体に記録する記録手段と、前記取得された撮影画像の画質を評価する第１の評価手段と、前記記録媒体に記録された撮影画像のうち、前記第１の評価手段により算出された評価結果が所定の条件を満たす撮影画像の画質を補正する補正手段と、前記補正手段により画質が補正された撮影画像の画質を評価する第２の評価手段と、前記第２の評価手段による評価結果が所定の条件を満たす撮影画像を前記記録媒体から削除する候補として決定する決定手段と、を備えることを特徴としている。
【０００６】
請求項２に記載の発明は、請求項１に記載の画像撮影装置において、前記記録媒体の残り容量が予め設定された容量以下であるか否かを判定する判定手段を備え、前記決定手段は、前記判定手段により、前記記録媒体の残り容量が予め設定された容量以下であると判定された場合に、前記第２の評価手段による評価結果が所定の条件を満たす撮影画像を前記記録媒体から削除する候補として決定することを特徴としている。
【０００７】
請求項３に記載の発明は、請求項１又は２に記載の画像撮影装置において、前記第１の評価手段及び第２の評価手段は、撮影画像の輝度又は濃度のヒストグラム分布に基づいて、当該撮影画像の画質を評価することを特徴としている。
【０００８】
請求項４に記載の発明は、請求項１又は２に記載の画像撮影装置において、前記第１の評価手段及び第２の評価手段は、撮影画像の周波数特性に基づいて、当該撮影画像の画質を評価することを特徴としている。
【０００９】
請求項５に記載の発明は、請求項１又は２に記載の画像撮影装置において、前記第１の評価手段及び第２の評価手段は、撮影画像のブレ又は焦点ボケの度合いに基づいて当該撮影画像の画質を評価することを特徴としている。
【００１０】
請求項６に記載の発明は、請求項１〜５の何れか一項に記載の画像撮影装置において、前記補正手段は、撮影画像に対してガンマ補正を施すことを特徴としている。
【００１１】
請求項７に記載の発明は、請求項１〜５の何れか一項に記載の画像撮影装置において、前記補正手段は、撮影画像の輝度又は濃度のヒストグラム分布を所定の変換式により変換することにより、当該撮影画像を補正することを特徴としている。
【００１２】
請求項８に記載の発明は、請求項１〜５の何れか一項に記載の画像撮影装置において、前記補正手段は、撮影画像の輝度又は濃度のヒストグラム分布の平坦化を行うことにより、当該撮影画像を補正することを特徴としている。
【００１３】
請求項９に記載の発明は、請求項１〜５の何れか一項に記載の画像撮影装置において、前記補正手段は、画像を平滑化するための空間フィルタを用いて撮影画像を補正することを特徴としている。
【００１４】
請求項１０に記載の発明は、請求項１〜５の何れか一項に記載の画像撮影装置において、前記補正手段は、所定の微分フィルタを用いて撮影画像を鮮鋭化する処理を行うことを特徴としている。
【００１５】
請求項１１に記載の発明は、請求項１〜５の何れか一項に記載の画像撮影装置において、前記補正手段は、撮影画像を周波数変換し、周波数変換された撮影画像に対して、鮮鋭化処理又はぼかし処理を施すことを特徴としている。
【００１６】
請求項１２に記載の発明は、請求項１〜５の何れか一項に記載の画像撮影装置において、前記補正手段は、点像分布関数又は線像分布関数を用いて、撮影画像の直線ブレ又は焦点ボケを補正することを特徴としている。
【００１７】
請求項１３に記載の発明は、請求項１〜１２の何れか一項に記載の画像撮影装置において、画像劣化が生じる方式で見本画像を撮影する劣化画像撮像手段と、前記劣化画像撮像手段により得られた画像の劣化成分を算出する算出手段と、画像劣化の度合い別に、前記算出手段により算出された画像の劣化成分の情報を記憶する記憶手段と、を備え、前記撮像手段により取得された撮影画像の評価結果が所定の条件を満たす場合、前記補正手段は、当該撮影画像に対応する劣化成分を前記記憶手段から読み出し、その読み出された劣化成分に基づいて当該撮影画像を補正することを特徴としている。
【発明の効果】
【００１８】
本発明によれば、撮影画像を自動的に評価し、評価結果が所定の条件を満たす画像（例えば、評価値が所定値以下の画像）を自動的に補正し、補正処理後に再度、当該画像を評価し、その評価結果が所定の条件を満たす画像（例えば、評価値が所定値以下の画像）のみを削除可能としたことにより、必要以上に画像を削除することなくなるとともに、記録媒体の空き領域を確実に確保することができる。
【００１９】
また、予め、画像劣化が生じる方式で見本画像を撮影して、画像劣化の度合い別に劣化成分を算出して記憶し、実際に撮影された画像の評価結果が所定の条件を満たす場合（例えば、評価値が所定値以下である場合）、記憶手段に記憶された劣化成分に基づいて画像補正を行うことにより、個々の画像撮影装置の性能に合わせた画像補正を容易に行うことができる。また、画像補正の際に劣化成分を算出する必要がないため、画像補正処理を高速に行うことができる。
【発明を実施するための最良の形態】
【００２０】
以下、図面を参照して、本発明の実施形態を詳細に説明する。
まず、本実施形態における構成について説明する。
【００２１】
図１に、本発明の画像撮影装置を適用した実施形態に係るデジタルカメラ１の回路構成を示す。デジタルカメラ１は、図１に示すように、モータ（Ｍ）２１、レンズ光学系２２、ＣＣＤ２３、ＴＧ（タイミング発生器）２４、垂直ドライバ２５、Ｓ／Ｈ（サンプルホールド回路）２６、Ａ／Ｄ変換器２７、カラープロセス回路２８、ＤＭＡ（Direct Memory Access）コントローラ２９、ＤＲＡＭ（Dynamic RAM） I／Ｆ（ＤＲＡＭインターフェース）３０、ＤＲＡＭ３１、制御部３２、ＶＲＡＭ（Video RAM）コントローラ３３、ＶＲＡＭ３４、デジタルビデオエンコーダ３５、キー入力部３６、ＪＰＥＧ（Joint Photograph coding Experts Group）回路３７、フラッシュメモリ３８、通信部３９、手ブレ検出センサ４０、表示部１３により構成される。
【００２２】
モータ（Ｍ）２１は、制御部３２から入力される制御信号に従って駆動し、レンズ光学系２２の絞り位置を移動させる。
【００２３】
ＣＣＤ２３は、撮像素子が平面状（２次元）に配列された構造を有し、光入力を電気信号に変換して蓄積する光電変換部、蓄積された電荷を読み出す走査部、電気信号として出力する出力部により構成され、撮影モードでのモニタリング状態において、タイミング発生器（ＴＧ）２４、垂直ドライバ２５によって駆動される。
【００２４】
ＣＣＤ２３の出力部から出力された電気信号は、アナログ値の信号の状態でＲＧＢの各原色成分毎に適宜ゲイン調整され、Ｓ／Ｈ（サンプルホールド回路）２６は、このゲイン調整された信号をサンプルホールドする。Ａ／Ｄ変換器２７は、Ｓ／Ｈ（サンプルホールド回路）２６から出力されたアナログ信号をデジタル信号に変換し、カラープロセス回路２８に出力する。
【００２５】
カラープロセス回路２８は、Ａ／Ｄ変換器２７から出力されたデジタル信号に対して画素補間処理及びγ補正処理を含むカラープロセス処理を施し、デジタル値の輝度信号Ｙ及び色差信号Ｃｂ、Ｃｒを生成し、ＤＭＡコントローラ２９に出力する。
【００２６】
ＤＭＡコントローラ２９は、カラープロセス回路２８から出力された輝度信号Ｙ及び色差信号Ｃｂ、Ｃｒを、同じくカラープロセス回路２８から出力された複合同期信号、メモリ書込みイネーブル信号及びクロック信号を用いて一度ＤＭＡコントローラ２９内部のバッファに書込み、ＤＲＡＭ Ｉ／Ｆ３０を介してバッファメモリとして使用されるＤＲＡＭ３１にＤＭＡ転送を行う。
【００２７】
制御部３２は、ＣＰＵ（Central Processing Unit）、ＣＰＵで実行される動作プログラムを固定的に記憶したＲＯＭ（Read Only Memory）、ワークメモリとして使用されるＲＡＭ（Random Access Memory）等により構成され、デジタルカメラ１の各部の動作を制御する。以下、制御部３２による制御動作について説明する。
【００２８】
制御部３２は、撮影された画像の画質を評価する画像評価処理を実行する。また、制御部３２は、撮影モードが指定されると、フラッシュメモリ３８（記録媒体）の残り容量が予め設定された容量以下であるか否かを判定し、記録媒体の残り容量が予め設定され容量以下であると判定された場合、記録媒体に記録された画像のうち、画像評価処理により算出された評価値が所定値以下である撮影画像の画質を補正する補正処理を実行する。そして、制御部３２は、画質補正が施された画像の評価処理を再度行い、この画像評価処理により算出された評価値が所定値以下の画像を、記録媒体から削除する候補として決定する。そして、制御部３２は、削除候補の画像の中から、キー入力部３６の操作により指定された画像を記録媒体から削除する処理を行う。以下の、実施形態において、記録媒体に記録済みのデータを消去する場合だけでなく、記録済みのデータに新しいデータを上書きする場合も、記録済みのデータが「削除」されたものとみなすことにする。
【００２９】
画像評価の方法としては、撮影画像の輝度又は濃度のヒストグラム分布から評価する方法と（画像評価方法１：図７及び図８参照）、撮影画像の周波数特性に基づいて評価する方法と（画像評価方法２：図９〜図１１参照）、撮影画像のブレ又は焦点ボケの度合いに基づいて評価する方法（画像評価方法３：図１２〜図１６）等がある。本実施形態の画像評価処理では、これら複数の画像評価方法の全てが適用される場合を示すが、これら画像評価方法１〜３のうちの何れか一つ又は２つ以上を組み合わせてもよい。
【００３０】
撮影画像の画質を補正する方法としては、撮影画像にガンマ補正を施す方法と（画像補正方法１：図１８参照）、撮影画像の輝度又は濃度のヒストグラム分布を所定の変換式による変換する方法と（画像補正方法２：図１９及び図２０参照）、撮影画像の輝度又は濃度のヒストグラム分布を平坦化する方法と（画像補正方法３：図２１及び図２２参照）、撮影画像を平滑化する方法と（画像補正方法４：図２３参照）、微分フィルタを用いて撮影画像を鮮鋭化する方法と（画像補正方法５：図２４及び図２５参照）、周波数変換された撮影画像に鮮鋭化処理を施す方法と（画像補正方法６）、ＰＳＦ（Point Spread Function）、ＬＰＦ（Line Spread Function）等を用いて撮影画像の直線ブレ又は焦点ボケを補正する方法（画像補正方法７）等がある。本実施形態では、画像補正方法１〜７のうちの何れか一つ又は２つ以上を組み合わせた画像補正方法を適用することができる。
【００３１】
ＶＲＡＭコントローラ３３は、制御部３２の制御により、ＤＲＡＭ３１に一時記憶されている画像データ（輝度信号Ｙ及び色差信号Ｃｂ、Ｃｒ）を、ＤＲＡＭ３１からＤＲＡＭ Ｉ／Ｆ３０を介して読出し、ＶＲＡＭ３４に書込む。また、再生モード時、制御部３２により伸長された画像データをＶＲＡＭ３４に展開して記憶させる。
【００３２】
デジタルビデオエンコーダ３５は、ＶＲＡＭ３４に記憶されている画像データを、ＶＲＡＭ３４からＶＲＡＭコントローラ３３を介して定期的に読出し、これらのデータを元にビデオ信号を生成して表示部１３に出力する。
【００３３】
キー入力部３６は、電源キー、シャッターキー、モードスイッチ、メニューキー及び十字キー等の各種キーにより構成され、それらのキー操作に応じた信号を制御部３２に出力する。キー入力部３６のシャッターキーは、２段階のストロークで動作し、一般的に「半押し」と表現されている第１段階の操作状態でＡＥ（自動露光）処理やＡＦ（オートフォーカス）処理を始めとする撮影の準備を行い、一般的に「全押し」と表現されている、より強く押圧操作した第２段階の操作状態で撮影を実行する。
【００３４】
ＪＰＥＧ回路３７は、ＣＣＤ２３から取り込まれ、ＤＲＡＭ３１に一時的に記憶されている画像信号を、ＡＤＣＴ（Adaptive Discrete Cosine Transform：適応離散コサイン変換）、エントロピー符号化方式であるハフマン符号化などの処理によりデータ圧縮する。
【００３５】
フラッシュメモリ３８は、デジタルカメラ１の記録媒体として着脱自在に装着されるメモリーカード内に封入された不揮発性メモリであり、ＪＥＰＧ回路３７により符号化された画像信号が記録される。なお、フラッシュメモリ３８は、デジタルカメラ１に内蔵されていてもよい。
【００３６】
通信部３９は、外部機器とケーブルを介したデータ通信を行うためのインターフェースを備える。例えば、通信部３９は、ＩＥＥＥ１３９４やＵＳＢ（Universal Serial Bus）等により通信を行うためのインターフェースを備える。なお、外部機器と無線通信手段により通信可能としてもよい。
【００３７】
手ブレ検出センサ４０は、加速度センサ等であり、撮影時の手ブレ量（ブレ方向（θ）、ブレの距離（Ｌ）等、後述）を検出する。表示部１３は、撮影モード時にはモニタ表示部（電子ファインダ）として機能し、デジタルビデオエンコーダ３５からのビデオ信号に基づいた表示を行うことで、その時点でＶＲＡＭコントローラ３３から取込んでいる画像情報に基づくスルー画像をリアルタイムに表示する。
【００３８】
次に、本実施形態における動作について説明する。
まず、図２のフローチャートを参照して、デジタルカメラ１において実行される画像撮影・削除処理について説明する。
【００３９】
まず、表示部１３に、図４（ａ）に示すようなメニュー画面が表示される（ステップＳ１）。図４（ａ）に示すメニュー画面には、撮影モードを指定するための「撮影モード」ボタン、撮影画像の再生を指定するための「再生モード」ボタン、撮影画像の削除を指定するための「画像の削除」ボタン、その他の各種の設定を指定するための「各種の設定」ボタン等が表示されている。
【００４０】
メニュー画面において、「画像の削除」ボタンが指定されると（ステップＳ３；ＹＥＳ）、図４（ｂ）及び図５（ａ）に示すような削除モード画面が表示され、キー入力部３６の操作により、画像の自動削除、画質補正等の各種の設定処理が行われ（ステップＳ４）、設定処理後、ステップＳ１に戻る。メニュー画面において、「撮影モード」ボタン及び「画像の削除」ボタン以外のボタンが指定された場合（ステップＳ３；ＮＯ）、その指定されたモードに対応する処理が行われ（ステップＳ５）、ステップＳ１に戻る。
【００４１】
削除モード画面において、「自動削除の設定」ボタンが指定されると、図５（ｂ）に示すような自動削除設定画面が表示される。自動削除設定画面には、図５（ｂ）に示すように、設定項目として、撮影画像の自動削除機能を有効（ＯＮ）にするか無効（ＯＦＦ）にするかを設定するための「自動削除機能」欄と、画像の削除方法を設定するための「削除方法」欄、画像の評価方法を設定するための「評価方法」欄と、各評価機能を有効（ＯＮ）にするか無効（ＯＦＦ）にするかを設定するための「評価項目」欄（例えば、「露出条件評価」項目、「コントラスト評価」項目、「ボケ・ブレ評価」項目）が表示されている。また、各評価項目欄には、評価項目の基準値を設定するための「基準値」欄が表示されている。ここで、評価項目の基準値とは、補正対象、削除対象の基準となる数値であり、この基準値に基づいて画像の評価値が決定される。
【００４２】
図５（ａ）に示す削除モード画面において、「画質の補正」ボタンが指定されると、図５（ｃ）に示すような補正項目画面が表示される。補正項目画面には、図５（ｃ）に示すように、画質の補正項目として、「露出（濃度レベル）の補正」ボタン、「トーンカーブ（コントラスト）修正」ボタン、「輪郭エッジの強調（シャープ）」ボタン、「ボケ・ブレの補正」ボタン、「ノイズの除去」ボタン、「色調（明度/色相/彩度）の修正」ボタンが表示されている。
【００４３】
図５（ａ）に示す削除モード画面において、「画質の自動補正の設定」ボタンが指定されると、図５（ｅ）に示すような画質の自動補正設定画面が表示される。図５（ｅ）に示す自動補正設定画面には、撮影画像の自動補正機能、露出条件の補正機能、輪郭エッジを強調する機能、ボケ、ブレの補正機能、ノイズの除去機能の各機能を有効（ＯＮ）にするか無効（ＯＦＦ）にするかを設定するための「自動補正機能」欄、「露出条件の補正」欄、「輪郭エッジ強調」欄、「ボケ・ブレの補正」欄、「ノイズの除去」欄が表示されている。
【００４４】
図５（ｃ）に示す補正項目画面において、「露出（濃度レベル）の補正」ボタンが指定されると、図５（ｄ）に示すような露出（濃度レベル）の補正画面が表示され、キー入力部３６の操作により、ＲＧＢの色毎に、撮影画像の濃度レベルが設定される。また、図５（ｃ）に示す補正項目画面において、「トーンカーブ（コントラスト）修正」ボタンが指定されると、図５（ｆ）に示すようなトーンカーブ修正画面が表示され、キー入力部３６の操作により、撮影画像のコントラストが設定される。
【００４５】
図４（ａ）に示すメニュー画面において、「撮影モード」ボタンが指定された場合（ステップＳ２；ＹＥＳ）、デジタルカメラ１のレンズ位置等の初期設定が行われ（ステップＳ６）、記録媒体（フラッシュメモリ３８）の残り容量が設定容量以下であるか否かが判定される（ステップＳ７）。
【００４６】
ステップＳ７において、記録媒体の残り容量が設定容量以下であると判定された場合（ステップＳ７；ＹＥＳ）、各撮影画像の評価値に基づいて記録媒体から画像データを削除する画像削除処理が行われる（ステップＳ８）。画像削除処理の詳細は、後に図３を参照して説明する。
【００４７】
ステップＳ７において、記録媒体の残り容量が設定容量より大きいと判定された場合（ステップＳ７；ＮＯ）又はステップＳ８の画像削除処理の終了後、測光処理、露出条件の設定等が行われる（ステップＳ９）。次いで、ズーム処理が行われ（ステップＳ１０）、オートフォーカス処理が行われる（ステップＳ１１）。
【００４８】
表示部１３に、図６（ａ）に示すような撮影スルー画像が表示され（ステップＳ１２）、キー入力部３６のシャッターキーの押下により撮影が指示されると（ステップＳ１３；ＹＥＳ）、撮影処理が行われ（ステップＳ１４）、撮影画像データが取得される。撮影処理後、撮影処理で取得された撮影画像データの画像評価処理が行われる（ステップＳ１５）。画像評価処理の詳細は、後に図７〜図１６を参照して説明する。シャッターキーが押されなかった場合には（ステップＳ１３；ＮＯ）、その他（撮影モード以外）のキー入力処理、表示処理等が行われ（ステップＳ２５）、本画像撮影・削除処理が終了する。
【００４９】
撮影画像データの画像評価処理が終了すると、撮影画像データとその評価値が対応付けて記録媒体に記録され（ステップＳ１６）、撮影画像データ及びその評価値が表示部１３に表示される（ステップＳ１７）。なお、ステップＳ１７において、図６（ｂ）に示すように、画像データのサイズ、撮影条件、画質評価値等の画像の詳細データを表示するようにしてもよい。また、画質評価値の各項目や総合評価値は、５点法や１０点法又は％に換算して表示するようにしてもよい。
【００５０】
次いで、撮影画像の評価値が所定値以下であるか否かが判定される（ステップＳ１８）。ここで、所定値とは、図５（ｂ）の自動削除設定画面において設定された基準値に基づいて決定され、例えば、各評価項目において設定された基準値の合計値等である。ステップＳ１８において、撮影画像の評価値が所定値より大きいと判定された場合（ステップＳ１８；ＮＯ）、その他（撮影モード以外）のキー入力処理、表示処理等が行われ（ステップＳ２５）、本画像撮影・削除処理が終了する。
【００５１】
ステップＳ１８において、撮影画像の評価値が所定値以下であると判定された場合（ステップＳ１８；ＹＥＳ）、表示部１３には、図６（ｃ）に示すような自動補正確認プロンプト画面が表示される。自動補正確認プロンプト画面には、図６（ｃ）に示すように、該当する撮影画像と、「この画像を自動補正してもよろしいですか？（Ｙ/Ｎ）」というメッセージ及び評価値が表示される。
【００５２】
図６（ｃ）に示す自動補正確認プロンプト画面において、画像の自動補正が指定されなかった場合（ステップＳ１９；ＮＯ）、ステップＳ２５に移行する。自動補正確認プロンプト画面において、画像の自動補正が指定された場合（ステップＳ１９；ＹＥＳ）、該当する撮影画像の画質の自動補正処理が行われる（ステップＳ２０）。画質の自動補正処理の詳細は、後に図１８〜図２５を参照して説明する。
【００５３】
ステップＳ２０における画質の自動補正処理が終了すると、補正済みの撮影画像に対する画像評価処理が行われる（ステップＳ２１）。画像評価処理の詳細は、後に図７〜図１６を参照して説明する。次いで、表示部１３に、図６（ｄ）に示すような補正済み画像の更新確認プロンプト画面が表示される。補正済み画像の更新確認プロンプト画面には、図６（ｄ）に示すように、該当する撮影画像と、「画質を自動補正しました。この画像に更新しますか？（Ｙ/Ｎ）」というメッセージ及び評価値が表示される。
【００５４】
補正済み画像の更新確認プロンプト画面において、補正済み画像の更新が指定されなかった場合（ステップＳ２３；ＮＯ）、ステップＳ２５に移行する。補正済み画像の更新確認プロンプト画面において、補正済み画像の更新が指定された場合（ステップＳ２３；ＹＥＳ）、その補正済み画像が符号化され、評価値とともに記録媒体に記録される（ステップＳ２４）。そして、その他（撮影モード以外）のキー入力処理、表示処理等が行われ（ステップＳ２５）、本画像撮影・削除処理が終了する。
【００５５】
なお、ステップＳ１５及びＳ２１における画像評価処理は、１枚撮影される度に行うようにしてもよいし、一連の撮影処理の終了後に全画像分の処理をまとめて行うようにしてもよい。特に、連写撮影や動画撮影では、撮影性能を保つため、後者の方法を用いるのが好ましい。
【００５６】
次に、図３のフローチャートを参照して、図２に示した画像削除処理の詳細について説明する。
【００５７】
まず、画像データの自動削除設定がＯＮであるか否か（画像データの自動削除機能が有効であるか否か）が判定される（ステップＳ３０）。ステップＳ３０において、自動削除設定がＯＦＦであると判定された場合（ステップＳ３０；ＮＯ）、表示部１３に、記録媒体の残り容量がなくなっている旨が表示されるとともに（ステップＳ３１）、マニュアル操作での画像削除を確認するための確認プロンプト画面が表示される（ステップＳ３２）。
【００５８】
表示部１３に表示された確認プロンプト画面上において、マニュアル操作での画像削除が指定されると（ステップＳ３３；ＹＥＳ）、マニュアル操作による画像の削除モードに移行し（ステップＳ３４）、本画像削除処理が終了する。ステップＳ３４では、ユーザは、記録媒体の容量を確認しながら画像を削除する。マニュアル操作での画像削除が指定されなかった場合（ステップＳ３３；ＮＯ）、本画像削除処理が終了する。
【００５９】
ステップＳ３０において、自動削除設定がＯＮであると判定された場合（ステップＳ３０；ＹＥＳ）、表示部１３に、図６（ｅ）に示すように、「メモリーの残りがありません！削除する候補を検索します。」というメッセージが表示され、記録媒体に記録済みの画像は画像評価処理済であるか否かが判定される（ステップＳ３５）。
【００６０】
ステップＳ３５において、記録済み画像が画像評価処理済でないと判定された場合（ステップＳ３５；ＮＯ）、記録済み画像の画像評価処理が行われ（ステップＳ３６）、後述のステップＳ３７に移行する。ステップＳ３６の画像評価処理の詳細は、後に図７〜図１６を参照して説明する。
【００６１】
ステップＳ３５において、記録済み画像が画像評価処理済であると判定された場合（ステップＳ３５；ＹＥＳ）又はステップＳ３６の画像評価処理が終了すると、記録媒体に記録された画像のうち、評価値が所定値以下の画像が検索され（ステップＳ３７）、評価値が最も低い順に、補正対象の画像が選択される（ステップＳ３８）。
【００６２】
補正対象の撮影画像が選択されると、表示部１３には、図６（ｃ）に示すような自動補正確認プロンプト画面が表示される（ステップＳ３９）。自動補正確認プロンプト画面には、図６（ｃ）に示すように、選択された撮影画像と、「この画像を自動補正してもよろしいですか？（Ｙ/Ｎ）」というメッセージ及び評価値が表示される。
【００６３】
図６（ｃ）に示す自動補正確認プロンプト画面において、画像の自動補正が指定されなかった場合（ステップＳ４０；ＮＯ）、後述のステップＳ４５に移行する。自動補正確認プロンプト画面において、画像の自動補正が指定された場合（ステップＳ４０；ＹＥＳ）、該当する撮影画像の画質の自動補正処理が行われる（ステップＳ４１）。画質の自動補正処理の詳細は、後に図１８〜図２５を参照して説明する。
【００６４】
画質の自動補正処理が終了すると、補正済み画像を再度評価する画像評価処理が行われる（ステップＳ４２）。ステップＳ４２の画像評価処理の詳細は、後に図７〜図１６を参照して説明する。次いで、表示部１３に、補正済み画像及びその評価値が表示され（ステップＳ４３）、その補正済み画像の評価値が所定値以下であるか否かが判定される（ステップＳ４４）。ここで、所定値とは、図５（ｂ）の自動削除設定画面において設定された基準値に基づいて決定され、例えば、各評価項目において設定された基準値の合計値等である。
【００６５】
ステップＳ４４において、補正済み画像の評価値が所定値より大きいと判定された場合（ステップＳ４４；ＮＯ）、次の補正対象の候補が検索され（ステップＳ４９）、ステップＳ３８に戻る。ステップＳ４４において、補正済み画像の評価値が所定値以下であると判定された場合（ステップＳ４４；ＹＥＳ）、その補正済み画像が記録媒体から削除する候補として決定され、表示部１３に、図６（ｆ）に示すように、削除確認プロンプト画面が表示される（ステップＳ４５）。削除確認プロンプト画面には、該当する画像と、「この画像を削除してもよろしいですか？（Ｙ/Ｎ）」というメッセージ及び評価値が表示される。
【００６６】
削除確認プロンプト画面において、該当する画像の削除が指定されなかった場合（ステップＳ４６；ＮＯ）、記録媒体に記録済みの画像の中から、次の補正対象の候補が検索され（ステップＳ４９）、ステップＳ３８に戻る。削除確認プロンプト画面において、該当する画像の削除が指定された場合（ステップＳ４６；ＹＥＳ）、該当画像が記録媒体から削除され（ステップＳ４７）、現在の記録媒体の残り容量が新規撮影で得られる画像データを保存するのに十分であるか否かが判定される（ステップＳ４８）。
【００６７】
ステップＳ４８において、記録媒体の残り容量が新規撮影の画像データを保存するのに十分であると判定された場合（ステップＳ４８；ＹＥＳ）、本画像削除処理が終了する。ステップＳ４８において、記録媒体の残り容量が新規撮影の画像データを保存するのに十分でないと判定された場合（ステップＳ４８；ＮＯ）、記録媒体に記録済みの画像の中から、次の補正対象の候補が検索され（ステップＳ４９）、ステップＳ３８に戻る。
【００６８】
なお、図３のフローチャートでは、表示部１３に、図６（ｆ）に示すような削除確認プロンプト画面を表示して、削除する画像をユーザに確認させ、ユーザから指定された画像を削除するようにしたが、削除候補として決定された画像を自動的に削除するようにしてもよい。
【００６９】
次に、図２及び図３の画像評価処理で適用される画像評価方法の詳細を説明する。
【００７０】
<画像評価方法１：輝度ヒストグラム分布に基づく評価>
まず、図７及び図８を参照して、撮影画像の輝度又は濃度のヒストグラム分布から画質を評価する方法について説明する。
【００７１】
２次元画像データの各点（ｉ、ｊ）の輝度値をｆ（ｉ、ｊ）＝ｘとし、輝度値別に画素数を集計することにより、図７に示すような輝度分布図（輝度のヒストグラム分布）Ｐ（ｘ）が作成される。一般に、輝度のヒストグラム分布が中央からはずれて右端側（輝度値大：濃度のヒストグラム分布では左端側）に偏ると、露出オーバーで、高輝度領域が飽和したり、いわゆる「白飛び」した画像になる。逆に、輝度ヒストグラム分布が左端側（濃度ヒストグラム分布では右端側）に偏ると、露出アンダーで、低輝度領域が飽和したり、「黒つぶれ」した画像になる。また、分布が中央に集中していても、その分布幅が狭くて一様に広く分布していない場合には、ダイナミックレンジが狭く、中間諧調や細かな濃淡に欠けた画像となる。これらの何れも、撮影時の露出条件が悪く、写りの良くない画像の典型である。撮影画像の輝度又は濃度分布特性の指標値や統計量から、露出条件の良し悪しを評価できる。最近のデジタルカメラの多くは、この輝度値を集計演算して、ヒストグラム（頻度分布図）をグラフ表示したりできる機種が多いので、この機能を流用すれば、比較的容易に画像評価を行うことができる。
【００７２】
本実施形態では、撮影画像の輝度又は濃度のヒストグラム分布から、輝度又は濃度の最頻値（最も頻度数が多い輝度値又は濃度値）、平均値、中央値（分布の頻度数の順位が中央にある輝度値又は濃度値）等の分布の代表値、分布の最大値、最小値、分散、標準偏差、歪み度、尖り度等を算出し、算出結果から、画像評価のための評価値が決定される。
【００７３】
輝度又は濃度のヒストグラム分布の平均値ｘM、分散（Variance）σ2、標準偏差（Standard Deviation）σ、歪み度（Skewness）Ｓkew、尖り度（Kurtosis）Ｋurtは、画素ｉの輝度値をｘi、評価対象の画像の総画素数をＮとすると、それぞれ、式（１）〜（５）のように表される。
【数１】
2006050494.tif 000003
【００７４】
ここで、分布の歪み度Ｓkewは、分布が左右対称かどうかの非対称性を示す。分布が左右対称であれば０、右に大きなはずれ値があれば（露出オーバーの場合）正、左に大きなはずれ値があれば（露出アンダーの場合）負となる。また、分布の尖り度Ｋurtは、分布が正規分布ならば０、正規分布より広がりが大きければ（ダイナミックレンジが広い場合）正、正規分布より広がりが小さければ（ダイナミックレンジが狭い場合）負となる。
【００７５】
このように、輝度又は濃度のヒストグラム分布の統計量に基づいて、分布の幅や広がり、その偏り、分布曲線の形状や特徴を判別し、画質の評価が可能となり、これらを元に画質の評価値が決定される。
【００７６】
撮影画像がカラー画像である場合、図８に示すように、ＲＧＢの各色毎に、輝度又は濃度のヒストグラム分布ＰR（ｘ）、ＰG（ｘ）、ＰB（ｘ）を作成し、色毎に、分布の平均値、分散、標準偏差、歪み度、尖り度、変動係数ＣＶ等の統計量を算出し、算出された色毎の統計量の違いや比較によって、画質の評価値を決定するようにしてもよい。ここで、変動係数ＣＶは、式（６）のように定義される。
【数２】
2006050494.tif 000004
【００７７】
<画像評価方法２：周波数特性に基づく評価>
次に、図９〜図１１を参照して、撮影画像の周波数特性から画質を評価する方法について説明する。
【００７８】
撮影画像（２次元画像データ）の各点の輝度値又は濃度値ｆ（ｘ，ｙ）を周波数変換したＦ（ｕ，ｖ）を算出し、周波数毎の振幅値や係数から、画像の周波数特性、画像全体に占める高周波成分や低周波成分の比率等から、画像のボケの程度や、鮮明度（シャープネス）、コントラスト等を評価することができる。周波数変換としては、フーリエ変換、離散フーリエ変換（ＤＦＴ：Discrete Fourier Transform）、離散コサイン変換（ＤＣＴ：Discrete Cosine Transform）、それらを高速計算する高速フーリエ変換（ＦＦＴ：Fast Fourier Transform）アルゴリズム等を用いることができる。最近のデジタルカメラは、ＪＰＥＧ等の圧縮符号化のためにＤＣＴ用のＤＳＰ（Digital Signal Processor）やソフトウェアを内蔵しているものが多いので、この機能を流用すれば、比較的容易に画像評価を行うことができる。
【００７９】
２次元画像データの輝度値又は濃度値ｆ（ｘ，ｙ）の離散フーリエ変換式は、式（７）のように表される。
【数３】
2006050494.tif 000005
【００８０】
また、２次元画像データの輝度値又は濃度値ｆ（ｘ，ｙ）の離散コサイン変換式は、式（８）のように表される。
【数４】
2006050494.tif 000006
【００８１】
このように、２次元画像ｆ（ｘ、ｙ）を周波数空間に変換して、周波数成分Ｆ（ｕ、ｖ）を算出し、全周波数成分に占める高周波成分の比率、全周波数成分に占める低周波成分の比率に基づいて画像の評価値を決定することができる。また、図９に示すように、予め、撮影画像を順次走査して１次元信号に並び替えてから周波数変換を施して、評価値を算出してもよいし、図１０に示すように、撮影画像を周波数変換したＦ（ｕ、ｖ）の振幅値を、ジグザグ走査などで１次元に並び替えてから、評価値を算出するようにしてもよい。
【００８２】
本実施形態では、全周波数成分に占める高周波成分の比率を画質評価値とし、全周波数成分に占める低周波成分の比率をボケ評価値とする。即ち、画質評価値、ボケ評価値は、それぞれ、式（９）、式（１０）のように表すことができる。
画質評価値＝（高周波成分の振幅値合計）／（全周波数成分の振幅値合計）  （９）
ボケ評価値＝（全周波数成分の振幅値合計−高周波成分の振幅値合計）／（全周波数成分の振幅値合計）    （１０）
【００８３】
また、図１１に示すように、撮影画像の輝度値又は濃度値ｆ（ｘ、ｙ）の周波数分布（振幅（度数）−周波数の関係）を求め、画像評価方法１と同様に、この周波数分布の平均値、分散、標準偏差等の統計量を算出し、算出された統計量に基づいて画質の評価値を決定するようにしてもよい。周波数成分ｉの振幅値をＦ[ｉ]、周波数成分の合計をＮとすると、この分布の平均値ｆM、分散σ2、標準偏差σは、それぞれ、以下の式（１１）〜式（１３）のように表すことができる。
【数５】
2006050494.tif 000007
【００８４】
<画像評価方法３：画像のブレの具合又はボケの度合いに基づく評価>
次に、図１２〜図１６を参照して、撮影画像のブレの具合又はボケの度合いに基づいて画質を評価する方法について説明する。
【００８５】
図１２に、画像の変換や復元等に用いるＰＳＦ（Point Spread Function：点像分布関数、点広がり関数）法等を利用して、ブレやボケを判別評価する方法について示す。一般に、ＰＳＦ法による画像復元では、ブレやボケのない本来の画像i（ｘ，ｙ）と、ブレやボケを生じさせる成分ｐ（ｘ，ｙ）により、劣化した劣化画像ｇ（ｘ，ｙ）を表す。即ち、劣化画像ｇ（ｘ，ｙ）は、式（１４）のように表される。
ｇ（ｘ，ｙ）＝ｐ（ｘ，ｙ）＊ｉ（ｘ，ｙ）      （１４）
ここで、＊は、コンボリューション（畳み込み積分）演算を表す。
【００８６】
ｐ（ｘ，ｙ）を、ＰＳＦ、ＬＳＦ（Line Spread Function：線像分布関数）として、劣化させたｐ（ｘ，ｙ）を算出することができれば、本来の画像ｉ（ｘ，ｙ）を、デコンボリューション（逆畳み込み積分）演算により求めることができる。以下、本来の画像ｉ（ｘ，ｙ）の算出方法を説明する。
【００８７】
まず、式（１４）に示す劣化画像ｇ（ｘ、ｙ）を、フーリエ変換等により周波数軸に変換する。周波数変換された劣化画像をＧ（ｕ、ｖ）とすると、Ｇ（ｕ、ｖ）は、式（１５−１）のように表される。
【数６】
2006050494.tif 000008
【００８８】
また、ｐ（ｘ、ｙ）、ｉ（ｘ、ｙ）を周波数変換したものを、それぞれＰ（ｕ、ｖ）、Ｉ（ｕ、ｖ）とすると、Ｇ（ｕ、ｖ）は、式（１５−２）のように表すこともできる。
Ｇ（ｕ、ｖ）＝Ｐ（ｕ、ｖ）・Ｉ（ｕ、ｖ）      （１５−２）
ここで、コンボリューションは、周波数軸上では、フーリエ変換同士の掛け算となる。周波数領域で、式（１５−２）のＰ（ｕ、ｖ）が算出できれば、逆フィルタとして１／Ｐ（ｕ、ｖ）を計算し、本来の画像ｉ（ｘ、ｙ）を復元することができる。即ち、Ｉ（ｕ、ｖ）＝Ｇ（ｕ、ｖ）／Ｐ（ｕ、ｖ）（フーリエ変換同士の割り算）を算出し、これを逆フーリエ変換してｉ（ｘ、ｙ）を求めればよい。式（１６）にＩ（ｕ、ｖ）の逆フーリエ変換式を示す。
【数７】
2006050494.tif 000009
【００８９】
本実施形態では、このＰＳＦ法におけるＰＳＦやＰ（ｕ、ｖ）の算出方法を、一般のアマチュアの撮影画像で、露出設定の失敗と並んで多い、「焦点ボケ」画像や「直線ブレ」画像の判別や評価に適用する。一般に、直線ブレした画像をフーリエ変換したＧ（ｕ、ｖ）の振幅値の分布は、図１２（ｂ）に示すような直線ブレのパターンとなり、中央付近の高周波成分がブレ角度に従って数本の細長い筋や傾斜した楕円形状に抜けて歪んだパターンが現れる。また、焦点ボケした画像をフーリエ変換したＧ（ｕ、ｖ）は、図１２（ｄ）に示すような焦点ボケの画像パターンとなり、中央部の高周波成分が同心円状に抜けて歪んだパターンとなる。また、直線ブレと焦点ボケが共に生じると、図１２（ｃ）に示すような、両方の特徴を合せもったパターンとなる。
【００９０】
本実施形態では、直線ブレの評価値が、Ｇ（ｕ、ｖ）の振幅値が周期1/Ｌでゼロ交差する方向のブレ方向角度（θ）と、ブレの距離又はその画素数（Ｌ）に基づいて決定されるものとする。また、焦点ボケの評価値が、Ｇ（ｕ，ｖ）の振幅値が周期１．０１π/ｒでゼロ交差する同心円状の焦点の広がり半径又はその画素数（ｒ）に基づいて決定されるものとする。ブレ方向角度（θ）、ブレの距離（Ｌ）、焦点の広がり半径（ｒ）を算出するには、ＰＳＦ法、多階調Ｈｏｕｇｈ変換法、逆フィルタ法等を適用できるが、以下では、ＰＳＦ法と多階調Ｈｏｕｇｈ変換を用いた算出方法について説明する。
【００９１】
まず、図１３（ａ）に示すように、直線ブレ等の劣化画像をフーリエ変換したＧ（ｕ，ｖ）のフーリエ・スペクトルの中央部の円内を、（複数点を通る直線の推定や検出等に広く用いられる）Ｈｏｕｇｈ変換により、θ―ρ空間上の多階調の２次元画像Ｈ（θ、ρ）に変換し、Ｈ（θ、ρ）値の角度θ毎に切り出したＨθ（ρ）波形を比較し、各角度θ（０°≦θ＜180°）におけるエントロピーＥ（θ）を求める。エントロピーＥ（θ）は、式（１７）のように表される。
【数８】
2006050494.tif 000010
式（１７）のＨsum(θ)は、劣化画像Ｇ（ｕ，ｖ）のフーリエ・スペクトルが示す円内の画素値の合計である。式（１７）に示すエントロピーＥ（θ）が最小となる角度θを、直線ブレ方向（θ）とみなすことができる。また、図１３（ｂ）に示すように、直線ブレ方向（θ）のＨθ（ρ）波形の極小点の周期（ＴＬ）を求めると、ブレ距離（Ｌ）＝1／ＴＬが求まる。
【００９２】
焦点ボケの広がり半径（ｒ）を算出するには、まず、図１４に示すように、Ｈ（θ、ρ）をθ軸方向に０°〜１７９°まで平均化したＣ（ρ）波形を作成する。Ｃ（ρ）は、式（１８）のように表される。
【数９】
2006050494.tif 000011
Ｃ（ρ）波形の極小点の周期（Ｔｒ）から、焦点ボケの広がり半径（ｒ）は、半径（ｒ）＝１．０１×π／Ｔｒより求まる。
【００９３】
図１５及び図１６に、上記の方法で用いたＨｏｕｇｈ変換（ハフ変換）の原理を示す。Ｈｏｕｇｈ変換では、ρ＝ｘcosθ＋ｙsinθの変換式により、一般の２次元座標（ｘ−ｙ平面上の点）はθ‐ρ平面上の曲線に変換され、ｘ−ｙ平面上の直線は、θ‐ρ平面上の点に逆変換できる。例えば、ｘ−ｙ平面上の点Ａ、Ｂ、Ｃを通る直線Ｌを求める場合には、図１５に示すように、点Ａ、Ｂ、ＣをＨｏｕｇｈ変換した、θ‐ρ平面上の曲線Ａ’、Ｂ’、Ｃ’が交差する交点Ｐ（θL、ρL）を求める。そして、その交点Ｐ（θL、ρL）をｘ−ｙ平面上にＨｏｕｇｈ逆変換し、ゼロ点から下ろした垂線の長さがρLで、垂線の角度がθLとなる直線Ｌが、求める直線となる。
【００９４】
Ｈｏｕｇｈ変換の配列Ｈ（θ、ρ）の明度（多階調を有する輝度又は濃度）をＨ（θ、ρ）＝ｆ（ｘ、ｙ）として多階調の２次元画像として定義し、図１５に示すような複数曲線の交点を求める場合には、図１６に示すように、明度（階調値）を加算して、階調値が最も大きくなる点（図１６の斜線部分）を交点として、近似計算で迅速に求めることができる。
【００９５】
以上のように、画像を劣化させたＰＳＦ関数を算出する方法、ゼロクロス法、多階調Ｈｏｕｇｈ変換等により、劣化画像ｇ（ｘ、ｙ）を周波数軸上にフーリエ変換したＧ（ｕ，ｖ）から、直線ブレ画像の方向（角度θ）、直線ブレの距離（画素数Ｌ）、焦点ボケ画像の焦点広がり半径（画素数ｒ）を算出し、これらθ、Ｌ、ｒに基づいて画像の評価値が決定される。θ、Ｌ、ｒの値が大きくなるほど評価値は小さくなる。
【００９６】
図１７に、本実施形態で適用される画質補正処理の概念図を示す。
画質補正処理は、デジタルカメラ１に入力された撮影画像の各画素の輝度値（又は濃度値）f(i,j)或いはその周辺の画素配列に、所定の変換表などを参照しながら、画質補正演算を施し、補正後の出力画像ｇ（i,j）或いは所定の特徴量を出力するものである。
【００９７】
変換表としては、例えば、輝度変換などの入出力特性の変換式、変換用のＬＵＴ（Look Up Table）（変換表）、加重マトリックス（行列）、空間フィルタ係数などの演算子（オペレータ）等を利用すればよく、画質補正演算には、点処理演算、注目画素と周囲８点等を対象にする局所処理演算、画像全体に演算を施す大局処理演算を適用することができる。
【００９８】
点処理演算としては、濃度ヒストグラム変換処理、２値化処理、座標変換処理、画像間演算がある。局所処理演算としては、ノイズ除去演算、微分演算、積和演算、フィルタ演算、エッジ抽出処理、膨張収縮処理、細線化処理などがある。大局処理演算としては、フーリエ変換等の直交変換演算、テクスチャー解析処理等を適用することができる。
【００９９】
次に、図２及び図３の画質の自動補正処理で適用される画像補正方法１〜７について説明する。
【０１００】
<画像補正方法１：ガンマ補正>
まず、図１８を参照して、ガンマ補正処理について説明する。
ガンマ（γ）補正は、ＣＣＤ（Charge Coupled Device）等の撮像素子の入出力特性曲線を、変換テーブルを用いて補正するもので、入力画像の輝度又は濃度をf(i,j)＝Ｌ、γ補正後の出力画像の輝度又は濃度をg(i,j)＝Ｌ'とすると、γ補正の変換式は、式（１９）のように定義される。
【数１０】
2006050494.tif 000012
図１８の変換テーブルαは、式（１９）においてγ＝０．５、（定数）＝（２５５）-2とした場合の入力Ｌと出力Ｌ'の関係を示すテーブルである。
【０１０１】
<画像補正方法２：輝度又は濃度のヒストグラム分布の変換による補正>
次に、図１９及び図２０を参照して、撮影画像の輝度又は濃度のヒストグラム分布を変換することによる補正について説明する。
【０１０２】
入力された撮影画像の輝度をf(i,j)＝ｘ、変換後の画像の輝度をg(i,j)＝ｙとすると、図１９に示すように、撮影画像の輝度のヒストグラム分布Ｐ（ｘ）を、線形又は非線形の輝度変換式ｙ＝Ｑ（ｘ）によって変換すると、Ｐ（ｘ）とは異なる特性の輝度ヒストグラム分布Ｐ（ｙ）が得られる。輝度ヒストグラム分布Ｐ（ｘ）における輝度ｘの最小値をａ、最大値をｂ、変換後の輝度ヒストグラム分布Ｐ（ｙ）における輝度ｙの最小値をｕ、最大値をｖとすると、図１９に示した輝度変換式ｙ＝Ｑ（ｘ）は、式（２０）のように表される。
【数１１】
2006050494.tif 000013
ｖ−ｕ＞ｂ−ａである場合、ダイナミックレンジが狭くコントラストの悪い入力画像の輝度ヒストグラム分布Ｐ（ｘ）を輝度変換式（２０）により変換することにより、ダイナミックレンジが広がり、コントラストや中間諧調の濃淡表現を改善することができる。なお、式（２０）における定数ｖ、ｕ、ａ、ｂの値を変えれば、ダイナミックレンジの収縮や、輝度ヒストグラム分布の左右シフト（全体の輝度を上げる、下げる）等にも応用できる。
【０１０３】
図２０に、その他の輝度変換式の例を示す。
図２０（ａ）は、画像の輝度を、（輝度が低いものほど）全体的に暗くする補正処理を行う場合の非線形の輝度変換式の一例であり、式（２１）のように表される。
ｙ＝ｖ（ｘ/ｂ）2      （２１）
【０１０４】
図２０（ｂ）及び（ｃ）は、画像の輝度を、（輝度が高いものほど）全体的に明るくする補正処理を行う場合の非線形の輝度変換式の例であり、それぞれ、式（２２）、式（２３）のように表される。
ｙ＝−ｖ｛（ｘ−ｂ）/ｂ｝2＋ｖ      （２２）
ｙ/ｖ＝log(1+μ・ｘ/ｂ)/log(1+μ)      （２３）
【０１０５】
また、図２０（ｄ）及び（ｅ）は、輝度が高い領域の輝度を上げ、輝度が低い領域の輝度を下げる場合の輝度変換式の例である。図２０（ｄ）に示すＳ字状の輝度変換式は、式（２４）のように表される。
【数１２】
2006050494.tif 000014
図２０（ｅ）に示す輝度変換式は、式（２５）のように表される。
【数１３】
2006050494.tif 000015
図２０（ｄ）又は（ｅ）に示すような輝度変換式を用いると、階調表現のトーンカーブ（色調曲線）やコントラストを改善することができる。
【０１０６】
このように、線形や非線形の輝度又は濃度の変換式により、濃度や輝度のヒストグラム分布を変換することにより、露出条件の設定が悪くて「白飛び」や「黒つぶれ」した画像、コントラストや階調表現の悪い画像等のダイナミックレンジの拡大や、コントラストの改善、画面の明暗、中間諧調の濃淡表現、トーンカーブの改善など、画質の補正や加工を行うことができる。
【０１０７】
また、ＲＧＢ各色別の輝度ヒストグラム分布毎に上記の補正処理を施して、ヒストグラム分布を改善したり、ＲＧＢ各色の輝度の加重バランスなどを補正するようにしてもよい。例えば、水中撮影写真等によくあるように、Ｂ（青）色成分の露出オーバー気味で輝度は高いが、Ｒ（赤）色成分が（水に吸収されて）減じて露出アンダーになってしまった画像では、Ｂ（青）色成分のヒストグラムを少し左側にシフトして弱め、逆にＲ（赤）色成分のヒストグラムを少し右側にシフトして強めるなど、ＲＧＢ各色のヒストグラムの偏りを調整して、全体の露出バランスを補正することができる。
【０１０８】
〈画像補正方法３：輝度又は濃度のヒストグラム分布の平坦化による補正〉
次に、図２１及び図２２を参照して、撮影画像の輝度又は濃度のヒストグラム分布の平坦化による補正方法について説明する。
【０１０９】
図２１に示すように、入力画像の輝度又は濃度のヒストグラム分布Ｐ（ｘ）の偏った分布（図２１（ａ））を、Ｐ（ｘ）の積分値を階調変換式として用いることにより、一様にフラットなヒストグラム分布Ｐ（ｙ）（図２１（ｂ））が得られる。
【０１１０】
図２１に示すヒストグラム分布の平坦化処理について、図２２のフローチャートを参照して説明する。
【０１１１】
補正対象の画像が入力されると（ステップＳ６０）、入力画像の横画素数×縦画素数が全画素数Ｎとして設定される（ステップＳ６１）。次いで、補正後の出力画像の階調数ｍ（０〜ｍ−１）が決定され（ステップＳ６２）、平均出現頻度数（１階調あたりの画素数）がＮ／ｍに設定される（ステップＳ６３）。
【０１１２】
次いで、入力画像のヒストグラムが計算される（ステップＳ６４）。次いで、階調輝度（又は濃度）をカウントするためのカウンタ値ｋ（０〜ｍ−１）が０に設定され（ステップＳ６５）、階調ｋの合計数をカウントするためのカウンタ値Ｈist(k)が０に設定される（ステップＳ６６）。
【０１１３】
次いで、入力画像の階調輝度（又は濃度）が低い方から順に画素が取り出され、出力画像の該当画素の階調がｋに設定される（ステップＳ６７）。そして、階調ｋの合計数Ｈist(k)がインクリメントされ（ステップＳ６８）、現在の合計数Ｈist(k)が平均出現頻度数Ｎ／ｍより大きいか否かが判定される（ステップＳ６９）。
【０１１４】
ステップＳ６９において、Ｈist(k)がＮ／ｍ以下であると判定された場合（ステップＳ６９；ＮＯ）、ステップＳ６７に戻り、ステップＳ６７及びＳ６８の処理が繰り返される。ステップＳ６９において、Ｈist(k)がＮ／ｍより大きいと判定された場合（ステップＳ６９；ＹＥＳ）、階調ｋがインクリメントされ（ステップＳ７０）、現在の階調ｋが、階調の最大値ｍ−１より大きいか否かが判定される（ステップＳ７１）。
【０１１５】
ステップＳ７１において、階調ｋがｍ−１以下であると判定された場合（ステップＳ７１；ＮＯ）、ステップＳ６６に戻り、ステップＳ６６〜Ｓ７０の処理が繰り返される。ステップＳ７１において、階調ｋがｍ−１より大きいと判定された場合（ステップＳ７１；ＹＥＳ）、階調変換後の画像が出力され（ステップＳ７２）、本ヒストグラム平坦化処理が終了する。
【０１１６】
〈画像補正方法４：画像の平滑化による補正〉
次に、図２３を参照して、撮影画像の平滑化による補正方法について説明する。
【０１１７】
入力画像の着目画素とその周囲からなる９画素に対し、図２３（ａ）又は（ｂ）に示すような平均値（ＭＥＡＮ）フィルタＡ（Ｋ,Ｌ）による空間フィルタ演算を施すことにより、画像全体が平滑化される。入力画像をｆ（i,j）、フィルタ演算後の出力画像をｇ(i,j)とすると、平均値フィルタＡ（Ｋ,Ｌ）を用いた空間フィルタ演算は、式（２６）のように表される。
【数１４】
2006050494.tif 000016
このようなフィルタ演算によって、入力画像のごま塩ノイズが除去され、明瞭度が落ち、ぼかす効果がある。明瞭度をあまり落とさずに、ノイズを除去するためには、メディアン（中央値）フィルタを用いた処理が有効である。メディアンフィルタは、着目画素とその周囲の９画素のうち、輝度が大きい方から５番目の中央値を、着目画素の新しい輝度値とするフィルタである。
【０１１８】
〈画像補正方法５：微分フィルタを用いた鮮鋭化による補正〉
次に、図２４及び図２５を参照して、撮影画像の輪郭やエッジ成分を強調して、画像のボケを減じて、鮮鋭化する補正処理について説明する。
【０１１９】
図２４（ａ）のように、階調が変化する部分のエッジがボケた画像f（i,j）を一次微分することによって、画像f（i,j）の勾配（Gradient）を求めると、図２４（ｂ）のように、階調が変化する勾配部分が抽出され、凹凸が得られる。画像f（i,j）の一次微分は、式（２７）のように表される。
Δxf(i,j)＝f(i+1,j)−f(i-1,j)、    Δyf(i,j)＝f(i,j+1)−f(i,j-1)    （２７）
ここで、Δxは、ｘ方向（横方向）の一次微分であり、Δyは、ｙ方向（縦方向）の一次微分を表す。また、方向をもたない勾配は、式（２８）、式（２９）、式（３０）のように表される。
【数１５】
2006050494.tif 000017
【０１２０】
画像f（i,j）の一次微分を更に微分する二次微分では、図２４（ｃ）に示すように、左右逆向きのＮ字状の波形となり、高周波成分が強調された画像が得られる。画像f（i,j）の二次微分（Laplaciant）は、式（３１）又は式（３２）のように表される。
【数１６】
2006050494.tif 000018
【０１２１】
式（２７）に示す一次微分は、図２５（ａ）に示すようなPrewittフィルタ、図２５（ｂ）に示すようなSobelフィルタ、Kirschフィルタ、Robertsフィルタ等が用いられる。図２５において、入力画像にPrewittフィルタを作用させると、Ａ−Ａ'線上の入力画像は、図２５（ｄ）に示すような出力画像となり、入力画像にSobelフィルタを作用させると、Ａ−Ａ'線上の入力画像は、図２５（ｅ）に示すような出力画像となる。また、原画像（入力画像）から、二次微分が施された画像を差し引く演算を行うフィルタの一例を図２５（ｃ）に示す。図２５（ｃ）に示すようなフィルタを用いると、図２５（ｆ）に示すように、エッジ部分の高周波成分が強調された画像が得られる。即ち、このフィルタは、画像のボケたエッジや輪郭などを強調し、画像の鮮明度（シャープネス）を増加し、鮮鋭化する補正に利用することができる。
【０１２２】
図２５（ａ）〜（ｃ）に示すような微分フィルタをＷ（Ｋ,Ｌ）とすると、微分フィルタＷ（Ｋ,Ｌ）を用いたフィルタ演算は、式（３３）のように表される。
【数１７】
2006050494.tif 000019
例えば、図２５（ｃ）に示したフィルタＷ（Ｋ,Ｌ）による演算は、式（３３−１）のように表される。
【数１８】
2006050494.tif 000020
【０１２３】
〈画像補正方法６：周波数変換された画像に鮮鋭化処理を施す補正〉
次に、周波数変換された撮影画像に鮮鋭化処理を施すことによる補正処理について説明する。
【０１２４】
画像補正方法６は、２次元の離散フーリエ変換（ＤＦＴ）、離散コサイン変換（ＤＣＴ）等の周波数変換により周波数空間上に変換した画像データの低周波成分或いは高周波成分に対して、鮮鋭化処理、ぼかし処理、輪郭強調処理等のフィルタ処理を施すものである。２次元画像ｆ（ｘ,ｙ）の離散フーリエ変換の変換式、離散コサイン変換の変換式は、それぞれ、画像評価方法２で示した式（７）、式（８）と同一である。
【０１２５】
周波数変換されたＦ（ｕ，ｖ）の直流成分及び低周波成分を残し（又は、増強し）、高周波成分を０にすると、低周波通過フィルタ（ＬＰＦ）と同等の処理となり、エッジや輪郭を弱め、画像をぼかす補正となる。また、周波数変換されたＦ（ｕ，ｖ）の直流成分及び高周波成分を残し（又は増強し）、低周波成分を０とすると、高周波通過フィルタ（ＨＰＦ）と同等の処理となり、エッジや輪郭を強調する補正となる。また、周波数変換されたＦ（ｕ，ｖ）の高周波成分の値を増大させるか、Ｆ（ｕ，ｖ）に一次微分フィルタ又は二次微分フィルタ（ラプラシアンフィルタ）を作用させた場合も、エッジや輪郭を強調する補正となる。周波数空間上で上述のような処理が施された画像を周波数逆変換（離散フーリエ逆変換、離散コサイン逆変換等）すると、鮮鋭化処理、ぼかし処理等が施された画像を得ることができる。
【０１２６】
〈画像補正方法７：ＰＳＦ法を用いた補正〉
次に、ＰＳＦ法等を用いて、手ブレなどによる直線ブレ画像、ピンボケ等の焦点ボケ画像の補正処理について説明する。
【０１２７】
直線ブレや焦点ボケのない本来の画像i(x,y)が、ブレやボケなどの劣化成分ｐ(x,y)により劣化した画像ｇ(x,y)をフーリエ変換したＧ（ｕ，ｖ）を求め、その周波数空間上のパターンから、画像評価方法３に示したように、直線ブレの方向（角度θ）、直線ブレの距離（画素数Ｌ）、焦点ボケの広がり半径（画素数ｒ）が算出される。
【０１２８】
θ方向に周期１/Ｌでゼロ交差する直線ブレが生じている場合、式（１５−２）のＰ（ｕ、ｖ）は、下記の式（３４）のように表される。
Ｐ（ｕ、ｖ）＝sin(πｆＬ)／πｆＬ        （３４）
ここで、ｆ＝ｕcosθ＋ｖsinθである。そして、Ｇ（ｕ、ｖ）と、式（３４）のＰ（ｕ、ｖ）から、ブレの無い本来の画像のフーリエ変換画像Ｉ（ｕ、ｖ）＝Ｇ（ｕ、ｖ）／Ｐ（ｕ、ｖ）を求め、このＩ（ｕ、ｖ）を式（１６）に示すように逆フーリエ変換して本来の画像ｉ（ｘ、ｙ）を求めることができる。
【０１２９】
周期１．０１π/ｒで同心円状のゼロ交差をする焦点ボケが生じている場合、式（１５−２）のＰ（ｕ、ｖ）は、下記の式（３５）のように表される。
Ｐ（ｕ、ｖ）＝２・Ｊ１（ｒ・Ｒ）／ｒ・Ｒ      （３５）
ここで、式（３５）のＪ１は、１次の第１種ベッセル関数である。そして、上述と同様に、受信された画像データをフーリエ変換したＧ（ｕ、ｖ）と、式（３５）のＰ（ｕ、ｖ）から、焦点ボケの無い本来の画像のフーリエ変換画像Ｉ（ｕ、ｖ）＝Ｇ（ｕ、ｖ）／Ｐ（ｕ、ｖ）を求め、このＩ（ｕ、ｖ）を式（１６）に示すように逆フーリエ変換して本来の画像ｉ（ｘ、ｙ）を求めることができる。
【０１３０】
なお、一般的なＰＳＦ、ＬＳＦ、逆フィルタ法、最小２乗法等の方法を用いてｐ（ｘ，ｙ）やＰ（ｕ，ｖ）を算出してもよいし、ウィナー・フィルタ等のその他の方法で逆フィルタ（インバースフィルタ）を求めたり、デコンボリューション演算を行って、画像の復元やブレ、ボケの補正を行うようにしてもよい。
【０１３１】
以上のように、本実施形態のデジタルカメラ１によれば、記録媒体の残り容量が予め設定された容量以下である場合、撮影画像を自動的に評価して、評価値が所定値以下の画像を自動的に補正し、補正処理後であっても画像の評価値が所定値以下である画像のみを削除するようにしたことにより、必要以上に画像を削除することなくなるとともに、記録媒体の空き領域を確実に確保することができる。
【０１３２】
（変形例）
次に、図２６〜図２８を参照して、上述の実施形態の変形例について説明する。本変形例では、デジタルカメラ１の生産工程で、ブレ、ボケ等の劣化成分を計測してデジタルカメラ１内に設定しておき、この設定された劣化成分に基づいて撮影画像の補正処理を行う例について説明する。なお、本変形例では、デジタルカメラ１の生産工程で設定される画像の劣化成分が、直線ブレ、焦点ボケである場合を示すが、他の劣化成分についても本変形例を適用することができる。
【０１３３】
まず、図２６のフローチャートを参照して、デジタルカメラ１の生産工程において実行されるブレ劣化成分設定処理について説明する。
【０１３４】
まず、撮影対象である見本画像ｉ(x,y)が設定されるとともに（ステップＳ８０）、直線ブレ方向角度θが０に設定され（ステップＳ８１）、更に、直線ブレ距離Ｌが０に設定される（ステップＳ８２）。次いで、ブレ方向角度がθ、ブレ距離がＬの直線ブレを生じさせて見本画像が撮影され、撮影画像ｇ(x,y)が取得される（ステップＳ８３）。ステップＳ８３では、ロボットハンドや、専用の実験装置等を用いて、ブレが生じるように撮影される。
【０１３５】
次いで、撮影画像ｇ(x,y)がフーリエ変換され、周波数空間での撮影画像Ｇ(u,v)が求められる（ステップＳ８４）。次いで、見本画像ｉ(x,y)をフーリエ変換したＩ(u,v)と、撮影画像をフーリエ変換したＧ(u,v)により、劣化成分p(x,y)をフーリエ変換したＰ(u,v)＝Ｇ(u,v)/Ｉ(u,v)が算出され、この算出されたＰ(u,v)が、ブレ方向角度がθ、ブレ距離がＬの劣化成分Ｐ１(θ,Ｌ)として、制御部３２のメモリ（ＲＯＭ）に記憶される（ステップＳ８５）。
【０１３６】
次いで、ブレ距離Ｌがインクリメントされ（ステップＳ８６）、現在のブレ距離Ｌが予め設定された最大値Ｌmax以上であるか否かが判定される（ステップＳ８７）。ステップＳ８７において、ブレ距離ＬがＬmax未満であると判定された場合（ステップＳ８７；ＮＯ）、ステップＳ８３に戻り、ステップＳ８３〜Ｓ８６の処理が繰り返される。
【０１３７】
ステップＳ８７において、ブレ距離ＬがＬmax以上であると判定された場合（ステップＳ８７；ＹＥＳ）、ブレ方向角度θがインクリメントされ（ステップＳ８８）、現在のブレ方向角度θが予め設定された最大値θmax以上であるか否かが判定される（ステップＳ８９）。
【０１３８】
ステップＳ８９において、ブレ方向角度θがθmax未満であると判定された場合（ステップＳ８９；ＮＯ）、ステップＳ８３に戻り、ステップＳ８３〜Ｓ８８の処理が繰り返される。ステップＳ８９において、ブレ方向角度θがθmax以上であると判定された場合（ステップＳ８９；ＹＥＳ）、本ブレ劣化成分設定処理が終了する。
【０１３９】
次に、図２７のフローチャートを参照して、デジタルカメラ１の生産工程において実行される焦点ボケ劣化成分設定処理について説明する。
【０１４０】
まず、撮影対象である見本画像ｉ(x,y)が設定されるとともに（ステップＳ１００）、焦点ボケ広がり半径ｒが０に設定される（ステップＳ１０１）。次いで、広がり半径がｒの焦点ボケを生じさせて見本画像が撮影され、撮影画像ｇ(x,y)が取得される（ステップＳ１０２）。ステップＳ１０２では、ロボットハンドや、専用の実験装置等を用いて、焦点ボケが生じるように撮影される。
【０１４１】
次いで、撮影画像ｇ(x,y)がフーリエ変換され、周波数空間での撮影画像Ｇ(u,v)が求められる（ステップＳ１０３）。次いで、見本画像ｉ(x,y)をフーリエ変換したＩ(u,v)と、撮影画像をフーリエ変換したＧ(u,v)により、劣化成分p(x,y)をフーリエ変換したＰ(u,v)＝Ｇ(u,v)/Ｉ(u,v)が算出され、この算出されたＰ(u,v)が、焦点ボケ広がり半径ｒの劣化成分Ｐ２(ｒ)として、制御部３２のメモリ（ＲＯＭ）に記憶される（ステップＳ１０４）。
【０１４２】
次いで、広がり半径ｒがインクリメントされ（ステップＳ１０５）、現在の広がり半径ｒが予め設定された最大値ｒmax以上であるか否かが判定される（ステップＳ１０６）。ステップＳ１０６において、広がり半径ｒがｒmax未満であると判定された場合（ステップＳ１０６；ＮＯ）、ステップＳ１０２に戻り、ステップＳ１０２〜Ｓ１０５の処理が繰り返される。ステップＳ１０６において、広がり半径ｒがｒmax以上であると判定された場合（ステップＳ１０６；ＹＥＳ）、本ブレ劣化成分設定処理が終了する。
【０１４３】
次に、図２８のフローチャートを参照して、出荷後のデジタルカメラ１において実行される画像撮影・補正処理について説明する。
【０１４４】
キー入力部３６の操作により撮影モードが指定されると、測光処理、露出条件の設定等が行われる（ステップＳ１１０）。測光処理後、キー入力部３６のシャッターキーの押下により撮影が指示されると、撮影処理が行われ、撮影画像ｇ(x,y)が取得される（ステップＳ１１１）。取得された撮影画像ｇ(x,y)は、記録媒体（フラッシュメモリ３８）に記録される。
【０１４５】
ステップＳ１１１の撮影処理において、手ブレ検出センサ４０によりブレが検出された場合、検出されたブレの方向θ、ブレ距離Ｌが制御部３２のＲＡＭに記憶される（ステップＳ１１２）。また、撮影画像ｇ(x,y)から焦点ボケが検出された場合、焦点ボケの広がり半径ｒが制御部３２のメモリ（ＲＡＭ）に記憶される（ステップＳ１１２）。
【０１４６】
次いで、撮影画像ｇ(x,y)がフーリエ変換され、周波数空間での撮影画像Ｇ(u,v)が求められる（ステップＳ１１３）。次いで、ステップＳ１１２において記憶されたブレ距離Ｌが予め設定された値（設定値）以上であるか否か（即ち、評価値が所定値以下であるか否か）が判定される（ステップＳ１１４）。
【０１４７】
ステップＳ１１４において、ブレ距離Ｌが設定値未満である（即ち、評価値が所定値より大きい）と判定された場合（ステップＳ１１４；ＮＯ）、後述のステップＳ１１９に移行する。ステップＳ１１４において、ブレ距離Ｌが設定値以上である（即ち、評価値が所定値以下である）と判定された場合（ステップＳ１１４；ＹＥＳ）、ステップＳ１１２において記憶されたブレ方向θ及びブレ距離Ｌに対応する劣化成分Ｐ１(θ,Ｌ)が、制御部３２のメモリ（ＲＯＭ）から読み出され（ステップＳ１１５）、Ｇ(u,v)及びＰ１(θ,Ｌ)から、ブレのない本来の画像ｉ(x,y)をフーリエ変換したＩ(u,v)＝Ｇ(u,v)/Ｐ１(θ,Ｌ)が算出される（ステップＳ１１６）。
【０１４８】
次いで、ステップＳ１１６で算出されたＩ(u,v)がフーリエ逆変換されることにより、復元画像ｉ(x,y)が算出され（ステップＳ１１７）、記録媒体に記録されたｇ(x,y)が復元画像ｉ(x,y)に更新される（ステップＳ１１８）。
【０１４９】
次いで、ステップＳ１１２において記憶された焦点ボケの広がり半径ｒが予め設定された値（設定値）以上であるか否か（即ち、評価値が所定値以下であるか否か）が判定される（ステップＳ１１９）。ステップＳ１１９において、焦点ボケの広がり半径ｒが設定値未満である（即ち、評価値が所定値より大きい）と判定された場合（ステップＳ１１９；ＮＯ）、本画像撮影・補正処理が終了する。
【０１５０】
ステップＳ１１９において、焦点ボケの広がり半径ｒが設定値以上である（即ち、評価値が所定値以下である）と判定された場合（ステップＳ１１９；ＹＥＳ）、ステップＳ１１２において記憶された焦点ボケの広がり半径ｒに対応する劣化成分Ｐ２(ｒ)が、制御部３２のメモリ（ＲＯＭ）から読み出され（ステップＳ１２０）、Ｇ(u,v)及びＰ２(ｒ)から、焦点ボケのない本来の画像ｉ(x,y)をフーリエ変換したＩ(u,v)＝Ｇ(u,v)/Ｐ２(ｒ)が算出される（ステップＳ１２１）。
【０１５１】
次いで、ステップＳ１２１で算出されたＩ(u,v)がフーリエ逆変換されることにより、復元画像ｉ(x,y)が算出され（ステップＳ１２２）、記録媒体に記録されたｇ(x,y)が復元画像ｉ(x,y)に更新され（ステップＳ１２３）、本画像撮影・補正処理が終了する。
【０１５２】
以上のように、デジタルカメラ１の生産工程において、予め、ブレやボケを生じさせるように見本画像を撮影して、ブレやボケの度合い別に劣化成分Ｐ(u,v)を算出してメモリに設定し、出荷後のデジタルカメラ１で撮影された画像を、メモリに設定された劣化成分Ｐ(u,v)に基づいて補正することにより、個々のカメラのレンズの性能に合わせた画像補正を容易に行うことができる。また、出荷後において劣化成分Ｐ(u,v)を算出する必要がないため、画像補正処理を高速に行うことができる。
【０１５３】
最近は、オートフォーカス機能が搭載されたデジタルカメラが普及しているが、カメラの機種毎に、オートフォーカスの精度や特性が異なっている。また、露出の条件や周囲の明るさ、被写体の配置等によって、常に所望の被写体に焦点が合った画像を撮影できるとは限らず、また、絞りの度合いや焦点距離によって、被写界深度も変化したりする。従って、図２６〜図２８に示すように、生産出荷前に、画像劣化の度合い別に劣化成分を算出してメモリに設定しておき、出荷後に撮影された画像を、メモリに設定された劣化成分に基づいて補正することにより、個々の機種でのレンズ特性やフォーカス特性の違いにも対応した補正を行うことができる。
【０１５４】
なお、本実施形態及び変形例における記述内容は、本発明の趣旨を逸脱しない範囲で適宜変更可能である。
【図面の簡単な説明】
【０１５５】
【図１】本発明の画像撮影装置を適用したデジタルカメラの回路構成を示すブロック図。
【図２】本実施形態のデジタルカメラにおいて実行される画像撮影・削除処理を示すフローチャート。
【図３】図２における画像削除処理の詳細を示すフローチャート。
【図４】デジタルカメラの表示部における画面表示例を示す図。
【図５】デジタルカメラの表示部における画面表示例を示す図。
【図６】デジタルカメラの表示部における画面表示例を示す図。
【図７】画像の輝度ヒストグラム分布から画像を評価する方法の一例を示す図。
【図８】画像の輝度ヒストグラム分布から画像を評価する方法の一例を示す図。
【図９】画像の周波数特性により画像を評価する方法の一例を示す図。
【図１０】画像の周波数特性により画像を評価する方法の一例を示す図。
【図１１】画像の周波数特性により画像を評価する方法の一例を示す図。
【図１２】画像の直線ブレ、焦点ボケの度合いの測定による画像評価方法の一例を示す図。
【図１３】直線ブレの方向の角度（θ）の算出方法と、直線ブレの距離（Ｌ）の算出方法を示す図。
【図１４】焦点ボケの広がり半径（ｒ）の算出方法を示す図。
【図１５】ハフ（Ｈｏｕｇｈ）変換を説明するための図。
【図１６】多階調ハフ変換を説明するための図。
【図１７】画質補正処理の概念図。
【図１８】ガンマ補正による画像補正方法を示す図。
【図１９】輝度又は濃度のヒストグラム分布の変換による画像補正方法を示す図。
【図２０】輝度又は濃度のヒストグラム分布の変換による画像補正方法を示す図。
【図２１】輝度又は濃度のヒストグラムの平坦化による画像補正方法を示す図。
【図２２】ヒストグラム平坦化処理を示すフローチャート。
【図２３】画像の平滑化による画像補正方法を示す図。
【図２４】微分フィルタを用いた鮮鋭化による画像補正方法を示す図。
【図２５】微分フィルタを用いた鮮鋭化による画像補正方法を示す図。
【図２６】デジタルカメラの生産工程に実行されるブレ劣化成分設定処理を示すフローチャート。
【図２７】デジタルカメラの生産工程に実行される焦点ボケ劣化成分設定処理を示すフローチャート。
【図２８】画像撮影・補正処理を示すフローチャート。
【符号の説明】
【０１５６】
１  デジタルカメラ（画像撮影装置）
２  外部機器
３  画像処理装置
４  表示装置
５  印刷装置
１３  表示部
２１  モータ（Ｍ）
２２  レンズ光学系
２３  ＣＣＤ
２４  ＴＧ（タイミング発生器）
２５  垂直ドライバ
２６  Ｓ／Ｈ（サンプルホールド回路）
２７  Ａ／Ｄ変換器
２８  カラープロセス回路
２９  ＤＭＡコントローラ
３０  ＤＲＡＭＩ／Ｆ
３１  ＤＲＡＭ
３２  制御部
３３  ＶＲＡＭコントローラ
３４  ＶＲＡＭ
３５  デジタルビデオエンコーダ
３６  キー入力部
３７  ＪＰＥＧ回路
３８  フラッシュメモリ
３９  通信部
４０  手ブレ検出センサ
【図１】
2006050494.tif 000021
【図２】
2006050494.tif 000022
【図３】
2006050494.tif 000023
【図４】
2006050494.tif 000024
【図５】
2006050494.tif 000025
【図６】
2006050494.tif 000026
【図７】
2006050494.tif 000027
【図８】
2006050494.tif 000028
【図９】
2006050494.tif 000029
【図１０】
2006050494.tif 000030
【図１１】
2006050494.tif 000031
【図１２】
2006050494.tif 000032
【図１３】
2006050494.tif 000033
【図１４】
2006050494.tif 000034
【図１５】
2006050494.tif 000035
【図１６】
2006050494.tif 000036
【図１７】
2006050494.tif 000037
【図１８】
2006050494.tif 000038
【図１９】
2006050494.tif 000039
【図２０】
2006050494.tif 000040
【図２１】
2006050494.tif 000041
【図２２】
2006050494.tif 000042
【図２３】
2006050494.tif 000043
【図２４】
2006050494.tif 000044
【図２５】
2006050494.tif 000045
【図２６】
2006050494.tif 000046
【図２７】
2006050494.tif 000047
【図２８】
2006050494.tif 000048