"撮像装置によって撮像された画像における、前記撮像装置の撮像光学系の収差による画質の劣化を回復するための画像処理装置であって、前記画像に対して前記撮像光学系の倍率色収差の補正を行う第1の倍率色収差補正手段と、前記倍率色収差が補正された画像について、色に関する特徴量を取得する第1の取得手段と、前記撮像光学系の収差を示す関数に基づいた画像回復フィルタを前記倍率色収差が補正された画像に適用する画像回復処理手段と、前記画像回復フィルタが適用された画像に対し、前記色に関する特徴量を取得する第2の取得手段と、前記第1の取得手段と前記第2の取得手段が取得した前記色に関する特徴量の変化に応じて、前記画像回復フィルタが適用された画像の画素値を補正する画素値補正手段とを有し、前記画像回復フィルタは、倍率色収差の補正を行わないフィルタであることを特徴とする画像処理装置。前記画像がRGB成分から構成され、前記色に関する特徴量が、G成分の値に対するR成分もしくはB成分の値の差分であり、前記画素値補正手段は画素のR成分もしくはB成分の値を補正することを特徴とする請求項1に記載の画像処理装置。前記画像がRGB成分から構成され、前記色に関する特徴量が、RGB成分から生成した輝度成分の値に対するR成分もしくはB成分の値の差分であり、前記画素値補正手段は画素のR成分もしくはB成分の値を補正することを特徴とする請求項1に記載の画像処理装置。さらに、前記画像回復フィルタが適用された画像から倍率色収差を検出して補正する第2の倍率色収差補正手段を有し、前記第2の取得手段が、前記第2の倍率色収差補正手段によって倍率色収差が補正された画像に対して前記色に関する特徴量を取得することを特徴とする請求項1乃至請求項3のいずれか1項に記載の画像処理装置。前記画像回復フィルタは、前記撮像光学系の点像分布関数の色成分間の差異を、ある色成分の点像分布関数の位置に他の色成分の点像分布関数を平行移動させることによって低減した点像分布関数に基づいて生成されたものであることを特徴とする請求項1乃至請求項4のいずれか1項に記載の画像処理装置。前記画像回復フィルタは、前記撮像光学系の光学伝達関数の位相の周波数特性から直線成分を除去した前記光学伝達関数に基づいて生成されたものであることを特徴とする請求項1乃至請求項4のいずれか1項に記載の画像処理装置。被写体の光学像を結像する撮像光学系と、前記光学像を撮像する撮像素子と、請求項1乃至請求項6のいずれか1項に記載の画像処理装置とを有することを特徴とする撮像装置。撮像装置によって撮像された画像における、前記撮像装置の撮像光学系の収差による画質の劣化を回復するための、画像処理装置によって実行される画像処理方法であって、前記画像処理装置の倍率色収差補正手段が、前記画像に対して前記撮像光学系の倍率色収差の補正を行う工程と、前記画像処理装置の第1の取得手段が、前記倍率色収差が補正された画像について、色に関する特徴量を取得する第1の取得工程と、前記画像処理装置の画像回復処理手段が、前記撮像光学系の収差を示す関数に基づいた画像回復フィルタを前記倍率色収差が補正された画像に適用する工程と、前記画像処理装置の第2の取得手段が、前記画像回復フィルタが適用された画像に対し、前記色に関する特徴量を取得する第2の取得工程と、前記画像処理装置の前記第1の取得工程と前記第2の取得工程で取得された前記色に関する特徴量の変化に応じて、画素値補正手段が、前記画像回復フィルタが適用された画像の画素値を補正する工程とを有し、前記画像回復フィルタは、倍率色収差の補正を行わないフィルタであることを特徴とする画像処理方法。コンピュータを、請求項1乃至請求項6のいずれか1項に記載の画像処理装置の各手段として機能させるためのプログラム。このように、撮像画像に対して撮像光学系の点像分布関数を用いた画像回復処理を施すことにより、諸収差を補正して画質を向上させる技術が知られている。しかしながら、実際の撮像では、入力画像の撮像状態と、適用する画像回復フィルタの状態が最適には一致しない場合がある。例えば、立体被写体の撮像画像がその一例である。撮像装置は、オートフォーカス機能やマニュアルフォーカスにより被写体空間の1つの面に焦点を合わせて撮像するため、合焦面に位置する被写体は比較的先鋭に撮像される。しかし、他の被写体同一被写体の、合焦面から異なる位置にある部分を含むは合焦面からの距離に応じたぼけ量で撮像される。そして、被写体距離に関する情報が、合焦面までの距離のみの場合、この被写体距離と画角に最適な画像回復フィルタを選択若しくは生成して使用することになる。その結果、合焦物体に対しては最適な画像回復フィルタが適用されるため、望ましい回復結果を得ることができるが、非合焦物体については画像回復フィルタが最適でないため、多少の回復効果はあるものの、ぼけを解消することはできない。一方で、写真において被写体のぼけは、被写体の立体感や注視物体を表現するための手法として用いられている。例えば、被写界深度の浅い望遠レンズを用いて主被写体にピントを合わせ、背景を意図的にぼかす写真表現などである。このような写真表現の存在を考慮すれば、合焦被写体はより先鋭化され、非合焦被写体にはぼけが残存する上述の画像回復処理は適切とも考えられる。しかし、合焦距離の被写体に最適であり、非合焦距離の被写体には最適でない画像回復フィルタを用いて画像回復処理した場合、非合焦距離の被写体に色付きが発生する場合がある。ここでの色付きとは、画像の各色成分のぼけ方の関係が画像回復処理の前後で異なることにより、画像回復処理後の画像の例えば非合焦距離の被写体非合焦被写体のエッジ部に、被写体に無い色偽色が発生することである。さらに、このような色付きは立体被写体の撮像に限らず発生する場合がある。本質的には合焦しているか否かに関わらず、画像を撮像した際の収差の状態と、適用する画像回復フィルタが補正する収差の状態とが異なる場合に色付きが発生する。このような色付きを低減する方法として、画像回復処理前の画像の色情報に基づいて、画像回復処理後の画像の色を補正する方法がある。つまり画像の各画素において画像回復処理による色の変化を判定し、画像回復処理による色付きを低減させるという方法である。しかし、画像回復処理前の撮像画像にはすでに、軸上色収差、色の球面収差、色のコマ収差が原因となる色のにじみや、倍率色収差が原因となる色ズレなど、撮像光学系の諸収差による色付きが存在している。画像回復処理前の画像にこのような色付きが存在する場合、画像回復処理前の画像の色情報に基づく方法では画像回復処理後の画像における色付きを十分に抑制できないことがある。画質を向上させるための画像回復処理において生じる色付きは、無視できない画質劣化であり、十分な抑制が必要である。特許文献1では、撮像画像の非合焦範囲に対して画像回復処理を行うことで、撮像光学系の浅い被写界深度を補おうとするものである。特許文献1記載の画像回復処理では、非業小範囲における先鋭度の向上は実現できたとしても、画像回復処理前の撮像画像に色付きがある場合、画像回復処理後の画像に発生する色付きを十分抑制することはできない。本発明はこのような従来技術の課題に鑑みてなされたものであり、画像回復処理後の画像における色付きを抑制することの可能な画像処理装置および画像処理方法を提供することを目的とする。上述の目的は、撮像装置によって撮像された画像における、撮像装置の撮像光学系の収差による画質の劣化を回復するための画像処理装置であって、画像に対して撮像光学系の倍率色収差の補正を行う第1の倍率色収差補正手段と、倍率色収差が補正された画像について、色に関する特徴量を取得する第1の取得手段と、撮像光学系の収差を示す関数に基づいた画像回復フィルタを倍率色収差が補正された画像に適用する画像回復処理手段と、画像回復フィルタが適用された画像に対し、色に関する特徴量を取得する第2の取得手段と、第1の取得手段と第2の取得手段が取得した色に関する特徴量の変化に応じて、画像回復フィルタが適用された画像の画素値を補正する画素値補正手段とを有し、画像回復フィルタは、倍率色収差の補正を行わないフィルタであることを特徴とする画像処理装置によって達成される。このような構成により本発明によれば、画像回復処理後の画像における色付きを抑制することの可能な画像処理装置および画像処理方法を提供することができる。以下、添付図面を参照して、本発明をその例示的な実施形態に基づいて詳細に説明する。はじめに、本発明の実施形態における画像回復処理の概要と、倍率色収差の検出及び補正処理について説明する。画像回復処理基本的には、上述した従来の画像回復処理と同様、画像回復フィルタを用いる。以下の式は、画像回復フィルタとして使用可能な一例としてのウィナーフィルタの特性を示す。ここで、Mu,vはウィナーフィルタの周波数特性、|Hu,v|は光学伝達関数の絶対値MTF、SNRは画像信号とノイズ信号との強度比である。上式からわかるように、この方法は周波数ごとに、MTFが小さいほど回復ゲイン回復度合を抑制し、MTFが大きいほど回復ゲインを強くするものである。一般的に撮像光学系のMTFは低周波側が大きく高周波側が小さくなるため、実質的に画像の高周波側の回復ゲインを抑制している。上述の通り、画像回復フィルタは2次元画像に適用する場合、画素に対応するタップを有する2次元フィルタとして構成され、タップ数は撮像光学系の収差特性や要求される回復精度に応じて決定することができる。図2aは、例として、水平及び垂直方向に11タップを有する11×11タップの2次元フィルタで構成された画像回復フィルタを模式的に示している。図2aでは各タップの値係数値の記載を省略しているが、画像回復フィルタの1水平ラインのタップ係数の分布例を図2bに示す。画像回復フィルタのタップの値係数値の分布が、撮像光学系の収差によって空間的な広がりを持つ点像分布関数を、理想的には元の1点に戻す回折の影響を無視した場合役割を果たしている。画像回復フィルタのタップ係数値は、画像回復処理の工程で実行される画素ごとのコンボリューション処理畳み込み積分、積和に用いられる。コンボリューション処理では、処理対象の画素の信号値に画像回復処理を適用するために、処理対象の画素を画像回復フィルタの中心タップと一致させる。そして、画像回復フィルタの各タップ位置において、対応する画素信号値とタップ係数値との積を求め、積の総和で処理対象の画素の信号値を置き換える。画像回復処理の実空間と周波数空間での特性を図3を用いて説明する。図3において、aは画像回復処理前の点像分布関数を、bは回復後の点像分布関数をそれぞれ示している。また、図3cの401は画像回復処理前のMTF、402は画像回復処理後のMTFを示し、図3dの403は画像回復処理前のPTF、404は画像回復処理後のPTFを示している。図3aからわかるように、画像回復処理前の点像分布関数は非対称な広がりを有し、この非対称性により、対応するPTFは零ではない値を有する図3d403。画像回復処理は、図3cの402に示すようにMTFを増幅し、図3dの404に示すようにPTFを零に補正するため、画像回復処理後の点像分布関数は図3bに示すように対称で先鋭になる。画像回復フィルタは、撮像光学系の光学伝達関数の逆関数1/Hに基づいて設計した関数を逆フーリエ変換して得ることができる。画像回復フィルタには上述のウィナーフィルタを用いることができるが、他のフィルタを用いてもよい。ウィナーフィルタを用いる場合、光学伝達関数の逆関数を逆フーリエ変換することで、コンボリューション処理に用いる実空間の画像回復フィルタを作成することができる。次に、倍率色収差を考慮した画像回復フィルタの生成について説明する。図4aの502は、501の結像位置を基準として平行移動した位置に同形状の点像分布関数として結像した2つの色成分を表している。このときのPTFを図4bに示す。図4bの601、602は図4aの501、502と対応しており、602は位相が傾きをもった直線となっている。画像回復フィルタは、このような直線的なPTFも図3dの403のような非直線のPTFも、もととなる光学伝達関数に含まれているため図4bの601のように補正することができる。また、予めPTFから直線成分を除去した光学伝達関数をもとに画像回復フィルタを生成することで、平行移動としての倍率色収差の成分を補正しない画像回復フィルタを生成することができる。一方、実空間で倍率色収差成分を除去する場合、図4aの502を501の位置に平行移動して差異を低減した点像分布関数を生成し、フーリエ変換することで、同様に倍率色収差の成分を除去した光学伝達関数を生成することができる。実際には501と502の点像分布関数の形状は異なるため、例えば重心を一致させる方法や501と502の差分の二乗平均を最小化するなどして両者を位置合わせすることができる。この光学伝達関数をもとに画像回復フィルタを生成することで、平行移動としての倍率色収差の成分を補正しない画像回復フィルタを生成することができる。倍率色収差の成分を補正しない画像回復フィルタを適用することで、倍率色収差を画像回復処理と独立して処理することが可能になる。つまり、色成分ごとの非対称収差の補正および先鋭化と、倍率色収差の平行移動の成分である色ズレの補正とを分離して処理することが可能になる。なお、光学伝達関数は同一の撮像状態ズーム位置や絞りの開口径であっても撮像光学系の像高画像の位置に応じて変化するので、画像回復フィルタは像高に応じて変更して使用することが望ましい。以上が、画像回復処理の概要である。倍率色収差の検出と補正倍率色収差の補正量は、上述のようにPTFの直線成分や点像分布関数の形状など、撮像光学系の設計値から算出することができる。しかし、撮像光学系の製造ばらつきや撮像時の光源分光の変動により倍率色収差の画像上での現れ方が変化する場合があり、このような変化に対応する方法として画像から倍率色収差を検出する方法がある。画像から倍率色収差を検出して補正する方法に関して図5を用いて説明する。後述するように、本実施形態に係る撮像装置図1においては、画像処理部104の有する画像回復処理部111が以下の処理を実行することができるが、ここでは一般的な手順として説明する。まず、S201で撮像画像を入力画像として取得する。次に、S202では、倍率色収差による色ズレが顕著に現れるエッジ部分を画像から検出する。エッジの検出にはY輝度プレーンを使用する。光学中心からの動径方向に大きく画素値の変化するエッジに限定して検出することにより、精度の高い色ズレ量取得が可能になる。また、Yプレーンにおいて、倍率色収差による色ズレはにじみとなって現れるので、画素値の単調増加もしくは単調減少が複数画素連続するような、ある程度幅のあるエッジを検出するのが良い。次に、S203では、S202で検出された各エッジにおいて色ズレ量を取得する。ここで扱う色ズレの方向は、光学中心と各エッジの位置関係により上/下方向、左/右方向、右斜め上/左斜め下方向、左斜め上/右斜め下方向の何れかを適用することで処理の簡略化を図る。各エッジにおける色ズレ量の取得には、色成分間の相関を用いる。例えば、色成分間の差分絶対値の和を判定して色ズレ量を取得することができる。Gプレーンに対しRプレーンもしくはBプレーンを色ズレ方向に移動させながら、検出されたエッジ付近の画素において、色成分間の差分絶対値和が最小となる場所を探索する。ここで検出された差分絶対値和が最小となった位置と、移動前の位置とから、Gプレーンに対するRプレーンもしくはBプレーンの色ズレ量を取得することができる。色ズレ量は、Gプレーンに対しRプレーンもしくはBプレーンが光学中心方向へずれている場合に負の値とし、Gプレーンに対しRプレーンもしくはBプレーンが光学中心とは逆方向へずれている場合に正の値とする。S204では、S202で検出された各エッジの像高と、S203で取得した各エッジの色ズレ量から、像高と色ズレの関係を求めることで補正データを作成する。ここでいう像高とは、光学中心に相当する画素以降、単に光学中心というからの距離である。以下、補正データの作成手順について具体的に説明する。1S202で検出されたエッジの像高をL、S203で取得された色ズレ量をDとした場合、下記のように像高に対する色ズレ率Mを求める。M=D/L式92図6aに示すように画像を像高に応じた複数の領域ここではh1~h8に分割し、エッジが属する領域を選別する。3画像内で検出された各エッジについて12を行い、像高に応じて分割した領域ごとに色ズレ率Mを集計し、各領域での色ズレ率Mの平均値を各領域の色ズレ率として求める。4図6bに示すように、像高と色ズレ率とから、像高lと色ズレ率Mとの関係を表す高次の多項近似式Flを算出し、これを補正データとする。図6bは以下の式10のように3次の多項式により補正データを算出した例である。なお式10においてa、b、cは係数を示している。M=Fl=al3bl2cl式10なお、エッジの検出と色ズレ量取得は、画像内の全エッジに対して行っても良い。しかし、像高に応じて分割した各領域で、ある閾値以上の数の色ズレ率が集計された段階でエッジ検出と色ズレ量取得を終了するなど、信頼度を保ちながら処理の効率化を図ることもできる。また、像高に応じて分割した複数の領域のうち、該当するエッジが見つかった領域のみを高次の多項近似式算出に使用することで、該当するエッジが見つからない領域があった場合でも補正データの作成が可能となる。次に、S205では、S204で作成された補正データを使用して色ズレを補正する。補正は対象色成分の画素位置をシフトすることで行い、具体的には例えば以下に示す方法で補正することができる。まず、補正したいプレーンRプレーン、Bプレーンの画素X,Yにおいて、画素X,Yの像高Lと、補正データとから、画素X,Yにおける色ズレ率Mを求める。なお、光学中心に対応する画素が0,0の座標系とする。M=FL式11次に色ズレ補正により生成する画素の座標X1,Y1を求める。X1=M×X式12Y1=M×Y式13補正したいプレーンにおいて上記座標X1,Y1に相当する画素値を、一般的な補間処理により生成し、画素X,Yの画素値とする。これらを全画素について行うことで、色ズレ補正を行う。以上が、画像から倍率色収差補正データを作成し、補正するまでの処理概要である。撮像装置の構成図1は、本発明の実施形態にかかる画像処理装置の一例としての撮像装置の構成例を示すブロック図である。本実施形態では撮像光学系および撮像素子を有する撮像装置に本発明を適用した例を説明するが、本発明において撮像光学系や撮像素子といった撮像画像を生成するための構成要素は必須でない。不図示の被写体の光学像を撮像光学系101で撮像素子102に結像する。撮像素子102で結像光が電気信号に変換され、A/Dコンバータ103でデジタル信号に変換されて、画像処理部104に入力される。画像処理部104は画像回復処理部111と、画像回復処理以外の画像処理を行う、他の画像処理部112から構成される。画像回復処理部111は、まず、状態検知部107から撮像装置の撮像状態の情報を得る。状態検知部107は撮像装置の撮像状態の情報をシステムコントローラ110から直接得ても良いし、例えば撮像光学系101に関する撮像状態の情報は撮像光学系制御部106から得てもよい。次に画像回復処理部111は、撮像状態に応じた画像回復フィルタを例えば記憶部108から選択し、画像回復処理部111は画像処理部104に入力された画像に対して倍率色収差補正処理および画像回復フィルタ適用処理を行う。画像回復処理の詳細については後述する。記憶部108で保持するデータは画像回復フィルタではなく、画像回復フィルタを生成するために必要な情報例えば点像分布関数や光学伝達関数に関する情報でもよい。この場合、画像回復処理部111は撮像状態に応じた点像分布関数又は光学伝達関数に関する情報を記憶部108から選択して撮像状態に応じた画像回復フィルタを生成し、この画像回復フィルタを用いて画像回復処理を行う。また、他の画像処理部112は画像回復処理後の画像に対し、ガンマ補正やカラーバランス調整など、所定の画像処理を行ないJPEG等の画像ファイルを生成する。画像処理部104で処理された出力画像は、システムコントローラ110が画像記録媒体109に所定のフォーマットで保存する。また、表示部105には、画像回復処理後の画像に表示用の所定の処理を行った画像を表示しても良いし、画像回復処理を行わない、又は簡易的な回復処理を行った画像を表示しても良い。一連の制御はシステムコントローラ110で行われ、撮像光学系101の機械的な駆動絞り101a、フォーカスレンズ101b、光学ズームの駆動などはシステムコントローラ110の指示により撮像光学系制御部106で行う。システムコントローラ110は例えばCPUやMPUのようなプログラマブルなコントローラであり、記憶部108に記憶されたプログラムを実行することで撮像装置の全体的な動作を実現する。なお、画像処理部104はハードウェアで実現してもよいが、少なくとも一部をシステムコントローラ110がソフトウェア的に実現することも可能である。絞り101aは、Fナンバーの撮像状態設定として開口径が制御される。フォーカスレンズ101bは、被写体距離に応じてピント調整を行うために不図示のオートフォーカスAF機構や手動のマニュアルフォーカス機構によりレンズの位置が制御される。撮像光学系101にはローパスフィルタや赤外線カットフィルタ等の光学素子を入れても構わない。ただし、ローパスフィルタ等の光学伝達関数の特性に影響を与える素子を用いる場合には、画像回復フィルタを作成する時点でその光学素子による光学伝達関数の変化を考慮する。赤外カットフィルタにおいても、分光波長の点像分布関数の積分値であるRGBチャンネルの各点像分布関数、特にRチャンネルの点像分布関数に影響するため、画像回復フィルタを作成する時点で赤外カットフィルタによる点像分布関数の変化を考慮する。また、撮像光学系101は撮像装置の一部として図示しているが、レンズ交換式カメラのような交換可能な構成であっても良い。図7に画像回復処理部111での画像回復処理を説明するためのフローチャートを示す。S301で画像回復処理部111は、状態検知部107から実際の撮像状態の情報を取得する。撮像状態とは、例えばズーム位置、絞り径、被写体距離等の撮像条件である。なお、これら撮像状態の情報は、画像ファイルのヘッダ部分などに書き込まれていることが多く、撮像済みの画像に対して画像回復処理を行う場合には、画像ファイルから取得してもよい。S302で第1の倍率色収差補正手段としての画像回復処理部111は、倍率色収差の補正を行う。ここで行う倍率色収差の補正は上述の通り画素位置をシフトすることによって行うものであり、光学系の諸収差を補正する画像回復フィルタは適用しない。倍率色収差の補正量は、撮像光学系の設計値から算出したものを使用してもよいし、上述したように画像から倍率色収差を検出して補正量を決定してもよい。S303で画像回復処理部111は、撮像条件に適した画像回復フィルタを選択する。このとき、選択された画像回復フィルタを必要に応じて補正しても構わない。これは、予め記憶部108に用意しておく画像回復フィルタのデータ数を低減するために、離散的な撮像状態のデータを用意しておいて、画像回復処理を実行する際に、実際の撮像状態に応じて画像回復フィルタを補正することである。また、画像回復フィルタの選択ではなく、画像回復フィルタを生成するために必要な光学伝達関数に関する情報から、撮像状態に応じた画像回復フィルタを生成してもよい。また、本実施形態で用いる画像回復フィルタは、倍率色収差補正のための平行移動成分が除去されている。つまり、画像回復フィルタは倍率色収差補正成分を含んでおらず、収差の非対称性および広がりを補正して画像の先鋭化を行うが、倍率色収差の平行移動成分に関しては補正しない。上述の通り、このような画像回復フィルタを用いることにより、倍率色収差を画像回復処理と独立して処理することが可能になる。つまり、色成分ごとの非対称収差の補正および先鋭化と、倍率色収差の平行移動の成分である色ズレの補正とを分離して処理することが可能になる。S304で第1の取得手段としての画像回復処理部111は、S302によって倍率色収差が補正された画像の色に関する特徴量の一例として、各画素における色差情報を取得する。ここでは、入力画像はRGB成分から構成され、色差を基準色Gに対する信号値の差分と定義する。従って、色差は以下の式により計算される。ここで、R1、G1、B1は、画像回復処理前のRGBの各信号値、x,yは画像上の座標値、C1r、C1bは、以下の式で定義されるR成分、B成分のG成分に対する色差である。C1rx,y=R1x,y-G1x,y式14C1bx,y=B1x,y-G1x,y式15なお、ここで、入力画像が例えばベイヤ配列の原色カラーフィルタを持つ撮像素子によって撮像された所謂RAW画像である場合、各画素はR,G,B成分のうち1つの値しか有していない。そのため、R画素又はB画素におけるG成分の値をG画素の値を適応補間して生成することにより、C1r及びC1bを求める。一方、既に色補間処理済みで、各画素がRGB成分の値を有している場合には、各画素においてC1r及びC1bを求める。また、各画素がRGBの成分の値を有している場合には、各画素においてRGBの成分の値からY輝度の成分の値を求め、R成分、B成分の輝度成分に対する色差をC1r及びC1bとして求めるようにしてもよい。次にS305で画像回復処理部111は、S303で選択もしくは生成した画像回復フィルタを用いて、入力画像の各画素に対して上述したコンボリューション処理を行う。これにより、撮像光学系101の収差による点像分布関数の非対称性を補正し、画像のぼけ成分を除去もしくは低減することができる。次のS306で第2の取得手段としての画像回復処理部111は、S305での画像回復処理後の画像の色に関する特徴量の一例として、各画素についてS304と同様に色差情報を取得する。ここで、R2、G2、B2は、画像回復処理後のRGBの各信号値、x,yは画像上の座標値、C2r、C2bは、以下の式で定義されるR成分、B成分のG成分に対する色差である。C2rx,y=R2x,y-G2x,y式16C2bx,y=B2x,y-G2x,y式17なお、ここでも、C1r及びC1bを求める場合と同様、入力画像が色補間前のRAW画像である場合、R画素又はB画素におけるG成分の値をG画素の値を適応補間して生成することにより、C2r及びC2bを求める。一方、既に色補間処理済みで、各画素がRGB成分の値を有している場合には、各画素においてC2r及びC2bを求める。次にS307で画素値補正手段としての画像回復処理部111は、S304とS306で取得した、画像回復処理前後の各画素の色差情報C1r、C1b、C2r、C2bに基づき、色付き抑制のための画素値調整が必要か否かを決定する。この決定は、以下に示すように、色差が増大しているか否か、色差の符号が反転しているか否かに基づいて行うことができる。画像回復処理部111は、RGBの画素のうちR画素とB画素又は、各画素におけるR成分とB成分に対してそれぞれ以下の判定を行い、いずれかを満たす場合はその色成分の画素に対し、色付き抑制のための画素値調整が必要と決定する。Rに対する判定判定条件1:C1rx,yとC2rx,yが同符号、且つ|C1rx,y|<|C2rx,y|判定条件2:C1rx,yとC2rx,yが異符号Bに対する判定判定条件1:C1bx,yとC2bx,yが同符号、且つ|C1rx,y|<|C2rx,y|判定条件2:C1bx,yとC2bx,yが異符号この判定結果により、画素値調整が必要と決定された画素について、画素値補正手段としての画像回復処理部111はS308において色付き抑制のための画素値調整を適用する。画素値調整は例えば以下に示すように、色差が増大した場合は画像回復前の色差を使用し、色差の符号が反転した場合は色差を0とするものであってよい。Rに対する画素値調整判定条件1を満たす場合色差増大;R2x,y=G2x,yC1rx,y式18判定条件2を満たす場合色差反転:R2x,y=G2x,y式19Bに対する画素値調整判定条件1を満たす場合色差増大:B2x,y=G2x,yC1rx,y式20判定条件2を満たす場合色差反転:B2x,y=G2x,y式21このように、図7のフローチャートに従った処理により、画像処理部104画像回復処理部111における画像回復処理が実行される。なお、撮像光学系101の光学伝達関数は1つの撮像状態においても撮像光学系の画角像高に応じて変化するので、ここで説明した画像回復処理を、像高に応じたて分割した領域図6のh1~h8ごとに変更して行うことが望ましい。具体的には、画像回復フィルタを画像上にコンボリューション処理をしながら走査させ、領域毎にフィルタを順次変更すれば良い。すなわちS303の処理を各対象画素に対して実行する。なお、本実施形態では理解及び説明を容易にするため、画像回復フィルタの適用=画像回復処理として説明したが、画像回復処理の一部として、他の補正処理、例えば歪曲補正処理や周辺光量補正処理やノイズ低減処理等を含めてもよい。ここで、S302、S304、S305、S306、S308で行なわれる倍率色収差補正、画像回復処理前後の色差情報取得、画像回復フィルタ適用の作用および効果の例を図9及び図10を用いて説明する。図9a~cは、図8に示したような、周辺部に点光源が存在する被写界を撮像した入力における、点光源とその近傍の水平方向における画素値の例を示す模式図である。図9a~cはいずれも上段にR成分とG成分の画素値を、下段には上段で示したR成分とG成分の差分、つまり色差を示している。図9bは、非合焦時の状態を示している。合焦時の状態を示す図9aと比べると、倍率色収差による色ズレは同様に発生しているが、G成分、R成分ともにぼけ量が大きく、またG成分とR成分は同じぼけ方をしている。倍率色収差補正後はR成分とG成分の位置が揃い、G成分とR成分のぼけ方が同じであるため、この段階で色にじみは発生していない。画像回復フィルタ適用後は、G成分、R成分ともにぼけが低減されているが、G成分とR成分のぼけの残り方が異なるため、色にじみが発生している。この色にじみが本実施形態において抑制の対象となる色付きである。図9bに対応する画像回復フィルタ適用前後の色差情報の関係を図10bに示す。画像回復フィルタの適用前である倍率色収差補正後の色差は全領域で0となっているが、画像回復フィルタ適用後はエッジ部分で色差が発生している。従って、S307における判定処理では色差増大により画素値調整が必要と決定され、S308で画像回復処理部111によって色付き抑制のための画素値調整が適用される。S308での画素値調整処理では、画像回復フィルタの適用前の色差である倍率色収差補正後の色差が使用されるため、画素値調整後の色差は全領域で0となり、色付きが適切に抑制された画像が得られる。このように、本実施形態では、非合焦時部分に対する画像回復フィルタの適用によって発生する色付きを効果的に抑制することができる。図9aで示した合焦時と図9bで示した非合焦時のそれぞれで、画像回復フィルタ適用前後のG成分とR成分のぼけ方の違いについて図11を用いて説明する。図11a~dは撮像光学系101のMTFと画像回復フィルタ適用後のMTFの空間周波数特性を示す模式図である。また、図11eは画像回復フィルタによるMTFの増減率回復ゲインを示す模式図である。図11aは合焦時の撮像光学系101のMTFであり、R成分はG成分よりも劣化している。従って図9aの入力画像ではG成分とR成分のぼけ方が異なり、点光源部分の画素値の変化はR成分の方がG成分よりも広がった状態になっている。図11bは合焦時の画像回復フィルタ適用後のMTFであり、G成分とR成分ともにMTFが適切に補正され、G成分とR成分は一致している。これは図11aのMTFに図11eに示す回復ゲインが乗算された結果である。画像回復フィルタは合焦物体に最適な画像回復処理が行なえるように設計されているため、合焦時は適切な画像回復処理がされる。従って図9aの画像回復フィルタ適用後ではG成分とR成分のぼけが適切に除去され、色にじみも除去されている。一方、図11cは非合焦時の撮像光学系101のMTFであり、図11aの合焦時と比較するとG成分もR成分も劣化した状態で一致している。従って図9bの入力画像ではG成分とR成分のぼけ方が一致しているが図9aで示した合焦時よりぼけが大きくなっている状態である。図11dは非合焦時の画像回復フィルタ適用後のMTFであり、G成分とR成分ともにMTFが補正されているが、G成分はR成分よりも劣化した状態になっている。これは図11cのMTFに図11eに示す回復ゲインが乗算された結果である。画像回復フィルタは合焦物体に最適な画像回復処理が行なえるように設計されているため、非合焦時は期待した画像回復結果にならない場合がある。従って、図9bの画像回復フィルタ適用後ではG成分とR成分のぼけが低減しているが、G成分の方がR成分よりも広がった状態になっている。本実施形態と異なり、画像回復フィルタに倍率色収差補正成分を含めた場合の例を図9cに示す。図9cは、図9bと同様の非合焦時において、倍率色収差補正成分を含めた画像回復フィルタを適用した場合の例を示している。入力画像は図9bと同じであり、また倍率色収差補正を行わないが、画像回復フィルタで倍率色収差補正も同時に行うため、画像回復フィルタ適用後の画像も図9bと同じである。画像回復フィルタ適用前後の色差情報の関係を図10cに示す。画像回復フィルタの適用前の入力画像の色差は色ズレにより左側のエッジ部分で正の値、右側のエッジ部分で負の値として存在する。また、画像回復フィルタ適用後は両エッジ部分で負の値の色差が発生している。従って、S307における判定処理では、左側のエッジ部では色差反転しているとの判定により画素値調整が必要と決定され、S308で色付き抑制のための画素値調整が行なわれるので色付きが抑制される。しかし、右側のエッジは色差の増大も反転もないことから、画素値調整は不要と決定され、S308での画素値調整は行われない。そのため、右側のエッジでは色付きが抑制されないまま画像回復処理結果として出力される。このように、画像回復フィルタに倍率色収差補正成分を含めた場合には画像回復フィルタの適用によって生じた色付きが十分に抑制されず、画質が劣化する。以上説明したように、本実施形態によれば、撮像光学系の収差を補正するための画像回復フィルタとして、倍率色収差補正を行わない画像回復フィルタを用いるとともに、倍率色収差は画像回復フィルタを適用する前に、画像回復フィルタとは別の方法で補正する。そのため、実際の撮像状態での収差の状態と、画像回復フィルタが想定している収差の状態が異なるときに発生する、画像回復フィルタに起因する色付きを、倍率色収差に起因する色付きと分離して補正することが可能になり、色付きを効果的に抑制可能となる。つまり、倍率色収差の色ズレによる色付きと、ぼけの除去時に発生する色付きとを切り分けることができ、ぼけの除去時に発生する色付きのみを抑制することが可能になる。このように、本実施形態によれば、画像回復処理全体として色付きを抑えた画像回復処理が可能となる。第2の実施形態次に、本発明の第2の実施形態について説明する。なお、本実施形態も第1の実施形態対と同様、図1に示した撮像装置に適用した場合について説明する。第1の実施形態とは画像処理部104画像回復処理部111における画像回復処理のみ異なるため、図12を用いて本実施形態における画像回復処理について説明する。なお、図12におけるS401からS405までの処理は、図7のS301からS305までの処理と同様であるため、説明は割愛する。本実施形態では、S405で画像回復フィルタを適用した後に、S406で画像から倍率色収差を検出して生成した補正量を使用した倍率色収差の補正を行うことを特徴とする。図13及び図14は、本実施形態においてS402、S404、S405、S406、S407、S408で行なわれる倍率色収差補正、画像回復処理前後の色差情報取得、画像回復フィルタ適用の作用及び効果の例を示す模式図である。図13a~dは、図8に示したような、周辺部に点光源が存在する被写界を撮像した入力における、点光源とその近傍のメリジオナル方向画像中心から放射方向における画素値の例を示す模式図である。図13a~dはいずれも上段にR成分とG成分の画素値を、下段には上段で示したR成分とG成分の差分、つまり色差を示している。図9とは着目している画素の方向が異なるため、非対称収差となり、左右のエッジでぼけ方が異なっている。倍率色収差を補正しない画像回復フィルタを適用した後に、S406で画像から倍率色収差を検出して補正することの効果について図13aを用いて説明する。図13aは、合焦時の状態を示している。入力画像には倍率色収差によりG成分とR成分の色ズレが発生し、さらにG成分とR成分のぼけ方が異なっており、左右のエッジの収差が異なるためにぼけ方が異なっている。画像回復フィルタ適用前に行うS402での倍率色収差補正の結果を、図13aの、倍率色収差補正後1で示している。S402の倍率色収差補正では、左右のエッジの収差の状態が異なったまま、色ズレができるだけ目立たないような倍率色収差補正を行う。倍率色収差の補正量は、撮像光学系101の設計値から算出したものを使用してもよいし、画像から倍率色収差を検出して補正量を決定してもよい。画像回復フィルタの適用により、G成分、R成分のエッジは非対称収差が補正され且つ先鋭化されているが、G成分、R成分間の色ズレが残留している。この、残留している色ズレの原因は、撮像光学系101の製造上のばらつきや、光源分光の変動などにより、画像撮像時の収差の状態と、画像回復フィルタが対象としている収差の状態とが異なることが挙げられる。また、画像回復フィルタ適用前に行うS402での倍率色収差補正において、画像から倍率色収差を検出して補正量を決定する場合、非対称性の残留したエッジで倍率色収差量の検出を行っている。そのため、画像回復処理後のエッジの状態に適切な倍率色収差の補正量と異なる補正量が得られる場合があり、その場合には色ズレが残留する。本実施形態はこのような理由によって画像回復フィルタの適用後に残留する色ズレに対処する。S406において第2の倍率色収差補正手段としての画像回復処理部111は、画像回復フィルタ適用後の画像から倍率色収差を検出して補正を行う。S406での倍率色収差補正処理後の状態を図13aの倍率色収差補正後2で示している。倍率色収差補正後2では残留した色ズレが適切に補正された状態であり、入力画像に対してはエッジは非対称収差が補正され且つ先鋭化されG成分とR成分が一致することで色にじみも除去されている。次に行なわれるS407からS409の処理は、第1の実施形態で説明した図7のS306からS308の処理と同様であるため、説明は割愛する。図13aに示した合焦時に、S404、S407で取得された色差情報を図14aに示す。この場合、画像回復フィルタの適用前の倍率色収差補正後1での色差はエッジ部分で値を持っているが、画像回復フィルタの適用後の倍率色収差補正後2では全領域で0となっている。従って、S408における判定処理では、色付き抑制のための画素値調整は不要と判断され、倍率色収差補正後2の画像がそのまま画像回復処理結果となる。図13aと同じ入力画像に対し、本実施形態と異なり画像回復フィルタ適用後にS406の倍率色収差補正を行なわない場合の例を図13bに示す。この場合、画像回復フィルタ適用後の色ズレが残留している状態でS407の色差情報取得処理が行われる。図14bに示すように、画像回復処理前後で左側のエッジでは色差の反転、右側のエッジでは色差の増大がそれぞれ発生しているため、いずれもS409で色付き抑制のための画素値調整が行なわれる。この場合の出力画像は図13bの画素値調整後の状態である。S409における色付き抑制処理によって色付きは抑制されているが、S406の倍率色収差補正を行なった場合の図12aと比べると、色付きが残留していることがわかる。このように、画像回復フィルタ適用後、色付き抑制処理の要否判定前の画像で倍率色収差を検出し、補正を行うことで、色付きをより適切に抑制することができる。以上、本発明の第2の実施形態について説明した。このように、画像回復フィルタを適用した後に、画像から倍率色収差を検出して補正する場合でも本発明は有効である。図13cに、図13aに対応する非合焦時の状態を示す。入力画像には倍率色収差によりG成分とR成分の色ズレが発生し、さらにG成分とR成分のぼけ方は同じであるが、左右のエッジの収差が異なるために左右のエッジでぼけ方が異なっている状態である。S402で倍率色収差補正した後の倍率色収差補正後1では、R成分とG成分の位置が揃い、G成分とR成分のぼけ方が同じであるため、この段階で色にじみは発生していない。画像回復フィルタ適用後は、G成分、R成分ともにぼけが低減されているが、G成分とR成分のぼけの残り方が異なるために色にじみが発生している。この色にじみは左右のエッジの状態が異なり、色ズレとして扱うことができる。S406において、画像回復フィルタ適用後の画像から倍率色収差を検出し、補正した結果を示す倍率色収差補正後2では、左右のエッジの色にじみ方が同じになっている。この場合の画像回復処理前後の色差情報の関係を図14cに示す。画像回復フィルタの適用前の倍率色収差補正後1での色差は全領域で0となっているが、画像回復フィルタの適用後の倍率色収差補正後2はエッジ部分で色差が発生している。従って、S408における判定処理では色付き抑制のための画素値調整が必要と判定され、S409で色付き抑制のための画素値調整が行なわれる。S409では、画像回復フィルタの適用前の倍率色収差補正後1の色差が使用されるため、画素値調整後の色差は全領域で0となり、色差が適切に抑制された画像が画像回復処理結果として得られる。本実施形態と異なり、画像回復フィルタに倍率色収差補正成分を含めた場合の例を図13dに示す。図13dにおいて、入力画像は図13cと同じである。S402での倍率色収差補正を行わずに画像回復フィルタで倍率色収差補正も同時に行うため、画像回復フィルタ適用後の画像と、画像回復フィルタ適用後の画像から倍率色収差を検出して補正した倍率色収差補正後2の画像は、図13cと同様である。この場合の画像回復処理前後の色差情報の関係を図14dに示す。左右のエッジとも一部で色差の反転があるため、S408における判定処理では色付き抑制のための画素値調整が必要と判定され、S409で色付き抑制のための画素値調整が行なわれる。この場合の出力画像は図13dの画素値調整後の状態である。色付き抑制処理によって色付きは抑制されているが、本実施形態による図13cと比べると色付きが残っていることが確認される。以上説明したように、本実施形態によれば、倍率色収差を補正しない画像回復フィルタ適用後、色付き抑制のための画素値調整の有無を判定する前に、画像回復フィルタ適用後の画像から倍率色収差を検出して補正する。これにより、一層良好な色付き抑制効果を実現することが可能になる。その他の実施形態上述の実施形態においては、画像回復処理前後における画像の色に関する特徴量として、色差を用いた場合を説明した。しかし、色差に限らず、画像回復フィルタの適用によって生じた色付きを補正するために用いることが可能な任意の色に関する特徴量を用いることができる。例えば、画像回復処理前後における画像にて、Gプレーンに対しRプレーンもしくはBプレーンをメリジオナル方向に移動させながら、Yプレーンで検出されたエッジ付近の画素において、色成分間の差分絶対値和が最小となる場所を探索する。ここで検出された差分絶対値和が最小となった位置と、移動前の位置とから、Gプレーンに対するRプレーンもしくはBプレーンの色ズレ方向と色ズレ量を取得し、この色ズレ方向と色ズレ量を色に関する特徴量とすることができる。画像回復処理前後における画像の色ズレ方向が反転していれば、色付き抑制のための画素値調整を行い、色ズレ方向が反転していなくとも画像回復処理後の色ズレ量ほうが大きいのであれば、色付き抑制のための画素値調整を行うようにすればよい。また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェアプログラムを、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータまたはCPUやMPU等がプログラムを読み出して実行する処理である。"