(57)【要約】【課題】既知の反射率を有する被写体の撮影や、照明の分光分布の測定を行うことなく、被写体に対する照明光源の分光分布情報を簡易に取得する。【解決手段】マルチバンド撮影を行う撮像部１０によって、被写体１２をストロボ光源１３下で撮影した画像データと、照明光源１４下で撮影した画像データについて、画像差分演算部１１１でこれらの画像差分を算出する。分光分布差分演算部１１３は、分光分布保持データベース１１２に保持された各種の分光分布情報について、ストロボ設定部１１６で設定されたストロボ光源１３の分光分布情報との差分である分光分布差分を算出する。そして差分比較部１１４で、画像差分との一致度が最も高い分光分布差分を選択し、分光分布推定部１１５で、該選択された分光分布差分に対応する分光分布情報を、照明光源１４に対する分光分布情報として推定する。【選択図】  図２
2011022044.tif 000002
【特許請求の範囲】
【請求項１】
被写体に対する照明光源の分光分布情報を推定する光源推定装置であって、
複数の分光特性によるマルチバンド撮影を行ってマルチバンド画像情報を得る撮像手段と、
複数種類の光源の分光分布情報を保持する分光分布保持手段と、
ストロボ光源の分光分布情報を設定するストロボ設定手段と、
前記撮像手段によって、前記被写体を前記ストロボ光源の下でマルチバンド撮影したマルチバンド画像情報と、前記被写体を前記照明光源の下でマルチバンド撮影したマルチバンド画像情報と、の差分であるマルチバンド画像差分を算出する画像差分算出手段と、
前記分光分布保持手段に保持された複数の分光分布情報のそれぞれについて、前記ストロボ光源の分光分布情報との差分である分光分布差分を算出する分光分布差分算出手段と、
該算出された複数の分光分布差分のそれぞれについて、前記マルチバンド画像差分と比較した一致度を算出し、該一致度が最も高い分光分布差分を選択する差分比較手段と、
前記差分比較手段で選択された分光分布差分に対応する分光分布情報を、前記照明光源に対する分光分布情報として推定する分光分布推定手段と、
を有することを特徴とする光源推定装置。
【請求項２】
前記ストロボ設定手段は、前記ストロボ光源の種類に応じた分光分布情報を前記分光分布保持手段より選出して設定することを特徴とする請求項１に記載の光源推定装置。
【請求項３】
前記ストロボ設定手段は、ユーザ指示に応じた種類の前記ストロボ光源の分光分布情報を設定することを特徴とする請求項２に記載の光源推定装置。
【請求項４】
前記分光分布保持手段は、前記複数種類の光源の分光分布情報として、それぞれの光源の分光分布に前記撮像手段の分光特性を乗じた値を全波長域において積分した値を保持することを特徴とする請求項１乃至３のいずれか１項に記載の光源推定装置。
【請求項５】
前記差分比較手段は、前記一致度として、前記分光分布差分と前記マルチバンド画像差分とにおけるバンド毎の差分を全バンドについて合計した値を算出することを特徴とする請求項１に記載の光源推定装置。
【請求項６】
さらに、前記差分比較手段で算出された最も高い一致度が予め定められたレベルに満たない場合に、前記分光分布保持手段に保持された複数の光源の組み合わせによる混合光源の複数組に対し、それぞれの混合分光分布情報を生成する混合分光分布の生成手段を有し、
前記分光分布差分算出手段は、前記混合分光分布情報のそれぞれについて、前記ストロボ光源の分光分布情報との差分である混合分光分布差分を算出し、
前記差分比較手段は、前記混合分光分布差分のそれぞれについて前記マルチバンド画像差分と比較した一致度を算出して該一致度が最も高い混合分光分布差分を選択し、
前記分光分布推定手段は、前記差分比較手段で選択された混合分光分布差分に対応する混合分光分布情報を、前記照明光源に対する分光分布情報として推定する
ことを特徴とする請求項１に記載の光源推定装置。
【請求項７】
前記混合分光分布の生成手段は、前記混合光源に対し、それぞれの光源に対する分光分布情報を予め定められたゲインおよび混合比に基づいて合成することによって、前記混合分光分布情報を算出することを特徴とする請求項６に記載の光源推定装置。
【請求項８】
被写体に対する照明光源の分光分布情報を推定する光源推定方法であって、
ストロボ光源の分光分布情報を設定するストロボ設定ステップと、
複数の分光特性によるマルチバンド撮影を行ってマルチバンド画像情報を得る撮像装置によって、前記被写体を前記ストロボ光源の下でマルチバンド撮影したマルチバンド画像情報と、前記被写体を前記照明光源の下でマルチバンド撮影したマルチバンド画像情報と、の差分であるマルチバンド画像差分を算出する画像差分算出ステップと、
分光分布データベースに保持された複数種類の光源の分光分布情報のそれぞれについて、前記ストロボ光源の分光分布情報との差分である分光分布差分を算出する分光分布差分算出ステップと、
該算出された複数の分光分布差分のそれぞれについて、前記マルチバンド画像差分と比較した一致度を算出し、該一致度が最も高い分光分布差分を選択する差分比較ステップと、
該選択された分光分布差分に対応する分光分布情報を、前記照明光源に対する分光分布情報として推定する分光分布推定ステップと、
を有することを特徴とする光源推定方法。
【請求項９】
コンピュータで実行されることにより、該コンピュータを請求項１乃至７のいずれか１項に記載の光源推定装置の各手段として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、撮影画像から光源の分光分布を推定する光源推定装置および光源推定方法に関する。
【背景技術】
【０００２】
撮影対象の分光波形を実用上十分な精度で復元可能とするマルチバンドカメラが知られている。マルチバンドカメラとしては、４種類以上の複数フィルタを用いるものが一般的であり、特殊な６バンドフィルタ配列とデモザイクによってワンショットでマルチバンド情報を取得するカメラも実用化されている。このようなマルチバンドカメラは、撮像されたマルチバンド画像に基づき、物体の分光反射率を正確に予測するために利用される。
【０００３】
マルチバンドカメラによって実際に取得されるマルチバンド情報は分光放射輝度であるため、該分光放射輝度を、光源の分光分布と物体の分光反射率データに分離する必要がある。そのため従来では、以下のような方法によって物体の分光反射率を推定していた。
【０００４】
まず、本撮影を行うに先立って、シーン中に分光反射率が既知である被写体を入れて撮影を行い、該撮影によって得られたマルチバンド情報を該被写体の分光反射率で除算することによって、光源の分光分布データを得ておく。そして本撮影においては、得られたマルチバンド情報を光源の分光分布データで除算することによって、該撮影された被写体の分光反射率が得られる。また、撮影時の光源の分光分布を測定器等により入手し、この測定データに基づいて、被写体の分光反射率を求めることも可能である。
【０００５】
このように、分光分布が既知である光源を用いたマルチバンド撮影によって被写体の分光分布を推定する方法として、特許文献１に記載された方法が知られている。この方法においては、液晶チューナブルフィルタと既知の分光分布を持つプローブ光源を利用している。すなわち、プローブ光源と照明光源を照射したマルチバンド情報と、照明光源のみを照射したマルチバンド情報との差分を、プローブ光の分光分布で除算することによって、被写体の分光分布を推定している。
【０００６】
以上説明したように従来のマルチバンド撮影では、被写体の分光反射率を推定する際に、光源の分光分布を利用していた。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００１−０８３００９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかしながら、上記特許文献１による分光分布の推定方法によれば、プローブ光源毎に撮影を行う必要があり、最大で１６回もの撮影を必要とする。そのため、動被写体についての撮影では、被写体の位置がずれてしまうため、正確なマルチバンド情報を取得することができない。また、撮影時間内に照明の条件が変動する場合にも、正確なマルチバンド情報は取得できない。更に、撮像バンド数を削減して撮像回数を削減した場合には、分光反射率の推定精度が著しく低下してしまう。また、使用する液晶チューナブルフィルタや各波長のプローブ光が安定するまでには所定の時間を要するため、これら機器に精通していない一般のユーザでは容易に操作することはできなかった。
【０００９】
上述したように、一般のユーザでは、マルチバンドカメラによる撮影時に、光源の分光分布の測定用に測定器を用意して操作することは困難であり、さらに、撮影時に既知の分光反射率の被写体を撮影することも困難である。このような場合には光源の分光分布が得られないため、被写体の分光情報を推定することができない。そこで、測定器の扱いに精通していない一般ユーザでも、光源の分光分布を取得できるような技術が求められる。
【００１０】
本発明は上述した問題を解決するためになされたものであり、被写体に対する光源の分光分布を簡易に取得可能とする光源推定装置および光源推定方法を提供することを目的とする。
【課題を解決するための手段】
【００１１】
上記目的を達成するための一手段として、本発明の光源推定装置は以下の構成を備える。
【００１２】
すなわち、被写体に対する照明光源の分光分布情報を推定する光源推定装置であって、複数の分光特性によるマルチバンド撮影を行ってマルチバンド画像情報を得る撮像手段と、複数種類の光源の分光分布情報を保持する分光分布保持手段と、ストロボ光源の分光分布情報を設定するストロボ設定手段と、前記撮像手段によって、前記被写体を前記ストロボ光源の下でマルチバンド撮影したマルチバンド画像情報と、前記被写体を前記照明光源の下でマルチバンド撮影したマルチバンド画像情報と、の差分であるマルチバンド画像差分を算出する画像差分算出手段と、前記分光分布保持手段に保持された複数の分光分布情報のそれぞれについて、前記ストロボ光源の分光分布情報との差分である分光分布差分を算出する分光分布差分算出手段と、該算出された複数の分光分布差分のそれぞれについて、前記マルチバンド画像差分と比較した一致度を算出し、該一致度が最も高い分光分布差分を選択する差分比較手段と、前記差分比較手段で選択された分光分布差分に対応する分光分布情報を、前記照明光源に対する分光分布情報として推定する分光分布推定手段と、を有することを特徴とする。
【発明の効果】
【００１３】
上記構成からなる本発明によれば、被写体に対する光源の分光分布を簡易に取得することができる。
【図面の簡単な説明】
【００１４】
【図１】本発明に係る一実施形態における光源推定装置の構成を示すブロック図、
【図２】本実施形態における光源推定処理を示すフローチャート、
【図３】本実施形態におけるストロボ設定ＵＩを説明する図、
【図４】本実施形態における分光分布差分算出処理を示すフローチャート、
【図５】本実施形態における分光分布データベースの模式図、
【図６】本実施形態における光源比較処理を示すフローチャート、
【図７】第２実施形態における光源推定装置の構成を示すブロック図、
【図８】第２実施形態における光源推定処理を示すフローチャート、
【図９】第２実施形態における仮想ＭＩＸ光源比較処理を示すフローチャート、である。
【発明を実施するための形態】
【００１５】
以下、添付の図面を参照して、本発明をその好適な実施形態に基づいて詳細に説明する。なお、以下の実施形態において示す構成は一例に過ぎず、本発明は図示された構成に限定されるものではない。
【００１６】
＜第１実施形態＞
●システム構成
図１は、本実施形態における光源推定装置の構成を示すブロック図である。同図において、１０は撮像部、１１は画像処理部、１２は被写体、１３はストロボ光源、１４は照明光源、である。撮像部１０は、例えば６種類の分光特性によるマルチバンド撮影を行ってマルチバンド画像情報を得る、６バンドカメラである。本実施形態では、ストロボ光源１３および照明光源１４でそれぞれ照射された被写体１２に対し、撮像部１０のマルチバンド撮影によって得られたマルチバンド画像情報（ストロボ画像情報および照明画像情報）を画像処理部１１に転送する。画像処理部１１ではこれらのマルチバンド画像情報から照明光源の分光分布を推定する。
【００１７】
画像処理部１１は概略として以下の構成を有する。まず、ストロボ画像情報と照明画像情報の差分（マルチバンド画像差分）を算出する画像差分演算部１１１と、複数種類の光源に対する分光分布情報を格納する、分光分布保持手段としての分光分布データベース１１２を有する。そして、ストロボ分光分布と分光分布データベース１１２内の照明分光分布の差分を算出する分光分布差分演算部１１３と、マルチバンド画像差分と分光分布差分を比較する差分比較部１１４を有する。さらに、差分比較結果に基づいて照明光源１４の分光分布を推定する分光分布推定部１１５と、ストロボ光源１３の分光分布情報をユーザ指示に応じて設定するストロボ設定部１１６、を有する。
【００１８】
なお、ストロボ光源１３としては、大型フラッシュ、小型フラッシュ、マクロ用フラッシュ、内蔵フラッシュのいずれでもかまわないが、広い波長域にスペクトルをもつブロードな特性の光源であることが望ましい。本実施形態では、ストロボ光源１３として小型フラッシュのキセノンガス、グリップオンタイプを想定している。分光分布データベース１１２には、ストロボ光源１３として使用可能な様々なストロボ光源の分光分布をはじめ、標準光源、各種蛍光灯、特殊光源の分光反射率、等が格納されている。
【００１９】
●光源推定処理
上記構成における光源推定処理について、図２のフローチャートを用いて説明する。まずＳ２０１において、ストロボ設定部１１６で制御されるＵＩにより、使用するストロボ光源１３の種類を設定する。このストロボ設定の詳細については後述するが、この処理によってストロボ光源１３の分光分布が設定される。その後、Ｓ２０２で撮像部１０がストロボ光源１３を用いて、被写体１２のマルチバンド撮影を行うことによって、ストロボ画素値Ｐs_n（ｎはバンド１〜６を示す）を得る。そしてＳ２０３では、撮像部１０が照明光源１４の下で被写体１２のマルチバンド撮影を行うことによって、照明下画素値Ｐx_nを得る。
【００２０】
するとＳ２０４でマルチバンド画像差分算出処理を行う。すなわち、画像差分演算部１１１において、Ｓ２０２で撮影したストロボ画素値Ｐs_nとＳ２０３で撮影した照明下画素値Ｐx_nとの差分（マルチバンド画像差分Ｄp）を、以下の式(1)より求める。
【００２１】
Ｄp＝Σ(Ｐs_n−Ｐx_n)      ・・・(1)
次にＳ２０５では分光分布差分算出処理を行う。すなわち、分光分布差分演算部１１３において、ストロボ光源１３の分光分布Ｉs_nと、複数の分光分布Ｉi_n（ｉは分光分布番号）に対する差分（分光分布差分Ｄs_i）を、以下の式(2)より算出する。なお、これらの分光分布Ｉs_n，Ｉi_nは、分光分布データベース１１２に予め保持されている。この分光分布差分算出処理の詳細については後述する。
【００２２】
Ｄs_i＝Σ(Ｉs_n−Ｉi_n)      ・・・(2)
次にＳ２０６で差分比較部１１４において、Ｓ２０４で求めたマルチバンド画像差分ＤpとＳ２０５で求めた分光分布差分Ｄs_iとを比較して、差分評価値Ｐiを求める。この比較処理の詳細については後述する。
【００２３】
そしてＳ２０７で分光分布推定部１１５において、Ｓ２０６で得られた差分評価値Ｐiに基づいて、分光分布データベース１１２に保持されている複数の分光分布のうち、照明光源１４に最も近いと判定される光源を選出する。このように選出された光源が、照明光源１４として推定されることになる。
【００２４】
●ストロボ設定ＵＩ
図３（ａ）に、上記Ｓ２０１においてストロボ情報を指定するためのＵＩ例を示す。同図において、２１０はストロボ形態選択コンボボックスであり、ユーザが使用するストロボ光源１３について、外付タイプと内蔵タイプのいずれかを選択すると、該選択された項目が反転表示される。２１１はストロボ種類選択コンボボックスであり、ユーザがストロボ光源１３の型番を選択すると、選択した項目が反転表示される。２１２は分光分布表示ウィンドウであり、設定されたストロボ光源１３に該当する分光分布が、分光分布データベース１１２から選択されて表示される。このように、選択したストロボの分光分布がＵＩ表示されることにより、該ストロボが照明光源１４を推定するのに適しているかをユーザが判断する一助となる。
【００２５】
２１３は選択ボタンであり、ユーザによって押下されると、各コンボボックスの設定内容に応じたストロボの分光分布情報がメモリに保持される。２１４はキャンセルボタンであり、ユーザによって押下されるとメモリの解放等を行って設定ＵＩのダイアログウィンドウを終了する。
【００２６】
以下、上記ストロボ設定ＵＩにおけるストロボ種類設定動作について、図３（ｂ）の状態遷移図を用いて説明する。まずステート２２０では、初期設定の読み込みや、図２ＢのＵＩを表示する等の初期化動作を行い、ステート２２１へ移行する。
【００２７】
ステート２２１では、図３（ａ）のＵＩにおけるユーザ操作待ち状態となる。ここで、ストロボ形態選択コンボボックス２１０にて項目選択が行われるとステート２２２へ移行する。同じく、ストロボ種類選択コンボボックス２１１にて項目選択が行われるとステート２２３へ移行する。また、選択ボタン２１３が押下されるとステート２２５へ移行し、キャンセルボタン２１４が押下されるとステート２２６へ移行する。
【００２８】
ステート２２２では、ストロボ形態選択コンボボックス２１０にて選択された項目を反転表示し、ステート２２３へ移行する。ステート２２３では、ストロボ種類選択コンボボックス２１１にて選択された項目を反転表示し、ステート２２４へ移行する。
【００２９】
ステート２２４では、ステート２２２およびステート２２３にて選択された項目、すなわち設定されたストロボ情報に従って、対応する分光分布を分光分布データベース１１２から選択して描画ウィンドウ２１２に描画し、ステート２２１へ移行する。
【００３０】
ステート２２５では、設定したストロボの分光分布をメモリ（不図示）に格納し、ステート２２６へ移行する。そしてステート２２６では、終了に関する動作を行う。
【００３１】
なお、本実施形態ではストロボの形態や種類を図３（ａ）のＵＩにより取得する例を示したが、取得方法はこの例に限定されない。例えば、通信ケーブルを介して外部より取得するようにしても良い。
【００３２】
●分光分布差分算出処理
以下、上記Ｓ２０５における、ストロボ光源１３の分光分布Ｉs_nと分光分布データベース１１２に保持された複数の分光分布Ｉi_nとの差分である、分光分布差分Ｄs_iの算出処理について、図４のフローチャートを用いて説明する。
【００３３】
まずＳ３０１において、分光分布データベース１１２に保持された光源の番号を表す変数ｉを１に設定し、メモリ確保等の初期化処理を行う。次にＳ３０２において、ストロボ光源１３の分光分布として以下の式(3)より、ストロボ分光分布Ｌs(λ)とフィルタ特性Ｔn(λ)の乗算値の全波長域における積分値（６バンド積分値）を算出する。
【００３４】
Ｉs_n＝Σ{Ｌs(λ)×Ｔn(λ)}      ・・・(3)
次にＳ３０３では、光源ｉの分光分布として以下の式(4)より、光源ｉの分光分布Ｌi(λ)とフィルタ特性Ｔn(λ)の乗算値の全波長域における積分値（６バンド積分値）を算出する。なお、光源ｉの分光分布Ｌi(λ)は分光分布データベース１１２に格納されており、例えば図５に示すように、３８０〜７８０nmの１０nm間隔にて、各波長域のスペクトル情報が光源番号の順にラスタライズされているとする。
【００３５】
Ｉi_n＝Σ{Ｌi(λ)×Ｔn(λ)}      ・・・(4)
次にＳ３０４では、上記式(2)により、Ｓ３０２による積分値Ｉs_nとＳ３０３による積分値Ｉi_nの差分Ｄs_i（６バンド差分）を算出し、メモリに格納する。
【００３６】
そしてＳ３０５において、分光分布データベース１１２内の全ての光源について処理を行ったか否かを判定し、未処理の光源が残っていればｉに１を加算してＳ３０３へ戻り、残っていなければ終了に関する処理を行う。
【００３７】
●光源比較処理
以下、上記Ｓ２０６における、マルチバンド画像差分Ｄpと分光分布差分Ｄs_iによる光源比較処理について、図６のフローチャートを用いて説明する。
【００３８】
まずＳ５０１において、分光分布データベース１１２に保持された光源の番号を表す変数ｉに１を設定し、後述する差分評価値の最小値を格納するための変数ＭＩＮに所定の定数ｋを設定して、メモリ確保等の初期化処理を行う。次にＳ５０２において、以下の式(5)により、差分評価値Ｐiを算出する。式(5)に示す差分評価値Ｐiは、バンド毎の差分として示されるマルチバンド画像差分Ｄpと分光分布差分Ｄs_iの一致度（差分ノルム）を、全バンドについて合計した値として算出される。差分評価値Ｐiは、その値が小さいほど、マルチバンド画像差分Ｄpと分光分布差分Ｄs_iの一致度が高いことを示す。
【００３９】
Ｐi＝Ｄp−Ｄs_i      ・・・(5)
次にＳ５０３において、差分評価値ＰiがＭＩＮより小さければＳ５０４へ進み、ＭＩＮの値を差分評価値Ｐiで置き換え、この際の光源ＩＤ、すなわち分光分布データベース１１２において変数ｉに対応する光源のＩＤをメモリに格納した後、Ｓ５０５に進む。一方、Ｓ５０３で差分評価値ＰiがＭＩＮ以上であれば、そのままＳ５０５へ進む。Ｓ５０５では、分光分布データベース１１２内の全ての光源について処理を行ったか否かを判定し、未処理の光源が残っていればｉに１を加算してＳ５０２へ戻り、残っていなければ終了に関する処理を行う。
【００４０】
図６のフローチャートに示す処理が終了した時点で、Ｓ５０４にて最終的に設定された最小の差分評価値Ｐiに対応する光源ＩＤが、Ｓ２０７で照明光源１４として選出される。
【００４１】
このように本実施形態においては、既知のストロボ光源に対する差分であるマルチバンド画像差分Ｄpと、同じく既知のストロボ光源に対する差分である分光分布差分Ｄs_iについて、その一致度（ノルム）を検証する。そして、分光分布データベース１１２内で該一致度が最も高い光源が、照明光源１４として選出される。
【００４２】
以上説明した様に本実施形態によれば、ストロボ光と照明光による実際のマルチバンド撮影情報と、ストロボおよび所定光源の既知の分光分布との一致度に基づき、照明光源の分光分布を推定する。これにより、既知の反射率を持つ被写体を撮影したり、測定器による分光分布測定を行うことなく、照明光源の分光分布を推定することができる。
【００４３】
更に、様々な光源情報を有する分光分布データベースを参照するため、様々な光源に対して高精度な推定を行うことが可能であり、全波長での撮像を行わなくとも、１〜数回の撮影で、照明の分光分布を簡易に推定することができる。
【００４４】
なお、本実施形態では一致度として６バンドの差分ノルムを用いる例を示したが、この差分に代えて、６バンドから算出される色差を用いることも可能である。更に、６バンドの差分ノルムと６バンドから算出される色差との組み合わせからなる差分評価値等、様々な評価値を用いることが可能である。
【００４５】
また、データベース内のストロボ光源については、各機種ごとに分光分布を格納する例を示したが、個体差を考慮して測定値に基づく分光分布をストロボ個体ごとに格納することも可能である。また、経時変化を考慮して分光分布情報を更新することも有効である。
【００４６】
＜第２実施形態＞
以下、本発明に係る第２実施形態について説明する。上述した第１実施形態では、照明が単一光源であると仮定して分光分布を推定する例を示したが、第２実施形態ではさらに、複数の光源、すなわち混合光源（以下、ＭＩＸ光源）の分光分布も推定可能とする例を示す。以下では、特に第１実施形態と異なる部分のみについて説明する。
【００４７】
図７は、第２実施形態における光源推定装置の構成を示すブロック図である。同図においては画像処理部６の構成が、上述した第１実施形態での図１に示した画像処理部１１と異なる。画像処理部６は、概略として以下の構成を有する。まず、上述した第１実施形態と同様に、画像差分演算部６０１と分光分布データベース６０２、分光分布差分演算部６０３、差分比較部６０４、分光分布推定部６０５、およびストロボ設定部６０６を有する。そして第２実施形態における特徴的な構成として、ＭＩＸ光源推定を行うか否かの判定処理を行うＭＩＸ光源推定判定部６０７と、ＭＩＸ光源の分光分布を生成するＭＩＸ光源分光分布生成部６０８を有する。
【００４８】
●光源推定処理
以下、第２実施形態における光源推定処理について、図８のフローチャートを用いて説明する。まずＳ７０１において、ストロボ設定部６０６で制御されるＵＩ等によって、使用するストロボ光源１３の情報を設定する。その後、Ｓ７０２で撮像部１０がストロボ光源１３を用いて、被写体のマルチバンド撮影を行う。そしてＳ７０３では、撮像部１０が照明光源１４の下で被写体のマルチバンド撮影を行う。するとＳ７０４で画像差分演算部６０１において、Ｓ７０２で撮影した画素値とＳ７０３で撮影した画素値との差分（マルチバンド画像差分Ｄp）を、上記式(1)より求める。次にＳ７０５では分光分布差分演算部６０３において、予め保持しているストロボの分光分布と、分光分布データベース６０２に保持されている複数の分光分布に対する差分（分光分布差分Ｄs_i）を、上記式(2)より算出する。次にＳ７０６で差分比較部６０４において、Ｓ７０４で求めたマルチバンド画像差分ＤpとＳ７０５で求めた分光分布差分Ｄs_iとを比較して、差分評価値Ｐiを求める。ここまでは、上述した第１実施形態に示した図２のＳ２０１〜Ｓ２０６と同様である。
【００４９】
第２実施形態では次にＳ７０７でＭＩＸ光源推定判定部６０７において、Ｓ７０６で得られた差分評価値Ｐiに基づいて、ＭＩＸ光源推定処理の可否を判定する。すなわち、最小の差分評価値Ｐiが所定の値よりも大きい場合に、マルチバンド画像差分Ｄpと分光分布差分Ｄs_iとにおいて満足する一致度が得られなかったとして、Ｓ７０８に進む。Ｓ７０８では仮想ＭＩＸ光源についての比較処理を行うが、その詳細については後述する。一方、Ｓ７０７で満足する一致度が得られた場合には、Ｓ７０８の仮想ＭＩＸ光源比較処理を行わずに、そのままＳ７０９へ進む。
【００５０】
そしてＳ７０９では分光分布推定部６０５において、Ｓ７０６またはＳ７０８で得られた差分評価値Ｐiに基づいて、分光分布データベース６０２に保持されている複数の分光分布から、照明光源１４に最も近いと判定される光源を選出する。すなわち第２実施形態では、差分評価値Ｐiに対応する単独光源またはＭＩＸ光源に係る光源が、照明光源１４として選出される。
【００５１】
●仮想ＭＩＸ光源比較処理
以下、上記Ｓ７０８における仮想ＭＩＸ光源比較処理について、図９のフローチャートを用いて説明する。まずＳ８０１でＭＩＸ光源分光分布生成部６０８において、分光分布データベース６０２に保持された複数の分光分布の組み合わせによるＭＩＸ光源の複数組について、以下の式(6)により、仮想ＭＩＸ光源の分光分布（混合分光分布情報）を生成する。式(6)は、ｊ番目の組み合わせとしてのｉ1光源とｉ2光源のＭＩＸ光源による分光分布モデルを示し、α，βはゲイン、kは混合比である。すなわち仮想ＭＩＸ光源の分光分布は、ＭＩＸ光源を構成する各光源に対する分光分布情報を、ゲインα，βおよび混合比ｋに基づいて合成することによって算出される。なお、ここでＭＩＸ光源としては、分光分布データベース６０２に保持された全ての分光分布の組み合わせとしても良いし、所定の分光分布の組み合わせのみに対応させても良い。
【００５２】
Ｉmix_j_n＝Σ{(α×ｋ×Ｌi1(λ)＋β×(1-ｋ)×Ｌi2(λ))×Ｔn(λ)}      ・・・(6)
次にＳ８０２で分光分布差分演算部６０３において、予め保持しているストロボの分光分布Ｉs_nと、Ｓ８０１で生成された複数の仮想ＭＩＸ光源の分光分布Ｉmix_j_nの差分（混合分光分布差分）Ｄs_jを、以下の式(7)より算出する。
【００５３】
Ｄs_j＝Σ(Ｉs_n−Ｉmix_j_n)      ・・・(7)
そしてＳ８０３において差分比較部６０４が、上述した第１実施形態の図６と同様の方法で、Ｓ７０４で求めたマルチバンド画像差分Ｄpと、Ｓ８０２で求めた仮想ＭＩＸ光源の分光分布差分Ｄs_jの全てを比較し、新たな差分評価値Ｐiを得る。
【００５４】
以上説明した様に第２実施形態によれば、ストロボ光と照明光による実際の撮影情報と、ストロボと仮想ＭＩＸ光源の分光分布とに基づき、照明光源の分光分布を推定することができる。これにより、データベースに保持された所定光源からは十分な推定が行えない場合にも、それらの組み合わせからなる仮想ＭＩＸ光に対する推定が可能となるため、より詳細な推定を行うことができる。
【００５５】
なお、第２実施形態では、差分一致度が所定のレベルを満たしていない場合にＭＩＸ光源推定を行う例を示したが、２種類以上の光源に対して同等の最良一致度をもつ光源があった場合に、該ＭＩＸ光源推定処理を実施するようにしても良い。
【００５６】
また、データベース内の全ての光源の組み合わせから仮想ＭＩＸ光を生成する例を示したが、一致度が最も高い光源と２番目に高い光源に限定して仮想ＭＩＸ光を生成することで、処理時間を短縮することができる。
【００５７】
＜その他の実施形態＞
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（またはＣＰＵやＭＰＵ等）がプログラムを読み出して実行する処理である。
【図１】
2011022044.tif 000003
【図２】
2011022044.tif 000004
【図３】
2011022044.tif 000005
【図４】
2011022044.tif 000006
【図５】
2011022044.tif 000007
【図６】
2011022044.tif 000008
【図７】
2011022044.tif 000009
【図８】
2011022044.tif 000010
【図９】
2011022044.tif 000011