"撮影画像の画像回復処理を行う画像処理方法であって、前記撮影画像の撮影条件に応じた光学伝達関数を取得する工程と、前記光学伝達関数の絶対値が所定の閾値以下になる閾値周波数を取得する工程と、前記光学伝達関数を用いて、前記閾値周波数に応じたゲイン特性を有する画像回復フィルタを生成する工程と、前記画像回復フィルタを用いて前記撮影画像の画像回復処理を行う工程と、を有することを特徴とする画像処理方法。前記閾値周波数は、前記光学伝達関数において前記絶対値が前記所定の閾値以下となる周波数のうち最も低い周波数であることを特徴とする請求項1に記載の画像処理方法。前記閾値周波数は、撮影時に検出されることを特徴とする請求項1または2に記載の画像処理装置。前記閾値周波数は、撮影時のFナンバーに基づく周波数であることを特徴とする請求項1または2に記載の画像処理方法。前記画像回復フィルタの前記ゲイン特性は、前記閾値周波数が低くなるにつれて最大ゲインが低下するように変化することを特徴とする請求項1乃至4のいずれか1項に記載の画像処理方法。前記画像回復フィルタの前記ゲイン特性は、前記閾値周波数が低くなるにつれて前記画像回復フィルタのゲインを低周波側に高い重みを付けるように変化することを特徴とする請求項1乃至4のいずれか1項に記載の画像処理方法。前記閾値周波数は、前記撮影画像の位置に応じて異なり、前記画像回復フィルタの前記ゲイン特性は、前記撮影画像の位置に応じて異なることを特徴とする請求項1乃至6のいずれか1項に記載の画像処理方法。撮影画像の画像回復処理を行う画像処理装置であって、前記撮影画像の撮影条件に応じた光学伝達関数を取得する光学伝達関数取得手段と、前記光学伝達関数の絶対値が所定の閾値以下になる閾値周波数を取得する閾値周波数取得手段と、前記光学伝達関数を用いて、前記閾値周波数に応じたゲイン特性を有する画像回復フィルタを生成する生成手段と、前記画像回復フィルタを用いて前記撮影画像の画像回復処理を行う処理手段と、を有することを特徴とする画像処理装置。撮影画像の画像回復処理を行う撮像装置であって、撮像光学系と、前記撮像光学系を介して得られた被写体像を光電変換して撮影画像を生成する撮像素子と、前記撮影画像の撮影条件に応じた光学伝達関数を取得する光学伝達関数取得手段と、前記光学伝達関数の絶対値が所定の閾値以下になる閾値周波数を取得する閾値周波数取得手段と、前記光学伝達関数を用いて、前記閾値周波数に応じたゲイン特性を有する画像回復フィルタを生成する生成手段と、前記画像回復フィルタを用いて前記撮影画像の画像回復処理を行う処理手段と、を有することを特徴とする撮像装置。コンピュータに撮影画像の画像回復処理を実行させる画像処理プログラムであって、前記撮影画像の撮影条件に応じた光学伝達関数を取得する工程と、前記光学伝達関数の絶対値が所定の閾値以下になる閾値周波数を取得する工程と、前記光学伝達関数を用いて、前記閾値周波数に応じたゲイン特性を有する画像回復フィルタを生成する工程と、前記画像回復フィルタを用いて前記撮影画像の画像回復処理を行う工程と、を有することを特徴とする画像処理プログラム。しかしながら、撮影条件によっては、撮影画像の少なくとも一部で、撮像素子のナイキスト周波数の帯域内で光学伝達関数の振幅成分MTFが零になる場合がある。以下、これを零落ちと呼び、零落ちする周波数を零落ち周波数という。零落ちの原因としては、収差や回折、手ぶれなどが挙げられる。零落ち周波数の情報は伝達撮像されないため、零落ち周波数での画像回復処理を行うことができない。また、零落ち周波数以外の周波数のみで画像回復処理を行うと、画像のエッジ部にリンギングが発生してしまう。この結果、高品質な回復画像を得ることができない。特許文献1の構成において、回復度合を調節することは可能である。しかし特許文献1の構成は、ユーザによる回復度合の変更が可能であるに過ぎず、光学系光学伝達関数の周波数特性に応じた適切なパラメータを自動的に決定するものではない。また、零落ち周波数に応じて回復度合を調節するものではないため、零落ちが存在する光学系による差撮影画像を良好に回復することはできない。そこで本発明は、良好な画像回復処理が可能な画像処理方法、画像処理装置、撮像装置、および、画像処理プログラムを提供する。本発明の一側面としての画像処理方法は、撮影画像の画像回復処理を行う画像処理方法であって、前記撮影画像の撮影条件に応じた光学伝達関数を取得する工程と、前記光学伝達関数の絶対値が所定の閾値以下になる閾値周波数を取得する工程と、前記光学伝達関数を用いて、前記閾値周波数に応じたゲイン特性を有する画像回復フィルタを生成する工程と、前記画像回復フィルタを用いて前記撮影画像の画像回復処理を行う工程とを有する。本発明の他の側面としての画像処理装置は、撮影画像の画像回復処理を行う画像処理装置であって、前記撮影画像の撮影条件に応じた光学伝達関数を取得する光学伝達関数取得手段と、前記光学伝達関数の絶対値が所定の閾値以下になる閾値周波数を取得する閾値周波数取得手段と、前記光学伝達関数を用いて、前記閾値周波数に応じたゲイン特性を有する画像回復フィルタを生成する生成手段と、前記画像回復フィルタを用いて前記撮影画像の画像回復処理を行う処理手段とを有する。本発明の他の側面としての撮像装置は、撮影画像の画像回復処理を行う撮像装置であって、撮像光学系と、前記撮像光学系を介して得られた被写体像を光電変換して撮影画像を生成する撮像素子と、前記撮影画像の撮影条件に応じた光学伝達関数を取得する光学伝達関数取得手段と、前記光学伝達関数の絶対値が所定の閾値以下になる閾値周波数を取得する閾値周波数取得手段と、前記光学伝達関数を用いて、前記閾値周波数に応じたゲイン特性を有する画像回復フィルタを生成する生成手段と、前記画像回復フィルタを用いて前記撮影画像の画像回復処理を行う処理手段とを有する。本発明の他の側面としての画像処理プログラムは、コンピュータに撮影画像の画像回復処理を実行させる画像処理プログラムであって、前記撮影画像の撮影条件に応じた光学伝達関数を取得する工程と、前記光学伝達関数の絶対値が所定の閾値以下になる閾値周波数を取得する工程と、前記光学伝達関数を用いて、前記閾値周波数に応じたゲイン特性を有する画像回復フィルタを生成する工程と、前記画像回復フィルタを用いて前記撮影画像の画像回復処理を行う工程とを有する。本発明の他の目的及び特徴は、以下の実施例において説明される。本発明によれば、良好な画像回復処理が可能な画像処理方法、画像処理装置、撮像装置、および、画像処理プログラムを提供することができる。以下、本発明の実施例について、図面を参照しながら詳細に説明する。各図において、同一の部材については同一の参照番号を付し、重複する説明は省略する。まず、本実施例で説明される用語の定義および画像回復処理画像処理方法について説明する。ここで説明される画像処理方法は、後述の各実施例において適宜用いられる。[入力画像]入力画像は、撮像光学系を介して撮像素子で受光することで得られたデジタル画像撮影画像であり、レンズと各種の光学フィルタ類を含む撮像光学系の収差による光学伝達関数OTFにより劣化している。撮像光学系は、レンズだけでなく曲率を有するミラー反射面を用いて構成することもできる。入力画像の色成分は、例えばRGB色成分の情報を有する。色成分としては、これ以外にもLCHで表現される明度、色相、彩度や、YCbCrで表現される輝度、色差信号など一般に用いられている色空間を選択して用いることができる。その他の色空間として、XYZ、Lab、Yuv、JChを用いることが可能である。更に、色温度を用いてもよい。入力画像や出力画像には、レンズの焦点距離、絞り値、撮影距離などの撮影条件や、この画像を補正するための各種の補正情報を付帯することができる。撮像装置から別の画像処理装置に画像を受け渡して補正処理を行う場合、上述のように撮影画像に撮影条件や補正に関する情報を付帯することが好ましい。撮影条件や補正に関する情報の他の受け渡し方法として、撮像装置と画像処理装置を直接または間接的に接続して受け渡すようにしてもよい。[画像回復処理]続いて、画像回復処理の概要について説明する。撮影画像劣化画像をgx,y、もとの画像をfx,y、光学伝達関数OTFのフーリエペアである点像分布関数PSFをhx,yとしたとき、以下の式1が成立する。gx,y=hx,yfx,y...1ここで、はコンボリューション畳み込み積分、積和、x,yは撮影画像上の座標である。また、式1をフーリエ変換して周波数面での表示形式に変換すると、周波数ごとの積で表される式2が得られる。Gu,v=Hu,vFu,v...2ここで、Hは点像分布関数PSFhをフーリエ変換することにより得られた光学伝達関数OTFであり、G,Fはそれぞれ劣化した画像g、もとの画像fをフーリエ変換して得られた関数である。u,vは2次元周波数面での座標、すなわち周波数である。撮影された劣化画像gからもとの画像fを得るには、以下の式3のように両辺を光学伝達関数Hで除算すればよい。Gu,v/Hu,v=Fu,v...3そして、Fu,v、すなわちGu,v/Hu,vを逆フーリエ変換して実面に戻すことにより、もとの画像fx,yが回復画像として得られる。H-1を逆フーリエ変換したものをRとすると、以下の式4のように実面での画像に対するコンボリューション処理を行うことで、同様にもとの画像fx,yを得ることができる。gx,yRx,y=fx,y...4ここで、Rx,yは画像回復フィルタと呼ばれる。画像が2次元画像である場合、一般的に、画像回復フィルタRも画像の各画素に対応したタップセルを有する2次元フィルタとなる。また、画像回復フィルタRのタップ数セルの数は、一般的に多いほど回復精度が向上する。このため、要求画質、画像処理能力、収差の特性等に応じて実現可能なタップ数が設定される。画像回復フィルタRは、少なくとも収差の特性を反映している必要があるため、従来の水平垂直各3タップ程度のエッジ強調フィルタハイパスフィルタなどとは異なる。画像回復フィルタRは光学伝達関数OTFに基づいて設定されるため、振幅成分および位相成分の劣化の両方を高精度に補正することができる。また、実際の画像にはノイズ成分が含まれるため、上記のように光学伝達関数OTFの逆数をとって作成した画像回復フィルタRを用いると、劣化画像の回復とともにノイズ成分が大幅に増幅されてしまう。これは、画像の振幅成分にノイズの振幅が付加されている状態に対して、光学系のMTF振幅成分を全周波数に渡って1に戻すようにMTFを持ち上げるためである。光学系による振幅劣化であるMTFは1に戻るが、同時にノイズのパワースペクトルも持ち上がってしまい、結果的にMTFを持ち上げる度合回復ゲインに応じてノイズが増幅されてしまう。したがって、ノイズが含まれる場合には、鑑賞用画像としては良好な画像は得られない。このことは、以下の式5-1、5-2で表される。Gu,v=Hu,vFu,vNu,v...5-1Gu,v/Hu,v=Fu,vNu,v/Hu,v...5-2ここで、Nはノイズ成分である。ノイズ成分が含まれる画像に関しては、例えば以下の式6で表されるウィナーフィルタのように、画像信号とノイズ信号の強度比SNRに応じて回復度合を制御する方法がある。ここで、Mu,vはウィナーフィルタの周波数特性、|Hu,v|は光学伝達関数OTFの絶対値振幅成分MTFである。この方法では、周波数ごとに、振幅成分MTFが小さいほど回復ゲイン回復度合を小さくし、振幅成分MTFが大きいほど回復ゲインを大きくする。一般的に、撮像光学系の振幅成分MTFは低周波側が高く高周波側が低くなるため、この方法では、実質的に画像の高周波側の回復ゲインを低減することになる。続いて、図2および図3を参照して、画像回復フィルタについて説明する。画像回復フィルタは、撮像光学系の収差特性や要求される回復精度に応じてそのタップ数が決定される。図2の画像回復フィルタは、一例として、11×11タップの2次元フィルタである。また図2では、各タップ内の値係数を省略しているが、この画像回復フィルタの一断面を図3に示す。画像回復フィルタの各タップの値係数値の分布は、収差により空間的に広がった信号値PSFを、理想的には元の1点に戻す機能を有する。画像回復フィルタの各タップは、画像の各画素に対応して画像回復処理の工程でコンボリューション処理畳み込み積分、積和される。コンボリューション処理では、所定の画素の信号値を改善するために、その画素を画像回復フィルタの中心と一致させる。そして、画像と画像回復フィルタの対応画素ごとに画像の信号値とフィルタの係数値の積をとり、その総和を中心画素の信号値として置き換える。続いて、図4および図5を参照して、画像回復の実空間と周波数空間での特性について説明する。図4は、点像分布関数PSFの説明図であり、図4aは画像回復前の点像分布関数PSF、図4bは画像回復後の点像分布関数PSFを示している。図5は、光学伝達関数OTFの振幅成分MTF図5aと位相成分PTF図5bの説明図である。図5a中の破線Aは画像回復前のMTF、一点鎖線Bは画像回復後のMTFを示す。また図5b中の破線Aは画像回復前のPTF、一点鎖線Bは画像回復後のPTFを示す。図4aに示されるように、画像回復前の点像分布関数PSFは、非対称な広がりを有し、この非対称性により位相成分PTFは周波数に対して非直線的な値を有する。画像回復処理は、振幅成分MTFを増幅し、位相成分PTFがゼロになるように補正するため、画像回復後の点像分布関数PSFは対称で先鋭な形状になる。このように画像回復フィルタは、撮像光学系の光学伝達関数OTFの逆関数に基づいて設計された関数を逆フーリエ変換して得ることができる。本実施例で用いられる画像回復フィルタは適宜変更可能であり、例えば上述のようなウィナーフィルタを用いることができる。ウィナーフィルタを用いる場合、式6を逆フーリエ変換することで、実際に画像に畳み込む実空間の画像回復フィルタを作成することが可能である。また、光学伝達関数OTFは1つの撮影状態においても撮像光学系の像高画像の位置に応じて変化する。このため、画像回復フィルタは像高に応じて変更して用いられる。次に、図1を参照して、本発明の実施例1における画像処理方法について説明する。図1は、本実施例における画像処理方法画像処理プログラムのフローチャートである。図1のフローチャートは、後述の画像処理装置の指令に基づいて実行される。まずステップS11において、撮影画像を入力画像として取得する。撮影画像は、撮像装置と画像処理装置を有線または無線で接続して取得することができる。また撮影画像は、記憶媒体を介して取得することもできる。続いてステップS12において、撮影画像の撮影条件を取得する。ステップS12は、画像処理装置の撮影条件取得手段により実行される。撮影条件は、焦点距離、絞り値、および、撮影距離等である。また、レンズがカメラ本体に交換可能に装着される撮像装置の場合、撮像条件は更にレンズIDやカメラIDを含む。撮影条件に関する情報は、撮像装置から直接取得することができる。またこの情報は、画像に付帯された情報から取得することもできる。次にステップS13において、撮影条件に応じた光学伝達関数OTFを取得する。ステップS13は、画像処理装置の光学伝達関数取得手段により実行される。この光学伝達関数OTFは、撮影条件に適した光学伝達関数OTFであり、予め保持された複数の光学伝達関数OTFの中から選択することができる。また、光学伝達関数OTFを生成するための関数とその生成に用いられる係数群を予め保持しておき、撮影条件に適した光学伝達関数OTFを本ステップにて新たに生成することもできる。更に、絞り、撮影距離、および、ズームレンズの焦点距離などの撮影条件が特定の撮影条件の場合、予め保持されている撮影条件に対応する光学伝達関数OTFから補間処理により生成することもできる。この場合、保持する画像回復フィルタのデータ量を低減することが可能である。補間処理としては、例えばバイリニア補間線形補間やバイキュービック補間等が用いられるが、これに限定されるものではない。続いてステップS14において、撮像光学系の光学伝達関数OTFの絶対値振幅成分MTFが所定の閾値以下となる閾値周波数を取得する。ステップS14は、画像処理装置の閾値周波数取得手段により実行される。閾値周波数は、撮影時に検出することができるが、これに限定されるものではなく、例えば撮影時のFナンバーに基づいて取得することもできる。閾値周波数の詳細は後述する。続いてステップS15において、撮像光学系の光学伝達関数OTFを用いて、閾値周波数に応じたゲイン特性を有する画像回復フィルタを生成する。ステップS15は、画像処理装置の生成手段により実行される。ステップS15にて生成された画像回復フィルタは、閾値周波数という観点から光学伝達関数OTFの特性を考慮して得られた画像回復フィルタである。画像回復フィルタの生成の詳細については後述する。なお、撮像光学系の光学伝達関数OTFを元に、閾値周波数に応じたゲイン特性を有する第2の光学伝達関数OTFを生成し、第2の光学伝達関数OTFに基づいて画像回復フィルタを生成してもよい。このようにして生成された画像回復フィルタも、閾値周波数に応じたゲイン特性を有する画像回復フィルタとなる。続いてステップS16において、ステップS15で生成された画像回復フィルタを用いて撮影画像の画像回復処理を行う。すなわち、撮影画像に画像回復フィルタをコンボリューション畳み込みすることで、撮影画像の画像回復処理を行う。ステップS16は、画像処理装置の処理手段により実行される。そしてステップS17において、ステップS16での画像回復処理の結果に基づいて回復画像出力画像を出力し、図1のフローは終了する。次に、図6を参照して、図1中のステップS14における閾値周波数を取得する方法について詳述する。図6は、本実施例における閾値周波数の説明図であり、光学伝達関数の振幅成分MTF絶対値と空間周波数との関係周波数特性を示している。ここでは、横軸空間周波数の1を撮像素子の画素ピッチで決定されるナイキスト周波数として正規化している。図6aは、互いに異なる3つの周波数特性A、B、Cを有する光学伝達関数の振幅成分MTFを示している。図6aでは、周波数特性A、Bのそれぞれの場合の振幅成分MTFが閾値t1となるときの周波数fa、fbを閾値周波数として設定する。周波数特性Cについては、閾値t1以下となる周波数閾値周波数は存在しない。一方、図6bは、互いに異なる2つの周波数特性D、Eを有する光学伝達関数の振幅成分MTFを示している。図6bでは、周波数特性D、Eのそれぞれの場合の振幅成分MTFが閾値t2となるときの周波数fd、feを閾値周波数として設定する。このように、光学伝達関数の振幅成分MTFの所定の閾値tt1、t2は目的に応じて適切な値に設定することが可能である。例えば、閾値周波数を零落ち周波数として取得する場合、振幅成分MTFが0になるときt=0の空間周波数を閾値周波数に設定すればよい。ただし、実用的には振幅成分MTFが0でない場合でも十分小さな値として扱うことが可能である。このとき、例えばt=0.03MTFが3の場合の空間周波数を閾値周波数として設定することができる。また図6bでは、t=0.50MTFが50の場合の空間周波数を閾値周波数として設定している。この場合、零落ち周波数の取得ではなく、基本性能を評価している意味合いが強い。また、偽解像がある場合など振幅成分MTFが0近傍で波打つような場合、t=0になる周波数が複数存在することになる。本実施例では、このような複数の周波数振幅成分MTFが所定の閾値以下となる複数の周波数のうち最も低い周波数を閾値周波数として設定することが好ましい。閾値周波数の具体的な取得方法としては、2次元のMTFデータを1つずつ所定の閾値と比較し、所定の閾値以下で且つ暫定の閾値周波数以下の場合にはその周波数を新たに暫定の閾値周波数として記憶する。これを繰り返すことで、所定の閾値以下となる最も低い周波数を閾値周波数として取得することができる。2次元のMTFデータは、保持している2次元のOTFデータの絶対値をとることで生成することができる。また変形例として、アジムス方向も記憶情報に加えることで、アジムス方向ごとに閾値周波数を取得することもできる。さらに、MTFデータのうち、特定のアジムス方向について1次元データとして閾値周波数を取得することもできる。閾値周波数は、全アジムス方向を対象に、所定の閾値以下となる最も低い周波数としてもよく、また、詳細にアジムス方向ごとに取得してもよい。なお、閾値周波数零落ち周波数は、図1中のステップS12にて取得した撮影時のFナンバーを用いて、回折限界の理論式に基づいて算出することもできる。また、光学伝達関数OTFがアジムス方向によって異なる場合、閾値周波数をアジムス方向に依存した特性を有する関数やルックアップテーブルのデータとすることもできる。または、全アジムス方向のうち、最も低い周波数を代表値として閾値周波数に設定することもできる。次に図7を参照して、図1中のステップS15で生成される画像回復フィルタの回復度合について説明する。画像回復フィルタを設計するための関数として、例えば上述の式6が用いられる。図7は、本実施例における画像回復フィルタの回復度合GAINを設定する際の説明図であり、図1中のステップS14で得られた閾値周波数から回復度合を決定するための関数を示している。図7中の2つの曲線2、3は、基本の回復度合の強弱を示す例である。このような回復度合の強弱は、プリセット値として記憶部に複数保持しておく。即ち、記憶部には、図7に示すような、閾値周波数と回復度合の対応を示す情報が格納されている。なお、プリセット値を所定の範囲においてユーザが任意に設定できるよう構成してもよい。図7中の曲線2が適用される場合、図6aの周波数特性A、B、Cを例とすると、周波数特性A、Bにおける回復度合は、それぞれの閾値周波数fa、fbに応じた回復度合Ga、Gbに設定される。周波数特性Cは零落ちが無い閾値周波数が存在しないため、その回復度合は曲線2の最大値である2に設定される。曲線3が適用される場合も同様に、それぞれの閾値周波数に応じた回復度合いに設定される。このような関数は、ルックアップテーブルで保持することもできる。そして、式6のSNRの項をパラメータとして用いて、画像回復フィルタの周波数特性が所望の回復度合となるように設計する。本実施例において、画像回復フィルタの回復度合GAINを表す値として、周波数特性の絶対値回復ゲインの最大値最大ゲインを用いる。図8は、本実施例における回復ゲインの説明図であり、最大ゲインがGa、Gb、2となるように設計された画像回復フィルタのそれぞれのゲイン特性f1、f2、f3を示している。ゲイン特性とは、空間周波数に対する回復度合を示す特性である。図8は、一例として、図7の曲線2を適用した場合のゲイン特性を示しており、ゲイン特性f1、f2、f3は、閾値周波数がfa、fbの場合、および、閾値周波数が存在しない場合にそれぞれ対応している。図7および図8に示されるように、閾値周波数fa、fbのそれぞれの場合のゲイン特性f1、f2で比較すると、閾値周波数fbのときの最大ゲインGbよりも閾値周波数faのときの最大ゲインGaのほうが小さい。すなわち、画像回復フィルタのゲイン特性は、閾値周波数が低くなるにつれて最大ゲインが低下するように変化する。また図8に示されるように、閾値周波数fa、fbのそれぞれの場合のゲイン特性f1、f2で比較すると、ゲイン特性f1のゲイン回復度合のほうがゲイン特性f2のゲインよりも低周波側に重み付けが行われている。すなわち、画像回復フィルタのゲイン特性は、閾値周波数が低くなるにつれて画像回復フィルタのゲインを低周波側に高い重みを付けるように変化する。また閾値周波数は、撮影画像の位置に応じて異なる。このため、画像回復フィルタのゲイン特性は、撮影画像の位置に応じて異なっている。図9は、画像回復前後の光学伝達関数の振幅成分MTFの説明図である。図9aは、図6a中の周波数特性A、B、Cに対する本実施例の画像回復処理を行った結果である。画像回復前の周波数特性A、B、Cに対して、図8に示されるようなゲイン特性f1、f2、f3を有する画像回復フィルタを適用することにより、画像回復後の周波数特性A1、B1、C1にそれぞれ変化する。一方、図9bは、光学伝達関数の振幅成分MTFの周波数特性閾値周波数を考慮せずに画像回復処理を行った結果を示している。画像回復前の周波数特性A、B、Cに対して、画像回復後の周波数特性A2、B2、C2にそれぞれ変化する。図9aでは、閾値周波数に応じた適切なゲイン特性を有する画像回復フィルタを用いて画像回復処理を行っているためリンギングの発生は低減する。一方、図9bではリンギングが発生し、高品質な回復画像を得ることができない。さらに、複数の閾値周波数を同時に用いて第一の閾値周波数および第二の閾値周波数により複合的に光学伝達関数の周波数特性ゲイン特性を分析することも可能である。複数の閾値周波数を用いる例として、式7を示す。式7は、式6に対してCu,vという周波数の重み付けを行って得られた式である。例えば、第一の閾値をt=0.03として第一の閾値周波数を取得し、第二の閾値をt=0.50として第二の閾値周波数を取得する。そして、第一の閾値周波数よりも低周波側のみを増幅するようにCu,vを設定することで、零落ち周波数に増幅ゲインが適用されないようにし、第二の閾値周波数で最大ゲインを設定して回復度合を決定することができる。すなわち、回復ゲインの周波数帯域の重み付けおよび回復度合の両方を、取得した2つの閾値周波数から自動的に決定することが可能になる。次に、図10を参照して、本実施例における撮像装置について説明する。図10は、本実施例における撮像装置200の構成図である。撮像装置200には、撮影画像の画像回復処理上述の画像処理方法を行う画像処理プログラムがインストールされており、この画像回復処理は撮像装置200の内部の画像処理部204画像処理装置により実行される。撮像装置200は、撮像光学系201レンズおよび撮像装置本体カメラ本体を備えて構成されている。撮像光学系201は、絞り201aおよびフォーカスレンズ201bを備え、撮像装置本体カメラ本体と一体的に構成されている。ただし本実施例はこれに限定されるものではなく、撮像光学系201が撮像装置本体に対して交換可能に装着される撮像装置にも適用可能である。撮像素子202は、撮像光学系201を介して得られた被写体像結像光を光電変換して撮影画像を生成する。すなわち被写体像は、撮像素子202により光電変換が行われてアナログ信号電気信号に変換される。そして、このアナログ信号はA/Dコンバータ203によりデジタル信号に変換され、このデジタル信号は画像処理部204に入力される。画像処理部204画像処理装置は、このデジタル信号に対して所定の処理を行うとともに、上述の画像回復処理を行う。まず画像処理部204撮影条件取得手段は、状態検知部207から撮像装置の撮像条件情報を取得する。撮像条件情報とは、絞り、撮影距離、または、ズームレンズの焦点距離等に関する情報である。状態検知部207は、システムコントローラ210から直接に撮像条件情報を取得することができるが、これに限定されるものではない。例えば撮像光学系201に関する撮像条件情報は、撮像光学系制御部206から取得することもできる。本実施例の画像回復処理の処理フロー画像処理方法は、図1を参照して説明したとおりである。光学伝達関数OTFまたは光学伝達関数OTFの生成に必要な係数データは、記憶部208に保持されている。画像処理部204で処理した出力画像は、画像記録媒体209に所定のフォーマットで保存される。表示部205には、本実施例の画像回復処理を行った画像に表示用の所定の処理を行った画像が表示される。ただしこれに限定されるものではなく、高速表示のために簡易処理を行った画像を表示部205に表示するように構成してもよい。本実施例における一連の制御はシステムコントローラ210により行われ、撮像光学系201の機械的な駆動はシステムコントローラ210の指示に基づいて撮像光学系制御部206により行われる。撮像光学系制御部206は、Fナンバーの撮影状態設定として、絞り201aの開口径を制御する。また撮像光学系制御部206は、被写体距離に応じてピント調整を行うため、不図示のオートフォーカスAF機構や手動のマニュアルフォーカス機構により、フォーカスレンズ201bの位置を制御する。なお、絞り201aの開口径制御やマニュアルフォーカスなどの機能は、撮像装置200の仕様に応じて実行しなくてもよい。撮像光学系201には、ローパスフィルタや赤外線カットフィルタ等の光学素子を入れても構わないが、ローパスフィルタ等の光学伝達関数OTFの特性に影響を与える素子を用いる場合、画像回復フィルタを作成する時点での考慮が必要になる場合がある。赤外カットフィルタに関しても、分光波長の点像分布関数PSFの積分値であるRGBチャンネルの各PSF、特にRチャンネルのPSFに影響するため、画像回復フィルタを作成する時点での考慮が必要になる場合がある。次に、図11を参照して、本発明の実施例2における画像処理装置および画像処理システムについて説明する。図11は、本実施例における画像処理システム300の構成図である。なお、本実施例の画像回復処理の処理フロー画像処理方法は、図1を参照して説明した実施例1と同様であるため、その説明は省略する。図11において、画像処理装置301は、本実施例の画像処理方法をコンピュータに実行させるための画像処理ソフトウェア306を搭載したコンピュータ機器である。撮像機器302は、カメラ、顕微鏡、内視鏡、または、スキャナなどの撮像装置である。記憶媒体303は、半導体メモリ、ハードディスク、または、ネットワーク上のサーバなど、撮影画像を記憶した記憶手段である。画像処理装置301は、撮像機器302または記憶媒体303から撮影画像データを取得し、所定の画像処理を行った画像データを出力機器305、撮像機器302、記憶媒体303のいずれか一つまたは複数に出力する。また、その出力先を画像処理装置301に内蔵された記憶部に保存することもできる。出力機器305は、例えばプリンタである。画像処理装置301にはモニタである表示機器304が接続されている。このため、ユーザは表示機器304を通して画像処理作業を行うとともに、補正された画像を評価することができる。画像処理ソフトウェア306は、本実施例の画像回復処理画像処理方法を行うほか、必要に応じて現像やその他の画像処理を行う。なお、本実施例における画像処理を行うためのデータの内容や機器間での受け渡しなどに関する情報補正情報については、個々の画像データに付帯させることが好ましい。必要な補正情報を画像データに付帯させることで、本実施例の画像処理装置を搭載した機器であれば、適切に本実施例の補正処理を行うことが可能である。上記各実施例によれば、撮影条件での光学伝達関数の特性を考慮して適切な回復度合ゲイン特性の回復フィルタを生成して用いることができるため、リンギング等の影響を低減して高品質な画像を得ることが可能である。従って、各実施例によれば、良好な画像回復処理が可能な画像処理方法、画像処理装置、撮像装置、および、画像処理プログラムを提供することができる。以上、本発明の好ましい実施形態について説明したが、本発明はこれらの実施形態に限定されず、その要旨の範囲内で種々の変形及び変更が可能である。"