"撮像装置により生成された入力画像に対して画像回復処理を行う画像処理方法であって、前記入力画像および該入力画像を生成した際の前記撮像装置における撮影条件に関する情報を得るステップと、前記撮影条件に関する情報に応じた前記撮像装置の光学伝達関数を得るステップと、前記光学伝達関数を用いて得られる指標値が所定値となる特定周波数をアジムス方向ごとに求めるステップと、前記アジムス方向ごとに求められた前記特定周波数よりも低周波数側の帯域と高周波数側の帯域とを分ける窓関数を生成するステップと、前記窓関数および前記光学伝達関数を用いて、前記入力画像のうち前記低周波数側の帯域に対しては前記画像回復処理を行い、前記高周波数側の帯域に対しては前記画像回復処理を制限する画像回復フィルタを生成するステップと、前記画像回復フィルタを用いて前記画像回復処理を行うステップとを有することを特徴とする画像処理方法。前記特定周波数は、前記指標値が前記所定値となる周波数のうち最も低い周波数であることを特徴とする請求項1に記載の画像処理方法。前記窓関数は、前記アジムス方向に応じて変化する前記特定周波数に対応する非回転対称な関数であることを特徴とする請求項1又は2に記載の画像処理方法。前記窓関数を生成するステップにおいて、前記光学伝達関数の絶対値を用いて前記窓関数を生成することを特徴とする請求項1から3のいずれか一項に記載の画像処理方法。前記特定周波数を求めるステップにおいて、前記入力画像上での位置ごとに前記特定周波数を求め、前記画像回復フィルタを生成するステップにおいて、前記位置ごとに前記画像回復フィルタの特性を変更することを特徴とする請求項1から4のいずれか一項に記載の画像処理方法。撮像装置により生成された入力画像に対して画像回復処理を行う画像処理装置であって、前記入力画像および該入力画像を生成した際の前記撮像装置における撮影条件に関する情報を得る入力部と、前記撮影条件に関する情報に応じた前記撮像装置の光学伝達関数を得て、該光学伝達関数を用いて得られる指標値が所定値となる特定周波数をアジムス方向ごとに求める演算部と、前記アジムス方向ごとに求められた前記特定周波数よりも低周波数側の帯域と高周波数側の帯域とを分ける窓関数を生成する窓関数生成部と、前記窓関数および前記光学伝達関数を用いて、前記入力画像のうち前記低周波数側の帯域に対しては前記画像回復処理を行い、前記高周波数側の帯域に対しては前記画像回復処理を制限する画像回復フィルタを生成するフィルタ生成部と、前記画像回復フィルタを用いて前記画像回復処理を行う処理部とを有することを特徴とする画像処理装置。撮影光学系により形成された被写体像を光電変換する撮像素子を含み、該撮像素子からの出力を用いて撮影画像を生成する撮像系と、前記撮影画像である入力画像に対して画像回復処理を行う請求項6に記載の画像処理装置とを有することを特徴とする撮像装置。コンピュータに、撮像装置により生成された入力画像に対する画像回復処理を行わせる画像処理プログラムであって、前記入力画像および該入力画像を生成した際の前記撮像装置における撮影条件に関する情報を得るステップと、前記撮影条件に関する情報に応じた前記撮像装置の光学伝達関数を得るステップと、前記光学伝達関数を用いて得られる指標値が所定値となる特定周波数をアジムス方向ごとに求めるステップと、前記アジムス方向ごとに求められた前記特定周波数よりも低周波数側の帯域と高周波数側の帯域とを分ける窓関数を生成するステップと、前記窓関数および前記光学伝達関数を用いて、前記入力画像のうち前記低周波数側の帯域に対しては前記画像回復処理を行い、前記高周波数側の帯域に対しては前記画像回復処理を制限する画像回復フィルタを生成するステップと、前記画像回復フィルタを用いて前記画像回復処理を行うステップとを有することを特徴とする画像処理プログラム。しかしながら、光学系のレンズの特性や撮影条件に依っては、撮像素子のナイキスト周波数の帯域内でMTFが零又はほぼ零になる以下、これを零落ちといい、零落ちが生じる周波数を零落ち周波数と称する場合がある。零落ちの原因としては、収差や回折が挙げられる。また、撮像装置の使用者の手振れによっても零落ちが発生する。零落ちの例を図12aに示す。横軸は空間周波数を、縦軸はMTFを示す。図中の下向き矢印の周波数で零落ちが生じており、この周波数が零落ち周波数である。画像回復処理に一般的に用いられるウィナーフィルタWienerfilterを用いた場合、画像回復フィルタの周波数特性の絶対値回復ゲイン特性は、図12bに示すようになる。図12bにおいて、横軸は空間周波数を、縦軸は回復ゲイン値を示している。ウィナーフィルタの詳細については後述する。この場合、図中の下向き矢印よりも低周波数側の帯域は増幅されているが、該矢印よりも高周波数側の帯域に対してはローパスフィルタ効果が得られることになる。このときの実空間での画像回復フィルタのプロファイル、すなわちフィルタ内の位置タップごとの係数値は図12cに示すようになる。この画像回復フィルタを用いた画像回復処理を行うと、入力画像における高周波数側のディテールがぼけてしまう。また、図13を用いて別の例を示す。図13aの横軸は、サンプリング周波数ナイキスト周波数の2倍までの周波数を示している。図13bは回復ゲイン特性を示し、図13cは実空間での画像回復フィルタのプロファイルを示す。図13cでは、画像回復フィルタの係数値が狭い位置範囲内の複数箇所で大きく変化している。ここで、光学系の製造誤差や入力画像中の輝度飽和等の原因で、入力画像のPSFが光学系のPSFと異なる場合には、回復された画像回復画像にリンギング等のアーティファクトが発生する場合がある。このような場合に、図13cに示す画像回復フィルタはPSFの不一致に対して敏感となり、この結果、アーティファクトが発生し易くなる。このように、高精細な回復画像を得るためには、零落ちに対する対策がきわめて重要である。特許文献1では、画像回復の度合いを調節可能にしているが、これは使用者による画像回復の度合いの変更を可能にするものであって、零落ち周波数に応じて画像回復フィルタの特性を変更するものではない。このため、特許文献1にて開示された方法では、零落ちによるアーティファクトを抑制するために画像回復の度合いを下げることはできるが、高精細な回復画像を得ることはできない。本発明は、MTFの零落ちが存在するような場合でも、高精細な回復画像を得ることができるようにした画像処理技術を提供する。本発明の画像処理方法は、撮像装置により生成された入力画像に対して画像回復処理を行う。該画像処理方法は、入力画像および該入力画像を生成した際の撮像装置における撮影条件に関する情報を得るステップと、撮影条件に関する情報に応じた撮像装置の光学伝達関数を得るステップと、光学伝達関数を用いて得られる指標値が所定値となる特定周波数をアジムス方向ごとに求めるステップと、アジムス方向ごとに求められた特定周波数よりも低周波数側の帯域と高周波数側の帯域とを分ける窓関数を生成するステップと、窓関数および光学伝達関数を用いて、入力画像のうち上記低周波数側の帯域に対しては画像回復処理を行い、上記高周波数側の帯域に対しては画像回復処理を制限する画像回復フィルタを生成するステップと、画像回復フィルタを用いて画像回復処理を行うステップとを有することを特徴とする。また、本発明の他の一側面としての画像処理装置は、撮像装置により生成された入力画像に対して画像回復処理を行う。該画像処理装置は、入力画像および該入力画像を生成した際の撮像装置における撮影条件に関する情報を得る入力部と、撮影条件に関する情報に応じた撮像装置の光学伝達関数を得て、該光学伝達関数を用いて得られる指標値が所定値となる特定周波数をアジムス方向ごとに求める演算部と、アジムス方向ごとに求められた特定周波数よりも低周波数側の帯域と高周波数側の帯域とを分ける窓関数を生成する窓関数生成部と、窓関数および光学伝達関数を用いて、入力画像のうち上記低周波数側の帯域に対しては画像回復処理を行い、上記高周波数側の帯域に対しては画像回復処理を制限する画像回復フィルタを生成するフィルタ生成部と、画像回復フィルタを用いて画像回復処理を行う処理部とを有することを特徴とする。なお、撮影光学系により形成された被写体像を光電変換する撮像素子を含み、該撮像素子からの出力を用いて撮影画像を生成する撮像系と、撮影画像である入力画像に対して画像回復処理を行う上記画像処理装置とを有する撮像装置も、本発明の他の一側面を構成する。さらに、本発明の画像処理プログラムは、コンピュータに、撮像装置により生成された入力画像に対する画像回復処理を行わせる、該画像処理プログラムは、入力画像および該入力画像を生成した際の前記撮像装置における撮影条件に関する情報を得るステップと、撮影条件に関する情報に応じた撮像装置の光学伝達関数を得るステップと、光学伝達関数を用いて得られる指標値が所定値となる特定周波数をアジムス方向ごとに求めるステップと、アジムス方向ごとに求められた特定周波数よりも低周波数側の帯域と高周波数側の帯域とを分ける窓関数を生成するステップと、窓関数および光学伝達関数を用いて、入力画像のうち上記低周波数側の帯域に対しては画像回復処理を行い、上記高周波数側の帯域に対しては画像回復処理を制限する画像回復フィルタを生成するステップと、画像回復フィルタを用いて画像回復処理を行うステップとを有することを特徴とする。本発明によれば、入力画像のうち特定周波数例えば、零落ち周波数よりも低周波数側の帯域に対しては画像回復処理を行い、該特定周波数よりも高周波数側の帯域に対しては画像回復処理を制限する画像回復フィルタを生成する。これにより、零落ちが存在するような場合でも、零落ち周波数より高周波数側にて画像のぼけが生じない高精細な回復画像を得ることができる。以下、本発明の実施例について図面を参照しながら説明する。まず、以下に説明する実施例にて用いる用語の定義と画像回復処理について説明する。入力画像入力画像は、撮像装置において撮影光学系により形成された被写体像を光電変換した撮像素子からの出力を用いて生成されたデジタル画像である。このデジタル画像は、レンズや光学フィルタ等の光学素子を含む撮影光学系の収差を含む光学伝達関数OTFにより劣化した画像である。撮像素子は、CMOSやCCD等の光電変換素子により構成される。撮影光学系は、曲率を有するミラー反射面を含んでもよい。また、撮影光学系は、撮像装置に対して着脱交換が可能であってもよい。撮像装置において、撮像素子および該撮像素子の出力を用いてデジタル画像入力画像を生成する信号処理回路により撮像系が構成される。入力画像の色成分は、例えばRGB色成分の情報を有している。色成分の扱いとしては、これ以外にもLCHで表現される明度、色相および彩度や、YCbCrで表現される輝度および色差信号等、一般に用いられている色空間を選択して用いることができる。その他の色空間としては、例えば、XYZ,Lab,Yuv,JChを用いることが可能であり、さらに色温度を用いることも可能である。入力画像や回復画像出力画像には、入力画像を生成した際の撮像装置における撮影光学系の焦点距離や絞り値、さらに撮影距離等の撮影条件に関する情報以下、撮影条件情報というを付帯することができる。また、入力画像を補正するための各種の補正情報も付帯することができる。撮像装置から、これとは別に設けられた画像処理装置に入力画像を出力し、該画像処理装置にて画像回復処理を行う場合には、入力画像に撮影条件情報や補正情報を付帯することが好ましい。撮影条件情報や補正情報は、入力画像に付帯する以外に、撮像装置から画像処理装置に直接または間接的に通信により受け渡すこともできる。画像回復処理撮像装置により生成された入力画像劣化画像をgx,yとし、元の画像劣化していない画像をfx,yとし、光学伝達関数OTFのフーリエペアである点像分布関数PSFをhx,yとする場合、以下の式が成り立つ。ただし、はコンボリューション畳み込み積分または積和を示し、x,yは入力画像上の座標位置を示す。gx,y=hx,yfx,yこの式をフーリエ変換して周波数面での表示形式に変換すると、以下の式のように周波数ごとの積の形式になる。Hは点像分布関数PSFhをフーリエ変換したものであり、光学伝達関数OTFに相当する。G,Fはそれぞれ、g,fをフーリエ変換したものである。u,vは2次元周波数面での座標、すなわち周波数を示す。Gu,v=Hu,vFu,v撮像装置により生成された劣化画像から元の画像を得るためには、以下のように、上記式の両辺をHで除算すればよい。Gu,v/Hu,v=Fu,vFu,v、すなわちGu,v/Hu,vを逆フーリエ変換して実面に戻すことで、元の画像fx,yである回復画像が得られる。ここで、H-1を逆フーリエ変換したものをRとすると、以下の式のように実面での画像に対するコンボリューション処理を行うことで、同様に元の画像fx,yである回復画像を得ることができる。gx,yRx,y=fx,yRx,yが画像回復フィルタである。入力画像が2次元であるとき、一般に画像回復フィルタも該2次元画像の各画素に対応したタップセルを有する2次元フィルタとなる。また、一般に画像回復フィルタのタップ数セル数が多いほど画像回復精度が向上するため、出力画像としての要求画質、画像処理装置としての画像処理能力、撮影光学系の収差の特性等に応じて実現可能なタップ数を設定する。画像回復フィルタは、少なくとも収差の特性を反映している必要があるため、従来の水平垂直各3タップ程度のエッジ強調フィルタハイパスフィルタ等とは全く異なる。また、画像回復フィルタは、光学伝達関数OTFに基づいて生成されるため、劣化画像入力画像における振幅成分と位相成分の劣化をともに高精度に補正することができる。また、実際の入力画像にはノイズ成分が含まれる。このため、上記のように光学伝達関数OTFの完全な逆数をとって作成した画像回復フィルタを用いると、劣化画像が回復されるだけでなくノイズ成分が大幅に増幅されてしまう。これは、入力画像の振幅成分にノイズの振幅が付加されている状態に対して撮影光学系のMTF振幅成分を全周波数にわたって1に戻すようにMTFを持ち上げるためである。撮影光学系による振幅劣化であるMTFは1に戻るが、同時にノイズ成分のパワースペクトルも持ち上がってしまい、結果的にMTFを持ち上げる度合い、すなわち回復ゲインに応じてノイズが増幅されてしまう。したがって、ノイズ成分がある入力画像からは鑑賞用画像として良好な回復画像が得られない。これを式で示すと以下のように表せる。Nはノイズ成分を表す。Gu,v=Hu,vFu,vNu,vGu,v/Hu,v=Fu,vNu,v/Hu,vこの点については、例えば、式1に示すウィナーフィルタのように画像信号とノイズ信号の強度比SNRに応じて回復度合いを制御する方法が知られている。Mu,vはウィナーフィルタの周波数特性を示し、|Hu,v|は光学伝達関数OTFの絶対値MTFを示す。この方法は、周波数ごとに、MTFが小さいほど回復ゲインを抑制し、MTFが大きいほど回復ゲインを強くするものである。一般に、撮影光学系のMTFは、低周波数側が高く、高周波数側が低くなるため、実質的に画像信号の高周波数側の回復ゲインを抑制する方法となる。画像回復フィルタを図2および図3を用いて説明する。画像回復フィルタは、撮影光学系の収差特性や要求される画像回復精度に応じてタップ数が決められる。図2では、例として、11×11タップの2次元画像回復フィルタを示している。図2では、各タップ内の値係数値を省略しているが、この画像回復フィルタの1つの断面を図3に示す。画像回復フィルタの各タップ内の値の分布が、収差によって空間的に広がった信号値PSFを、理想的には元の1点に戻す役割を果たす。画像回復処理では、画像回復フィルタの各タップの値が、入力画像における各タップに対応する各画素に対してコンボリューション畳み込み積分や積和ともいうされる。コンボリューションの処理では、ある画素の信号値を改善するために、その画素を画像回復フィルタの中心と一致させる。そして、入力画像と画像回復フィルタの対応画素ごとに入力画像の信号値と画像回復フィルタのタップの値係数値との積をとり、その総和を中心画素の信号値として置き換える。画像回復処理の実空間と周波数空間での特性を図4および図5を用いて説明する。図4のaは画像回復前のPSFを示し、bは画像回復後のPSFを示している。また、図5のMのaは画像回復前のMTFを示し、Mのbは画像回復後のMTFを示している。さらに、図5のPのaは画像回復前のPTFを示し、Pのbは画像回復後のPTFを示している。画像回復前のPSFは非対称な広がりを持っており、この非対称性によりPTFは周波数に対して非直線的な値を持つ。画像回復処理は、MTFを増幅し、PTFを零に補正するため、画像回復後のPSFは対称で、かつ鮮鋭になる。画像回復フィルタの作成法については、撮影光学系の光学伝達関数OTFの逆関数に基づいて設計した関数を逆フーリエ変換して得ることができる。例えば、ウィナーフィルタを用いる場合、式1を逆フーリエ変換することで、実際に入力画像に畳み込む実空間の画像回復フィルタを作成することができる。また、光学伝達関数OTFは、同じ撮影条件であっても撮影光学系の像高画像上での位置に応じて変化するので、画像回復フィルタは像高に応じて変更して使用される。以下、本発明の具体的な実施例について説明する。図1には、本発明の実施例1である画像処理方法の手順を示している。この画像処理方法は、画像処理装置に含まれるCPU等により構成されるコンピュータがコンピュータプログラムとしての画像処理プログラムに従って実行する。ステップS11では、コンピュータは、撮像装置から入力画像以下、撮影画像ともいうを取得する。撮像装置からの撮影画像の取得は、該撮像装置と画像処理装置とを有線または無線による通信を介して行ってもよいし、半導体メモリや光ディスク等の記憶媒体を介して行ってもよい。次にステップS12では、コンピュータは、撮像装置が撮影画像を生成した際の撮影条件に関する情報以下、撮影条件情報というを取得する。撮影条件は、撮像光学系の焦点距離や絞り値のほか、撮影距離や撮像装置の識別情報カメラID等を含む。また、撮影光学系の交換が可能な撮像装置においては、撮影条件に、該撮影光学系交換レンズの識別情報レンズIDを含む。撮影条件情報は、撮影画像に付帯された情報として取得してもよいし、有線または無線による通信や記憶媒体を介して取得してもよい。次にステップS13では、コントローラは、撮影条件に対応する適した光学伝達関数を取得する。この光学伝達関数は、画像処理装置内又は該装置外のメモリに予め保持された複数の光学伝達関数から撮影条件に対応するものを選択して取得してもよい。また、光学伝達関数を生成するための関数を予め記憶しておき、その変数に撮影条件を代入して生成算出することで取得してもよい。なお、離散的に選択された焦点距離、絞り値および撮影距離等の撮影条件に対する光学伝達関数を予めメモリに記憶しておいてもよい。この場合において、実際の撮影条件が記憶された撮影条件と異なるときは、記憶された撮像条件のうち実際の撮像条件に近い複数の撮像条件に対応する光学伝達関数からの補間演算により実際の撮影条件に対応する光学伝達関数を生成できる。これにより、予めメモリに記憶させる光学伝達関数のデータ量を削減することができる。補間処理の方法としては、バイリニア補間線形補間やバイキュービック補間等が挙げられるが、どのような方法であってもよい。次にステップS14では、コンピュータは、ステップS13で取得した光学伝達関数を用いて得られる指標値が所定値としての閾値となる特定周波数を求める。本実施例では、光学伝達関数OTFの絶対値MTFを指標値として用いるが、光学伝達関数を用いて得られる指標値であれば、他の指標値を用いてもよい。詳細については後述する。次にステップS15では、コンピュータは、特定周波数を零とする窓関数を生成する。窓関数の詳細については後述する。次にステップS16では、コンピュータは、光学伝達関数と窓関数とを用いて画像回復フィルタを生成する。この詳細についても後述する。次にステップS17では、コンピュータは、撮影画像に画像回復フィルタをコンボリューションして撮影画像に対する画像回復処理を行う。そして、ステップS18にて回復画像を出力する。ここで、ステップS14での特定周波数の求め方について詳しく説明する。図6には、MTFの周波数特性を示している。ここで示すMTFは、横軸の1を撮像素子の画素ピッチで決まるナイキスト周波数として正規化したものである。図6中のa,b,cは、撮影光学系において光学特性が互いに異なる3つのアジムス方向でのMTFの例を示している。a,bのMTFが閾値t1となる周波数fa,fbがそれぞれの特定周波数である。cのMTFは閾値以下とならないので、特定周波数は存在しない。閾値t1は、目的に応じて適宜設定することが可能である。例えば、特定周波数を、前述した零落ち周波数として使用する場合には、t1=0とすればよい。ただし、実用的にはMTFが完全に0にならなくとも十分小さな値として0と等価とみなせるので、t1=0.03MTFが3等と0よりも大きな値として、MTFがこの値となる周波数を零落ち周波数としての特定周波数としてもよい。偽解像がある場合等のように周波数に対してMTFが0の近傍で増加したり減少したりする場合には、MTFが閾値に対して高い側から低い側に該閾値を通過するように変化する周波数が複数存在することになる。この場合は、該複数の周波数のうち最も低い周波数を特定周波数として用いることが好ましい。さらに、図6に示すようにMTFがアジムス方向によって異なる場合には、特定周波数もアジムス方向によって異なるように設定することが望ましい。この場合、アジムス方向によって異なる特定周波数を関数データとして保持したり、ルックアップテーブルデータとして保持したりすることができる。また、アジムス方向によって異なる特定周波数のうち1つ例えば最も低い特定周波数を代表特定周波数として用いることもできる。次に、ステップS16での画像回復フィルタを生成する処理について詳しく説明する。ここでは、画像回復フィルタを設計するための基本となる関数として、例えば前述した式1により表されるウィナーフィルタを用いる。式1中のMu,vと窓関数Wu,vとを用いて、式2に従って画像回復フィルタの周波数特性MLu,vを生成する。MLu,vは、窓関数により画像回復処理を行う周波数帯域以下、画像回復帯域ともいうを制限された画像回復フィルタである。言い換えれば、窓関数により、撮影画像のうち特定周波数よりも低周波数側の帯域に対しては画像回復処理を行い、特定周波数よりも高周波数側の帯域に対しては画像回復処理を制限するように生成された画像回復フィルタである。このような画像回復フィルタを図9に示す。図9aには、図12aと同じMTFを示している。横軸は空間周波数を、縦軸はMTFを示す。また、図9bには、特定周波数零落ち周波数に対応して画像回復帯域を制限した画像回復フィルタの周波数特性を示している。横軸は空間周波数を、縦軸は回復ゲイン値を示している。図12bとの比較から分かるように、図9bの画像回復フィルタでは、高周波数側がローパスフィルタ化していないため、図9cに示す画像回復フィルタのプロファイル横軸はフィルタでのタップ位置を、縦軸はタップの係数値を示すが尖鋭である。このため、この画像回復フィルタを用いて画像回復処理を行うと、撮影画像においてMTFが低いもののディテールとして存在する画像部分をそのままにして画像回復処理によりぼかすことなく、低周波数側の画像部分を良好に回復させることができる。また、図10aには、図13aと同じMTFを示している。横軸は空間周波数を、縦軸はMTFを示す。さらに、図10bには、特定周波数零落ち周波数に対応して画像回復帯域を制限した画像回復フィルタの周波数特性を示している。横軸は空間周波数を、縦軸は回復ゲイン値を示している。また、図10cには、画像回復フィルタのプロファイルを示している。図13bとの比較において、図10bの画像回復フィルタでは、係数値が狭い位置範囲内の複数箇所で大きく変化するようなことがない。このため、この画像回復フィルタは、零落ちしていない周波数帯域の画像部分を良好に回復させながらも、PSFの不一致に対して耐性のある画像回復フィルタとなる。なお、撮像光学系のOTFの特性は、像高、すなわち撮影画像上での位置ごとに異なるため、零落ち周波数も像高に依存する。このため、撮影画像上での位置ごとに特定周波数を求め、該位置ごとに生成する画像回復フィルタの特性を変更するようにしてもよい。特定周波数が異なるという点以外は、本実施例で説明しているフローチャートに従って処理を行えばよい。次に、ステップS15にて生成する窓関数について詳しく説明する。図7には、互いに直交する2つのアジムス方向のMTFを示している。図8aおよびcには、図7のMTFの分布を上から二次元画像面に垂直な方向から見て示している。図中の破線は、図7に示した閾値t1により決定された、アジムス方向に応じて変化する特定周波数を示す。また、実線は、その内側においては値が0より大きい窓関数の値が、0となる周波数である。ここいう窓関数の値としての0は、1よりも十分小さく、0とみなせる値を含む意味である。図8bには、図8aに示した窓関数のプロファイルを示している。破線は図8aにおける垂直方向のプロファイルを、実線は同水平方向のプロファイルを示している。図8bに示すように、アジムス方向に応じたMTFの零落ち特性に応じて窓関数を生成することが好ましい。図8bの窓関数は、アジムス方向に応じて変化する特定周波数零落ち周波数に対応する非回転対称な関数である。別の例として、図8cおよびdには、アジムス方向に応じて変化する零落ち周波数特定周波数のうち最小値を特定周波数代表特定周波数として生成された、回転対称な窓関数を示している。この場合、MTFが零落ちしていないにもかかわらず、窓関数によって帯域を制限されているアジムス方向については画像回復帯域が狭くなる。なお、これらのように、アジムス方向に応じて変化する実際のMTFOTFの絶対値から得た特定周波数を用いて窓関数を生成することにより、撮影光学系の光学伝達関数OTFに精度良く適合した窓関数、ひいては画像回復フィルタを生成することができる。なお、窓関数としては、一般的に知られているガウス窓やハン窓等の窓関数を、光学伝達関数から得られたアジムス方向ごとの特定周波数特性に応じて非対称に変形して用いてもよいし、新たに窓関数を生成してもよい。さらに、アジムス方向に応じて変化する特定周波数に対して精度良く適合する窓関数を生成する方法として、各像高画像上での位置でのMTFに基づいて窓関数Wu,vを生成することができる。この例として、式3が挙げられる。ただし、Nは実数のパラメータである。画像回復処理では、画像回復フィルタを生成するために像高ごとのOTFを用いるため、このように像高ごとのMTFを用いた窓関数を生成することにより、新たにデータ量が増大することはない。さらに、特定周波数のアジムス特性を精度良く関数で再現するのが難しい場合にも、MTFは対応することができる。これは窓関数で抑制すべき周波数とMTFとを対応させることができるためである。図1のフローチャートに戻り、ステップS17では、コンピュータは、ステップS16で生成した画像回復フィルタを用いて撮影画像にコンボリューションを行い、さらに所定の現像処理を行うことで、出力画像としての回復画像を、ステップS18にて出力する。本実施例によれば、MTFの零落ちが存在するような入力画像に対しても、零落ち周波数より高周波数側での画像のぼけやリンギング等の弊害の発生を抑制しつつ、高精細な回復画像を得ることができる。次に、本発明の実施例2である撮像装置について、図11を用いて説明する。該撮像装置は、実施例1にて説明した画像処理方法による画像回復処理を行う画像処理装置を搭載している。撮影光学系201は、不図示の被写体からの光に被写体像を形成させる。撮像素子202は、該被写体像を光電変換して電気信号であるアナログ信号を出力する。撮像素子202から出力されたアナログ信号は、A/Dコンバータ203によってデジタル信号に変換され、画像処理部204に入力される。画像処置部204は、入力されたデジタル信号に対して各種画像処理を行うことにより撮影画像入力画像を生成する。そして、該撮影画像に対して実施例1で説明した画像回復処理を行う。撮像素子202から画像処理部204のうち撮影画像を生成する部分までの系が撮像系に相当する。また、画像処理部204のうち画像回復処理を行う部分が、入力部、演算部、窓関数生成部、フィルタ生成部および処理部を含む画像処理装置に相当する。画像処理部204に画像回復処理を行わせるために、状態検知部207は、撮影画像を生成する際における撮像装置の撮像条件情報を得る。撮像条件情報については、実施例1にて説明した通りである。状態検知部207は、撮像条件情報をシステムコントローラ210から得てもよい。また、撮影光学系に関する撮像条件情報については、撮影光学系201の絞り201aの動作やズームレンズ又はフォーカスレンズ等のレンズ201bの移動を制御する撮影光学系制御部206から得てもよい。そして、画像処理部204は、図1のフローチャートにより説明した画像回復処理を実行する。光学伝達関数OTFまたは光学伝達関数OTFを生成するための係数データは、予め記憶部208に保持されている。画像処理部204は、画像回復処理により生成した出力画像である回復画像を半導体メモリや光ディスク等の画像記録媒体209に出力して記録させたり、表示部205に出力して表示させたりする。以上説明した一連の動作は、システムコントローラ210により制御される。なお、撮影光学系201にはローパスフィルタや赤外線カットフィルタ等の光学素子を含ませてもよい。ただし、ローパスフィルタ等、撮影光学系201の光学伝達関数OTFに影響を与える光学素子を用いる場合には、該光学素子に関する画像回復フィルタを作成する時点で該光学素子に関する考慮が必要になる。赤外カットフィルタにおいても、分光波長の点像分布関数PSFの積分値であるRGBチャンネルのそれぞれのPSF、特にRチャンネルのPSFに影響するため、画像回復フィルタを作成する時点での考慮が必要になる。なお、前述したように、撮影光学系201は、撮像装置に対して交換可能なものであってもよい。実施例2では、画像処理方法を使用する画像処理装置を搭載した撮像装置について説明したが、本発明の画像処理方法は、パーソナルコンピュータにインストールされる画像処理プログラムによっても実施することができる。この場合、パーソナルコンピュータが本発明の画像処理装置に相当する。パーソナルコンピュータは、撮像装置により生成された画像回復処理前の撮影画像入力画像を、該撮像装置から有線/無線通信を介して、または他のパーソナルコンピュータからインターネット等の回線を介して取得する。撮影画像が記録された画像記録媒体を介して該撮影画像を取得してもよい。そして、撮影画像を取得したパーソナルコンピュータは、画像処理プログラムによって画像回復処理を行い、その結果得られた回復画像を出力する。以上説明した各実施例は代表的な例にすぎず、本発明の実施に際しては、各実施例に対して種々の変形や変更が可能である。"