(57)【特許請求の範囲】
【請求項１】
シーンのマルチスペクトル画像から観測スペクトル情報を選択し、
照明の基底ベクトルと表面反射率の基底ベクトルを用いて、照明の分光分布の制約範囲と物体の表面反射率の制約範囲を生成し、
前記観測スペクトル情報と、前記照明の分光分布の制約範囲と、前記物体の表面反射率の制約範囲とを用いて、前記照明の基底ベクトルの重み係数と前記表面反射率の基底ベクトルの重み係数とをパラメタとする観測方程式を生成し、
前記観測方程式から前記シーンの照明の分光分布と物体の表面反射率とを算出する
分光画像処理方法。
【請求項２】
前記照明の分光分布の制約範囲と前記物体の表面反射率の制約範囲のもと、前記観測方程式の重み係数に初期パラメタを与えて、観測スペクトルの推定値を算出し、
前記観測スペクトル情報と前記観測スペクトルの推定値の誤差に基づき、前記重み係数を更新し、
前記誤差が所定内となったときの前記観測スペクトルの推定値を構成する照明の分光分布と物体の表面反射率とを、前記シーンにおける照明の分光分布と物体の表面反射率として算出する
請求項１に記載の分光画像処理方法。
【請求項３】
前記観測スペクトル情報として、前記マルチスペクトル画像から線形独立なスペクトル情報を選択する
請求項１又は請求項２に記載の分光画像処理方法。
【請求項４】
前記照明の基底ベクトルとして、日射スペクトルの直達成分ベクトルと散乱成分ベクトルとを用いる
特徴とする請求項１から請求項３のいずれかに記載の分光画像処理方法。
【請求項５】
照明色が同一とみなせる範囲で異なるシーンを複数回撮影し、前記観測方程式を生成するに必要な数の前記線形独立なスペクトル情報を得る
請求項３に記載の分光画像処理方法。
【請求項６】
シーンのマルチスペクトル画像から観測スペクトル情報を選択するスペクトル情報選択手段と、
照明の基底ベクトルと表面反射率の基底ベクトルを用いて、照明の分光分布の制約範囲と物体の表面反射率の制約範囲を生成する制約範囲生成手段と、
前記観測スペクトル情報と、前記照明の分光分布の制約範囲と、前記物体の表面反射率の制約範囲とを用いて、前記照明の基底ベクトルの重み係数と前記表面反射率の基底ベクトルの重み係数とをパラメタとする観測方程式を生成する観測方程式生成手段と、
前記観測方程式から前記シーンの照明の分光分布と物体の表面反射率とを算出する算出手段と
を有する分光画像処理装置。
【請求項７】
前記算出手段は、
前記照明の分光分布の制約範囲と前記物体の表面反射率の制約範囲のもと、前記観測方程式の重み係数に初期パラメタを与えて、観測スペクトルの推定値を算出する手段と、
前記観測スペクトル情報と前記観測スペクトルの推定値の誤差に基づき、前記重み係数を更新する手段と、
前記誤差が所定内となったときの前記観測スペクトルの推定値を構成する照明の分光分布と物体の表面反射率とを前記シーンにおける照明の分光分布と物体の表面反射率として算出する手段と
を有する請求項６に記載の分光画像処理装置。
【請求項８】
前記スペクトル情報選択手段は、前記マルチスペクトル画像から線形独立なスペクトル情報を選択する
請求項６又は請求項７に記載の分光画像処理装置。
【請求項９】
前記スペクトル情報選択手段は、前記照明の基底ベクトルとして、日射スペクトルの直達成分ベクトルと散乱成分ベクトルとを用いる
特徴とする請求項６から請求項８のいずれかに記載の分光画像処理装置。
【請求項１０】
コンピュータに、
シーンのマルチスペクトル画像から観測スペクトル情報を選択する処理と、
照明の基底ベクトルと表面反射率の基底ベクトルを用いて、照明の分光分布の制約範囲と物体の表面反射率の制約範囲を生成する処理と、
前記観測スペクトル情報と、前記照明の分光分布の制約範囲と、前記物体の表面反射率の制約範囲とを用いて、前記照明の基底ベクトルの重み係数と前記表面反射率の基底ベクトルの重み係数とをパラメタとする観測方程式を生成する処理と、
前記観測方程式から前記シーンの照明の分光分布と物体の表面反射率とを算出する処理と
を実行させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、分光画像処理方法、分光画像処理装置およびプログラムに関する。
【背景技術】
【０００２】
撮影装置を用いて撮影された画像の色情報は、撮影時の環境照明光が物体の表面で反射され、撮影装置に入射した光をもとに生成される。そのため、表面反射率が同一である物体を撮影しても、環境照明光が異なると撮影画像の色情報も変化する。多数の波長の光が記録されたマルチスペクトル画像を用いて、撮影シーン中の物体を観測スペクトルの形状を手掛かりに認識する技術において、同一の物体であるにも関わらず、環境照明光の変化によって観測スペクトルが異なると、物体の識別が困難になる。この問題を解決するには、観測スペクトルから、環境照明光の影響を除去し、物体の表面反射率を得る必要がある。すなわち、観測スペクトルから照明と物体の分光特性を分離、復元する必要がある。【０００３】
画像などの、シーンから観測される色情報を手掛かりに、撮影シーン中の照明の分光分布や物体の表面反射率を、自動で推定する方法として、シーン中の物体色または照明色に関する仮定を置き、その仮定を手掛かりとして、照明と物体の分光特性を推定する方法が、既に提案されている。【０００４】
例えば、特許文献１に記載されている方法では、環境照明光が白色光であることと、反射率の平均値が灰色となることと、人の視覚特性としてエッジ領域の情報をもとに色を充填することを仮定している。この方法では、照明の分光分布と物体の表面反射率が予め保持されたそれぞれの主成分ベクトルと平均ベクトルの線型和となり、かつ上記仮定を満たすために必要なエネルギーを定義し、このエネルギーを最適化することで、照明の分光分布と物体の表面反射率を推定する。【０００５】
また、特許文献２に記載されている方法では、撮影画像中において肌色もしくは、または、グレーの物体の占める面積が広いことを仮定し、照明の範囲を黒体軌跡上の点に対応付け、色温度を変化させながら、仮定を満たすのに必要なエネルギーを算出、最適化することで色温度を推定する。この方法で推定された色温度をもとに照明を推定することが可能となる。【０００６】
上記関連技術による照明・反射率推定方法の構成図を図１１に示す。図１１の構成図は、特許文献１および特許文献２をもとに生成したブロック図である。上記関連技術は、色情報取得手段１と照明・反射率主成分ベクトル保存メモリ２、分光推定手段３とを備える。【０００７】
色情報取得手段１は、シーンから色情報を取得し出力する。照明・反射率主成分ベクトル保存メモリ２は、照明の分光分布および物体の表面反射率の制約範囲を示すため、主成分ベクトルと平均ベクトルを保持する。分光推定手段３は、照明の分光分布や物体の表面反射率をそれぞれの主成分ベクトルと平均ベクトルの線型和としつつ、予め設定された物体色もしくは照明色の仮定を満たすために必要なエネルギーを、色情報を用いて算出し、算出されたエネルギーを最適化する。更に、算出されたエネルギーが最適となるときの照明の分光分布と物体の表面反射率を推定値として出力する。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開平０７−０６６９８６号公報
【特許文献２】特開２００３−２０９８５６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００９】
画像として観測された色情報からシーン中の照明色と物体色を分離し、それぞれの分光特性を復元することは不良設定問題となる。この問題を解決するために従来技術では、シーン中の物体色または照明色に関する仮定を設定することで、照明の分光分布と物体の表面反射率を復元する。そのため、設定した仮定を満たすような限られたシーンについては、精度よく照明の分光分布と物体の表面反射率を推定することができるが、仮定と一致しないシーンについては、照明および物体の分光特性の推定精度が悪化するという課題があった。【００１０】
本発明の目的は、上記課題に鑑みて発明されたものであって、その目的は、シーン中の物体色および照明色に関する仮定を用いることなく観測スペクトルから物体および照明の分光特性を精度よく復元する分光画像処理方法、分光画像処理装置およびプログラムを提供することにある。
【課題を解決するための手段】
【００１１】
本発明は、シーンのマルチスペクトル画像から観測スペクトル情報を選択し、照明の基底ベクトルと表面反射率の基底ベクトルを用いて、照明の分光分布の制約範囲と物体の表面反射率の制約範囲を生成し、前記観測スペクトル情報と、前記照明の分光分布の制約範囲と、前記物体の表面反射率の制約範囲とを用いて、前記照明の基底ベクトルの重み係数と前記表面反射率の基底ベクトルの重み係数とをパラメタとする観測方程式を生成し、前記観測方程式から前記シーンの照明の分光分布と物体の表面反射率とを算出する分光画像処理方法である。【００１２】
本発明は、シーンのマルチスペクトル画像から観測スペクトル情報を選択するスペクトル情報選択手段と、照明の基底ベクトルと表面反射率の基底ベクトルを用いて、照明の分光分布の制約範囲と物体の表面反射率の制約範囲を生成する制約範囲生成手段と、前記観測スペクトル情報と、前記照明の分光分布の制約範囲と、前記物体の表面反射率の制約範囲とを用いて、前記照明の基底ベクトルの重み係数と前記表面反射率の基底ベクトルの重み係数とをパラメタとする観測方程式を生成する観測方程式生成手段と、前記観測方程式から前記シーンの照明の分光分布と物体の表面反射率とを算出する算出手段とを有する分光画像処理装置である。【００１３】
本発明は、コンピュータに、シーンのマルチスペクトル画像から観測スペクトル情報を選択する処理と、照明の基底ベクトルと表面反射率の基底ベクトルを用いて、照明の分光分布の制約範囲と物体の表面反射率の制約範囲を生成する処理と、前記観測スペクトル情報と、前記照明の分光分布の制約範囲と、前記物体の表面反射率の制約範囲とを用いて、前記照明の基底ベクトルの重み係数と前記表面反射率の基底ベクトルの重み係数とをパラメタとする観測方程式を生成する処理と、前記観測方程式から前記シーンの照明の分光分布と物体の表面反射率とを算出する処理とを実行させるプログラムである。
【発明の効果】
【００１４】
本発明によれば、観測されたシーンに依らず、得られた分光スペクトル情報を手掛かりに、照明の分光分布および物体の表面反射率を安定して精度よく推定することができる。
【図面の簡単な説明】
【００１５】
【図１】図１は本発明による第１の実施の形態の概要の構成図である。
【図２】図２は本発明による第１の実施の形態の構成図である。
【図３】図３は本発明による第１の実施の形態のフローチャートである。
【図４】図４は本発明による第２の実施の形態の構成図である。
【図５】図５は本発明による第２の実施の形態のフローチャートである。
【図６】図６は本発明による第３の実施の形態の構成図である。
【図７】図７は本発明による第３の実施の形態のフローチャートである。
【図８】図８は照明分光分布の基底ベクトルの例である。
【図９】図９はマクベスカラーチャート２４色の直流成分および主成分ベクトルから作られる表面反射率の基底ベクトルの例である。
【図１０】図１０は照明分光分布の制約範囲のイメージを表した図である。
【図１１】図１１は関連する技術の照明・反射率推定方法の構成図である。
【発明を実施するための形態】
【００１６】
本発明は、シーンから観測された色情報を手掛かりに、シーン中の照明と物体の分光特性を推定する分光画像処理技術に関し、特に、高い波長分解能を持つカメラで観測された分光スペクトルを用いて、シーン中の物体色または照明色に関する仮定を置くことなく、照明の分光分布と物体の表面反射率を推定する分光画像処理技術に関する。
【００１７】
以下、本発明の実施の形態について図面を参照して詳細に説明する。
【００１８】
[第１の実施の形態]
図１は、本発明による第１の実施の形態の概要を示すブロック図である。図１に示す第１の実施の形態は、マルチスペクトル情報取得手段１１と、照明・反射率制約範囲保持メモリ１２と、色仮定なし分光推定手段１３を備える。
【００１９】
本発明による第１の実施の形態を構成する各手段の概要を説明する。
【００２０】
マルチスペクトル情報取得手段１１は、シーンから取得したマルチスペクトル情報をもとに、最適化に十分な数のスペクトル情報を抽出し、出力する。
【００２１】
照明・反射率制約範囲保持メモリ１２は、照明の分光分布と物体の表面反射率の制約範囲を保持する。
【００２２】
色仮定なし分光推定手段１３では、マルチスペクトル情報取得手段１１で算出された最適化に十分な数のスペクトル情報をもとに、照明・反射率制約範囲保持メモリ１２に保持された照明の分光分布と物体の表面反射率の制約範囲を満たすように、観測方程式を生成し、この観測方程式を解くことで、照明の分光分布と物体の表面反射率とを算出する。そして、色仮定なし分光推定手段１３は、算出した照明の分光分布と物体の表面反射率とを出力する。
【００２３】
本発明による第１の実施の形態は、上述した関連技術の構成と比較して、マルチスペクトル情報取得手段１１と、照明・反射率制約範囲保持メモリ１２と、色仮定なし分光推定手段１３が異なる。以下、マルチスペクトル情報取得手段１１と、照明・反射率制約範囲保持メモリ１２と、色仮定なし分光推定手段１３の構成の詳細を説明する。
【００２４】
マルチスペクトル情報取得手段１１は、図１１における色情報取得手段１と同様に、シーンからマルチスペクトル情報を得る。ただし、色情報取得手段１では、出力する色情報の数を考慮しないのに対し、マルチスペクトル情報取得手段１１は、後段の色仮定なし分光推定手段１３での処理を考慮し、最適化に必要な数のスペクトル情報を出力する点で異なる。ここでNは観測される分光スペクトルのバンド数、Mは観測される分光スペクトルの数、Pallは推定する全パラメタ数、PIは一つの照明の分光分布を再現するのに必要なパラメタ数、PRは一つの物体の表面反射率を再現するのに必要なパラメタ数を示すとすると、例えば、式（１）を満たす十分な数の線形独立なスペクトル情報を出力する。
【００２５】
【数１】
0006015964.tif 000002
図２は、本発明による第１の実施の形態を示すブロック図である。マルチスペクトル情報取得手段１１の内部の構成を説明する。図２に示すマルチスペクトル情報取得手段１１は、マルチスペクトル画像撮影手段１１１とスペクトル選択手段１１２とを備える。
【００２６】
マルチスペクトル画像撮影手段１１１は、シーンのマルチスペクトル情報が記録されたマルチスペクトル画像を出力する。
【００２７】
スペクトル選択手段１１２は、マルチスペクトル画像撮影手段１１１によって得たマルチスペクトル画像のデータの中から、例えば、式（１）の条件を満たす、線形独立な関係にある複数のスペクトル情報を選択し出力する。
【００２８】
図１１における分光推定手段３では、予め設定した物体色または照明色に関する仮定に基づいて、シーン中の照明の分光分布および物体の表面反射率を推定する。具体的には、設定した物体色または照明色に関する仮定を満たすのに必要なエネルギーを定義し、これを最適化することにより、シーン中の照明の分光分布および物体の表面反射率を推定する。例えば、「物体色の平均が灰色である」（グレイワールド仮説）、「画像中の最も明るい物体は白色である」、「画像中に肌の色が含まれる」、「照明は白色である」などの仮定を満たすために必要なエネルギーを算出し、算出されたエネルギーが最適となるときの照明の分光分布および物体の表面反射率を推定値として出力する。
【００２９】
一方、第１の実施の形態における色仮定なし分光推定手段１３は、物体色または照明色に関する仮定を用いることなく、照明の分光分布および物体の表面反射率を算出する。ここで、照明の分光分布と物体の表面反射率の制約範囲を用いて照明の分光分布および物体の表面反射率を高精度に推定する。これら分光特性の推定には、式（１）を満たす十分な数の線形独立なスペクトル情報と、照明の分光分布および物体の表面反射率のモデル式から観測方程式を生成し、この観測方程式を解くことで照明の分光分布および物体の表面反射率を推定値として出力する。
【００３０】
色仮定なし分光推定手段１３の内部の構成を説明する。図２に示す色仮定なし分光推定手段１３は、観測方程式生成手段１３１と方程式計算手段１３２とを備える。
【００３１】
観測方程式生成手段１３１では、N次元のベクトルである、照明の分光分布Iおよび物体の表面反射率R のモデルを、照明の分光分布の基底ベクトルIpIbasisと物体の表面反射率の基底ベクトルRpRbasisそれぞれに対し、可変な重み係数apI,bpRを用いて式（２）のように生成し、モデル式（２）における基底ベクトルの各波長成分IpIbasis (λn)、RpRbasis (λn)と、式（１）を満たす最適化に十分な数の線型独立なスペクトル情報Em(λn)（m=1,…,M, n=1,…,N）が与えられるとき、式（３）のような観測方程式を生成する。式（１）の条件から方程式の数が未知数の数以上となり、かつ、線形独立なスペクトル情報で構成されるため、照明の分光分布I(λn)と物体の表面反射率R(λn)の最適解が得られることになる。
【００３２】
【数２】
0006015964.tif 000003
【００３３】
【数３】
0006015964.tif 000004
方程式計算手段１３２は、式（３）に示す観測方程式を解く。
【００３４】
第１の実施の形態の照明・反射率制約範囲保存メモリ１２は、照明分光分布や物体表面反射率の基底ベクトルを保持し、これら基底ベクトルを用いて、様々な制約方法によって制約された範囲に分布する照明分光分布および物体表面反射率を出力する。例えば、照明の分光分布の場合には、図８に示すように、日射モデルから計算される日射スペクトルの成分である、直達成分I1basis、散乱成分I2basisを基底ベクトルに用いることが考えられる。この日射モデルでは日時や場所などの容易に取得可能な照明条件によって、照明の分光分布の分布範囲が制約することが可能である。また、非特許文献１（新編色彩科学ハンドブック第2版 p69-72）に示すCIE昼光における、自然昼光の相対分光分布の平均ベクトルおよび主成分ベクトルを、照明の分光分布の基底ベクトルに用いることが考えられる。この主成分ベクトルの数を制限することや、非特許文献１に示される方法を用いて、主成分ベクトルの係数を相関色温度に基づいて計算することで、照明分光分布の分布範囲を制約できる。物体の表面反射率の場合には、被写体として計測される物体の反射率データベースの平均ベクトルおよび主成分ベクトルを基底ベクトルとして用いることが考えられ、この主成分ベクトルの数を制限することで物体の反射率の分布範囲を制約できる。例えば、図９に示すマクベスカラーチャートから生成した、平均ベクトルおよび主成分ベクトルを基底ベクトルとして用いることが考えられる。
【００３５】
照明・反射率制約範囲保存メモリ１２の内部の構成を説明する。図２に示す照明・反射率制約範囲保存メモリ１２は、日射スペクトル成分計算手段１２１と表面反射率基底ベクトル保存メモリ１２２と有限次元線型和生成手段１２３とを備える。
【００３６】
日射スペクトル成分計算手段１２１は、日射モデルなどを用いて算出される日射スペクトル成分などを算出し、各成分を、PI（< N）個の基底ベクトルIpIbasis ( pI =1,…, PI)として保持する。
【００３７】
表面反射率基底ベクトル保存メモリ１２２は、PR（<N）個の、物体表面反射率の基底ベクトルRpRbasis  ( pR =1,…, PR)を保持する。
【００３８】
有限次元線型和生成手段１２３では、照明分光分布の基底ベクトルとして保持された日射スペクトル成分IpIbasisと物体表面反射率の基底ベクトルRpRbasisそれぞれに対し、可変な重み係数apI,bpRを用いて、式（２）に示すように、線型和をとることで、照明の分光分布I、物体の表面反射率RをN次元のベクトルとして生成する。
【００３９】
図１０は照明分光分布の制約範囲のイメージを表した図である。例えば、Nが３、PIが２のとき、図１０に示すように、照明の分光分布の波長毎の強度を各次元の座標の値とし、一つの照明の分光分布が座標上の一点として表わされる照明の分光分布が取り得る分布空間では、基底ベクトルの線型和として表わされる照明の分光分布は平面となる。
【００４０】
すなわち、これは、照明の分光分布の取り得る範囲が制約されることを意味する。同様にして、式（２）の方法で、照明の分光分布I、物体の表面反射率Rを生成することで、照明の分光分布I、物体の表面反射率R がN次元の空間において取り得る範囲についても制約される。このようにして照明の分光分布Iと物体の表面反射率Rの制約範囲を算出する。
【００４１】
図３は、本発明を実施するための第１の実施の形態における分光画像処理方法の動作を示すフローチャートである。
【００４２】
マルチスペクトル画像撮影手段１１１は、シーンからマルチスペクトル情報を獲得し、マルチスペクトル画像を出力する（ステップＳ１０１）。
【００４３】
スペクトル選択手段１１２は、マルチスペクトル画像からスペクトルを選択し、十分な数の線型独立なスペクトル情報を出力する（ステップＳ１０２）。
【００４４】
日射スペクトル成分計算手段１２１は、自動で取得可能な照明条件を自動で取得する（ステップＳ１０３）。
【００４５】
日射スペクトル成分計算手段１２１は、日射スペクトル成分を計算し、出力する（ステップＳ１０４）。
【００４６】
限次元線型和生成手段１２３は、日射スペクトル成分と表面反射率基底ベクトル保存メモリ１２２に保持してある表面反射率の基底ベクトルをもとに、有照明の分光分布と物体の表面反射率の制約範囲をそれぞれ生成する（ステップＳ１０５）。
【００４７】
観測方程式生成手段１３１は、十分な数の線型独立なスペクトル情報と照明分光分布と物体表面反射率の制約を用いて、最適解の得られる観測方程式を生成する（ステップＳ１０６）。
【００４８】
方程式計算手段１３２は、最適解の得られる観測方程式を解き、物体の表面反射率および照明の分光分布を出力する（ステップＳ１０７）。
【００４９】
第１の実施の形態は、マルチスペクトル情報取得手段１１によってシーンから十分な数の線型独立なスペクトル情報を獲得し、照明・反射率制約範囲保存メモリ１２に保持される照明の分光分布および物体の表面反射率の制約範囲を得ることで、色仮定なし分光推定手段１３を用いて物体色や照明色の仮定を用いずに、照明の分光分布と物体の表面反射率を推定することを可能としている。そのため、特定の仮定を満たすようなシーンでなくとも、安定して照明の分光分布および物体の表面反射率を精度よく推定することが可能となる。
【００５０】
第１の実施の形態における分光画像処理方法は、日射スペクトル成分計算手段１２１を用いて日射スペクトル成分を計算することで、照明において狭い制約範囲を出力することを可能とする。照明において狭い範囲の制約を得ることで、観測方程式生成手段１３１は、推定するパラメタの数が減り計算の容易な観測方程式の生成が可能となり、その結果、方程式計算手段１３２は、高精度に推定された物体表面反射率および照明分光分布を出力することが可能となる。
【００５１】
[第２の実施の形態]
図４は、本発明を実施するための第２の実施の形態を示すブロック図である。図４にも示されるように、本発明を実施するための第２の実施の形態における分光画像処理方法では、色仮定なし分光推定手段２３が第１の実施の形態における色仮定なし分光推定手段１３と異なり、その他の構成要素については第１の実施の形態と同様である。第１の実施の形態と同様の構成要素については、図２と同一の符号を付し、詳細な説明を省略する。
【００５２】
第２の実施の形態では、色仮定なし分光推定手段２３は、照明・反射率推定値算出手段２３１と観測スペクトル推定値算出手段２３２と誤差算出手段２３３と誤差最適化手段２３４を備える。
【００５３】
照明・反射率推定値算出手段２３１は、式（２）のような照明の分光分布I(λn)と、物体の表面反射率R(λn)の制約範囲のもと、照明の分光分布および物体の表面反射率の重み係数に初期パラメタを与えることで、照明の分光特性および物体の表面反射率の推定値を算出する。
【００５４】
推定観測スペクトル推定値算出手段２３２は、照明・反射率推定値算出手段２３１で算出された照明の分光特性および物体の表面反射率の推定値を式（３）の右辺に代入し観測スペクトルの推定値を算出する。
【００５５】
誤差算出手段２３３は、観測方程式（３）の誤差を考慮し、推定観測スペクトル推定値算出手段２３２で算出された観測スペクトル推定値（式（３）の右辺）と、最適化に十分な数の色情報として与えられた複数のマルチスペクトルから得られる観測値Em(λn)を用いて、誤差を算出する。照明の分光分布、物体の表面反射率、観測スペクトルの計測誤差をそれぞれ、εI(λn), εR(λn), εL(λn)とすると、観測方程式（３）は式（４）のように表せる。これらの計測誤差の最適化を近似する方法として、例えば、式（５）に示すような最小二乗誤差エネルギーを定義し、このエネルギーを最小化する方法がある。
【００５６】
【数４】
0006015964.tif 000005
【００５７】
【数５】
0006015964.tif 000006
【００５８】
誤差最適化手段２３４は、例えばLevenberg-Marquardt法などの非線型最適化方法や、粗密探索法を用いて、誤差算出手段２３３で算出される誤差を最適化する。具体的には、誤差が最小となるように、照明・反射率推定値算出手段２３１に与えるパラメタの更新を繰り返し、最終的に最適なパラメタから算出される物体の表面反射率および照明の分光分布を最適解として出力する。
【００５９】
色仮定なし分光推定手段１３では、誤差最適化手段２３４で算出された照明の分光分布と物体の表面反射率の推定値を出力する。
【００６０】
図５は、本発明を実施するための第２の実施の形態における分光画像処理方法動作の一例を示すフローチャートである。第１の実施の形態と同様の動作については図３と同一の符号を付し、詳細な説明を省略する。
【００６１】
第２の実施の形態は、以下の点で第１の実施の形態と異なる。
【００６２】
照明・反射率推定値算出手段２３１は初期パラメタと照明の分光分布と物体の表面反射率の制約範囲をもとに、照明の分光分布と物体の表面反射率の推定値を算出する（ステップＳ２０６）。
【００６３】
観測スペクトル推定値算出手段２３２は照明の分光分布と物体の表面反射率の推定値をもとに、観測スペクトルの推定値を算出する（ステップＳ２０７）。
【００６４】
誤差算出手段２３３は、観測スペクトルの推定値と、観測された十分な数の線型独立なスペクトル情報の誤差を算出する（ステップＳ２０８）。
【００６５】
誤差最適化手段２３４は、誤差を用いて、誤差最適化手段ごとに設定された値を算出し、算出した値が予め規定された値より小さければステップＳ１０８に、小さくなければステップＳ２１０の動作に移動する（ステップＳ２０９）。
【００６６】
誤差最適化手段２３４は、誤差を最小化するように照明・反射率推定値算出手段２３１に与えるパラメタを更新して、Ｓ２０６の動作に戻る（ステップＳ２１０）。
【００６７】
[第３の実施の形態]
図６は、本発明を実施するための第３の実施の形態を示すブロック図である。図６にも示されるように、本発明を実施するための第３の実施の形態における分光画像処理方法では、マルチスペクトル情報取得手段２１が第１の実施の形態のマルチスペクトル情報取得手段１１と異なり、その他の構成要素については第１の実施の形態と同様である。第１の実施の形態と同様の構成要素については、図２と同一の符号を付し、詳細な説明を省略する。
【００６８】
第３の実施の形態では、マルチスペクトル情報取得手段２１は、マルチスペクトル複数回数計測手段２１１を備える。
【００６９】
マルチスペクトル複数回数計測手段２１１は、観測されたシーン中に十分な数の線形独立なスペクトル情報が含まれない場合に対応できるように、照明色が同一とみなせる範囲で、式（１）に示される条件を満たす十分な数の線形独立なスペクトル情報に達するまで異なるシーンを撮影し、後段の最適化に十分な数の線形独立な関係にある複数のスペクトル情報を出力する。
【００７０】
即ち、複数回数計測手段２１１は、例えば、シーン中に十分な線形独立な観測スペクトルが存在しないような均一シーンの撮影において、十分な線形独立な観測スペクトルが得られるまでシーン撮影を繰り返す（動画撮影も含む）際に使用することができる。
【００７１】
図７は、本発明を実施するための第３の実施の形態における分光画像処理方法動作の一例を示すフローチャートである。第１の実施の形態と同様の動作については図３と同一の符号を付し、詳細な説明を省略する。
【００７２】
第３の実施の形態は、以下の点で第１の実施の形態と異なる。
【００７３】
マルチスペクトル複数回数計測手段２１１は、シーンからマルチスペクトル情報を計測する（ステップＳ３０１）。
【００７４】
マルチスペクトル複数回数計測手段２１１は、シーンから計測されたマルチスペクトル情報が、十分な数の線型独立なマルチスペクトル情報である場合はステップＳ１０６に、そうでない場合はステップＳ３０１に移動する（ステップＳ３０２）。
【００７５】
第３の実施の形態における分光画像処理方法は、マルチスペクトル複数回計測手段２１１を備え、式（１）に示される条件を満たす十分な数の線形独立なスペクトル情報に達するまで異なるシーンを撮影することで、観測されたシーン中に十分な数の線形独立なスペクトル情報が含まれない場合でも、十分な数の線形独立なスペクトル情報を得るといった効果が得られる。
【００７６】
尚、上述した説明からも明らかなように、各部をハードウェアで構成することも可能であるが、コンピュータプログラムにより実現することも可能である。この場合、プログラムメモリに格納されているプログラムで動作するプロセッサによって、上述した各実施の形態と同様の機能、動作を実現させる。また、上述した実施の形態の一部の機能のみをコンピュータプログラムにより実現することも可能である。
【００７７】
また、上記の実施形態の一部又は全部は、以下の付記のようにも記載されうるが、以下には限られない。
【００７８】
（付記１）  シーンのマルチスペクトル画像から観測スペクトル情報を選択し、
照明の基底ベクトルと表面反射率の基底ベクトルを用いて、照明の分光分布の制約範囲と物体の表面反射率の制約範囲を生成し、
前記観測スペクトル情報と、前記照明の分光分布の制約範囲と、前記物体の表面反射率の制約範囲とを用いて、前記照明の基底ベクトルの重み係数と前記表面反射率の基底ベクトルの重み係数とをパラメタとする観測方程式を生成し、
前記観測方程式から前記シーンの照明の分光分布と物体の表面反射率とを算出する
分光画像処理方法。
【００７９】
（付記２）  前記照明の分光分布の制約範囲と前記物体の表面反射率の制約範囲のもと、前記観測方程式の重み係数に初期パラメタを与えて、前記観測スペクトルの推定値を算出し、
前記観測スペクトル情報と前記観測スペクトルの推定値の誤差に基づき、前記重み係数を更新し、
前記誤差が所定内となったときの前記観測スペクトルの推定値を構成する照明の分光分布と物体の表面反射率とを、前記シーンにおける照明の分光分布と物体の表面反射率として算出する
付記１に記載の分光画像処理方法。
【００８０】
（付記３）  前記観測スペクトル情報として、前記マルチスペクトル画像から線形独立なスペクトル情報を選択する
付記１又は付記２に記載の分光画像処理方法。
【００８１】
（付記４）  前記照明の基底ベクトルとして、日射スペクトルの直達成分ベクトルと散乱成分ベクトルとを用いる
特徴とする付記１から付記３のいずれかに記載の分光画像処理方法。
【００８２】
（付記５）  前記線形独立なスペクトル情報が前記照明の基底ベクトルの重み係数の個数と前記表面反射率の基底ベクトルの重み係数の個数の総和と同数かそれより多く観測できるまで複数回計測する
付記３又は付記４に記載の分光画像処理方法。
【００８３】
（付記６）  シーンのマルチスペクトル画像から観測スペクトル情報を選択するスペクトル情報選択手段と、
照明の基底ベクトルと表面反射率の基底ベクトルを用いて、照明の分光分布の制約範囲と物体の表面反射率の制約範囲を生成する制約範囲生成手段と、
前記観測スペクトル情報と、前記照明の分光分布の制約範囲と、前記物体の表面反射率の制約範囲とを用いて、前記照明の基底ベクトルの重み係数と前記表面反射率の基底ベクトルの重み係数とをパラメタとする観測方程式を生成する観測方程式生成手段と、
前記観測方程式から前記シーンの照明の分光分布と物体の表面反射率とを算出する算出手段と
を有する分光画像処理装置。
【００８４】
（付記７）  前記算出手段は、
前記照明の分光分布の制約範囲と前記物体の表面反射率の制約範囲のもと、前記観測方程式の重み係数に初期パラメタを与えて、前記観測スペクトルの推定値を算出する手段と、
前記観測スペクトル情報と前記観測スペクトルの推定値の誤差に基づき、前記重み係数を更新する手段と、
前記誤差が所定内となったときの前記観測スペクトルの推定値を構成する照明の分光分布と物体の表面反射率とを前記シーンにおける照明の分光分布と物体の表面反射率として算出する手段と
を有する付記６に記載の分光画像処理装置。
【００８５】
（付記８）  前記スペクトル情報選択手段は、前記マルチスペクトル画像から線形独立なスペクトル情報を選択する
付記６又は付記７に記載の分光画像処理装置。
【００８６】
（付記９）  前記スペクトル情報選択手段は、前記照明の基底ベクトルとして、日射スペクトルの直達成分ベクトルと散乱成分ベクトルとを用いる
特徴とする付記６から付記８のいずれかに記載の分光画像処理装置。
【００８７】
（付記１０）  前記スペクトル情報選択手段は、前記線形独立なスペクトル情報が前記照明の基底ベクトルの重み係数の個数と前記表面反射率の基底ベクトルの重み係数の個数の総和と同数かそれより多く観測できるまで複数回計測する
付記８又は付記９に記載の分光画像処理装置。
【００８８】
（付記１１）  コンピュータに、
シーンのマルチスペクトル画像から観測スペクトル情報を選択する処理と、
照明の基底ベクトルと表面反射率の基底ベクトルを用いて、照明の分光分布の制約範囲と物体の表面反射率の制約範囲を生成する処理と、
前記観測スペクトル情報と、前記照明の分光分布の制約範囲と、前記物体の表面反射率の制約範囲とを用いて、前記照明の基底ベクトルの重み係数と前記表面反射率の基底ベクトルの重み係数とをパラメタとする観測方程式を生成する処理と、
前記観測方程式から前記シーンの照明の分光分布と物体の表面反射率とを算出する処理と
を実行させるプログラム。
【００８９】
以上、これまで述べてきた各実施の形態は、本発明の好的な実施形態であり、上記実施の形態のみに本発明の範囲を限定するものではなく、本発明の要旨を逸脱しない範囲に置いて種々の変更を施した形態での実施が可能である。
【００９０】
本出願は、２０１１年１２月２８日に出願された日本出願特願２０１１−２８８８８６号を基礎とする優先権を主張し、その開示の全てをここに取り込む。
【符号の説明】
【００９１】
１          色情報取得手段
２          分光推定手段
３          照明・反射率主成分ベクトル保存メモリ
１１        マルチスペクトル情報取得手段
１２        色仮定なし分光推定手段
１３        照明・反射率制約範囲保存メモリ
２１        マルチスペクトル情報取得手段
２３        色仮定なし分光推定手段
１１１      マルチスペクトル画像撮影手段
１１２      スペクトル選択手段
１２１      日射スペクトル成分計算手段
１２２      表面反射率基底ベクトル保存メモリ
１２３      有限次元線型和生成手段
１３１      観測方程式生成手段
１３２      方程式計算手段
２１１      マルチスペクトル複数回計測手段
２３１      照明・反射率推定値算出手段
２３２      観測スペクトル推定値算出手段
２３３      誤差算出手段
２３４      誤差最適化手段
【図１】
0006015964.tif 000007
【図２】
0006015964.tif 000008
【図３】
0006015964.tif 000009
【図４】
0006015964.tif 000010
【図５】
0006015964.tif 000011
【図６】
0006015964.tif 000012
【図７】
0006015964.tif 000013
【図８】
0006015964.tif 000014
【図９】
0006015964.tif 000015
【図１０】
0006015964.tif 000016
【図１１】
0006015964.tif 000017