(57)【特許請求の範囲】
【請求項１】
共通する被写体を有し相異なる照明光下で撮像された複数の画像から、照明光色と共通被写体の色とを計算する色計算装置であって、
各画像における前記共通被写体の色から、前記共通被写体の色として存在しうる範囲である物体色範囲を求める物体色範囲計算部と、
前記物体色範囲計算部によって得られた、相異なる複数の照明光下の画像における前記共通被写体の物体色範囲を用いて、前記共通被写体の色を計算する物体色計算部と、
各画像から求められた前記共通被写体の物体色範囲を表す曲線上で、前記物体色計算部によって得られた前記共通被写体の色に最も近い点から、各画像の照明光色を計算する照明光色計算部と
を備える色計算装置。
【請求項２】
前記物体色計算部は、前記共通被写体を有する相異なる照明光下の画像において前記物体色範囲計算部から得られた物体色範囲に共通して含まれる領域の中心を前記共通被写体の色とする、請求項１に記載の色計算装置。
【請求項３】
前記物体色計算部は、前記共通被写体を有する相異なる照明光下の画像において前記物体色範囲計算部から得られた物体色範囲に対して、任意の２種類以上の物体色範囲で囲まれる領域の中心を前記共通被写体の色とする、請求項１に記載の色計算装置。
【請求項４】
前記物体色計算部は、前記共通被写体を有する相異なる照明光下の画像において前記物体色範囲計算部から得られる複数の物体色範囲の中から２種類ずつの組を作成し、各組から得られる共通領域の平均を前記共通被写体の色とする、請求項１に記載の色計算装置。
【請求項５】
前記照明光色計算部は、前記共通被写体が複数存在する場合に、前記物体色計算部によって得られた各共通物体の物体色計算結果から得られる照明光色の平均を照明光色とする、請求項１〜４のいずれか一つに記載の色計算装置。
【請求項６】
共通する被写体を有し相異なる照明光下で撮像された複数の画像から、照明光色と共通被写体の色とを計算する色計算装置が行う色計算方法であって、
各画像における前記共通被写体の色から、前記共通被写体の色として存在しうる範囲である物体色範囲を求める物体色範囲計算ステップと、
前記物体色範囲計算ステップによって得られた、相異なる複数の照明光下の画像における前記共通被写体の物体色範囲を用いて、前記共通被写体の色を計算する物体色計算ステップと、
各画像から求められた前記共通被写体の物体色範囲を表す曲線上で、前記物体色計算ステップによって得られた前記共通被写体の色に最も近い点から、各画像の照明光色を計算する照明光色計算ステップと
を有する色計算方法。
【請求項７】
共通する被写体を有し相異なる照明光下で撮像された複数の画像から、照明光色と共通被写体の色とを計算するコンピュータに対し、
各画像における前記共通被写体の色から、前記共通被写体の色として存在しうる範囲である物体色範囲を求める物体色範囲計算ステップと、
前記物体色範囲計算ステップによって得られた、相異なる複数の照明光下の画像における前記共通被写体の物体色範囲を用いて、前記共通被写体の色を計算する物体色計算ステップと、
各画像から求められた前記共通被写体の物体色範囲を表す曲線上で、前記物体色計算ステップによって得られた前記共通被写体の色に最も近い点から、各画像の照明光色を計算する照明光色計算ステップと
を実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像に写っている環境の照明光や物体の色を推定する技術に関する。
【背景技術】
【０００２】
画像における画素値（Ｒ，Ｇ，Ｂ）は、以下の式１に示されるように、照明光（Ｅ（λ）：分光分布）と物体の色（ρ（λ）：表面反射率）の両方の情報を含んでいる。そのため、入力機器によって入力された画像の画素値は、入力機器の分光感度特性と、分光分布と表面反射率との積分で表される。ＲＧＢそれぞれの分光感度特性は、［ｒ］（λ）、［ｇ］（λ）、［ｂ］（λ）として表す。なお、［ｒ］は、ｒの上に横棒（バー）が付されていることを表す。［ｇ］、［ｂ］も同様である。また、Ｒ、Ｇ、Ｂおよびρの右下に付与したインデックスｉは、それぞれシーン中のｉ番目の物体の画素値、表面反射率であることを示す。
【数１】
0005739790.tif 000002
【０００３】
照明光や物体の色を計算する関連技術では、照明光や物体の色に関する性質に基づいて計算する手法がある。例えば、画像中の鏡面反射成分を利用する手法や、画像中の無彩色領域を利用する手法や、シーン中の物体の色に関する仮定を利用する手法等がある。鏡面反射成分を利用する技術（非特許文献１）では、誘電体からの反射光成分のうち鏡面反射成分が照明光の色を表すことを利用して照明光の色を推定する。無彩色領域を利用する手法（特許文献１）は、無彩色の領域からの反射光の色が照明光の色を反映していることを利用している。具体的には、事前に様々な照明光の下での白色（紙、色票）の画像から無彩色領域の色変化を求めておき、画像中の無彩色領域に近い色を、事前に求めた無彩色領域の色の中から探索し、その平均を照明光の色として推定する。また、シーン中の物体の色の平均が灰色（無彩色）であると仮定することによって、シーン中の反射光の色を照明光色として推定する手法（非特許文献２）もある。
【０００４】
上述の各関連技術は、１枚の画像のみを用いることを前提としている。しかし、複数枚の画像を用いる手法もある。このような手法は、同一照明光下の異なるシーンの画像を複数枚用いる手法と、同一シーンの相異なる照明光下での画像を複数枚用いる手法とに分類される。前者は、１枚の画像を用いる手法をそのまま適用することにより実現できる。後者の具体例としては、非特許文献２に記載の技術を拡張した手法（特許文献２参照）がある。この技術は、各画像に共通する被写体が複数あることを前提としている。この技術では、各被写体の平均の色が、照明光の変化によらず黒体放射軌跡（照明光の色の範囲とみなせる）の近傍にある場合に、選択した被写体の平均が灰色になると判断する。そして、その際の平均値を照明光の色として推定する。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００２−２９０９８８号公報
【特許文献２】特開平９−３０４１８０号公報
【非特許文献】
【０００６】
【非特許文献１】G.J.Klinker, S.A.Shafer, T.Kanade, ”Using a color reflection model to a separate highlights from object color”, Proceedings of the First International Conference on Computer Vision, pp.145-150 (1987).
【非特許文献２】R.Gershon and A.D.Jepson,”The Computation of Color Constant Descriptors in Chromatic Images”, Color Research and Application, Vol.14, No.6, pp.325-334 (1989).
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、上述した関連手法には以下のような課題があった。
鏡面反射成分を利用する手法では、画像中に鏡面反射成分があることを前提としている。そのため、画像中に鏡面反射成分が無い場合には適用できないという問題がある。また、無彩色領域を利用する手法では、画像中に鏡面反射は不要であるが、推定精度が無彩色領域の有無に依存するという問題がある。この手法では、さらに、彩度の低い領域の色を照明光の色として誤推定するという問題がある。灰色仮説を利用する手法では、仮説が成立しない場合に、推定精度が画像中の色分布の影響を強く受けるという問題がある。
【０００８】
また、画像を複数枚用いる場合であっても、各画像に共有する被写体の色分布によっては、それらの平均の色が黒体放射軌跡の近傍に位置しないと照明光の色を正しく推定することができないという問題がある。
上記事情に鑑み、本発明は、共通の被写体を含む相異なる照明光下の画像から、共通被写体の色又は照明光の色を計算する技術を提供することを目的としている。
【課題を解決するための手段】
【０００９】
本発明の一態様は、共通する被写体を有し相異なる照明光下で撮像された複数の画像から、照明光色と共通被写体の色とを計算する色計算装置であって、各画像における前記共通被写体の色から、前記共通被写体の色として存在しうる範囲を求める物体色範囲計算部と、前記物体色範囲計算部によって得られた、相異なる複数の照明光下の画像における前記共通被写体の物体色範囲を用いて、前記共通被写体の色を計算する物体色計算部と、前記物体色計算部によって得られた前記共通被写体の色を用いて、相異なる照明光下の画像の照明光色を計算する照明光色計算部とを備える色計算装置である。
【００１０】
本発明の一態様は、上記の色計算装置であって、前記物体色計算部は、前記共通被写体を有する相異なる照明光下の画像において前記物体色範囲計算部から得られた物体色範囲に共通して含まれる領域の中心を前記共通被写体の色とする。
【００１１】
本発明の一態様は、上記の色計算装置であって、前記物体色計算部は、前記共通被写体を有する相異なる照明光下の画像において前記物体色範囲計算部から得られた物体色範囲に対して、任意の２種類以上の物体色範囲で囲まれる領域の中心を前記共通被写体の色とする。
【００１２】
本発明の一態様は、上記の色計算装置であって、前記物体色計算部は、前記共通被写体を有する相異なる照明光下の画像において前記物体色範囲計算部から得られる複数の物体色範囲の中から２種類ずつの組を作成し、各組から得られる共通領域の平均を前記共通被写体の色とする。
【００１３】
本発明の一態様は、上記の色計算装置であって、前記照明光色計算部は、前記共通被写体が複数存在する場合に、前記物体色計算部によって得られた各共通物体の物体色計算結果から得られる照明光色の平均を照明光色とする。
【００１４】
本発明の一態様は、共通する被写体を有し相異なる照明光下で撮像された複数の画像から、照明光色と共通被写体の色とを計算する色計算装置が行う色計算方法であって、各画像における前記共通被写体の色から、前記共通被写体の色として存在しうる範囲を求める物体色範囲計算ステップと、前記物体色範囲計算ステップによって得られた、相異なる複数の照明光下の画像における前記共通被写体の物体色範囲を用いて、前記共通被写体の色を計算する物体色計算ステップと、前記物体色計算ステップによって得られた前記共通被写体の色を用いて、相異なる照明光下の画像の照明光色を計算する照明光色計算ステップとを有する色計算方法である。
【００１５】
本発明の一態様は、上記の色計算装置であって、共通する被写体を有し相異なる照明光下で撮像された複数の画像から、照明光色と共通被写体の色とを計算するコンピュータに対し、各画像における前記共通被写体の色から、前記共通被写体の色として存在しうる範囲を求める物体色範囲計算ステップと、前記物体色範囲計算ステップによって得られた、相異なる複数の照明光下の画像における前記共通被写体の物体色範囲を用いて、前記共通被写体の色を計算する物体色計算ステップと、前記物体色計算ステップによって得られた前記共通被写体の色を用いて、相異なる照明光下の画像の照明光色を計算する照明光色計算ステップとを実行させるためのプログラムである。
【発明の効果】
【００１６】
本発明により、共通の被写体を含む相異なる照明光下の画像から、共通被写体の色又は照明光の色を計算することが可能となる。
【図面の簡単な説明】
【００１７】
【図１】本発明の一実施形態である色計算装置の機能構成を示すブロック図である。
【図２】色計算装置１０の処理の流れを示すフローチャートである。
【図３】物体色を求める処理の概略を示す図である。
【図４】共通被写体の色から照明光色を求める処理の概略を示す図である。
【図５】色計算装置１０の処理の流れの変形例を示すフローチャートである。
【発明を実施するための形態】
【００１８】
図１は、本発明の一実施形態である色計算装置の機能構成を示すブロック図である。色計算装置は、バスで接続されたＣＰＵ（Central Processing Unit）やメモリや補助記憶装置などを備える。色計算装置は、色計算プログラムを実行することによって、画像入力部１０１、データ蓄積部１０２、色情報変換部１０３、物体色範囲計算部１０４、物体色計算部１０５、照明光色計算部１０６を備える装置として機能する。なお、色計算装置の各機能の全て又は一部は、ＡＳＩＣ（Application Specific Integrated Circuit）やＰＬＤ（Programmable Logic Device）やＦＰＧＡ（Field Programmable Gate Array）等のハードウェアを用いて実現されても良い。色計算プログラムは、コンピュータ読み取り可能な記録媒体に記録されても良い。コンピュータ読み取り可能な記録媒体とは、例えばフレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置である。色計算プログラムは、電気通信回線を介して伝送されても良い。
【００１９】
画像入力部１０１は、照明光２２から出力された照明で照らされた被写体２１が撮像されたカラー画像（以下、単に「画像」という。）を色計算装置１０に入力する。画像入力部１０１は、画像を色計算装置１０に入力可能であればどのような構成であっても良い。例えば、画像入力部１０１は、撮像素子を備え、画像を撮像しても良い。例えば、画像入力部１０１は、スキャナとして構成され、画像をスキャンして色計算装置１０に入力しても良い。例えば、画像入力部１０１は、カメラによって撮像された画像を色計算装置１０に入力するためのインターフェースとして構成されても良い。
【００２０】
以下の説明では、画像は、Ｒ（赤）、Ｇ（緑）、Ｂ（青）の三成分から構成される。なお、画像は、ＨＳＶやＹＣｂＣｒ、ＣＩＥＬＡＢ等色空間の値を画素値とする画像であって良い。
被写体２１は、複数枚の画像に共通する１種類以上の物体である。照明光２２は、被写体２１を照射する光源であり、複数存在しても良い。
【００２１】
データ蓄積部１０２は、磁気ハードディスク装置や半導体記憶装置などの記憶装置を用いて構成される。データ蓄積部１０２は、画像入力部１０１によって入力された画像を記憶する。データ蓄積部１０２は、色情報変換部１０３、物体色範囲計算部１０４、物体色計算部１０５、照明光色計算部１０６の各処理で得られた画像や途中の処理結果を記憶しても良い。この場合、データ蓄積部１０２は、コンピュータにおけるメモリに相当する。
【００２２】
色情報変換部１０３は、画像入力部１０１によって入力された画像の色情報を、画像を撮像した装置の機器特性を排除した色情報に変換する。色情報変換部１０３は、変換後の画像をデータ蓄積部１０２に書き込む。色情報の変換に用いられるパラメータは、画像入力部１０１に応じて予め設定されている。以下の説明における色情報は、ＣＩＥ（国際照明委員会）で規定されたＸＹＺ三刺激値とする。ただし、色情報は、機器特性に依存しない色空間で表現されたものであればどのようなものでも構わない。
【００２３】
物体色範囲計算部１０４は、色情報変換部１０３によって色情報が変換された画像を用いて処理を行う。データ蓄積部１０２に格納されている画像の共通被写体領域の色情報に対し、被写体の色として存在しうる色範囲を求める。そして、物体色範囲計算部１０４は、求めた色範囲の情報（以下、「色範囲情報」という。）をデータ蓄積部１０２に書き込む。
【００２４】
物体色計算部１０５は、複数の色範囲情報に基づいて被写体２１の物体色を求め、求めた物体色の情報をデータ蓄積部１０２に書き込む。この処理において、物体色計算部１０５が用いる複数の色範囲情報は、同一の被写体２１に関する色範囲情報である。物体色計算部１０５が求める物体色は、色空間における座標値である。
照明光色計算部１０６は、物体色計算部１０５によって求められた物体色に基づいて照明光色を求める。照明光色計算部１０６は、求めた物体色をデータ蓄積部１０２に書き込む。
【００２５】
図２は、色計算装置１０の処理の流れを示すフローチャートである。まず、画像入力部１０１が画像を色計算装置１０に入力する（ステップＳ２０１）。データ蓄積部１０２は、入力された画像を格納する。色計算装置１０は、この処理を複数回繰り返すことによって同一の被写体２１が撮像された複数の画像をデータ蓄積部１０２に格納する。各画像は、それぞれ異なる照明光２２の下で撮像された画像である。具体的例は以下の通りである。屋外に設置された定点観測カメラや監視カメラ等で撮像された映像を、ある時間間隔でフレームとして抜き出して静止画として生成された画像。異なる日時に撮像した建設風景写真。なお、複数の画像それぞれにおいて、厳密に撮像位置やアングルが一致していなくともよい。
【００２６】
次に、色情報変換部１０３は、画像の色情報を変換し、機器特性を排除した画像に変換する（ステップＳ２０２）。機器特性とは、上述した式１に示す分光感度特性に対応する。ただし、このような特性に関する情報は入手できない場合も多い。その場合、式２及び式３に示す関係式に基づいてγ係数と３×３の行列（αｉ，ｊ）を求め、機器特性として用いても良い。色情報変換部１０３は、以下の式３を用いてＲＧＢ画素値からＸＹＺ三刺激値を求めることにより、色情報への変換を行う。
【数２】
0005739790.tif 000003
【数３】
0005739790.tif 000004
【００２７】
次に、物体色範囲計算部１０４は、共通被写体領域における色情報を用いて、共通被写体の色としてとりうる範囲を求める（ステップＳ２０３〜Ｓ２０６）。ｉ番目の照明光下における共通被写体領域Ｐの色情報（ＸＹＺ三刺激値）は以下の式４として表される。
【数４】
0005739790.tif 000005
【００２８】
式４の左辺は、ＸＹＺ三刺激値、即ち、共通被写体領域の色を入力機器に依存しない色情報に変換した結果である。式４の右辺のうち、第１項はＣＩＥ（国際照明委員会）で規定されたＸＹＺ三刺激値に対応した等色関数を表し、第２項はｉ番目の照明光の波長毎の成分を表し、第３項は共通被写体の表面反射率の波長毎の成分を表す。また、Ｎは可視光領域（例えば、４００〜７００ｎｍ）の分割数である。
【００２９】
式４において、照明光と、共通被写体領域の表面反射率とは、いずれも未知の値である。しかし、物体色範囲計算部１０４は、自然光の分光特性が黒体放射に近いことを利用して共通被写体の色範囲を求める。なお、黒体放射は、プランクの放射則にしたがって色温度Ｔと波長λの関数で式５のように表現される。
【数５】
0005739790.tif 000006
【００３０】
色温度Ｔ毎の分光特性（式５）を式４にあてはめることによって、照明光が色温度Ｔの黒体放射に近い分光特性を持つ場合の共通被写体の表面反射率（ρ（λ））が得られる。照明光の分光特性を既知とし、式４の右辺第１項及び第２項の積を３行Ｎ列からなる行列Ａとみなす。そして、擬似逆行列を両辺に作用させることによって、表面反射率（式６の左辺）が得られる。ただし、左上にｔを有するＡは、行列Ａの転置を表す。このような計算を様々な色温度Ｔの黒体放射に適用する。
【数６】
0005739790.tif 000007
【００３１】
上記の表面反射率に対し、等色関数を用いて色空間上の座標値に変換する。ここでは、色空間としてｘｙ色度空間の場合を説明するが、どのような色空間であっても構わない。
【数７】
0005739790.tif 000008
【００３２】
式７における左辺は、式６で得られた共通被写体の表面反射率のＸＹＺ三刺激値であり、これを用いてｘｙ色度を計算する。
【数８】
0005739790.tif 000009
【００３３】
このように、照明光が黒体放射の分光特性で近似できると仮定することによって、未知の照明光下の画像における共通被写体領域の色から、その被写体の色の範囲を求めることができる。物体色範囲計算部１０４は、このような処理を、同一被写体が撮像された複数の画像に対して繰り返し実行する。上述したように、物体色範囲計算部１０４の処理に用いられる各画像は、異なる照明光２２の下で撮像された画像である。色範囲は、色空間上の複数の点の集合であっても良いし、色空間上の点の集合の範囲（始点と終点）であっても良いし、これら点の集合を関数近似したものでも良い。以下の説明では、共通被写体の色の範囲を二次関数（ｙ＝ａ（ｘ＾２）＋ｂｘ＋ｃ）の形式で表現する。ｘおよびｙは、それぞれｘｙ色度の座標値である。また、“＾２”の記載は、その左にある値（上記二次関数の場合はｘ）の二乗を表す。
【００３４】
次に、物体色計算部１０５は、相異なる照明光２２の下で撮像された複数の画像から得られた共通被写体の色の範囲を用いて、この被写体（共通被写体）の色を求める（ステップＳ２０７）。図３は、物体色を求める処理の概略を示す図である。図３Ａは、３種類の照明光２２の下でそれぞれ撮像された同一の被写体領域の色から得られた被写体の色範囲をＬ１，Ｌ２，Ｌ３の曲線として示した図である。図示されるように、処理に求める色範囲の曲線Ｌ１〜Ｌ３の全てが１点で交わっている場合は、物体色計算部１０５は交点を共通被写体の色（ｘｙ色度）として求める。
【００３５】
しかし、処理に求める色範囲の曲線の全てが１点で交わるとは限らない。図３Ｂは、４種類の照明光２２の下でそれぞれ撮像された同一の被写体領域の色から得られた被写体の色範囲をＬ４〜Ｌ７の曲線として示した図である。図３Ｂに示すように、処理に求める色範囲の４つの曲線（Ｌ４〜Ｌ７）の全てが１点で交わらない場合には、物体色計算部１０５は、以下のように共通被写体の色を求める。なお、図３Ｂ〜図３Ｅでは、処理に求める色範囲の曲線が４つである場合を例として説明するが、処理に求める色範囲の曲線の数は４つに限定される必要は無い。また、データ蓄積部１０２に格納される全ての色範囲の曲線が処理に用いられても良いし、データ蓄積部１０２に格納される一部の色範囲の曲線が処理に用いられても良い。
【００３６】
図３Ｃは、４つの曲線（Ｌ４〜Ｌ７）で囲まれる領域Ａを示す図である。物体色計算部１０５は、領域Ａの重心（中心）を共通被写体の色（ｘｙ色度）として求めても良い。図３Ｄは、任意の２曲線で囲まれる領域の論理和で表される領域（領域Ｂ）を示す図である。物体色計算部１０５は、領域Ｂの重心（中心）を共通被写体の色（ｘｙ色度）として求めても良い。図３Ｅは、任意の２曲線の交点（以下、「色交点」といい、図中“ｘ”で示す。）の例を示す図である。物体色計算部１０５は、色交点の平均の位置を共通被写体の色（ｘｙ色度）として求めても良い。
【００３７】
照明光色計算部１０６は、ステップＳ２０７で得られた共通被写体の色に基づいて、各画像の照明光色を求める（ステップＳ２０８〜Ｓ２１１）。図４は、共通被写体の色から照明光色を求める処理の概略を示す図である。図４の曲線Ｌ８は、ある照明下の画像で得られた共通被写体の色範囲を示す曲線を表す。また、図４の点Ｐは、ステップＳ２０７で得られた共通被写体の色を表す。照明光色計算部１０６は、曲線Ｌ８上の点のうち最も点Ｐに近い点を、照明光の色（図中、×で示す）として求める。
【００３８】
本発明の色計算装置１０によれば、共通の被写体を含む複数種の照明光下の画像から、共通被写体の色を計算することが可能となる。また、共通被写体の色のみならず、照明光の色を計算することも可能である。そのため、任意の照明光下での画像への変換や、他の照明光環境で取得した画像との合成等も容易に行うことが可能となる。
【００３９】
＜変形例＞
（（式４の変形例））
式４のように直接的に表面反射率を求める代わりに、基底関数の係数を求めることによって間接的に表面反射率を求めても良い。基底関数は、事前に様々な物体の表面反射率を主成分分析して得られた基底であり、これらの線形結合で表面反射率を表現することができる。基底関数は、例えば以下の文献に開示されている。
【００４０】
J. Cohen, "Dependency of the spectral reflectance curves of the Munsell color chips", Psychonomic Science, vol.1, pp.369-370 (1964).
基底関数を用いる場合には、式４は以下の式９のように表される。
【数９】
0005739790.tif 000010
【００４１】
上式において、ｓｋｌはＭ種類（ｋ＝１，２，…，Ｍ）の基底関数における波長毎の成分（ｌ＝１，２，…，Ｎ）である。αｋは、基底関数の係数である。ｓｋｌとαｋとの積は、表面反射率を表す。式９の場合、右辺第１項から第３項までの積を３行Ｍ列からなる行列Ａ’とみなし、式４の場合と同等に擬似逆行列を両辺に作用させることで基底関数の係数が得られる。ここで得られた係数と基底関数から表面反射率を求め、式７及び式８にしたがって共通被写体の色を求めることができる。
【００４２】
（（フローチャートの変形例））
図５は、色計算装置１０の処理の流れの変形例を示すフローチャートである。変形例の処理は、相異なる照明光下の画像に共通する被写体が複数存在する場合に用いられる処理である。なお、図５において、図２と同じ処理には同じ符号を付す。
まず、画像入力部１０１は、相異なるｉ種類の照明光下でｊ種類の被写体を含む画像を色計算装置１０に入力する（ステップＳ２０１）。色情報変換部１０３は、各画像の色情報を変換する（ステップＳ２０２）。
【００４３】
次に、物体色範囲計算部１０４は、１つめの共通被写体の領域に対して（ステップＳ２１２）、相異なる照明光下の画像から被写体の色範囲を求める（ステップＳ２０３〜Ｓ２０７）。物体色計算部１０５は、照明光条件のみが異なる画像から得られた物体色範囲を用いて、１つめの共通被写体の色を計算する（ステップＳ２０７）。照明光色計算部１０６は、１つめの共通被写体の色に基づいて照明光色を計算する（ステップＳ２０８〜Ｓ２１１）。ここまでの処理は、ステップＳ２１２において変数ｊに１を代入して初期化する処理を除いて、上述した図２の処理フローと同じである。
次に、ステップＳ２０３〜Ｓ２１１までの処理を、２つめ以降の共通被写体に対して行う（ステップＳ２１３、Ｓ２１４）。照明光色計算部１０６は、同一照明光下の異なる被写体領域から得られた照明光の色を平均し、照明光の色を求める（ステップＳ２１５）。
【００４４】
以上、この発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。
【符号の説明】
【００４５】
１０…色計算装置，  １０１…画像入力部，  １０２…データ蓄積部，  １０３…色情報変換部，  １０４…物体色範囲計算部，  １０５…物体色計算部，  １０６…照明光色計算部，  ２１…被写体，  ２２…照明光
【図１】
0005739790.tif 000011
【図２】
0005739790.tif 000012
【図３】
0005739790.tif 000013
【図４】
0005739790.tif 000014
【図５】
0005739790.tif 000015