(57)【要約】【課題】少ない印刷／測定によって正確に印刷結果を予測する。【解決手段】印刷において記録媒体に付着される複数の色材の使用量を軸とする色材量空間において、前記軸の上に予め選択された複数の一次色に基づいて前記色材量空間に配置される複数の直交格子点について、色材打込量制限に対する大小関係および各格子点間での大小関係を維持しつつ、各色材の使用量を単独で増減させたグラデーションに基づいて印刷したときの参照測定値と、前記複数の格子点に対応する前記参照測定値を前記予測手段に入力したときの予測値と、の誤差が極小化する位置へ変位させた格子点を選択し、各格子点の色材量セットで印刷を行った結果を測定した測定値に基づいて任意の色材使用量で印刷したときの印刷結果を予測する。【選択図】図１４
2010147783.tif 000002
【特許請求の範囲】
【請求項１】
印刷において記録媒体に付着される複数の色材の使用量を軸とする色材量空間の中から複数の格子点を選択し、各格子点の色材量セットで印刷を行った結果を測定した測定値に基づいて任意の色材使用量で印刷したときの印刷結果を予測する予測手段を用いて行う印刷結果予測方法であって、
前記複数の格子点は、前記軸の上に予め選択された複数の一次色に基づいて前記色材量空間に配置される複数の直交格子点について、色材打込量制限に対する大小関係および各格子点間での大小関係を維持しつつ、各色材の使用量を単独で増減させたグラデーションに基づいて印刷したときの参照測定値と、前記複数の格子点に対応する前記参照測定値を前記予測手段に入力したときの予測値と、の誤差が極小化する位置へ変位させた格子点であることを特徴とする印刷結果予測方法。
【請求項２】
前記グラデーションは、色彩値の変化が一定となるように前記色材の使用量を増減させてある請求項１に記載の印刷結果予測方法。
【請求項３】
前記軸の上に予め選択される複数の一次色には、前記色材の使用量の最大値と最小値とを含んでおり、これらの一次色に基づいて前記色材量空間に配置される複数の格子点を固定して前記誤差を極小化させる請求項１または請求項２に記載の印刷結果予測方法。
【請求項４】
各々ｍ個の一次色が前記軸上に選択された２つの色材量軸を含んで形成される色材量平面において、前記複数の直交格子点のうち、一方の軸において原点からｎ番目の一次色と他方の軸において原点からｍ−ｎ番目の一次色との直交格子点と、一方の軸において原点からｍ−ｎ番目の一次色と他方の軸において原点からｎ番目の一次色との直交格子点と、が前記色材打込量制限の内外に分かれないように選択される請求項１〜請求項３の何れか１項に記載の印刷結果予測方法。（ｍとｎは正の整数）
【請求項５】
前記色材打込量制限内には前記複数の直交格子点のうち、所定数が含まれる請求項１〜または請求項４のいずれか１項に記載の印刷結果予測方法。
【請求項６】
請求項１から請求項５のいずれか１項に記載の印刷結果予測方法によって予測された予測結果に基づいて色変換を行うことにより印刷装置に印刷を実行させることを特徴とする印刷制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は印刷結果予測方法および印刷制御装置に関し、特に色材を記録媒体に付着させて再現される色を予測する印刷結果予測方法および印刷制御装置に関する。
【背景技術】
【０００２】
従来、印刷装置が印刷可能な色の一部についてのみ印刷と測色を行い、その測定結果に基づいて印刷装置が印刷可能な色全体についての測定結果を予測することが行われている。例えば、セル分割ユール・ニールセン分光ノイゲバウアモデル（予測モデルの一種）に基づいて、任意のインク量での分光反射率を予測する手法が提案されている。かかる構成においては、インク量空間を分光反射率が既知の格子点（インク量（被覆率）空間において均等に分布する格子点）によって囲まれた複数のセルに分割し、セルを構成する各格子点の分光反射率に基づいてセル内の任意のインク量の分光反射率を予測する。インク量空間において比較的近い位置に存在する格子点の分光反射率に基づいて任意のインク量の分光反射率を予測することができるため、良好な予測精度を実現することができる。
【特許文献１】国際公開ＷＯ２００５／０４３８８４号のパンフレット
【発明の開示】
【発明が解決しようとする課題】
【０００３】
近年、家庭用プリンタの高画質化にともなって使用できるインクの種類が増加しており、上述したセル分割ユール・ニールセン分光ノイゲバウアモデルにおけるインク量空間も多次元化してきている。これに伴い、上述した予測モデルに利用される格子点の個数は、インク種の増加に応じてべき乗的に増加する。上述した予測モデルに利用される格子点についてはカラーパッチを印刷して測色する必要があるが、大量のカラーパッチを印刷／測色するには、膨大な時間を要するため、例えば変換プロファイルの作成等の作業に要する期間が長期化するという問題があった。これに対して、上述した予測モデルに利用する格子点の数を減少させることにより、変換プロファイルの作成等の作業に要する時間や手間を抑制することができるが、良好な予測精度を確保することができなくなるという問題があった。
上記課題を解決するために、本発明は、印刷結果を予測するモデルにおいて印刷／測定を必要とするカラーパッチ数を減らしつつ、高い精度の印刷結果を予測可能にすることを目的とする。
【課題を解決するための手段】
【０００４】
本発明の印刷結果予測方法においては、以下に説明する予測手段を用いて色の予測を行う。当該予測手段においては、まず印刷において記録媒体に付着させる複数の色材の使用量を軸とする色材量空間のなかから複数の格子点を選択する。そして、当該予測手段は、選択した各格子点の色材量セットで印刷を行ったときの測定値を取得し、当該取得した各測定値に基づいて任意の色材量セットで印刷したときの印刷結果を予測する。すなわち、前記予測手段は、一部の前記色材量セットについて予め測定値を得ておき、当該測定値に基づいて任意の前記色材量セットで印刷したときの印刷結果を予測する。
【０００５】
以上のような予測手段を用いた本発明の予測方法において、まず各色材の使用量を単独で増減させたグラデーションに基づいて印刷したときの参照測定値を取得する。そして、各色材の使用量の軸上に、複数の一次色を選択する。この一次色に基づいて前記色材量空間に複数の格子点が配置される。このように配置された複数の格子点は、前記参照測定値と、前記複数の格子点に対応する前記参照測定値を前記予測手段に入力したときの予測値と、の誤差が極小化する位置へ変位される。その際、色材打込量制限に対する大小関係と、各格子点間での大小関係とを維持するように、制約条件を課した上で変位を行う。従って、仮設定した時点での意図を反映しつつ実測値と実測値との誤差を少なくし、精度の高い印刷結果の予測が可能となる。
【０００６】
また、前記軸の上に予め選択される複数の一次色には、前記色材の使用量の最大値と最小値とを含んでおり、これらの一次色に基づいて前記色材量空間に配置される複数の格子点を固定して前記誤差を極小化させるようにしてもよい。前記予測の精度を向上させるに当り、色材量軸の全範囲にわたって予測精度を向上させるためには、予測元となる格子点が色材量軸の全範囲に分布している必要がある。すなわち、色材量軸の端部に格子点を確実に配置することにより、予測精度が向上する。
【０００７】
また、各々ｍ個の一次色が前記軸上に選択された２つの色材量軸を含んで形成される色材量平面において、前記複数の直交格子点のうち、一方の軸において原点からｎ番目の一次色と他方の軸において原点からｍ−ｎ番目の一次色との直交格子点と、一方の軸において原点からｍ−ｎ番目の一次色と他方の軸において原点からｎ番目の一次色との直交格子点と、が前記色材打込量制限の内外に分かれないように選択されるように構成してもよい。（ｍとｎは正の整数）このように原点側からの順番と最大値側からの順番とが転置した関係にある直交格子点間で、色材打込量制限に対する大小関係が異なる場合は、前記偏差を極小化した結果、インク量空間において直交格子点が偏在してしまう可能性がある。直交格子点が偏在すると予測精度が低下するため、このような状況を排除することによって、予測精度を向上することができる。
【０００８】
また、前記総色材打込量制限内には前記複数の直交格子点のうち、所定数が含まれるようにしてもよい。すなわち、所望の予測制度を得るために必要な格子点数を予め決定しておき、必要最小限の直交格子点で印刷結果を予測するようにすれば、上述した工数を最小限に止める事ができる。
【０００９】
さらに、以上説明した印刷結果予測方法を利用して色変換プロファイルを作成することも可能である。すなわち、本発明の印刷結果予測方法で予測した印刷結果とインク量セットとの対応関係を規定したルックアップテーブル等の色変換プロファイルを作成することができる。この場合も、色変換プロファイルの作成に要する作業を軽減することができる。なお、色変換プロファイルを使用して、さらに別の変換プロファイルを作成するようにしてもよい。さらに、当該色変換プロファイルを利用して色変換を実行する印刷制御装置において本発明を実現することも可能である。さらに、以上のようにして作成された色変換プロファイルを用いて印刷を行う印刷制御装置によれば、画像データの意味する色と該画像データに基づいて印刷した印刷結果の色との誤差が、所定の色材の使用量が単調に変化するときに、複数のピークを持ち、少なくとも一部のピークは、互いに略一定となる印刷を実現することができる。
【００１０】
なお、本発明の技術的思想は、方法のみならず、当該方法を実行するコンピュータ等のハードウェアや、上記方法に沿った処理手順をコンピュータ等のハードウェアに実行させるプログラムとしても具体的に実現可能なことは言うまでもない。また、本発明の方法は、単体として存在するものに限られず、ある方法の一部として組み込まれる場合もある。例えば、本発明の手法により印刷結果を予測する方法の工程を一部に備えた印刷方法や画像処理方法においても本発明が実現できることはいうまでもない。また、本発明を構成する各工程が複数の実体的な装置において分散して実行される工程によって実現されるものであってもよい。例えば、本発明の一部の工程がパーソナルコンピュータにて実現され、他の工程が印刷装置にて実現されるものであってもよい。むろん、本発明の各工程がネットワークを介して分散していてもよい。
【発明を実施するための最良の形態】
【００１１】
次に、本発明の実施の形態を以下の順序で説明する。
Ａ．印刷結果予測装置：
Ａ−１．分光反射率と色の予測：
Ａ−２．分光反射率データベースの作成：
Ｂ．色変換ルックアップテーブルの作成：
Ｃ．印刷制御装置（処理）：
Ｄ．変形例：
【００１２】
Ａ．印刷結果予測装置
図１は、本実施形態にかかる印刷結果予測装置のハードウェア構成を示している。同図において、印刷結果予測装置は主にコンピュータ１０によって構成されており、コンピュータ１０はＣＰＵ１１とＲＡＭ１２とＲＯＭ１３とハードディスクドライブ（ＨＤＤ）１４と汎用インターフェイス（ＧＩＦ）１５とビデオインターフェイス（ＶＩＦ）１６と入力インターフェイス（ＩＩＦ）１７とバス１８とから構成されている。バス１８は、コンピュータ１０を構成する各要素１１〜１７の間でのデータ通信を実現するものであり、図示しないチップセット等によって通信が制御されている。ＨＤＤ１４には、オペレーティングシステム（ＯＳ）を含む各種プログラムを実行するためのプログラムデータ１４ａが記憶されており、当該プログラムデータ１４ａをＲＡＭ１２に展開しながらＣＰＵ１１が当該プログラムデータ１４ａに準じた演算を実行する。ＧＩＦ １５は、例えばＵＳＢ規格に準じたインターフェイスを提供するものであり、外部のプリンタ２０と分光反射率計３０をコンピュータ１０に接続させている。ＶＩＦ １６はコンピュータ１０を外部のディスプレイ４０に接続し、ディスプレイ４０に画像を表示するためのインターフェイスを提供する。ＩＩＦ １７はコンピュータ１０を外部のキーボード５０ａとマウス５０ｂに接続し、キーボード５０ａとマウス５０ｂからの入力信号をコンピュータ１０が取得するためのインターフェイスを提供する。
【００１３】
図２は、コンピュータ１０にて実行されるプログラムのソフトウェア構成を示している。同図において、コンピュータ１０では、おもにＯＳ_Ｐ１とアプリケーション（ＡＰＬ）Ｐ２とプリンタドライバ（ＰＤＶ）Ｐ３と分光反射率計ドライバ（ＭＤＶ）Ｐ４とディスプレイドライバ（ＤＤＶ）Ｐ５が実行されている。ＡＰＬ_Ｐ２は、印刷結果としての分光反射率および色を予測するためのアプリケーションプログラムである。ＭＤＶ_Ｐ４は、ＡＰＬ_Ｐ２の要求に応じて分光反射率計３０を制御し、当該制御によって得られた分光反射率のデータをＡＰＬ_Ｐ２に出力する。ＰＤＶ_Ｐ３は、ＡＰＬ_Ｐ２の要求に応じてプリンタ２０を制御し、当該制御によって各種パッチを印刷させる。以下、フローチャートを用いて各プログラムＰ１〜Ｐ４が実行する処理を詳細に説明する。
【００１４】
Ａ−１．分光反射率と色の予測
図２に示すように、ＡＰＬ_Ｐ２は、大きく分光反射率予測部（ＲＰＭ）Ｐ２ａと、色予測部（ＣＰＭ）Ｐ２ｂと、データベース作成部（ＤＢＧ）Ｐ２ｃと、プロファイル作成部（ＰＦＧ）Ｐ２ｄとから構成されている。ＲＰＭ_Ｐ２ａは、予めＤＢＧ_Ｐ２ｃがＨＤＤ１４に作成した分光反射率データベースＲＤＢを使用して任意のインク量セットで印刷したときの分光反射率Ｒ（λ）を予測する。ＣＰＭ_Ｐ２ｂは、ＲＰＭ_Ｐ２ａが予測した分光反射率Ｒ（λ）に基づいて所定光源における色を予測する。なお、ＲＰＭ_Ｐ２ａとＣＰＭ_Ｐ２ｂ、および、ＲＰＭ_Ｐ２ａとＣＰＭ_Ｐ２ｂを実行させるハードウェア資源が、本発明の予測手段に相当する。
【００１５】
図３は本実施形態にかかる印刷結果予測処理の流れを示し、図４は分光反射率データベースＲＤＢを示している。ステップＳ１００においては、ＲＰＭ_Ｐ２ａがＨＤＤ１４から分光反射率データベースＲＤＢを取得する。分光反射率データベースＲＤＢは、本実施形態のプリンタ２０が印刷に使用するＣＭＹＫｌｃｌｍインクの各インク量ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlmの組合せであるインク量セットφと、当該インク量セットφに基づいてプリンタ２０が印刷を実際に行ったときの分光反射率Ｒ（λ）との対応関係を記述したテーブルである。なお、インク量セットφは６次元のインク量空間における任意の座標によって表現することができ、その組合せの数は膨大なものとなる。
【００１６】
図５は、インク量空間を示している。なお、図示の便宜上、ＹＫｌｃｌｍインクのインク量ｄy，ｄk，ｄlc，ｄlmを一定としたＣＭインク量平面を図示している。本実施形態においては、各インク量ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlmが８ビット（２５６階調）のデジタル階調で表され、階調が小さい方がインク量の少量側に対応し、階調が大きい方がインク量の多量側に対応している。このように、各インク量ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlmが８ビットで表現した場合には、インク量セットφの組合せは、２５６6個となり、膨大な数となる。従って、インク量セットφのうち一部の代表的なものを格子点（○で図示）として選択し、当該格子点についてのインク量セットφと分光反射率Ｒ（λ）の対応関係のみを分光反射率データベースＲＤＢに記憶させることとしている。分光反射率データベースＲＤＢの作成処理の詳細については後述する。
【００１７】
分光反射率データベースＲＤＢが取得できると、ステップＳ１１０にて、ＲＰＭ_Ｐ２ａは分光反射率Ｒ（λ）と色を予測したいインク量セットφを取得する。予測対象のインク量セットφは、例えばユーザから入力されてもよいし、後述するプロファイル作成部（ＰＦＧ）Ｐ２ｄから入力されてもよい。ここでは、任意のインク量セットφが取得される。ステップＳ１２０においては、ＲＰＭ_Ｐ２ａは、予測したいインク量セットφを、以下に説明する分光プリンティングモデルを適用することによって分光反射率Ｒ（λ）を予測する。
【００１８】
図６は、本実施形態のプリンタ２０の印刷方式を模式的に示している。同図において、プリンタ２０は、ＣＭＹＫｌｃｌｍのインクごとに複数のノズル２１ａ，２１ａ・・・を備えた印刷ヘッド２１を備えており、ノズル２１ａ，２１ａ・・・が吐出するＣＭＹＫｌｃｌｍのインクごとのインク量を上述したインク量セットφによって指定された量とする制御が印刷データに基づいて行われる。各ノズル２１ａ，２１ａ・・・が吐出したインク滴は印刷用紙上において微細なドットとなり、多数のドットの集まりによってインク量セットφに応じたインク被覆率の印刷画像が印刷用紙上に形成されることとなる。
【００１９】
ＲＰＭ_Ｐ２ａが使用する予測モデル（分光プリンティングモデル）は、本実施形態のプリンタ２０で使用され得る任意のインク量セットφで印刷を行った場合の分光反射率Ｒ（λ）を予測分光反射率Ｒs（λ）として予測するための予測モデルである。分光プリンティングモデルにおいては、インク量空間における複数の格子点について実際にカラーパッチを印刷し、その分光反射率Ｒ（λ）を分光反射率計によって測定することにより得られた分光反射率データベースＲＤＢを予め取得しておく（ステップＳ１００）。そして、この分光反射率データベースＲＤＢを使用したセル分割ユール・ニールセン分光ノイゲバウアモデル(Cellular Yule-Nielsen Spectral Neugebauer Model)による予測を行うことにより、正確に任意のインク量セットφで印刷を行った場合の分光反射率Ｒ（λ）を予測する。
【００２０】
以下、セル分割ユール・ニールセン分光ノイゲバウアモデルによる予測を説明する。イセル分割ユール・ニールセン分光ノイゲバウアモデルは、よく知られた分光ノイゲバウアモデルとユール・ニールセンモデルとに基づいている。なお、以下の説明では、説明の簡略化のためＣＭＹの３種類のインクを用いた場合のモデルについて説明するが、同様のモデルを本実施形態のＣＭＹＫｌｃｌｍを含む任意のインクセットを用いたモデルに拡張することは容易である。また、セル分割ユール・ニールセン分光ノイゲバウアモデルについては、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（λ）は、以下の（１）式で与えられる。
【数１】
2010147783.tif 000003
【００２２】
ここで、ａiはｉ番目の領域の面積率であり、Ｒi（λ）はｉ番目の領域の分光反射率である。添え字ｉは、インクの無い領域（ｗ）と、シアンインクのみの領域（ｃ）と、マゼンタインクのみの領域（ｍ）と、イエローインクのみの領域（ｙ）と、マゼンタインクとイエローインクが吐出される領域（ｒ）と、イエローインクとシアンインクが吐出される領域（ｇ）と、シアンインクとマゼンタインクが吐出される領域（ｂ）と、ＣＭＹの３つのインクが吐出される領域（ｋ）をそれぞれ意味している。また、ｆc，ｆm，ｆyは、ＣＭＹ各インクを１種類のみ吐出したときにそのインクで覆われる面積の割合（「インク被覆率(Ink area coverage)」と呼ぶ）である。
【００２３】
インク被覆率ｆc，ｆm，ｆyは、図７（Ｂ）に示すマーレイ・デービスモデルで与えられる。マーレイ・デービスモデルでは、例えばシアンインクのインク被覆率ｆcは、シアンのインク量ｄcの非線形関数であり、例えば１次元ルックアップテーブルによってインク量ｄcをインク被覆率ｆcに換算することができる。インク被覆率ｆc，ｆm，ｆyがインク量ｄc，ｄm，ｄyの非線形関数となる理由は、単位面積に少量のインクが吐出された場合にはインクが十分に広がるが、多量のインクが吐出された場合にはインクが重なり合うためにインクで覆われる面積があまり増加しないためである。他の種類のＭＹインクについても同様である。
【００２４】
分光反射率に関するユール・ニールセンモデルを適用すると、前記（１）式は以下の（２ａ）式または（２ｂ）式に書き換えられる。
【数２】
2010147783.tif 000004
ここで、ｎはユール・ニールセンのｎ値といい、「ドットゲイン」という現象を補正するものである。このｎ値は用紙特性や解像度により決定される１以上の所定の係数であり、予め適切な値を求めておく。例えばｎ＝１０に設定することができる。前記の（２ａ）式および（２ｂ）式は、ユール・ニールセン分光ノイゲバウアモデル(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とインク量ｄcとの関係を示している。ここでは、１種類のインクのインク量の範囲０〜ｄcmaxも３つの区間に分割されており、区間毎に０から１まで単調に増加する非線形の曲線によってセル分割モデルにて使用する仮想的なインク被覆率ｆcが求められる。他のインクについても同様にインク被覆率ｆm，ｆyが求められる。
【００２９】
図８（Ｃ）は、図８（Ａ）の中央のセルＣ５内にある任意のインク量セット（ｄc，ｄm）にて印刷を行った場合の予測分光反射率Ｒs（λ）の算出方法を示している。インク量セット（ｄc，ｄm）にて印刷を行った場合の分光反射率Ｒ（λ）は、以下の（３）式で与えられる。
【数３】
2010147783.tif 000005
ここで、（３）式におけるインク被覆率ｆc，ｆmは図８（Ｂ）のグラフで与えられる値である。また、セルＣ５を囲む４つの格子点に対応する分光反射率Ｒ（λ）11，（λ）12，（λ）21，（λ）22は分光反射率データベースＲＤＢを参照することにより取得することができる。これにより、（３）式の右辺を構成するすべての値を確定することができ、その計算結果として任意のインク量セットφ（ｄc，ｄm）にて印刷を行った場合の予測分光反射率Ｒs（λ）を算出することができる。波長λを可視波長域にて順次シフトさせていくことにより、可視波長域における予測分光反射率Ｒs（λ）を得ることができる。このように、インク量空間を各格子点によって囲まれた複数のセルに分割すれば、分割しない場合に比べて予測分光反射率Ｒs（λ）をより精度良く算出することができる。
【００３０】
以上のようにして予測対象のインク量セットφに基づいて、プリンタ２０が印刷した場合の予測分光反射率Ｒs（λ）が予測できると、ステップＳ１３０において、ＣＰＭ_Ｐ２ｂが当該予測分光反射率Ｒs（λ）を取得し、当該予測分光反射率Ｒs（λ）に基づいて色彩値を予測する。
【００３１】
図９は、本発明のＣＰＭ_Ｐ２ｂが予測分光反射率Ｒs（λ）に基づいて色彩値を予測する様子を模式的に示している。同図において、ＲＰＭ_Ｐ２ａが予測した予測分光反射率Ｒs（λ）の各波長λにおいて所望の光源の分光エネルギーのスペクトルを乗算することにより、印刷物からの反射光のスペクトルを予測する。さらに、反射光のスペクトルに対して所望の観察条件での感度関数ｘ（λ），ｙ（λ），ｚ（λ）を畳み込み、正規化をすることにより、三刺激値ＸＹＺを算出する。本実施形態においては、特に示さない限りＣＩＥ１９３１  ２°観測者の観察条件で三刺激値ＸＹＺを算出するものとする。本実施形態では光源として、ＣＩＥ標準のＤ６５光源を入力する。さらに、ＣＰＭ_Ｐ２ｂは、三刺激値ＸＹＺにＣＩＥ標準の変換式を適用することにより、ＣＩＥＬＡＢ色空間のＬ*ａ*ｂ*値を算出する。このように、本発明の予測手段としてのＲＰＭ_Ｐ２ａとＣＰＭ_Ｐ２ｂが予測を行うことにより、任意のインク量セットφにて印刷を行った場合の印刷結果として予測分光反射率Ｒs（λ）および色彩値を得ることができる。
【００３２】
Ａ−２．分光反射率データベースの作成
図１０は、以上説明したＲＰＭ_Ｐ２ａが使用する分光反射率データベースＲＤＢの作成処理の流れを示している。図２に示すように、データベース作成部（ＤＢ作成部）Ｐ２ｃは、グラデーションパッチデータ作成部（ＧＰ作成部）Ｐ２ｃ１と、格子点設定部Ｐ２ｃ２と格子点最適化部Ｐ２ｃ３とデータベース作成部（ＤＢ作成部）Ｐ２ｃ４とカラーパッチデータ作成部（ＣＰ作成部）Ｐ２ｃ５とから構成されている。これらの各モジュールが、図１０に示す分光反射率データベースＲＤＢの作成処理を実行する。
【００３３】
ステップＳ２００においては、ＤＢＧ_Ｐ２ｃが分光反射率データベースＲＤＢの作成条件を取得する。ここでは、印刷に使用する機種、メディア（印刷用紙等の記録媒体）、インクセット、解像度、総インク打込量制限（ＴＤＬ）、観察光源等を取得する。以下に説明する本実施形態では、インクセットとしてＣＭＹＫｌｃｌｍが取得され、印刷用紙として光沢紙が取得され、観察光源としてＤ６５光源が取得されているものを例にとって説明する。
【００３４】
ステップＳ２１０においては、ＧＰ作成部Ｐ２ｃ１が、各インクについてインク量を等量ずつ変化させたグラデーションパッチ（以下、ＧＰ１と記載する。）をプリンタ２０に印刷させるための処理を実行する。まず、ＧＰ作成部Ｐ２ｃ１は、ＧＰ１を印刷するためのインク量セットφ１を生成する。
図１１は、ＧＰ１を印刷するためのインク量セットφ１の一部を示している。同図において、ＣインクについてのＧＰ１を印刷するためのインク量セットφ１が示されている。ここでは、Ｃインクを除くＭＹＫｌｃｌｍインク量ｄM，ｄy，ｄk，ｄlc，ｄlmを０とし、Ｃインク量ｄCを０〜２５５まで等量ずつ変化させたインク量セットφ１が生成されている。同様に他のＭＹＫｌｃｌｍインクについても、ＭＹＫｌｃｌｍインク量ｄM，ｄy，ｄk，ｄlc，ｄlmを単独で０〜２５５の間で変化させたインク量セットφ１を生成する。これにより、２５６×６個のグラデーションパッチを印刷するためのインク量セットが生成されることとなる。
すなわち、インク量空間の各インク量軸上のインク量セットφが、ＧＰ１を印刷するためのインク量セットφとして生成される。なお、本実施形態では、Ｃインク量ｄCを０〜２５５の間で変化させる間隔は、１としているが、この間隔は、２以上であってもよい。ただし、後述の等色差グラデーションパッチ（ＧＰ２）を作成する都合上、ＧＰ１は密な間隔で作成されることが望ましい。
【００３５】
図１２は、ＧＰ１の一例を示している。各インクについてグラデーションパッチを印刷するためのインク量セットφ１は、それぞれ矩形状のグラデーションパッチが印刷できるように、プリンタ２０に出力する印刷データにおいてレイアウトされる。図１２の例では、各インクについてインク量ｄC，ｄM，ｄy，ｄk，ｄlc，ｄlmの少インク量側から多インク量側にグラデーションパッチが隙間なく並ぶようにレイアウトされている。ＧＰ１を印刷するための印刷データはプリンタドライバ（ＰＤＶ）Ｐ３に出力され、ＰＤＶ_Ｐ３が順次ハーフトーン処理やラスタライズ処理を行って、プリンタ２０にグラデーションパッチを印刷させる。なお、印刷用紙はステップＳ１００にて指定された光沢紙とする。
【００３６】
ステップＳ２２０においては、ＧＰ作成部Ｐ２ｃ１が分光反射率計ドライバ（ＭＤＶ）Ｐ４を制御して、ステップＳ２１０で印刷されたＧＰ１の分光反射率Ｒ１（λ）の測定を分光反射率計３０に実行させる。ここでは、すべてのＧＰ１の分光反射率Ｒ１（λ）を測定し、その分光反射率Ｒ１（λ）をＭＤＶ_Ｐ４を介してＧＰ作成部Ｐ２ｃ１が取得する。
【００３７】
ステップＳ２３０においては、ＧＰ作成部Ｐ２ｃ１が、各インクについてグラデーションパッチ間の色差が均等になるように変化させた等色差グラデーションパッチ（以下、ＧＰ２と記載する。）をプリンタ２０に印刷させるための処理を実行する。まず、取得した分光反射率Ｒ１（λ）をＣＰＭ Ｐ２ｂに出力し、図９に示す手順に従って各グラデーションパッチの色彩値を算出する。すなわち、光沢紙に印刷されたＧＰ１にＤ６５光源を照射したときの各波長における反射光の分光エネルギーとして、分光反射率Ｒ（λ）と分光エネルギーＰ（λ）を各波長に掛け合わせた値を算出し、反射光の分光エネルギーのスペクトルに対して人間の分光感度特性に応じた等色関数ｘ（λ），ｙ（λ），ｚ（λ）をそれぞれ畳み込み積分し、係数ｋによって正規化することにより、３刺激値Ｘ，Ｙ，Ｚを得る。以上を数式で表すと下記の（４）式となる。
【数４】
2010147783.tif 000006
そして３刺激値Ｘ，Ｙ，Ｚを公知の変換式によって変換することにより、ターゲットＴＧにＤ６５光源を照射したときの色を示すＬ*ａ*ｂ*値を得ることができる。ＧＰ作成部Ｐ２ｃ１は、得られたＬ*ａ*ｂ*値に基づいて、各インクについて隣接するパッチとの色差が等しくなるようにＧＰ２を印刷するためのインク量セットφ２を生成する。すなわち後述の格子点となりうるインク量セットとして、等色差となるインク量セットφ２を用意することにより、効率的に分光反射率データベースに登録すべきインク量セットを探索できるようになる。
【００３８】
ステップＳ２４０においては、ＧＰ作成部Ｐ２ｃ１が、ＧＰ１を印刷するための印刷データはプリンタドライバ（ＰＤＶ）Ｐ３に出力し、ＰＤＶ_Ｐ３が順次ハーフトーン処理やラスタライズ処理を行って、プリンタ２０にグラデーションパッチを印刷させる。なお、印刷用紙はステップＳ１００にて指定された光沢紙とする。
【００３９】
ステップＳ２５０においては、ＧＰ作成部Ｐ２ｃ１が分光反射率計ドライバ（ＭＤＶ）Ｐ４を制御して、ステップＳ２４０で印刷されたＧＰ１の分光反射率Ｒ２（λ）の測定を分光反射率計３０に実行させる。ここでは、すべてのＧＰ１の分光反射率Ｒ２（λ）を測定し、その分光反射率Ｒ２（λ）をＭＤＶ_Ｐ４を介してＧＰ作成部Ｐ２ｃ１が取得する。
【００４０】
ステップＳ２６０においては、ＧＰ作成部Ｐ２ｃ１が、取得した分光反射率Ｒ２（λ）と、ＧＰ２を印刷するためのインク量セットφ２とを対応付けたデータをグラデーション測定データとしてＨＤＤ１４等の記憶媒体に記憶させる。以下、ＧＰ２を測定して得られた分光反射率Ｒ２（λ）を参照分光反射率Ｒr（λ）と表記する。
【００４１】
図１３は、ＧＰ２を印刷するためのインク量セットφ２の一部を示している。同図において、ＣインクについてのＧＰ２を印刷するためのインク量セットφ２が示されている。ここでは、Ｃインクを除くＭＹＫｌｃｌｍインク量ｄM，ｄy，ｄk，ｄlc，ｄlmを０とし、Ｃインク量ｄCを０〜２５５まで等量ずつ変化させたインク量セットφ２が生成されている。同様に他のＭＹＫｌｃｌｍインクについても、ＭＹＫｌｃｌｍインク量ｄM，ｄy，ｄk，ｄlc，ｄlmを単独で０〜２５５の間で変化させたインク量セットφ２を生成する。これにより、２５６×６個のグラデーションパッチを印刷するためのインク量セットが生成されることとなる。ＧＰ２においては、等色差でインク量セットが生成されているので、Ｃインク量ｄCを０〜２５５の間で変化させる間隔は、図１３に示すようにランダムになっている。
【００４２】
ステップＳ２７０においては、ＧＰ作成部Ｐ２ｃ１が、ＧＰ２の色彩値を算出する。すなわち、ＨＤＤ１４からグラデーション測定データを取得し、図９に示した手順で演算を行って、各参照分光反射率Ｒr（λ）の物体に所定光源（本実施形態ではＤ６５光源）を照射したときの色彩値（Ｌ*ａ*ｂ*値）を順次算出する。順次計算された色彩値はグラデーション測定データに参照色彩値として追記される。以上のように用意されたグラデーション測定データを利用して、ステップＳ２８０の格子点最適化処理が実行される。
【００４３】
図１４は、格子点最適化処理のフローチャートである。同図に示す格子点最適化処理が開始されると、ステップＳ４１０において、格子点設定部Ｐ２ｃ２が、インクセットの中から２色の組合せを選択する。例えば、ＣＭＹＫｌｃｌｍの６色のインクセットを有するプリンタであれば、６Ｃ２＝１５通りの組合せの中から２色の組合せを選択する。ステップＳ４１０〜Ｓ４４０は、これら１５通り全ての組合せが選択されるまで繰り返し実行される。なお、以下の説明ではＣとＭの２色が選択された場合を例にとって説明を行う。
【００４４】
ステップＳ４２０においては、格子点設定部Ｐ２ｃ２が、ＣＭインク量平面に格子点を仮設定する。格子点は、上述したグラデーション測定データに登録されたインク量セットの中から選択される。まず、ＣＭインク量平面におけるＴＤＬ内に含まれる格子点数を設定する。この格子点数は、諸条件（機種、メディア、解像度、ＴＤＬ等）に基づいて経験的に決定されるものであり、例えば、予めＨＤＤ１４等に記憶されている数値を取得したり、ユーザからの入力設定を受付けたりしてその数を格子点数として設定する。
【００４５】
図１５は、ＣＭインク量平面に格子点を仮設定した状態を示す図である。同図において、○で図示するポイントが仮設定された格子点を示している。同図においては、ＣＭインク量平面においてインク量の最少量側および最多量側の両端（０，２５５）を含めて各インク量軸に５個の格子点Ｃ０〜Ｃ４，Ｍ０〜Ｍ４が設定されるようになっており、５×５個の直交格子が形成されている。このときＣＭインク量平面においてＴＤＬ内に含まれる格子点数は、格子点の配置位置とＴＤＬとに応じて９〜２４個となる。図１５においては、１２個の格子点がＴＤＬ内に含まれているが、このＴＤＬが異なるインク量平面だったり、Ｃ１〜Ｃ３の位置を変えたりすることにより、ＴＤＬ内に含まれる格子点数が増減することが分かる。つまり、Ｃ１〜Ｃ３，Ｍ１〜Ｍ３の位置を変更することにより、ＴＤＬ内に上記設定された数の格子点が含まれるように格子点の仮設定を実行する。
【００４６】
格子点の仮設定は自動で行ってもよいし手動で行ってもよい。自動で行う場合は、例えば、ｄｍ＝ｄｃに対して線対称に格子点を配置しつつ、格子点位置をずらしながらＴＤＬ内に入る格子点数をカウントし、ＴＤＬ内の格子点が上記設定された数の格子点となるように調節を行うことが考えられる。また、手動で行う場合は、図１５に示すようなＵＩを視認しつつマウスポインタ等を用いて格子点位置をユーザが設定し、ＴＤＬ内の格子点数を調整することが考えられる。なお、格子点は、偏りなく配置されることが好ましく、例えば、Ｃインク量軸では高階調側に格子点が多く配置されているのに、Ｍインク量軸においては低階調側に格子点が多く配置されるといった状況が発生しないようにする。そのためにも、自動的に格子点を仮設定する場合は、上述したように上記インク量平面において対称的に格子点が配置されるようにする。
【００４７】
なお、各格子点は、ＣＭインク量ｄ C，ｄMのみならず、他のＹＫｌｃｌｍインク量ｄy，ｄk，ｄlc，ｄlmについても同様に仮設定される。従って、仮設定の段階で、格子点の数は、５6個となり、各インク量軸上の格子点の数はそれぞれ５個ずつとなる。そして、ＣＥＥ_Ｐ２ｃ４が、当該グラデーション測定データに記述されたインク量セットφ２のうち、仮設定した格子点と一致するものに格子点フラグをセットする。上述したように、グラデーション測定データに記述されたインク量セットφ２はインク量空間の各軸上のものであり、格子点は各インク量軸上に５個ずつ仮設定されているため、各インクのグラデーションについて５個ずつインク量セットφが格子点と一致することとなる。
【００４８】
ステップＳ４３０においては、格子点最適化部Ｐ２ｃ３が、仮設定された格子点に対して制約条件を設定する。この制約条件は、ＴＤＬに隣接する格子点に対応するインク量セットのインク量合計とＴＤＬとの関係を定める制約条件と、各インク量軸における格子点の並び順および各インク量軸においてインク量が最大の格子点と最小の格子点位置を固定する制約条件である。まず、ＴＤＬとの関係において定められる制約条件は、図１５に示す例においては、下記（５）式のようになる。
【数５】
2010147783.tif 000007
【００４９】
すなわち、上記ステップＳ４２０において仮設定された段階でＴＤＬよりもインク量が少ない格子点については、格子点最適化後においてもＴＤＬよりもインク量が少なくなるようにし、仮設定された段階でＴＤＬよりもインク量が多い格子点については、格子点最適化後においてもその関係を維持させるようにＣ１〜Ｃ４，Ｍ１〜Ｍ４に制約を課す。ＴＤＬは、プリンタ２０や印刷用紙やインクの物性等によって定められており、当該制限値を超えるような合計値では、例えばインクにじみなどの不具合が生じ、正常な印刷が不可能となるからである。このようにＴＤＬ近傍の格子点に着目して、ＴＤＬとの大小関係を不等式で表した制約条件が一つ目の制約条件である。
【００５０】
なお、この制約条件にもとづいて、インク量軸上に仮設定された格子点が、適切であるか否かを判断するための条件として以下のものが考えられる。図１６に不適切に仮設定された格子点の例を示した。同図においては、Ｃ４＋Ｍ２がＴＤＬの内側に位置するのに対し、Ｃ２＋Ｍ４がＴＤＬの外側に位置している。このようにＣインク量軸とＭインク量軸とで格子点における添字を転置した関係にある格子点のインク量を足し合わせた値が、ＴＤＬの内外に分かれてしまうように格子点を仮設定すると、格子点の配置が偏る可能性がある。従って、各々ｍ個の格子点をインク量軸上に設定された２つのインク量軸を含んで形成されるインク量平面において、一方の軸においてｎ番目の格子点と他方の軸においてｍ−ｎ番目の格子点とのインク量を足し合わせた総インク打込量と、一方の軸においてｍ−ｎ番目の格子点と他方の軸においてｎ番目の格子点とのインク量を足し合わせた総インク打込量と、がＴＤＬの内外に分かれないように仮設定を行うことが好ましい。
【００５１】
また、各インク量軸における格子点の並び順及び各インク量軸においてインク量が最大の格子点と最小の格子点位置を固定する制約条件は下記（６）式のようになる。
【数６】
2010147783.tif 000008
すなわち、各一次色のインク量軸上において仮設定された時点での格子点の並び順が格子点最適化後においても維持されるように各格子点の大小関係を規定し、さらに最小の格子点と最大の格子点についてはその位置が変更されないように固定する。
【００５２】
ステップＳ４４０においては、格子点最適化部Ｐ２ｃ３が、全ての２色の組合せについて、格子点数の設定と制約条件の設定とが終了したか否かを判断する。すなわち、上述した６Ｃ２＝１５通りの２色の組合せの全てについて上述した格子点数の設定、格子点の仮設定、制約条件の設定が終了するまでステップＳ４１０〜Ｓ４４０を繰り返す。全ての組合せについて処理が終了するとステップＳ４６０に進む。
【００５３】
ステップＳ４６０においては、格子点最適化部Ｐ２ｃ３が、インク量セットφ２の中から、制約条件を満たす格子点を選択して格子点とする。すなわち、ステップＳ４２０にてＣ１〜Ｃ３，Ｍ１〜Ｍ３，Ｙ１〜Ｙ３，Ｋ１〜Ｋ３，ｌｃ１〜ｌｃ３，ｌｍ１〜ｌｍ３に対応して付けた格子点フラグを解除し、Ｃ１〜Ｃ３，Ｍ１〜Ｍ３，Ｙ１〜Ｙ３，Ｋ１〜Ｋ３，ｌｃ１〜ｌｃ３，ｌｍ１〜ｌｍ３のそれぞれに対してステップＳ４３０において設定された制約条件を満たすようなインク量セットφ２の組合せを選択し、選択されたインク量セットφ２に格子点フラグを付ける。なお、Ｃ０，Ｃ４，Ｍ０，Ｍ４，Ｙ０，Ｙ４，Ｋ０，Ｋ４，ｌｃ０，ｌｃ４，ｌｍ０，ｌｍ４等のように各インク量軸において端に位置する格子点については、格子点フラグの変更は行わない。
【００５４】
具体的には、上記等色差グラデーションパッチに設定された階調値の組合せから階調値０と２５５とを除いた全ての階調値から、ＣＭＹＫｌｃｌｍ各色から３つの階調値をそれぞれ選択し、上記（３）式を満たすように上記Ｃ１〜Ｃ３，Ｍ１〜Ｍ３，Ｙ１〜Ｙ３，Ｋ１〜Ｋ３，ｌｃ１〜ｌｃ３，ｌｍ１〜ｌｍ３に設定し、設定された組み合わせが上記（５）、（６）式を満たすか否かを判断し、両式を満たすインク量セットに格子点フラグを付す。
【００５５】
ステップＳ４７０においては、格子点最適化部Ｐ２ｃ３が、目的関数ΔＥを算出する。すなわち、設定された格子点に基づいて上記分光プリンティングモデルにてインク量セットφ２の分光反射率を予測し、予測された分光反射率から求められる色彩値と、各インク量セットφ２の実際の色彩値との差分の可算和を求める。
【００５６】
図１７は、格子点フラグがセットされたグラデーション測定データの一部を示している。同図において、グラデーション測定データに記述されたインク量セットφ２のうち、格子点と一致するもののみに格子点フラグがセットされている。ステップＳ４６０においては、格子点最適化部Ｐ２ｃ３が、上述した分光反射率予測部（ＲＰＭ）Ｐ２ａに対して、グラデーション測定データに記述されたインク量セットφ２のうち、格子点フラグがセットされていないインク量セットφ２を順次出力する。それと同時に、格子点最適化部Ｐ２ｃ３は、ＲＰＭ_Ｐ２ａに対して、グラデーション測定データのうち格子点フラグがセットされたものを分光反射率データベースＲＤＢとして使用するように指示する。
【００５７】
これを受けて、ＲＰＭ_Ｐ２ａは、入力された各インク量セットφに対して上述した分光プリンティングモデルを適用して、予測分光反射率Ｒs（λ）を算出する。グラデーション測定データのうち格子点フラグがセットされたものを分光反射率データベースＲＤＢとして使用するため、ステップＳ２７０では当該格子点フラグがセットされたインク量セットφ２が分光プリンティングモデルにおける格子点（各セルの頂点）として使用されることとなる。当該格子点は、インク量軸上のものに限定されるが、予測対象のインク量セットφ２もインク量軸上のものに限定されるため予測を行うことができる。すなわち、上記の（３）式において、Ｃインク量軸上のインク量セットφ２および格子点を考えた場合、Ｍインクのインク被覆率ｆm＝０となり、上記の（３）式は下記の（７）式に書き換えられる。
【数７】
2010147783.tif 000009
【００５８】
上記（７）式では、インク量軸上の格子点の分光反射率Ｒ11（λ），Ｒ21（λ）のみに基づいて、当該格子点の間のインク量軸上の任意のインク量セットφ２についての予測分光反射率Ｒs（λ）を算出することができる。なお、図８（Ｃ）に図示したセルはインク量軸上の線分になると考えることができ、矩形状のセルの頂点を示す格子点は線分の端点になると考えることができる。むろん、Ｃインク以外のインクについても、インク量軸上の格子点の分光反射率Ｒ（λ）に基づいて、当該インク量軸上の任意のインク量セットφ２についての予測分光反射率Ｒs（λ）を算出することができる。以上のようにしてインク量軸上の各インク量セットφ２について予測分光反射率Ｒs（λ）が算出できると、当該予測分光反射率Ｒs（λ）をグラデーション測定データに追記する。
【００５９】
そして格子点最適化部Ｐ２ｃ３が、グラデーション測定データに追記した予測分光反射率Ｒs（λ）を色予測部（ＣＰＭ）Ｐ２ｂに順次出力する。これにより、ＣＰＭ_Ｐ２ｂは、図９に示した手順を行って、各予測分光反射率Ｒs（λ）に所定光源（ステップＳ１００で指定されたＤ６５光源）を照射したときの色彩値（Ｌ*ａ*ｂ*値）を順次算出させる。順次計算された予測色彩値は、ＣＥＥ_Ｐ２ｃ４によってグラデーション測定データに予測色彩値として追記される。
【００６０】
次に格子点最適化部Ｐ２ｃ３は、上述した参照色彩値と予測色彩値との色差をグラデーション測定データの各インク量セットφ２について算出する。色差は、例えばＣＩＥ２０００色差式によって算出される。算出された色差は、グラデーション測定データに追記される。なお、格子点フラグがセットされた格子点については、色差を予め０とする。上記の（４）式に格子点のインク量セットφ２を代入しても、格子点の分光反射率Ｒ（λ）がそのまま予測分光反射率Ｒs（λ）として算出され、色差は０となるからである。そして、各インク量セットφ２について算出された色差の和を目的関数ΔＥとして設定する。すなわち目的関数ΔＥは、設定された各格子点についての実測色彩値と予測色彩値の差分の和で表される。各格子点のΔＥの総和は下記（８）式で表される。
【数８】
2010147783.tif 000010
ここで算出される目的関数ΔＥが、本発明の誤差に相当する。
【００６１】
ステップＳ４７５においては、ステップＳ４７０において算出されたΔＥが、これまでで算出された中で最小のΔＥであるか否か、すなわちｍｉｎΔＥであるか否かを判断する。ＲＡＭ１２にはこれまでに算出された中で最も小さいｍｉｎΔＥとそのｍｉｎΔＥに対応するインク量セットの組合せを示す情報とが記憶されている。従って、ＲＡＭ１２に記憶されたｍｉｎΔＥとステップＳ４７０で算出されたΔＥとを比較し、算出されたΔＥがｍｉｎΔＥよりも小さい場合はこの算出されたΔＥを新たなｍｉｎΔＥとしてＲＡＭ１２に記憶してステップＳ４８０に進む。その際、ΔＥの算出元となったインク量セットの組合せを示す情報も更新される。算出されたΔＥがｍｉｎΔＥよりも大きい場合は、ＲＡＭ１２のｍｉｎΔＥとこれに対応するインク量セットを更新せずにステップＳ４８０に進む。
【００６２】
ステップＳ４８０においては、上記制約条件の範囲で格子点として設定可能な全てのインク量セットについて、Ｃ１〜Ｃ３，Ｍ１〜Ｍ３，Ｙ１〜Ｙ３，Ｋ１〜Ｋ３，ｌｃ１〜ｌｃ３，ｌｍ１〜ｌｍ３へ設定完了したか否かを判断する。設定完了していれば、ステップＳ４９０に進み、設定完了していなければステップＳ４６０に戻って新たなインク量セットに対して格子点フラグを設定する。
【００６３】
ステップＳ４９０においては、ＤＢ作成部Ｄ２ｃ４が、ＲＡＭ１２に記憶されたインク量セットの組合せについて、上記グラデーション測定データに格子点フラグをつけて、図１０のステップＳ２９０に処理に進む。
【００６４】
ステップＳ２９０においては、ＣＰデータ作成部Ｐ２ｃ５が、格子点フラグを付されたインク量軸上の格子点に基づいて、カラーパッチとして印刷すべきインク量セットを設定する。そのために、まず、格子点フラグを付されたインク量軸上の格子点を組み合わせることにより、インク量空間全体に直交格子上の格子点を生成する。さらに、生成された格子点のうち、各インク量ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlmの合計値（インクデューティ値）が所定の制限値内にあるものとこの所定の制限値から所定量だけ大きい範囲に入っているものとを除いて、実際に使用する格子点から除外し、除外されなかった格子点をカラーパッチ印刷用の格子点として設定する。
【００６５】
図１８は、ＣＭインク量平面およびｌｃｌｍインク量平面にて生成された格子点を示している。同図には、ＣＭインク量平面における制限値を図示しており、当該制限値よりインク量ｄc，ｄm（ｄy，ｄk，ｄlc，ｄlmがある一定値のとき）の合計が小さい格子点（○で図示）と当該制限値のすぐ外にある格子点（★で図示）とを除いて、カラーパッチの対象外の格子点（●で図示）で示されている。このような制限値は、プリンタ２０や印刷用紙やインクの物性等によって定められており、当該制限値を超えるような合計値では、例えばインクにじみなどの不具合が生じ、正常な印刷が不可能となる。ただし、制限値のすぐ外にある格子点については制限値付近の予測精度向上のためにカラーパッチ印刷用の格子点に採用する。図示の便宜上、ＣＭインク量平面における制限値の判定を例示したが実際には、ＣＭＹＫｌｃｌｍインク量ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlmの６次元空間で制限値の判定を行う。
【００６６】
ステップＳ３００においては、ＣＰデータ作成部Ｐ２ｃ５が、ステップＳ２９０にて除外されなかった格子点についてのカラーパッチをプリンタ２０に印刷させるための処理を実行する。ここでは、少なくとも各格子点が示すインク量セットφを、カラーパッチを印刷するためのインク量セットφとする。
【００６７】
ステップＳ３１０においては、ＣＰデータ作成部Ｐ２ｃ５と分光反射率計ドライバ（ＭＤＶ）Ｐ４がカラーパッチの分光反射率Ｒ（λ）の測定を分光反射率計３０に実行させる。ここでは、すべてのカラーパッチの分光反射率Ｒ（λ）を測定し、その分光反射率Ｒ（λ）をＭＤＶ_Ｐ４を介してデータベース記憶部（ＤＢＲ）Ｐ２ｃ９が取得する。
【００６８】
そして、ステップＳ３２０では、測定した分光反射率Ｒ（λ）と、グラデーションパッチを印刷するためのインク量セットφとを対応付けたデータを生成し、分光反射率データベースＲＤＢとしてＨＤＤ１４に記憶させる。記憶した分光反射率データベースＲＤＢは、Ａ−１．節で説明した分光反射率と色の予測において使用されることとなる。
【００６９】
Ｂ．色変換ルックアップテーブルの作成
図１９は、色変換ルックアップテーブル作成処理の流れを示している。図２に示すように、ＡＰＬ_Ｐ２のプロファイル作成部（ＬＴＭ）Ｐ２ｄは、ＬＵＴ用格子点生成部（ＬＧＧ）Ｐ２ｄ１とＬＵＴ生成部（ＬＴＧ）Ｐ２ｄ２とから構成されている。ステップＳ５００においては、ＬＧＧ_Ｐ２ｄ１がインク量空間においてＬＵＴ用格子点を生成する。ＬＵＴ用格子点を生成するには、プリンタ２０の色再現ガマットの全体を網羅するようにＬＵＴ用格子点を生成させるとともに、例えば粒状性や色恒常性や平滑性やトータルインクデューティが総合的に良好となる空間域からＬＵＴ用格子点を選択するのが望ましい。例えば、国際公開ＷＯ２００５／０４３８８４号のパンフレットに開示された手法によって好ましいＬＵＴ用格子点を選択することができる。
【００７０】
以上のようにしてＬＵＴ用格子点が用意できると、ステップＳ５１０において、ＬＴＧ_Ｐ２ｄ２は、ＬＵＴ用格子点が示すインク量セットφを分光反射率予測部（ＲＰＭ）Ｐ２ａに対して順次出力し、ＬＵＴ用格子点が示すインク量セットφで印刷した場合の予測分光反射率Ｒs（λ）を算出させる。このとき、上述した手法によって作成された分光反射率データベースＲＤＢが使用されることとなる。さらに、ステップＳ５２０においては、予測分光反射率Ｒs（λ）を色予測部（ＣＰＭ）Ｐ２ｂに順次出力する。これにより、ＣＰＭ_Ｐ２ｂは、図９に示した手順を行って、各予測分光反射率Ｒs（λ）に所定光源（本実施形態ではＤ６５光源）を照射したときの色彩値（Ｌ*ａ*ｂ*値）を順次算出させる。ステップＳ５３０では、ＬＴＧ_Ｐ２ｄ２が、順次計算された色彩値とＬＵＴ用格子点が示すインク量セットφとの対応関係を規定した色変換ルックアップテーブルを作成する。この色変換ルックアップテーブルによれば、絶対色空間としてのＣＩＥＬＡＢ色空間と、プリンタ２０固有のインク量空間との対応関係を得ることができる。
【００７１】
Ｃ．印刷制御装置（処理）
図２０は、印刷制御処理の流れを示している。図２に示すように、プリンタドライバ（ＰＤＶ）Ｐ３は、レンダラ（ＲＮＰ）Ｐ３ａと色変換部（ＣＣＰ）Ｐ３ｂとハーフトーン部（ＨＴＰ）Ｐ３ｃとラスタライザ（ＲＬＰ）Ｐ３ｄとから構成されている。ステップＳ６００においては、ＲＮＰ_Ｐ３ａが例えば描画コマンド等から構成された印刷指示データに応じて、ｓＲＧＢ色空間のＲＧＢ値で各画素の色が表現された印刷イメージを生成する。ステップＳ６１０においては、ＣＣＰ_Ｐ３ｂが当該印刷イメージに対して色変換処理を実行する。ここでは、まずｓＲＧＢ色空間とＣＩＥＬＡＢ色空間との対応関係を規定したｓＲＧＢプロファイルを使用して、印刷イメージを各画素がＣＩＥＬＡＢ色空間で表現されたデータに変換する。さらに、各画素がＣＩＥＬＡＢ色空間で表現された印刷イメージを、上述した色変換ルックアップテーブルを参照して、インク量空間のデータに変換する。その際に、線形補間等の補間演算を実行する。また、ＣＩＥＬＡＢ色空間において、ガマットマッピングを行う。ステップＳ６２０では、ＨＴＰ_Ｐ３ｃがインク量空間で各画素が表された印刷イメージを取得し、誤差拡散法やディザ法によってハーフトーンデータに変換する。さらに、ステップＳ６３０では、ＲＬＰ_Ｐ３ｄがハーフトーンデータを取得し、当該ハーフトーンデータをプリンタ２０における複数のノズル２１ａ，２１ａ・・・への吐出可否を指示するためのデータに変換し、プリンタ２０に出力する。
【００７２】
Ｄ．変形例
ところで、上述した（５）式や（６）式の制約条件に代えて、下記（９）式の制約条件を課すことによっても前記実施形態と同様の効果を得ることができる。無論、前記（９）式の制約条件を、前記（５）式や前記（６）式の制約と併用しても構わない。
【数９】
2010147783.tif 000011
上記（９）式において、Ｔ１〜Ｔ４は定数であり、ステップ４２０でインク量軸上に格子点が仮設定された時点の各格子点に対応するインク量セットに基づいて決定される。すなわち、各インク量軸上に同数の格子点を配置し、各軸上において原点からの順番が等しい格子点のインク量の和が、それぞれＴ１〜Ｔ４となるように制限する。
【００７３】
図２１は、前記（９）式で表される制約条件を説明する図である。同図は、説明の便宜上、ＣＭインク量平面を例にとって示してある。例えば、各色２５６階調を取りうる場合に、格子点の仮設定時にＣ０＝０，Ｍ０＝０に設定すると、Ｃ０とＭ０の和であるＴ０は０であるので、制約条件を満たすインク量セットは、Ｃ０＝０、Ｍ０＝０に固定される。また、格子点の仮設定時にＣ４＝２５５，Ｍ４＝２５５に設定すると、Ｃ４とＭ４の和であるＴ４は５１０であるので、制約条件を満たすインク量セットは、Ｃ４＝２５５、Ｍ４＝２５５に固定される。また、前記（９）式の制約を課すことにより、例えばＣ１＋Ｍ１で規定される格子点はその和がＴ１となるような値（図２１においてＴ１の二点鎖線上の値）しか取れなくなる。このような制約下において、いずれかのインク量軸において高階調側に格子点を偏在させつつ他のインク量軸において低階調側に格子点を偏在させるような非対称的な配置を取ると、前記目的関数ΔＥが上昇してしまうことを考えるとＣ１＋Ｍ１やＣ２＋Ｍ２やＣ３＋Ｍ３がｄｃ＝ｄｍの近辺に配置されることになる。すなわち、前記（９）式によって格子点に課す制約は、前記（６）式と略同等の作用を奏すると言える。
【００７４】
さらに、図２１に示すように、ＴＤＬは傾きが「−１」に近く、前記（９）式の制約下においてＣ１＋Ｍ１やＣ２＋Ｍ２やＣ３＋Ｍ３が変動可能な方向とＴＤＬの方向とは略一致している。従って、仮設定された時点でＴＤＬよりも総インク打込量が多い格子点は、前記（９）式の制約条件の下で目的関数ΔＥを極小化した時点でもＴＤＬの外側に位置し、仮設定された時点でＴＤＬよりも総インク打込量が少ない格子点は目的関数ΔＥが極小化した時点でもＴＤＬの内側に位置する可能性が高い。すなわち、前記（９）式によって格子点に課す制約は、前記（５）式と略同等の作用を奏すると言える。
【図面の簡単な説明】
【００７５】
【図１】印刷結果予測装置のハードウェア構成を示すブロック図である。
【図２】印刷結果予測装置のソフトウェア構成を示すブロック図である。
【図３】印刷結果予測処理のフローチャートである。
【図４】分光反射率データベースの一例を示す図である。
【図５】ＣＭインク量平面における格子点を示す図である。
【図６】プリンタによる印刷を模式的に説明する図である。
【図７】分光ノイゲバウアモデルを示す図である。
【図８】セル分割ユール・ニールセン分光ノイゲバウアモデルを示す図である。
【図９】分光反射率から色を特定する様子を示す図である。
【図１０】分光反射率データベース作成処理のフローチャートである。
【図１１】グラデーションパッチφ１を印刷するためのインク量セットを示す図である。
【図１２】グラデーションパッチφ１を示す図である。
【図１３】グラデーションパッチφ２を印刷するためのインク量セットを示す図である。
【図１４】格子点最適化処理のフローチャートである。
【図１５】ＣＭインク量平面に格子点を仮設定した状態を示す図である。
【図１６】仮設定された格子点として不適切な例を示す図である。
【図１７】グラデーション測定データを示す図である。
【図１８】ＣＭインク量平面において生成された格子点を示す図である。
【図１９】色変換ルックアップテーブル作成処理のフローチャートである。
【図２０】印刷制御処理のフローチャートである。
【図２１】変形例にかかる制約条件を説明する図である。
【符号の説明】
【００７６】
１０…コンピュータ、１１…ＣＰＵ、１２…ＲＡＭ、１３…ＲＯＭ、１４…ＨＤＤ、１５…ＧＩＦ、１６…ＶＩＦ、１７…ＩＩＦ、１８…バス、Ｐ１…ＯＳ、Ｐ２…ＡＰＬ、Ｐ２ａ…ＲＰＭ、Ｐ２ｂ…ＣＰＭ、Ｐ２ｃ…データベース作成部、Ｐ２ｃ１…ＧＰ作成部、Ｐ２ｃ２…格子点設定部、Ｐ２ｃ３…格子点最適化部、Ｐ２ｃ４…ＤＢ作成部、Ｐ２ｃ５…ＣＰ作成部、Ｐ２ｄ…ＰＦＧ、Ｐ２ｄ１…ＬＧＧ、Ｐ２ｄ２…ＬＴＧ、Ｐ３…ＰＤＶ、Ｐ３ａ…ＲＮＰ、Ｐ３ｂ…ＣＣＰ、Ｐ３ｃ…ＨＴＰ、Ｐ３ｄ…ＲＬＰ、Ｐ４…ＭＤＶ、Ｐ５…ＤＤＶ 
【図１】
2010147783.tif 000012
【図２】
2010147783.tif 000013
【図３】
2010147783.tif 000014
【図４】
2010147783.tif 000015
【図５】
2010147783.tif 000016
【図６】
2010147783.tif 000017
【図７】
2010147783.tif 000018
【図８】
2010147783.tif 000019
【図９】
2010147783.tif 000020
【図１０】
2010147783.tif 000021
【図１１】
2010147783.tif 000022
【図１３】
2010147783.tif 000023
【図１４】
2010147783.tif 000024
【図１５】
2010147783.tif 000025
【図１６】
2010147783.tif 000026
【図１７】
2010147783.tif 000027
【図１８】
2010147783.tif 000028
【図１９】
2010147783.tif 000029
【図２０】
2010147783.tif 000030
【図２１】
2010147783.tif 000031
【図１２】
2010147783.tif 000032