(57)【要約】【課題】観察環境に応じたカラーマッチングを行うには光源種を推定する必要があるが、そのために測光器を用いた測光を行うことは煩雑であった。【解決手段】まず、推定対象となる光源種Ａ〜Ｅごとに２つのカラーパッチからなるペアパッチを有する光源推定用画像９０２を印刷し、該光源推定用画像を光学フィルタ９０３を介して観察する。そして、２つのカラーパッチが等色として観察されるペアパッチを検出し、その対応する光源種Ｂを観察環境の光源種として推定する。このペアパッチは、推定対象となる光源種（光源Ａ〜Ｅ）の分光分布と光学フィルタ９０２の分光透過率に応じて、該光源種による観察環境下で等色として観察されるようなメタメリックペアである。【選択図】  図１０Ｂ
2009071550.tif 000002
【特許請求の範囲】
【請求項１】
観察環境の光源種を推定する光源推定方法であって、
推定対象となる光源種ごとに２つのカラーパッチからなるペアパッチを有する光源推定用画像を出力するパッチ出力ステップと、
前記光源推定用画像を光学フィルタを介して観察することにより、前記２つのカラーパッチが等色として観察されるペアパッチを検出し、該検出されたペアパッチの対応する光源種を観察環境の光源種として推定する推定ステップと、を有し、
前記ペアパッチは、前記推定対象となる光源種の分光分布と前記光学フィルタの分光透過率に応じて、該光源種による観察環境下で等色として観察されるように生成されていることを特徴とする光源推定方法。
【請求項２】
前記推定ステップにおいては、前記ペアパッチにおける少なくとも一方のカラーパッチを前記光学フィルタを介して観察することを特徴とする請求項１に記載の光源推定方法。
【請求項３】
前記推定ステップにおいては、前記ペアパッチにおけるそれぞれのカラーパッチを互いに異なる特性を有する前記光学フィルタを介して観察することを特徴とする請求項１に記載の光源推定方法。
【請求項４】
前記ペアパッチにおける一方のカラーパッチは、推定対象となる全ての光源種に対して同一色であることを特徴とする請求項３に記載の光源推定方法。
【請求項５】
前記推定ステップにおいては、推定対象となる全ての光源種に対して同一の前記光学フィルタを用いることを特徴とする請求項１乃至４のいずれか１項に記載の光源推定方法。
【請求項６】
前記推定ステップにおいては、推定対象となる光源種に応じて、前記光学フィルタを切り替えることを特徴とする請求項１乃至４のいずれか１項に記載の光源推定方法。
【請求項７】
前記パッチ出力ステップにおいては、前記光源推定用画像を記録用紙上に印刷出力することを特徴とする請求項１乃至６のいずれか１項に記載の光源推定方法。
【請求項８】
光学フィルタを介して観察されることによって観察環境の光源種を推定するために印刷出力される、複数のペアパッチからなる光源推定用画像を生成する画像生成方法であって、
前記光源推定用画像の印刷に用いられる色材の分光反射率を取得する分光反射率取得ステップと、
前記色材の分光反射率に対して前記光学フィルタの分光透過率を乗じて補正する補正ステップと、
推定対象となる全ての光源種のそれぞれについて、その分光分布と、前記色材の補正後の分光反射率に基づいて、当該光源種による観察環境下で等色として観察されるように、前記ペアパッチの色を算出するパッチ色の算出ステップと、
を有することを特徴とする画像生成方法。
【請求項９】
前記補正ステップにおいては、前記光学フィルタに対応すると予め定められた色材の分光反射率に対して、前記光学フィルタの分光透過率を乗じることを特徴とする請求項８に記載の画像生成方法。
【請求項１０】
前記色材は、少なくとも３種の色材であることを特徴とする請求項８または９に記載の画像生成方法。
【請求項１１】
前記色材は、ＣＭＹＫ印刷用とＲＧＢ印刷用であることを特徴とする請求項１０に記載の画像生成方法。
【請求項１２】
光学フィルタを介して観察されることによって観察環境の光源種を推定するために印刷出力される、複数のペアパッチからなる光源推定用画像を生成する画像生成装置であって、
前記光源推定用画像の印刷に用いられる色材の分光反射率を取得する分光反射率取得手段と、
前記色材の分光反射率に対して前記光学フィルタの分光透過率を乗じて補正する補正手段と、
推定対象となる全ての光源種のそれぞれについて、その分光分布と前記色材の分光反射率、および前記光学フィルタの分光透過率に基づいて、当該光源種による観察環境下で等色として観察されるように、前記ペアパッチの色を算出するパッチ色の算出手段と、
を有することを特徴とする画像生成装置。
【請求項１３】
コンピュータで実行されることにより、該コンピュータに請求項８乃至１１のいずれか１項に記載の画像生成方法を実行させることを特徴とするコンピュータプログラム。
【請求項１４】
請求項１３に記載のコンピュータプログラムを記録したことを特徴とするコンピュータ読取可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は光源推定方法ならびに画像生成方法および画像生成装置に関し、特にメタメリックペアによる等色判定を行う光源推定方法ならびに画像生成方法および画像生成装置に関する。
【背景技術】
【０００２】
近年、デスクトップ・パブリッシング（ＤＴＰ）技術やコンピュータグラフィクス（ＣＧ）技術の発展に伴い、出力物の色を管理するカラーマネージメント技術の重要性が高まりつつある。例えばＣＧでは、モニタ上のカラー画像を観察しつつＣＧの作成／編集作業等を行った後に、カラープリンタで入稿原稿の確認を行う場合がある。この場合、モニタ上のカラー画像とプリンタ出力画像とが知覚的に一致していることが望ましい。他方でＤＴＰに於いては、入稿先での本印刷を、入稿前にカラープリンタで出力して確認する。この場合、入稿先での本印刷画像とプリンタ出力画像とが知覚的に一致していることが望ましい。何れの場合に於いても、この要求を実現する為にはカラーマネージメント技術が必須となる。
【０００３】
このようなカラーマネージメントを実現する技術としては、International Color Consortium（ＩＣＣ）により定められたＩＣＣカラープロファイルを用いるカラーマネージメントシステム（ＣＭＳ）が一般的である。ＣＭＳでは、まずカラーマッチングを行う為のデバイスインディペンデントなProfile Connection Space（ＰＣＳ）を定義する。その上で、デバイス色空間からＰＣＳへの色変換を規定するソース側プロファイル、ならびにＰＣＳからデバイス色空間への色変換を規定するディスティネーション側プロファイルを用いて、カラーマネージメントを実現する。尚、ＰＣＳはハブ色空間と呼ばれることもある。色処理は、これらの２種のカラープロファイルに基づき、次の様な変換処理より実施される。すなわち、まずソース側プロファイルにより入力された画像の入力側デバイスに適合したデバイス色空間での色信号値を、ＰＣＳでの色信号値に変換する。その後、さらにディスティネーション側プロファイルにより、出力側デバイスに適合したデバイス色空間での色信号に変換する。
【０００４】
このようにＣＭＳによれば、ＣＧにて使用されるモニタ−プリンタシステム、あるいはＤＴＰにて使用されるプルーフシステムなど幅広く柔軟に対応できる。例えば上述したＣＧやＤＴＰでは、以下のようにソース側およびディスティネーション側のプロファイルを指定することにより、所望の画像とプリンタ出力画像との知覚的一致を図ることが可能である。すなわちＣＧでは、ソース側プロファイルにモニタの特性を記したプカラーロファイルを指定すると共に、ディスティネーション側プロファイルにプリンタの特性を記したカラープロファイルを指定する。またＤＴＰでは、ソース側プロファイルに入稿先の印刷機の特性を記したプカラーロファイルを指定すると共に、ディスティネーション側プロファイルにプリンタの特性を記したカラープロファイルを指定する。
【０００５】
但し、より正確に知覚的一致を図るためには、画像の観察環境として光源の特性を正しく把握して、プロファイルを適切に作成しなければならない。これは、印刷物を観察した際の色は、光源の光源色と印刷物の物体色とから、次式の様に求まる為である。
【０００６】
2009071550.tif 000003
【０００７】
ここでＲ(λ)は印刷物の分光反射率、Ｓ(λ)は光源の分光分布、ｘ(λ)、ｙ(λ)、ｚ(λ)は２度視野のＸＹＺ等色関数であり、上式から明らかな様に、光源の分光分布により色は変化する。従って、近似的にでもより正確に光源の特性を推定する必要がある。
【０００８】
さらに昨今ではより進んだカラーマネージメント技術として、International Commission on Illumination（ＣＩＥ）で標準化された、CIECAM02等のカラーアピアランスモデル（ＣＡＭ）が提案されている。このモデルでは、心理物理量を求める為に観察環境の情報が入力値として必要であり、今後はこのような背景からも、光源の特性を推定する必要が増してくると考えられる。
【０００９】
光源の特性を推定する手段としては、センサーを備えた測光器を用いるものが一般的である。更に、特性を推定する際には光源の波形情報が非常に重要である。例えば、同じ５０００Ｋの色温度を有するＦ１０光源とＦ８光源とで印刷物の色の見えを比較しても、波形が異なることに起因して色の見えが異なることが度々ある。従って、測光器としては分光分布を測定できるものが望ましい。
【発明の開示】
【発明が解決しようとする課題】
【００１０】
しかしながら、上述したような従来の測光器を用いる場合、一般にパーソナルコンピュータと接続して測光値を記録する必要があり、システム構成としての手軽さに欠ける。また測光器の中には、ＸＹＺ値やＸＹ値しか出力できないものもある。これは、このような型の測光器では、赤色／緑色／青色を有する光学フィルタ等を用いて大まかに数色程度への分光を行って検知することに起因する。このような測光器は、回折格子等を用いる分光測光器と比較して安価に構成できるという利点があるが、上述した様に分光情報を出力することができないため、光源の波形特性まで推定することは困難である。
【００１１】
光源の波形特性までを簡易に推定する手段としては、メタメリックペアを用いる方法が考えられる。しかしながら、メタメリックペアの生成はインクの特性に依存するため、等色時と非等色時との色差を大きくすることができず、条件等色しているかの否かの判定が困難となる場合がある。この例を図２６に示す。同図はあるメタメリックペアの分光反射率の比較例を示したものであるが、視覚の感度が高い波長域では、分光反射率の関数形状（波形）が比較的似た形状となっており、判別が困難であることが分かる。
【００１２】
本発明は上述した課題を解決するためになされたものであり、観察環境の光源種を、測光器を用いることなく容易に推定可能な光源推定方法ならびに画像生成方法および画像生成装置を提供することを目的とする。
【課題を解決するための手段】
【００１３】
上記目的を達成するための一手法として、本発明の光源推定方法は以下の工程を備える。
【００１４】
すなわち、観察環境の光源種を推定する光源推定方法であって、推定対象となる光源種ごとに２つのカラーパッチからなるペアパッチを有する光源推定用画像を出力するパッチ出力ステップと、前記光源推定用画像を光学フィルタを介して観察することにより、前記２つのカラーパッチが等色として観察されるペアパッチを検出し、該検出されたペアパッチの対応する光源種を観察環境の光源種として推定する推定ステップと、を有し、前記ペアパッチは、前記推定対象となる光源種の分光分布と前記光学フィルタの分光透過率に応じて、該光源種による観察環境下で等色として観察されるように生成されていることを特徴とする。
【００１５】
また、本発明の画像生成方法は以下の工程を備える。
【００１６】
すなわち、光学フィルタを介して観察されることによって観察環境の光源種を推定するために印刷出力される、複数のペアパッチからなる光源推定用画像を生成する画像生成方法であって、前記光源推定用画像の印刷に用いられる色材の分光反射率を取得する分光反射率取得ステップと、前記色材の分光反射率に対して前記光学フィルタの分光透過率を乗じて補正する補正ステップと、推定対象となる全ての光源種のそれぞれについて、その分光分布と、前記色材の補正後の分光反射率に基づいて、当該光源種による観察環境下で等色として観察されるように、前記ペアパッチの色を算出するパッチ色の算出ステップと、を有することを特徴とする。
【発明の効果】
【００１７】
以上の構成からなる本発明によれば、光源推定用画像を光学フィルタを介して観察することにより、観察環境の光源種を、測光器を用いることなく容易に推定することができる。
【発明を実施するための最良の形態】
【００１８】
以下、添付の図面を参照して、本発明をその好適な実施形態に基づいて詳細に説明する。なお、以下の実施形態において示す構成は一例に過ぎず、本発明は図示された構成に限定されるものではない。
【００１９】
＜第１実施形態＞
●装置構成
図１は、本実施形態における画像処理装置のシステム構成を示すブロック図である。同図において、１０１はＣＰＵ、１０２はメインメモリ、１０３はＳＣＳＩインタフェース、１０４はネットワークインタフェース、１０５はＨＤＤである。また、１０６はグラフィックアクセラレータ、１０７はカラーモニタ、１０８はＵＳＢコントローラ、１０９はカラープリンタ、１１０はキーボード／マウスコントローラ、１１１はキーボード、１１２はマウスである。そして、１１３はローカルエリアネットワーク、１１４はＰＣＩバスである。
【００２０】
以下、上記構成からなる本実施形態の画像処理装置における画像出力動作について説明する。
【００２１】
まず、ＨＤＤ１０５に格納されている画像処理アプリケーションならびにソフトウェア・ラスタイメージプロセッサ（以下、Ｓｏｆｔ・Ｒｉｐ）が、ＣＰＵ１０１からの指令により起動される。続いて、画像処理アプリケーション内の処理にしたがって、ＨＤＤ１０５に格納されている画像データが、ＣＰＵ１０１からの指令に基づきＳＣＳＩインタフェース１０３を介して、ＰＣＩバス１１４経由でメインメモリ１０２に転送される。また、ＬＡＮに接続されているサーバに格納されている画像データあるいはインターネット上の画像データは、ＣＰＵ１０１からの指令によりネットワークインタフェース１０４を介して、ＰＣＩバス１１４経由でメインメモリ１０２に転送される。
【００２２】
メインメモリ１０２に保持されている画像データは、ＣＰＵ１０１からの指令によりＰＣＩバス１１４経由でグラフィックアクセラレータ１０６に転送される。グラフィックアクセラレータ１０６は、画像データをＤ／Ａ変換した後にディスプレイケーブルを通じてカラーモニタ１０７に送信し、カラーモニタ１０７上に画像データが表示される。ここで、ユーザが画像処理アプリケーションにてメインメモリ１０２に保持されている画像をプリンタ１０９から出力するよう指令すると、画像処理アプリケーションはＳｏｆｔ・Ｒｉｐを起動し、画像をＳｏｆｔ・Ｒｉｐへ転送する。Ｓｏｆｔ・Ｒｉｐは、後述する画像処理動作に基づき、ＲＧＢ２４ビット画像をＣＭＹＫ３２ビット画像に変換した後、ＵＳＢコントローラ１０８を介して前記ＣＭＹＫ画像データをプリンタ１０９へ送信する。以上の一連の動作の結果として、プリンタ１０９よりＣＭＹＫ画像が印字される。
【００２３】
●画像出力処理（ＲＧＢ→ＣＭＹＫ）
上記構成における画像処理アプリケーション動作のうち、Ｓｏｆｔ・ＲｉｐによるＲＧＢ２４ビット画像からＣＭＹＫ３２ビット画像への変換について、図２のフローチャートを用いて説明する。
【００２４】
まずステップＳ２０１にて、適当なソースカラープロファイルならびにディスティネーションカラープロファイルをＨＤＤ１０５より読み込む。ここで、ソースカラープロファイルは、デバイスディペンデントな色空間（ＲＧＢ色空間、ＣＭＹＫ色空間等）からデバイスインディペンデントな色空間（ハブ色空間）への色変換特性が記述されたファイルである。ハブ色空間としては一般に、Ｌ*ａ*ｂ*色空間、ＸＹＺ色空間等が用いられる。また変換特性は、例えばＲＧＢからＸＹＺへの３×３変換マトリクス、あるいは変換ルックアップテーブル（ＬＵＴ）として記述される。他方、ディスティネーションカラープロファイルは、ハブ色空間からデバイスディペンデントな色空間（ＲＧＢ色空間、ＣＭＹＫ色空間等）への色変換特性が記述されたファイルである。この変換特性は、例えばＸＹＺとＣＭＹＫとの対応関係を示すＬＵＴとして記述される。
【００２５】
次にステップＳ２０２にてＳｏｆｔ・Ｒｉｐは、ソースカラープロファイルに記述された変換特性に基づき、画像処理アプリケーションより転送されたＲＧＢ２４ビット画像をＬ*ａ*ｂ*２４ビット画像へ変換する。変換後の画像はメインメモリ１０２もしくはＨＤＤ１０５に格納する。続いてステップＳ２０３にてＳｏｆｔ・Ｒｉｐは、ディスティネーションカラープロファイルに記述された変換特性に基づき、Ｌ*ａ*ｂ*２４ビット画像をＣＭＹＫ３２ビット画像へ変換し、メインメモリ１０２もしくはＨＤＤ１０５に格納する。そしてステップＳ２０４にて、メインメモリ１０２もしくはＨＤＤ１０５に格納されたＣＭＹＫ３２ビット画像を、プリンタドライバに転送する。
【００２６】
●デスティネーションプロファイル作成処理
上記構成における画像処理装置から画像を出力するに際しては、予め、適当なデスティネーションプロファイルを作成しておく必要がある。以下、本実施形態におけるデスティネーションプロファイルの作成処理について、図３のフローチャートを用いて説明する。この処理は、ＨＤＤ１０５に格納されているプロファイル作成ソフトウェアが、ＣＰＵ１０１からの指令により起動されことによって開始される。
【００２７】
まずステップＳ３０１にて、本実施形態の特徴である光源推定用のパッチ画像をプリンタ１０９より記録用紙上に印刷出力する。本実施形態では、このパッチ出力処理によって印刷された光源推定用パッチ画像をユーザが観察することによって、現在の光源を推定する。なお、この光源推定用パッチの生成方法の詳細については、図１１のフローチャートを用いて後述する。
【００２８】
ここで、印刷された光源推定用パッチ画像の例およびその観察方法について、図９および図１０Ａ，図１０Ｂを用いて詳細に説明する。
【００２９】
図９において、９０１は光源、９０２は光源推定用パッチ画像、９０３は光学フィルタ、９０４は観察者である。光源推定用パッチ画像９０２には、図１０Ａに示すように２行のパッチが印刷されており、パッチ数は判別対象の光源数×２である。本実施形態では、判別対象の光源数を５種として説明する。光学フィルタ９０３は図１０Ｂに示すように、光源推定用パッチ画像の下のパッチ行のみを覆うように配置される。この様な配置により、光源推定用パッチ画像９０２における上の行のパッチは、光源からの光線を通常の様に反射して観察者９０４に観察される。これに対し、下の行のパッチは、一度光学フィルタ９０３を透過した光線を反射し、さらにこの反射光がもう一度光学フィルタ９０３を透過した光線が、観察者９０４に観察される。
【００３０】
ここで図１０Ａに示すパッチ画像は、図１０Ｂにおける光源Ｂのように、観察環境における光源（この場合光源Ｂ）に相当する上下（ペア）のパッチが等色、もしくは近しい色として観察されるように生成されている。すなわち、図１０Ａに示すパッチ画像の下段に対して光学フィルタ９０３を配置することにより、ペアのパッチについて、等色時と非等色時との色差を大きくすることができ、ユーザは最も互いの色が近く見えたペアを容易に選択することができる。
【００３１】
図２に戻り、次にステップＳ３０２では、ユーザに対し、ステップＳ３０１で印刷された光源推定用パッチ画像の観察結果に基づき、適切な光源を選択するように要求する。図４に、この際のユーザインタフェース例を示す。図４に示すユーザインタフェース４０１によればすなわち、光学フィルタ９０３をパッチ下段に載置した図１０Ｂにおいて最も近い色に見えたペアパッチに相当する光源を、選択リスト４０２から選択するように要求する。図１０Ｂに示す例では光源Ｂが等色しているため、ユーザは選択リスト４０２から光源Ｂを選択し、ＯＫボタン４０３を押下する。これにより、本実施形態のデスティネーションプロファイルを作成する際に参照すべき光源が選択される。
【００３２】
次にステップＳ３０３では、ステップＳ３０２でユーザにより選択された光源に基づき、その光源下で観察する場合の印刷色再現特性データを算出する。この処理により、ＣＭＹＫ値の観察環境下における、Ｌ*ａ*ｂ*値との対応関係が算出される。なお、この算出手順の詳細については、図５のフローチャートを用いて後述する。
【００３３】
そしてステップＳ３０４では、ステップＳ３０３で算出された色再現特性データに基づいてＩＣＣプロファイルを生成する。ここで、ＩＣＣプロファイルはＬＵＴに基づくデータ構造を持ち、所定のＬ*ａ*ｂ*値に対して出力すべきＣＭＹＫ値を記憶するものである。このＩＣＣプロファイル算出手順の詳細については、図７のフローチャートを用いて後述する。
【００３４】
このようにして作成されたＩＣＣプロファイルが、現在の観察環境下における光源を反映した最適なデスティネーションプロファイルである。
【００３５】
●印刷色再現特性データ算出処理
以下、上述したステップＳ３０３における印刷色再現特性データの算出処理の詳細について、図５のフローチャートを用いて説明する。
【００３６】
まずステップＳ５０１にて、ＣＭＹＫ値に基づく印刷パッチの分光反射率を示す分光反射率データを取得する。この分光反射率データは、予め所定のパッチを測色することによって得られているとする。
【００３７】
そしてステップＳ５０２において、ステップＳ３０２で選択された光源について、その相対分光分布Ｓ(λ)をＨＤＤ１０５より読み出すことにより取得する。例えば、選択された光源がＦ８光源であった場合には、相対分光分布Ｓ(λ)は、図１３にて実線で表した関数形状として得られる。
【００３８】
ステップＳ５０３では、ステップＳ５０２で取得した相対分光分布Ｓ(λ)に基づいて、完全白色板に相当するＸＹＺ値であるＸw，Ｙw，Ｚwを下式により算出する。ここで、完全白色板の分光反射率は全波長において１としているため、分光反射率の関数は式中に表記されていない。また、ｘ(λ)，ｙ(λ)，ｚ(λ)は２度視野ＸＹＺ等色関数であり、図６に示す関数形状を呈する。
【００３９】
2009071550.tif 000004
【００４０】
ステップＳ５０４では、所定のＣＭＹＫパッチ色に対応する分光反射率Ｒi(λ)を取得する。ここでインデックスｉは、パッチの順序に対応するものである。
【００４１】
そしてステップＳ５０５において、取得した分光反射率Ｒi(λ)、相対分光分布Ｓ(λ)、並びにＸＹＺ等色関数に基づき、選択された光源下で所定のパッチを観察した際のＸＹＺ値を下式により算出する。
【００４２】
2009071550.tif 000005
【００４３】
ステップＳ５０６では、ステップＳ５０５で算出したＸＹＺ値と、完全白色版のＸＹＺ値とに基づき、下式の様にＬ*ａ*ｂ*値を算出する。尚、下式において、ｘr＝Ｘi／Ｘw、ｙr＝Ｙi／Ｙw、ｚr＝Ｚi／Ｚwである。
【００４４】
2009071550.tif 000006
【００４５】
さらに、算出したＬ*ａ*ｂ*値を、所定のパッチに対する値Ｌi，ｂi，ａiとして記憶する。
【００４６】
ステップＳ５０７では、総てのＣＭＹＫパッチ色に対してＬ*ａ*ｂ*値算出を行ったか否かを判定し、行っていればステップＳ５０８へ進み、行っていなければステップＳ５０４へ戻り、次のＣＭＹＫパッチ色に対する処理を行う。
【００４７】
ステップＳ５０８では、上記ステップＳ５０４〜Ｓ５０７で算出したＬ*ａ*ｂ*値を、選択された光源における印刷の色再現特性データとして、メインメモリ１０２に記憶する。
【００４８】
以上の処理により、選択された光源下においてプリンタ印刷を観察する場合の、色再現特性データが算出される。
【００４９】
●ＩＣＣプロファイル作成処理
以下、上述したステップＳ３０４におけるＩＣＣプロファイル生成処理について、図７のフローチャートを用いて詳細に説明する。
【００５０】
まずステップＳ７０１において、ステップＳ３０３で生成した色再現特性データをメインメモリ１０２より取得する。
【００５１】
そしてステップＳ７０２では、所定の順序に基づき、ＩＣＣプロファイルのＬＵＴのテーブル座標に相当するＬ*ａ*ｂ*値を取得する。
【００５２】
ステップＳ７０３では、ステップＳ７０２で取得したＬ*ａ*ｂ*値で示される色に対してデバイス色域の内外判定を行い、色域外であれば色域表面へ写像する。この写像は図８に示す様に、Ｌ*軸上の一点Ｆを焦点とする方式の写像で実現され、デバイス色域８０１の外に位置する色Ｃは、夫々デバイス色域表面の色Ｃmへと写像される。尚、Ｌ*ａ*ｂ*値で示される色がデバイス色域内に位置する場合には、値の変換は行わない。また、デバイス色域は、ステップＳ７０１で取得した色再現特性データに基づき、Convex Hull技法等を用いて算出される。
【００５３】
ステップＳ７０４では、ステップＳ７０１で取得した色再現特性データに基づき、ステップＳ７０５で取得したＬ*ａ*ｂ*値をＣＭＹＫ値に変換し、メインメモリ１０２に記憶する。この変換には例えば、ＣＭＹＫ値からＬ*ａ*ｂ*値への変換に基づいて、所望のＬ*ａ*ｂ*値を出力するＣＭＹＫ値を探索する、という探索技法を用いれば良い。
【００５４】
ステップＳ７０５では、ＬＵＴの全てのテーブル座標に対してＣＭＹＫ値算出を行ったか否かを判定し、行っていればステップＳ７０６へ進み、行っていなければステップＳ７０２へ戻って次のテーブル座標の処理を行う。
【００５５】
ステップＳ７０６では、上記ステップＳ７０２〜Ｓ７０５で算出したＣＭＹＫ値に基づいてＩＣＣプロファイルフォーマットに従ったファイルを生成し、ディスティネーションプロファイルとしてＨＤＤ１０５に記憶する。
【００５６】
●光源推定用パッチ画像作成処理
以下、上述したステップＳ３０１で印刷される光源推定用のパッチ画像について、その作成処理を図１１のフローチャート等を用いて詳細に説明する。
【００５７】
まずステップＳ１１０１では、パッチ画像作成のみに使用する、２種の色材料ＬＵＴを生成する。一方は、Ｒ／Ｇ／Ｂの３種のインクを使用するＬＵＴ（以下、ＲＧＢ・ＬＵＴ）、他方はＣ／Ｍ／Ｙ／Ｋの４種のインクを使用するＬＵＴ（以下、ＣＭＹＫ・ＬＵＴ）である。両ＬＵＴともにＣ'Ｍ'Ｙ'値に対応するテーブルであり、テーブルを参照するためのＣ'Ｍ'Ｙ'値の範囲は[０，２５５]である。また、テーブルに記述されているＲＧＢ各インクの値の範囲は同様に[０，２５５]であり、ＣＭＹＫ各インクの値の範囲も同様に[０，２５５]である。ここで、テーブルを参照する際のＣ'Ｍ'Ｙ'のステップを以下に示す。
【００５８】
Ｃ'＝{0,32,64,96,128,160,192,224,255}
Ｍ'＝{0,32,64,96,128,160,192,224,255}
Ｙ'＝{0,32,64,96,128,160,192,224,255}
ＲＧＢ・ＬＵＴに記述されるＲＧＢインク値は、Ｃ'Ｍ'Ｙ'値より以下のように算出される。
【００５９】
Ｒ＝１−Ｃ'，Ｇ＝１−Ｍ'，Ｂ＝１−Ｙ'
また、ＣＭＹＫ・ＬＵＴに記述されるＣＭＹＫインク値は、Ｃ'Ｍ'Ｙ'値より以下のように算出される。
【００６０】
Ｋ＝Min(Ｃ'，Ｍ'，Ｙ')、
Ｃ＝Ｃ'−Ｋ，Ｍ＝Ｍ'−Ｋ，Ｙ＝Ｙ'−Ｋ
ステップＳ１１０２では、作成した２種のＬＵＴそれぞれについて、テーブルに記載された色材量に基づく７２９色のパッチを印刷し、分光反射率を測定してメインメモリ１０２に記憶する。ここで、ＣＭＹＫ・ＬＵＴに対する分光反射率をＲ1i(λ)、ＲＧＢ・ＬＵＴに対する分光反射率をＲ2i(λ)と表す。なお、インデックスｉは夫々のパッチ順序に対応するものである。ここで、この分光反射率取得処理において取得された分光反射率の例を図１２Ａ，図１２Ｂに示す。図１２Ａは、ＣＭＹＫ・ＬＵＴにおけるグレイ色の分光反射率特性であり、図１２Ｂは、ＲＧＢ・ＬＵＴにおけるグレイ色の分光反射率特性である。
【００６１】
ステップＳ１１０３では、ＲＧＢ・ＬＵＴに対応する分光反射率Ｒ2i(λ)を、光学フィルタ９０３の分光透過率Ｔ(λ)の２乗を乗じることにより補正したＲ2'i(λ)を算出する。ここで分光透過率の２乗を乗ずるのは、光源からの光線が光学フィルタ９０３を透過した後、印刷物で反射された後で更にもう一度、光学フィルタ９０３を透過することを考慮するためである。
【００６２】
なお、本実施形態の光学フィルタ９０３は予め与えられたものであるため、分光透過率Ｔ(λ)もこの透過率乗算処理の時点で既に定まった値を有する。使用する光学フィルタ９０３を選択する際には、その分光透過率特性について次の注意が必要である。
【００６３】
本実施形態での光学フィルタ９０３に対する要求を簡単に説明すると、分光分布の関数形状が異なる波長域、すなわち分光分布の特徴を示す波長域において、透過率を遮断するような特性を有することが求められる。これは次の理由による。本実施形態の光学フィルタ９０３は、光源が変化した際のカラーパッチ色の変化を、光学フィルタ９０３を介さない場合のそれと差を持たせるようにする、といった効果を担う。これにより、分光分布の特徴を容易に検出することができる。さらに光学フィルタ９０３によって遮断する波長域は、視覚が感度を有する波長域でなければならないという条件を持つ。
【００６４】
ここで参考として、図１３Ａに、判別対象の光源のうち、Ｆ７，Ｆ８，Ｆ１０の三種の標準光源の相対分光分布を、それぞれ点線、実線、一点破線で示す。また図１３Ｂに、図１３Ａに示したＦ８光源とＦ７光源、および図６に示した等色関数とを同時に比較するグラフを示す。これらの関係を考慮して本実施形態では、光学フィルタ９０３として図１４に示す分光透過率特性を有するフィルタを使用する。光学フィルタ９０３が図１４に示す特性を有することにより、視覚が感度を有する波長域の全域において透過率が低くなるために分光分布の差が小さくなり、光源が変化した際のパッチ色の僅かな変化を際立たせることができる。これにより、ペアパッチが等色であるか否かを容易に判別することができる。
【００６５】
但し本発明における光学フィルタとしては、図１４のように急峻な特性を持つものに限らず、一般のカラーフィルタの様にブロードな特性を持つものも、光源を変化させた場合の色の変化に差をつける効果があるため、当然ながら適用可能である。このブロードな光学フィルタの適用可能性については、通常印刷でのメタメリズムを考え合わせれば容易に推測できる。
【００６６】
図１１に戻り、ステップＳ１１０４では、判別対象であるＮ種の光源の相対分光分布Ｓi(λ)をＨＤＤ１０５より取得する。なお、相対分光分布Ｓi(λ)の一例は、図１３Ａに示した通りである。
【００６７】
ステップＳ１１０５では、分光反射率Ｒ1i(λ)と補正後の分光反射率Ｒ2'i(λ)、並びに相対分光分布Ｓj(λ)とから、ｊ番目の光源下で観察した際のｉ番パッチのＸＹＺ値を、下式により算出する。まず、ＣＭＹＫ・ＬＵＴに対するＸＹＺ値であるＸ1ji，Ｙ1ji，Ｚ1jiは以下のようになる。
【００６８】
2009071550.tif 000007
【００６９】
また、ＲＧＢ・ＬＵＴに対するＸＹＺ値であるＸ2ji，Ｙ2ji，Ｚ2jiは以下のようになる。
【００７０】
2009071550.tif 000008
【００７１】
以下、ｊ番目の光源に基づきＣＭＹＫ・ＬＵＴをＸＹＺ値に変換したものをＬＵＴ1j、ＲＧＢ・ＬＵＴをＸＹＺ値に変換したものをＬＵＴ2jと称する。算出したＬＵＴは、メインメモリ１０２に記憶する。
【００７２】
ステップＳ１１０６では、ステップＳ１１０５で作成した２×Ｎ種のＸＹＺ値ＬＵＴに基づき、メタメリックペアの印刷色を算出する。この算出手順の詳細については、図１５のフローチャートを用いて後述する。
【００７３】
上述した図３のステップＳ３０１では、以上のように生成されたメタメリックペアの印刷色に基づいてパッチを印刷することにより、本実施形態における光源推定用パッチ画像９０２が得られる。
【００７４】
●メタメリックペア算出処理
以下、上述したステップＳ１１０６におけるメタメリックペアの算出処理について、図１５のフローチャート等を用いて詳細に説明する。
【００７５】
まずステップＳ１５０１において、ステップＳ１５０２〜Ｓ１５０８でのメタメリックペア算出対象となるひとつの光源ｊを選択し、更に、指定された光源ｊに対応するＬＵＴ1jとＬＵＴ2jを、メインメモリ１０２から取得する。なお、ここでの光源ｊの選択方法については特に限定しない。
【００７６】
ステップＳ１５０２では、選択された光源ｊにおけるメタメリックペアの一方であるＣＭＹＫ印刷用のパッチ色について、そのＣ'Ｍ'Ｙ'値であるＣ'Ｍ'Ｙ'1jを更新する。なお、Ｃ'Ｍ'Ｙ'値の更新方法については特に限定しない。
【００７７】
次にステップＳ１５０３では、選択された光源ｊに対応するＬＵＴ1jとＬＵＴ2jに基づき、光源ｊの下でＣ'Ｍ'Ｙ'1jで印刷されたＣＭＹＫパッチに等色する、ＲＧＢ印刷用のパッチのＣ'Ｍ'Ｙ'値であるＣ'Ｍ'Ｙ'2jを算出する。この算出は例えば、補間技法と探索技法とを組み合わせることによって行うことができる。
【００７８】
ステップＳ１５０４では、ステップＳ１５０３で算出したメタメリックペアを光源ｊとは異なる他の光源ｋで観察した際の平均色差ΔＥを、ＬＵＴ1kとＬＵＴ2kに基づいて以下の手順により算出する。すなわち、まずＣＭＹＫ印刷パッチを光源ｋの下で観察した際の色を、ＬＵＴ1kとＣ'Ｍ'Ｙ'1jとから補間により算出する。そして、ＲＧＢ印刷パッチを光源ｋの下で観察した際の色を、ＬＵＴ2kとＣ'Ｍ'Ｙ'2jとから補間により算出する。続いて、得られた２つの色の色差ΔＥk（k≠j）を、ＣＩＥΔＥ９４色差式に基づきＮ−１個算出する。そして、算出したＮ−１個の色差の平均をΔＥとする。但し、ΔＥ９４色差式を適用する際の白色基準としては、ペーパーホワイトのＸＹＺ値を用いるとする。
【００７９】
ステップＳ１５０５では、ステップＳ１５０４で算出した平均色差ΔＥと、既に設定されている最大平均色差ΔＥmaxとの比較を行い、ΔＥ＞ΔＥmaxであればステップＳ１５０６へ進み、そうでなければステップＳ１５０７へ進む。但し、初めてステップＳ１５０５を処理する場合には、最大平均色差ΔＥmaxが未定であるため、無条件でステップＳ１５０６へ進む。
【００８０】
ステップＳ１５０６では、ステップＳ１５０３で算出したメタメリックペアの値Ｃ'Ｍ'Ｙ'1j並びにＣ'Ｍ'Ｙ'2jを、光源判別用のパッチ候補として記憶すると共に、ΔＥmaxをΔＥで更新する。これにより、光源ｊについて、他の光源との平均色差ΔＥが最大となるようなパッチ色を得ることができる。
【００８１】
そしてステップＳ１５０７では、Ｃ'Ｍ'Ｙ'1jの更新が完了したか否かの判定を行い、完了していればステップＳ１５０８へ進み、完了していなければステップＳ１５０２に戻って次のパッチ色についての処理を開始する。
【００８２】
ステップＳ１５０８では、メタメリックペアの候補値、すなわちステップＳ１５０６で最終的に更新された光源判別用の２つのパッチ候補Ｃ'Ｍ'Ｙ'1j並びにＣ'Ｍ'Ｙ'2jを、光源ｊ判別用のパッチのデバイスカラーとして記憶する。
【００８３】
ステップＳ１５０９では、Ｎ種の光源の全てに対し、メタメリックペア印刷用のデバイスカラー算出が終了したか否かを判定する。終了していればステップＳ１５１０へ進み、終了していなければステップＳ１５０１へ戻って次の光源に対する処理を開始する。
【００８４】
ステップＳ１５１０では、上述したステップＳ１５０１〜Ｓ１５０９の処理ループによって算出したＮ組のメタメリックペア印刷用のＣ'Ｍ'Ｙ'デバイスカラー値を、メインメモリ１０２に記憶する。これにより、本実施形態の光源推定用パッチ画像におけるメタメリックペアの算出処理が終了する。
【００８５】
以上説明したように本実施形態によれば、メタメリックペアからなるパッチ画像について、光学フィルタを用いることにより等色／非等色の色差を容易に検出可能とし、最適な光源を容易に推定することができる。したがって、測光センサ等の機器を用いることなく、観察環境を適切に反映したプロファイルを作成することが可能となり、最適なカラーマッチングを行うことができる。
【００８６】
＜第２実施形態＞
以下、本発明に係る第２実施形態について説明する。第２実施形態における画像処理装置のシステム構成は、上述した第１実施形態と同様であるため、説明を省略する。第２実施形態においても上述した第１実施形態と同様に、メタメリックペアからなる光源推定用のパッチ画像を観察することによって、現在の光源に最適なデスティネーションプロファイルを作成する。このディスティネーションプロファイルの作成手順は、上述した第１実施形態で示した図３のフローチャートと同様である。
【００８７】
ただし第２実施形態では、メタメリックペアの一方のカラーパッチを、判別対象の光源に拠らず同一の色で印刷することを特徴とする。このようにメタメリックペアの一方のカラーパッチを同一色とすることで、最も等色するメタメリックペアを同一色に対する比較で判定することができるため、上述した第１実施形態と比較して、最も近い色となるメタメリックペアを判定し易くなる。しかしながら他方で、メタメリックペアの一方のカラーパッチが同一色であるという拘束条件が入るため、非等色時の色差が第１実施形態と比較して小さくなる。
【００８８】
第２実施形態においてはさらに、光学フィルタを２種使用することを特徴とする。第１実施形態のように１種類の光学フィルタを用いる場合、光源を推定するのに満足な、あるいは推定し易い十分な透過率を持つフィルタが市販されていない場合がある。第２実施形態では光学フィルタを２種使用することによって、適当な光学フィルタが入手できないという不都合を回避することが可能である。しかしながら他方では、機材を余分に使用するためにコストや煩雑さが増える惧れがある。
【００８９】
●光源推定用パッチ画像の観察方法
以下、第２実施形態において印刷される光源推定用パッチ画像の例およびその観察方法について、図１６および図１７Ａ，図１７Ｂを用いて詳細に説明する。
【００９０】
図１６において、１６０１は光源、１６０２は光源推定用パッチ画像、１６０３は光学フィルタＡ、１６０４は光学フィルタＢ、１６０５は観察者である。光源推定用パッチ画像１６０２には、図１７Ａに示すように２行のパッチが印刷されており、パッチ数は判別対象の光源数＋１である。第２実施形態では、判別対象の光源数が５種であり、全パッチ数が６である例を示す。図１７Ｂに示すように、光学フィルタＡ１６０３は光源推定用パッチ画像の上のパッチ行のみを覆うように配置され、光学フィルタＢ１６０４は光源推定用パッチ画像の下のパッチ行のみを覆うように配置される。この様な配置により、光源推定用パッチ画像１６０２における上の行のパッチは、一度光学フィルタＡ１６０３を透過した光線を反射し、さらにこの反射光がもう一度光学フィルタＡ１６０３を透過した光線が、観察者１６０５に観察される。また下の行のパッチは、一度光学フィルタＢ１６０４を透過した光線を反射し、さらにこの反射光がもう一度光学フィルタＢ１６０４を透過した光線が、観察者１６０５に観察される。
【００９１】
ここで図１７Ａに示すパッチ画像は、図１７Ｂにおける光源Ｂのように、観察環境における光源（この場合光源Ｂ）に相当する上下のペアパッチが等色、もしくは近しい色として観察されるように生成されている。すなわち、図１７Ａに示すパッチ画像の上段に対して光学フィルタＡ１６０３を、下段に対して光学フィルタＢ１６０４を配置することにより、ペアのパッチについて、等色時と非等色時との色差を大きくすることができる。したがってユーザは、最も互いの色が近く見えたペアを容易に選択することができる。
【００９２】
ここで図１９Ａ，図１９Ｂに、第２実施形態において使用される光学フィルタＡ１６０３，Ｂ１６０４の特性の一例を示す。図１９Ａ，図１９Ｂによれば、これら２つの光学フィルタを利用することにより、上述した第１実施形態で図１４に示したのとほぼ同様の特性を得られることが分かる。
【００９３】
●光源推定用パッチ画像作成処理
以下、図３のステップＳ３０１で印刷される、第２実施形態の光源推定用パッチ画像について、その作成処理を図１８のフローチャート等を用いて詳細に説明する。
【００９４】
まずステップＳ１８０１では、パッチ画像作成のみに使用する、２種の色材料ＬＵＴを生成する。一方は、ＲＧＢの３種のインクを使用するＲＧＢ・ＬＵＴ、他方はＣＭＹＫの４種のインクを使用するＣＭＹＫ・ＬＵＴである。両ＬＵＴともにＣ'Ｍ'Ｙ'値に対応するテーブルであり、テーブルを参照するためのＣ'Ｍ'Ｙ'値の範囲は[０，２５５]である。また、テーブルに記述されているＲＧＢ各インクの値の範囲は同様に[０，２５５]であり、ＣＭＹＫ各インクの値の範囲も同様に[０，２５５]である。ここで、テーブルを参照する際のＣ'Ｍ'Ｙ'のステップを以下に示す。
【００９５】
Ｃ'＝{0,32,64,96,128,160,192,224,255}
Ｍ'＝{0,32,64,96,128,160,192,224,255}
Ｙ'＝{0,32,64,96,128,160,192,224,255}
ＲＧＢ・ＬＵＴに記述されるＲＧＢインク値は、Ｃ'Ｍ'Ｙ'値より以下のように算出される。
【００９６】
Ｒ＝１−Ｃ'，Ｇ＝１−Ｍ'，Ｂ＝１−Ｙ'
また、ＣＭＹＫ・ＬＵＴに記述されるＣＭＹＫインク値は、Ｃ'Ｍ'Ｙ'値より以下のように算出される。
【００９７】
Ｋ＝Min(Ｃ'，Ｍ'，Ｙ')、
Ｃ＝Ｃ'−Ｋ，Ｍ＝Ｍ'−Ｋ，Ｙ＝Ｙ'−Ｋ
ステップＳ１８０２では、作成した２種のＬＵＴそれぞれについて、テーブルに記載された色材量に基づく７２９色のパッチを印刷し、分光反射率を測定してメインメモリ１０２に記憶する。ここで、ＣＭＹＫ・ＬＵＴに対する分光反射率をＲ1i(λ)、ＲＧＢ・ＬＵＴに対する分光反射率をＲ2i(λ)と表す。なお、インデックスｉは夫々のパッチ順序に対応するものである。
【００９８】
ステップＳ１８０３では、ＣＭＹＫ・ＬＵＴに対応する分光反射率Ｒ1i(λ)を、光学フィルタＡ１６０３の分光透過率Ｔ1(λ)の２乗を乗じることにより補正したＲ1'i(λ)を算出する。またステップＳ１８０４では、ＲＧＢ・ＬＵＴに対応する分光反射率Ｒ2i(λ)を、光学フィルタＢ１６０４の分光透過率Ｔ2(λ)の２乗を乗じることにより補正したＲ2'i(λ)を算出する。なお、２種の光学フィルタは予め与えられたものである。
【００９９】
ステップＳ１８０５では、判別対象であるＮ種の光源の相対分光分布Ｓi(λ)をＨＤＤ１０５より取得する。なお、相対分光分布Ｓi(λ)の一例は、図１３Ａに示した通りである。
【０１００】
ステップＳ１８０６では、補正後の分光反射率Ｒ1'i(λ)とＲ2'i(λ)、並びに相対分光分布Ｓj(λ)とから、ｊ番目の光源下で観察した際のｉ番パッチのＸＹＺ値を、下式により算出する。まず、ＣＭＹＫ・ＬＵＴに対するＸＹＺ値であるＸ1ji，Ｙ1ji，Ｚ1jiは以下のようになる。
【０１０１】
2009071550.tif 000009
【０１０２】
また、ＲＧＢ・ＬＵＴに対するＸＹＺ値であるＸ2ji，Ｙ2ji，Ｚ2jiは以下のようになる。
【０１０３】
2009071550.tif 000010
【０１０４】
以下、ｊ番目の光源に基づきＣＭＹＫ・ＬＵＴをＸＹＺ値に変換したものをＬＵＴ1j、ＲＧＢ・ＬＵＴをＸＹＺ値に変換したものをＬＵＴ2jと称する。算出したＬＵＴは、メインメモリ１０２に記憶する。
【０１０５】
ステップＳ１８０７では、ステップＳ１８０６で作成した２×Ｎ種のＸＹＺ値ＬＵＴに基づき、メタメリックペアの印刷色を算出する。この算出手順の詳細については、図２０のフローチャートを用いて後述する。
【０１０６】
第２実施形態では、図３のステップＳ３０１において、以上のように生成されたメタメリックペアの印刷色に基づいてパッチを印刷することにより、光源推定用パッチ画像１６０２が得られる。
【０１０７】
●メタメリックペア算出処理
以下、上述したステップＳ１８０７におけるメタメリックペアの算出処理について、図２０のフローチャート等を用いて詳細に説明する。
【０１０８】
まずステップＳ２００１において、メタメリックペアの一方であるＣＭＹＫ印刷のパッチ色について、そのＣ'Ｍ'Ｙ'値であるＣ'Ｍ'Ｙ'1を更新する。
【０１０９】
次にステップＳ２００２では、メタメリックペア算出対象となるひとつの光源ｊを選択し、更に、指定された光源ｊに対応するＬＵＴ1jとＬＵＴ2jを、メインメモリ１０２から取得する。
【０１１０】
次にステップＳ２００３では、選択された光源ｊに対応するＬＵＴ1jとＬＵＴ2jに基づき、光源ｊの下でＣ'Ｍ'Ｙ'1で印刷されたＣＭＹＫパッチに等色する、ＲＧＢ印刷用のパッチのＣ'Ｍ'Ｙ'値であるＣ'Ｍ'Ｙ'2jを算出する。
【０１１１】
ステップＳ２００４では、ステップＳ２００３で算出したメタメリックペアを光源ｊとは異なる他の光源ｋで観察した際の平均色差ΔＥを、ＬＵＴ1kとＬＵＴ2kに基づいて以下の手順により算出する。すなわち、まずＣＭＹＫ印刷パッチを光源ｋの下で観察した際の色を、ＬＵＴ1kとＣ'Ｍ'Ｙ'1とから補間により算出する。そして、ＲＧＢ印刷パッチを光源ｋの下で観察した際の色を、ＬＵＴ2kとＣ'Ｍ'Ｙ'2jとから補間により算出する。続いて、得られた２つの色の色差ΔＥk（k≠j）を、ＣＩＥΔＥ９４色差式に基づきＮ−１個算出する。そして、算出したＮ−１個の色差の平均をΔＥmean_jとし、最小色差をΔＥmin_jとする。但し、ΔＥ９４色差式を適用する際の白色基準としては、ペーパーホワイトのＸＹＺ値を用いるとする。
【０１１２】
ステップＳ２００５では、Ｎ種の光源の全てに対し、メタメリックペア算出が終了したか否かを判定する。終了していればステップＳ２００６へ進み、終了していなければステップＳ２００２へ戻って次の光源に対する処理を開始する。
【０１１３】
ステップＳ２００６では、ステップＳ２００３〜Ｓ２００５の処理で算出したＮ個の最小色差ΔＥmin_jをチェックし、全てのΔＥmin_jが１以上であればステップＳ２００７へ進む。一方、１未満のものが１つでもあれば、当該Ｃ'Ｍ'Ｙ'1に対して算出されたメタメリックペアは、全ての光源に対して十分な色差が確保できないため、このパッチを非採用としてステップＳ２００９へ進む。
【０１１４】
ステップＳ２００７では、ステップＳ２００４で算出した平均色差ΔＥmean_jの平均ΔＥを算出し、既に設定されている最大平均色差ΔＥmaxとの比較を行う。該比較の結果、ΔＥ＞ΔＥmaxであればステップＳ２００８へ進み、そうでなければステップＳ２００９へ進む。但し、初めてステップＳ２００７を処理する場合には、最大平均色差ΔＥmaxが未定であるため、無条件でステップＳ２００８へ進む。
【０１１５】
ステップＳ２００８では、ステップＳ２００２〜Ｓ２００５の処理で算出したメタメリックペアの値Ｃ'Ｍ'Ｙ'1並びにＣ'Ｍ'Ｙ'2jを、光源ｊ判別用のパッチのデバイスカラーとして記憶すると共に、ΔＥmaxをΔＥで更新する。
【０１１６】
そしてステップＳ２００９では、Ｃ'Ｍ'Ｙ'1の更新が完了したか否かの判定を行い、完了していればステップＳ２０１０へ進み、完了していなければステップＳ２００１に戻って次のパッチ色についての処理を開始する。
【０１１７】
ステップＳ２０１０では、上述したステップＳ２００１〜Ｓ２００９の処理ループによって算出したＮ組のメタメリックペア印刷用であるＮ＋１色のＣ'Ｍ'Ｙ'デバイスカラー値を、メインメモリ１０２に記憶する。これにより、第２実施形態の光源推定用パッチ画像におけるメタメリックペアの算出処理が終了する。
【０１１８】
以上説明したように第２実施形態によれば、メタメリックペアの一方のカラーパッチを、判別対象の光源に拠らず同一色とすることにより、最も等色するメタメリックペアを同一色に対して比較・判定すればよいため、より容易な判定が可能となる。さらに、光学フィルタを２種使用することにより、光学フィルタを容易に入手することができる。
【０１１９】
＜第３実施形態＞
以下、本発明に係る第３実施形態について説明する。第３実施形態における画像処理装置のシステム構成は、上述した第１実施形態と同様であるため、説明を省略する。第３実施形態においても上述した第１実施形態と同様に、メタメリックペアからなる光源推定用のパッチ画像を観察することによって、現在の光源に最適なデスティネーションプロファイルを作成する。このディスティネーションプロファイルの作成手順は、上述した第１実施形態で示した図３のフローチャートと同様である。
【０１２０】
ただし第３実施形態では、メタメリックペアの一方のカラーパッチを、判別対象の光源に拠らず同一の色で印刷するとともに、光源に応じた２種の光学フィルタを用いることを特徴とする。このように、推定すべき光源に応じて光学フィルタを変更することにより、非等色時の色差を大きくすることができ、より判別し易くすることが可能である。しかしながら他方では、機材を余分に使用するためにコストや煩雑さが増える惧れがある。
【０１２１】
また第３実施形態では、安価なプリンタへの適用を考慮して、特にＣＭＹプリンタにおける光源推定用パッチ作成例を示す。但し、インクに起因するメタメリズムの効果を期待できないため、非等色時の色差を大きく取れない惧れがある。第３実施形態において２種の光学フィルタを光源に応じて変更することにより、この点が補完される。
【０１２２】
さらに第３実施形態では、光学フィルタを紙に密着させないようにして、カラーパッチの観察を行う。これは、例えば干渉フィルタ等、光学フィルタの種類に拠っては、その特性から紙に密着させることが不適な場合があるためである。
【０１２３】
●光源推定用パッチ画像の観察方法
以下、第３実施形態において印刷される光源推定用パッチ画像の例およびその観察方法について、図２１および図２２Ａ，図２２Ｂを用いて詳細に説明する。
【０１２４】
図２１において、２１０１は光源、２１０２は光源推定用パッチ画像、２１０３は光学フィルタＡ、２１０４は光学フィルタＢ、２１０５は観察者である。なお、光学フィルタＡ２１０３と光学フィルタＢ２１０４は、観察するカラーパッチに応じて使用を切り替える。光源推定用パッチ画像２１０２には、図２２Ａに示すように２行のパッチが印刷されており、パッチ数は判別対象の光源数＋１である。第３実施形態では、判別対象の光源数が５種であり、全パッチ数が６である例を示す。図２２Ｂに示すように、光学フィルタＡ２１０３は光源推定用パッチ画像の下のパッチ行における左２つを観察する際に用いられ、光学フィルタＢ２１０４は同じく下のパッチ行における右３つを観察する際に用いられる。なお、図２２Ｂでは便宜的にパッチの直上に光学フィルタを設置してあるかのように図示しているが、実際には図２１に示す様に、光学フィルタはパッチ画像に触れないように設置する。この様な配置により、光源推定用パッチ画像２１０２における上の行のパッチは、光源からの光線を通常の様に反射して観察者２１０５に観察される。これに対し、下の行のパッチは、光源からの光線を通常のように反射し、その反射光が光学フィルタを透過した光線が、観察者２１０５に観察される。
【０１２５】
ここで図２２Ａに示すパッチ画像は、図２２Ｂにおける光源Ｂのように、観察環境における光源（この場合光源Ｂ）に相当する上下のペアパッチが等色、もしくは近しい色として観察されるように生成されている。
【０１２６】
●光源推定用パッチ画像作成処理
以下、図３のステップＳ３０１で印刷される、第３実施形態の光源推定用パッチ画像について、その作成処理を図２３のフローチャート等を用いて詳細に説明する。
【０１２７】
まずステップＳ２３０１では、パッチ画像作成のみに使用する、ＣＭＹの３種のインクを使用するＬＵＴ（以下、ＣＭＹ・ＬＵＴ）を生成する。このＣＭＹ・ＬＵＴは、Ｃ'Ｍ'Ｙ'値に対応するテーブルであり、テーブルを参照するためのＣ'Ｍ'Ｙ'値の範囲は[０，２５５]である。また、テーブルに記述されているＣＭＹ各インクの値の範囲も同様に[０，２５５]である。ここで、テーブルを参照する際のＣ'Ｍ'Ｙ'のステップを以下に示す。
【０１２８】
Ｃ'＝{0,32,64,96,128,160,192,224,255}
Ｍ'＝{0,32,64,96,128,160,192,224,255}
Ｙ'＝{0,32,64,96,128,160,192,224,255}
ここで、ＣＭＹ・ＬＵＴに記述されるＣＭＹインク値は、以下に示すようにＣ'Ｍ'Ｙ'値と等しい。
【０１２９】
Ｃ＝Ｃ'，Ｍ＝Ｍ'，Ｙ＝Ｙ'
ステップＳ２３０２では、作成したＣＭＹ・ＬＵＴについて、テーブルに記載された色材量に基づく７２９色のパッチを印刷し、分光反射率Ｒi(λ)を測定してメインメモリ１０２に記憶する。なお、インデックスｉはパッチ順序に対応するものである。
【０１３０】
ステップＳ２３０３では、分光反射率Ｒi(λ)に光学フィルタＡ２１０３の分光透過率ＴA(λ)を乗じたＲAi(λ)と、光学フィルタＢ２１０４の分光透過率ＴB(λ)を乗じたＲBi(λ)を算出する。ここで分光透過率を乗ずるのは、光源からの光線は印刷物で反射された後で更に光学フィルタを透過することを考慮するためである。
【０１３１】
なお、本実施形態で使用する２種の光学フィルタは予め与えられたものであるため、それぞれの分光透過率もこの算出の時点で既に定まった値を有する。使用する光学フィルタを選択する際には、その分光透過率特性について次の注意が必要である。
【０１３２】
すなわち光学フィルタとしては、上述した第１実施形態と同様に、分光分布の関数形状が異なる波長域において、透過光を遮断するものが効果的である。更に第３実施形態では、推定対象の光源毎に光学フィルタを選択することができる。ここで、例えば図１３Ａに示したＦ１０光源の様な三波長型の光源では、他の光源に対して波長のピークが際立って高いという特徴を有すると共に、色温度が異なる光源でもピークの高さが大きく異なる。そこで、特に三波長型の光源を推定するためのメタメリックペアに対して非条件等色時の色差を大きくするには、ピークを有する波長域で透過光を遮断するものが効果的である。そこで第３実施形態では、光学フィルタＡ２１０３，光学フィルタＢ２１０４として、それぞれ図２４Ａ，図２４Ｂに示す特性を有するフィルタを使用する。特に、図２４Ａの特性を有する光学フィルタＡ２１０３は三波長型以外の光源の推定用として利用し、図２４Ｂの特性を有する光学フィルタＢ２１０４は三波長型の光源の推定用として利用する。
【０１３３】
図２３に戻り、ステップＳ２３０４では、判別対象であるＮ種の光源の相対分光分布Ｓi(λ)をＨＤＤ１０５より取得する。
【０１３４】
ステップＳ２３０５では、分光反射率Ｒi(λ)と相対分光分布Ｓj(λ)とから、ｊ番目の光源下で観察した際のｉ番パッチのＸＹＺ値を、下式により算出する。
【０１３５】
2009071550.tif 000011
【０１３６】
以下、このステップで作成したｊ番目の光源に対するＸＹＺ値によるＬＵＴをＬＵＴ1jと称する。
【０１３７】
ステップＳ２３０６では、分光反射率ＲAi(λ)と相対分光分布Ｓj(λ)とから、j番目の光源下で観察した際のｉ番パッチのＸＹＺ値を、下式により算出する。
【０１３８】
2009071550.tif 000012
【０１３９】
以下、このステップで生成したj番目の光源に対するＸＹＺ値によるＬＵＴをＬＵＴAjと称する。
【０１４０】
ステップＳ２３０７では、分光反射率ＲBi(λ)と相対分光分布Ｓj(λ)とから、j番目の光源下で観察した際のｉ番パッチのＸＹＺ値を、下式により算出する。
【０１４１】
2009071550.tif 000013
【０１４２】
以下、このステップで生成したj番目の光源に対するＸＹＺ値によるＬＵＴをＬＵＴBjと称する。
【０１４３】
ステップＳ２３０８では、以上のように作成した３×Ｎ種のＸＹＺ値ＬＵＴに基づき、メタメリックペアの印刷色を算出する。この算出手順の詳細については、図２５のフローチャートを用いて後述する。
【０１４４】
第３実施形態では、図３のステップＳ３０１において、以上のように生成されたメタメリックペアの印刷色に基づいてパッチを印刷することにより、光源推定用パッチ画像２１０２が得られる。
【０１４５】
●メタメリックペア算出処理
以下、上述したステップＳ２３０８におけるメタメリックペアの算出処理について、図２５のフローチャート等を用いて詳細に説明する。
【０１４６】
まずステップＳ２５０１において、メタメリックペアの一方である、光学フィルタを介さずに観察されるパッチ色について、そのＣ'Ｍ'Ｙ'値であるＣ'Ｍ'Ｙ'1を更新する。
【０１４７】
次にステップＳ２５０２では、メタメリックペア算出対象となる光源ｊを選択する。
【０１４８】
次にステップＳ２５０３では、選択された光源ｊが、光学フィルタＡ２１０３と光学フィルタＢ２１０４のいずれに対応する光源であるかに応じて、パッチのＣ'Ｍ'Ｙ'値を算出する。ここで、光学フィルタＡ２１０３に対応する光源とはすなわち、図２２Ｂの例では光源Ａあるいは光源Ｂであり、光学フィルタＢ２１０４に対応する光源とは、図２２Ｂの例では光源Ｃ，光源Ｄ，光源Ｅのいずれかである。光源ｊが光学フィルタＡ２１０３に対応する場合、光源ｊに対応するＬＵＴjとＬＵＴAJに基づき、光源ｊの下でＣ'Ｍ'Ｙ'1で印刷されたＣＭＹパッチに等色するパッチのＣ'Ｍ'Ｙ'値であるＣ'Ｍ'Ｙ'2jを算出する。一方、光源ｊが光学フィルタＢ２１０４に対応する場合、光源ｊに対応するＬＵＴjとＬＵＴBJに基づき、光源ｊの下でＣ'Ｍ'Ｙ'1で印刷されたＣＭＹパッチに等色するパッチのＣ'Ｍ'Ｙ'値であるＣ'Ｍ'Ｙ'2jを算出する。
【０１４９】
ステップＳ２５０４では、ステップＳ２５０３で算出したメタメリックペアを光源ｊとは異なる他の光源ｋで観察した際の平均色差ΔＥを、以下の手順により算出する。
【０１５０】
すなわち、光源ｊが光学フィルタＡ２１０３に対応する場合には、まず該光学フィルタを介さずに観察した際の色を、ＬＵＴkとＣ'Ｍ'Ｙ'1とから補間ににより算出する。そして、光学フィルタＡ２１０３を介して観察した際の色を、ＬＵＴAkとＣ'Ｍ'Ｙ'2jとから補間により算出する。また、光源ｊが光学フィルタＢ２１０４に対応する場合には、まず該光学フィルタを介さずに観察した際の色を、ＬＵＴkとＣ'Ｍ'Ｙ'1とから補間により算出する。そして、光学フィルタＢ２１０４を介して観察した際の色を、ＬＵＴBkとＣ'Ｍ'Ｙ'2jとから補間により算出する。
【０１５１】
そして、得られた２つの色の色差ΔＥk（k≠j）を、ＣＩＥΔＥ９４色差式に基づきＮ−１個算出する。そして、算出したＮ−１個の色差の平均をΔＥmean_jとし、最小色差をΔＥmin_jとする。但し、ΔＥ９４色差式を適用する際の白色基準としては、ペーパーホワイトのＸＹＺ値を用いるとする。
【０１５２】
ステップＳ２５０５では、Ｎ種の光源の全てに対し、メタメリックペア算出が終了したか否かを判定する。終了していればステップＳ２５０６へ進み、終了していなければステップＳ２５０２へ戻って次の光源に対する処理を開始する。
【０１５３】
ステップＳ２５０６では、ステップＳ２５０３〜Ｓ２５０５の処理で算出したＮ個の最小色差ΔＥmin_jをチェックし、全てのΔＥmin_jが１以上であればステップＳ２５０７へ進み、１未満のものが１つでもあれば、ステップＳ２５０９へ進む。
【０１５４】
ステップＳ２５０７では、ステップＳ２５０４で算出した平均色差ΔＥmean_jの平均ΔＥを算出し、既に設定されている最大平均色差ΔＥmaxとの比較を行う。該比較の結果、ΔＥ＞ΔＥmaxであればステップＳ２５０８へ進み、そうでなければステップＳ２５０９へ進む。但し、初めてステップＳ２５０７を処理する場合には、最大平均色差ΔＥmaxが未定であるため、無条件でステップＳ２５０８へ進む。
【０１５５】
ステップＳ２５０８では、ステップＳ２５０２〜Ｓ２５０５の処理で算出したメタメリックペアの値Ｃ'Ｍ'Ｙ'1並びにＣ'Ｍ'Ｙ'2jを、光源ｊ判別用のパッチのデバイスカラーとして記憶すると共に、ΔＥmaxをΔＥで更新する。
【０１５６】
そしてステップＳ２５０９では、Ｃ'Ｍ'Ｙ'1の更新が完了したか否かの判定を行い、完了していればステップＳ２５１０へ進み、完了していなければステップＳ２５０１に戻って次のパッチ色についての処理を開始する。
【０１５７】
ステップＳ２５１０では、上述したステップＳ２５０１〜Ｓ２５０９の処理ループによって算出したＮ組のメタメリックペア印刷用であるＮ＋１色のＣ'Ｍ'Ｙ'デバイスカラー値を、メインメモリ１０２に記憶する。これにより、第３実施形態の光源推定用パッチ画像におけるメタメリックペアの算出処理が終了する。
【０１５８】
以上説明したように第３実施形態によれば、上述した第１および第２実施形態に対してさらに、光源に応じて光学フィルタを変更することにより、非等色時の色差を確実に大きくすることができ、より容易な判定が可能となる。
【０１５９】
以上説明した第１乃至第３実施形態においては、光学フィルタとして急峻な特性を有する例を示したが、第１実施形態でも説明したように、ブロードな特性を有するカラーフィルタを用いることも可能である。
【０１６０】
また、以上説明した第１乃至第３実施形態においては、使用する光学フィルタは多くとも２種類であった。しかしながら、光学フィルタの数は制限されるものでは無く、例えばカラーパッチ数分の光学フィルタを用いることも可能である。
【０１６１】
また、光源推定用パッチ画像におけるメタメリックペアを直接、視覚的に観察する例を示したが、これを撮像装置によって撮像し、該撮像画像を観察しても、本発明と同様の効果が得られる。
【０１６２】
また、光源推定用パッチ画像におけるメタメリックペアの配置については、特に限定しない。上述した各実施形態では直線状に配置する例を示したが、例えばこれを円環状に配置しても構わない。
【０１６３】
＜他の実施形態＞
以上、実施形態例を詳述したが、本発明は例えば、システム、装置、方法、プログラム若しくは記録媒体(記憶媒体)等としての実施態様をとることが可能である。具体的には、複数の機器（例えば、ホストコンピュータ、インタフェース機器、撮像装置、webアプリケーション等）から構成されるシステムに適用しても良いし、また、一つの機器からなる装置に適用しても良い。
【０１６４】
尚本発明は、前述した実施形態の機能を実現するソフトウェアプログラムを、システムあるいは装置に直接あるいは遠隔から供給し、そのシステムあるいは装置のコンピュータが該供給されたプログラムコードを読み出して実行することによっても達成される。なお、この場合のプログラムとは、コンピュータ読取可能であり、実施形態において図に示したフローチャートに対応したプログラムである。
【０１６５】
従って、本発明の機能処理をコンピュータで実現するために、該コンピュータにインストールされるプログラムコード自体も本発明を実現するものである。つまり、本発明は、本発明の機能処理を実現するためのコンピュータプログラム自体も含まれる。
【０１６６】
その場合、プログラムの機能を有していれば、オブジェクトコード、インタプリタにより実行されるプログラム、OSに供給するスクリプトデータ等の形態であっても良い。
【０１６７】
プログラムを供給するための記録媒体としては、以下に示す媒体がある。例えば、フロッピー（登録商標）ディスク、ハードディスク、光ディスク、光磁気ディスク、MO、CD-ROM、CD-R、CD-RW、磁気テープ、不揮発性のメモリカード、ROM、DVD(DVD-ROM，DVD-R)などである。
【０１６８】
プログラムの供給方法としては、以下に示す方法も可能である。すなわち、クライアントコンピュータのブラウザからインターネットのホームページに接続し、そこから本発明のコンピュータプログラムそのもの(又は圧縮され自動インストール機能を含むファイル)をハードディスク等の記録媒体にダウンロードする。また、本発明のプログラムを構成するプログラムコードを複数のファイルに分割し、それぞれのファイルを異なるホームページからダウンロードすることによっても実現可能である。つまり、本発明の機能処理をコンピュータで実現するためのプログラムファイルを複数のユーザに対してダウンロードさせるWWWサーバも、本発明に含まれるものである。
【０１６９】
また、本発明のプログラムを暗号化してCD-ROM等の記憶媒体に格納してユーザに配布し、所定の条件をクリアしたユーザに対し、インターネットを介してホームページから暗号化を解く鍵情報をダウンロードさせることも可能である。すなわち該ユーザは、その鍵情報を使用することによって暗号化されたプログラムを実行し、コンピュータにインストールさせることができる。
【０１７０】
また、コンピュータが、読み出したプログラムを実行することによって、前述した実施形態の機能が実現される。さらに、そのプログラムの指示に基づき、コンピュータ上で稼動しているOSなどが、実際の処理の一部または全部を行い、その処理によっても前述した実施形態の機能が実現され得る。
【０１７１】
さらに、記録媒体から読み出されたプログラムが、コンピュータに挿入された機能拡張ボードやコンピュータに接続された機能拡張ユニットに備わるメモリに書き込まれた後、実行されることによっても、前述した実施形態の機能が実現される。すなわち、該プログラムの指示に基づき、その機能拡張ボードや機能拡張ユニットに備わるCPUなどが実際の処理の一部または全部を行うことが可能である。
【図面の簡単な説明】
【０１７２】
【図１】本発明に係る一実施形態における画像処理装置のシステム構成を示すブロック図である。である光源推定方法並びに該推定用カラーパッチ作成装置の実施形態である、
【図２】本実施形態におけるＳｏｆｔ・ＲＩＰによる画像出力処理を示すフローチャートである。
【図３】本実施形態におけるデスティネーションプロファイルの作成処理を示すフローチャートである。
【図４】本実施形態における光源選択用のユーザインタフェース例を示す図である。
【図５】本実施形態における印刷色再現特性データの算出処理を示すフローチャートである。
【図６】ＸＹＺ等色関数の関数形状を示す図である。
【図７】本実施形態におけるＩＣＣプロファイルの生成処理を示すフローチャートである。
【図８】本実施形態における色域写像方式を説明する図である。
【図９】本実施形態における光源推定用パッチ画像を観察するための構成を示す図である。
【図１０Ａ】本実施形態における光源推定用パッチ画像例を示す図である。
【図１０Ｂ】本実施形態における光源推定用パッチ画像を光学フィルタを介して観察した例を示す図である。
【図１１】本実施形態における光源推定用パッチ画像の作成処理を示すフローチャートである。
【図１２Ａ】ＣＭＹＫ印刷色の分光反射率の一例を示す図である。
【図１２Ｂ】ＲＧＢ印刷色の分光反射率の一例を示す図である。
【図１３Ａ】光源の相対分光分布例を示す図である。
【図１３Ｂ】光源の相対分光分布と等色関数を比較する図である。
【図１４】本実施形態における光学フィルタの分光透過率を示す図である。
【図１５】本実施形態におけるメタメリックペア算出処理を示すフローチャートである。
【図１６】第２実施形態における光源推定用パッチ画像を観察するための構成を示す図である。
【図１７Ａ】第２実施形態における光源推定用パッチ画像例を示す図である。
【図１７Ｂ】第２実施形態における光源推定用パッチ画像を光学フィルタを介して観察した例を示す図である。
【図１８】第２実施形態における光源推定用パッチ画像の作成処理を示すフローチャートである。
【図１９Ａ】第２実施形態における光学フィルタＡの分光透過率を示す図である。
【図１９Ｂ】第２実施形態における光学フィルタＢの分光透過率を示す図である。
【図２０】第２実施形態におけるメタメリックペア算出処理を示すフローチャートである。
【図２１】第３実施形態における光源推定用パッチ画像を観察するための構成を示す図である。
【図２２Ａ】第３実施形態における光源推定用パッチ画像例を示す図である。
【図２２Ｂ】第３実施形態における光源推定用パッチ画像を光学フィルタを介して観察した例を示す図である。
【図２３】第３実施形態における光源推定用パッチ画像の作成処理を示すフローチャートである。
【図２４Ａ】第３実施形態における光学フィルタＡの分光透過率を示す図である。
【図２４Ｂ】第３実施形態における光学フィルタＢの分光透過率を示す図である。
【図２５】第３実施形態におけるメタメリックペア算出処理を示すフローチャートである。
【図２６】従来のインク特性に依存して生成したメタメリックペアの分光反射率の例を示す図である。
【図１】
2009071550.tif 000014
【図２】
2009071550.tif 000015
【図３】
2009071550.tif 000016
【図４】
2009071550.tif 000017
【図５】
2009071550.tif 000018
【図６】
2009071550.tif 000019
【図７】
2009071550.tif 000020
【図８】
2009071550.tif 000021
【図９】
2009071550.tif 000022
【図１０Ａ】
2009071550.tif 000023
【図１０Ｂ】
2009071550.tif 000024
【図１１】
2009071550.tif 000025
【図１２Ａ】
2009071550.tif 000026
【図１２Ｂ】
2009071550.tif 000027
【図１３Ａ】
2009071550.tif 000028
【図１３Ｂ】
2009071550.tif 000029
【図１４】
2009071550.tif 000030
【図１５】
2009071550.tif 000031
【図１６】
2009071550.tif 000032
【図１７Ａ】
2009071550.tif 000033
【図１７Ｂ】
2009071550.tif 000034
【図１８】
2009071550.tif 000035
【図１９Ａ】
2009071550.tif 000036
【図１９Ｂ】
2009071550.tif 000037
【図２０】
2009071550.tif 000038
【図２１】
2009071550.tif 000039
【図２２Ａ】
2009071550.tif 000040
【図２２Ｂ】
2009071550.tif 000041
【図２３】
2009071550.tif 000042
【図２４Ａ】
2009071550.tif 000043
【図２４Ｂ】
2009071550.tif 000044
【図２５】
2009071550.tif 000045
【図２６】
2009071550.tif 000046