"光学系を用いた撮像により得られた画像を処理する画像処理方法であって、画像を取得するステップと、前記光学系の収差情報に基づいて、フィルタ値が2次元分布を持つように作成された画像回復フィルタを用いて前記画像に対する画像回復処理を行う画像回復ステップとを有し、該画像回復ステップにおいて、前記画像回復フィルタとして、少なくとも2つの色において互いに異なる前記フィルタ値の2次元分布を有する複数の画像回復フィルタを用い、前記複数の画像回復フィルタは、前記フィルタ値の絶対値が最大であるセルの位置と該画像回復フィルタの中心セルの位置とのずれ量が前記光学系の倍率色収差量に応じて互いに異なり、前記画像のぼけ成分とともに倍率色収差成分を低減するフィルタであることを特徴とする画像処理方法。前記画像回復フィルタは、光が前記光学系に入射してから前記撮像により入力画像が取得されるまでの光学伝達関数の逆関数に基づいて生成された関数を逆フーリエ変換することにより作成されることを特徴とする請求項1に記載の画像処理方法。前記画像回復ステップにおいて、前記撮像に際しての撮像状態を示す情報を取得し、該撮像状態情報に応じた前記画像回復フィルタを選択又は作成して用いることを特徴とする請求項1又は2に記載の画像処理方法。光学系を用いた撮像により得られた画像を処理する画像処理装置であって、前記光学系の収差情報に基づいて、フィルタ値が2次元分布を持つように作成された画像回復フィルタを記憶する記憶手段と、前記画像回復フィルタを用いて前記画像に対する画像回復処理を行う画像回復手段とを有し、該画像回復手段は、前記画像回復フィルタとして、少なくとも2つの色において互いに異なる前記フィルタ値の2次元分布を有する複数の画像回復フィルタを用い、前記複数の画像回復フィルタは、前記フィルタ値の絶対値が最大であるセルの位置と該画像回復フィルタの中心セルの位置とのずれ量が互いに異なり、前記画像のぼけ成分とともに倍率色収差成分を低減するフィルタであることを特徴とする画像処理装置。光学系により形成された被写体像を光電変換して画像を取得する撮像系と、該画像を処理する請求項4に記載の画像処理装置とを有することを特徴とする撮像装置。光学系を用いた撮像により得られた画像を処理する画像処理装置に以下のステップを実行させるプログラムであって、画像を取得するステップと、前記光学系の収差情報に基づいて、フィルタ値が2次元分布を持つように作成された画像回復フィルタを用いて前記画像に対する画像回復処理を行う画像回復ステップとを有し、該画像回復ステップにおいて、前記画像回復フィルタとして、少なくとも2つの色において互いに異なる前記フィルタ値の2次元分布を有する複数の画像回復フィルタを用い、前記複数の画像回復フィルタは、前記フィルタ値の絶対値が最大であるセルの位置と該画像回復フィルタの中心セルの位置とのずれ量が前記光学系の倍率色収差量に応じて互いに異なり、前記画像のぼけ成分とともに倍率色収差成分を低減するフィルタであることを特徴とするプログラム。光学系の諸収差により劣化した画像を良好に補正して高画質な画像を得るためには、上述したぼけ成分と色ずれ成分を低減する処理を行う必要がある。しかしながら、これらの処理を行うソフトウェアやそれを実装した撮像装置では、ユーザーにストレスを感じさせないために処理速度を高速にしなければならない。また、特許文献1にて開示された画像処理方法では、光学系の色による撮像倍率の相違により発生する倍率色収差をR,G,Bの色成分ごとに拡大縮小を行うことで、色ずれ成分のある程度の補正が可能とも考えられる。しかしながら、実際の画像における色ずれ成分の発生原因を光学的に考察すると、この方法では、以下のように十分な補正効果が得られない。光学系は、被写体から入射した光をレンズにより屈折させて像面に結像させる。この際、色波長によってレンズの屈折率が異なるため、色ごとの焦点距離に差が生じる。光学系は、一般に、波長分散と形状を考慮した複数のレンズを用いて色収差を抑制するように設計されるが、光学系の小型化やコスト低減の要求により、色収差がある程度残留する。一方、被写体像を光電変換する撮像素子を含む撮像系では、被写体像をカラーフィルタを通してR,G,B信号に分解し、これらR,G,B信号の混合によってカラー画像を生成する。したがって、R,G,Bの各チャンネルの信号値は、撮像面に到達した400nm~700nmの波長の光をカラーフィルタの分光透過率特性に応じて積分した強度情報となっている。このため、G信号の中でも、G用カラーフィルタの分光透過率特性に応じた透過波長帯域幅の間で焦点距離の異なった波長光が積分されている。このことは、R信号及びB信号でも同様である。すなわち、倍率色収差は、R,G,Bチャンネルごとに画像上の結像位置がずれるとともに、各チャンネルの像をぼかしてしまう。したがって、特許文献1の画像処理方法によってR,G,B信号の位置ずれを幾何変換により補正しても十分な色ずれ成分の補正はできない。また、特許文献1の画像処理方法では、色ずれ成分以外のぼけ成分色にじみ成分を含むの補正を行うことはできない。特許文献1にて開示された色ずれ成分の補正方法と、従来の画像回復フィルタを用いた画像回復処理とを順次行うことで、色ずれ成分だけでなくぼけ成分を低減することは可能である。ただし、処理工程が増え、処理速度が低下する。本発明は、画像のぼけ成分と倍率色収差成分を高速に補正することができるようにした画像処理方法、画像処理装置、撮像装置、及びプログラムを提供する。本発明の一側面としての画像処理方法は、光学系を用いた撮像により得られた画像を処理する画像処理方法であって、画像を取得するステップと、光学系の収差情報に基づいて、フィルタ値が2次元分布を持つように作成された画像回復フィルタを用いて該画像に対する画像回復処理を行う画像回復ステップとを有する。画像回復ステップにおいて、画像回復フィルタとして、少なくとも2つの色において互いに異なるフィルタ値の2次元分布を有する複数の画像回復フィルタを用いる。そして、複数の画像回復フィルタは、フィルタ値の絶対値が最大であるセルの位置と該画像回復フィルタの中心セルの位置とのずれ量が光学系の倍率色収差量に応じて互いに異なり、画像のぼけ成分とともに倍率色収差成分を低減するフィルタであることを特徴とする。また、本発明の他の一側面としての画像処理装置は、光学系を用いた撮像により得られた画像を処理する画像処理装置であって、光学系の収差情報に基づいて、フィルタ値が2次元分布を持つように作成された画像回復フィルタを記憶する記憶手段と、画像回復フィルタを用いて該画像に対する画像回復処理を行う画像回復手段とを有する。画像回復手段は、画像回復フィルタとして、少なくとも2つの色において互いに異なるフィルタ値の2次元分布を有する複数の画像回復フィルタを用いる。そして、複数の画像回復フィルタは、フィルタ値の絶対値が最大であるセルの位置と該画像回復フィルタの中心セルの位置とのずれ量が互いに異なり、画像のぼけ成分とともに倍率色収差成分を低減するフィルタであることを特徴とする。また、本発明の他の一側面としてのプログラムは、光学系を用いた撮像により得られた画像を処理する画像処理装置に以下のステップを実行させるプログラムであって、画像を取得するステップと、光学系の収差情報に基づいて、フィルタ値が2次元分布を持つように作成された画像回復フィルタを用いて該画像に対する画像回復処理を行う画像回復ステップとを有する。画像回復ステップにおいて、画像回復フィルタとして、少なくとも2つの色において互いに異なるフィルタ値の2次元分布を有する複数の画像回復フィルタを用いる。そして、複数の画像回復フィルタは、フィルタ値の絶対値が最大であるセルの位置と該画像回復フィルタの中心セルの位置とのずれ量が光学系の倍率色収差量に応じて互いに異なり、画像のぼけ成分とともに倍率色収差成分を低減するフィルタであることを特徴とする。なお、光学系により形成された被写体像を光電変換して画像を取得する撮像系と、該画像を処理する上記画像処理装置とを有する撮像装置も本発明の他の一側面を構成する。本発明によれば、画像回復フィルタがぼけ成分とともに倍率色収差成分を低減するように作成されているので、画像回復処理によってぼけ成分と倍率色収差成分とを同時に低減することができる。したがって、ぼけ成分と倍率色収差成分とを高速な処理によって低減することができる。図1には、本発明の実施例である画像処理装置を備えた画像処理方法を使用するデジタルカメラやビデオカメラ等の撮像装置の構成を示している。不図示の被写体からの光束は、撮像光学系101によって、CCDセンサやCMOSセンサ等により構成される撮像素子102上に結像する。撮像光学系101は、不図示の変倍レンズ、絞り101a及びフォーカスレンズ101bを含む。変倍レンズを光軸方向に移動させることで、撮像光学系101の焦点距離を変更するズームが可能である。また、絞り101aは、絞り開口径を増減させて、撮像素子102に到達する光量を調節する。フォーカスレンズ101bは、被写体距離に応じてピント調整を行うために、不図示のオートフォーカスAF機構やマニュアルフォーカス機構によって光軸方向の位置が制御される。撮像素子102上に形成された被写体像は、該撮像素子102により電気信号に変換される。撮像素子102からのアナログ出力信号は、A/Dコンバータ103によりデジタル撮像信号に変換され、画像処理部104に入力される。画像処理部104は、入力されたデジタル撮像信号に対して各種処理を行うことで、カラー入力画像を生成する画像生成部104aを有する。また、画像処理部104は、この入力画像に対して画像回復処理を行う画像回復処理部画像回復手段104bを有する。撮像素子102から画像生成部104aまでが撮像系に相当する。画像回復処理部104bは、状態検知部107から撮像光学系101の状態以下、撮像状態というの情報を得る。撮像状態とは、例えば、撮像光学系101の焦点距離ズーム位置、絞り開口径絞り値、Fナンバー、フォーカスレンズ位置被写体距離である。なお、状態検知部107は、システムコントローラ110から撮像状態の情報を得てもよいし、撮像光学系101を制御する撮像光学系制御部106から得てもよい。そして、画像回復処理部104bは、撮像状態に応じた画像回復フィルタを記憶部記憶手段108から選択し、入力画像に対して画像回復処理を行う。状態検知部107、画像回復処理部104b及び記憶部108により、撮像装置内で画像処理装置が構成される。図2には、画像回復処理部104b以下の説明では、画像処理部104というで行われる画像回復に関する処理画像処理方法のフローチャートを示している。画像処理部104は、画像処理用コンピュータにより構成され、コンピュータプログラムに従って該処理を実行する。ステップS1では、画像処理部104は、撮像素子102からの出力信号に基づいて生成された画像以下、入力画像というを画像回復処理の対象として取得する。次に、ステップS2では、画像処理部104は、状態検知部107から撮像状態情報を取得する。ここでは、撮像状態を、ズーム位置、絞り開口径、及び被写体距離の3つとする。次に、ステップS3では、画像処理部104は、記憶部108に記憶された画像回復フィルタの中から、ステップS2で取得した撮像状態に対応する画像回復フィルタを選択する。ここで、記憶部108には、画像回復フィルタの数データ数を低減するために、離散的に選択された撮像状態に対する画像回復フィルタのみが記憶格納されている。このため、ステップS2で取得した撮像状態に対応する又は該撮像状態にきわめて近い撮像状態に対応する画像回復フィルタが記憶部108に格納されていない場合は、その撮像状態にできるだけ近い画像回復フィルタを選択する。そして、その画像回復フィルタを、以下のステップS4~S6で、ステップS2で取得した撮像状態に最適化するように補正することで、実際に使用する画像回復フィルタを作成する。図3には、記憶部108に記憶された、離散的に選択された撮像状態に対する画像回復フィルタを模式的に示す。記憶部108に記憶された画像回復フィルタは、前述したように、ズーム位置状態A、絞り開口径状態B及び被写体距離状態Cの3つの撮像状態を軸とした撮像状態空間中に離散的に配置されている。撮像状態空間中の各点黒丸の座標が、記憶部108に記憶されている画像回復フィルタを示す。なお、図3では、画像回復フィルタを各撮像状態に対して直交した線上の格子点に配置しているが、画像回復フィルタを格子点から外して配置しても構わない。また、撮像状態の種類は、ズーム位置、絞り開口径及び被写体距離に限らず、その数も3つでなくてもよく、4つ以上の撮像状態による4次元以上の撮像状態空間を構成してその中に画像回復フィルタを離散的に配置してもよい。このように離散的に配置された画像回復フィルタから、実際に使用する画像回復フィルタを作成する方法については後述する。ここで、画像回復フィルタの例を図4に示す。画像回復フィルタでは、撮像光学系101の収差量に応じてセルタップ数、すなわちカーネルサイズが決められる。図4に示す画像回復フィルタは、11×11セルを有する2次元フィルタである。各セルが画像の1画素に対応する。画像回復フィルタを、100以上のセルに分割した2次元フィルタとすることで、撮像光学系101による球面収差、コマ収差、軸上色収差、軸外色フレア等の結像位置から大きく広がる収差に対しても良好な画像回復結果を得ることができる。また、このような実空間での画像回復フィルタを、入力画像に対してコンボリューション処理することで、フーリエ変換を行うことなく画像を回復することができる。また、各セルは、撮像光学系101の球面収差、コマ収差、軸上色収差、軸外色フレア等の収差情報に応じて、図5に示すような値を有するように設定される。図5には、画像回復フィルタの1つの断面でのセル値フィルタ値を示している。画像回復フィルタは、フィルタ値が2次元分布を持つように作成されている。また、図5の例では、画像の中心から外れた位置での、コマ収差のような非対称性をもった収差を回復するフィルタとなっている。画像回復フィルタは、撮像光学系101等の光学伝達関数OTFを計算若しくは計測し、その逆関数を逆フーリエ変換することで作成される。一般的にはノイズの影響を考慮する必要があるため、ウィーナーフィルタや関連する回復フィルタの作成方法を選択して用いることができる。また、ここにいう光学伝達関数は、撮像光学系101のみならず、撮像素子102の出力信号から画像処理部104にて入力画像が生成されるまでの、光学伝達関数を劣化させる要因を含めることが望ましい。すなわち、画像回復フィルタを、光が撮像光学系101に入射してから撮像素子102による撮像によって入力画像が取得されるまでの光学伝達関数の逆関数に基づいて生成された関数を逆フーリエ変換することにより作成するようにするとよい。撮像光学系101以外で光学伝達関数を劣化させる要因としては以下のようなものがある。例えば、撮像素子102の前面に配置されるローパスフィルタ図示せずは、光学伝達関数の周波数特性に対して高周波成分を抑制する。また、同様に撮像素子102の前面に配置される赤外線カットフィルタは、分光波長の点像分布関数PSFの積分値であるRGBチャンネルの各PSF、特にRチャンネルのPSFに影響する。さらに、撮像素子102の画素開口の形状や開口率も周波数特性に影響する。他にも、被写体を照らす光源の分光特性や、各種波長フィルタの分光特性も光学伝達関数を劣化させる要因として挙げられる。したがって、これらを考慮した広義の光学伝達関数に基づいて画像回復フィルタを作成することが望ましい。なお、撮像光学系101は、撮像装置の一部として設けられてもよいし、撮像装置に対して着脱可能な交換式のものであってもよい。また、本実施例では、入力画像がRGB形式のカラー画像であるため、R,G,Bの複数の色成分に対応した複数ここでは3つの画像回復フィルタを作成する。撮像光学系101には色収差があり、色成分ごとにぼけ方が異なるため、色成分ごとに最適な画像回復フィルタを得るためには、その特性を色収差に基づいて異ならせるべきだからである。つまり、図5に示したセル値の2次元分布が、色成分ごとに異なる3つの画像回復フィルタを作成する。なお、図4には、11×11セルを有する画像回復フィルタを示したが、画像回復フィルタのセル数は、収差量に応じて任意に設定すればよい。例えば、収差量が大きい場合には、30×30セルや50×50セルのようにセル数を増加させてもよい。また、画像回復フィルタの縦横のセル数セル配列に関しては、図4に示すような正方配列である必要はなく、コンボリューション処理時にセル配列を考慮するようにすれば、任意に変更することができる。画像回復フィルタの具体的な選択及び作成補正方法について説明する。図3において、大きな白丸で示した撮像状態が、ステップS2で取得した実際の撮像状態であるとする。実際の撮像状態又はそのきわめて近くに、記憶部108に記憶格納された画像回復フィルタ以下、格納フィルタというが存在する場合には、その格納フィルタを選択して画像回復処理に用いる。実際の撮像状態又はそのきわめて近くに、格納フィルタが存在しない場合には、以下の方法により画像回復フィルタを選択又は作成補正する。まず、前述したステップS3において、画像処理部104は、実際の撮像状態と複数の格納フィルタに対応する撮像状態との間の撮像状態空間内での距離をそれぞれ算出する。そして、算出した距離のうち最も短い距離にある撮像状態に対応する格納フィルタを選択する。このような格納フィルタを選択することで、実際の撮像状態と該格納フィルタに対応する撮像状態との相違量以下、状態相違量というが最も小さくなる。このため、該格納フィルタに対する補正量を少なくすることができ、実際の撮像状態に対応する画像回復フィルタにより近い画像回復フィルタを作成することができる。図3において、小さな白丸で示した撮像状態に対応する格納フィルタが選択されたものとする。なお、図3には、撮像状態空間を概念的に示しているが、実際の各格納フィルタのデータには、座標値の情報が必要である。このため、格納フィルタ自体に座標値情報を持たせたり、予めアドレス座標が決まった多次元の配列空間に各格納フィルタのデータを配置したりすればよい。次に、ステップS4では、画像処理部104は、ステップS3で選択した格納フィルタに対応する撮像状態と実際の撮像状態との状態相違量ΔA,ΔB,ΔCを算出する。そして、ステップS5では、これらの状態相違量ΔA,ΔB,ΔCに基づいて、状態補正係数を算出する。さらに、画像処理部104は、ステップS6において、該状態補正係数を用いて、ステップS3で選択した格納フィルタを補正する。これにより、実際の撮像状態に対応した画像回復フィルタを作成することができる。次に、ステップS7では、画像処理部104は、選択又は作成された画像回復フィルタを用いて、ステップS1で取得した入力画像に対して画像回復処理を行う。すなわち、上記画像回復フィルタを用いて入力画像に対するコンボリューション処理を行うことで、撮像光学系101の収差に起因する画像のぼけ成分及び色にじみ成分を低減又は除去することができる。なお、ステップS3~ステップS7での処理は、撮像装置とは別の装置パーソナルコンピュータ等で行ってもよい。この場合、該別の装置で得られた画像回復フィルタを、撮像装置の記憶部108に格納インストールする処理を行えばよい。図1において、画像処理部104で画像回復処理が行われた出力画像は、半導体メモリや光ディスク等の画像記録媒体109に所定のフォーマットで保存される。この出力画像は、画像回復処理によって鮮鋭化された高画質の画像である。また、出力画像は、表示部105に表示されたり、撮像装置の外部プリンタやデータベースに出力されたりする。上述した撮像素子102での光電変換、画像処理部104での画像処理、画像記録媒体109への記録、及び表示部105への画像表示といった一連の動作の制御は、システムコントローラ110によって行われる。また、撮像光学系101のズーム駆動やフォーカス駆動は、システムコントローラ110の指示により撮像光学系制御部106が制御する。なお、撮像光学系101の中心像高の横倍率が最大像高の横倍率よりも1.25倍以上大きくなるように撮像光学系101を設計するとよい。これにより、撮像光学系101のうち最も被写体側に位置する前玉レンズの径を小さくしたり、像面湾曲の補正を行ったりすることができる。本実施例では、画像回復処理によって画像のぼけ成分と色にじみ成分だけでなく、色ずれ成分も低減補正する。このため、画像回復処理において、ぼけ成分と色にじみ成分と色ずれ成分を補正できるよう、R,G,Bの色成分ごとの複数3つの画像回復フィルタを、前述したように撮像状態に応じて選択したり選択した画像回復フィルタを補正して作成したりする。色成分ごとの画像回復フィルタの例を図6に示す。図6においてR,G,Bで示す画像回復フィルタは、それぞれR,G,B用の画像回復フィルタである。これらの画像回復フィルタは、ぼけ成分と色にじみ成分を補正するフィルタであるが、倍率色収差に対応する色ずれ成分の補正機能は有していない。また、R′,G′,B′で示す画像回復フィルタも、それぞれR,G,B用の画像回復フィルタである。ただし、これらの画像回復フィルタは、ぼけ成分と色にじみ成分の補正機能に加えて、色ずれ成分の補正を行う機能も有する。具体的には、該3つの画像回復フィルタのうち少なくとも1つここではR′とB′において、セル値フィルタ値の絶対値が最大であるセル以下、最大値セルというの位置が、その画像回復フィルタの中心セルの位置からずれている。各画像回復フィルタ内で黒いセルが最大値セルであり、×を付したセルが中心セルである。R,G,Bで示す画像回復フィルタは、色成分ごとに光学伝達関数OTFが異なることからフィルタ特性は異なるが、最大値セルの位置は互いに同じ中心セルの位置である。一方、R′,G′,B′で示す画像回復フィルタは、フィルタ特性が異なっているのに加えて、最大値セルの位置の中心セルに対するずれ量が互いに異なっている。これら画像回復フィルタ間での最大値セルの位置の相対ずれ量は、撮像光学系101の倍率色収差量に応じた量である。図7には、R′で示す画像回復フィルタの1つの断面でのセル値フィルタ値を示している。R′,G′,B′で示す画像回復フィルタにおいて、フィルタ内のセル値の分布を色成分ごとに異ならせることで色ごとのぼけ成分及び色にじみ成分を補正できる。さらに、R′,G′,B′で示す画像回復フィルタにおいて、上記分布の位置をずらす最大値セルの位置と中心セルの位置とのずれ量を異ならせる。ただし、本実施例ではG′でのずれ量は0ことで色ずれ成分も補正できる。なお、図6では、11×11セルの画像回復フィルタを示したが、収差量が大きい場合には、30×30セルや50×50セルのように、収差量に応じてセル数を増加させてもよい。また、画像回復フィルタ内のセル値の分布を移動させたときに、元のセル領域から分布がはみ出してしまう場合には、画像回復フィルタのセル数を増加させることで対応すればよい。次に、R′,G′,B′で示す画像回復フィルタの効果を、図8を用いて説明する。図8において、aは撮像光学系101の収差の影響がなければ矩形のプロファイルを持つはずの被写体像が撮像光学系101の収差によってぼけ像として劣化している入力画像を示している。実線、破線及び一点鎖線はそれぞれ、G成分、R成分及びB成分を表しており、aでは、G成分、R成分及びB成分でぼけの広がりが互いに異なる上に、ピーク位置が互いにずれている。各色成分でのプロファイルの劣化がぼけ成分であり、色成分間でのぼけ量の違いによる色づきが色にじみ成分であり、ピーク位置のずれが色ずれ成分である。bは、図7にR,G,Bで示す画像回復フィルタを用いて、aに示す入力画像に対する画像回復処理を行った場合の回復画像を示している。この回復画像において、R,G,Bの被写体像が本来の矩形プロファイルを持つようにぼけ成分色にじみ成分を含むは補正されたが、色ずれ成分は補正されていない。さらに、色成分ごとに画像の倍率を拡大縮小して色ずれを補正する方法により、cに示すような画像を得ることは可能である。一方、R′,G′,B′で示す画像回復フィルタを用いてaに示す入力画像に対して画像回復処理を行うと、c′に示すように、色成分ごとのぼけ成分色にじみ成分を含むが補正されるとともに、色ずれ成分も同時に補正される。このように、本実施例によれば、従来では互いに別の処理として行う必要があったぼけ成分の補正処理画像回復処理と色ずれ成分の補正処理とを、1つ一度の画像回復処理の中で行うことができる。したがって、ぼけ成分と色ずれ成分倍率色収差成分とを高速な処理によって低減することができ、良好な出力画像を得ることができる。前述した光学伝達関数OTFは、同じ撮像状態であっても画角像高に応じて変化する。このため、画像回復処理は、画像における像高に応じて分割された領域ごとに画像回復フィルタを変更して行うことが望ましい。画像回復フィルタを画像上で走査しながらコンボリューション処理を行い、分割領域ごとに画像回復フィルタを順次変更すればよい。また、画像回復処理の対象とする入力画像における画素の信号値に応じて画像回復フィルタを補正することで、入力画像の暗い領域で顕著となるノイズの増加を抑制することができる。入力画像の暗い領域では、明るい領域に比べて、画像信号に対してノイズ成分が相対的に大きいため、画像回復処理によってノイズを強調してしまうおそれがある。このため、暗い領域に対しては、画像回復フィルタの各セルの値を補正して画像回復効果を低下させることで、ノイズの増大を抑制することができる。信号値がある閾値以下の画素又は領域については画像回復処理を行わないようにしてもよい。また、画像回復処理後の回復画像に補正誤差がある場合には、この補正誤差を補正するための別の画像処理を行ってもよい。これにより、実際の撮像状態と画像回復処理に用いた画像回復フィルタに対応する撮像状態との誤差に起因する回復画像における残留収差成分をさらに低減することができる。また、入力画像内に飽和輝度領域が存在する場合に、その周辺部では収差成分が顕著に現れる。しかし、画像回復処理を行う際に、輝度値が飽和している画素については本来の輝度値が分からず、良好な補正が行えないおそれがある。この場合には、画像回復処理後の回復画像から残留している収差成分を検出して適応的に補正処理を行うことで、さらなる高画質化が達成できる。この適応的な補正処理は、本実施例での画像回復処理によってぼけ成分が低減された画像に最終的に残留した収差を補正するので、画像回復処理よりも後の段階で行うことが望ましい。なお、本実施例では、画像がRGB形式の場合について説明したが、YCC形式等、他の色形式の画像にも本発明を適用することができる。また、本実施例では、本発明の画像処理方法を使用する画像処理装置を搭載した撮像装置について説明したが、本発明の画像処理方法は、パーソナルコンピュータにインストールされる画像処理プログラムによっても実施することができる。この場合、パーソナルコンピュータが本発明の画像処理装置に相当する。パーソナルコンピュータは、撮像装置により生成された画像回復処理前の画像入力画像を取り込み取得し、画像処理プログラムによって画像回復処理を行って、その結果得られた画像を出力する。以上説明した各実施例は代表的な例にすぎず、本発明の実施に際しては、各実施例に対して種々の変形や変更が可能である。"