"撮影画像中で輝度値が所定の値より大きな値となる位置である輝度飽和位置を検出する輝度飽和位置検出部と、前記輝度飽和位置の周辺部の輝度傾斜を検出する輝度傾斜検出部と、前記輝度飽和位置の画像と、前記輝度飽和位置に対応するPSFPointSpreadFunction画像と、前記輝度傾斜とに基づいて、前記輝度傾斜が大きいほど輝度値が大きくなるように、前記輝度飽和位置に撮像された被写体によって形成される撮像面上の輝度分布を推定する光源像推定部と、前記輝度分布を用いて、前記撮影画像から不要光の輝度値を減算する不要光減算部とを備える画像処理装置。前記光源像推定部は、前記輝度飽和位置の画像と前記輝度飽和位置に対応するPSF画像との畳込み積分を行うことにより光源像モデルを作成する光源像モデル作成部と、前記輝度傾斜が大きいほど輝度値が大きくなるように前記光源像モデルの輝度値を調整することにより、前記撮像面上の輝度分布を推定する光源像ゲイン調整部とを備える請求項1に記載の画像処理装置。前記光源像ゲイン調整部は、撮影画像の輝度傾斜と撮像面上の輝度最大値との予め定められた関係を利用して、検出された前記輝度傾斜に対応する前記撮像面上の輝度最大値を推定し、推定した前記輝度最大値を用いて前記光源像モデルの輝度値を調整する請求項2に記載の画像処理装置。前記撮影画像は、回折光学素子を含む光学系を用いて撮影された画像であり、前記輝度傾斜検出部は、前記輝度飽和位置の周辺部のうち前記光学系の光軸と反対側の輝度傾斜を検出する請求項1~3のいずれか1項に記載の画像処理装置。前記画像処理装置は、集積回路として構成されている請求項1~4のいずれか1項に記載の画像処理装置。請求項1~5のいずれか1項に記載の画像処理装置と、光学系及び撮像素子を有し、前記撮影画像を出力する撮像部とを備える撮像装置。撮影画像中で輝度値が所定の値より大きな値となる位置である輝度飽和位置を検出する輝度飽和位置検出ステップと、前記輝度飽和位置の周辺部の輝度傾斜を検出する輝度傾斜検出ステップと、前記輝度飽和位置の画像と、前記輝度飽和位置に対応するPSFPointSpreadFunction画像と、前記輝度傾斜とに基づいて、前記輝度傾斜が大きいほど輝度値が大きくなるように、前記輝度飽和位置に撮像された被写体によって形成される撮像面上の輝度分布を推定する光源像推定ステップと、前記輝度分布を用いて、前記撮影画像から不要光の輝度値を減算する不要光減算ステップとを含む画像処理方法。請求項7に記載の画像処理方法をコンピュータに実行させるためのプログラム。特許文献1のように、不要光として、設計次数以外の不要回折光0次光や2次光のみに着目して画像処理を行った場合、次数に依存せずに発生する不要光を撮影画像から低減することができない。さらに、PSF画像と撮影画像の畳込み積分により不要光を推定する場合、被写体が明るいことにより撮影画像の輝度が飽和していると、輝度飽和位置では誤った被写体の明るさ輝度から不要光を推定することになる。そのため、不要光を正確に推定することが難しい。被写体の明るさが、撮影画像の飽和レベルより明るくなるほど、推定した不要光成分の大きさが実際より小さくなる。また、次数に依存せずに発生する不要光は、被写体が明るい程より顕著となる。つまり、被写体が明るく、撮影画像の輝度が飽和する場合には、適切に不要光を低減することが困難である。また、特許文献2では、露光時間の異なる複数の画像を撮影する必要があるため、撮影するタイミングによって被写体が動いて撮影画像が変動する動画撮影の場合などにおいて、不要光を正確に推定することが困難となる場合がある。また、複数の画像を演算する必要があるので、演算コストも増大する。そこで、本発明は、上記課題を解決するためになされたものであって、明るい被写体が撮影された1枚の撮影画像を用いて、当該撮影画像よりも不要光が適切に低減された出力画像を生成することができる画像処理装置、撮像装置、及び画像処理方法を提供することを目的とする。上記の目的を達成するために、本発明の一態様に係る画像処理装置は、撮影画像中で輝度値が所定の値より大きな値となる位置である輝度飽和位置を検出する輝度飽和位置検出部と、前記輝度飽和位置の周辺部の輝度傾斜を検出する輝度傾斜検出部と、前記輝度飽和位置の画像と、前記輝度飽和位置に対応するPSFPointSpreadFunction画像と、前記輝度傾斜とに基づいて、前記輝度傾斜が大きいほど輝度値が大きくなるように、前記輝度飽和位置に撮像された被写体によって形成される撮像面上の輝度分布を推定する光源像推定部と、前記輝度分布を用いて、前記撮影画像から不要光の輝度値を減算する不要光減算部とを備える。また、本発明の一態様に係る撮像装置は、上記画像処理装置と、光学系及び撮像素子を有し、前記撮影画像を出力する撮像部とを備える。なお、本発明は、このような画像処理装置として実現することができるだけでなく、このような画像処理装置が備える特徴的な構成部の動作をステップとする画像処理方法として実現することができる。また、本発明は、画像処理方法に含まれる各ステップをコンピュータに実行させるプログラムとして実現することもできる。そして、そのようなプログラムは、CD-ROMCompactDiscReadOnlyMemory等の非一時的な記録媒体あるいはインターネット等の伝送媒体を介して配信することができるのは言うまでもない。本発明によれば、明るい被写体が撮影された1枚の撮影画像を用いて、当該撮影画像よりも不要光が適切に低減された出力画像を生成することが可能となる。本発明の一実施形態に係る画像処理装置は、撮影画像中で輝度値が所定の値より大きな値となる位置である輝度飽和位置を検出する輝度飽和位置検出部と、前記輝度飽和位置の周辺部の輝度傾斜を検出する輝度傾斜検出部と、前記輝度飽和位置の画像と、前記輝度飽和位置に対応するPSF画像と、前記輝度傾斜とに基づいて、前記輝度傾斜が大きいほど輝度値が大きくなるように、前記輝度飽和位置に撮像された被写体によって形成される撮像面上の輝度分布を推定する光源像推定部と、前記輝度分布を用いて、前記撮影画像から不要光の輝度値を減算する不要光減算部とを備える。この構成によれば、撮像画像の輝度傾斜が大きいほど輝度値が大きくなるように、輝度飽和位置に撮像された被写体によって形成される撮像面上の輝度分布を推定することができる。つまり、輝度が飽和していることにより撮影画像から被写体の正しい輝度分布が得られない場合であっても、輝度傾斜を用いて撮像面上の輝度分布を正確に推定することができる。このように推定された撮像面上の輝度分布を用いて撮影画像から不要光の輝度値を減算することにより、撮影画像よりも不要光が適切に低減された出力画像を1枚の撮影画像を用いて生成することが可能となる。つまり、明るい被写体が撮影された撮影画像において被写体像の輝度が飽和している場合であっても、撮影画像から不要光を適切に低減することが可能となる。さらに、この構成によれば、PSF画像に基づいて推定された撮像面上の輝度分布を用いて、撮影画像から不要光の輝度値を減算することができる。つまり、設計次数以外の不要回折光のみに着目して画像処理を行うわけではないので、次数に依存せずに発生する不要光を撮影画像から低減することも可能となる。また、本発明の別の実施形態に係る画像処理装置において、前記光源像推定部は、前記輝度飽和位置の画像と前記輝度飽和位置に対応するPSF画像との畳込み積分を行うことにより光源像モデルを作成する光源像モデル作成部と、前記輝度傾斜が大きいほど輝度値が大きくなるように前記光源像モデルの輝度を調整することにより、前記撮像面上の輝度分布を推定する光源像ゲイン調整部とを備える。この構成によれば、輝度飽和位置の画像とPSF画像との畳込み積分を行うことにより光源像モデルを作成することができる。このように作成された光源像モデルの輝度値を、輝度傾斜を用いて調整することにより、撮像面上の輝度分布を正確に推定することが可能となる。また、本発明の別の実施形態に係る画像処理装置において、前記光源像ゲイン調整部は、撮影画像の輝度傾斜と撮像面上の輝度最大値との予め定められた関係を利用して、検出された前記輝度傾斜に対応する前記撮像面上の輝度最大値を推定し、推定した前記輝度最大値を用いて前記光源像モデルの輝度値を調整する。この構成によれば、撮影画像の輝度傾斜と撮像面上の輝度最大値との予め定められた関係を利用して輝度最大値を正確に推定することができる。このように推定された輝度最大値を用いて光源像モデルの輝度を調整することにより、被写体によって形成される撮像面上の輝度分布を正確に推定することができる。したがって、撮影画像よりも不要光が適切に低減された出力画像を生成することが可能となる。また、本発明の別の実施形態に係る画像処理装置において、前記撮影画像は、回折光学素子を含む光学系を用いて撮影された画像であり、前記輝度傾斜検出部は、前記輝度飽和位置の周辺部のうち前記光学系の光軸と反対側の輝度傾斜を検出する。この構成によれば、輝度飽和位置の周辺部のうち光学系の光軸と反対側の輝度傾斜を検出することができる。回折光学素子を含む光学系を用いて撮影された画像では、輝度飽和位置の光軸側に大きな不要光が発生する。つまり、光軸と反対側の輝度傾斜は、光軸側の輝度傾斜よりも不要光の影響が少ない。したがって、光軸と反対側の輝度傾斜を検出することにより、検出される輝度傾斜に対する不要光の影響を抑制することができる。その結果、撮像面上の輝度分布をさらに正確に推定することが可能となる。また、本発明の別の実施形態に係る画像処理装置は、集積回路として構成されてもよい。以下、本発明の実施の形態の詳細を説明する前に、明るい被写体を撮影した場合に、被写体像が飽和し不要光が視認可能になる状態について、回折光学素子を含む光学系を使用して撮影された撮影画像を用いて説明する。図1は、本発明の実施の形態における光学系の構成例を模式的に示している。光学系200は、負のパワーを持つレンズ201と、正のパワーを持つ回折レンズ202とを備える。光学系200の光軸210は、撮像素子208の撮像面209と交差している。回折レンズ202は、回折光学素子に相当する。回折レンズ202は、互いに異なる材料からなる第1部材203及び第2部材204を備える。第1部材203の一方の面は、非球面形状に形成されている。また、第1部材203の他方の面には、回折格子206が光軸を中心に輪帯状に形成されている。回折格子206の表面は、第2部材204によって非球面形状に被覆されている。このような光学系200を介して、被写体の像が撮像素子208の撮像面209に結像される。撮像面に結像された被写体の像は、撮像素子208に撮影画像として取り込まれる。撮像素子208は、CCDChargeCoupledDeviceやCMOSComplementaryMetalOxideSemiconductorなどにより構成される。絞り211は、撮像素子208の撮像面209に入射する光線を調整する。回折格子206の格子厚dは、式1により求められる。ここで、n1は第1部材203の屈折率を表わし、n2は第2部材204の屈折率を表わす。λは波長である。ここでは、光学系200は、撮像用途の光学系である。そのため、λは、400nmから700nm程度の可視域の波長範囲の値をとる。mは回折次数である。ここでは、m=1としている。つまり、回折格子206の格子厚dは、1次回折光の回折効率が高くなるように設計されている。この可視域の波長範囲においてdがほぼ一定値となる屈折率n1及びn2を持つ第1部材203及び第2部材204を組合せることにより、回折レンズ202は、可視域の波長範囲全領域で高い1次回折効率を実現することができることが知られている例えば参考文献特許第4077508号公報を参照。本発明の実施の形態においては、n1及びn2がこのような組合せとなる第1部材203及び第2部材204が用いられる。図2Aに、図1の光学系200の画角45度付近でのPSF画像を示す。図2Bに、PSF画像の不要光成分の分布が見やすいように、図2AのPSF画像の明るさを50倍したPSF画像を示す。図2A及び図2Bでは画像の左方向が光軸方向である。つまり、画像の左方向に、光軸と撮像面とが交差する位置に対応する画像位置が存在する。なお、PSFとは、光学系の点光源に対する応答を表す関数である。また、PSF画像とは、PSFを表す画像である。つまり、PSF画像とは、点光源が撮影された画像に相当する。図3Aは、図2AのPSF画像の輝度最大位置周辺部の画像水平方向の輝度推移を表わす。図3Aにおいて、縦軸は輝度値を表わし、横軸は画像位置を表わしている。図3Bは、図3Aの縦軸のスケールを拡大したときの輝度推移を表わす。図2A~図3Bに示すように、回折光学素子を含む光学系を用いて撮影された撮影画像では、大きな不要光が発生する。なお、不要光とは、画像に写ってしまう不要な光である。つまり、不要光とは、本来、画像に写ることが望ましくない光である。この不要光によって画質が劣化する。不要光には、0次、あるいは2次の回折光などの設計次数以外の回折光以下、不要回折光ともいうだけでなく、設計次数である1次回折光の不要光も含まれる。この1次回折光の不要光は、回折格子の各輪帯の影響により発生する不要光である。つまり、設計次数である1次回折光においても非球面レンズに比較して大きな不要光が原理的に発生する。特に、回折格子の格子厚dが大きい場合や、明るい被写体を撮影する場合には、図3Bに示すように、PSF画像の最大輝度の位置より光軸側の位置に、この大きな不要光が顕著に現れる。例えば、明るい被写体が撮影された撮影画像において被写体像の位置の輝度が飽和したときに、被写体像の位置の光軸側の位置に大きな不要光が発生する。不要光は、撮影画像に悪影響を及ぼす。以下に、このような不要光が発生したときの撮影画像の一例を、図4~図6Bを用いて説明する。図4は、本発明の実施の形態における被写体を示す。図5A及び図5Bは、本発明の実施の形態における撮影画像を示す。具体的には、図5A及び図5Bは、図1の光学系で図4に示す明るい被写体を撮影した場合の撮影画像を示す。本実施の形態では、図5A及び図5Bは、画角45度付近に撮像される蛍光灯の撮影画像を示す。蛍光灯は、図4に示すような長方形の輝度が均一な光源と仮定している。そして、図4に示す被写体の画像と図2AのPSF画像とを畳込み積分することにより、図5A及び図5Bに示す実際の撮像画像が模擬的に作成される。図5Aは、蛍光灯の最大輝度値が0.7となる場合の撮影画像を示す。図5Bは、蛍光灯の最大輝度値が1.3となる場合の撮影画像を示す。なお、図5A及び図5Bでは、撮影画像の輝度飽和レベルは1.0としている。図5Bでは、撮像面209上で輝度値が1.0を超える位置に対応する画像位置で、撮影画像の輝度が飽和している。図6A及び図6Bは、図5A及び図5Bそれぞれの蛍光灯の撮影画像の中心付近における画像水平方向の輝度推移を示す。図6A及び図6Bにおいて、縦軸は輝度値を表わし、横軸は画像位置を表わす。図5A~図6Bから分かるように、蛍光灯の最大輝度値が大きくなるほど、蛍光灯の光軸方向画像の左方向に二重像のように写り込む不要光の輝度値も大きくなる。つまり、図5Aでは目立たない不要光が、図5Bでは視認できるレベルの明るさになっている。このことから、蛍光灯の最大輝度値が大きくなるほど不要光の輝度値が大きくなり、撮影画像に悪影響を及ぼすことが分かる。図5Bのように、撮影画像の輝度が飽和している画像位置以下、輝度飽和位置というにおいては、撮像素子208の撮像面209における実際の輝度分布が損なわれる。つまり、撮影画像の輝度飽和位置では、撮像面209上における実際の輝度分布とは異なる輝度分布が形成される。そのため、特許文献1の手法を用いた場合には、輝度飽和位置から得られる誤った被写体の明るさから不要光が推定される。その結果、特許文献1の手法では、不要光を適切に推定することが難しく、撮影画像から不要光を低減させることが困難である。そこで、以下に説明する本発明の実施の形態では、図5B及び図6Bに示したような、撮影画像の輝度が飽和している場合に発生する不要光を低減することができる画像処理装置及び撮像装置について説明する。なお、以下で説明する実施の形態は、いずれも本発明の好ましい一具体例を示す。つまり、以下の実施の形態で示される数値、形状、材料、構成要素、構成要素の配置及び接続形態、ステップ、ステップの順序などは、本発明の一例であり、本発明を限定する主旨ではない。本発明は、請求の範囲の記載に基づいて特定される。したがって、以下の実施の形態における構成要素のうち、本発明の最上位概念を示す独立請求項に記載されていない構成要素は、本発明の課題を達成するために必ずしも必要ではないが、より好ましい形態を構成する構成要素として説明される。実施の形態以下、本発明の実施の形態について図面を参照しながら説明する。図7Aは、本発明の実施の形態における撮像装置100の構成を示すブロックずである。また、図7Bは、本発明の実施の形態における光源像推定部の構成の一例を示すブロック図である。図7Aにおいて、撮像装置100は、撮像部101と画像処理装置102とを備える。撮像部101は、図1に示す光学系200及び撮像素子208を備える。撮像部101は、被写体を撮像し、撮影画像Itx,yを出力する。ここで、x及びyは、それぞれ画像水平方向及び垂直方向の画像位置を表わす。なお、画像位置とは、撮影画像上の位置である。例えば、画像位置は、撮影画像を構成する各画素の位置を示す。画像処理装置102は、撮影画像から不要光の成分を低減し、不要光による画質の劣化が少ない出力画像を出力する。具体的には、画像処理装置102は、撮影画像中の被写体像のうち所定の輝度閾値Is以上の被写体像を光源像とみなし、その光源によって形成される撮像面209上の輝度分布を推定する。そして、画像処理装置102は、推定された撮像面209上の輝度分布に基づいて撮影画像から不要光の輝度値を減算することにより、撮影画像よりも不要光による画質の劣化が少ない出力画像を生成する。なお、不要光の輝度値とは、不要光によって形成される像不要光像の輝度値である。なお、本実施の形態では、画像処理装置102が撮像装置100に備えられる場合について説明するが、画像処理装置102は、必ずしも撮像装置100に備えられる必要はない。例えば、画像処理装置102は、撮像部101を備える撮像装置から撮影画像を取得してもよい。以下に、画像処理装置102の詳細について説明する。図7Aに示すように、画像処理装置102は、輝度飽和位置検出部110と、輝度傾斜検出部112と、光源像推定部105と、不要光減算部106とを備える。輝度飽和位置検出部110は、撮影画像中で輝度値が所定の値より大きな値となる位置である輝度飽和位置を検出する。つまり、輝度飽和位置検出部110は、撮影画像において、輝度値が閾値より大きい画像位置を輝度飽和位置として検出する。輝度傾斜検出部112は、輝度飽和位置の周辺部の輝度傾斜を検出する。輝度傾斜とは、撮影画像における輝度値の空間的な変化の割合を示す。光源像推定部105は、輝度飽和位置の画像と、輝度飽和位置に対応するPSF画像と、輝度傾斜とに基づいて、輝度傾斜が大きいほど輝度値が大きくなるように、撮像面209上の光源像を推定する。光源像とは、光源によって形成される輝度分布である。また、光源とは、輝度飽和位置に撮像された被写体に相当する。本実施の形態では、図7Bに示すように、光源像推定部105は、光源像モデル作成部111と、畳込み積分部115及びPSF抽出部116を有する光源像ゲイン調整部113とを備える。なお、図7Bに示す光源像推定部105の構成は一例であり、光源像推定部105は、必ずしもこのように構成される必要はない。不要光減算部106は、推定された撮像面209上の光源像を用いて、撮影画像から不要光の輝度値を減算する。つまり、不要光減算部106は、推定された撮像面209上の光源像から得られる不要光成分を撮影画像から減算することにより、撮影画像よりも不要光による画像劣化が少ない出力画像を生成する。次に、輝度飽和位置検出部110と、光源像推定部105に含まれる光源像モデル作成部111との処理を詳細に説明する。輝度飽和位置検出部110は、撮像部101が出力する撮影画像において輝度閾値Isより大きな輝度値を持つ画像位置を、輝度飽和位置として検出する。そして、輝度飽和位置検出部110は、検出した輝度飽和位置を示すデータをメモリ等に記憶する。輝度閾値Isは、例えば、撮影画像において表現可能な輝度値の範囲が0~1.0であれば、0.98程度に設定される。輝度閾値Isは、撮像素子208の撮像特性などに応じて設定されればよい。つまり、輝度飽和位置検出部110は、輝度が飽和している画像位置に撮像された被写体を光源と仮定して検出する。実際の撮影環境においては、蛍光灯やランプなど、他の被写体に比べて非常に明るい被写体の輝度が飽和することが多い。したがって、輝度飽和位置に撮像された被写体を光源と仮定することは妥当性を損なわない。図8は、図5Bの撮影画像から検出された輝度飽和位置を示す。次に、光源像モデル作成部111は、輝度飽和位置の画像と輝度飽和位置に対応するPSF画像との畳込み積分を行うことにより光源像モデルを作成する。光源像モデル作成部111は、畳込み積分部115とPSF抽出部116とを備える。畳込み積分部115は、輝度飽和位置の輝度値を一定値Icと設定し、かつ、その他の画像位置の輝度値を0と設定した光源形状画像Ifx,yを作成する。そして、畳込み積分部115は、PSF抽出部116によって抽出された画像位置毎のPSF画像Ipsfx,yと光源形状画像Ifx,yとの畳込み積分を式2に示すように行うことにより、光源像モデルImx,yを作成する。ここで、u及びvは、それぞれ画像水平方向及び垂直方向の画像位置を表わす。なお、畳込み積分部115は、Imx,yを畳込み積分により求めたが、必ずしもこのように求める必要はない。例えば、畳込み積分部115は、以下のようにImx,yを求めてもよい。まず、畳込み積分部115は、Ifx,y及びIpsfx,yをそれぞれFFTFastFourierTransformなどを用いてフーリエ変換する。そして、畳込み積分部115は、フーリエ変換後のデータを周波数領域上で乗算する。最後に、畳込み積分部115は、乗算後のデータを逆フーリエ変換することにより光源像モデルImx,yを算出する。なお、畳込み積分部115は、演算量などを考慮して空間領域で演算を行うか周波数領域で演算を行うかを決めればよい。PSF抽出部116は、メモリ等に予め保存された複数のPSF画像の中から各画像位置に対応するPSF画像を抽出する。なお、PSF画像Ipsfx,yが画像位置毎に予め保存された場合、膨大なメモリ量が必要となる。そこで、PSF抽出部116は、画像位置の集合であるブロック例えば、64×64画素ごとに保存されたPSF画像の中から、各画像位置が属するブロックのPSF画像を、各画像位置に対応するPSF画像として抽出してもよい。この場合、例えば、正方形のブロック、長方形のブロック、もしくは光軸を中心とした輪帯状のブロックごとにPSF画像が保存されればよい。これにより、PSF画像を保存するためのメモリ量を削減することができる。なお、ブロックサイズを大きくするほど、PSF画像の保存に必要なメモリ量は削減できるが、実際のPSF画像と抽出されたPSF画像とのかい離が大きい画像位置が増加してしまう。その結果、光源像の推定精度が低下してしまう。したがって、撮像装置に搭載するメモリ量と光源像の推定精度とのバランスを考慮して、ブロックサイズが決定されることが好ましい。図9Aは、図2Aに示したPSF画像及び図8の光源形状画像から式2を用いて作成された光源像モデルを示す。図9Bは、図9Aの光源像モデルの輝度推移を示す。なお、光源像モデルは、ここでは輝度最大値が1.0となるように正規化されている。正規化は、実装時のシステム仕様に応じて実施されればよい。次に、輝度傾斜検出部112の処理について説明する。輝度傾斜検出部112は、輝度飽和位置検出部110が検出した輝度飽和位置の周辺部における撮影画像の輝度傾斜Dmを検出する。具体的には、まず、輝度傾斜検出部112は、隣り合う画像位置の輝度値の差分の絶対値を微分値として微分演算により算出する。図10は、図6Bの輝度推移から微分演算により算出した、輝度飽和位置の周辺部における輝度値の微分値の推移を示す。ここでは、輝度傾斜検出部112は、輝度飽和位置の光軸側の微分値の最大値を輝度傾斜Dmとして検出する。なお、輝度傾斜Dmは、必ずしも輝度飽和位置の周辺部の微分値の最大値である必要はない。例えば、輝度傾斜Dmは、輝度飽和位置の周辺部の微分値の平均値であってもよい。次に、光源像推定部105に含まれる光源像ゲイン調整部113の処理を説明する。図11は、輝度傾斜Dmと撮像面209での光源像の実際の輝度最大値との関係を示す。図11において、第1データ301は、図4に示した蛍光灯を撮影したときの、輝度傾斜Dmと撮像面209での実際の輝度最大値との関係を示す。第1データ301から、輝度傾斜Dmと輝度最大値とは、ほぼ比例の関係にあることが分かる。また、図11において、第2データ302は、図4に示した蛍光灯の幅を約1/2倍にした蛍光灯を撮影したときの、輝度傾斜Dmと輝度最大値との関係を示す。第2データ302から、蛍光灯の幅を1/2倍にしても、輝度傾斜Dmと輝度最大値との関係は大きく変化しないことがわかる。すなわち、第1データ301と第2データ302との間で、同一の輝度傾斜Dmに対する輝度最大値は1割程度しか変化していない。つまり、光源形状の変化に対して輝度傾斜Dmがあまり変動しないことが分かる。したがって、光源像ゲイン調整部113は、輝度飽和位置の周辺部の輝度傾斜Dmを検出することにより、光源形状に関わらず、光源の撮像面209での実際の輝度最大値を推定することが可能となる。このように、撮影画像の一部領域において輝度が飽和している場合でも、飽和している領域の周辺部の輝度傾斜は、飽和している領域に撮像された光源の撮像面209上の実際の輝度最大値に比例している。さらに、被写体の形状が変わっても飽和している領域の周辺部の輝度傾斜の大きさは、大きく変化しない。したがって、光源像ゲイン調整部113は、輝度傾斜Dmと輝度最大値との関係を予め記憶しておけば、輝度傾斜Dmから輝度最大値を精度よく推定することができる。そして、このように推定された輝度最大値を用いて光源像モデルの輝度値を調整することにより、画像処理装置102は、撮影画像から不要光を低減することができる。具体的には、光源像ゲイン調整部113は、式3を利用して、検出した輝度傾斜Dmから撮像面209での輝度最大値Imaxを算出する。なお、この式3が、撮影画像の輝度傾斜と撮像面209上の輝度最大値との予め定められた関係に相当する。本実施の形態においては、式3におけるパラメータA及びBを、A=2.6153、B=0とした。これらのパラメータA及びBの値は、輝度傾斜Dmから得られる輝度最大値Imaxが、実際の輝度最大値図11の第1データ301より1割程度小さくなるように設定されている。なお、輝度傾斜Dmから得られる輝度最大値Imaxが、実際の輝度最大値図11の第1データ301と同一となるようにパラメータA及びBの値を設定してもよい。しかし、図11の第2データ302のように、光源形状が変わった場合に、輝度傾斜Dmに対する実際の輝度最大値も若干変動する可能性がある。そこで、上記のように、パラメータA及びBの値は、輝度最大値Imaxが第1データ301によって示される実際の輝度最大値より若干小さめとなるように設定されることが好ましい。これにより、光源像ゲイン調整部113は、光源形状が変化したときに輝度最大値を実際より大きく推定してしまうことを抑制することができる。その結果、不要光減算部106は、撮影画像から過剰に不要光成分を減算することを抑制することができる。より具体的には、本実施の形態では、最小二乗法などを用いて図11の第1データ301から得られる近似直線よりも、式3のパラメータBの値が小さめに設定されている。これにより、輝度最大値Imaxが実際の輝度最大値より大きく推定されることが抑制される。光源像ゲイン調整部113は、式3に基づき算出した輝度最大値Imaxを光源像モデルImx,yに乗じることにより、推定光源像Iex,yを作成する。つまり、光源像ゲイン調整部113は、輝度飽和位置の輝度値が輝度最大値Imaxと一致するように、光源像モデルImx,yの輝度値を調整することにより、撮像面209上の光源像を推定する。なお、ここでは、輝度傾斜Dmと輝度最大値Imaxとの予め定められた関係は、1次関数で表わされたが、必要に応じて2次関数など適当な多項式で表わされてもよい。また、輝度傾斜Dmと輝度最大値Imaxとの関係は、輝度傾斜Dmに輝度最大値Imaxが対応付けられたテーブルで表わされてもよい。この場合、光源像ゲイン調整部113は、そのテーブルを参照して、輝度傾斜Dmに対応する輝度最大値Imaxを推定すればよい。なお、輝度傾斜Dmと輝度最大値Imaxとの関係が数式として記憶された場合の方が、テーブルとして記憶される場合よりも、メモリ量を低減することができる。また、輝度傾斜Dmと輝度最大値Imaxとの予め定められた関係は、画像位置ごとに記憶されていることが好ましい。これにより、光源像ゲイン調整部113は、撮影画像中の輝度飽和位置に応じて輝度最大値をさらに正確に推定することが可能となる。なお、輝度傾斜Dmと輝度最大値Imaxとの予め定められた関係は、画像位置の集合であるブロック例えば64×64画素ごとに記憶されていてもよい。この場合、輝度傾斜Dmと輝度最大値Imaxとの予め定められた関係が画像位置ごとに記憶されている場合よりも、メモリ量を低減することができる。以上のように、光源像ゲイン調整部113は、輝度傾斜と輝度最大値との予め定められた関係を利用して、検出された輝度傾斜に対応する撮像面209上の輝度最大値を推定する。そして、光源像ゲイン調整部113は、推定した輝度最大値を用いて光源像モデルの輝度値を調整する。このように輝度値が調整された推定光源像Iex,yは、撮影画像中で輝度が飽和することにより失われた、光源によって形成された撮像面209上の輝度分布を復元したものとみなすことができる。次に、不要光減算部106の処理について説明する。不要光減算部106は、輝度飽和位置以外の画像位置において、撮影画像Itx,yから推定光源像Iex,yを減算することにより、撮影画像よりも不要光が低減された出力画像を生成する。つまり、不要光減算部106は、推定光源像のうち輝度飽和位置以外の画像位置の輝度値を、撮影画像の対応する画像位置の輝度値から減算する。図12Aは、本発明の実施の形態における出力画像の一例を示す。図12Bは、図12Aに示す出力画像の中心付近の画像水平方向の輝度推移を示す。図12Bにおいて、実線311は出力画像の輝度推移を表わし、点線312は撮影画像Itx,yの輝度推移を表わし、破線313は推定光源像Iex,yから得られる不要光の輝度推移撮影画像からの減算量を表わす。減算量は、Itx,yの輝度が飽和している位置では0となっている。図12Aより、出力画像では撮影画像の不要光が大幅に低減されているのが分かる。次に、上記と被写体の形状が異なる場合について説明する。図13Aは、本発明の実施の形態における被写体の他の一例を示す。具体的には、図13Aは、図4で示した被写体よりも幅が約1/2の被写体を示す。より具体的には、図13Aは、図4で示した蛍光灯の幅を約1/2倍にした光源を示す。図13Bは、本発明の実施の形態における撮影画像の他の一例を示す。具体的には、図13Bは、図13Aに示す被写体を撮影したときの撮影画像Itx,yを示す。図13Bの撮影画像では、図5Bと同様に、被写体が撮像された画像位置の輝度が飽和している。図13Cは、本発明の実施の形態における出力画像の他の一例を示す。具体的には、図13Cは、図13Bに示す撮影画像から不要光を低減した後の出力画像を示す。図13Dは、図13Cの出力画像の画像中央部付近の画像水平方向の輝度推移を示す。図13Dにおいて、縦軸が輝度値を表わし、横軸が画像位置を表わす。図12Bと同様に、実線321は出力画像の輝度推移を表わし、点線322は撮影画像Itx,yの輝度推移を表わし、破線323は推定光源像Iex,yから得られる不要光の輝度推移撮影画像からの減算量を表わす。図14Aは、本発明の実施の形態における被写体の他の一例を示す。具体的には、図14Aは、ハロゲン電球などの円形の光源を示す。図14Bは、本発明の実施の形態における撮影画像の他の一例を示す。具体的には、図14Bは、図14Aに示す被写体を撮影したときの撮影画像Itx,yを示す。図14Bの撮影画像では、図5Bと同様に、被写体が撮像された画像位置の輝度が飽和している。図14Cは、本発明の実施の形態における出力画像の他の一例を示す。具体的には、図14Cは、図14Bに示す撮影画像から不要光を低減した後の出力画像を示す。図14Dは、図14Cの出力画像の画像中央部付近の画像水平方向の輝度推移を示す。図14Dにおいて、縦軸が輝度値を表わし、横軸が画像位置を表わす。図12Bと同様に、実線331は出力画像の輝度推移を表わし、点線332は撮影画像Itx,yの輝度推移を表わし、破線333は推定光源像Iex,yから得られる不要光の輝度推移撮影画像からの減算量を表わす。なお、図13C又は図14Cの出力画像を生成する際に、光源像ゲイン調整部113は、図12Aの出力画像を生成したときと同一のパラメータA及びBの値が設定された式3を用いて、輝度最大値を推定している。図13C又は図14Cの出力画像から分かるように、被写体の形状が変わっても、出力画像では撮影画像の不要光が大幅に低減されている。次に、以上のように構成された撮像装置100の各種動作について説明する。図15Aは、本発明の実施の形態における画像処理を示すフローチャートである。図15Aにおいて、撮像部101は、すでに、光学系200を介して撮像された撮影画像を作成しているとする。また、図15Bは、本発明の実施の形態における光源像推定処理の一例を示すフローチャートである。まず、輝度飽和位置検出部110は、撮影画像中の輝度飽和位置を検出するS102。輝度傾斜検出部112は、ステップS102で検出した輝度飽和位置の周辺部の輝度傾斜Dmを、微分演算を行うことにより検出するS103。光源像推定部105は、検出された輝度傾斜Dmに基づいて撮像面209上の光源像を推定するS104。本実施の形態におけるステップS104の処理の詳細の一例は、図15Bに示すとおりである。なお、ステップS104の処理は、必ずしも図15Bように行われる必要はない。つまり、ステップS104では、輝度飽和位置の画像と、輝度飽和位置に対応するPSF画像と、輝度傾斜とに基づいて、輝度傾斜が大きいほど輝度値が大きくなるように、輝度飽和位置に撮像された被写体によって形成される撮像面上の輝度分布が推定されればよい。以下に、図15Bに示すフローチャートについて説明する。まず、PSF抽出部116は、ステップS102で検出した輝度飽和位置に対応するPSF画像をメモリ等にあらかじめ保存された複数のPSF画像の中から抽出するS105。なお、PSF抽出部116は、必要に応じてメモリに保存していたPSF画像の回転処理を行うことにより、輝度飽和位置に対応するPSF画像を作成してもよい。あるいは、PSF抽出部116は、ブロック毎に同一のPSF画像をあらかじめ保存しておき、画素飽和位置が属するブロックのPSF画像を読み出してきてもよい。続いて、畳込み積分部115は、ステップS102で検出した輝度飽和位置に基づき作成した光源形状画像と、ステップS103で抽出したPSF画像の2つ画像の畳込み積分を行うことにより、光源像モデルを作成するS106。なお、畳込み積分部115は、畳込み積分の代わりに、2つの画像をそれぞれFFTなどによりフーリエ変換した後に乗算し、逆フーリエ変換することにより、光源像モデルを作成してもよい。次に、光源像ゲイン調整部113は、輝度傾斜が大きいほど輝度値が大きくなるように光源像モデルの輝度値を調整することにより、光源像の輝度分布を推定するS107。具体的には、光源像ゲイン調整部113は、撮影画像の輝度傾斜と撮像面上の輝度最大値との予め定められた関係を利用して、検出された輝度傾斜Dmに対応する撮像面上の輝度最大値Dmaxを推定する。そして、光源像ゲイン調整部113は、光源像モデルと輝度最大値Dmaxとに基づき、撮像面209上での光源の輝度分布を光源像として推定する。以上のように、ステップS104において光源像が推定された後、不要光減算部106は、撮影画像の輝度飽和位置以外の画像位置において、撮影画像の輝度値から、ステップS106で推定した撮像面209での光源像の輝度値を減算することにより、不要光が低減された出力画像を生成するS108。以上、本実施の形態によれば、明るい被写体を撮影した撮影画像に視認可能な不要光が発生する場合であっても、撮影画像の輝度飽和位置の画像とPSF画像とから光源像モデルを推定し、また、輝度飽和位置の周辺の輝度傾斜から推定された輝度最大値を用いて光源像モデルの輝度値を調整することにより、1回の撮影で得られた撮影画像を用いて当該撮影画像から適切に不要光を低減することが可能となる。つまり、本実施の形態における画像処理装置102によれば、撮像画像の輝度傾斜が大きいほど輝度値が大きくなるように、輝度飽和位置に撮像された被写体によって形成される撮像面209上の輝度分布を推定することができる。つまり、輝度が飽和していることにより撮影画像から被写体の正しい輝度分布が得られない場合であっても、輝度傾斜を用いて撮像面209上の輝度分布を正確に推定することができる。このように推定された撮像面209上の輝度分布を用いて撮影画像から不要光の輝度値を減算することにより、1枚の撮影画像を用いて当該撮影画像よりも不要光が適切に低減された出力画像を生成することが可能となる。さらに、画像処理装置102によれば、PSF画像に基づいて推定された撮像面209上の輝度分布を用いて、撮影画像から不要光の輝度値を減算することができる。つまり、設計次数以外の不要回折光のみに着目して画像処理を行うわけではないので、次数に依存せずに発生する不要光を撮影画像から低減することも可能となる。またさらに、画像処理装置102によれば、撮影画像の輝度傾斜と撮像面209上の輝度最大値との予め定められた関係を利用して輝度最大値を正確に推定することができる。このように推定された輝度最大値を用いて光源像モデルの輝度値を調整することにより、被写体によって形成される撮像面209上の輝度分布を正確に推定することができる。したがって、撮影画像よりも不要光が適切に低減された出力画像を生成することが可能となる。以上、本発明の一態様に係る撮像装置100及び画像処理装置102について、実施の形態に基づいて説明したが、本発明は、これらの実施の形態に限定されるものではない。本発明の趣旨を逸脱しない限り、当業者が思いつく各種変形を本実施の形態に施したものも、本発明の範囲内に含まれる。例えば、上記実施の形態において、輝度傾斜検出部112は、図10に示すように、輝度飽和位置の周辺部のうち光軸210側の輝度傾斜を検出していたが、光軸210とは反対側の輝度傾斜を検出してもよい。例えば、図10では、輝度傾斜検出部112は、輝度飽和位置の右側における微分値の最大値を輝度傾斜として検出してもよい。回折レンズ202を含む光学系200を用いて撮影された画像では、輝度飽和位置の光軸210側に大きな不要光が発生する。つまり、光軸210と反対側の輝度傾斜は、光軸210側の輝度傾斜よりも不要光の影響が少ない。したがって、輝度傾斜検出部112は、光軸210と反対側の輝度傾斜を検出することにより、検出される輝度傾斜に対する不要光の影響を抑制することができる。その結果、光源像推定部105は、撮像面209上の光源像をさらに正確に推定することが可能となる。また、上記実施の形態において、光源像推定部105は、輝度傾斜に基づいて光源像モデルの輝度値を調整していたが、必ずしも光源像モデルの輝度値を調整する必要はない。例えば、光源像推定部105は、輝度傾斜に基づいて、PSF画像の輝度値、又は撮影画像中の輝度飽和位置の画像の輝度値を調整してもよい。つまり、光源像推定部105は、必ずしも図7Bのように構成される必要はない。このような場合であっても、PSF画像の輝度値、又は撮影画像中の輝度飽和位置の画像の輝度値を、輝度傾斜が大きいほど大きくなるように調整することにより、光源像推定部105は、光源像モデルの輝度値を調整する場合と同様に光源像を正確に推定することができる。また、上記実施の形態において、光源像ゲイン調整部113は、輝度傾斜から輝度最大値を推定していたが、輝度最大値の代わりに、光源像モデルなどの輝度値を調整するための係数を推定してもよい。なお、上記実施の形態では、白黒画像の不要光低減について説明したが、例えばベイヤー配列のカラーフィルタを用いた撮像素子や、プリズムを用いた赤、緑、青用の3枚構成の撮像素子などを用いて撮影されたカラー用途の撮影画像の場合であっても、本発明を適用することができる。具体的には、画像処理装置102は、赤、緑、青それぞれの画像について、上記と同様に不要光を低減した画像を作成する。そして、画像処理装置102は、不要光を低減した赤、緑、青の画像を合成してカラー画像を作成すればよい。この時、画像処理装置102は、赤、緑、青それぞれ色毎に異なるPSF画像を用いて光源像を推定してもよい。また、画像処理装置102は、例えば緑など一色のみの光源像モデルを計算し、その光源像モデルを用いて輝度飽和位置周辺の輝度傾斜に基づき、赤、緑、青の各色の不要光を低減してもよい。このとき、画像処理装置102は、各色に応じて光源像の輝度最大値の比をあらかじめ決めておき、その比に基づいて各色で不要光の低減量を異なるように設定してもよい。なお、上記実施の形態では、回折光学素子を含む光学系を用いて撮影された撮影画像から不要光を低減する例を示した。しかし、回折光学素子を含まない光学系を用いて撮影された撮影画像であっても、画像処理装置102は、上記実施の形態と同様に撮影画像から不要光を低減することができる。例えば、光学系の収差が大きいなどの理由により撮影画像に不要光が発生する場合についても、飽和している画像領域の周辺部の輝度傾斜は、飽和している光源像の輝度最大値に比例し、かつ被写体が変わっても飽和している画像領域の周辺部の輝度傾斜の大きさは大きく変化しない。したがって、回折光学素子を含まない光学系を用いて撮影された撮影画像に対しても本発明が適用できることは言うまでもない。この場合、撮影画像の輝度傾斜と撮像面209上の輝度最大値との予め定められた関係を示す式3のパラメータA及びBは、光学系に応じて適当に設定されればよい。例えば、撮影画像の輝度傾斜と撮像面209上の輝度最大値との組み合わせを予め複数測定し、測定された複数の組み合わせにおいて最小二乗法などを用いて算出される近似直線に基づいて、パラメータA及びBが設定されればよい。また、上記実施の形態における画像処理装置102が備える構成要素の一部または全部は、1個のシステムLSILargeScaleIntegration:大規模集積回路から構成されているとしてもよい。例えば、画像処理装置102は、輝度飽和位置検出部110と、輝度傾斜検出部112と、光源像推定部105と、不要光減算部106とを備えるシステムLSIから構成されてもよい。システムLSIは、複数の構成部を1個のチップ上に集積して製造された超多機能LSIであり、具体的には、マイクロプロセッサ、ROMReadOnlyMemory、RAMRandomAccessMemoryなどを含んで構成されるコンピュータシステムである。前記ROMには、コンピュータプログラムが記憶されている。前記マイクロプロセッサが、前記コンピュータプログラムに従って動作することにより、システムLSIは、その機能を達成する。なお、ここでは、システムLSIとしたが、集積度の違いにより、IC、LSI、スーパーLSI、ウルトラLSIと呼称されることもある。また、集積回路化の手法はLSIに限るものではなく、専用回路又は汎用プロセッサで実現してもよい。LSI製造後に、プログラムすることが可能なFPGAFieldProgrammableGateArray、あるいはLSI内部の回路セルの接続や設定を再構成可能なリコンフィギュラブルプロセッサを利用してもよい。さらには、半導体技術の進歩又は派生する別技術によりLSIに置き換わる集積回路化の技術が登場すれば、当然、その技術を用いて機能ブロックの集積化を行ってもよい。バイオ技術の適用等が可能性としてありえる。また、本発明は、このような特徴的な処理部を備える画像処理装置として実現することができるだけでなく、画像処理装置に含まれる特徴的な処理部をステップとする画像処理方法として実現することもできる。また、画像処理方法に含まれる特徴的な各ステップをコンピュータに実行させるコンピュータプログラムとして実現することもできる。そして、そのようなコンピュータプログラムを、CD-ROM等のコンピュータ読取可能な非一時的な記録媒体あるいはインターネット等の通信ネットワークを介して流通させることができるのは、言うまでもない。"