(57)【要約】【課題】実際に印刷を行う印刷装置で使用される色材群とは異なる色材群の混色を予測し、該予測された混色を該印刷装置にて印刷して視認可能にする。【解決手段】前記第１の色材群と異なる第２の色材群を利用して作成される混色の分光反射率を、混色元となる前記第２の色材群における各色材の分光反射率と前記混色における前記第２の色材群の各色材の使用割合とに基づいて予測する分光反射率予測手段と、前記混色の分光反射率に近似する分光反射率を前記記録媒体上に再現する前記色材量セットを予測する色材量セット予測手段と、当該予測した色材量セットに基づいて前記印刷装置に印刷を実行させる印刷制御手段と、を備えさせる。【選択図】図８
2009177789.tif 000002
【特許請求の範囲】
【請求項１】
第１の色材群の色材を記録媒体に付着させる印刷を印刷装置に実行させるにあたり、前記第１の色材群における色材の量の組み合わせである色材量セットに基づいた印刷を前記印刷装置に実行させる印刷制御装置であって、
前記色材量セットとインデックスとの対応関係を規定することにより作成されたルックアップテーブルを参照することにより、指定されたインデックスに対応する前記色材量セットを前記印刷装置に指定して印刷させる印刷手段を備え、
前記ルックアップテーブルにおいては、
前記第１の色材群と異なる第２の色材群を利用して作成される混色について、混色元となる前記第２の色材群における各色材の分光反射率と前記混色における前記第２の色材群の各色材の使用割合とに基づいて所定の予測モデルで予測された分光反射率に近似する分光反射率を前記記録媒体上に再現するように予測された前記色材量セットと、前記混色を特定するインデックスと、の対応関係が規定されていることを特徴とする印刷制御装置。
【請求項２】
前記色材量セット予測手段は、前記第１の色材群における各色材の分光反射率と前記第１の色材群における各色材の使用割合とに基づいて分光反射率予測手段によって予測された結果が、前記混色の分光反射率に近似する分光反射率となるように、前記第１の色材群における各色材の使用割合を変化させて繰り返し分光反射率予測手段に予測させることで、前記色材量セットの予測を行うことを特徴とする請求項１に記載の印刷制御装置。
【請求項３】
前記色材量セットの予測は、波長によって異なる重みを加味しつつ前記混色の分光反射率に対する近似性を評価する評価値に基づいて行うことを特徴とする請求項１又は２に記載の印刷制御装置。
【請求項４】
前記重みは、人間の目の分光感度特性に基づいて設定されることを特徴とする請求項３に記載の印刷制御装置。
【請求項５】
前記重みは、前記ターゲット分光反射率に基づいて設定されることを特徴とする請求項３に記載の印刷制御装置。
【請求項６】
前記重みは、所定の光源の分光エネルギー分布に基づいて設定されることを特徴とする請求項３に記載の印刷制御装置。
【請求項７】
第１の色材群を記録媒体に付着させる印刷を実行する印刷装置と、前記第１の色材群における色材の使用量の組み合わせである色材量セットを前記印刷装置に指定し、当該色材量セットに基づいた印刷を実行させる印刷制御装置とからなる印刷システムであって、
前記印刷制御装置は、
前記色材量セットとインデックスとの対応関係を規定することにより作成されたルックアップテーブルを参照することにより、指定されたインデックスに対応する前記色材量セットを前記印刷装置に指定して印刷させる印刷手段を備え、
前記ルックアップテーブルにおいては、
前記第１の色材群と異なる第２の色材群を利用して作成される混色について、混色元となる前記第２の色材群における各色材の分光反射率と前記混色における前記第２の色材群の各色材の使用割合とに基づいて所定の予測モデルで予測された分光反射率に近似する分光反射率を前記記録媒体上に再現するように予測された前記色材量セットと、前記混色を特定するインデックスと、の対応関係が規定されており、
前記印刷装置は、
前記色材量セットに基づいて印刷を実行する印刷実行手段を具備することを特徴とする印刷システム。
【請求項８】
第１の色材群の色材を記録媒体に付着させる印刷を印刷装置に実行させるにあたり、前記第１の色材群における色材の量の組み合わせである色材量セットに基づいた印刷を前記印刷装置に実行させる機能をコンピュータに実現させるための印刷制御プログラムであって、
前記色材量セットとインデックスとの対応関係を規定することにより作成されたルックアップテーブルを参照することにより、指定されたインデックスに対応する前記色材量セットを前記印刷装置に指定して印刷させる印刷機能をコンピュータに実現させており、
前記ルックアップテーブルにおいては、
前記第１の色材群と異なる第２の色材群を利用して作成される混色について、混色元となる前記第２の色材群における各色材の分光反射率と前記混色における前記第２の色材群の各色材の使用割合とに基づいて所定の予測モデルで予測された分光反射率に近似する分光反射率を前記記録媒体上に再現するように予測された前記色材量セットと、前記混色を特定するインデックスと、の対応関係が規定されていることを特徴とする印刷制御プログラム。 
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、印刷制御装置、印刷システムおよび印刷制御プログラムに関し、特に混色の分光反射率を再現させるための印刷制御装置、印刷システムおよび印刷制御プログラムに関する。
【背景技術】
【０００２】
分光的な再現性に注目した印刷方法が提案されている（特許文献１参照。）。この文献においては、ターゲット画像に分光的かつ測色的に一致するような印刷を行うために、プリンティングモデルを使用し、ターゲットの分光反射率（ターゲットスペクトル）にフィッティングするようにプリンタ色（ＣＭＹＫＯＧ）の組み合わせを最適化している。このようにすれば、当該プリンタ色（ＣＭＹＫＯＧ）に基づく印刷を行うことにより、分光的にターゲット画像が再現でき、結果として測色的にも再現性の高い印刷結果を得ることができる。
【特許文献１】特表２００５−５０８１２５号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
ところで、印刷業界には、正式な印刷を行う前に印刷結果を確認するためのプルーフ等、他の印刷装置が無い状況においても該他の印刷装置による印刷結果を作成する必要性や要望がある。また、既に作成済みの色（印刷済みの色、自然界に存在する色、絵画や文化財や文書等に既に着色されている色、絵の具の色等）を混ぜ合わせた色を、実際に印刷して目で見て確認したいという要望もある。印刷装置同士であれば、予め印刷装置同士の色空間を対応付けるＬＵＴ等を作成しておき、作成したＬＵＴに基づく印刷を行ったりすることで実現可能ではある。ただし、大きなサイズのＬＵＴを記憶させるとメモリサイズを圧迫することにもなるし、印刷結果を再現したい他の印刷装置が特定されていない場合はこのようなＬＵＴを持つこと自体が現実的でなくなる。また、印刷結果以外で作成済みの色を混ぜ合わせた結果を予測するのは、そもそもＬＵＴ等では実現不可能である。
【０００４】
本発明は、前記課題にかんがみてなされたもので、実際に印刷を行う印刷装置で使用される色材群とは異なる色材群の混色を予測し、該予測された混色を該印刷装置にて印刷して視認可能にする印刷制御装置、印刷システムおよび印刷制御プログラムを提供する。
【課題を解決するための手段】
【０００５】
前記課題を解決するために、本発明の印刷制御装置は、印刷手段と分光反射率予測手段と色材量セット予測手段と混色印刷手段とを備える構成としてある。印刷手段は、前記色材量セットとインデックスとの対応関係を規定することにより作成されたルックアップテーブルを参照することにより、指定されたインデックスに対応する前記色材量セットを前記印刷装置に指定して印刷させる。この印刷において参照されるルックアップテーブルにおいては、第２の色材群を利用して作成される混色を前記記録媒体上に再現するように予測した色材量セットが、前記混色を特定するインデックスに対応付けて規定されている。前記第２の色材群は、前記第１の色材群と異なる色材群である。前記混色を前記第１の色材群を用いて前記記録媒体上に再現するための色材量セットは、混色元となる前記第２の色材群における各色材の分光反射率と前記混色における前記第２の色材群の各色材の使用割合とに基づいて前記混色の分光反射率を所定の予測モデルで予測し、該混色の分光反射率に近似する分光反射率を前記記録媒体上に再現するように予測されたものである。前記混色の分光反射率に近似する分光反射率を前記記録媒体上に再現する前記色材量セットによれば、光源が変動しても混色の元となった色材を実際に混合して作成される色と同様の色を示す印刷結果を得ることができる。
【０００６】
なお、前記分光反射率取得手段は、実際に前記色材についての分光反射率測定を行うことにより前記分光反射率を取得してもよいし、ユーザー等から前記色材の分光反射率の入力を受付けてもよい。また、前記印刷装置は少なくとも複数の前記色材を前記記録媒体に付着させることができればよく、インクジェットプリンタやレーザープリンタや昇華型プリンタ等の種々の印刷装置が採用可能である。
【０００７】
また、前述した色材量セットを予測するにあたり、前記色材量セット予測手段は、前記第１の色材群における各色材の分光反射率と前記第１の色材群における各色材の使用割合とに基づいて分光反射率予測手段によって予測された結果が、前記混色の分光反射率に近似する分光反射率となるように、前記第１の色材群における各色材の使用割合を変化させて繰り返し分光反射率予測手段に予測させることで、前記色材量セットの予測を行う構成としてもよい。このように構成すると、第１の色材群を利用して第２の色材群の混色結果を再現する際に、最も適した色材量セットで印刷を行うことが出来るようになる。
【０００８】
また、予測された色材量セットの近似度合を評価する、波長によって異なる重みを加味しつつ前記混色の分光反射率に対する近似性を評価する評価値に基づいて行うことにより実現可能である。前記重みの好適な一例として、人間の目の分光感度特性に基づいて設定されるものを適用してもよい。このようにすることにより、人間の分光感度が敏感な波長について優先的に分光反射率を近似させることができ、視覚的に再現精度のよい印刷結果を得ることができる。より具体的な例として、前記重みが３刺激値に対応する各等色関数の線形結合に基づいて設定されるようにしてもよい。このようにすれば、３刺激値に対応する各等色関数に対応する波長域を総合的に重視することが可能な前記重みを設定することができる。
【０００９】
さらに、前記重みが前記混色の分光反射率に基づいて設定されてもよい。例えば、前記混色の分光反射率が強いスペクトルを有している波長域については、最終的に視覚に大きく影響すると考えられるため、当該波長域を優先させて近似させることが望ましい。また、前記重みが所定の光源の分光エネルギー分布に基づいて設定されてもよい。所定の光源の分光エネルギー分布に応じて前記重みを設定することにより、光源が強いスペクトルを有している波長域について優先させて近似させることでき、当該光源における視覚的再現性を向上させることができる。また、単一光源のみならず複数光源での再現性を総合的に考慮して、前記重みを複数光源の分光エネルギーの線形結合に基づいて設定されるようにしてもよい。
【００１０】
さらに、本発明の技術的思想は、具体的な印刷制御装置にて具現化されるのみならず、その方法としても具現化することができる。すなわち、上述した印刷制御装置が行う各手段に対応する工程を有する方法としても本発明を特定することができる。むろん、上述した印刷制御装置がプログラムを読み込んで上述した各手段を実現する場合には、当該各手段に対応する機能を実行させるプログラムや当該プログラムを記録した各種記録媒体においても本発明の技術的思想が具現化できることは言うまでもない。なお、本発明の印刷制御装置は、単一の装置のみならず、複数の装置によって分散して存在可能であることはいうまでもない。例えば、印刷制御装置が有する各手段が、パーソナルコンピュータ上で実行されるプリンタドライバと、プリンタの双方において分散することも可能である。また、プリンタ等の印刷装置に本発明の印刷制御装置の各手段を包含させることも可能である。
【発明を実施するための最良の形態】
【００１１】
以下、下記の順序に従って本発明の実施形態を説明する。
１．印刷制御装置の構成：
２．印刷データ生成処理：
３．混色の印刷データ生成処理：
４．印刷制御処理：
４−１．１Ｄ−ＬＵＴ作成処理：
４−２．印刷制御データ生成処理：
５．分光プリンティングモデル：
６．変形例：
６−１：変形例１：
６−２：変形例２：
６−３：変形例３：
６−４：変形例４：
６−５：変形例５：
６−６：変形例６：
６−７：変形例７：
【００１２】
１．印刷制御装置の構成：
図１は、本発明の一実施形態にかかる印刷制御装置のハードウェア構成を示している。同図において、印刷制御装置は主にコンピュータ１０によって構成されており、コンピュータ１０はＣＰＵ１１とＲＡＭ１２とＲＯＭ１３とハードディスクドライブ（ＨＤＤ）１４と汎用インターフェイス（ＧＩＦ）１５とビデオインターフェイス（ＶＩＦ）１６と入力インターフェイス（ＩＩＦ）１７とバス１８とから構成されている。バス１８は、コンピュータ１０を構成する各要素１１〜１７の間でのデータ通信を実現するものであり、図示しないチップセット等によって通信が制御されている。ＨＤＤ１４には、オペレーティングシステム（ＯＳ）を含む各種プログラムを実行するためのプログラムデータ１４ａが記憶されており、当該プログラムデータ１４ａをＲＡＭ１２に展開しながらＣＰＵ１１が当該プログラムデータ１４ａに従って演算を実行する。ＧＩＦ １５は、例えばＵＳＢ規格に準じたインターフェイスを提供するものであり、外部のプリンタ２０と分光反射率計３０をコンピュータ１０に接続させている。ＶＩＦ １６はコンピュータ１０を外部のディスプレイ４０に接続し、ディスプレイ４０に画像を表示するためのインターフェイスを提供する。ＩＩＦ １７はコンピュータ１０を外部のキーボード５０ａとマウス５０ｂに接続し、キーボード５０ａとマウス５０ｂからの入力信号をコンピュータ１０が取得するためのインターフェイスを提供する。
【００１３】
図２は、コンピュータ１０にて実行されるプログラムのソフトウェア構成を概略的なデータの流れとともに示している。同図において、コンピュータ１０では、おもにＯＳ Ｐ１と見本印刷アプリケーション（ＡＰＬ）Ｐ２と１Ｄ−ＬＵＴ生成アプリケーション（ＬＵＧ）Ｐ３ａとプリンタドライバ（ＰＤＶ）Ｐ３ｂと分光反射率計ドライバ（ＭＤＶ）Ｐ４とディスプレイドライバ（ＤＤＶ）Ｐ５が実行されている。ＯＳ Ｐ１は、各プログラムが使用可能なＡＰＩのひとつとして画像機器インターフェイス（ＧＤＩ） Ｐ１ａとスプーラＰ１ｂを提供しており、ＡＰＬ Ｐ２の要求に応じてＧＤＩ Ｐ１ａが呼び出され、さらにＧＤＩ Ｐ１ａの要求に応じてＰＤＶ Ｐ３ｂやＤＤＶ Ｐ５が呼び出される。ＧＤＩ Ｐ１ａはコンピュータ１０がプリンタ２０やディスプレイ４０等の画像出力装置における画像出力を制御するための汎用的な仕組みを提供し、一方のＰＤＶ Ｐ３ｂやＤＤＶ Ｐ５はプリンタ２０やディスプレイ４０の機種固有の処理等を提供する。また、スプーラＰ１ｂは、ＡＰＬ Ｐ２やＰＤＶ Ｐ３ｂやプリンタ２０の間に介在し、ジョブのコントロール等を実行する。ＡＰＬ Ｐ２は、見本チャートＳＣを印刷するためのアプリケーションプログラムであり、ＲＧＢビットマップ形式の印刷データＰＤを生成し、ＧＤＩ Ｐ１ａに対して当該印刷データＰＤを出力する。また、印刷データＰＤを生成するにあたっては、ＭＤＶ Ｐ４からターゲットの測色データＭＤを取得する。ＭＤＶ Ｐ４は、ＡＰＬ Ｐ２の要求に応じて分光反射率計３０を制御し、当該制御によって得られた分光反射率データＲＤをＡＰＬ Ｐ２に出力する。
【００１４】
ＡＰＬ Ｐ２が生成した印刷データＰＤはＧＤＩ Ｐ１ａやスプーラＰ１ｂを経由してＰＤＶ Ｐ３ｂに出力され、ＰＤＶ Ｐ３ｂが印刷データＰＤに基づいてプリンタ２０に出力可能な印刷制御データＣＤを生成する処理を実行する。ＰＤＶ Ｐ３ｂが生成した印刷制御データＣＤはＯＳ Ｐ１が提供するスプーラＰ１ｂを介してプリンタ２０に出力され、プリンタ２０が当該印刷制御データＣＤに基づく動作を行うことにより見本チャートＳＣを印刷用紙上に印刷させる。以上においては、全体の処理の流れ概略的に説明したが、以下、フローチャートを用いて各プログラムＰ１〜Ｐ４が実行する処理を詳細に説明する。
【００１５】
２．印刷データ生成処理
図３は、ＡＰＬ Ｐ２が実行する印刷データ生成処理の流れを示している。図２に示すようにＡＰＬ Ｐ２はＵＩモジュール（ＵＩＭ）Ｐ２ａと測定制御モジュール（ＭＣＭ）Ｐ２ｂと印刷データ生成モジュール（ＰＤＧ）Ｐ２ｃと混色の印刷データ生成モジュールＰ２ｄとから構成されており、これらの各モジュールＰ２ａ，Ｐ２ｂ，Ｐ２ｃが図３に示す各ステップを実行する。なお、混色の印刷データ生成モジュールＰ２ｄについては後述の３節や６−７節で行われる混色の作成に使用されるモジュールである。ステップＳ１００においては、ＵＩＭ Ｐ２ａがＧＤＩ Ｐ１ａおよびＤＤＶ Ｐ５を介して、見本チャートＳＣを印刷させる旨の印刷指示を受付けるためのＵＩ画面Ｘを表示させる。ＵＩ画面Ｘにおいては、見本チャートＳＣのテンプレートを示す表示が設けられている。
【００１６】
図４は、ＵＩ画面Ｘの一例を示している。同図において、前記テンプレートＴＰが表示されており、当該テンプレートＴＰにはカラーパッチをレイアウトするための１２個の枠ＦＬ１〜ＦＬ１２が設けられている。ＵＩ画面Ｘには各枠ＦＬ１〜ＦＬ１２をマウス５０ｂのクリックによって選択することが可能となっており、枠ＦＬ１〜ＦＬ１２をクリックすると分光反射率測定を開始させるか否かを指示するための選択ウィンドウＷが表示される。また、ＵＩ画面Ｘにおいては、見本チャートＳＣの印刷を実行させるか否かを指示するためのボタンＢ１も設けられている。ステップＳ１１０においては、ＵＩＭ Ｐ２ａがマウス５０ｂによる各枠ＦＬ１〜ＦＬ１２のクリックを検出し、検出された場合にはステップＳ１２０にて分光反射率測定を開始させるか否かを指示するための選択ウィンドウＷを表示させる。ステップＳ１３０においては選択ウィンドウＷにおけるマウス５０ｂのクリックを検出し、キャンセルがクリックされた場合にはステップＳ１１０に戻る。一方、分光反射率測定実行がクリックされた場合には、ステップＳ１４０においてＭＣＭ Ｐ２ｂがＭＤＶ Ｐ４を介して分光反射率計３０にターゲットＴＧの分光反射率Ｒ（λ）としてターゲット分光反射率Ｒt（λ）の測定を実行させ、当該ターゲット分光反射率Ｒt（λ）を格納した分光反射率データＲＤを取得する。
【００１７】
ステップＳ１４０におけるターゲット分光反射率Ｒt（λ）の測定が完了すると、最も標準的な光源であるＤ６５光源を照射したときのＣＩＥＬＡＢ色空間における色彩値（Ｌ*ａ*ｂ*値）を算出する。そして、当該Ｌ*ａ*ｂ*値を所定のＲＧＢプロファイルを使用してＲＧＢ値に変換し、当該ＲＧＢ値を表示用ＲＧＢ値として取得する。なお、ＲＧＢプロファイルは絶対色空間としてのＣＩＥＬＡＢ色空間と本実施形態のＲＧＢ色空間との等色関係を規定したプロファイルであり、例えばＩＣＣプロファイルを使用することができる。
【００１８】
図５は、ステップＳ１４０において、分光反射率データＲＤから表示用ＲＧＢ値を算出する様子を模式的に示している。ターゲットＴＧについてターゲット分光反射率Ｒt（λ）を測定した結果、図示するようターゲット分光反射率Ｒt（λ）の分布を示す分光反射率データＲＤが得られたとする。なお、ターゲットＴＧは、分光的な再現の目標とする物体表面を意味し、例えば他の印刷装置や塗装装置等で形成した人工的物体表面や、自然物の表面等が該当する。一方、Ｄ６５光源は図示するような可視波長域において不均一な分光エネルギーＰ（λ）の分布を有しており、ターゲットＴＧにＤ６５光源を照射したときの各波長の反射光の分光エネルギーは、ターゲット分光反射率Ｒt（λ）と分光エネルギーＰ（λ）を各波長について掛け合わせた値となる。さらに、反射光の分光エネルギーのスペクトルに対して人間の分光感度特性に応じた等色関数ｘ（λ），ｙ（λ），ｚ（λ）をそれぞれ畳み込み積分し、係数ｋによって正規化することにより、３刺激値Ｘ，Ｙ，Ｚを得る。以上を数式で表すと下記の（１）式となる。
【数１】
2009177789.tif 000003
【００１９】
３刺激値Ｘ，Ｙ，Ｚを所定の変換式によって変換することにより、ターゲットＴＧにＤ６５光源を照射したときの色を示すＬ*ａ*ｂ*値を得ることができ、さらにＲＧＢプロファイルを使用することにより、表示用ＲＧＢ値を得ることができる。ステップＳ１４５においては、テンプレートＴＰにおいてクリックされた枠ＦＬ１〜ＦＬ１２を前記表示用ＲＧＢ値で塗りつぶした表示に更新する。これにより、標準的な光源であるＤ６５光源でのターゲットＴＧの色をＵＩ画面にて感覚的に把握することが可能となる。ステップＳ１４５が完了すると、ステップＳ１５０において固有のインデックスを生成するとともに、当該インデックスと、ステップＳ１１０にてクリックされた枠ＦＬ１〜ＦＬ１２の位置情報と表示用ＲＧＢ値を分光反射率データＲＤに対応付けてＲＡＭ１２に記憶する。ステップＳ１５０が完了すると、ステップＳ１１０に戻り、ステップＳ１２０〜Ｓ１５０を繰り返して実行する。これにより、他の枠ＦＬ１〜ＦＬ１２を選択し、他の枠ＦＬ１〜ＦＬ１２について他のターゲットＴＧのターゲット分光反射率Ｒt（λ）の測定を行うことができる。
【００２０】
本実施形態においては、それぞれ異なる１２種類のターゲットＴＧ１〜ＴＧ１２が用意されており、ターゲットＴＧ１〜ＴＧ１２のそれぞれについてのターゲット分光反射率Ｒt（λ）が分光反射率測データＲＤとして取得されるものとする。従って、ステップＳ１５０においては、各枠ＦＬ１〜ＦＬ１２についての分光反射率測データＲＤと固有のインデックスを対応付けたデータが順次ＲＡＭに記憶されていくこととなる。なお、インデックスはそれぞれの値が固有となるように生成されればよく、インクリメントによって生成してもよいし、重複しない乱数によって生成してもよい。
【００２１】
ステップＳ１１０において、各枠ＦＬ１〜ＦＬ１２のクリックが検出されない場合には、ステップＳ１６０にて見本チャートＳＣの印刷を実行させる旨のボタンＢ１のクリックを検出し、検出されない場合にはステップＳ１１０に戻る。一方、見本チャートＳＣの印刷を実行させる旨のボタンＢ１のクリックが検出された場合には、ステップＳ１７０にてＰＤＧ Ｐ２ｃが印刷データＰＤを生成する。
【００２２】
図６は、印刷データＰＤの構成を模式的に示している。同図において、印刷データＰＤはドットマトリクス状に配列した多数の画素によって構成されており、各画素が４バイト（８ビット×４）の情報を有している。印刷データＰＤは、図４に示したテンプレートＴＰと同様の画像を示しており、テンプレートＴＰの各枠ＦＬ１〜ＦＬ１２に対応する領域以外の画素は、テンプレートＴＰに対応する色のＲＧＢ値を有している。ＲＧＢ各チャネルの階調値はそれぞれ８ビット（２５６階調）によって表現され、上述した４バイトのうち３バイトがＲＧＢ値を格納するために使用される。例えば、テンプレートＴＰの各枠ＦＬ１〜ＦＬ１２以外の色が（Ｒ，Ｇ，Ｂ）＝（１２８，１２８，１２８）の一様な中間グレーで表される場合、印刷データＰＤにおける各枠ＦＬ１〜ＦＬ１２に対応する領域以外の画素は（Ｒ，Ｇ，Ｂ）＝（１２８，１２８，１２８）の色情報を有することとなる。なお、残りの１バイトは使用されない。
【００２３】
一方、テンプレートＴＰの各枠ＦＬ１〜ＦＬ１２に対応する画素も４バイトの情報を有しており、通常、ＲＧＢ値が格納される３バイトを使用してインデックスを格納する。このインデックスは、ステップＳ１５０にて各枠ＦＬ１〜ＦＬ１２ごとに生成した固有のものであり、ＰＤＧ Ｐ２ｃはインデックスをＲＡＭ１２から取得し、各枠ＦＬ１〜ＦＬ１２に対応する画素に対応するインデックスを格納する。このようにインデックスをＲＧＢ値の代わりに格納した各枠ＦＬ１〜ＦＬ１２に対応する画素については、残りの１バイトを使用してインデックスが格納された旨のフラグを立てる。これにより、各画素がＲＧＢ値を格納しているか、インデックスを格納しているかを判別することができる。本実施形態では、インデックスを格納するために３バイトを使用することができるため、３バイト以下の情報量で表現できるインデックスをステップＳ１５０で生成しておく必要がある。以上のようにしてビットマップ形式の印刷データＰＤが生成できると、ステップＳ１８０において、ＰＤＧ Ｐ２ｃがインデックステーブルＩＤＢを生成する。
【００２４】
図７は、インデックステーブルＩＤＢの一例を示している。同図において、各枠ＦＬ１〜ＦＬ１２に対応して生成された固有のインデックスのそれぞれに対して、測定によって得られたターゲット分光反射率Ｒt（λ）およびＤ６５光源におけるＬ*ａ*ｂ*値に対応する表示用ＲＧＢ値が格納されている。インデックステーブルＩＤＢの生成が完了すると、印刷データＰＤはＧＤＩ Ｐ１ａやスプーラＰ１ｂを経由してＰＤＶ Ｐ３ｂに出力される。印刷データＰＤは、外形上、通常のＲＧＢビットマップ形式と変わらないため、ＯＳ Ｐ１が提供するＧＤＩ Ｐ１ａやスプーラＰ１ｂにおいても通常の印刷ジョブと同様に処理することができる。一方、インデックステーブルＩＤＢは、直接、ＰＤＶ Ｐ３ｂに出力される。なお、本実施形態においては、インデックステーブルＩＤＢを新たに生成するようにしたが、既存のインデックステーブルＩＤＢにインデックスとターゲット分光反射率Ｒt（λ）と表示用ＲＧＢ値の新たな対応関係を追記するようにしてもよい。また、以上の印刷データ生成処理と後述する印刷制御処理は必ずしも同一の装置において連続して実行する必要はなく、印刷データ生成処理と印刷制御処理を例えばＬＡＮやインターネット等の通信回線によって接続された複数のコンピュータ上にて個別に実行してもよい。
【００２５】
３．混色の印刷データ生成処理
さらにＡＰＬ Ｐ２は混色の印刷データ生成も可能になっている。図８は混色の印刷データ生成部Ｐ２ｂが主体となって実行する混色の印刷データ生成処理の流れを示している。ステップＳ４００において、ＵＩＭ Ｐ２ａがＧＤＩ Ｐ１ａおよびＤＤＶ Ｐ５を介して、混色の作成と作成した混色を印刷させる旨の印刷指示を受付けるためのＵＩ画面Ｙを表示させる。
【００２６】
図９は、ＵＩ画面Ｙの一例を示している。同図においてＵＩ画面Ｙには、混色の元になる色を指定する指定枠ＦＬ２１〜ＦＬ２４が設けられている。指定枠ＦＬ２１〜ＦＬ２４はマウス５０ｂのクリックによって選択可能になっており、指定枠ＦＬ２１〜ＦＬ２４をクリックするとカラーパレットＣＰが新たなウィンドウとしてポップアップ表示される。カラーパレットＣＰには、混色の元となる色見本ＣＬ１〜ＣＬ１６が一覧表示されており、マウス５０ｂで色見本ＣＬ１〜ＣＬ１６のいずれかをクリックするとクリックされた色が指定枠に指定される。また、ＵＩ画面Ｙには、混色の割合を指定するスライダＳＬと、混色の予測を指示するボタンＢ２も設けられている。図９においては、一例として４つの指定枠を正方形の各角に表示してあり、４色までの混色の作成が可能なように記載してある。各指定枠の間にはスライダＳＬが表示されており、スライダＳＬの位置で各指定枠の色の使用割合（混色において各指定枠の色が使用される比率）が決定される。無論、２色以上であれば任意の色数の混色が可能である。
【００２７】
混色の印刷データ生成処理が開始されると、ステップＳ４００において、混色の指定を行うＵＩ画面Ｙを表示する。続くステップＳ４１０においては、混色の元となる色の指定を受付ける。具体的にはＵＩＭ Ｐ２ａがマウス５０ｂによる指定枠ＦＬ２１〜ＦＬ２４のいずれかのクリックを検出し、検出された場合はステップＳ４２０に進んでカラーパレットＣＰのウィンドウをポップアップ表示させる。そしてカラーパレットＣＰの色見本ＣＬ１〜ＣＬ１６の何れかのクリックを検出し、検出された色見本を受付け、受付けた色見本をステップＳ４００でクリックされた指定枠に対して設定してステップＳ４３０に進む。色見本が設定された指定枠には、この色見本が表示される。一方、ステップＳ４１０においてマウス５０ｂのクリックを検出しない場合は、ステップＳ４４０に進む。
【００２８】
続くステップＳ４３０においては、使用割合の指定を受付ける。具体的には、ＵＩＭ Ｐ２ａがマウス５０ｂによるスライダＳＬのドラッグアンドドロップを検出し、検出された場合にはドラッグアンドドロップ中のマウス５０ｂの動きに追随してスライダＳＬを移動させる。なお、より詳細な使用割合を指定するために、各見本色の使用割合を数値で入力する等、使用割合の指定方法は様々なものが採用可能である。
【００２９】
続くステップＳ４４０においては、指定枠ＦＬ２１〜ＦＬ２４のうち２枠以上に色見本が設定されているか否かを判断する。２つ以上の指定枠に色見本が設定されていない場合は、混色が不可能なのでステップＳ４１０に戻り、色見本が２枠以上に設定済みであればステップＳ４５０に進む。
【００３０】
ステップＳ４５０においては、混色の予測を開始するか否かを判断する。具体的には、ＵＩＭ Ｐ２ａがマウス５０ｂによるボタンＢ２のクリックを検出し、検出した場合にはステップＳ４６０に進んで混色の分光反射率の計算を実行する。マウス５０ｂによるボタンＢ２のクリックを検出しない場合は、ステップＳ４１０に戻る。
【００３１】
ステップＳ４６０においては、指定枠に指定された色見本をスライダＳＬにて指定された使用割合にて混合した混色の分光反射率の予測を行う。混色の分光反射率は、使用割合に応じて各色見本の分光反射率に対し重み付けを行った線型結合で計算したり、後述の４節に記載の分光プリンティングモデルにおけるノイゲバウアモデル等で計算したりすることができる。例えば、４色の色見本の混色を線型結合で計算する場合は、分光反射率Ｒ１（λ）の色見本と、分光反射率Ｒ２（λ）の色見本と、分光反射率Ｒ３（λ）の色見本と、分光反射率Ｒ４（λ）の色見本と、を、ｆ１：ｆ２：ｆ３：ｆ４（ｆ１＋Ｆ２＋ｆ３＋ｆ４＝１、０≦ｆ１≦１、０≦ｆ２≦１、０≦ｆ３≦１、０≦ｆ４≦１）、の割合で混合して作成される混色の分光反射率Ｒｍｉｘ（λ）は、Ｒｍｉｘ（λ）＝ｆ１×Ｒ１（λ）＋ｆ２×Ｒ２（λ）＋ｆ３×Ｒ３（λ）＋ｆ４×Ｒ４（λ）、の式にて計算できる。
【００３２】
また、同様に４色の色見本の混色を、分光プリンティングモデルのセル分割ユール・ニールセン分光ノイゲバウアモデルを利用して計算する場合は、４節の分光プリンティングモデルにおけるインクセット（ＣＭＹやＣＭＹＫｌｃｌｍ等）の代わりに指定枠ＦＬ２１〜ＦＬ２４に指定された色を使用し、インク量セットの代わりにスライダＳＬにて指定された使用割合を使用することになる。
【００３３】
なお、混色の予測は、前述のように複数の色見本の混色を予測するのみならず、他のプリンタにおける印刷結果を予測することも可能である。すなわち、実際に印刷実行するプリンタ２０以外の他のプリンタＰにおけるインクセットの組合せにより作成される色の分光反射率を予測するのである。より具体的には、プリンタＰのインクセットに基づいて後述の４節のように分光プリンティングモデルによる分光反射率データベースを構築し、プリンタＰの任意のインク量セットが入力された場合の分光反射率を予測するのである。このように他のプリンタＰにて使用されるインクセットや、混色の元となる色として前述したカラーパレットＣＰからＦＬにて受付けられた各色等が第２の色材群に相当する。一方、プリンタ２０のインクセットは、第１の色材群に相当する。
【００３４】
ステップ４６０における分光反射率Ｒｍｉｘ（λ）の算出が完了すると、ステップ４７０において、固有のインデックスを作成するとともに、当該インデックスを分光反射率データＲＤに対応付けてＲＡＭ１２に記憶する。また、分光反射率データＲＤに対して、最も標準的な光源であるＤ６５光源を照射したときのＣＩＥＬＡＢ色空間における色彩値（Ｌ*ａ*ｂ*値）を算出する。そして、当該Ｌ*ａ*ｂ*値を所定のＲＧＢプロファイルを使用してＲＧＢ値に変換し、当該ＲＧＢ値を表示用ＲＧＢ値として測色データＭＤと対応付けてＲＡＭ１２に記憶する。なお、ＲＧＢプロファイルは絶対色空間としてのＣＩＥＬＡＢ色空間と本実施形態のＲＧＢ色空間との等色関係を規定したプロファイルであり、例えばＩＣＣプロファイルを使用することができる。この分光反射率データＲＤから表示用ＲＧＢ値を算出する様子は前述の印刷データ生成処理と同様なので、説明を省略する。
【００３５】
続くステップＳ４８０においては、混色の印刷を実行させる旨のボタンＢ３のクリックを検出し、検出されない場合にはステップＳ４１０に戻る。一方、混色の印刷を実行させる旨のボタンＢ３のクリックが検出された場合には、ステップＳ４９０にてＰＤＧ Ｐ２ｃが印刷データＰＤを生成する。作成される印刷データは前述の印刷データ生成処理と同様である。印刷データＰＤが生成されると、ステップＳ５００において、ＰＤＧ Ｐ２ｃがインデックステーブルＩＤＢを生成する。生成されるインデックステーブルＩＤＢについても前述の印刷データ生成処理と同様である。インデックステーブルＩＤＢの生成が完了すると、印刷データＰＤはＧＤＩ Ｐ１ａやスプーラＰ１ｂを経由してＰＤＶ Ｐ３ｂに出力される。一方、インデックステーブルＩＤＢは、直接、ＰＤＶ Ｐ３ｂに出力される。
【００３６】
４．印刷制御処理
図１０は、ＬＵＧ Ｐ３ａとＰＤＶ Ｐ３ｂが実行する印刷制御処理の全体的な流れを示している。１Ｄ−ＬＵＴ生成処理（ステップＳ２００）をＬＵＧ Ｐ３ａが担当し、一方の印刷制御データ生成処理（ステップＳ３００）をＰＤＶ Ｐ３ｂが担当する。１Ｄ−ＬＵＴ生成処理は印刷制御データ生成処理に先行して行われてもよいし、１Ｄ−ＬＵＴ生成処理と印刷制御データ生成処理を並行して行うようにしてもよい。
【００３７】
４−１．１Ｄ−ＬＵＴ生成処理
図１１は、１Ｄ−ＬＵＴ生成処理の流れを示している。図２に示すようにＬＵＧ Ｐ３ａは、インク量セット算出モジュール（ＩＣＭ）Ｐ３ａ１と分光反射率予測モジュール（ＲＰＭ）Ｐ３ａ２と評価値算出モジュール（ＥＣＭ）Ｐ３ａ３とＬＵＴ出力モジュール（ＬＯＭ） Ｐ３ａ４とから構成されている。ステップＳ２１０においては、ＩＣＭ Ｐ３ａ１がインデックステーブルＩＤＢを取得する。ステップＳ２２０においては、インデックステーブルＩＤＢから一つのインデックスを選択し、当該インデックスに対応付けられている分光反射率データＲＤを取得する。ステップＳ２３０においては、ＩＣＭ Ｐ３ａ１が前記分光反射率データＲＤが示すターゲット分光反射率Ｒt（λ）や混色分光反射率Ｒｍｉｘ（λ）と同様の分光反射率Ｒ（λ）が再現可能なインク量セットを算出する処理を行う。その際に、上述したＲＰＭ Ｐ３ａ２とＥＣＭ Ｐ３ａ３を使用する。
【００３８】
図１２は、前記分光反射率データＲＤが示すターゲット分光反射率Ｒt（λ）や混色分光反射率Ｒｍｉｘ（λ）と同様の分光反射率Ｒ（λ）が再現可能なインク量セットを算出する処理の流れを模式的に示している。ＲＰＭ Ｐ３ａ２は、ＩＣＭ Ｐ３ａ１からのインク量セットφの入力に応じて、当該インク量セットφに基づいてプリンタ２０が所定の印刷用紙にインクを吐出させたときの分光反射率Ｒ（λ）を予測し、当該分光反射率Ｒ（λ）を予測分光反射率Ｒs（λ）としてＥＣＭ Ｐ３ａ３に出力する。
【００３９】
ＥＣＭ Ｐ３ａ３は、分光反射率データＲＤが示すターゲット分光反射率Ｒt（λ）や混色分光反射率Ｒｍｉｘ（λ）と予測分光反射率Ｒs（λ）の差分Ｄ（λ）を各波長λについて算出し、各波長λごとに重みが課せられた重み関数ｗ（λ）を当該差分Ｄ（λ）に乗算する。この値の二乗平均の平方根を評価値Ｅ（φ）として算出する。以上の計算を数式で表すと下記の（２）式のように表すことができる。
【数２】
2009177789.tif 000004
【００４０】
前記の（２）式において、Ｎは波長λの有限区分数を意味する。前記の（２）式において、評価値Ｅ（φ）が小さければ小さいほど、各波長λにおけるターゲット分光反射率Ｒt（λ）や混色分光反射率Ｒｍｉｘ（λ）と予測分光反射率Ｒs（λ）の差が少ないということができる。すなわち、評価値Ｅ（φ）が小さければ小さいほど、入力したインク量セットφによってプリンタ２０が印刷したときに記録媒体上にて再現される分光反射率Ｒ（λ）と、対応するターゲットＴＧから得られたターゲット分光反射率Ｒt（λ）や混色分光反射率Ｒｍｉｘ（λ）が近似しているということができる。さらに、前記の（１）式によれば、光源の変動に応じてインク量セットφによってプリンタ２０が印刷したときの記録媒体と対応するターゲットＴＧが示す絶対的な色彩値は双方とも変動するものの、分光反射率Ｒ（λ）が近似すれば光源の変動に拘わらず相対的には同じ色に知覚されるということができる。従って、評価値（φ）が小さくなるインク量セットφによれば、あらゆる光源においてターゲットＴＧと同じ色に知覚される印刷結果を得ることができるということができる。
【００４１】
また、本実施形態において、重み関数ｗ（λ）は下記の（３）式のものを使用する。
【数３】
2009177789.tif 000005
前記の（３）式においては、等色関数ｘ（λ），ｙ（λ），ｚ（λ）を加算することにより、重み関数ｗ（λ）が定義されている。なお、前記の（３）式の右辺全体に所定の係数を乗算して、重み関数ｗ（λ）の値の範囲を正規化してもよい。前記の（１）式によれば、等色関数ｘ（λ），ｙ（λ），ｚ（λ）が大きい波長域ほど、色彩値（Ｌ*ａ*ｂ*値）に大きく影響するということができる。従って、等色関数ｘ（λ），ｙ（λ），ｚ（λ）を加算した重み関数ｗ（λ）を使用すれば、色への影響が大きい波長域を重視した二乗誤差が評価可能な評価値Ｅ（φ）を得ることができる。例えば、人間の目に知覚されない近紫外波長域においてはｗ（λ）が０となり、当該波長域における差分Ｄ（λ）は評価値Ｅ（φ）の増大に寄与しないこととなる。
【００４２】
すなわち、必ずしも全可視波長域においてターゲット分光反射率Ｒt（λ）や混色分光反射率Ｒｍｉｘ（λ）と予測分光反射率Ｒs（λ）との差が小さくなくても、人間の目に特に強く知覚される波長域においてターゲット分光反射率Ｒt（λ）や混色分光反射率Ｒｍｉｘ（λ）と予測分光反射率Ｒs（λ）とが似ていれば小さい値の評価値Ｅ（φ）を得ることができ、人間の目に知覚に即した分光反射率Ｒ（λ）の近似性の指標として評価値Ｅ（φ）を使用することができる。算出された評価値Ｅ（φ）はＩＣＭ Ｐ３ａ１に返される。すなわち、ＩＣＭＰ ３ａ１が任意のインク量セットφをＲＰＭ Ｐ３ａ２とＥＣＭ Ｐ３ａ３に出力することにより、最終的に評価値Ｅ（φ）がＩＣＭ Ｐ３ａ１に返される構成となっている。ＩＣＭ Ｐ３ａ１は任意のインク量セットφに対応して評価値Ｅ（φ）を得ることを繰り返し実行することにより、目的関数としての評価値Ｅ（φ）が極小化するようなインク量セットφの最適解を算出する。この最適解を算出する手法としては、様々な最適化手法を用いることができるが、例えば勾配法といった非線形最適化手法を用いることが望ましい。
【００４３】
図１３は、ステップＳ２３０においてインク量セットφが最適化されていく様子を模式的に示している。同図において、インク量セットφが最適化されていくにつれて、インク量セットφで印刷を行った場合の予測分光反射率Ｒs（λ）がターゲット分光反射率Ｒt（λ）や混色分光反射率Ｒｍｉｘ（λ）に近づいていく。また、重み関数ｗ（λ）を使用することにより、等色関数ｘ（λ），ｙ（λ），ｚ（λ）が大きい波長域ほど、予測分光反射率Ｒs（λ）のターゲット分光反射率Ｒt（λ）や混色分光反射率Ｒｍｉｘ（λ）への拘束が強くなっており、予測分光反射率Ｒs（λ）がターゲット分光反射率Ｒt（λ）や混色分光反射率Ｒｍｉｘ（λ）の差が小さくなっている。このように、等色関数ｘ（λ），ｙ（λ），ｚ（λ）が大きく、視覚に大きく影響する波長域について優先的に予測分光反射率Ｒs（λ）をターゲットＴＧのターゲット分光反射率Ｒt（λ）や混色分光反射率Ｒｍｉｘ（λ）に拘束することができるため、任意の光源を照射したときの見た目が近くなるようなインク量セットφを算出することができる。以上により、いずれの光源においてもターゲットＴＧと似たような見た目をプリンタ２０にて再現させることが可能なインク量セットφを算出することができる。なお、最適化の終了条件は、インク量セットφ更新の繰り返し回数としてもよいし、評価値Ｅ（φ）の閾値としてもよい。
【００４４】
以上のようにして、ステップＳ２３０においてＩＣＭ Ｐ３ａ１がターゲットＴＧと同様の分光反射率Ｒ（λ）が再現可能なインク量セットφを算出すると、ステップＳ２４０においてインデックステーブルＩＤＢに記述されたインデックスのすべてがステップＳ２２０にて選択されたか否かを判定し、すべて選択されていない場合にはステップＳ２２０に戻り、次のインデックスを選択する。このようにすることにより、すべてのインデックスについてターゲットＴＧと同様の色が再現可能なインク量セットφを算出することができる。すなわち、印刷データ生成処理（図２）のステップＳ１４０において測色を行ったすべてのターゲットＴＧ１〜ＴＧ１２についてターゲットＴＧ１〜ＴＧ１２と同様の分光反射率Ｒ（λ）が再現可能なインク量セットφを算出することができる。ステップＳ２４０において、すべてのインデックスについて最適なインク量セットφが算出されたことが判定されると、ステップＳ２５０において、ＬＯＭ Ｐ３ａ４が１Ｄ−ＬＵＴを生成し、当該１Ｄ−ＬＵＴをＣＤＧ Ｐ３ｂに出力する。
【００４５】
図１４は、１Ｄ−ＬＵＴの一例を示している。同図において、各インデックスに対応して最適なインク量セットφが格納されている。すなわち、各ターゲットＴＧ１〜ＴＧ１２について、各ターゲットＴＧ１〜ＴＧ１２と似たような見た目をプリンタ２０にて再現させることが可能なインク量セットφを記述した１Ｄ−ＬＵＴを用意することができる。１Ｄ−ＬＵＴをＣＤＧ Ｐ３ｂに出力すると、１Ｄ−ＬＵＴ生成処理が完了し、次の印刷制御データ生成処理（ステップＳ３００）を実行させる。
【００４６】
以上説明したように、分光反射率予測の手法は、混色の分光反射率を予測する際の手法と、該混色の分光反射率を実現するインク量セットを計算する手法とで、同じ手法が利用可能である。すなわち、混色をプリンタ２０にて印刷するための色材料セットを予測するにあたり、第２の色材群の各色材の分光反射率と使用割合とに基づいて混色の分光反射率Ｒｍｉｘ（λ）を算出する手法と、プリンタ２０の各インクの分光反射率と各インクの使用割合と基づいて予測分光反射率ＲＳ（λ）を算出する手法と、に同じ手法を利用しつつ、予測分光反射率ＲＳ（λ）が混色分光反射率Ｒｍｉｘ（λ）に近似する分光反射率となるように、プリンタ２０における各インクの使用割合を変化させて繰り返し分光反射率予測を行うのである。なお、予測分光反射率ＲＳ（λ）と混色分光反射率Ｒｍｉｘ（λ）とが近似するか否かは、前述の評価関数等で評価すればよい。このように分光反射率の予測手法を共通化することにより、各分光反射率の予測の間で不整合が生じにくくなり、混色の予測がより適切なものになる。また、予測手法に使用するアルゴリズムが共通化されることで、プログラムサイズのコンパクト化にも寄与する。
【００４７】
４−２．印刷制御データ生成処理
図１５は、印刷制御データ生成処理の流れを示している。図２に示すようにＣＤＧ Ｐ３ｂは、モード判別モジュール（ＭＩＭ）Ｐ３ｂ１とインデックス分版モジュール（ＩＳＭ）Ｐ３ｂ２とＲＧＢ分版モジュール（ＣＳＭ）Ｐ３ｂ３とハーフトーンモジュール（ＨＴＭ）Ｐ３ｂ４とラスタ化モジュール（ＲＴＭ）Ｐ３ｂ５とから構成されている。ステップＳ３１０においては、モード判別モジュール（ＭＩＭ）Ｐ３ｂ１が印刷データＰＤを取得する。ステップＳ３２０において、ＭＩＭ Ｐ３ｂ１は印刷データＰＤから一の画素を選択する。ステップＳ３３０において、ＭＩＭ Ｐ３ｂ１は当該選択した画素にインデックスが格納された旨のフラグが立っているか否かを判定する。当該フラグが立っていないと判定された場合には、ステップＳ３４０にてＣＳＭ Ｐ３ｂ３が３Ｄ−ＬＵＴを参照して、当該画素についての色変換（分版）を実行する。
【００４８】
図１６は、３Ｄ−ＬＵＴの一例を示している。同図において、３Ｄ−ＬＵＴはＲＧＢ値とインク量セットφ（ｄＣ，ｄＭ，ｄＹ，ｄＫ，ｄｌｃ，ｄｌｍ）との対応関係が色空間における複数の代表的な座標について記述されたテーブルであり、ＣＳＭ Ｐ３ｂ３は３Ｄ−ＬＵＴを参照して当該画素が有するＲＧＢ値に対応するインク量セットφを取得する。その際に、３Ｄ−ＬＵＴに直接記述されていないＲＧＢ値については補間演算を行うことにより、対応するインク量セットφを取得する。なお、３Ｄ−ＬＵＴの作成方法として、特開２００６−８２４６０号公報等を採用することができる。当該公報においては、特定光源における色の再現性や、再現色の階調性や、粒状性や、再現色の光源非依存性や、ガマットや、インクデューティが総合的に良好となる３Ｄ−ＬＵＴが作成される。
【００４９】
一方、ステップＳ３３０において、選択した画素にインデックスが格納された旨のフラグが立っていると判定された場合には、ステップＳ３５０にてＩＳＭ Ｐ３ｂ２が１Ｄ−ＬＵＴを参照して、当該画素についての色変換（分版）を実行する。すなわち、インデックスが格納された旨のフラグが立っている画素から、インデックスを取得するとともに、１Ｄ−ＬＵＴにて当該インデックスに対応付けられているインク量セットφを取得する。ステップＳ３４０またはステップＳ３５０のいずれかにおいて、当該画素についてのインク量セットφが取得できると、ステップＳ３６０においてすべての画素についてインク量セットφが取得できたか否かを判定する。ここでインク量セットφが未取得の画素が残っている場合には、ステップＳ３２０に戻り次の画素を選択する。
【００５０】
以上の処理を繰り返して実行することにより、すべての画素についてインク量セットφを取得することができる。すべての画素についてインク量セットφが取得できると、すべての画素がインク量セットφで表現された印刷データＰＤに変換されたこととなる。以上のように各画素について１Ｄ−ＬＵＴと３Ｄ−ＬＵＴのいずれを使用するかを判定することにより、インデックスが格納された枠Ｆ１〜Ｆ１２に対応する画素については、各光源において各ターゲットＴＧ１〜ＴＧ１２に近い色が再現可能なインク量セットφを取得することができるとともに、ＲＧＢ値が格納された画素については３Ｄ−ＬＵＴの作成指針（例えば、粒状性を重視する等。）に基づく色再現が可能なインク量セットφを取得させることができる。
【００５１】
ステップＳ３７０においては、各画素がインク量セットφで表現された印刷データＰＤをＨＴＭ Ｐ３ｂ４が取得し、ハーフトーン処理を実行する。ＨＴＭ Ｐ３ｂ４はハーフトーン処理をするにあたっては公知のディザ法や誤差拡散法等を使用することができる。ハーフトーン処理が完了した印刷データＰＤにおいては、各画素が各インクを吐出させるか否を示す吐出信号を有している。ステップＳ３８０では、ハーフトーン処理が完了した印刷データＰＤをＲＴＭ Ｐ３ｂ５が取得し、当該印刷データＰＤにおける吐出信号をプリンタ２０が有する印刷ヘッドの各走査パスおよび各ノズルに割り振る処理を実行する。以上によりプリンタ２０に出力可能な印刷制御データＣＤが生成でき、プリンタ２０の制御に必要な信号を添付した印刷制御データＣＤをスプーラＰ１ｂおよびプリンタ２０に出力する。これによりプリンタ２０が印刷用紙上にインクを吐出して、見本チャートＳＣを形成する。
【００５２】
以上のようにして印刷用紙上に形成された見本チャートＳＣの枠ＦＬ１〜ＦＬ１２に対応する領域において、各ターゲットＴＧ１〜ＴＧ１２のターゲット分光反射率Ｒt（λ）を再現したり、混色の分光反射率Ｒｍｉｘ（λ）を再現したりすることができる。すなわち、枠ＦＬ１〜ＦＬ１２に対応する領域は、ターゲットＴＧ１〜ＴＧ１２の複数光源下での色に追従するようなインク量セットφで印刷されているため、各光源下においてターゲットＴＧ１〜ＴＧ１２と似たような色を再現することができる。例えば、見本チャートＳＣを室内で視認したときの各枠ＦＬ１〜ＦＬ１２に対応する領域の色はターゲットＴＧ１〜ＴＧ１２を室内で視認したときの色を再現するし、見本チャートＳＣを室外で視認したときの各枠ＦＬ１〜ＦＬ１２に対応する領域の色もターゲットＴＧ１〜ＴＧ１２を室外で視認したときの色を再現することができる。
【００５３】
なお、究極的には、ターゲットＴＧ１〜ＴＧ１２と完全に同じ分光反射率Ｒ（λ）を有する見本チャートＳＣや予測された混色の分光反射率と完全に同じ分光反射率を再現すれば、いかなる光源においてもターゲットＴＧ１〜ＴＧ１２や実際に混色された印刷結果と同様の色を再現することができる。しかしながら、プリンタ２０が使用可能なインク（色材の種類）はＣＭＹＫｌｃｌｍに限られているため、完全に同じ分光反射率Ｒ（λ）が再現可能なインク量セットφを求めることは実質的に不可能である。また、知覚色に影響しない波長域について完全な分光反射率Ｒ（λ）が再現可能なインク量セットφを求めても、視覚的な再現精度の実現においては無駄となる。これに対して、本発明では、等色関数ｘ（λ），ｙ（λ），ｚ（λ）に基づく重み付けを行った評価値Ｅ（φ）を利用してターゲット分光反射率Ｒt（λ）や予測された混色への近似性を評価しているため、視覚的に十分な精度が達成可能なインク量セットφを求めることができる。
【００５４】
一方、印刷用紙上に形成された見本チャートＳＣの枠ＦＬ１〜ＦＬ１２に対応する領域や混色の印刷される領域においては、上述した３Ｄ−ＬＵＴに基づいたインク量セットφによって印刷がなされることとなる。そのため、当該領域についての印刷パフォーマンスは３Ｄ−ＬＵＴに基づくものとなる。上述したとおり本実施形態においては枠ＦＬ１〜ＦＬ１２や混色が印刷される領域以外の領域は中間グレーの一様な画像を示すが、当該領域において３Ｄ−ＬＵＴが目標とする印刷パフォーマンスを満足させることができる。すなわち、再現色の階調性や、粒状性や、再現色の光源非依存性や、ガマットや、インクデューティが総合的に良好となる印刷を実現することができる。
【００５５】
５．分光プリンティングモデル
図１７は、本実施形態のプリンタ２０の印刷方式を模式的に示している。同図において、プリンタ２０は、ＣＭＹＫｌｃｌｍのインクごとに複数のノズル２１ａ，２１ａ・・・を備えた印刷ヘッド２１を備えており、ノズル２１ａ，２１ａ・・・が吐出するＣＭＹＫｌｃｌｍのインクごとのインク量を上述したインク量セットφ（ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlm）によって指定された量とする制御が印刷制御データＣＤに基づいて行われる。各ノズル２１ａ，２１ａ・・・が吐出したインク滴は印刷用紙上において微細なドットとなり、多数のドットの集まりによってインク量セットφ（ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlm）に応じたインク被覆率の印刷画像が印刷用紙上に形成されることとなる。
【００５６】
ＲＰＭ Ｐ３ａ２が使用する予測モデル（分光プリンティングモデル）は、本実施形態のプリンタ２０で使用され得る任意のインク量セットφ（ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlm）で印刷を行った場合の分光反射率Ｒ（λ）を予測分光反射率Ｒs（λ）として予測するための予測モデルである。分光プリンティングモデルにおいては、インク量空間における複数の代表点について実際にカラーパッチを印刷し、その分光反射率Ｒ（λ）を分光反射率計によって測定することにより得られた分光反射率データベースＲＤＢを用意する。そして、この分光反射率データベースＲＤＢを使用したセル分割ユール・ニールセン分光ノイゲバウアモデル(Cellular Yule-Nielsen Spectral Neugebauer Model)による予測を行うことにより、正確に任意のインク量セットφ（ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlm）で印刷を行った場合の分光反射率Ｒ（λ）を予測する。
【００５７】
図１８は、分光反射率データベースＲＤＢを示している。同図に示すように分光反射率データベースＲＤＢはインク量空間（本実施形態では６次元であるが、図の簡略化のためＣＭ面のみ図示。）における複数の格子点のインク量セットφ（ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlm）について実際に印刷／測定をして得られた分光反射率Ｒ（λ）が記述されたルックアップテーブルとなっている。例えば、各インク量軸を分割する５グリッドの格子点を発生させる。ここでは５13個もの格子点が発生し、膨大な量のカラーパッチの印刷／測定をすることが必要となるが、実際にはプリンタ２０にて同時に搭載可能なインク数や同時に吐出可能なインクデューティの制限があるため、印刷／測定をする格子点の数は絞られることとなる。
【００５８】
さらに、一部の格子点のみ実際に印刷／測定をし、他の格子点については実際に印刷／測定を行った格子点の分光反射率Ｒ（λ）に基づいて分光反射率Ｒ（λ）を予測することにより、実際に印刷／測定を行うカラーパッチの個数を低減させてもよい。分光反射率データベースＲＤＢは、プリンタ２０が印刷可能な印刷用紙ごとに用意されている必要がある。厳密には、分光反射率Ｒ（λ）は印刷用紙上に形成されたインク膜（ドット）による分光透過率と印刷用紙の反射率によって決まるものであり、印刷用紙の表面物性（ドット形状が依存）や反射率の影響を大きく受けるからである。次に、分光反射率データベースＲＤＢを使用したセル分割ユール・ニールセン分光ノイゲバウアモデルによる予測を説明する。
【００５９】
ＲＰＭ Ｐ３ａ２は、ＩＣＭ Ｐ３ａ１の要請に応じて分光反射率データベースＲＤＢを使用したセル分割ユール・ニールセン分光ノイゲバウアモデルによる予測を実行する。この予測にあたっては、ＩＣＭ Ｐ３ａ１から予測条件を取得し、この予測条件を設定する。具体的には、印刷用紙やインク量セットφを印刷条件として設定する。例えば、光沢紙を印刷用紙として予測を行う場合には、光沢紙にカラーパッチを印刷することにより作成した分光反射率データベースＲＤＢが設定される。
【００６０】
分光反射率データベースＲＤＢの設定ができると、ＩＣＭ Ｐ３ａ１から入力されたインク量セットφ（ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlm）を分光プリンティングモデルに適用する。セル分割ユール・ニールセン分光ノイゲバウアモデルは、よく知られた分光ノイゲバウアモデルとユール・ニールセンモデルとに基づいている。なお、以下の説明では、説明の簡略化のためＣＭＹの３種類のインクを用いた場合のモデルについて説明するが、同様のモデルを本実施形態のＣＭＹＫｌｃｌｍを含む任意のインクセットを用いたモデルに拡張することは容易である。また、セル分割ユール・ニールセン分光ノイゲバウアモデルについては、Color Ｒes Appl 25, 4-19, 2000およびＲ Balasubramanian, Optimization of the spectral Neugebauer model for printer characterization, J. Electronic Imaging 8(2), 156-166 (1999)を参照。
【００６１】
図１９は、分光ノイゲバウアモデルを示す図である。分光ノイゲバウアモデルでは、任意のインク量セット（ｄc，ｄm，ｄy）で印刷したときの印刷物の予測分光反射率Ｒs（λ）は、以下の（４）式で与えられる。
【数４】
2009177789.tif 000006
【００６２】
ここで、ａiはｉ番目の領域の面積率であり、Ｒi（λ）はｉ番目の領域の分光反射率である。添え字ｉは、インクの無い領域（ｗ）と、シアンインクのみの領域（ｃ）と、マゼンタインクのみの領域（ｍ）と、イエローインクのみの領域（ｙ）と、マゼンタインクとイエローインクが吐出される領域（ｒ）と、イエローインクとシアンインクが吐出される領域（ｇ）と、シアンインクとマゼンタインクが吐出される領域（ｂ）と、ＣＭＹの３つのインクが吐出される領域（ｋ）をそれぞれ意味している。また、ｆc，ｆm，ｆyは、ＣＭＹ各インクを１種類のみ吐出したときにそのインクで覆われる面積の割合（「インク被覆率(Ink area coverage)」と呼ぶ）である。
【００６３】
インク被覆率ｆc，ｆm，ｆyは、図１９（Ｂ）に示すマーレイ・デービスモデルで与えられる。マーレイ・デービスモデルでは、例えばシアンインクのインク被覆率ｆcは、シアンのインク量ｄcの非線形関数であり、例えば１次元ルックアップテーブルによってインク量ｄcをインク被覆率ｆcに換算することができる。インク被覆率ｆc，ｆm，ｆyがインク量ｄc，ｄm，ｄyの非線形関数となる理由は、単位面積に少量のインクが吐出された場合にはインクが十分に広がるが、多量のインクが吐出された場合にはインクが重なり合うためにインクで覆われる面積があまり増加しないためである。他の種類のＭＹインクについても同様である。
【００６４】
分光反射率に関するユール・ニールセンモデルを適用すると、前記（４）式は以下の（５ａ）式または（５ｂ）式に書き換えられる。
【数５】
2009177789.tif 000007
ここで、ｎは１以上の所定の係数であり、例えばｎ＝１０に設定することができる。前記の（５ａ）式および（５ｂ）式は、ユール・ニールセン分光ノイゲバウアモデル(Yule-Nielsen Spectral Neugebauer Model)を表す式である。
【００６５】
本実施形態で採用するセル分割ユール・ニールセン分光ノイゲバウアモデル(Cellular Yule-Nielsen Spectral Neugebauer Model)は、上述したユール・ニールセン分光ノイゲバウアモデルのインク量空間を複数のセルに分割したものである。
【００６６】
図２０（Ａ）は、セル分割ユール・ニールセン分光ノイゲバウアモデルにおけるセル分割の例を示している。ここでは、説明の簡略化のために、ＣＭインクのインク量ｄc，ｄmの２つの軸を含む２次元インク量空間でのセル分割を描いている。なお、インク被覆率ｆc，ｆmは上述したマーレイ・デービスモデルにてインク量ｄc，ｄmと一意の関係にあるため、インク被覆率ｆc，ｆmを示す軸と考えることもできる。白丸は、セル分割のグリッド点（「格子点」と呼ぶ）であり、２次元のインク量（被覆率）空間が９つのセルＣ１〜Ｃ９に分割されている。各格子点に対応するインク量セット（ｄc，ｄm）は、分光反射率データベースＲＤＢに規定された格子点に対応するインク量セットとされている。すなわち、上述した分光反射率データベースＲＤＢを参照することにより、各格子点の分光反射率Ｒ（λ）を得ることができる。従って、各格子点の分光反射率Ｒ（λ）00，Ｒ（λ）10，Ｒ（λ）20・・・Ｒ（λ）33は、分光反射率データベースＲＤＢから取得することができる。
【００６７】
実際には、本実施形態ではセル分割もＣＭＹＫｌｃｌｍの６次元インク量空間で行うとともに、各格子点の座標も６次元のインク量セットφ（ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlm）によって表される。そして、各格子点のインク量セットφ（ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlm）に対応する格子点の分光反射率Ｒ（λ）が分光反射率データベースＲＤＢ（例えば光沢紙のもの）から取得されることとなる。
【００６８】
図２０（Ｂ）は、セル分割モデルにて使用するインク被覆率ｆcとインク量ｄcとの関係を示している。ここでは、１種類のインクのインク量の範囲０〜ｄcmaxも３つの区間に分割されており、各区間毎に０から１まで単調に増加する非線形の曲線によってセル分割モデルにて使用する仮想的なインク被覆率ｆcが求められる。他のインクについても同様にインク被覆率ｆm，ｆyが求められる。
【００６９】
図２０（Ｃ）は、図２０（Ａ）の中央のセルＣ５内にある任意のインク量セット（ｄc，ｄm）にて印刷を行った場合の予測分光反射率Ｒs（λ）の算出方法を示している。インク量セット（ｄc，ｄm）にて印刷を行った場合の分光反射率Ｒ（λ）は、以下の（６）式で与えられる。
【数６】
2009177789.tif 000008
ここで、（６）式におけるインク被覆率ｆc，ｆmは図２０（Ｂ）のグラフで与えられる値である。また、セルＣ５を囲む４つの格子点に対応する分光反射率Ｒ（λ）11，（λ）12，（λ）21，（λ）22は分光反射率データベースＲＤＢを参照することにより取得することができる。これにより、（６）式の右辺を構成するすべての値を確定することができ、その計算結果として任意のインク量セットφ（ｄc，ｄm）にて印刷を行った場合の予測分光反射率Ｒs（λ）を算出することができる。波長λを可視波長域にて順次シフトさせていくことにより、可視波長域における予測分光反射率Ｒs（λ）を得ることができる。インク量空間を複数のセルに分割すれば、分割しない場合に比べて予測分光反射率Ｒs（λ）をより精度良く算出することができる。以上のようにして、ＲＰＭ Ｐ３ａ２がＩＣＭ Ｐ３ａ１の要請に応じて予測分光反射率Ｒs（λ）を予測することができる。
【００７０】
６．変形例
６−１：変形例１
以下の変形例においては、ターゲット分光反射率Ｒt（λ）を例にとって説明を行うが、混色の分光反射率Ｒｍｉｘ（λ）についても同様である。図２１は、変形例においてＥＣＭ Ｐ３ａ３が設定する重み関数ｗ（λ）を模式的に示している。同図において、ターゲットＴＧから得られたターゲット分光反射率Ｒt（λ）が示されており、当該ターゲット分光反射率Ｒt（λ）と各等色関数ｘ（λ），ｙ（λ），ｚ（λ）との相関係数ｃx，ｃy，ｃzがＥＣＭ Ｐ３ａ３によって算出されている。そして、下記の（７）式によって本変形にかかる重み関数ｗ（λ）を算出する。
【数７】
2009177789.tif 000009
前記の（７）式においては、ターゲットＴＧから得られたターゲット分光反射率Ｒt（λ）との相関が高い等色関数ｘ（λ），ｙ（λ），ｚ（λ）ほど線形結合の際の重みが大きくなるようにされている。以上のようにして得られた重み関数ｗ（λ）においては、ターゲットＴＧのターゲット分光反射率Ｒt（λ）が大きい波長域についての重みを強調することができる。従って、各光源下での反射光の分光エネルギーのスペクトルが強くなりがちな波長域を重視した評価値Ｅ（φ）を得ることができる。すなわち、特にターゲットＴＧのターゲット分光反射率Ｒt（λ）が大きい波長域については、ターゲットＴＧのターゲット分光反射率Ｒt（λ）と予測分光反射率Ｒs（λ）とのずれを許容しないようなインク量セットφの最適解を得ることができる。むろん、重み関数ｗ（λ）は各等色関数ｘ（λ），ｙ（λ），ｚ（λ）に由来しているため、人間の知覚に適合した評価値Ｅ（φ）を得ることができる。
【００７１】
６−２：変形例２
図２２は、別の変形例においてＥＣＭ Ｐ３ａ３が設定する重み関数ｗ（λ）を模式的に示している。同図において、ターゲットＴＧから得られたターゲット分光反射率Ｒt（λ）をそのまま重み関数ｗ（λ）として適用している。このようにすることによっても、特にターゲットＴＧのターゲット分光反射率Ｒt（λ）が大きい波長域についてターゲットＴＧの分光反射率Ｒ（λ）とターゲット分光反射率Ｒt（λ）とのずれを許容しないようなインク量セットφの最適解を得ることができる。
【００７２】
６−３：変形例３
図２３は、別の変形例においてＥＣＭ Ｐ３ａ３が設定する重み関数ｗ（λ）を模式的に示している。同図において、５種類（標準昼光系のＤ５０光源，Ｄ５５光源，Ｄ６５光源、白熱電球系のＡ光源、蛍光ランプ系のＦ１１光源）の各光源の分光エネルギーＰ D50（λ），Ｐ D55（λ），Ｐ D65（λ），ＰA（λ），Ｐ F11（λ）が示されている。本変形例においては、下記の（８）式によって、これらの分光エネルギーＰ D50（λ），Ｐ D55，Ｐ D65（λ），ＰA（λ），Ｐ F11（λ）を線形結合することにより重み関数ｗ（λ）を算出する。
【数８】
2009177789.tif 000010
【００７３】
前記の（８）式において、ｗ1〜ｗ5は各光源についての重みを設定する重み係数である。このように、光源の分光エネルギー分布Ｐ D50（λ），Ｐ D55（λ），Ｐ D65（λ），ＰA（λ），Ｐ F11（λ）に由来する重み関数ｗ（λ）を設定することにより、各光源下での反射光の分光エネルギーのスペクトルが強くなりがちな波長域を重視した評価値Ｅ（φ）を得ることができる。また、重み係数ｗ1〜ｗ5を調整することも可能である。例えば全光源における色の再現性をバランスよく確保したい場合にはｗ1＝ｗ2＝ｗ3＝ｗ4＝ｗ5とすればよいし、人工光源における色の再現性を重視したい場合にはｗ1，ｗ2，ｗ3＜ｗ4，ｗ5とすればよい。
【００７４】
６−４：変形例４
図２４は、変形例においてディスプレイ４０に表示されるＵＩ画面を示している。同図において、ＵＩ画面において複数のターゲット分光反射率Ｒt（λ）のグラフが表示されている。このようなＵＩ画面を表示させることにより、ユーザーがステップＳ１４０においてターゲットＴＧのターゲット分光反射率Ｒt（λ）を測定する代わりに、所望の波形のグラフをターゲットＴＧのターゲット分光反射率Ｒt（λ）として選択することができる。このようにすることにより、実際に分光反射率測定をしなくてもターゲット分光反射率Ｒt（λ）を設定することができる。むろん、グラフの波形をユーザーが直接編集できるようにしてもよい。例えば、新規の物体表面の開発を行う際に目標とするターゲット分光反射率Ｒt（λ）に編集しておけば、実際に物体表面を試作することなく、目標とするターゲット分光反射率Ｒt（λ）を有する見本チャートＳＣをプリンタ２０によって印刷させることができる。このような所望の波形のグラフは、混色の元となる見本色にも利用可能であり、所望のグラフにて表現される分光反射率を混色した色を予測可能となる。
【００７５】
６−５：変形例５
図２５は、変形例にかかる評価値（φ）を模式的に説明している。同図において、ターゲットＴＧのターゲット分光反射率Ｒt（λ）に対して前記５種類の光源を照射したときの色彩値（ターゲット色彩値）を上述した（１）式，図５によって算出する。一方、ＲＰＭ Ｐ３ａ２が予測した予測分光反射率Ｒs（λ）に対して前記５種類の光源を照射したときの色彩値（予測色彩値）も上述した（１）式（Ｒt（λ）をＲs（λ）に置き換えて使用），図５によって算出する。そして、各光源におけるターゲット色彩値と予測色彩値の色差ΔＥ（ΔＥ2000）をＣＩＥ  ＤＥ２０００の色差式に基づいて算出する。そして、各光源についての色差ΔＥをΔＥ D50，ΔＥ D55，ΔＥ D65，ΔＥA，ΔＥ F11とし、下記の（９）式によって評価値Ｅ（φ）を算出する。
【数９】
2009177789.tif 000011
前記の（２）式において、ｗ1〜ｗ5は各光源の重みを設定する重み係数であり、上述した変形例３の重み係数ｗ1〜ｗ5とほぼ同様の性質を有する。ここでも全光源における色の再現性をバランスよく確保したい場合にはｗ1＝ｗ2＝ｗ3＝ｗ4＝ｗ5とすればよいし、人工光源における色の再現性を重視したい場合にはｗ1，ｗ2，ｗ3＜ｗ4，ｗ5とすればよい。
【００７６】
６−６．変形例６
図２６〜図２７は、本発明の変形例にかかる印刷システムのソフトウェア構成を示している。図２４に示すように、上述した実施形態のＬＵＧ Ｐ３ａに相当する構成がＰＤＶ Ｐ３ｂの内部モジュール（１Ｄ−ＬＵＴ作成部）として備えられていてもよい。また、図２７に示すように、上述した実施形態のＬＵＧ Ｐ３ａに相当する構成が他のコンピュータ１１０において実行されてもよい。この場合、コンピュータ１０とコンピュータ１１０とが所定の通信インターフェイスＣＩＦによって接続され、コンピュータ１１０のＬＵＧ Ｐ３ａにて生成された１Ｄ―ＬＵＴが通信インターフェイスＣＩＦを介してコンピュータ１０に送信される。通信インターフェイスＣＩＦは、インターネットを介在するものであってもよい。その場合、コンピュータ１０はインターネット上のコンピュータ１１０から取得した１Ｄ−ＬＵＴを参照して色変換を行うことができる。さらに、プリンタ２０が各ソフトウェア構成Ｐ１〜Ｐ５を実行するようにしてもよい。むろん、ソフトウェア構成Ｐ１〜Ｐ５と同等の処理を実行するハードウェアがプリンタ２０に組み込まれる場合にも、本発明を実現することができる。
【００７７】
６−７．変形例７
前記実施形態では、複数色を指定してその混色を表示するＵＩについて説明したが、以上説明した混色の予測手法は、ユーザの意図する色（目標色）を再現するためにも利用可能である。例えば、ユーザが所望の色（目標色）をイメージしているにも関わらずその目標色がインクセットにて表現可能な色見本（代表色）の中に無く、目標色を有する印刷物等が手元に無いときに、目標色そのものではないが近い色の実物やデータを持っている場合や、近い色を指定するために必要な情報を持っている場合等に、これら近い色を指標（指針）色として目標色を探索することができる。
【００７８】
より具体的には、例えば、ユーザがイメージしている色（目標色ＴＣ）を再現するために混色する色見本の組み合わせを探索するにあたり、目標色ＴＣに類似する色（指標色ＩＣ）の指定をユーザから受付けて、まずは指標色ＩＣを再現するために必要な色材組合せとその混合割合とを求める。このようにして再現された指標色ＩＣは、目標色ＴＣに類似した色である。従って、指標色ＩＣにおける色材組合せの混色割合を変更することにより、より目標色ＴＣに近い色を再現するための混合割合を求めることが出来る可能性がある。そこで、指標色ＩＣの混色割合を変更する手段と、変更された混色割合の色を表示する手段と、をユーザに提供し、指標色ＩＣを起点としてユーザが混色割合を変更しつつ最もユーザの目標色ＴＣに近い色を再現できるようにする。すなわちユーザに混色割合を変更する自由度を与えることにより、ユーザの感性に基づいて目標色ＴＣにより近い色を再現する色材を探索可能にするのである。以下、指標色ＩＣとして画像データの一部の色を利用する場合を例にとって説明を行う。
【００７９】
図２８は、混色の印刷データ生成部Ｐ２ｂが主体となって実行する目標色探索処理のフローチャート、図２９は、指標色ＩＣの選択をユーザから受付けるためのユーザインターフェース（ＵＩ）の一例である。同図に示すＵＩ Ｗは、画像表示枠Ｗ１、指標色選択枠Ｗ２、混色割合指定部Ｗ３を備えており、混色割合指定部Ｗ３は第１近似色枠Ｗ３ａと第２近似色枠Ｗ３ｂと混色表示枠Ｗ３ｃと指標色表示枠Ｗ３ｄとを備える。画像表示枠Ｗ１は、ユーザから指定された画像データに基づく画像表示をするエリアである。指標色選択枠Ｗ２は、画像表示枠Ｗ１に表示された画像の中の一部を囲う選択枠であり、所定の操作入力が行われると選択枠内の色を指標色ＩＣとして設定する。指標色選択枠Ｗ２にて設定された指標色は、指標色表示枠Ｗ３ｃに表示される。なお、指標色選択枠Ｗ２の中の色が単一色でなくて複数色を含んでいる場合は、枠内の色を平均した色を指標色ＩＣとして設定する。
【００８０】
混色割合指定部Ｗ３は、指標色ＩＣを再現するための混色元となる色材の色である第１近似色と第２近似色とを表示する第１近似色枠Ｗ３ａと第２近似色枠Ｗ３ｂを所定間隔離して表示し、これら枠をバーで接続する。このバー上を移動可能なスライダとして、各近似色の混色割合に基づいた位置に混色表示枠Ｗ３ｃを表示する。この混色表示枠Ｗ３ｃには、その位置に対応する混色割合で各近似色を混是合わせた色が表示される。指標色表示枠Ｗ３ｄには、指標色選択枠Ｗ２で設定された指標色ＩＣが表示される。この指標色表示枠Ｗ３ｄと混色表示枠Ｗ３ｃは、並べて配置されており、混色と指標色ＩＣとを見比べながら混色割合を変更し、目標色ＴＣに近い色探索できるようになっている。なお、本変形例では、２色を混色して目標色を探索する例について説明するが、無論、３色以上の複数色を混色して目標色を探索してもよい。例えば３色を混色するのであれば、三角形の頂点に混色元となる各近似色を表示して、その三角形の内部に混色表示枠を表示する。そして、三角形の内部で混色表示枠を移動可能とし、前記実施形態の４色の混色と同様にその位置に応じた混色を算出して混色表示枠に表示することができる。
【００８１】
なお、前実施形態においては、ターゲット分光反射率Ｒt（λ）を実際に測定し、当該ターゲット分光反射率Ｒt（λ）とインデックスとを対応付けたインデックステーブルが用意されるようにしたが、予め多数のインデックスとターゲット分光反射率Ｒt（λ）とが登録されたインデックステーブルが用意されていてもよい。そこで、本変形例では、塗料メーカーが製造する各塗料について与えられたインデックスと、各塗料の塗布表面を測定したターゲット分光反射率Ｒt（λ）との対応関係を登録したインデックステーブルが予め用意されているものとする。当該インデックステーブルにおいては、上述した実施形態と同様に表示用ＲＧＢ値も登録されている。インデックステーブルが予め用意されている場合、ステップＳ１００では、見本チャートＳＣに再現したい塗料（インデックス）を選択するための処理がＡＰＬ Ｐ２によって行われる。
【００８２】
図２８の処理が開始されると、ステップＳ５０５において、図２９のＵＩ Ｗを表示する。そしてステップＳ５１０において、指標色ＩＣの選択を受付ける。より具体的には、例えば、ユーザがＵＩの画像表示枠Ｗ１にマウスのポインタを合わせてクリックする操作を行うと、画像表示枠Ｗ１に表示させる画像データの選択を受付ける画像選択ウィンドウ等が表示される。画像選択ウィンドウにおいて指標色ＩＣを含んだ画像データのファイルをユーザが選択すると、この選択された画像データを受付けて画像データに基づく画像を画像表示枠Ｗ１に表示する。次に、画像表示枠Ｗ１の内部に指標色表示枠Ｗ２が表示される。指標色表示枠Ｗ２は、カーソル移動キーの操作やマウスの動きと連動して画像表示枠Ｗ１の内部で移動する。ユーザは、指標色表示枠Ｗ２の内部にユーザのイメージに最も近い色が含まれるようにカーソル移動キーやマウスを操作し、決定キーやマウスの右クリックを行って指標色表示枠Ｗ２の色を指標色ＩＣに設定する。なお、指標色表示枠Ｗ２のサイズを変更できるようにしたり、指標色表示枠Ｗ２に囲まれた各画素の色を平均した色を別途表示する枠を設けたりすると、ユーザは適切な指標色ＩＣを選択しやすくなる。このようにして指標色表示枠Ｗ２を利用して指定された指標色ＩＣは、ＲＧＢ値（指標色ＲＧＢ値）として取得される。指標色の設定が終了すると、ステップＳ５１０は条件成立となり、ステップＳ５２０に進む。
【００８３】
ステップＳ５２０においては、指標色ＩＣを再現するための近似色を選定する処理を行う。本処理においては、まず、インデックステーブルに登録されている色の中から、指標色ＩＣを中心として所定範囲内にある色を抽出する。指標色ＩＣから大きく離れた色を利用して混色しても指標色を得られない可能性があるからである。インデックステーブルには前述したように表示用ＲＧＢ値が設定されている。そこで、例えば、指標色ＲＧＢ値をＰ１＝（Ｒ０，Ｇ０，Ｂ０）とすると、以下の条件を満たすＲＧＢ値をインデックステーブルから抽出する。
Ｒ０−Δｒ≦Ｒ≦Ｒ０＋Δｒ
Ｇ０−Δｇ≦Ｇ≦Ｇ０＋Δｇ
Ｂ０−Δｂ≦Ｂ≦Ｂ０＋Δｂ
上記条件式において、Δｒ、Δｇ、Δｂは、色空間におけるインデックスの密度に応じて決定される定数である。すなわち、前記条件式にて規定される範囲内に、所定数以上の表示用ＲＧＢ値が含まれるように決定されている。含まれるインデックス数が少なすぎると、含まれたインデックスの任意の２色を結ぶ線と指標色ＩＣとがかけ離れてしまう可能性があるからである。
【００８４】
なお、近似色の範囲指定を行うにあたり、明度・彩度・色相等の色属性のいずれかを優先するように範囲指定を行ってもよい。すなわち、ユーザが用意した指標色と想定している目標色との間で最も類似する色属性を指定する操作入力をユーザから受付けて、指定された属性を変更しないようなインデックスが近似色の候補に選定されるように範囲指定する。例えば、指標色と目標色との色合いが似ているが明るさや鮮やかさがずれているのであれば、色相に対する限定を強くして範囲指定する。
【００８５】
図３０は色属性に基づいて範囲指定を行う方法を説明する図である。同図に示す色見本の分布を得てその範囲指定を行うためには、まず、前述したインデックスの表示用ＲＧＢ値と指標色ＲＧＢ値とを、例えば、公知の変換式でＨＳＶ色空間（明度・彩度・色相を変数とする色空間）に変換する。そして、ＨＳＶ空間において、指標色ＲＧＢ値を変換したＨＳＶ値（Ｑ０）の色相角（Ｈ値）に対して±５度以内の断面扇形状の空間を特定する。すなわち、色相角が指標色に近似する空間が特定される。次に、ＨＳＶ空間において、指標色ＲＧＢ値を変換したＨＳＶ値（Ｑ０）の明度Ｖおよび彩度Ｓに対して±５以内となる明度Ｖおよび彩度Ｓとなる円環状空間を特定する。すなわち、明度Ｖおよび彩度Ｓが指標色に近似する空間が特定される。このように指定された範囲に含まれるインデックスは、色属性が指標色と近似するものとなる。
【００８６】
範囲指定が完了すると、抽出されたＲＧＢ値の中から、指標色を再現するのに最も適した近似色の組合せを選定する。この選定においては、例えば、抽出された色点から任意の２色（説明上、各色の座標をＸ，Ｙとする。）を選択し、線分ＸＹと指標色ＲＧＢ値Ａとのユークリッド距離ｄ０が最小値を取る組合せを第１近似色ＡＣ１と第２近似色ＡＣ２に選定する。このように第１近似色ＡＣ１と第２近似色ＡＣ２を結ぶ線分上に近い位置に指標色ＩＣを位置させることにより、指標色ＩＣにより近い色を再現可能な近似色の組合せを選定できる。
【００８７】
なお、選定に際して、距離ｄ０のみならず、色点Ｘと指標色ＲＧＢ値との距離ｄ１や色点Ｙと指標色ＲＧＢ値との距離ｄ２を考慮して近似色ＲＧＢ値を選定してもよい。その理由は、２色を結ぶ軸上に存在する色点であっても色空間で距離が離れていると２色を混色して再現される色よりも、軸がずれていても色点同士が近い２色を混色した方が指標色の再現性が高い場合があるからである。具体的には、距離ｄ０に重みｗ１を乗じ、距離ｄ１、ｄ２に重みｗ２を乗じて、これらの和を評価式ｗ１Ｚ＋ｗ２（ｄ１＋ｄ２）として、この式が極小化する組合せを近似色ＲＧＢ値に選定する等である。
【００８８】
ステップＳ５３０においては、以上のステップＳ５１０にて指定された指標色ＩＣとステップＳ５２０にて選定された近似色ＡＣ１，ＡＣ２を、図２９に示したＵＩの下部に表示する。すなわち、指標色表示枠Ｗ３ｄに指標色ＩＣを表示し、第１近似色枠に近似色ＡＣ１を表示し、第２近似色枠Ｗ３ｂに近似色ＡＣ２を表示する。そして、近似色ＡＣ１，ＡＣ２の混色結果を混色表示枠Ｗ３ｃに表示する。なお、初期状態でスライダに表示される混色は、近似色を半々で混色した結果を表示してもよいし、指標色ＲＧＢ値と各近似色ＲＧＢ値との距離に応じた割合で近似色を混色して表示してもよい。
【００８９】
ステップＳ５４０においては、混合割合の変更を指示する操作入力が行われたか否かを判断する。該操作入力が行われた場合はステップＳ５６０に進み、該操作入力が検出されなかった場合はステップＳ５８０に進む。図２９のＵＩにおいては、混合割合の変更は、スライダをバー上で移動させることにより行われるようになっており、第１近似色に近づくほど第１近似色の混合比率を高め、第２近似色に近づくほど第２近似色の混色比率を高める。例えば、［スライダと第１近似色の間隔］：［スライダと第２近似色の間隔］＝ａ：ｂ、であれば第１近似色の混色割合をｂ／（ａ＋ｂ）、第２近似色の混色割合をａ／（ａ＋ｂ）、とすることが考えられるが、無論これに限定されるものではない。
【００９０】
ステップＳ５６０においては、指定された混色割合の色をスライダに表示する。すなわち、第１近似色ＡＣ１と第２近似色ＡＣ２とをスライダにて指定された使用割合にて混合した混色の分光反射率Ｒ’ｍｉｘ（λ）の予測を行う。混色の分光反射率は、使用割合に応じて各色見本の分光反射率Ｒ’ｍｉｘ（λ）に対し重み付けを行った線型結合で計算したり、後述の４節に記載の分光プリンティングモデルにおけるノイゲバウアモデル等で計算したりすることができる。例えば、第１近似色ＡＣ１の分光反射率Ｒ１（λ）と、第２近似色ＡＣ２の分光反射率Ｒ２（λ）と、を、ｆ１：ｆ２（ｆ１＋ｆ２＝１、０≦ｆ１≦１、０≦ｆ２≦１）、の割合で混合して作成される混色の分光反射率Ｒ’ｍｉｘ（λ）は、Ｒ’ｍｉｘ（λ）＝ｆ１×Ｒ１（λ）＋ｆ２×Ｒ２（λ）、なる式にて計算できる。
【００９１】
また、同様に第１近似色と第２近似色の混色を、分光プリンティングモデルのセル分割ユール・ニールセン分光ノイゲバウアモデルを利用して計算する場合は、４節の分光プリンティングモデルにおけるインクセット（ＣＭＹやＣＭＹＫｌｃｌｍ等）の代わりに第１近似色と第２近似色を使用し、インク量セットの代わりにスライダにて指定された使用割合を使用することになる。
【００９２】
ステップ５６０における分光反射率Ｒ’ｍｉｘ（λ）の算出が完了すると、ステップ５７０において、固有のインデックスを作成するとともに、当該インデックスを分光反射率データＲＤに対応付けてＲＡＭ１２に記憶する。また、分光反射率データＲＤに対して、最も標準的な光源であるＤ６５光源を照射したときのＣＩＥＬＡＢ色空間における色彩値（Ｌ*ａ*ｂ*値）を算出する。そして、当該Ｌ*ａ*ｂ*値を所定のＲＧＢプロファイルを使用してＲＧＢ値に変換し、当該ＲＧＢ値を表示用ＲＧＢ値として測色データＭＤと対応付けてＲＡＭ１２に記憶する。なお、ＲＧＢプロファイルは絶対色空間としてのＣＩＥＬＡＢ色空間と本実施形態のＲＧＢ色空間との等色関係を規定したプロファイルであり、例えばＩＣＣプロファイルを使用することができる。この分光反射率データＲＤから表示用ＲＧＢ値を算出する様子は前述の印刷データ生成処理と同様なので、説明を省略する。
【００９３】
続くステップＳ５８０においては、混色の印刷を実行させる旨のボタンＷ４のクリックを検出し、検出されない場合にはステップＳ４１０に戻る。一方、混色の印刷を実行させる旨のボタンＷ４のクリックが検出された場合には、ステップＳ５９０にてＰＤＧ Ｐ２ｃが印刷データＰＤを生成する。作成される印刷データは前述の印刷データ生成処理と同様である。印刷データＰＤが生成されると、ステップＳ６００において、ＰＤＧ Ｐ２ｃがインデックステーブルＩＤＢを生成する。生成されるインデックステーブルＩＤＢについても前述の印刷データ生成処理と同様である。インデックステーブルＩＤＢの生成が完了すると、印刷データＰＤはＧＤＩ Ｐ１ａやスプーラＰ１ｂを経由してＰＤＶ Ｐ３ｂに出力される。一方、インデックステーブルＩＤＢは、直接、ＰＤＶ Ｐ３ｂに出力される。
【図面の簡単な説明】
【００９４】
【図１】印刷制御装置のハードウェア構成を示すブロック図である。
【図２】印刷制御装置のソフトウェア構成を示すブロック図である。
【図３】印刷データ生成処理の流れを示すフローチャートである。
【図４】ＵＩ画面Ｘの一例を示す図である。
【図５】分光反射率に基づいて色彩値を算出するための計算を説明する図である。
【図６】印刷データを示す図である。
【図７】インデックステーブルを示す図である。
【図８】混色の印刷データ生成処理の流れを示すフローチャートである。
【図９】ＵＩ画面Ｙの一例を示す図である。
【図１０】印刷制御処理の全体的な流れを示すフローチャートである。
【図１１】１Ｄ−ＬＵＴ生成処理の流れを示すフローチャートである。
【図１２】インク量セットを最適化する処理の流れを示す模式図である。
【図１３】インク量セットが最適化されていく様子を示す模式図である。
【図１４】１Ｄ−ＬＵＴを示す図である。
【図１５】印刷制御データ生成処理の流れを示すフローチャートである。
【図１６】３Ｄ−ＬＵＴを示す図である。
【図１７】プリンタの印刷方式を示す模式図である。
【図１８】分光反射率データベースを示す図である。
【図１９】分光ノイゲバウアモデルを示す図である。
【図２０】セル分割ユール・ニールセン分光ノイゲバウアモデルを示す図である。
【図２１】変形例にかかる重み関数を示す模式図である。
【図２２】変形例にかかる重み関数を示す模式図である。
【図２３】変形例にかかる重み関数を示す模式図である。
【図２４】変形例にかかるＵＩ画面を示す図である。
【図２５】変形例にかかる評価値を示す模式図である。
【図２６】変形例にかかる印刷システムのソフトウェア構成である。
【図２７】変形例にかかる印刷システムのソフトウェア構成である。
【図２８】変形例にかかる目標色探索処理のフローチャートである。
【図２９】変形例にかかるユーザインターフェース（ＵＩ）の一例である。
【図３０】変形例にかかる色属性に基づいて範囲指定を行う方法を説明する図である。
【符号の説明】
【００９５】
１０…コンピュータ、１１…ＣＰＵ、１２…ＲＡＭ、１３…ＲＯＭ、１４…ＨＤＤ、１５…ＧＩＦ、１６…ＶＩＦ、１７…ＩＩＦ、１８…バス、Ｐ１…ＯＳ、Ｐ１ａ…ＧＤＩ、Ｐ１ｂ…スプーラ、Ｐ２…ＡＰＬ、Ｐ２ａ…ＵＩＭ、Ｐ２ｂ…ＭＣＭ、Ｐ２ｃ…ＰＤＧ、Ｐ３ａ…ＬＵＧ、Ｐ３ｂ…ＰＤＶ、Ｐ３ａ１…ＩＣＭ、Ｐ３ａ２…ＲＰＭ、Ｐ３ａ３…ＥＣＭ、Ｐ３ａ４…ＬＯＭ、Ｐ４…ＭＤＶ、Ｐ５…ＤＤＶ。
【図１】
2009177789.tif 000012
【図２】
2009177789.tif 000013
【図３】
2009177789.tif 000014
【図４】
2009177789.tif 000015
【図５】
2009177789.tif 000016
【図６】
2009177789.tif 000017
【図７】
2009177789.tif 000018
【図８】
2009177789.tif 000019
【図９】
2009177789.tif 000020
【図１０】
2009177789.tif 000021
【図１１】
2009177789.tif 000022
【図１２】
2009177789.tif 000023
【図１３】
2009177789.tif 000024
【図１４】
2009177789.tif 000025
【図１５】
2009177789.tif 000026
【図１６】
2009177789.tif 000027
【図１７】
2009177789.tif 000028
【図１８】
2009177789.tif 000029
【図１９】
2009177789.tif 000030
【図２０】
2009177789.tif 000031
【図２１】
2009177789.tif 000032
【図２２】
2009177789.tif 000033
【図２３】
2009177789.tif 000034
【図２４】
2009177789.tif 000035
【図２５】
2009177789.tif 000036
【図２６】
2009177789.tif 000037
【図２７】
2009177789.tif 000038
【図２８】
2009177789.tif 000039
【図２９】
2009177789.tif 000040
【図３０】
2009177789.tif 000041