(57)【要約】【課題】撮影時に自動で適切な調整値を決定し、回復フィルタを生成可能な画像処理装置を提供すること。【解決手段】撮影条件に応じた光学伝達関数を取得する取得手段と、光学伝達関数の特性に基づいて、光学伝達関数のゲインを決定する決定手段と、光学伝達関数およびゲインを用いて、画像回復処理を行うための回復フィルタを生成する生成手段と、を有する。【選択図】図２
2015226086.tif 000002
【特許請求の範囲】
【請求項１】
撮影条件に応じた光学伝達関数を取得する取得手段と、
前記光学伝達関数の特性に基づいて、前記光学伝達関数のゲインを決定する決定手段と、
前記光学伝達関数および前記ゲインを用いて、画像回復処理を行うための回復フィルタを生成する生成手段と、を有することを特徴とする画像処理装置。
【請求項２】
前記決定手段は、前記ゲインが所定値となる周波数および前記ゲインの微分値が所定の条件を満たす周波数のうち少なくともいずれかと撮像装置のナイキスト周波数とを比較した結果を用いて前記ゲインを決定することを特徴とする請求項１に記載の画像処理装置。
【請求項３】
前記撮影条件は、撮像光学系と撮像装置の識別子、焦点距離、絞りおよび撮影距離の少なくともいずれかであることを特徴とする請求項１または２に記載の画像処理装置。
【請求項４】
前記回復フィルタを記憶する記憶手段を更に有し、
前記生成手段は、前記回復フィルタを生成した後、同一の撮影条件で前記回復フィルタを生成する際には、前記記憶手段から前記回復フィルタを取得することを特徴とする請求項１から３のいずれか１項に記載の画像処理装置。
【請求項５】
前記記憶手段は、所定の容量に収まるように、優先度に応じて前記回復フィルタを削除することを特徴とする請求項４に記載の画像処理装置。
【請求項６】
前記記憶手段は、前記回復フィルタの使用回数および前記回復フィルタのデータサイズの少なくともいずれかに基づいて、前記優先度を決定していることを特徴とする請求項５に記載の画像処理装置。
【請求項７】
撮像素子と、
前記撮像素子からの信号に基づいて画像処理を行う請求項１から６のいずれか１項に記載の画像処理装置と、を有することを特徴とする撮像装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像回復処理に用いる回復フィルタを生成可能な画像処理装置に関する。
【背景技術】
【０００２】
デジタルカメラ等の撮像装置により得られた画像は、撮像系の球面収差、コマ収差、像面湾曲および非点収差等によって起こるボケによって、画質が劣化する。これらの収差は、点像分布関数により表される。点像分布関数をフーリエ変換することにより得られる光学伝達関数は、収差の周波数空間における情報である。以下、光学伝達関数の情報を用いて画像の劣化を回復、復元する画像回復処理の方法について説明する。
【０００３】
劣化画像をｇ（ｘ，ｙ）、元画像をｆ（ｘ，ｙ）、劣化画像ｇ（ｘ，ｙ）を取得するために用いた撮像系の点像分布関数をｈ（ｘ，ｙ）としたとき、以下の式（１）が成り立つ。ただし、＊はコンボリューション（畳み込み積分、積和）を示し、（ｘ，ｙ）は実空間における画像の座標を示す。
【０００４】
ｇ（ｘ，ｙ）＝ｈ（ｘ，ｙ）＊ｆ（ｘ，ｙ）      （１）
式（１）をフーリエ変換して周波数空間での表示形式に変換すると、式（２）のように表される。
【０００５】
Ｇ（ｕ，ｖ）＝Ｈ（ｕ，ｖ）・Ｆ（ｕ，ｖ）      （２）
ここで、Ｈ（ｕ，ｖ）は点像分布関数ｈ（ｘ，ｙ）をフーリエ変換した光学伝達関数である。Ｇ（ｕ，ｖ），Ｆ（ｕ，ｖ）は、それぞれ劣化画像ｇ（ｘ，ｙ），元画像ｆ（ｘ，ｙ）をフーリエ変換したものである。（ｕ，ｖ）は２次元周波数空間での周波数（座標）を示す。
【０００６】
式（２）の両辺をＨ（ｕ，ｖ）で除算したものが式（３）である。式（３）の両辺を逆フーリエ変換したものが式（４）である。
【０００７】
Ｇ（ｕ，ｖ）／Ｈ（ｕ，ｖ）＝Ｆ（ｕ，ｖ）      （３）
ｇ（ｘ，ｙ）＊Ｒ（ｘ，ｙ）＝ｆ（ｘ，ｙ）      （４）
ここで、Ｒ（ｘ，ｙ）は、１／Ｈ（ｕ，ｖ）を逆フーリエ変換したものであり、回復フィルタである。回復フィルタとは、従来のシャープネスやメディアンフィルタ等の一般的なフィルタ処理と同様に、対象画素と、その周囲の画素に適用するフィルタである。回復フィルタは、光学伝達関数を元に生成するので、画像の領域ごとに異なる。そのため、画像の領域ごとに適切な補正が適切なゲインで行われる。
【０００８】
周辺光量補正や倍率色収差、歪曲の補正については、設計値の逆特性を補正値として適用することで、意図した補正が可能である。しかしながら、回復フィルタを単純な逆特性としてしまうと、特に元の応答性の悪い箇所で、部分的に非常に高いゲインのシャープネスをかけたような画像となってしまい、ノイズや黒沈み等の弊害が強調されてしまう。そのため、画像回復処理は、一般的に強度をユーザに調整させ、それを元に回復フィルタのゲインをコントロールしている。
【０００９】
回復フィルタの設計については、様々な式が考えられるが、概ね、式（４）を式（５）のように考え、式（５）のＲ’（ｘ，ｙ）について、式（６）のように考える。
【００１０】
ｇ（ｘ，ｙ）＊Ｒ’（ｘ，ｙ）＝ｆ（ｘ，ｙ）      （５）
Ｋ（ｘ，ｙ，Ｒ，α，β）＊Ｒ（ｘ，ｙ）＝Ｒ’（ｘ，ｙ）      （６）
式（６）のＫは、Ｒに応じた関数として定義されているが、実際には複雑な関数を定義する必要がある。場合によっては不定値が複数存在し、入力Ｒから一意に決めることができない場合がある。その際には、何かしらの目標値を定め、不定値に値を代入し、最適化処理を用いてＫを決定する。例えば、Ｒ’の絶対値を回復フィルタ全体で計算し、最大値をゲインとして定義すると、ユーザに調整値として目標のゲインを入力してもらい、そのゲインを満たすようなＫを最適化処理によって算出し、回復フィルタを生成することができる。ゲインが高いとノイズや黒沈みを強調してしまうため、どの程度までゲインを上げることができるかは撮影された画像によって異なる。高感度で撮影された画像はノイズも多く、ゲインをあまり上げることができないが、低感度で撮影された画像は高感度の場合よりもゲインを高く設定することができる。また、撮影時の光学伝達関数の特性によって、あまりにも悪い特性であればゲインを上げると弊害が目立ち、特性がよければほぼ完全な状態までゲインを強く上げることができる。このように、回復フィルタを生成する場合にゲインを入力してもらう必要があり、それらはあらかじめ一律で設定しておくことは難しい。特許文献１では、ユーザが調整値を入力した上で、その値が明らかに強い場合には、弊害を抑えるために目標とするゲインを低く設定している。
【先行技術文献】
【特許文献】
【００１１】
【特許文献１】特許第４０１９５０２号
【発明の概要】
【発明が解決しようとする課題】
【００１２】
しかしながら、ユーザに調整値を入力してもらうことが難しい場合も考えられる。例えば、撮像装置であるカメラ本体が撮影と同時に回復処理を行う場合には、撮影ごとにユーザに調整値を入力してもらうことは応答性を著しく損ね、シャッターチャンスを逃すことにもなりかねない。
【００１３】
このような課題を鑑みて、本発明は、撮影時に自動で適切な調整値を決定し、回復フィルタを生成可能な画像処理装置を提供することを目的とする。
【課題を解決するための手段】
【００１４】
本発明の一側面としての画像処理装置は、撮影条件に応じた光学伝達関数を取得する取得手段と、前記光学伝達関数の特性に基づいて、前記光学伝達関数のゲインを決定する決定手段と、前記光学伝達関数および前記ゲインを用いて、画像回復処理を行うための回復フィルタを生成する生成手段と、を有することを特徴とする。
【発明の効果】
【００１５】
本発明によれば、撮影時に自動で適切な調整値を決定し、回復フィルタを生成可能な画像処理装置を提供することができる。
【図面の簡単な説明】
【００１６】
【図１】本発明の実施形態に係る撮像装置の構成図である。
【図２】画像処理装置のブロック図である。
【図３】光学伝達関数取得部の処理フローを示した図である。
【図４】目標ゲイン決定部の処理フローを示した図である。
【図５】５０％ゲイン判定の処理フローを示した図である。
【図６】５０％ゲイン判定の際の模式図である。
【図７】０落ちゲイン判定の処理フローを示した図である。
【図８】０落ちゲイン判定の際の模式図である。
【図９】光学伝達関数の微分値判定の処理フローを示した図である。
【図１０】光学伝達関数の微分値判定の際の模式図である。
【図１１】フィルタ生成部の処理フローを示した図である。
【図１２】フィルタ生成フロー中のフィルタコピーに関する模式図である。
【発明を実施するための形態】
【００１７】
以下、本発明の実施例について、図面を参照しながら詳細に説明する。各図において、同一の部材については同一の参照番号を付し、重複する説明は省略する。
【００１８】
まず、図１を用いて、撮像装置１００の構成について説明する。図１は、本発明の実施形態に係る撮像装置１００の構成図である。撮像素子１２は撮影レンズ１０によって形成された光学像を電気信号に変換し、Ａ／Ｄ変換器１４は撮像素子１２からのアナログ信号をディジタル信号に変換する。画像処理回路１６は、Ａ／Ｄ変換器１４からのデータあるいはメモリ制御回路２０からのデータに対して所定のデモザイク処理や色変換処理を行う。メモリ制御回路２０は、Ａ／Ｄ変換器１４、画像処理回路１６を制御する。Ａ／Ｄ変換器１４のデータは、画像処理回路１６、メモリ制御回路２０を介して、あるいは直接メモリ制御回路２０を介して、画像表示メモリ２６あるいはメモリ２８に書き込まれる。画像表示メモリ２６に書き込まれた表示用の画像データは、Ｄ／Ａ変換器２２を介してＴＦＴ  ＬＣＤ等から成る画像表示部２４に表示される。画像表示部２４に撮像した画像データを逐次表示すれば、ライブビュー機能を実現することが可能である。メモリ２８は、所定枚数の静止画像や所定時間の動画像を格納するのに十分な記憶量を備えている。これにより、複数枚の静止画像を連続して撮影する連写撮影の場合にも、高速かつ大量の画像書き込みをメモリ２８に対して行うことが可能となる。また、メモリ２８は、撮像装置１００全体を制御するシステム制御回路３２の作業領域としても使用することが可能である。シャッタースイッチＳＷ１は、シャッターボタン（不図示）の操作途中でＯＮとなり、ＡＦ（オートフォーカス）処理、ＡＥ（自動露出）処理、ＡＷＢ（オートホワイトバランス）処理、ＥＦ（フラッシュプリ発光）処理等の動作開始を指示する。シャッタースイッチＳＷ２は、シャッターボタンの完全な押下でＯＮとなる。このとき、シャッタースイッチＳＷ２は、画像処理回路１６やメモリ制御回路２０での現像処理、メモリ２８から画像データを読み出し、圧縮処理後、記録媒体を含む記録部３０に画像データを書き込む記録処理等の動作開始を指示する。
【００１９】
図２は、画像処理装置として機能する画像処理回路１６のブロック図である。画像処理回路１６は、光学伝達関数取得手段部（取得手段、以下、取得部という）１０１、目標ゲイン決定部（決定手段）１０２およびフィルタ生成部（生成手段）１０３を有する。取得部１０１は、撮影条件に応じた光学伝達関数を取得する。目標ゲイン決定部１０２は、取得した光学伝達関数の特性から目標ゲインを決定する。フィルタ生成部１０３は、目標ゲインとなるように回復フィルタを生成する。
【００２０】
図３は、取得部１０１の処理フローを示した図である。ステップＳ３００で処理フローを開始すると、ステップＳ３０１で取得部１０１は撮影条件を取得する。撮影条件は、撮像光学系や撮像素子の識別子と、焦点距離、絞り、撮影距離等である。ステップＳ３０２で、取得した撮影条件、おもに、レンズ、焦点距離、撮影距離を用いて、全てのレンズの光学伝達関数を保持している光学伝達関数係数データベースから光学伝達関数の係数を取得する。なお、光学伝達関数は、各撮影条件の像高ごとに複数用意されており、またＲＧＢのチャンネルごとにも異なる場合がある。そのため、各撮影条件で使用する全ての光学伝達関数を一括で取得する必要がある。ステップＳ３０３で取得した光学伝達関数の係数をフィルタ生成部１０３に送信し、処理フローを終了する（ステップＳ３０４）。
【００２１】
図４は、目標ゲイン決定部１０２の処理フローを示した図である。ステップＳ４００で処理を開始すると、ステップＳ４０１で取得部１０１において取得した光学伝達関数を１つずつ取り出す。ステップＳ４０２で、ベース目標ゲインを取得する。ベース目標ゲインは、撮影条件によらず定数でもよいし、撮影時のＩＳＯ感度に応じた値でも構わない。例えば、ＩＳＯ感度に応じてベース目標ゲインを用意する場合、感度が１０００までは５、それ以上の感度であれば３という簡易的な値とすればよい。後述するステップＳ４０３〜Ｓ４０５で、光学伝達関数の特性を判定する。これらの判定方法は、１つだけ使用してもよいし、全て使用してもよい。ステップＳ４０３〜Ｓ４０５のそれぞれの判定で、ゲイン調整割合が算出可能なので、それらを全て目標ベースゲインに乗算し（ステップＳ４０６）、その値を目標ゲインとする（ステップＳ４０７）。ステップＳ４０８で、全ての光学伝達関数に上述した処理が行われたかどうかを判断する。ステップＳ４０８で処理が行われていないと判断された場合、ステップＳ４０１に戻り、処理が行われていない光学伝達関数の目標ゲインを算出する。ステップＳ４０８で処理が行われていると判断された場合、処理フローを終了する（ステップＳ４０９）。
【００２２】
図５は図４のステップＳ４０３の５０％ゲイン判定の処理フローを示した図であり、図６は５０％ゲイン判定の際の模式図である。この判定方法では、ゲインがＤＣの５０％となる周波数が、ナイキスト周波数に対してどの位置にあるかということを判定する。ステップＳ５００で処理を開始すると、ステップＳ５０１で光学伝達関数の中心、つまりＤＣのゲインを取得する。ステップＳ５０２で、ＤＣのゲインの５０％の値（６０２）を計算する。ステップＳ５０３で、光学伝達関数の中心から２次元の軸上方向、計４方向の値をそれぞれ抽出する。ステップＳ５０４で、ステップＳ５０３で抽出した４方向の値から１方向の値（６０１）ずつ取り出す。ステップＳ５０５で、ステップＳ５０２で求めた値（６０２）の周波数（６０４）を取得する。ステップＳ５０６で、ステップＳ５０５で求めた周波数（６０４）と機種によって異なる撮像装置１００のナイキスト周波数（６０３）とを比較する。ステップＳ５０５で求めた周波数（６０４）のほうが大きい場合、ゲインの割合は１とする（ステップＳ４０７）。ナイキスト周波数（６０３）と等しいか、ナイキスト周波数（６０３）のほうが大きい場合、ステップＳ５０８で、ステップＳ５０５で求めた周波数（６０４）とナイキスト周波数（６０３）の８割の周波数とを比較する。ステップＳ５０５で求めた周波数（６０４）のほうが大きい場合、ゲインの割合は０．８とする（ステップＳ４０９）。ナイキスト周波数（６０３）の８割の周波数と等しいか、ナイキスト周波数（６０３）の８割の周波数のほうが大きい場合、ステップＳ５１０で、ステップＳ５０５で求めた周波数（６０４）とナイキスト周波数（６０３）の５割の周波数とを比較する。ステップＳ５０５で求めた周波数（６０４）のほうが大きい場合、ゲインの割合は０．５とする（ステップＳ５１１）。ナイキスト周波数（６０３）の５割の周波数と等しいか、ナイキスト周波数（６０３）の５割の周波数のほうが大きい場合、ゲインの割合は０．３とする（ステップＳ５１２）。なお、判定の数値や割合は、仮の値である。また、ＤＣからナイキスト周波数を１００とした場合の割合を、そのままゲインの割合として用いてもよい。ステップＳ５１３で、４方向でゲインの割合が求められたかどうかを判断する。ステップＳ５１４で４方向のゲインの割合の平均値を出力し、処理フローを終了する（ステップＳ５１５）。なお、本実施形態では、ステップＳ５１４で、各方向の平均を算出しているが、最大値や最小値、最大と最小を省いた平均値等を算出してもよい。
【００２３】
図７は図４のステップＳ４０４の０落ちゲイン判定の処理フローを示した図であり、図８は０落ちゲイン判定の際の模式図である。この判定方法では、ゲインが０となる周波数が、ナイキスト周波数に対してどの位置にあるかということを判定する。ステップＳ７００で処理を開始すると、ステップＳ７０１で０落ちゲイン（８０２）を取得する。０落ちゲイン（８０２）を厳密に０としてしまうと、ノイズ等の影響もあり、正常に判定できない可能性がある。そこで、０落ちゲイン（８０２）をノイズとみなされない程度の所定値に設定する必要がある。本実施形態では、０落ちゲイン（８０２）をＤＣのゲインの３％の値に設定している。ステップＳ７０２で、光学伝達関数の中心から２次元の軸上方向、計４方向の値をそれぞれ抽出する。ステップＳ７０３で、ステップＳ７０２で抽出した４方向の値から１方向の値（８０１）ずつ取り出す。ステップＳ７０４で、ＤＣ側から見てステップＳ７０１で取得した０落ちゲイン（８０２）よりも低くなる最初の周波数（８０４）を取得する。ステップＳ７０５で、ステップＳ７０４で求めた周波数（８０４）と撮像装置１００のナイキスト周波数（８０３）とを比較する。ステップＳ７０４で求めた周波数（８０４）のほうが大きい場合、ゲインの割合は１とする（ステップＳ７０６）。ナイキスト周波数（８０３）と等しいか、ナイキスト周波数（８０３）のほうが大きい場合、ステップＳ７０７で、ステップＳ７０４で求めた周波数（８０４）とナイキスト周波数（８０３）の８割の周波数とを比較する。ステップＳ７０４で求めた周波数（８０４）のほうが大きい場合、ゲインの割合は０．８とする（ステップＳ７０８）。ナイキスト周波数（８０３）の８割の周波数と等しいか、ナイキスト周波数（８０３）の８割の周波数のほうが大きい場合、ステップＳ７０９で、ステップＳ７０４で求めた周波数（８０４）とナイキスト周波数（８０３）の５割の周波数とを比較する。ステップＳ７０４で求めた周波数（８０４）のほうが大きい場合、ゲインの割合は０．５とする（ステップＳ７１０）。ナイキスト周波数（８０３）の５割の周波数と等しいか、ナイキスト周波数（８０３）の５割の周波数のほうが大きい場合、ゲインの割合は０．３とする（ステップＳ７１１）。なお、判定の数値や割合は、仮の値である。また、ＤＣからナイキスト周波数を１００とした場合の割合を、そのままゲインの割合として用いてもよい。ステップＳ７１２で、４方向でゲインの割合が求められたかどうかを判断する。ステップＳ７１３で４方向のゲインの割合の平均値を出力し、処理フローを終了する（ステップＳ７１４）。なお、本実施形態では、ステップＳ７１３で、各方向の平均を算出しているが、最大値や最小値、最大と最小を省いた平均値等を算出してもよい。
【００２４】
図９は図４のステップＳ４０５の光学伝達関数の微分値判定の処理フローを示した図であり、図１０は光学伝達関数の微分値判定の際の模式図である。この判定方法では、ＤＣからナイキスト周波数まで走査し、隣あう２点を抽出し、２点間の傾きを用いて判定する。通常の光学伝達関数であれば、２点間の傾きは常にナイキスト周波数に向かって連続的に減少していく。仮に２点間の傾きが増加している部分があれば、ノイズ等の弊害が出る可能性がある。増加している部分がある場合でも、その位置がナイキスト周波数に対して高周波域であれば、低周波域であるよりも特性がよいと考えられる。ステップＳ９００で処理を開始すると、ステップＳ９０１で光学伝達関数の中心から２次元の軸上方向、計４方向の値をそれぞれ抽出する。ステップＳ９０２で、ステップＳ９０１で抽出した４方向の値から１方向の値（１００１）ずつ取り出す。ステップＳ９０３で、ＤＣ側から２点を順番に取り出し、２点間の傾きが初めて増加する周波数（１００３）を取得する。これは、光学伝達関数をグラフとして考えた場合の微分値を求め、その値が正（増加）となる周波数を求めることと同義である。ステップＳ９０５で、ステップＳ９０４で求めた周波数（１００３）と撮像装置１００のナイキスト周波数（１００２）とを比較する。ステップＳ９０４で求めた周波数（１００３）がナイキスト周波数（１００２）内に存在しない場合、ゲインの割合は１とする（ステップＳ９０５）。ステップＳ９０４で求めた周波数（１００３）がナイキスト周波数（１００２）内に存在し、ステップＳ９０６で、ステップＳ９０４で求めた周波数（１００３）とナイキスト周波数（１００２）の８割の周波数とを比較する。ステップＳ９０４で求めた周波数（１００２）のほうが大きい場合、ゲインの割合は０．５とする（ステップＳ９０７）。ナイキスト周波数（１００２）の８割の周波数と等しいか、ナイキスト周波数（１００２）の８割の周波数のほうが大きい場合、ゲインの割合は０．３とする（ステップＳ９０８）。なお、判定の数値や割合は、仮の値である。また、ＤＣからナイキスト周波数を１００とした場合の割合を、そのままゲインの割合として用いてもよい。ステップＳ９０９で、４方向でゲインの割合が求められたかどうかを判断する。ステップＳ９１０で４方向のゲインの割合の平均値を出力し、処理フローを終了する（ステップＳ９１１）。なお、本実施形態では、ステップＳ９１０で、各方向の平均を算出しているが、最大値や最小値、最大と最小を省いた平均値等を算出してもよい。
【００２５】
図１１はフィルタ生成部１０３の処理フローを示した図であり、図１２はフィルタ生成フロー中のフィルタコピーに関する模式図である。ステップＳ１１００で処理を開始すると、ステップＳ１１０１で取得部１０１から光学伝達関数の係数を取得する。なお、全画素の光学伝達関数のデータを取得してもよいが、現実的にはデータ量や処理速度との兼ね合いで、画像領域内部の離散点においてデータを保持し、実際の画像に適用する際に近傍の回復フィルタを補間して用いるほうがよい。処理速度の懸念がある場合には、全面ではなく、画面の一部分の回復フィルタを生成し、適用するときに補間して使う場合もある。ステップＳ１１０２で、画像の中心から周辺までの像高を１０分割した離散点を取得し、その１０像高分の光学伝達関数を生成する。ステップＳ１１０３で、その１０像高分の光学伝達関数を回転させ、全面の４分の１の領域について離散点の光学伝達関数を生成する。ステップＳ１１０４で、ナイキスト周波数で光学伝達関数を切り出す。ステップＳ１１０５で、撮像装置の特性である光学ローパスフィルタや開口劣化特性を掛け合わせる。ステップＳ１１０６で、目標ゲイン決定部１０２から回復フィルタの生成に必要な目標ゲインを取得する。ステップＳ１１０７で取得した目標ゲインを加味して、光学伝達関数の逆特性を求め、ステップＳ１１０８でその逆特性のデータをフーリエ変換することで、回復フィルタが求まる。光学伝達関数のデータは複素数データであるが、回復フィルタは実数の二次元フィルタとして生成される。図１２（ａ）に示すように、画面の４分の１の範囲で演算していた場合には、ステップＳ１１０９でさらに全面になるようにデータをコピーして、図１２（ｂ）に示すように、全面の回復フィルタを生成して処理フローを終了する（ステップＳ１１１０）。生成した回復フィルタは、画像の各画素に適用されることで画像回復処理が行われる。実際には、画像全体に適用する場合もあれば、拡大表示をするために画像の一部分に適用する場合もある。また、回復フィルタは、画像の全ての画素に対して生成せず、いくつかの離散点に対して生成することが一般的である。そのため、フィルタ生成部１０３で生成した回復フィルタから対象の画素位置における回復フィルタを補間等の方法を用いて算出し、その結果を適用してもよい。回復フィルタは、対象画素１つにつき、幅と高さを持つ二次元のデータ配列になっており、対象画素を含む領域に対して、回復フィルタの各値をそれぞれの画素について求める。
【００２６】
また、フィルタ生成部１０３において、毎回の撮影条件ごとに回復フィルタを作らずに、一度生成した回復フィルタを装置内の不図示のメモリにキャッシュしておいてもよい。そうすることで、再度同じ撮影条件で回復フィルタを生成する必要があった際に、回復フィルタを生成せずに、キャッシュをそのまま利用することが可能となる。このとき、装置内のメモリ（容量）には制限があるため、古いものからキャッシュを削除していくという基本的な仕組みに加え、より効率的になるように回復フィルタに優先度付けを行い、削除する回復フィルタを選択してもよい。優先度の付け方は複数考えられ、例えば、使用回数をインクリメントし保存しておき、使用回数の少ない順から削除する方法や、データサイズの大きさを考慮して削除する方法が考えられる。
【００２７】
以上、本発明の好ましい実施形態について説明したが、本発明はこれらの実施形態に限定されず、その要旨の範囲内で種々の変形及び変更が可能である。
【符号の説明】
【００２８】
１０１  光学伝達関数取得部（取得手段）
１０２  目標ゲイン決定部（決定手段）
１０３  フィルタ生成部（生成手段）
【図１】
2015226086.tif 000003
【図２】
2015226086.tif 000004
【図３】
2015226086.tif 000005
【図４】
2015226086.tif 000006
【図５】
2015226086.tif 000007
【図６】
2015226086.tif 000008
【図７】
2015226086.tif 000009
【図８】
2015226086.tif 000010
【図９】
2015226086.tif 000011
【図１０】
2015226086.tif 000012
【図１１】
2015226086.tif 000013
【図１２】
2015226086.tif 000014