(57)【要約】【課題】フレアの生源色に依らず適切なフレア補正処理を実行する。【解決手段】本発明の画像処理プログラムは、撮影光学系により取得されたカラーの処理対象画像を入力する入力手順と、前記処理対象画像におけるフレアの発生領域の色に基づき前記フレアの光源色が白色であったか否かを判別する判別手順（Ｓ９７）と、前記発生領域に対してフレア補正処理を施すフレア補正手順（Ｓ９３〜Ｓ９１０）とを含み、前記フレア補正手順では、前記判別手段による前記判別の結果に応じて前記フレア補正処理の内容を切り換える。【選択図】    図３
2014165833.tif 000002
【特許請求の範囲】
【請求項１】
撮影光学系により取得されたカラーの処理対象画像を入力する入力手段と、
前記処理対象画像におけるフレアの発生領域の色に基づき前記フレアの光源色が白色であったか否かを判別する判別手段と、
前記発生領域に対してフレア補正処理を施すフレア補正手段とを備え、
前記フレア補正手段は、
前記判別手段による前記判別の結果に応じて前記フレア補正処理の内容を切り換える
ことを特徴とする画像処理装置。
【請求項２】
請求項１に記載の画像処理装置において、
前記判別手段は、
前記撮影光学系の種類に基づき前記判別を行う
ことを特徴とする画像処理装置。
【請求項３】
請求項１又は請求項２に記載の画像処理装置において、
前記判別手段は、
前記発生領域の彩度に基づき前記判別を行う
ことを特徴とする画像処理装置。
【請求項４】
請求項３に記載の画像処理装置において、
前記判別手段は、
前記処理対象画像が色補正後の画像である場合には、前記発生領域の彩度に基づく代わりに、前記発生領域に色補正画素が存在するか否かに基づき前記判別を行う
ことを特徴とする画像処理装置。
【請求項５】
請求項１〜請求項４の何れか一項に記載の画像処理装置において、
前記判別手段は、
前記処理対象画像上の飽和画素領域を包含する円形領域を前記発生領域とみなす
ことを特徴とする画像処理装置。
【請求項６】
請求項５に記載の画像処理装置において、
前記判別手段は、
前記飽和画素領域のサイズに基づき前記発生領域のサイズを推定する
ことを特徴とする画像処理装置。
【請求項７】
請求項１〜請求項６の何れか一項に記載の画像処理装置において、
前記判別手段は、
前記発生領域から飽和画素領域を除外した領域のみを前記判別に使用する
ことを特徴とする画像処理装置。
【請求項８】
請求項１〜請求項７の何れか一項に記載の画像処理装置において、
前記フレア補正手段は、
前記発生領域を包含する処理領域を前記処理対象画像上に設定する領域設定手段と、
前記フレアの色付きを除去する除去処理を前記処理領域へ施す除去手段と、
前記フレアの光源色が非白色であった場合には、前記フレアの光源に固有の色を前記発生領域に復元させる復元処理を、前記除去処理後の前記処理領域へ施す復元手段と、
を備えたことを特徴とする画像処理装置。
【請求項９】
請求項８に記載の画像処理装置において、
前記復元処理は、
前記除去処理後の前記処理領域に対して前記除去処理前の前記処理領域を重み付け合成する処理である
ことを特徴とする画像処理装置。
【請求項１０】
請求項９に記載の画像処理装置において、
前記除去処理後の重みに対する前記除去処理前の重みの比は、
前記処理領域の位置毎に設定される
ことを特徴とする画像処理装置。
【請求項１１】
請求項１０に記載の画像処理装置において、
前記発生領域に近い位置の前記比よりも、前記発生領域から離れた位置の前記比の方が小さい
ことを特徴とする画像処理装置。
【請求項１２】
請求項８〜請求項１１の何れか一項に記載の画像処理装置において、
前記除去手段は、
前記発生領域の各色成分の中で基準となる基準色成分を設定する基準色設定手段と、
前記処理領域における基準色成分と着目色成分との間の相関を示す相関情報を算出する相関算出手段と、
前記フレアの色付き分布であるフレア色成分を、前記処理領域の前記相関情報と前記発生領域の基準色成分とに基づき推定するフレア色成分推定手段と、
を含むことを特徴とする画像処理装置。
【請求項１３】
撮影光学系が結像した被写体像を撮像するカラーの撮像素子と、
前記撮像装置が生成した画像を前記処理対象画像として処理する請求項１〜請求項１２の何れか一項に記載の画像処理装置と、
を備えたことを特徴とする撮像装置。
【請求項１４】
請求項１３に記載の撮像装置において、
回折光学素子を含む前記撮影光学系を備えた
ことを特徴とする撮像装置。
【請求項１５】
請求項１３に記載の撮像装置において、
回折光学素子を含む前記撮影光学系を装着可能である
ことを特徴とする撮像装置。
【請求項１６】
請求項１４又は請求項１５に記載の撮像装置において、
前記回折光学素子の白色光に対する応答特性は、
前記画像に発生するフレアの波長スペクトルが所定波長域にピークを持つように設定されている
ことを特徴とする撮像装置。
【請求項１７】
撮影光学系により取得されたカラーの処理対象画像を入力する入力手順と、
前記処理対象が像におけるフレアの発生領域の色に基づき前記フレアの光源色が白色であったか否かを判別する判別手順と、
前記発生領域に対してフレア補正処理を施すフレア補正手順とを含み、
前記フレア補正手順では、
前記判別手順による前記判別の結果に応じて前記フレア補正処理の内容を切り換える
ことを特徴とするコンピュータ実行可能な画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、撮影画像へフレア補正処理を施す画像処理装置、撮像装置、画像処理プログラムに関する。
【背景技術】
【０００２】
回折光学素子（ＤＯＥ：Diffractive Optical Element）を有する撮影光学系で撮影を行う際、極めて高輝度な物体（光源）が撮影画角内に存在していた場合には、撮影画像上の光源像の周辺に色滲みを伴うフレア（色フレア）が発生する。この色フレアは無色のフレアと比較すると著しく不自然である。
【０００３】
因みに、特許文献１には、フレア補正処理として、フレア発生領域の色相を背景の色相に合わせるように変換し、さらに平滑化フィルタをかけることにより、色フレアを目立たなくする電子カメラが開示されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許第４２５０５０６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、色フレアの発生源である光源の中には、固有の色を有したものも存在しており、そのような光源に起因した色フレアへ前述したフレア補正処理を施すと、色フレアの色滲みと共に光源本来の色までもが消えてしまう虞がある。
【０００６】
そこで本発明は、フレアの生源色に依らず適切なフレア補正処理を実行することのできる画像処理装置、撮像装置、及び画像処理プログラムを提供することを目的とする。
【課題を解決するための手段】
【０００７】
本発明の画像処理装置は、撮影光学系により取得されたカラーの処理対象画像を入力する入力手段と、前記処理対象画像におけるフレアの発生領域の色に基づき前記フレアの光源色が白色であったか否かを判別する判別手段と、前記発生領域に対してフレア補正処理を施すフレア補正手段とを備え、前記フレア補正手段は、前記判別手段による前記判別の結果に応じて前記フレア補正処理の内容を切り換える。
【０００８】
本発明の撮像装置は、撮影光学系が結像した被写体像を撮像するカラーの撮像素子と、前記撮像装置が生成した画像を前記処理対象画像として処理する本発明の画像処理装置とを備える。
【０００９】
本発明のコンピュータ実行可能な画像処理プログラムは、撮影光学系により取得されたカラーの処理対象画像を入力する入力手順と、前記処理対象が像におけるフレアの発生領域の色に基づき前記フレアの光源色が白色であったか否かを判別する判別手順と、前記発生領域に対してフレア補正処理を施すフレア補正手順とを含み、前記フレア補正手順では、前記判別手順による前記判別の結果に応じて前記フレア補正処理の内容を切り換える。
【発明の効果】
【００１０】
本発明によれば、フレアの光源色に依らず適切なフレア補正処理を実行することのできる画像処理装置、撮像装置、及び画像処理プログラムが実現する。
【図面の簡単な説明】
【００１１】
【図１】第１実施形態の画像処理装置の構成例を示すブロック図である。
【図２】現像処理の動作フローチャートである。
【図３】フレア補正処理の動作フローチャートである。
【図４】フレア色成分抽出処理の動作フローチャートである。
【図５】光源色判定処理の動作フローチャートである。
【図６】光源色復元処理の動作フローチャートである。
【図７】入力画像の例である。
【図８】検出された高輝度領域の例である。
【図９】フレアに起因した高輝度領域（飽和領域）の例である。
【図１０】ラベリングされた飽和領域の例である。
【図１１】飽和領域Ｓｉの半径ｒｉと、フレア領域Ａｉの半径との関係を示す模式図である。
【図１２】飽和領域Ｓｉとフレア領域Ａｉと処理領域Ａｉ’との関係を説明する図である。
【図１３】二値のマスク画像を説明する図である。
【図１４】グレースケールマスク画像を説明する図である。
【図１５】処理領域Ａｉ’のＲ成分（フレア色成分有り）とＧ成分との関係を示す図である。
【図１６】フレア色成分の抽出手順を説明する模式図である。
【図１７】処理領域Ａｉ’のＲ成分（フレア色成分無し）とＧ成分との関係を示す図である。
【図１８】（Ａ）は、彩度強調前の或る処理領域Ａｉ’の画素プロファイル（中心線上の画素プロファイル）であり、（Ｂ）は、彩度強調後の同じ処理領域Ａｉ’の画素プロファイルである。
【図１９】拡大前のグレースケールマスク画像を説明する図である。
【図２０】拡大後のグレースケールマスク画像を説明する図である。
【図２１】拡大前後のグレースケールマスク画像の画素プロファイル（中心線上の画素プロファイル）である。
【図２２】ステップＳ９９を説明する図である。
【図２３】第２実施形態の電子カメラの構成例を示すブロック図である。
【発明を実施するための形態】
【００１２】
［第１実施形態］
以下、本発明の第１実施形態を説明する。
【００１３】
図１は、本実施形態の画像処理装置の構成例を示すブロック図である。本実施形態の画像処理装置は、画像処理プログラムがインストールされたコンピュータ１１で構成される。
【００１４】
図１に示すコンピュータ１１は、データ読込部１２、記憶装置１３、ＣＰＵ１４、メモリ１５および入出力Ｉ／Ｆ１６、バス１７を有している。データ読込部１２、記憶装置１３、ＣＰＵ１４、メモリ１５および入出力Ｉ／Ｆ１６は、バス１７を介して相互に接続されている。さらに、コンピュータ１１には、入出力Ｉ／Ｆ１６を介して、入力デバイス１８（キーボード、ポインティングデバイスなど）とモニタ１９とがそれぞれ接続されている。なお、入出力Ｉ／Ｆ１６は、入力デバイス１８からの各種入力を受け付けるとともに、モニタ１９に対して表示用のデータを出力する。
【００１５】
データ読込部１２は、画像処理対象となる画像データを外部から読み込むときに用いられる。例えば、データ読込部１２は、データ読込部１２に挿入された記憶媒体から画像データを取得する読込デバイス（光ディスクドライブ、磁気ディスクドライブ、光磁気ディスクドライブなど）や、公知の通信規格に準拠して外部の装置と通信を行う通信デバイス（ＵＳＢインターフェース、ＬＡＮモジュール、無線ＬＡＮモジュールなど）で構成される。
【００１６】
記憶装置１３は、例えば、ハードディスクや、不揮発性の半導体メモリなどの記憶媒体で構成される。この記憶装置１３には、画像処理プログラムが記録されている。なお、記憶装置１３には、データ読込部１２から読み込まれた画像データを記憶しておくこともできる。
【００１７】
ＣＰＵ１４は、記憶装置１３に記憶された画像処理プログラムを実行し、コンピュータ１１の各部を統括的に制御するプロセッサである。ここで、本実施形態の画像処理プログラムには、カメラで撮影されたＲＡＷ画像に対する現像処理の機能が搭載されており、この現像処理を実行する際、ＣＰＵ１４は、各部を制御する制御部として機能する他に、フレア補正処理部２１としても機能する（フレア補正処理部２１の動作は後述する）。
【００１８】
メモリ１５は、画像処理プログラムでの各種演算結果を一時的に記憶する。このメモリ１５は、例えば揮発性のＳＤＲＡＭなどで構成される。
【００１９】
図２は、現像処理の動作フローチャートである。以下、図２の各ステップを順に説明する。
【００２０】
ステップＳ１：ＣＰＵ１４は、ユーザが現像処理の対象として指定したカラー画像（ＲＡＷ画像）の画像データを、メモリ１５上に読み込む（以下、これを「入力画像」と称す。）。この入力画像は、例えば、データ読込部１２の記憶媒体に書き込まれていた画像である。
【００２１】
ここで、入力画像には、色フレアの発生している可能性がある。色フレアの発生した画像は、ＤＯＥを含む撮影光学系で撮影された画像であって、例えば図７に示すとおりである（但し、図７はモノクロ画像なので、色フレアの色付きは現れていない。）。色フレアは通常、輪帯状又はそれに近いパターンをしているが、そのパターンは完全な輪帯状になるとは限らず、撮影条件（画像の撮影に使用された撮影光学系の結像特性、撮影時に画角内に存在していた光源の波長スペクトルなど）に依存する。また、色フレアの色づき分布（フレア色成分）も、撮影条件に依存する。
【００２２】
但し、本実施形態のフレア補正処理は、色フレアのパターン及び色付き分布が未知であっても実行可能であるので、入力画像には詳細な撮影条件データが付随していなくても構わない。よって、本実施形態では、撮影に使用された撮影光学系の種類情報が入力画像に付随していれば十分である。
【００２３】
ステップＳ２：ＣＰＵ１４は、メモリ１５上の入力画像に対して、フレア補正処理以外の現像処理を施す。ここでは簡単のため、フレア補正処理以外の現像処理は、色補間処理、彩度強調処理（色補正処理の一種）のみと仮定する。これによって、メモリ１５上の入力画像は、色補間後、かつ彩度強調後の状態となる。
【００２４】
ステップＳ３：ＣＰＵ１４は、入力画像に付随する、撮影光学系の種類情報に基づき、その撮影光学系にＤＯＥが含まれていたか否かを判別し、含まれていなかった場合には、入力画像に色フレアの発生している可能性は無いとみなしてステップＳ１０へ移行し、撮影光学系にＤＯＥが含まれていた場合には、入力画像に色フレアの発生している可能性があるとみなしてステップＳ４へ移行する。
【００２５】
ステップＳ４：ＣＰＵ１４は、入力画像のサイズ縮小版（縮小画像）を作成する。なお、縮小画像を作成するのは、後続するステップの演算負荷を軽減するためである。入力画像に対する縮小画像のサイズ比は、１／４、或いは、１／８などである。因みに、１／８よりもサイズ比を小さくすると、小さな光源を検出できなくなる虞があるので、サイズ比は１／８以上であることが望ましい。
【００２６】
ステップＳ５：ＣＰＵ１４は、縮小画像の各画素の輝度を所定値（ここでは飽和レベルとする。）と比較することにより、輝度が所定値を超過するような明るい領域（高輝度領域）を縮小画像から検出する（図８参照）。図８に示すとおり、検出された高輝度領域の中には、フレアに起因した高輝度領域だけでなく、フレアとは無関係な高輝度領域も混在している可能性がある。
【００２７】
なお、本ステップでは演算負荷を軽減するため、各画素の輝度の指標として、各画素の輝度値（Ｒ、Ｇ、Ｂの重み付け和）を使用する代わりに、各画素の各色成分の和（Ｒ、Ｇ、Ｂの和）を使用してもよい。或いは、各画素の輝度の指標として、各画素を代表する色成分の値（Ｇ強度）を使用してもよい。また、このような演算負荷の軽減方法は、画素の輝度を扱う他のステップでも採用することができる。
【００２８】
ステップＳ６：ＣＰＵ１４は、ステップＳ５で検出された高輝度領域の各々がフレアに起因したものであるか否かを判別し、フレアに起因したもの以外を排除する（図９参照。）。
【００２９】
なお、個々の高輝度領域がフレアに起因したものであるか否かの判定は、例えば次のとおり行われる。すなわち、ＣＰＵ１４は、高輝度領域と予め用意した円形パターンとを比較し、両者の相関度が一定以上である場合には、その高輝度領域をフレアに起因したものとみなし、相関度が一定未満である場合には、その高輝度領域をフレアに起因したものではないとみなす。なお、比較対象となる円形パターンとしては、フレアに起因した高輝度領域として平均的なものを使用すればよい。
【００３０】
或いは、個々の高輝度領域がフレアに起因したものであるか否かの判定は、次のとおり行われてもよい。すなわち、ＣＰＵ１４は、高輝度領域の中心から所定半径を有した円形領域を設定し、その円形領域の全域に高輝度領域が存在していた場合には、その高輝度領域をフレアに起因したものとみなし、円形領域の少なくとも１部が高輝度領域から外れていた場合には、その高輝度領域をフレアに起因したものではないとみなす。なお、円形領域の半径は、フレアに起因した高輝度領域として最も小さいものの半径と同程度に設定される。
【００３１】
したがって、本ステップで除外されなかった高輝度領域は、フレアの中心部における飽和領域に相当する。以下、本ステップで除外されなかった高輝度領域を、単に「飽和領域」と称す。
【００３２】
ステップＳ７：ＣＰＵ１４は、ステップＳ６で検出された１又は複数の飽和領域に対してラベリングを行い、通し番号（フレア番号ｉ＝１、２、…）を付与する（図１０参照）。これによって、入力画像に発生しているフレアの総数（すなわち、フレア番号ｉの最終値）が確定する。
【００３３】
ステップＳ８：ＣＰＵ１４は、フレアの総数がゼロである場合には、フレア補正処理は不要であるとみなしてステップＳ１０へ移行し、フレアの総数がゼロでない場合には、フレア補正処理を実行するべくステップＳ９へ移行する。
【００３４】
ステップＳ９：ＣＰＵ１４は、メモリ１５上の入力画像に対してフレア補正処理を施す。フレア補正処理の詳細は、後述する。
【００３５】
ステップＳ１０：ＣＰＵ１４は、メモリ１５上の入力画像を、ＪＰＥＧなどの所定の圧縮方式で圧縮してから、ユーザが指定した保存先（例えば、データ読込部１２に挿入された記憶媒体）へ保存し、現像処理のフローを終了する。
【００３６】
図３は、フレア補正処理の動作フローチャートである。以下、図３の各ステップを順に説明する。
【００３７】
ステップＳ９１：フレア補正処理部２１は、フレア番号ｉを初期値（１）に設定する。
【００３８】
ステップＳ９２：フレア補正処理部２１は、縮小画像上でフレア番号ｉに対応する飽和領域Ｓｉに着目すると、その飽和領域Ｓｉを包含する処理領域Ａｉ’を縮小画像上に設定する。この処理領域Ａｉ’は、飽和領域Ｓｉを中心としたフレアの発生領域（フレア領域）Ａｉを包含している必要がある。
【００３９】
ここで、図１１に示すとおり、飽和領域Ｓｉが大きいときほど、フレア領域Ａｉも大きくなる。但し、飽和領域Ｓｉの半径ｒｉが閾値ａ以下である場合（図１１（ａ）、（ｂ）、（ｃ））は、フレア領域Ａｉの半径は所定値βとなるが、飽和領域の半径ｒｉが閾値ａを超過する場合（図１１（ｄ）、（ｅ））は、フレア領域Ａｉの半径は半径ｒｉに応じた値α×ｒｉとなる。なぜなら、飽和領域Ｓｉの半径ｒｉが閾値ａ以下である場合は、光源の明るさが一定未満であり、低次の回折光しか発生しないのに対して、飽和領域Ｓｉの半径ｒｉが閾値ａを超過する場合は、光源の明るさが一定以上であり、低次の回折光に加えて高次の回折光も発生するからである。
【００４０】
そこで、本ステップにおけるフレア補正処理部２１は、処理領域Ａｉの設定を以下の手順（Ａ）〜（Ｅ）で行う。
【００４１】
（Ａ）フレア補正処理部２１は、縮小画像上で飽和領域Ｓｉの輝度重心に相当する座標ｇｉを算出し、座標ｇｉから飽和領域Ｓｉの輪郭までの平均距離を、飽和領域Ｓｉの半径ｒｉとして算出する。
【００４２】
（Ｂ）フレア補正処理部２１は、飽和領域Ｓｉの半径ｒｉが閾値ａ以下であるか否かを判別する。閾値ａ以下である場合は手順（Ｃ）へ移行し、そうでない場合は手順（Ｄ）へ移行する。
【００４３】
（Ｃ）フレア補正処理部２１は、図１２（ａ）に示すとおり、縮小画像上で座標ｇｉを中心とし、かつ半径が所定値βである円形領域を、フレア番号ｉに対応したフレア領域Ａｉに設定し（但しβ＞ｒｉ）、手順（Ｅ）へ移行する。
【００４４】
（Ｄ）フレア補正処理部２１は、図１２（ｂ）に示すとおり、縮小画像上で座標ｇｉを中心とし、かつ半径が（α×ｒｉ）である円形領域を、フレア番号ｉに対応したフレア領域Ａｉに設定する（但しα＞１）。
【００４５】
（Ｅ）フレア補正処理部２１は、図１２（ａ）又は図１２（ｂ）に示すとおり、縮小画像上でフレア領域Ａｉの全体を包含し、かつ一定のマージンγを持った包含領域を、フレア番号ｉに対応した処理領域Ａｉ’に設定する（但し、γ＞０）。
【００４６】
なお、本ステップでは、処理領域Ａｉ’のサイズをフレア領域Ａｉのサイズに応じて設定したが、処理領域Ａｉ’がフレア領域Ａｉを確実に包含できるのであれば、処理領域Ａｉ’のサイズを予め決められたサイズ（十分に大きなサイズ）としてもよい。また、処理領域Ａｉ’の形状は、円形でも矩形でも構わない。また、処理領域Ａｉ’の中心は、フレア領域Ａｉの中心から外れていても構わない。
【００４７】
また、本ステップでは、飽和領域Ｓｉの輝度重心を座標ｇｉとしたが、演算負荷を軽減するために、飽和領域Ｓｉの中心を座標ｇｉとしてもよい。
【００４８】
また、本ステップでは、飽和領域Ｓｉのサイズとして飽和領域Ｓｉの半径ｒｉを算出したが、飽和領域Ｓｉのサイズとして飽和領域Ｓｉの画素数を使用してもよい。その場合は、半径の閾値ａの代わりに画素数の閾値ａ’が使用される。
【００４９】
また、飽和領域Ｓｉのサイズとフレア領域Ａｉのサイズとの関係（図１１参照）は、ＤＯＥの種類（すなわち撮影光学系の種類）によって異なるため、本ステップのフレア補正処理部２１は、サイズの閾値を、撮影光学系の種類に応じて設定することが望ましい。
【００５０】
ステップＳ９３：フレア補正処理部２１は、フレア番号ｉに対応する処理領域Ａｉ’に対してフレア色成分抽出処理を施し、処理領域Ａｉ’に発生しているフレア色成分を抽出する。なお、本ステップで抽出されるフレア色成分は、２色のフレア色成分からなる。フレア色成分抽出処理の詳細は、後述する。
【００５１】
ステップＳ９４：フレア補正処理部２１は、ステップＳ９３で抽出されたフレア色成分（２色のフレア色成分）へサイズ拡大処理を施す。なお、前述したステップＳ４におけるサイズ比を１／Ｍとした場合、本ステップにおけるサイズ拡大処理ではフレア色成分（２色分のフレア色成分）のサイズをＭ倍にする。これによって、入力画像に対応するサイズのフレア色成分（２色のフレア色成分）が取得される。
【００５２】
ステップＳ９５：フレア補正処理部２１は、サイズ拡大後のフレア色成分（２色のフレア色成分）を、入力画像上で処理領域Ａｉ’に対応する領域から減算する。但し、減算に当たりＣＰＵ１４は、縮小画像上の座標ｇｉと前述した倍率Ｍとに基づき、座標ｇｉに対応する入力画像上の座標ｇｉ’を求め、その座標ｇｉ’に対してフレア色成分（２色のフレア色成分）の中心を一致させておく。これによって、入力画像上で処理領域Ａｉ’に対応する領域に発生していた色フレアのフレア色成分が除去される。つまり、その領域に発生していた色フレアは、無色のフレアとなる。
【００５３】
ステップＳ９６：フレア補正処理部２１は、縮小画像上のフレア領域Ａｉに対して光源色判定処理を施し、フレア領域Ａｉに発生しているフレアの光源色が白色であるか否かを判定する。光源色判定処理の詳細は、後述する。なお、ここでいう「白色」は、ヒトの目から見て白色に見える色のことであって、所謂白色光源の色（波長スペクトルがフラット）だけでなく、それに近い色をも含むものとする。
【００５４】
ステップＳ９７：フレア補正処理部２１は、ステップＳ９６で判定した光源色が白色であったか否かを判別し、非白色であった場合は、光源色復元処理を実行するべくステップＳ９８へ移行し、白色であった場合は、光源色復元処理をスキップするべくステップＳ９１１へ移行する。
【００５５】
ステップＳ９８：フレア補正処理部２１は、縮小画像上の処理領域Ａｉ’に関する光源色復元処理を実行することにより、復元後の処理領域Ａｉ’（＝復元画像ｉｍｇＣ）を取得する。光源色復元処理の詳細は、後述する。
【００５６】
ステップＳ９９：フレア補正処理部２１は、復元画像ｉｍｇＣのサイズ拡大版（＝復元画像ＩＭＧＣ）を作成する。復元画像ＩＭＧＣの作成方法の詳細は、後述する。
【００５７】
ステップＳ９１０：フレア補正処理部２１は、入力画像上で処理領域Ａｉ’に対応する領域を復元画像ＩＭＧＣに置換することにより、その領域における光源の色味を復元する。
【００５８】
ステップＳ９１１：フレア補正処理部２１は、フレア番号ｉが最終値に達したか否かを判別し、達していなかった場合はステップＳ９１２へ移行し、達していた場合はフレア補正処理のフローを終了する。
【００５９】
ステップＳ９１２：フレア補正処理部２１は、フレア番号ｉをインクリメントしてからステップＳ９２へ移行する。
【００６０】
図４は、フレア色成分抽出処理のフローチャートである。以下、各ステップを順に説明する。
【００６１】
ステップＳ９３１：フレア補正処理部２１は、フレア番号ｉに対応するフレア領域Ａｉの平均値を色成分毎に（ＲＧＢ別に）算出し、平均値が最小であった１つの色成分を、基準色成分に定める。以下、Ｇ成分が基準色成分に設定されたと仮定する。その場合、フレア色成分の推定は、Ｇ成分以外の色成分、すなわちＲ成分及びＢ成分の各々に関して行われる。
【００６２】
ステップＳ９３２：フレア補正処理部２１は、フレア番号ｉに対応した二値のマスク画像を作成する。二値のマスク画像は、図１３に示すとおり、処理領域Ａｉ’と同形同サイズのマスク画像であって、フレア領域Ａｉに対応する領域に開口部を設けたものである。さらに、フレア補正処理部２１は、その二値のマスク画像を、ガウシアンフィルタなどの平滑化フィルタで空間的に平滑化することにより、図１４に示すようなグレースケールマスク画像を作成する。このグレースケールマスク画像は、二値のマスク画像において、開口部と非開口部との境界のコントラストを意図的に低下させたものに相当する。このグレースケールマスク画像の単位面積当たりの平均値は、処理領域Ａｉ’の中央から周辺に向かって低下している。
【００６３】
なお、本ステップで使用される平滑化フィルタのサイズ（フィルタ径）は、フレア領域Ａｉのサイズが大きいときほど大きく設定されることが望ましい。
【００６４】
ステップＳ９３３：フレア補正処理部２１は、基準色成分（ここではＧ成分）以外の１つの色成分を着目色成分に設定する。以下、Ｒ成分が着目色成分に設定されたと仮定する。図１５は、或る処理領域Ａｉ’のＲ成分と、同じ処理領域Ａｉ’のＧ成分との関係を示す図（画素プロファイル）である。図１５に示すとおり、Ｒ成分においては、フレア領域Ａｉの内側かつ飽和領域Ｓｉの外側の領域に、フレア色成分が重畳されている。つまり、その領域のＲ成分は、同じ領域のＧ成分と比較して著しく高い値となっている。
【００６５】
本ステップのフレア補正処理部２１は、処理領域Ａｉ’のＲ成分を処理領域Ａｉ’のＧ成分で画素毎に除算することにより、Ｇ成分を基準としたＲ成分の相関色分布画像ＩＧＲを作成する。図１６（Ａ）の左側に示すのは、相関色分布画像ＩＧＲの模式図である。
【００６６】
ここで、相関色分布画像ＩＧＲの中央領域ａの多くの画素は、１に近い画素値を有しており、中央領域ａの周辺に位置するリング状領域ｂの多くの画素は、１よりも高い画素値を有しており、リング状領域ｂの外側に位置する最周辺領域ｃの多くの画素は、リング状領域ｂよりも低い画素値を有している。このうち、中央領域ａは、図１５に示した飽和領域Ｓｉに対応しており、リング状領域ｂは、図１５に示したフレア領域Ａｉの内側かつ飽和領域Ｓｉの外側の領域、すなわち、フレア色成分の重畳された領域に対応している。
【００６７】
ステップＳ９３４：フレア補正処理部２１は、相関色分布画像ＩＧＲのうち、周辺と比較して値の低い領域（中央領域ａ及び最周辺領域ｃ）のみを平滑化することにより、図１６（Ｂ）に示すような相関色分布画像ＩＧＲ”を取得する。なお、相関色分布画像ＩＧＲ”の取得は、図１６（Ａ）→図１６（Ｂ）に示す手順によって行うことができる。
【００６８】
すなわち、フレア補正処理部２１は、先ず、図１６（Ａ）の左側に示した相関色分布画像ＩＧＲの全体を、ガウシアンフィルタなどの平滑化フィルタで空間的に平滑化することにより、図１６（Ａ）の右側に示すような相関色分布画像ＩＧＲ’を取得する。この平滑化によると、周辺と比較して値の低い領域の値は、平滑化前よりも高めになり、周辺と比較して値の高い領域の値は、平滑化前よりも低めになる。続いて、フレア補正処理部２１は、平滑化前後の２つの相関色分布画像ＩＧＲ、ＩＧＲ’を画素毎に比較し、値の高かった方の画素を集めて１枚の合成画像を作成し、その合成画像を相関色分布画像ＩＧＲ”とする（図１６（Ｂ））。この相関色分布画像ＩＧＲ”は、元の相関色分布画像ＩＧＲにおいて、フレア色成分の重畳されていなかった領域（点線枠部分）のみを平滑化したものに相当する。
【００６９】
ステップＳ１８３５：フレア補正処理部２１は、相関色分布画像ＩＧＲ”に対して最小値フィルタ等による収縮処理を施すことにより、相関色分布画像ＩＧＲ”において周辺と比較して値の高い領域の面積を収縮し、図１６（Ｃ）に示すような相関色分布画像ＩＨを取得する。さらに、フレア補正処理部２１は、相関色分布画像ＩＨをガウシアンフィルタなどの平滑化フィルタで空間的に平滑化することにより、図１６（Ｄ）に示すような相関色分布画像ＩＨ’を取得する。この相関色分布画像ＩＨ’は、元の相関色分布画像ＩＧＲにおいてフレア色成分の重畳されていた領域を周辺の領域（点線枠部分）によって補間したものに相当する。
【００７０】
ステップＳ９３６：フレア補正処理部２１は、処理領域Ａｉ’の基準色成分（ここではＧ成分）に対して相関色分布画像ＩＨ’を画素毎に乗算することにより、フレア色成分の除去されたＲ成分ＩＲｎｅｗを算出する。これによって、処理領域Ａｉ’からＲ成分に関するフレア色成分が除去されたことになる。なお、図１７は、除去後の処理領域Ａｉ’の画素プロファイルである。
【００７１】
この図１７を図１５と比較すれば明らかなとおり、除去後には、フレア色成分の重畳されていた領域（フレア領域Ａｉの内側かつ飽和領域Ｓｉの外側の領域）のＲ強度は、同じ領域のＧ強度に近づけられているのがわかる。また、除去後には、フレア色成分の重畳されていた領域のＲ／Ｇ強度は、フレア領域Ａｉの外側の領域のＲ／Ｇ強度と、飽和領域Ｓｉの内側のＲ／Ｇ強度とによって補間されているのがわかる。また、除去後には、フレア色成分の重畳されていた領域のＲ強度分布形状は、同じ領域のＧ強度分布形状に対応した形状となっているのがわかる。
【００７２】
なお、本ステップのフレア補正処理部２１は、フレア色成分の除去されたＲ成分ＩＲｎｅｗを、元のＲ成分から差し引くことにより、Ｒ成分に関するフレア色成分ＩＦＲを算出しておく（フレア色成分の除去されたＲ成分ＩＲｎｅｗを、フレア色成分ＩＦＲに換算しておく。）。なぜなら、本実施形態では、最終画像を算出する前に、フレア色成分ＩＦＲに対して幾つかの処理（マスク処理やサイズ拡大処理など）を施す必要があるからである。
【００７３】
ステップＳ９３７：フレア補正処理部２１は、ステップＳ９３６で取得したフレア色成分ＩＦＲに対してステップＳ９３２で作成したグレースケールマスクを適用することにより、フレア色成分ＩＦＲの周辺部をぼかす。このように周辺部をぼかしておけば、処理領域Ａｉ’のＲ成分からフレア色成分ＩＦＲを除去（減算）した際に、処理領域Ａｉ’に不自然な階調段差が発生することを防げるからである。また、上述したとおり、グレースケールマスクの作成時に使用された平滑化フィルタのサイズは、フレア領域Ａｉのサイズに応じたサイズに設定されたので、マスク適用後のフレア色成分ＩＦＲの周辺部は、より自然なものとなる。これによって、Ｒ成分に関するフレア色成分ＩＦＲの推定は、完了である。
【００７４】
ステップＳ９３８：フレア補正処理部２１は、基準色成分以外の全ての色成分（ここではＲ成分、Ｂ成分）に関する推定が完了したか否かを判別し、完了していなければ未推定である色成分の推定を行うべくステップＳ９３３へ移行し、完了していればフレア色成分抽出処理のフローを終了する。なお、ここでは、Ｒ成分に関する推定しか説明しなかったが、Ｂ成分に関する推定は、Ｒ成分に関する推定と同様なので、説明を省略する。
【００７５】
図５は、光源色判定処理のフローチャートである。以下、各ステップを順に説明する。
【００７６】
ステップＳ９６１：フレア補正処理部２１は、縮小画像上のフレア領域Ａｉに着目し、そのフレア領域Ａｉから飽和領域Ｓｉを除外してなるリング状領域（図１１の斜線部）を、参照領域Ｒｉに設定する。
【００７７】
ステップＳ９６２：フレア補正処理部２１は、参照領域Ｒｉの平均色が無彩色（又は無彩色相当）であるか否かを判別し、無彩色（又は無彩色相当）でない場合には、ステップＳ９６３へ移行し、無彩色（又は無彩色相当）である場合には、ステップＳ９６６へ移行する。
【００７８】
ステップＳ９６３：フレア補正処理部２１は、参照領域Ｒｉの平均色が撮影光学系の固有フレア色であるか否かを、所定の閾範囲内に収まるか否かによって判別し、固有フレア色である場合にはステップＳ９６４へ移行し、固有フレア色でない場合にはステップＳ９６３へ移行する。以下、固有フレア色を詳しく説明する。
【００７９】
撮影光学系の中には、赤色のフレアを発生させ易いものや、青色のフレアを発生させ易いものや、紫色のフレアを発生させ易いものなどが存在するので、ここでは、或る撮影光学系が最も高頻度に発生させるフレアの色を、その撮影光学系の「固有フレア色」と称している。この固有フレア色は、撮影光学系に含まれるＤＯＥの種類（つまり撮影光学系の種類）によって決まる。なぜなら、ＤＯＥの白色光に対する応答特性は、撮影画像上に発生するフレアの波長スペクトルが予め決められた波長域にピーク（高強度）を持つように設定されている。或いは、ＤＯＥの白色光に対する応答特性は、撮影画像上に発生するフレアの波長スペクトルが予め決められた波長域にピーク（高強度）を持ち易いように設定されている。
【００８０】
なお、固有フレア色は撮影光学系の種類によって異なるため、本ステップのフレア補正処理部２１は、固有フレア色の閾値範囲を撮影光学系の種類に応じて設定する。
【００８１】
ステップＳ９６４：フレア補正処理部２１は、参照領域Ｒｉの平均色の彩度が閾値以上であるか否かを判別し、閾値以上である場合は、ステップＳ９６５へ移行し、閾値未満である場合は、ステップＳ９６６へ移行する。なお、本ステップのフレア補正処理部２１は、本ステップの彩度の閾値をステップＳ２における彩度強調処理の強度に応じて設定する。
【００８２】
なお、本ステップにおけるフレア補正処理部２１は、参照領域Ｒｉの平均色の彩度が閾値以上であるか否かを判別する代わりに、色補正済みの画素（色補正画素）が参照領域Ｒｉの中に存在するか否かを判別してもよい。以下、色補正画素を説明する。
【００８３】
図１８（Ａ）は、彩度強調前の或る処理領域Ａｉ’の画素プロファイル（中心線上の画素プロファイル）であり、図１８（Ｂ）は、彩度強調後の同じ処理領域Ａｉ’の画素プロファイルである。この処理領域Ａｉ’に施された彩度強調は、Ｇ成分及びＢ成分をＲ成分に比して相対的に強調するための彩度強調である。
【００８４】
図１８（Ａ）に示すとおり、彩度強調前の処理領域Ａｉ’には、全色成分がノイズレベルＴ１以上となった画素しか存在していないのに対し、図１８（Ｂ）に示すとおり、彩度強調後の処理領域Ａｉ’には、Ｒ成分がノイズレベルＴ１未満となった画素が発生している。しかも、そのような画素の少なくとも一部の画素は、Ｇ成分とＢ成分との少なくとも一方が飽和レベルＴ２に達している。このように、少なくとも１つの色成分が飽和レベルＴ２に達しており、かつ他の少なくとも１つの色成分がノイズレベルＴ１未満となっているような画素が、色補正画素である。このような色補正画素が発生するのは、参照領域Ｒｉである。
【００８５】
ステップＳ９６５：フレア補正処理部２１は、処理領域Ａｉ’に発生しているフレアの光源色は非白色であったと判定し、光源色判定処理のフローを終了する。
【００８６】
ステップＳ９６６：フレア補正処理部２１は、処理領域Ａｉ’に発生しているフレアの光源色は白色であったと判定し、光源色判定処理のフローを終了する。
【００８７】
図６は、光源色復元処理のフローチャートである。以下、各ステップを順に説明する。
【００８８】
ステップＳ９８１：フレア補正処理部２１は、ステップＳ９３２で作成したものと同様のグレースケールマスク画像（図１９参照）を取得する。
【００８９】
次に、フレア補正処理部２１は、このグレースケールマスク画像に対して正の所定値を乗算することにより、グレースケールマスク画像の画素値を向上させる。但し、フレア補正処理部２１は、乗算によって１を超過した画素値については１に置換することで、グレースケールマスク画像の画素値を０から１の範囲内に収める。その結果、グレースケールマスク画像の開口部の面積が拡大される（図２０参照）。拡大前におけるグレースケールマスク画像のプロファイル（開口率分布）は、図２１の点線で示すとおりであって、拡大後におけるグレースケールマスク画像のプロファイル（開口率分布）は、図２１の実線で示すとおりである。図２１に示すとおり、拡大後におけるグレースケールマスク画像の開口部は、フレア領域Ａｉとほぼ一致している。以下、拡大後におけるグレースケールマスク画像を、「グレースケールマスク画像ｉｍｇＦ」とおく。このグレースケールマスク画像ｉｍｇＦは、後述する重み付け合成（ステップＳ９８４）の重み分布として使用される。
【００９０】
ステップＳ９８２：フレア補正処理部２１は、フレア除去前の縮小画像ｉｍｇＡを用意する。フレア除去前の縮小画像ｉｍｇＡとは、縮小画像上の処理領域Ａｉ’のことである。
【００９１】
ステップＳ９８３：フレア補正処理部２１は、フレア除去前の縮小画像ｉｍｇＡからサイズ拡大前のフレア色成分を減算することにより、フレア除去後の縮小画像ｉｍｇＢを用意する（図２２（ａ））。サイズ拡大前のフレア色成分とは、ステップＳ９３において取得されたフレア色成分のことである。
【００９２】
ステップＳ９８４：フレア補正処理部２１は、フレア除去前の縮小画像ｉｍｇＡと、フレア除去後の縮小画像ｉｍｇＢとを重み付け合成することにより、光源の色味の復元された復元画像ｉｍｇＣを取得すると（図２２（ｂ））、光源色復元処理のフローを終了する。
【００９３】
ここで、本ステップにおける重み付け合成の重みは、画像上の位置毎（画素毎）に設定される。そして、縮小画像ｉｍｇＡの重み分布は、ステップＳ９８１で作成されたグレースケールマスク画像ｉｍｇＦの画素値分布と同じに設定される。つまり、本ステップの重み付け合成は、以下の式によって行われる。
【００９４】
【数１】
2014165833.tif 000003
【００９５】
なお、式（１）における（ｘ，ｙ）は、画素座標であり、ｒ，ｇ，ｂは、画素値のＲ成分、Ｇ成分、Ｂ成分である。
【００９６】
したがって、本ステップで取得される復元画像ｉｍｇＣにおいて、フレア領域に対応する中央領域は、フレア除去前の縮小画像Ａと同じ値を採り、フレア領域外に対応する周辺領域は、フレア除去後の縮小画像Ｂに近い値を採る。しかも、その周辺領域では、中央領域から離れるほど、フレア除去後の縮小画像Ｂが強く反映されている。
【００９７】
すなわち、復元画像ｉｍｇＣの中央領域では、光源の色味が復元され、周辺領域では、フレアの色付きが除去されている。しかも、その周辺領域では、中央部領域から離れるほど除去の程度が高くなっている。
【００９８】
但し、この復元画像ｉｍｇＣのサイズは小さいので、続くステップ９９では、復元画像ｉｍｇＣのサイズ拡大版（＝復元画像ＩＭＧＣ）を作成する必要がある。
【００９９】
以下、ステップＳ９９を詳細に説明する。ステップＳ９９は、以下の手順（ｉ）〜（ｉｖ）からなる。
【０１００】
（ｉ）フレア補正処理部２１は、図２２（ｃ）に示すとおり、フレア除去前の縮小画像ｉｍｇＡから、復元画像ｉｍｇＣを減算することにより、修正フレア色成分ｉｍｇＤを求める。この修正フレア色成分ｉｍｇＤは、真のフレア色成分を表す。すなわち、この修正フレア色成分ｉｍｇＤは、フレア色成分から光源の色味成分を除去したものを表す。
【０１０１】
（ｉｉ）フレア補正処理部２１は、図２２（ｄ）に示すとおり、修正フレア色成分ｉｍｇＤに対してサイズ拡大処理を施すことにより、サイズ拡大後の修正フレア色成分ＩＭＧＤを作成する。なお、前述したステップＳ４におけるサイズ比を１／Ｍとした場合、本ステップにおけるサイズ拡大処理では修正フレア色成分ｉｍｇＤのサイズをＭ倍にする。
【０１０２】
（ｉｉｉ）フレア補正処理部２１は、図２２（ｅ）に示すとおり、フレア除去前の非縮小画像ＩＭＧＡを用意する。フレア除去前の非縮小画像ＩＭＧＡは、入力画像上で処理領域Ａｉ’に対応する領域と同じものである。
【０１０３】
（ｉｖ）フレア補正処理部２１は、図２２（ｆ）に示すとおり、フレア除去前の非縮小画像ＩＭＧＡから、サイズ拡大後の修正フレア色成分ＩＭＧＤを減算することにより、サイズ拡大後の復元画像ＩＭＧＣを取得する（以上、ステップＳ９９の説明。）。
【０１０４】
［第１実施形態の効果］
以上、本実施形態の画像処理プログラムは、撮影光学系により取得されたカラーの処理対象画像を入力する入力手順（Ｓ１）と、処理対象画像におけるフレアの発生領域（Ａｉ）の色に基づきフレアの光源色が白色であったか否かを判別する判別手順（Ｓ９６、Ｓ９７）と、発生領域に対してフレア補正処理を施すフレア補正手順（Ｓ９３、Ｓ９４、Ｓ９５、Ｓ９８、Ｓ９９、Ｓ９１０）とを含み、フレア補正手順では、判別手順による判別の結果に応じてフレア補正処理の内容を切り換える。
【０１０５】
したがって、本実施形態の画像処理プログラムによれば、光源色が白色であるか否かに拘らず適切なフレア補正処理を実行することができる。
【０１０６】
また、本実施形態の判別手順では、撮影光学系の種類、発生領域の彩度の少なくとも一方に基づき判別を行う（Ｓ９６３、Ｓ９６４）。また、本実施形態の判別手順では、処理対象画像が色補正後の画像である場合には、発生領域の彩度に基づく代わりに、発生領域に色補正画素が存在するか否かに基づき判別を行う（Ｓ９６４）。よって、判別の精度は高い。
【０１０７】
また、本実施形態の判別手順では、処理対象画像上の飽和画素領域（Ｓｉ）を包含する円形領域を発生領域（Ａｉ）とみなす（Ｓ９２）。また、本実施形態の判別手段では、飽和画素領域のサイズに基づき発生領域のサイズを推定する（Ｓ９２）。よって、フレアの発生領域は、高精度に検知される。
【０１０８】
また、本実施形態の判別手順では、発生領域から飽和画素領域を除外した領域（Ｒｉ）のみを判別に使用する（Ｓ９６１）。よって、判別の効率は高い。
【０１０９】
また、本実施形態のフレア補正手順は、発生領域を包含する処理領域（Ａｉ’）を処理対象画像上に設定する領域設定手順（Ｓ９２）と、フレアの色付きを除去する除去処理を処理領域へ施す除去手順（Ｓ９３、Ｓ９４、Ｓ９５）と、フレアの光源色が非白色であった場合には、フレアの光源に固有の色を発生領域に復元させる復元処理を、除去処理後の処理領域へ施す復元手順（Ｓ９８、Ｓ９９、Ｓ９１０）とを含む。よって、光源色が非白色であった場合には、フレアの光源の色味が復元される。
【０１１０】
また、本実施形態の復元処理は、除去処理後の処理領域に対して除去処理前の処理領域を重み付け合成する処理である（Ｓ９８１、Ｓ９８２、Ｓ９８３、Ｓ９８４）。よって、光源固有の色を簡単に復元することができる。
【０１１１】
また、本実施形態の画像処理プログラムでは、除去処理後の重みに対する除去処理後の重みの比は、処理領域の位置毎に設定される。よって、復元をきめ細かに行うことができる。
【０１１２】
また、本実施形態の画像処理プログラムでは、発生領域に近い位置の比よりも、発生領域から離れた位置の比の方が小さい。よって、復元後の処理領域に不自然な階調段差が現れるのを防ぐことができる。
【０１１３】
また、本実施形態の除去手順は、発生領域の各色成分の中で基準となる基準色成分を設定する基準色設定手順（Ｓ９３１）と、処理領域における基準色成分と着目色成分との間の相関を示す相関情報を算出する相関算出手順（Ｓ９３３）と、フレアの色付き分布であるフレア色成分を、処理領域の相関情報と発生領域の基準色成分とに基づき推定するフレア色成分推定手順（Ｓ９３６）とを含む。よって、本実施形態の除去手順は、フレアの色付きを確実に除去するものの、フレアの光源の色味も一緒に除去する可能性がある。したがって、本実施形態の除去手順に本実施形態の復元手順を組み合わせることの効果は高い。
【０１１４】
［第１実施形態の補足］
なお、本実施形態では、フレア補正処理以外の現像処理として、色補間処理及び彩度強調処理を実行したが、彩度強調処理については省略しても構わない。但し、彩度強調処理を省略した場合は、ステップＳ９６４において、参照領域Ｒｉの彩度の高さを判別する代わりに色補正済みの画像（色補正画素）の有無を判別する、という代替手法を採用することはできない。
【０１１５】
また、本実施形態では、現像処理前の画像（ＲＡＷ画像）を入力画像としたが、現像処理後の画像（ＪＰＥＧ画像等）を入力画像としても構わない。この場合、彩度強調処理等の色補正処理が入力画像に対して実行済みであったならば、ステップＳ９６４において、参照領域Ｒｉの彩度の高さを判別する代わりに色補正済みの画素（色補正画素）の有無を判別する、という代替手法を採用することができる。
【０１１６】
また、上述したステップＳ９２の手順（Ｄ）では、飽和領域Ｓｉの半径ｒｉに所定値αを乗算したものをフレア領域Ａｉの半径とみなしたが、飽和領域Ｓｉの半径ｒｉに所定値σを加えたものをフレア領域Ａｉの半径とみなしてもよい。
【０１１７】
また、本実施形態のフレア補正処理部２１は、フレア色成分の推定を基準色成分以外の２つの色成分の各々について行ったが、２つの色成分のうち一方の色成分に関する推定を省略してもよい。Ｇ成分が基準色成分であると仮定したならば、例えば、フレア領域ＡｉにおけるＧ成分の平均値とＢ成分の平均値との差が所定範囲内であった場合、つまり、Ｇ成分とＢ成分とに差異が無かった場合には、Ｂ成分に関する推定を省略してもよい。
【０１１８】
また、本実施形態のフレア補正処理部２１は、基準色成分を自動的に設定したが、ユーザに指定させてもよい。
【０１１９】
また、本実施形態の現像処理では、演算負荷を軽減するために、入力画像をそのまま使用する代わりに入力画像のサイズ縮小版（縮小画像）を使用したが、演算負荷の軽減よりも精度向上を優先させたい場合には、入力画像をそのまま使用してもよい。
【０１２０】
また、本実施形態の現像処理では、入力画像におけるフレアの発生位置をＣＰＵ１４が自動的に検出したが、ユーザに指定させてもよい。
【０１２１】
また、本実施形態の現像処理では、フレア領域ＡｉをＣＰＵ１４が自動的に検出したが、ユーザに指定させてもよい。
【０１２２】
また、本実施形態の現像処理では、処理領域Ａｉ’をＣＰＵ１４が自動的に設定したが、ユーザに設定させてもよい。
【０１２３】
また、本実施形態のフレア補正処理部２１は、フレア領域Ａｉのフレア色成分を、処理領域Ａｉにおける基準色成分と着目色成分との間の相関と、フレア領域Ａｉの基準色成分とに基づき推定したが、フレア領域Ａｉの色分布に基づく他の手法により推定しても構わない。
【０１２４】
［第２実施形態］
以下、本実施形態の第２実施形態を説明する。
【０１２５】
図２３は、電子カメラの概略構成を示すブロック図である。図２３に示すとおり、電子カメラ１１１は、撮像光学系１１２と、レンズ駆動部１１３と、絞り１１４と、絞り駆動部１１５と、カラー撮像素子１１６と、ＡＦＥ１１７と、画像処理エンジン１１８と、第１メモリ１１９と、第２メモリ１２０と、メディアＩ／Ｆ１２１と、通信Ｉ／Ｆ１２２と、モニタ１２３と、レリーズ釦１２４とを有しており、この中でレンズ駆動部１１３、絞り駆動部１１５、ＡＦＥ１１７、第１メモリ１１９、第２メモリ１２０、メディアＩ／Ｆ１２１、通信Ｉ／Ｆ１２２、モニタ１２３、レリーズ釦１２４の各々は、画像処理エンジン１１８に接続されている。また、撮影光学系１１２は、例えば、ＤＯＥレンズを含んだ高機能な撮影光学系である。
【０１２６】
第１メモリ１１９は、揮発性の記憶媒体（ＳＤＲＡＭなど）で構成されており、画像処理エンジン１１８による画像処理の前工程や後工程で撮像画像を一時的に記憶する。一方、第２メモリ１２０は、フラッシュメモリ等の不揮発性の記憶媒体で構成されており、画像処理エンジン１１８によって実行されるプログラムを長期的に記憶する。
【０１２７】
この第１メモリ１１９には、画像処理エンジン１１８が実行すべき画像処理プログラムが格納されており、この画像処理プログラムには、第１実施形態と同様の現像処理の機能が搭載されている。現像処理を実行する際、画像処理エンジン１１８は、第１実施形態と同様、フレア補正処理部１２５として機能する。
【０１２８】
したがって、本実施形態の画像処理エンジン１１８は、電子カメラ１１１が撮影で取得したＲＡＷ画像や、メディアＩ／Ｆ１２１（又は通信I／F１２２）を介して読み込まれたＲＡＷ画像に対して、第１実施形態と同様の現像処理を施すことができる。
【０１２９】
なお、電子カメラ１１１が取得したＲＡＷ画像に対する現像処理のタイミングは、そのＲＡＷ画像が撮影された直後としてもよい。
【０１３０】
また、図２３において点線で示した部分（撮影光学系１１２を含むレンズユニット）は、電子カメラ１１１に対して交換可能であっても構わない。
【０１３１】
また、交換可能なレンズユニットで撮影されたＲＡＷ画像に対して現像処理を施す場合は、撮影光学系１１２の種類情報を、画像処理エンジン１１８がレンズユニットから読み込めばよい。
【０１３２】
また、本実施形態の現像処理も、第１実施形態の現像処理と同様に変形することが可能である。
【符号の説明】
【０１３３】
１１…コンピュータ、１２…データ読込部、１３…記憶装置、１４…ＣＰＵ、１５…メモリ、１６…入出力Ｉ／Ｆ、１７…バス
【図１】
2014165833.tif 000004
【図２】
2014165833.tif 000005
【図３】
2014165833.tif 000006
【図４】
2014165833.tif 000007
【図５】
2014165833.tif 000008
【図６】
2014165833.tif 000009
【図１１】
2014165833.tif 000010
【図１２】
2014165833.tif 000011
【図１５】
2014165833.tif 000012
【図１６】
2014165833.tif 000013
【図１７】
2014165833.tif 000014
【図１８】
2014165833.tif 000015
【図２１】
2014165833.tif 000016
【図２２】
2014165833.tif 000017
【図２３】
2014165833.tif 000018
【図７】
2014165833.tif 000019
【図８】
2014165833.tif 000020
【図９】
2014165833.tif 000021
【図１０】
2014165833.tif 000022
【図１３】
2014165833.tif 000023
【図１４】
2014165833.tif 000024
【図１９】
2014165833.tif 000025
【図２０】
2014165833.tif 000026