JP2011077867A

JP2011077867A - プロファイル作成方法、プログラム、印刷装置、および、表示制御装置

Info

Publication number: JP2011077867A
Application number: JP2009227787A
Authority: JP
Inventors: Nao Kaneko; 奈緒金子; Hisanori Nakajima; 久典中島
Original assignee: Seiko Epson Corp
Current assignee: Seiko Epson Corp
Priority date: 2009-09-30
Filing date: 2009-09-30
Publication date: 2011-04-14

Abstract

【課題】目的関数等を複数準備することなく複数の印刷媒体についてのプロファイルを作成する技術を提供する。
【解決手段】機器独立表色系の複数の格子点の分布を平滑化する平滑化処理と、平滑化された前記複数の格子点が示す色彩値を基本印刷媒体において再現させるインク量を、基本印刷媒体にインクを付着したときの画質を評価する目的関数を用いた最適化によって決定する最適化処理とを実行し、前記最適化処理における前記インク量の範囲を、前記基本印刷媒体に付着可能な前記インク量の制限値と、転用印刷媒体に付着可能な前記インク量の制限値との比である制限値比に基づいて制限し、前記平滑化処理と前記最適化処理で決定された前記インク量を前記制限値比に基づいて変換した前記インク量と、該インク量に対応する格子点の色彩値との対応関係を規定することにより前記転用印刷媒体のためのプロファイルを作成する。
【選択図】図２０

Description

この発明は、プロファイル作成方法、プログラム、印刷装置、および、表示制御装置に関する。

色変換プロファイルは、入力表色系と出力表色系との対応関係を示す情報であり、例えば、色変換ルックアップテーブルや色変換関数などの形式で利用される。色変換ルックアップテーブルの入力表色系の座標値は、入力表色系の色空間内の点の位置を表しており、出力表色系の座標値は出力表色系の色空間内の点の位置を表している。なお、本明細書において、任意の色空間内の点を「色点」または「格子点」とも呼ぶ。また、色変換ルックアップテーブルに登録されている入力値で表される色点および出力値で表される色点を、それぞれ「入力格子点」および「出力格子点」とも呼ぶ。
色変換ルックアップテーブルの入力格子点や出力格子点の配置を平滑なものにする技術として、例えば本願の出願人により開示された特許文献１に記載されたものがある。この平滑化では、Ｌａｂ表色系の格子点を移動させた後、目的関数を用いた最適化処理を利用して、移動後のＬ^*ａ^*ｂ^*格子点に対応する最適なインク量を決定している。この最適なインク量は、目的関数を最小とするようなインク量として決定される。

特開２００６−１９７０８０号公報

しかしながら、印刷媒体の種類ごとに同一のインク量に対する発色特性が異なるとともに、単位面積あたりに付着可能なインク量（デューティ制限値）も異なるため、複数の印刷媒体について色変換ルックアップテーブルを作成するためには印刷媒体ごとに目的関数を準備しておく必要があるという問題があった。さらに、目的関数には、各インク量のインクを印刷媒体に付着した場合に再現される色の色彩値に基づいて画質を評価する項が含まれるため、該色彩値を予測するための分光プリンティングモデルも印刷媒体ごとに準備しておかなければならないという問題があった。

本発明は、印刷媒体ごとに目的関数や分光プリンティングモデルを準備することなく複数の印刷媒体についてのプロファイルを作成する技術を提供することを目的とする。

本発明は、上述の課題の少なくとも一部を解決するためになされたものであり、以下の形態または適用例として実現することが可能である。

［適用例１］
プロファイル作成方法であって、
機器独立表色系の複数の格子点の分布を平滑化する平滑化処理と、平滑化された前記複数の格子点が示す色彩値を基本印刷媒体において再現させるインク量を、基本印刷媒体にインクを付着したときの画質を評価する目的関数を用いた最適化によって決定する最適化処理とを実行し、
前記最適化処理における前記インク量の範囲を、前記基本印刷媒体に付着可能な前記インク量の制限値と、前記基本印刷媒体とは異なる転用印刷媒体に付着可能な前記インク量の制限値との比である制限値比に基づいて制限し、
前記平滑化処理と前記最適化処理で決定された前記インク量を前記制限値比に基づいて変換した前記インク量と、該インク量に対応する格子点の色彩値との対応関係を規定することにより前記転用印刷媒体のためのプロファイルを作成することを特徴とするプロファイル作成方法。
この構成によれば、前記基本印刷媒体にインクを付着したときの画質を評価する目的関数を用いた最適化により、該基本印刷媒体とは異なる前記転用印刷媒体のためのプロファイルを作成することができる。前記最適化処理における前記インク量の範囲が、前記基本印刷媒体に付着可能な前記インク量の制限値と、前記転用印刷媒体に付着可能な前記インク量の制限値と制限値比に基づいて制限されるため、前記基本印刷媒体にインクを付着したときの画質を評価する目的関数が評価可能な前記インク量の範囲で最適化を行わせることができる。さらに、前記最適化処理の後の前記インク量の範囲を前記制限値比に基づいて復元することができる。

［適用例２］
前記制限値比は、前記基本印刷媒体に付着可能な前記インク量の制限値を、前記転用印刷媒体に付着可能な前記インク量の制限値で除算したものであり、
前記最適化処理における前記インク量の範囲は、前記転用印刷媒体に付着可能な前記インク量の制限値に前記制限値比を乗算した仮制限値によって制約されることを特徴とするプロファイル作成方法。
この構成によれば、前記基本印刷媒体に付着可能な前記インク量の制限値によって制限される前記インク量の範囲で最適化を行わせることができる。

［適用例３］
前記平滑化処理と前記最適化処理で決定された前記インク量を、前記制限値比によって除算した前記インク量と、該インク量に対応する前記機器独立表色系の色彩値との対応関係を規定することにより前記プロファイルを作成することを特徴とするプロファイル作成方法。
この構成によれば、最適化の後には、前記基本印刷媒体に付着可能な前記インク量の制限値によって制限される前記インク量の範囲を、前記転用印刷媒体に付着可能な前記インク量の制限値によって制限される前記インク量の範囲に復元することができる。

［適用例４］
前記基本印刷媒体に付着可能な前記インク量の制限値と、前記転用印刷媒体に付着可能な前記インク量の制限値は、前記基本印刷媒体と前記転用印刷媒体においてインクによる発色の変動が一定の基準よりも小さくなる前記インク量であることを特徴とするプロファイル作成方法。
この構成によれば、前記インク量の制限値を特定することができる。

［適用例５］
前記基本印刷媒体に付着可能な前記インク量の制限値と、前記転用印刷媒体に付着可能な前記インク量の制限値は、前記基本印刷媒体と前記転用印刷媒体においてインクによるにじみが観察される前記インク量であることを特徴とするプロファイル作成方法。
この構成によれば、前記インク量の制限値を特定することができる。

なお、本発明は、種々の形態で実現することが可能であり、例えば、プロファイルの作成方法および装置、そのための平滑化／最適化処理方法および装置、プロファイルを印刷装置に組み込む印刷装置の製造方法および製造システム、それらの方法または装置の機能を実現するためのコンピュータープログラム、そのコンピュータープログラムを記録した記録媒体等の形態で実現することができる。さらに、色変換プロファイルが組み込まれた印刷装置においても本発明の有用性を実現することができる。

本発明の一実施例におけるルックアップテーブル作成装置の構成を示すブロック図である。実施例の全体処理手順を示すフローチャートである。重み指定用ＵＩ画像を示す図である。メディアテーブルＭＴＢの一例を示す図である。設定テーブルＳＴＢの一例を示す図である。実施例のベースＬＵＴ作成処理手順を示すフローチャートである。図６のステップＳ１００〜Ｓ３００によってベース３Ｄ−ＬＵＴを作成する場合の処理内容を示す説明図である。入力表色系であるＲＧＢ表色系の色彩値とＬａｂ表色系の色彩値との対応関係を示す説明図である。図６のステップＳ１００〜Ｓ３００によってベース４Ｄ−ＬＵＴを作成する場合の処理内容を示す説明図である。ベースＬＵＴを用いた色補正ＬＵＴの作成方法を示す説明図である。実施例のスムージング処理に利用される力学モデルを示す説明図である。グレー軸格子点に対応する格子点がグレーターゲットに拘束される様子を示す図である。スムージング処理の典型的な処理手順を示すフローチャートである。図１３のステップＴ１００の詳細手順を示すフローチャートである。図１３のステップＳ１２０〜Ｓ１５０の処理内容を示す説明図である。最適化処理（図８のステップＴ１３０）の詳細手順を示すフローチャートである。本発明の一実施例におけるプリンターの構成を示すブロック図である。プリンターの印刷方式を説明する模式図である。プリンターのソフトウェア構成を示すブロック図である。メディア特性指定ＵＩ画像を示す図である。変換関数を設定する様子を示すグラフである。転用メディア色味と被転用メディア色味をａ^*ｂ^*平面にプロットしたグラフである。転用メディアのＬＵＴを作成する場合のグレーターゲットを示す図である。変形例のメディア特性指定ＵＩ画像を示す図である。

次に、本発明の実施の形態を以下の順序で説明する。
Ａ．装置構成と全体処理手順：
Ｂ．基本メディアのＬＵＴ作成手順：
Ｂ−１．全体手順：
Ｂ−２．力学モデル：
Ｂ−３．スムージング処理（平滑化／最適化処理）の処理手順：
Ｂ−４．最適化処理の内容：
Ｃ．印刷装置の構成：
Ｄ．転用メディアのＬＵＴ作成手順：
Ｅ．変形例：

Ａ．装置構成と全体処理手順：
図１は、本発明の一実施例におけるルックアップテーブル作成装置の構成を示すブロック図である。この装置は、ベースＬＵＴ作成モジュール１００と、色補正ＬＵＴ作成モジュール２００と、フォワードモデルコンバーター３００と、ＬＵＴ格納部４００と、ＬＵＴ作成条件設定モジュール７００と、を備えている。「ＬＵＴ」は、色変換プロファイルとしてのルックアップテーブルの略語である。これらのモジュール１００，２００やコンバーター３００の機能は、メモリに格納されたコンピュータープログラムをコンピューターが実行することによってそれぞれ実現される。また、ＬＵＴ格納部４００は、コンピューターが具備するハードディスク装置などの記録媒体によって実現される。また、コンピューターには表示装置（例えば、液晶ディスプレイ。）が接続されており、各処理に必要なＵＩ表示が該表示装置において行われる。さらに、コンピューターには入力装置（例えば、キーボード、マウス。）が接続されており、各処理に必要な情報が該入力装置を介して入力される。また、コンピューターにはプリンター１０と測色機が接続されている。

ベースＬＵＴ作成モジュール１００は、スムージング処理初期値設定モジュール１２０と、テーブル作成モジュール１４０とを有している。スムージング処理モジュール１３０は、色点移動モジュール１３２と、インク量最適化モジュール１３４と、画質評価指数コンバーター１３６とを有している。フォワードモデルコンバーター３００は、分光プリンティングモデルコンバーター３１０と、色算出部３２０とを備えている。ＬＵＴ作成条件設定モジュール７００は、インク量コンバーター７１０と、ＵＩモジュール７２０と、設定情報格納モジュール７３０とを備えている。設定情報格納モジュール７３０は、ハードディスク装置等に格納されたメディアテーブルＭＴＢと設定テーブルＳＴＢとを管理する。これらの各部の機能については後述する。

ＬＵＴ格納部４００は、インバースモデル初期ＬＵＴ４１０や、ベース３Ｄ−ＬＵＴ５１０，ベース４Ｄ−ＬＵＴ５２０，色補正３Ｄ−ＬＵＴ６１０，色補正４Ｄ−ＬＵＴ６２０などを格納するためのものである。ただし、インバースモデル初期ＬＵＴ４１０以外のＬＵＴは、ベースＬＵＴ作成モジュール１００や色補正ＬＵＴ作成モジュール２００によって作成されるものである。ベース３Ｄ−ＬＵＴ５１０は、ＲＧＢ表色系を入力とし、インク量を出力とする色変換ルックアップテーブルである。一方、ベース４Ｄ−ＬＵＴ５２０は、ＣＭＹＫ表色系を入力とし、インク量を出力とする色変換ルックアップテーブルである。なお、「３Ｄ」や「４Ｄ」は、入力値の数を意味している。これらのベースＬＵＴ５１０，５２０の入力表色系であるＲＧＢ表色系やＣＭＹＫ表色系は、いわゆる機器依存表色系では無く、特定のデバイスとは無関係に設定された仮想の表色系（あるいは抽象的な表色系）である。これらのベースＬＵＴ５１０，５２０は、例えば色補正ＬＵＴ６１０，６２０を作成する際に使用される。「ベースＬＵＴ」という名前は、色補正ＬＵＴを作成する基礎として用いられるからである。また、「ベースＬＵＴ」は本発明のプロファイル作成方法によって作成される「プロファイル」に相当する。色補正ＬＵＴ６１０，６２０は、標準的な機器依存表色系（例えばｓＲＧＢ表色系やＪＡＰＡＮＣＯＬＯＲ２００１表色系）を、特定のプリンターのインク量に変換するためのルックアップテーブルである。インバースモデル初期ＬＵＴ４１０については後述する。なお、本実施例では、シアン（Ｃ）、マゼンタ（Ｍ）、イエロー（Ｙ）、ブラック（Ｋ）の４種類のインクを利用可能なプリンターのためのＬＵＴを作成することとする。本実施例では、説明の簡略化のために４種類のインクを想定したが、他のインクについてのＬＵＴを作成する場合にも本発明を適用することができる。

図２は、本実施例の全体の処理手順を示すフローチャートである。ステップＳ０１では、ＬＵＴ作成条件設定モジュール７００のＵＩモジュール６１０が表示装置および入力装置を介してＬＵＴを作成したいメディアの指定を受け付ける。各メディアの一覧をリスト表示し、そのなかから所望のものを選択させるメディア選択ＵＩ画像を表示装置に表示させる。例えば、基本メディアとして基本光沢紙，基本マット紙，基本普通紙，基本プルーフ紙等を選択することができる。「基本メディア」とは、例えばプリンターメーカーが供給する印刷媒体を意味し、特に本発明においては分光プリンティングモデルコンバーター３１０と画質評価指数コンバーター１３６とが予め準備されている印刷媒体を意味する。基本メディアの特性は既知であり、各基本メディアにインクを付着させた場合の発色特性やデューティ制限値を特定するデータが予めメディアテーブルＭＴＢに格納されている。基本光沢紙，基本マット紙，基本普通紙，基本プルーフ紙は、それぞれ光沢紙系，マット紙系，普通紙系，プルーフ系の各系統に属する。

一方、前記メディア選択ＵＩ画像においては、基本メディアの他に、本発明の転用メディアとして、光沢紙系用紙，マット紙系用紙，普通紙系用紙，プルーフ系用紙，未分類用紙等を選択することができる。「転用メディア」とは、基本メディアのいずれとも同一でない印刷媒体を意味し、特に本発明においては分光プリンティングモデルコンバーター３１０と画質評価指数コンバーター１３６とが予め準備されていない印刷媒体を意味する。ユーザーがＬＵＴを作成しようとするメディアの系統を認知している場合には、該系統のメディアを選択することができる。メディアの系統が不明もしくは分類が困難である場合には、未分類用紙を選択することができる。なお、メディアをユーザーが指定するのではなく、メディアを測色することにより得られた色相等に基づいてメディアの系統を自動判別するようにしてもよい。ＬＵＴ作成対象のメディアが指定されると、メディアを特定する情報を設定テーブルＳＴＢに登録する（ステップＳ０２）。

ステップＳ０３において、ＬＵＴ作成条件設定モジュール７００は、指定されたメディアの系統に対応するデフォルト重みｗ_L*，ｗ_a*…をハードディスク装置等に格納されたメディアテーブルＭＴＢを参照することにより取得する。メディアの系統は、光沢紙系，マット紙系，普通紙系，プルーフ系，未分類系であり、各系統についてデフォルト重みｗ_L*，ｗ_a*がメディアテーブルＭＴＢに格納されている。ステップＳ０４において、ＵＩモジュール６１０は、重み指定用ＵＩ画像を表示装置に表示させ、入力装置により重みｗ_L*，ｗ_a*…の指定を受け付ける。

図３は重み指定用ＵＩ画像を示し、図４はメディアテーブルＭＴＢの一例を示している。重み指定用ＵＩ画像においては、後述する目的関数Ｅを構成する各項の重みｗ_L*，ｗ_a*…（０〜１００％）の指定を受け付ける。同図に示すように、粒状性と色恒常性とランニングコストとガマットと階調性の各設定項目についてスライダーバーが設けられており、ポインターが右側にスライドするほど好ましい設定がなされていることを意味する。また、スライダーバーの中央の位置が、重みｗ_L*，ｗ_a*…の中央値（５０％）に対応する。粒状性のポインターを右側にスライドさせるほど、重みｗ_GIの値が大きく設定される。色恒常性のポインターを右側にスライドさせるほど、重みｗ_CII(A)・・ｗ_CII(F12)の値が大きく設定される。なお、本実施例では、重みｗ_CII(A)・・ｗ_CII(F12)の値は互いに等しくする。むろん、光源の重要度に応じて重みｗ_CII(A)・・ｗ_CII(F12)の値を異ならせてもよい。ランニングコストのポインターを右側にスライドさせるほど、重みｗ_TIの値が大きく設定される。ガマットのポインターを右側にスライドさせるほど、重みｗ_GMIの値が大きく設定される。階調性のポインターを右側にスライドさせるほど、重みｗ_L*，ｗ_a*，ｗ_b*が大きく設定される。本実施例では、重みｗ_L*，ｗ_a*，ｗ_b*の値は互いに等しくする。なお、明度・彩度別に異なる重みｗ_L*，ｗ_a*，ｗ_b*が設定できるようにしてもよい。ポインターのスライダーバーにおける位置と各重みｗ_L*，ｗ_a*…の値の関係は、単調増加の関係にあればよく、線形関数や２次関数等の各種関数で規定することができる。

ステップＳ０４にてスライダーバーを最初に表示させる際の各ポインターの初期位置は、図４に示すメディアテーブルＭＴＢに規定されたデフォルト重みｗ_L*，ｗ_a*…に対応する位置とする。デフォルト重みｗ_L*，ｗ_a*…は、重みテーブルＷＴＢにおいてメディアの各系統について好ましい値が予め設定されている。図３に図示するように、普通紙系についてはランニングコストのデフォルト重みｗ_TIが中央値よりも大きく、ガマットのデフォルト重みｗ_GMIが中央値よりも小さくされており、それ以外は中央値とされている。マット紙系と未分類系については、すべてのデフォルト重みｗ_L*，ｗ_a*…が中央値とされている。光沢紙系については、粒状性と階調性とガマットの重みｗ_GI，ｗ_L*，ｗ_a*，ｗ_b*，ｗ_GMIが中央値よりも大きくされており、それ以外は中央値とされている。プルーフ紙系は、ガマットの重みｗ_GMIのみが中央値よりも大きくされており、それ以外は中央値とされている。ポインターの初期位置をユーザーが変更することなく決定ボタンをクリックすると、デフォルト重みｗ_L*，ｗ_a*…がそのまま設定されることとなる。

デフォルト重みｗ_L*，ｗ_a*…は、各系統のメディアの使用目的等を考慮して好適な値に設定されている、基本的には変更しなくてもよい。ユーザーが特に意図する場合には、初期位置からポインターを所望の位置にスライドさせることにより、所望の重みｗ_L*，ｗ_a*…を設定することができる。なお、各重みｗ_L*，ｗ_a*…は、相対的な大きさの差が意味をなし、全体的に一様に大小させることの意義は小さい。従って、ある項目のポインターを移動させたことにより、他の項目のポインターが逆方向に一様に移動するようにしてもよい。ステップＳ０５では、設定情報格納モジュール７３０が、決定ボタンがクリックされたときの各ポインターの位置に対応する重みｗ_L*，ｗ_a*…を設定テーブルＳＴＢに登録する。

ステップＳ０６では指定されたメディアが基本メディアであるか否かを判定し、基本メディアである場合には、ハードディスク装置等に格納されたメディアテーブルＭＴＢを参照して該基本メディアについてのデューティ制限値を取得する（ステップＳ０７）。本実施例では、ＣＭＹＫの４種類のインクを自然数の下付文字j（ｊ＝１〜４）によって区別し、メディアに付着させるインク量を個々のインクのインク量Ｉ₁〜Ｉ₄にベクトルＩ＝（Ｉ₁，Ｉ₂，Ｉ₃，Ｉ₄）によって表すこととする。なお、インク量Ｉ_j（後述するＩ_j(R,G,B)，ΔＩ_j，Ｉ_jr，ｈ_jも含む。）を下付文字ｊを付すことなく示す場合は、各インクのインク量Ｉ_jを各行要素として有する行列（ベクトル）を意味することとする。さらに、下付文字ｊ（ｊ＝５〜７）によってＣＭＹの３種類のインクを２種類ずつ混色したときの２次色のインク量を示すこととする。すなわち、Ｉ₅＝Ｉ₁＋Ｉ₂，Ｉ₆＝Ｉ₁＋Ｉ₃，Ｉ₇＝Ｉ₂＋Ｉ₃とする。２次色のインク量Ｉ₅〜Ｉ₇は、それぞれブルー（Ｂ）、レッド（Ｒ）、グリーン（Ｇ）の色相に対応する色をメディア上で再現させる。さらに、下付文字ｊ（ｊ＝８）によってＣＭＹＫの４種類のインクをすべて混色したときのインク量を示すこととする。すなわち、Ｉ₈＝Ｉ₁＋Ｉ₂＋Ｉ₃＋Ｉ₄とする。

本実施例では、各インクのインク量Ｉ_jは８ビットで表現されることとする。図４に示すように、デューティ制限値Ｄ_Ijは、各インク単独（１次色）と、２次色のインクの合計と、全インクの合計について記憶されている。デューティ制限値Ｄ_Ijは、各基本メディアに対して最大限付着させることができるインク量を意味し、例えばインクにじみが生じる下限値が設定される。インク滴のメディア上における物理的性質は、インクとメディアの組み合わせごとに異なっており、これらの組み合わせごとに異なるデューティ制限値Ｄ_Ijが設定されている。また、複数インクを混色する場合にも、単一インクとは異なる物理的性質を示すため、本実施例では１次色だけでなく、２次色（２インク混色）と全インク合計についてもデューティ制限値Ｄ_Ij（ｊ＝１〜８）が設定されている。基本メディアについてデューティ制限値Ｄ_Ijが取得できると、ステップＳ０８において、設定情報格納モジュール７３０が、該取得したデューティ制限値Ｄ_Ijを設定テーブルＳＴＢに格納させる。さらに、ステップＳ０８においては、インク量コンバーター７１０を無効とする旨の無効フラグを設定テーブルＳＴＢに格納させる。本明細書において、単にインク量Ｉ_jと表記した場合には下付文字ｊの範囲はｊ＝１〜４であり、デューティ制限値Ｄ_Ijと表記した場合には下付文字ｊの範囲はｊ＝１〜８であるとする。

図５は、設定テーブルＳＴＢの一例を示している。基本メディアが指定された場合には、設定テーブルＳＴＢに、上述した指定された基本メディアの種類と、重みｗ_L*，ｗ_a*…と、デューティ制限値Ｄ_Ijと、無効フラグと、グレーターゲットの色味（ａ_gt ^*，ｂ_gt ^*）が記憶される。基本メディアの場合、グレーターゲットの色味（ａ_gt ^*，ｂ_gt ^*）＝（０，０）とされる。一方、転用メディアが指定された場合には、以上とは異なる処理（ステップＳ０９〜）が実行されるが、ここではまず基本メディアが指定された場合について、ＬＵＴを作成する処理をひととおり説明することとする。

Ｂ．基本メディアのＬＵＴ作成手順
Ｂ−１．全体手順
図６は、実施例において基本メディアのベースＬＵＴを作成する手順を示すフローチャートである。図７（Ａ）〜（Ｃ）は、図６のステップＳ１００〜Ｓ３００によってベース３Ｄ−ＬＵＴを作成する場合の処理内容を示す説明図である。ステップＳ１００では、設定テーブルＳＴＢに記憶された情報に基づいて、フォワードモデルコンバーター３００（分光プリンティングモデルコンバーター３１０）とインバースモデル初期ＬＵＴ４１０と画質評価指数コンバーター１３６とが準備（起動）される。上述したように、基本メディアについては、該基本メディアのための分光プリンティングモデルコンバーター３１０と画質評価指数コンバーター１３６とが準備されているため、これらを起動して使用可能とすることができる。また、無効フラグが添付されているため、インク量コンバーター７１０は起動させない。ここで、「フォワードモデル」とは、インク量Ｉ_jを機器独立表色系の色彩値（測色値）に変換する変換モデルを意味し、「インバースモデル」とは、逆に、機器独立表色系の色彩値をインク量に変換する変換モデルを意味している。実施例では、機器独立表色系としてＣＩＥ−Ｌａｂ表色系を使用する。なお、以下では、ＣＩＥ−Ｌａｂ表色系の色彩値を、単に「Ｌ^*ａ^*ｂ^*値」または「Ｌａｂ値」とも呼ぶ。

図７（Ａ）に示すように、フォワードモデルコンバーター３００の前段を構成する分光プリンティングモデルコンバーター３１０は、複数（４）種類のインクのインク量Ｉ_jを基本メディアに付着させて形成したカラーパッチの分光反射率Ｒ（λ）に変換する。なお、本明細書において「カラーパッチ」という用語は、有彩色のパッチに限らず、無彩色のパッチも含む広い意味で使用される。また、「印刷する」とは、メディアにインク量に応じてインクを付着させることを指す。分光プリンティングモデルコンバーター３１０は、上述した４種類のインクのインク量Ｉ_jを入力としている。色算出部３２０は、分光反射率Ｒ（λ）からＬａｂ表色系の色彩値を算出する。この色彩値の算出には、予め選択された光源（例えば標準の光Ｄ５０）がカラーパッチの観察条件として使用される。なお、分光プリンティングモデルコンバーター３１０を作成する方法としては、例えば特表２００７−５１１１７５号公報に記載された方法を採用することが可能である。

インバースモデル初期ＬＵＴ４１０は、Ｌ^*ａ^*ｂ^*値を入力とし、インク量Ｉ_j出力とするルックアップテーブルである。この初期ＬＵＴ４１０は、例えば、Ｌ^*ａ^*ｂ^*空間を複数の小セルに区分し、各小セル毎に最適なインク量Ｉ_jを選択して登録したものである。この選択は、例えば、そのインク量Ｉ_jで基本メディアに印刷されるカラーパッチの画質を考慮して行われる。一般に、ある１つのＬ^*ａ^*ｂ^*値を再現するインク量Ｉ_jの組み合わせは多数存在する。そこで、初期ＬＵＴ４１０では、ほぼ同じＬ^*ａ^*ｂ^*値を再現する多数のインク量Ｉ_jの組み合わせの中から、画質等の所望の観点から最適なインク量を選択したものが登録されている。この初期ＬＵＴ４１０の入力値であるＬ^*ａ^*ｂ^*値は各小セルの代表値である。一方、出力値であるインク量Ｉ_jはそのセル内のいずれかのＬ^*ａ^*ｂ^*値を再現するものである。従って、この初期ＬＵＴ４１０では、入力値であるＬ^*ａ^*ｂ^*値と出力値であるインク量Ｉ_jとが厳密に対応したものとなっておらず、出力値のインク量をフォワードモデルコンバーター３００でＬ^*ａ^*ｂ^*値に変換すると、初期ＬＵＴ４１０の入力値とは多少異なる値が得られる。ただし、初期ＬＵＴ４１０として、入力値と出力値とが完全に対応するものを利用してもよい。また、初期ＬＵＴ４１０を用いずにベースＬＵＴを作成することも可能である。なお、小セル毎に最適なインク量を選択して初期ＬＵＴ４１０を作成する方法としては、例えば前記特表２００７−５１１１７５号公報に記載された方法を採用することが可能である。なお、前記特表２００７−５１１１７５号公報の方法では、対象の印刷媒体にカラーパッチを形成することにより、分光プリンティングモデルコンバーター３１０とインバースモデル初期ＬＵＴ４１０とが作成される。すなわち、基本メディアのベースＬＵＴを作成するためには、基本メディアにカラーパッチを形成することにより作成された分光プリンティングモデルコンバーター３１０とインバースモデル初期ＬＵＴ４１０とを準備することとなる。

図６のステップＳ２００では、ベースＬＵＴ作成のための初期入力値がユーザーによって設定される。図７（Ｂ）は、ベース３Ｄ−ＬＵＴ５１０の構成とその初期入力値設定の例を示している。ベース３Ｄ−ＬＵＴ５１０の入力値としては、ＲＧＢの各値として予め定められたほぼ等間隔の値が設定される。１組のＲＧＢ値はＲＧＢ色空間内の点を表していると考えられるので、１組のＲＧＢ値を「入力格子点」とも呼ぶ。ステップＳ２００においては、複数の入力格子点のうちから予め選択されたいくつかの少数の入力格子点に対するインク量Ｉ_jの初期値がユーザーによって入力される。本実施例では、ＲＧＢの各値を８ビットで表現した場合に、（Ｒ，Ｇ，Ｂ）＝（１６ｎ₁−１，１６ｎ₂−１，１６ｎ₃−１）を満足するすべて（１７³個）の入力格子点を選択することとする。なお、ｎ₁〜ｎ₃は、それぞれ０〜１６の整数であり、Ｒ，Ｇ，Ｂ＝−１のときはＲ，Ｇ，Ｂ＝０とする。この初期入力値が設定される入力格子点には、ＲＧＢ色空間における３次元色立体の頂点に相当する入力格子点が含まれることとなる。この３次元色立体の頂点では、ＲＧＢの各値がその定義範囲の最小値または最大値を取る。具体的には、（Ｒ，Ｇ，Ｂ）＝（０，０，０）、（０，０，２５５）、（０，２５５，０）、（２５５，０，０）、（０，２５５，２５５）、（２５５，０，２５５）、（２５５，２５５，０）、（２５５，２５５，２５５）である８つの入力格子点に関してインク量Ｉ_jの初期入力値が設定される。また、ｎ₁＝ｎ₂＝ｎ₃となる１７個の入力格子点（以下、グレー格子点と表記する。）は、ＲＧＢ色空間上のグレー軸上に存在することとなる。なお、（Ｒ，Ｇ，Ｂ）＝（２５５，２５５，２５５）の入力格子点に対するインク量Ｉ_jは、すべて０に設定される。他の入力格子点に対するインク量Ｉ_jの初期入力値は任意であり、例えば０に設定される。図７（Ｂ）の例では、（Ｒ，Ｇ，Ｂ）＝（０，０，３２）の入力格子点に対するインク量が０以外の値になっているが、これはこのＬＵＴ５１０が完成したときの値である。

図６のステップＳ３００では、スムージング処理モジュール１３０（図１）が、ステップＳ２００で設定された初期入力値に基づいてスムージング処理（平滑化／最適化処理）を実行する。図７（Ｃ）は、ステップＳ３００の処理内容を示している。図７（Ｃ）の左側には、スムージング処理前の状態における複数の色彩値の分布が２重丸と白丸とで示されている。これらの色彩値は、Ｌ^*ａ^*ｂ^*空間における３次元色立体ＣＳを構成している。各色彩値のＬ^*ａ^*ｂ^*座標値は、ベース３Ｄ−ＬＵＴ５１０の複数の入力格子点におけるインク量Ｉ_jを、フォワードモデルコンバーター３００（図７（Ａ））を用いてＬ^*ａ^*ｂ^*値に変換した値である。上述したように、ステップＳ２００では一部の少数の入力格子点についてのみインク量Ｉ_jの初期入力値が設定される。そこで、他の入力格子点に対するインク量の初期値は、初期入力値からスムージング処理初期値設定モジュール１２０（図１）によって設定される。この初期値設定方法については後述する。

Ｌａｂ表色系の３次元色立体ＣＳは、以下の８つの頂点（図７（Ｃ）の２重丸の点）を有している。
・点Ｐ_K：（Ｒ，Ｇ，Ｂ）＝（０，０，０）に対応する紙黒点。
・点Ｐ_W：（Ｒ，Ｇ，Ｂ）＝（２５５，２５５，２５５）に対応する紙白点。
・点Ｐ_C：（Ｒ，Ｇ，Ｂ）＝（０，２５５，２５５）に対応するシアン点。
・点Ｐ_M：（Ｒ，Ｇ，Ｂ）＝（２５５，０，２５５）に対応するマゼンタ点。
・点Ｐ_Y：（Ｒ，Ｇ，Ｂ）＝（２５５，２５５，０）に対応するイエロー点。
・点Ｐ_R：（Ｒ，Ｇ，Ｂ）＝（２５５，０，０）に対応するレッド点。
・点Ｐ_G：（Ｒ，Ｇ，Ｂ）＝（０，２５５，０）に対応するグリーン点。
・点Ｐ_B：（Ｒ，Ｇ，Ｂ）＝（０，０，２５５）に対応するブルー点。

図７（Ｃ）の右側は、スムージング処理後の格子点（色彩値）の分布を示している。スムージング処理は、Ｌ^*ａ^*ｂ^*空間における複数の格子点を移動させて、それらの格子点の分布を等間隔に近い平滑なものにする処理である。スムージング処理では、さらに、移動後の各格子点のＬ^*ａ^*ｂ^*値を再現するために最適なインク量Ｉ_jも決定される。この最適なインク量がベースＬＵＴ５１０の出力値として登録されると、ベースＬＵＴ５１０が完成する。

図８（Ａ）〜（Ｃ）は、入力表色系の格子点（すなわち入力格子点）とＬａｂ表色系の格子点との対応関係を示している。Ｌａｂ表色系の３次元色立体ＣＳの頂点は、ベースＬＵＴ５１０の入力表色系の３次元色立体の頂点と一対一に対応している。また、各頂点を結ぶ辺（稜線）も、両方の色立体で互いに対応しているものと考えることができる。スムージング処理前のＬａｂ表色系の各格子点の色彩値は、ベースＬＵＴ５１０の入力格子点にそれぞれ対応付けられており、従って、スムージング処理後のＬａｂ表色系の各格子点の色彩値もベースＬＵＴ５１０の入力格子点にそれぞれ対応付けられる。なお、ベースＬＵＴ５１０の入力格子点はスムージング処理によって変化しない。スムージング処理後のＬａｂ表色系の３次元色立体ＣＳは、ベースＬＵＴ５１０の出力表色系を構成するインクセットで再現可能な色域（ガマット）の全体に対応している。従って、ベースＬＵＴ５１０の入力表色系は、このインクセットで再現可能な色域の全体を表す表色系としての意義を有している。

ベースＬＵＴ５１０を作成する際に、Ｌ^*ａ^*ｂ^*空間においてスムージング処理を行う理由は以下の通りである。ベースＬＵＴ５１０では、なるべく大きな色域を再現できるように出力表色系のインク量Ｉ_jを設定したいという要望がある。特定のインクセットで再現可能な色域は、デューティ制限値Ｄ_Ij等に依存することとなる。一方、上述したフォワードモデルコンバーター３００は、これらの制限条件が考慮されておらず、再現可能な色域とは無関係に作成されている。そこで、スムージング処理の際にデューティ制限値Ｄ_Ij等の制限条件を考慮してＬ^*ａ^*ｂ^*空間内の色彩値の取り得る範囲を決定すれば、特定のインクセットで再現可能な色域を決定することが可能となる。なお、格子点の移動を行うアルゴリズムとしては、例えば、後述する力学モデルを使用したものが利用される。

図６のステップＳ４００では、スムージング処理の結果を用いて、テーブル作成モジュール１４０がベースＬＵＴ５１０を作成する。すなわち、テーブル作成モジュール１４０は、各入力格子点に対応付けられたＬａｂ表色系の格子点の色彩値を再現するための最適なインク量Ｉ_jをベースＬＵＴ５１０の出力値として登録する。なお、スムージング処理では、その計算負荷を軽減するために、ベースＬＵＴ５１０の入力格子点の一部のみに対応する格子点の色彩値のみを処理対象として選択することも可能である。例えば、ベースＬＵＴ５１０の入力格子点におけるＲＧＢ値の間隔が１６である場合に、スムージング処理の対象となる入力格子点におけるＲＧＢ値の間隔を３２に設定すれば、スムージング処理の負荷を半減することができる。この場合には、テーブル作成モジュール１４０は、スムージング処理結果を補間することによってベースＬＵＴ５１０のすべての入力格子点に対するインク量Ｉ_jを決定して登録する。

図９（Ａ）〜（Ｃ）は、図６のステップＳ１００〜Ｓ３００によってベース４Ｄ−ＬＵＴ５２０を作成する場合の処理内容を示す説明図である。図９（Ａ）は、図７（Ａ）と同じである。図９（Ｂ）に示すベース４Ｄ−ＬＵＴ５２０は、入力がＣＭＹＫ表色系である点が図７（Ｂ）に示したベース３Ｄ−ＬＵＴ５１０と異なっている。このベース４Ｄ−ＬＵＴ５２０の初期入力値としては、（Ｃ，Ｍ，Ｙ，Ｋ）＝（０，０，０，０），（０，０，２５５，０），（０，２５５，０，０），（０，２５５，２５５，０），（２５５，０，０，０），（２５５，０，２５５，０），（２５５，２５５，０，０），（２５５，２５５，２５５，０），（０，０，０，２５５），（０，０，２５５，２５５），（０，２５５，０，２５５），（０，２５５，２５５，２５５），（２５５，０，０，２５５），（２５５，０，２５５，２５５），（２５５，２５５，０，２５５），（２５５，２５５，２５５，２５５）である１６個の入力格子点に関してインク量の初期値が設定される。他の入力格子点に対するインク量の初期入力値は任意であり、例えば０に設定される。本実施例では、ベース４Ｄ−ＬＵＴ５２０を作成する場合においても、Ｃ＝Ｍ＝Ｙとなるグレー軸上のグレー格子点が入力格子点として１７個含まれることとする。

図９（Ｃ）は、スムージング処理の様子を示している。なお、Ｌ^*ａ^*ｂ^*空間においてベース４Ｄ−ＬＵＴ５２０に対応する色立体としては、図９（Ｃ）の右端に示すように、入力値のうちのＫ値のそれぞれの値に対して１つの３次元色立体ＣＳが存在する。この例では、Ｋ＝０の色立体とＫ＝３２の色立体とを含む複数の色立体ＣＳが図示されている。本明細書では、これらの個々の色立体ＣＳを「Ｋレイヤ」とも呼ぶ。この理由は、各色立体ＣＳが、ＣＭＹＫ値のうちのＫ値が一定でＣ，Ｍ，Ｙ値が可変である入力層に対応するものと考えることができるからである。複数の色立体ＣＳは、Ｋ値が大きいほど暗い色域を表現するものとなっている。これらの複数の色立体ＣＳは、入力表色系のＫ値が大きいほど濃ブラックインクＫのインク量Ｉ₄が多くなるように濃ブラックインクＫのインク量を決定することによって実現できる。上述したように、再現可能な色域はデューティ制限値Ｄ_Ij等によって制限される。このデューティ制限値Ｄ_Ijは、指定されたメディアの種類に依存する。一方、暗い色を再現する方法としては、濃ブラックインクＫなどの無彩色インクを用いる方法と、コンポジットブラックを用いる方法とがある。しかし、コンポジットブラックは合計インク量が多くなるので、濃ブラックインクＫに比べてデューティ制限値に抵触する可能性が高く、暗い色を再現するのには不利である。従って、入力表色系のＫ値が大きく濃ブラックインクＫのインク量Ｉ₄が多い色立体の方が、入力表色系のＫ値が小さく濃ブラックインクＫのインク量が少ない色立体に比べてより暗い色を再現することが可能となる。

図１０（Ａ），（Ｂ）は、ベースＬＵＴを用いた色補正ＬＵＴの作成方法を示す説明図である。図１０（Ａ）に示すように、ベース３Ｄ−ＬＵＴ５１０は、ＲＧＢ値をインク量Ｉ_jに変換する。変換後のインク量Ｉ_jは、フォワードモデルコンバーター３００によってＬ^*ａ^*ｂ^*値に変換される。一方、ｓＲＧＢ値は、既知の変換式に従ってＬ^*ａ^*ｂ^*値に変換される。この変換後のＬ^*ａ^*ｂ^*値は、その色域が、フォワードモデルコンバーター３００で変換されたＬ^*ａ^*ｂ^*値の色域と一致するようにガマットマッピングされる。一方、ベース３Ｄ−ＬＵＴ５１０とフォワードモデルコンバーター３００を通じて、ＲＧＢ値から変換したＬ^*ａ^*ｂ^*値を、逆方向ルックアップテーブルとして、逆変換ＬＵＴ５１１を作成する。ガマットマッピングされたＬ^*ａ^*ｂ^*値は、この逆変換ＬＵＴ５１１によってＲＧＢ値に変換される。このＲＧＢ値は、さらに、ベース３Ｄ−ＬＵＴ５１０によってインク量Ｉ_jに再度変換される。この最後のインク量Ｉ_jと最初のｓＲＧＢ値の対応関係をルックアップテーブルに登録することによって、色補正３Ｄ−ＬＵＴ６１０を作成することができる。この色補正３Ｄ−ＬＵＴ６１０は、ｓＲＧＢ表色系をインク表色系に変換する色変換テーブルである。

図１０（Ｂ）は、色補正４Ｄ−ＬＵＴ６２０を作成する方法を示している。図１０（Ａ）との違いは、３Ｄ−ＬＵＴ５１０およびその逆変換ＬＵＴ５１１の代わりに、４Ｄ−ＬＵＴ５２０およびその逆変換ＬＵＴ５２１を利用している点と、ｓＲＧＢ表色系をＬ^*ａ^*ｂ^*値に変換する既知変換式の代わりにＪＡＰＡＮＣＯＬＯＲ表色系（図中では「ｊＣＭＹＫ」と記したもの）をＬ^*ａ^*ｂ^*値に変換する既知変換式を使用している点である。良く知られているように、ＪＡＰＡＮＣＯＬＯＲは、ＣＭＹＫの４色で構成される表色系である。なお、図１０（Ｂ）の方法では、逆変換ＬＵＴ５２１において、Ｌ^*ａ^*ｂ^*値からＣＭＹＫ値に変換する際に、既知変換前の最初のｊＣＭＹＫ値のＫ値から、逆変換ＬＵＴ５２１のＫレイヤ（Ｋ値が一定を取る部分）が選択される。従って、色補正４Ｄ−ＬＵＴ６２０として、ベース４Ｄ−ＬＵＴ５２０のうちのＫレイヤにおける特性を反映したものを作成することが可能である。なお、設定情報格納モジュール７３０は、作成したベースＬＵＴ５１０，５２０や色補正ＬＵＴ６１０，６２０に設定テーブルＳＴＢを添付する。これにより、ベースＬＵＴ５１０，５２０や色補正ＬＵＴ６１０，６２０がどのようなメディアについて作成されたものか、どのような重みｗ_L*，ｗ_a*…で作成されたものであるかを識別することができる。

なお、通常は、ベースＬＵＴ５１０，５２０がプリンタドライバに実装されており、色補正ＬＵＴの作成処理以外の処理にも活用されているが、ここでは他の活用例の説明は省略する。以下では、実施例のスムージング処理（平滑化／最適化処理）に利用される力学モデルについて簡単に説明した後に、スムージング処理の処理手順、および、最適化処理の内容について順次説明する。

Ｂ−２．力学モデル：
図１１は、本実施例のスムージング処理（平滑化／最適化処理）に利用される力学モデルを示す説明図である。ここでは、Ｌ^*ａ^*ｂ^*色空間内に上述した入力格子点に対応する格子点（白丸および２重丸）が配列されている様子を示している。ただし、ここでは説明の便宜上、格子点の配置を２次元的に描いている。この力学モデルでは、着目格子点ｇに対して次式の仮想的な力Ｆp_gが係るものと仮定する。

ここで、Ｆ_gは着目格子点ｇが隣接格子点ｇｎ（ｎは１〜Ｎ）から受ける引力の合計値、Ｖ_gは着目格子点ｇの速度ベクトル、−ｋ_vＶ_gは速度に応じた抵抗力、Ｘ_gは着目格子点ｇの位置ベクトル、Ｘ_gnは隣接格子点ｇｎの位置ベクトル、ｋ_p，ｋ_gは係数である。係数ｋ_p，ｋ_gは予め一定の値に設定される。なお、文中では、ベクトルを示す矢印は省略される。

このモデルは、バネで互いに結ばれた質点の減衰振動モデルである。すなわち、着目格子点ｇに係る仮想合力Ｆp_gは、着目格子点ｇと隣接格子点ｇｎとの距離が大きいほど大きくなるバネ力Ｆ_gと、着目格子点ｇの速度が大きいほど大きくなる抵抗力−ｋ_vＶ_gとの合計値である。この力学モデルでは、各色彩値について、位置ベクトルＸ_gと速度ベクトルＦ_gの初期値を設定した後に、微小時間経過の位置ベクトルＸ_gと速度ベクトルＦ_gを順次算出してゆく。なお、複数の色彩値の速度ベクトルＶ_gの初期値は、例えば０に設定される。このような力学モデルを利用すれば、各色彩値が徐々に移動して、平滑な色彩値分布を得ることが可能である。

なお、各色彩値に係る力としては、バネ力Ｆ_gと抵抗力−ｋ_vＶ_g以外の力を用いても良い。例えば、本出願人により開示された特開２００６−１９７０８０号公報で説明されている他の種々の力をこの力学モデルで利用してもよい。また、力学モデルを適用して各色彩値を移動させる際に、特定の色彩値は、力学モデルによって移動しない拘束点として取り扱うことも可能である。本実施例では、上述した１７個のグレー軸格子点に対応する格子点の色彩値が、グレーターゲットの色味（ａ_gt ^*，ｂ_gt ^*）が示す色相方向にずれるように拘束する。基本メディアの場合には、グレーターゲットの色味（ａ_gt ^*，ｂ_gt ^*）＝（０，０）としていされているため、グレー軸格子点に対応する格子点の色彩値は、Ｌ^*ａ^*ｂ^*色空間のＬ^*軸上の位置を示すように拘束される。

図１２は、グレー軸格子点に対応する格子点の位置（色彩値）がグレーターゲットに拘束される様子を示している。同図に示すように、グレー軸格子点に対応する格子点（白丸）が、Ｌ^*ａ^*ｂ^*色空間のＬ^*軸上の紙黒点と紙白点とを結ぶ線分（グレーターゲット）を１６等分する位置を示すように拘束されている。なお、特開２００６−２１７１５０公報に開示された手法によって、特定の格子点をＬ^*ａ^*ｂ^*色空間の特定の位置に拘束することができる。このように、Ｌ^*ａ^*ｂ^*色空間のＬ^*軸上において均等に格子点を拘束することにより、スムージング処理後のベースＬＵＴ５１０，５２０を用いて色補正ＬＵＴ６１０，６２０を作成する際に行われるグレー軸まわりの補間演算の精度を向上させることができる。従って、グレー軸まわりの色再現性や階調性に優れた色補正ＬＵＴ６１０，６２０を作成することができる。

Ｂ−３．スムージング処理（平滑化／最適化処理）の処理手順：
図１３は、スムージング処理（図６のステップＳ３００）の典型的な処理手順を示すフローチャートである。ステップＴ１００では、スムージング処理初期値設定モジュール１２０（図１）が、スムージング処理の対象とする複数の格子点を初期設定する。

図１４は、ステップＴ１００の詳細手順を示すフローチャートである。ステップＴ１０２では、インク量の初期入力値（図７（Ｂ），図９（Ｂ））から、スムージング処理の対象となる各格子点の仮インク量が決定される。例えば、３Ｄ−ＬＵＴ用のスムージング処理では、次の（２）式、（３）式に従って、各入力格子点に対する仮インク量Ｉ_(R,G,B)が決定される。

ここで、Ｉ_(R,G,B)は、入力格子点のＲＧＢ値に対するインクセット全体のインク量Ｉ_j（図７の例では４種類のインクのインク量Ｉ_j）を表している。ＲＧＢ値が０または２５５を取る入力格子点に対するインク量は、図６のステップＳ２００においてユーザーによって予め入力された値である。前記（２）式および（３）式によれば、任意のＲＧＢ値における仮インク量Ｉ_(R,G,B)を求めることが可能である。

４Ｄ−ＬＵＴ用のスムージング処理では、次の（４）式、（５）式に従って、各入力格子点に対する仮インク量Ｉ_(C,M,Y,K)が決定される。

なお、（４）式からも理解できるように、４Ｄ−ＬＵＴ用のインク量の初期入力値は１６個存在するので、初期入力値の設定が煩雑である。そこで、例えば、インク量の初期入力値を設定する入力格子点を、Ｋ＝０の８個の頂点、すなわち、（Ｃ，Ｍ，Ｙ，Ｋ）＝（０，０，０，０），（０，０，２５５，０），（０，２５５，０，０），（０，２５５，２５５，０），（２５５，０，０，０），（２５５，０，２５５，０），（２５５，２５５，０，０），（２５５，２５５，２５５，０）の８個の頂点と、Ｋ＝２５５の１個の頂点、例えば、（Ｃ，Ｍ，Ｙ，Ｋ）＝（０，０，０，２５５）の頂点のみとし、Ｋ＝２５５の格子点のインク量を次の（６）式または（７）式で決定してもよい。

ここで、Ｉ_(C,M,Y,0)は、Ｋ＝０の８個の頂点におけるインク量の初期入力値から、前記（２）式と同様の式で算出されたインク量である。（６）式の関数ｆ_D1は値Ｉ_(C,M,Y,0)と値Ｉ_(0,0,0,255)の合計値がデューティ制限値Ｄ_I8をオーバーする場合に、値Ｉ_(C,M,Y,0)を減じることによって、インク量Ｉ_(C,M,Y,255)がデューティ制限値Ｄ_I8内に納まるようにする関数である。また、（７）式の関数ｆ_D2は、値Ｉ_(C,M,Y,0)と値Ｉ_(0,0,0,255)の合計値がデューティ制限値Ｄ_I8をオーバーする場合に、合計値（Ｉ_(C,M,Y,0)＋Ｉ_(0,0,0,255)）の全体を減じることによって、インク量Ｉ_(C,M,Y,255)がデューティ制限値Ｄ_I8内に納まるようにする関数である。

図１４のステップＴ１０４では、フォワードモデルコンバーター３００を用いて、仮インク量に対応する色彩値Ｌ^*ａ^*ｂ^*を求める。この演算は、以下の（８）式または（９）式で表すことができる。

ここで、Ｌ^* _(R,G,B)、ａ^* _(R,G,B)、ｂ^* _(R,G,B)、Ｌ^* _(C,M,Y,K)、ａ^* _(C,M,Y,K)、ｂ^* _(C,M,Y,K)は変換後の色彩値Ｌ^*ａ^*ｂ^*を示しており、関数ｆ_L*FM、ｆ_a*FM、ｆ_b*FMはフォワードモデルコンバーター３００による変換を意味している。なお、これらの式からも理解できるように、この変換後の色彩値Ｌ^*ａ^*ｂ^*は、ベースＬＵＴの入力値であるＲＧＢ値またはＣＭＹＫ値に対応付けられている。

図１４のステップＴ１０６では、ステップＴ１０４で得られた色彩値Ｌ^*ａ^*ｂ^*を、インバースモデル初期ＬＵＴ４１０を用いてインク量に再度変換する。ここで、インバースモデル初期ＬＵＴ４１０を用いてインク量に再度変換する理由は、インク量Ｉ_jの初期入力値や、ステップＴ１０２で決定された仮インク量が、Ｌ^*ａ^*ｂ^*値を再現するインク量として必ずしも好ましいインク量Ｉ_jでは無いからである。一方、インバースモデル初期ＬＵＴ４１０では、画質等を考慮した好ましいインク量が登録されているので、これを用いてＬ^*ａ^*ｂ^*値をインク量Ｉ_jに再度変換すれば、そのＬ^*ａ^*ｂ^*値を実現するための好ましいインク量Ｉ_jを初期値として得ることができる。ただし、ステップＴ１０６を省略してもよい。ステップＴ１０７では、上述したＬ^*軸上のグレーターゲットを設定する。

上述のステップＴ１００の処理の結果、スムージング処理の対象となる色彩値について、以下の初期値が決定される。
（１）ベースＬＵＴの入力格子点の値：（Ｒ，Ｇ，Ｂ）または（Ｃ，Ｍ，Ｙ，Ｋ）
（２）各入力格子点に対応するＬ^*ａ^*ｂ^*空間の格子点の初期座標値：（Ｌ^* _(R,G,B)，ａ^* _(R,G,B)，ｂ^* _(R,G,B)）または（Ｌ^* _(C,M,Y,K)，ａ^* _(C,M,Y,K)，ｂ^* _(C,M,Y,K)）
（３）各入力格子点に対応する初期インク量：Ｉ_(R,G,B)またはＩ_(C,M,Y,K)
以上の説明から理解できるように、初期値設定モジュール１２０は、代表的な入力格子点に関する入力初期値から他の入力格子点に関する初期値を設定する機能を有している。なお、初期値設定モジュール１２０は、スムージング処理モジュール１３０に含まれるものとしてもよい。

図１３のステップＴ１２０では、色点移動モジュール１３２が、上述した力学モデルに従ってＬ^*ａ^*ｂ^*空間内の色彩値を移動させる。

図１５（Ａ）〜（Ｄ）は、図１３のステップＴ１２０〜Ｔ１５０の処理内容を示す説明図である。図１５（Ａ）に示すように、スムージング処理前には、格子点の分布にはかなりの偏りがある。図１５（Ｂ）は、微小時間経過後の各格子点の位置を示している。この移動後の各色彩値のＬ^*ａ^*ｂ^*値を「ターゲット値（Ｌ_t ^*ａ_t ^*ｂ_t ^*値）」と呼ぶ。「ターゲット」という修飾語は、このＬ_t ^*ａ_t ^*ｂ_t ^*値が、以下で説明するインク量の最適値の探索処理の際の目標値として使用されるからである。

ステップＴ１３０では、インク量最適化モジュール１３４が、予め設定された目的関数Ｅを用いて、ターゲット値ＬＡＢ_tに対するインク量Ｉ_jの最適値を探索する（図１１（Ｃ）参照）。この目的関数Ｅを用いた最適化では、力学モデルで微小量だけ移動した後の色彩値の座標値ＬＡＢ_tに近いＬ^*ａ^*ｂ^*値を再現するインク量Ｉ_jの中で、複数のパラメータΔＬ^*，Δａ^*，ΔＧＩ，ΔＣＩＩ，ΔＴＩの２乗誤差の和がなるべく小さいインク量が最適なインク量Ｉ_jとして決定される。また、最適なインク量Ｉ_jの探索は、ステップＴ１００で設定された各入力格子点の初期インク量から開始される。従って、探索で得られるインク量Ｉ_jは、この初期インク量を修正した値となる。後で詳述するように、（ＥＱ１）式で与えられる目的関数Ｅは、（ＥＱ２）式のようなインク量ベクトルＩに関する２次形式の関数として書き表すことができる。インク量Ｉ_jの最適化は、このような２次形式の目的関数Ｅを用いて、２次計画法に従って実行される。なお、ステップＴ１３０の詳細手順や目的関数Ｅの内容については後述する。

図１３のステップＴ１４０では、ステップＴ１３０で探索されたインク量Ｉ_jに対応するＬ^*ａ^*ｂ^*値が、フォワードモデルコンバーター３００で再計算される（図１０（Ｄ）参照）。ここで、Ｌ^*ａ^*ｂ^*値を再計算する理由は、探索されたインク量Ｉ_jが目的関数Ｅを最小とするインク量なので、そのインク量Ｉ_jで再現されるＬ^*ａ^*ｂ^*値は、最適化処理のターゲット値ＬＡＢ_tから多少ずれているからである。こうして再計算されたＬ^*ａ^*ｂ^*値が、各格子点の移動後の座標値として使用される。

ステップＴ１５０では、各格子点の色彩値の移動量の平均値（ΔＬａｂ）_aveが、予め設定された閾値ε以下であるか否かが判定される。移動量の平均値（ΔＬａｂ）_aveが閾値εよりも大きい場合には、ステップＴ１２０に戻りステップＴ１２０〜Ｔ１５０のスムージング処理が継続される。一方、移動量の平均値（ΔＬａｂ）_aveが閾値ε以下の場合には、色彩値の分布が十分に平滑になっているので、スムージング処理が終了する。なお、閾値εは、予め適切な値が実験的に決定される。

このように、本実施例の典型的なスムージング処理（平滑化／最適化処理）では、力学モデルによって各格子点を微小時間間隔毎に移動させつつ、移動後の格子点に対応する最適なインク量Ｉ_jを最適化手法で探索する。そして、格子点の移動量が十分に小さくなるまでそれらの処理が継続される。この結果、図７（Ｃ）に示したように、スムージング処理によって、平滑な格子点分布を得ることが可能である。

Ｂ−４．最適化処理の内容：
最適化処理の目的関数Ｅ（図１５（Ｃ）参照）は、インク量の関数である色彩値（Ｌ^*ａ^*ｂ^*値）および画質評価指数に関するヤコビ行列Ｊを用いて表現することが可能である。各画質評価指数は、画質評価指数コンバーター１３６によって算出される。後述するように、各画質評価指数は、各インク量のインクを基本メディアに付着させた場合の画質を評価する指数である。ヤコビ行列Ｊは、例えば以下の（１０）式で表される。

（１０）式の右辺の第１行〜第３行は、色彩値Ｌ^*ａ^*ｂ^*を個々のインク量Ｉ_jで偏微分した値を示している。また、第４行以下は、１組のインク量Ｉ_j（ｊ＝１〜８）で印刷されるカラーパッチの画質を表す画質評価指数（粒状性指数ＧＩ（Graininess Index）と、色非恒常性指数ＣＩＩ（Color Inconstancy Index）と、ガマット評価指数ＧＭＩと、合計インク量ＴＩを個々のインク量Ｉ_jで偏微分した値を示している。なお、画質評価指数ＧＩ，ＣＩＩ，ＧＭＩ，ＴＩは、その値が小さいほど、インク量Ｉ_jで再現されるカラーパッチの画質が良い傾向にあることを示す指数である。

色彩値Ｌ^*ａ^*ｂ^*は、フォワードモデルコンバーター３００を用いて、以下の（１１）式でインク量Ｉ_jから変換される。

画質評価指数ＧＩ，ＣＩＩ，ＴＩ，ＧＭＩも、一般に１次色のインク量Ｉ_j（ｊ＝１〜４）の関数としてそれぞれ表現できる。

なお、（１３）式の色非恒常性指数ＣＩＩ_illの下付文字「ill」は、光源の種類を表している。上述した（１０）式では、光源の種類として、標準の光Ａと標準の光Ｆ１２とを用いている。なお、色非恒常性指数ＣＩＩの計算方法の例は後述するが、色非恒常性指数ＣＩＩとしては一種類または複数種類の任意の標準光源に関するものを利用することが可能である。

粒状性指数ＧＩは、各種の粒状性予測モデルを用いて算出可能であり、例えば以下の（１６）式で算出することができる。

ここで、ａＬは明度補正係数、ＷＳ（ｕ）はカラーパッチの印刷に利用されるハーフトーンデータが示す画像のウイナースペクトラム、ＶＴＦ（ｕ）は視覚の空間周波数特性、ｕは空間周波数である。ハーフトーンデータは、カラーパッチのインク量Ｉ_jからハーフトーン処理（プリンター１０が実行するハーフトーン処理と同一のものとする）によって決定される。前記（１５）式は一次元で表現しているが、空間周波数の関数として二次元画像の空間周波数を算出することは容易である。粒状性指数ＧＩの計算方法としては、例えば、本出願人により開示された特開２００６−１０３６４０号公報に記載された方法を利用することができる。特開２００６−１０３６４０号公報の方法では、印刷媒体にテストインク量Ｉ_jのインクを付着させて形成したカラーパッチを測定することにより得られた粒状性指数ＧＩに基づいて学習したニューラルネットワークによって任意のインク量Ｉ_jで印刷した場合の粒状性指数ＧＩを予測する。本実施例では、ニューラルネットワークが基本メディアに形成したカラーパッチの測定結果に基づいて学習されている。実体的には、画質評価指数コンバーター１３６がニューラルネットワークに任意のインク量Ｉ_jを入力することにより、該インク量Ｉ_jのインクを基本メディアに付着させた場合の粒状性指数ＧＩを算出する。

色非恒常性指数ＣＩＩは、例えば以下の（１７）式で与えられる。

ここで、ΔＬ^*は２つの異なる観察条件下（異なる光源下）におけるカラーパッチの明度差、ΔＣ^* _abは彩度差、ΔＨ^* _abは色相差を示す。色非恒常性指数ＣＩＩの計算時には、２つの異なる観察条件下でのＬ^*ａ^*ｂ^*値は、色順応変換（ＣＡＴ）を用いて標準観察条件（例えば標準の光Ｄ６５の観察下）に変換される。なお、観察条件下でのＬ^*ａ^*ｂ^*値は、上述したフォワードモデルコンバーター３００によって算出される。フォワードモデルコンバーター３００（分光プリンティングモデルコンバーター３１０）は基本メディアについて準備されたものであるため、色非恒常性指数ＣＩＩによれば各インク量Ｉ_jのインクを基本メディアに付着させた場合の色非恒常性を評価することができる。ＣＩＩについては、Billmeyer and Saltzman's Principles of Color Technology， 3rd edition， John Wiley & Sons， Inc， 2000， p.129， pp. 213-215を参照。

ガマット評価指数ＧＭＩは、フォワードモデルコンバーター３００によって得られる色彩値Ｌ^*ａ^*ｂ^*と目標色彩値Ｌ_GM ^*ａ_GM ^*ｂ_GM ^*との色差ΔＥ（ＣＩＥ１９７６）で与えられる。目標色彩値Ｌ_GM ^*ａ_GM ^*ｂ_GM ^*は、Ｌ^*ａ^*ｂ^*色空間の最外縁の色彩値とされる。ガマット評価指数ＧＭＩはすべての格子点について考慮されるのではなく、ガマットの頂点や稜線や外面上の格子点についてのみ考慮すればよい。また、各格子点に応じて目標色彩値Ｌ_GM ^*ａ_GM ^*ｂ_GM ^*は異なる。例えば、フォワードモデルコンバーター３００によって得られる色彩値Ｌ^*ａ^*ｂ^*と色相角が同じで、かつ、より高彩度（Ｌ^*ａ^*ｂ^*色空間の最外縁）の色彩値を目標色彩値Ｌ_GM ^*ａ_GM ^*ｂ_GM ^*とすれば高彩度側にガマットが広いか否かを評価することができる。なお、グレー軸格子点に対応する格子点についての目標色彩値Ｌ_GM ^*ａ_GM ^*ｂ_GM ^*を、グレーターゲットの色彩値とすることによっても、該格子点をグレーターゲットに拘束することもできる。

ヤコビ行列Ｊの複数の成分（「要素」とも呼ぶ）のうち、例えばＬ^*値に関する成分は、（１８）式で与えられる。

ここで、ｆ_L*FMは、フォワードモデルによるインク量ＩからＬ^*値への変換関数、Ｉ_rはインク量Ｉの現在値（平滑化／最適化処理前のインク量）、ｈ_jはｊ番目のインク量Ｉ_jの微小変動量である。Ｌ^*値について（１７）式を例示したが、ａ^*ｂ^*値についても同様である。Ｌ^*ａ^*ｂ^*値は、上述したフォワードモデルコンバーター３００（（１１）式）によって算出されるため、Ｌ^*ａ^*ｂ^*値は各インク量Ｉ_jのインクを基本メディアに付着させた場合の色彩値を意味する。前記（１４）、（１７）式に準じて、ヤコビ行列Ｊの最下行の要素を算出すると、ヤコビ行列Ｊの最下行の要素はすべて１となる。あるインクのインク量Ｉ_jが微小変動量ｈ_jだけ変動した場合の合計インク量ＴＩの変動量もｈ_jとなるからである。

最適化の目的関数Ｅは、例えば以下の（１９）式で与えられる。

ここで、右辺の各項の最初に記載されているｗ_L*，ｗ_a*等は、各項の重みである。各項の重みｗ_L*、ｗ_a*…は、ステップＳ０５において、ユーザーから指定され、設定テーブルＳＴＢに格納された重みｗ_L*，ｗ_a*…が使用される。特に、ユーザーがポインターの位置を初期位置から移動させなかった場合には、デフォルト重みｗ_L*，ｗ_a*…が使用される。従って、目的関数Ｅによって重要視される項目がメディアおよびユーザーの設定に依存することとなる。

（１９）式の右辺第１項ｗ_Ｌ*（ΔＬ^*−ΔＬ^* _t）²は、色彩値Ｌ^*の変動量ΔＬ^*，ΔＬ^* _tに関する２乗誤差である。これらの変動量ΔＬ^*，ΔＬ^* _tは、次の式で与えられる。

前記（２０）式の右辺における偏微分値はヤコビ行列（（１０）式）で与えられる値であり、Ｉ_jは最適化処理の結果として得られるインク量であり、Ｉ_jrは現在のインク量である。第１の変動量ΔＬ^*は、最適化処理によるインク量の変動量ΔＩ_jを、ヤコビ行列の成分である偏微分値で線形変換した量である。一方、第２の変動量ΔＬ^* _tは、ステップＴ１２０の平滑化処理で得られたターゲット値Ｌ^* _tと、現在インク量Ｉ_rで与えられる色彩値Ｌ^*（Ｉ_r）との差分である。なお、第２の変動量ΔＬ^* _tは、平滑化処理の前後におけるＬ^*値の差分と考えることが可能である。

前記（１９）式の右辺の第２項以降の各項も、前記（２０）式および（２１）式と同様の式で与えられる。すなわち、目的関数Ｅは、最適化処理によるインク量の変動量ΔＩ_jをヤコビ行列の成分で線形変換して得られる第１の変動量ΔＬ^*，Δａ^*，Δｂ^*，ΔＧＩ…と、パラメータＬ^*，ａ^*，ｂ^*，ＧＩ…に関する平滑化処理の前後における第２の変動量ΔＬ^* _t，Δａ^* _t，Δｂ^* _t，ΔＧＩ_t…と、の２乗誤差の和として与えられている。

ところで、第１の変動量ΔＬ^*，Δａ^*，Δｂ^*，ΔＧＩ…は、行列を用いて以下の（２２）式および（２３）式の形式に書き表すことが可能である。

前記（１９）式は、行列を用いて（２４）式のように表記できる。

ここで、Ｔは行列の転置を表している。行列Ｗ_Mはそれぞれ対角要素に重みを配置した対角行列（（２５）式参照）であり、行列ΔＭは各パラメータに関する目標変動量ベクトル（（２６）式参照）である。

（２６）式の右辺は、各パラメータＬ^*，ａ^*，ｂ^*，ＣＩＩ…（「要素」とも呼ぶ）に関するターゲット値と、現在のインク量Ｉ_rで与えられる各パラメータ値との差分である。各パラメータのターゲット値のうち、色彩値Ｌ^* _t，ａ^* _t，ｂ^* _tは平滑化処理（ステップＴ１２０）で決定される。画質評価指数のターゲット値と現在の画質評価指数から求められる目標変動量ΔＧＩ_t，ΔＣＩＩ_t，ΔＴＩ_t，ΔＧＭＩ_tについては、いくつかの決定方法がある。第１の方法は、目標変動量ΔＧＩ_t，ΔＣＩＩ_t，ΔＴＩ_t，ΔＧＭＩ_tとして所定の定数（例えばΔＧＩ_t＝−２，ΔＣＩＩ_t＝−１，ΔＴＩ_t＝−１，ΔＧＭＩ_t＝−１）を使用する方法である。なお、定数としてマイナスの値を使用する理由は、これらの画質評価指数は、より小さいほど高画質であることを示す指数だからである。また、粒状性指数ＧＩのターゲット値ＧＩ_tは、ゼロとすることも好ましい。第２の方法は、ターゲット値ＧＩ_t，ＣＩＩ_t，ＴＩ_t，ΔＧＭＩ_tを色彩値のターゲット値Ｌ^* _t，ａ^* _t，ｂ^* _tの関数として定義しておく方法である。以上のように、各パラメータのターゲット値は最適化処理前に決められているので、目標変動量ベクトルΔＭの各成分はすべて定数である。

前記（２４）式の右辺の各項のうち、第３項（Ｉ_r ^TＪ^T+ΔＭ^T）Ｗ_M（ＪＩ_r+ΔＭ）、は、最適化処理の結果として得られるインク量Ｉを含まないので定数である。一般に、最適化のための目的関数Ｅにおいて定数項は不要である。そこで、前記（２４）式から定数項を削除して全体に１／２を乗じると、次の（２７）式が得られる。

ここで、以下の（２８）式および（２９）式のように行列Ａおよびベクトルｇを定義すると、前記（２７）式は（３０）式のように書き表せる。

（３０）式で与えられる目的関数Ｅは、最適化で得られるインク量ベクトルＩに関する２次形式であることが理解できる。図１５（Ｃ）に示した（ＥＱ１）式と（ＥＱ２）式は、（１８）式と（２９）式とそれぞれ同じものである。

本実施例の最適化処理では、（３０）式のような２次形式の目的関数Ｅを用いるので、最適化手法として２次計画法を使用することが可能である。ここで、「２次計画法」とは、逐次２次計画法を含まない狭義の２次計画法を意味している。２次形式の目的関数を用いた２次計画法を利用すれば、準ニュートン法や逐次２次計画法などの他の非線形計画法を利用する場合に比べて、処理を大幅に高速化することが可能である。

ところで、本実施例における最適化処理によるインク量の探索は、以下の条件の下で実行される。
（最適化条件）目的関数Ｅを最小とする。
（制約条件）デューティ制限値を守る。

基本メディアの場合、デューティ制限値として、ステップＳ０８において設定テーブルＳＴＢに格納したデューティ制限値Ｄ_Ijがそのまま使用される。なお、画質評価指数コンバーター１３６とフォワードモデルコンバーター３００（分光プリンティングモデルコンバーター３１０）は、デューティ制限値Ｄ_Ijを満足するインク量Ｉ_jについて色彩値や画質評価指数ＧＩ等を予測することができる。

デューティ制限値に関する制約条件は、次の（３１）式で表すことができる。

ここで、ベクトルｂは、デューティ制限値の対象となるインク種類を識別するための係数であり、要素に０か１を持つベクトルである。例えば、１種類のインクに関するデューティ制限値の場合には、ベクトルｂの１個の要素のみが１となる。一方、全インクの合計インク量に関するデューティ制限値の場合には、ベクトルｂのすべての要素が１となる。（３１）式の右辺のＤ_Iは、個々のデューティ制限値Ｄ_Ijを要素とするベクトルである。（３１）式の右辺、左辺とも、ｊ＝１〜８であるとする。すなわち、デューティ制限値に関する制約条件を課す際には、２次色の合計インク量Ｉ₅〜Ｉ₇と全部の合計インク量Ｉ₈も考慮する。

各インク量Ｉ_j（ｊ＝１〜８）には、負でないという制約も存在する。この非負制限は、以下の（３２）式で表せる。

前記（３１）式と（３２）式とを合体すると、デューティ制限値は、次の（３３）式で与えられる。

この（３３）式で表される制約は、線形不等号制約である。一般に、２次計画法は線形制約の下で実行することが可能である。すなわち、本実施例における最適化処理では、（３３）式の制約の下で、前記（３０）式で与えられる２次形式の目的関数Ｅを用いた２次計画法を実行することによって、最適なインク量を探索する。この結果、この線形制約を厳密に満足しつつ、インク量探索を高速に実行することが可能である。

図１６は、最適化処理（図８のステップＴ１３０）の詳細手順を示すフローチャートである。ステップＴ１３２では、まず、前記（２６）式で与えられる目標変動量ΔＭを求める。この目標変動量ΔＭは、前述したように、ステップＴ１２０（平滑化処理）で得られたターゲット値Ｌ^* _t，ａ^* _t，ｂ^* _tと現在インク量Ｉ_rに基づいて決定される。

ステップＴ１３４では、前記（１０）式で与えられるヤコビ行列Ｊを算出する。なお、ヤコビ行列Ｊの各成分は、前記（１８）式で例示されるように、インク量の現在値Ｉ_r（平滑化／最適化前の値）に関して算出される値である。

ステップＴ１３６では、ヤコビ行列Ｊによる線形変換の結果ΔＬ^*，Δａ^*，Δｂ^*，ΔＧＩ…と、目標変動量ΔＭ（ΔＬ^* _t，Δａ^* _t，Δｂ^* _t，ΔＧＩ_t…）との差が最小になるように、インク量Ｉ_jの最適化を実行する。この最適化は、前記（３０）式で与えられる２次形式の目的関数Ｅを用いた２次計画法を実行することによって実現される。上述したようにステップＳ０８において設定テーブルＳＴＢに格納したデューティ制限値Ｄ_Ijによる制約の下でインク量Ｉ_jが最適化されるため、メディアに応じて異なる最適化の結果が得られることとなる。特にガマットの大きさは、デューティ制限値Ｄ_Ijに大きく依存し、各メディア間で異なったものとなる。むろん、画質評価指数コンバーター１３６とフォワードモデルコンバーター３００による色彩値や画質評価指数ＧＩ等の予測結果も各メディアに応じて異なるため、異なるインク量Ｉ_jの最適化の結果となる。

なお、図１４のフローチャートにおいて既に説明したように、ステップＴ１３０の最適化処理の後、収束が不十分と判断される場合には、平滑化処理（ステップＴ１２０）および最適化処理（ステップＴ１３０）が再度実行される。この際、平滑化／最適化処理の初期値としては、その前の平滑化／最適化処理で得られた値が利用される。なお、このような繰り返し処理は必須ではなく、少なくとも１回の平滑化／最適化処理を行えばよい。

Ｃ．印刷装置の構成：
図１７は、プリンター１０の構成を示している。同図において、プリンター１０はＣＰＵ１３とＲＡＭ１４とＲＯＭ１５とメモリーカードスロット２０とバス２１とＡＳＩＣ２２を備えている。ＲＯＭ１５に記憶されたプログラムデータ１５ａをＲＡＭ１４に展開しつつＣＰＵ１３がプログラムデータ１５ａにしたがった演算を行うことによりプリンター１０を制御するためのファームウェアＦＷが実行される。ファームウェアＦＷは、メモリーカードスロット２０に装着されたメモリーカードＭＣに記憶された印刷データＰＤに基づいて駆動データを生成する。ＡＳＩＣ２２は駆動データを取得し、紙送り機構１７やキャリッジモータ１８や印刷ヘッド１９の駆動信号を生成する。ＲＯＭ１５においては、色補正３Ｄ−ＬＵＴ６１０が記憶されている。色補正３Ｄ−ＬＵＴ６１０には設定テーブルＳＴＢが添付されており、設定テーブルＳＴＢに記述されたメディアごとに色補正３Ｄ−ＬＵＴ６１０が用意されている。プリンター１０はキャリッジ１１を備えており、キャリッジ１１はインクカートリッジ１２，１２・・・が取り付け可能なカートリッジホルダー１１ａ，１１ａ・・・を備えている。そして、各インクカートリッジ１２，１２・・・をカートリッジホルダー１１ａ，１１ａ・・・に挿入して固定することができる。

本実施形態においてカートリッジホルダー１１ａ，１１ａ・・・は４個備えられているものとし、カートリッジホルダー１１ａ，１１ａ・・・にシアン（Ｃ）、マゼンタ（Ｍ）、イエロー（Ｙ）、ブラック（Ｋ）のインクをそれぞれ充填した４種類のインクのインクカートリッジ１２，１２・・・が取り付けられる。キャリッジ１１は、インクカートリッジ１２，１２・・・から供給されるインクを多数のインクノズルＮｚから吐出する印刷ヘッド１９を備えている。

図１８は、本実施例のプリンター１０の印刷方式を模式的に示している。同図において、プリンター１０は、ＣＭＹＫのインクごとに複数のノズルＮｚ・・・を備えた印刷ヘッド１９を備えており、ノズルＮｚ・・・が吐出するインク量をＣＭＹＫごとに制御する処理が駆動データに基づいて行われる。各ノズルＮｚ・・・が吐出したインク滴は印刷用紙上において微細なドットとなり、多数のドットの集まりによって各インク量Ｉ_jに応じたインク被覆率の印刷画像が印刷用紙上に形成されることとなる。本明細書において、各インク量Ｉ_jに基づいて印刷することは、各インク量Ｉ_jに基づく量のインクがメディアに付着し、各インク量Ｉ_jに応じたインク被覆率の画像がメディア上に再現されることを意味する。本実施形態におけるプリンター１０はインクジェット方式のプリンターであるが、インクジェット方式の他にも種々のプリンターに対して本発明を適用可能である。

図１９は、ファームウェアＦＷのソフトウェア構成を示している。ファームウェアＦＷは、画像データ取得部ＦＷ１とレンダリング部ＦＷ２と色変換部ＦＷ３とハーフトーン部ＦＷ４とラスタライズ部ＦＷ５とから構成されている。画像データ取得部ＦＷ１は、メモリーカードＭＣに記憶された印刷データＰＤを印刷対象として取得する。印刷データＰＤは、文書データやグラッフィックデータであってもよいし、写真画像データであってもよい。レンダリング部ＦＷ２は、印刷データＰＤに基づいて印刷に使用する入力画像データＩＤを生成する。入力画像データＩＤは、印刷解像度（例えば２８８０×２８８０ｄｐｉ）に対応した画素数（印刷解像度×印刷実サイズ）の画素で構成されており、各画素が８ビット（０〜２５５）のＣＩＥ−ｓＲＧＢ色空間に準拠したＲＧＢ値で表現されている。

色変換部ＦＷ３は、入力画像データＩＤを取得し、該入力画像データＩＤを色変換する。具体的には、色変換部ＦＷ３は、色補正３Ｄ−ＬＵＴ６１０を参照しつつ補間演算を実行することにより、ＲＧＢ値を各インクのインク量Ｉ_jに変換する。メディアごとに色補正３Ｄ−ＬＵＴ６１０が用意されているが、対応する色補正３Ｄ−ＬＵＴ６１０が記憶されていないメディアがプリンター１０にセットされたり、指定されたりした場合には、該メディアに印刷を実行することができない。ハーフトーン部ＦＷ４は、色変換部ＦＷ３が出力した各インクのインク量Ｉ_jに基づくハーフトーン処理を実行する。ラスタライズ部ＦＷ５は、ハーフトーン処理後のハーフトーンデータの各画素（吐出可否）を印刷ヘッド１９の各主走査および各インクノズルＮｚに割り当て、駆動データを生成する。駆動データはＡＳＩＣ２２に出力され、ＡＳＩＣ２２が紙送り機構１７やキャリッジモータ１８や印刷ヘッド１９の駆動信号を生成する。本実施例では、プリンター１０上のファームウェアＦＷによって色変換処理が行われることとしたが、プリンター１０が接続されたコンピューター上で色変換処理が行われてもよい。すなわち、色補正３Ｄ−ＬＵＴ６１０はプリンターだけでなくプリンターを制御するコンピューター（印刷制御装置）に組み込まれてもよい。

上述したように、対応する色補正３Ｄ−ＬＵＴ６１０が記憶されていないメディアには印刷を実行することができない。色補正３Ｄ−ＬＵＴ６１０を作成する際に使用したフォワードモデルコンバーター３００や画質評価指数コンバーター１３６が、あるメディア（基本メディア）を対象として準備されたものであり、他のメディアについて色補正３Ｄ−ＬＵＴ６１０を使用しても所望の再現色や画質が実現できないからである。そのため、理想的には、あらゆる印刷媒体について色補正３Ｄ−ＬＵＴ６１０を用意しておく必要がある。すなわち、あらゆる印刷媒体について、上述した手順によって、色補正３Ｄ−ＬＵＴ６１０を作成する処理を行っておく必要がある。しかしながら、フォワードモデルコンバーター３００や画質評価指数コンバーター１３６を準備するためには各メディアに多数のカラーパッチを形成し、測定を行う必要があるため、あらゆるメディアについてフォワードモデルコンバーター３００や画質評価指数コンバーター１３６を準備するのは困難である。また、未知のメディアについてフォワードモデルコンバーター３００や画質評価指数コンバーター１３６を準備しておくのは、そもそも不可能である。そのため、本発明では、フォワードモデルコンバーター３００や画質評価指数コンバーター１３６をメディア毎に準備することなく、基本メディアとは異なる転用メディアについてのＬＵＴを作成する。以下、その手順について説明する。

Ｄ．転用メディアのＬＵＴ作成手順：
図２の全体処理のステップＳ０６において、転用メディアが指定されたと判定した場合、ステップＳ０９〜を実行することとなる。まず、ステップＳ０９において、ＬＵＴ作成条件設定モジュール７００は、被転用メディアを決定する。被転用メディアは、指定された転用メディアの系統と同じ系統の基本メディアとすることとする。指定された転用メディアの系統が未分類であった場合には、基本メディアのうち各インクの発色特性が標準的なものを被転用メディアとする。被転用メディアを特定する情報は、設定情報格納モジュール７３０によって設定テーブルＳＴＢに登録される。ステップＳ１０においては、設定情報格納モジュール７３０が、メディアテーブルＭＴＢに記憶された被転用メディアの発色特性データとデューティ制限値Ｄ_Ij（以下、基準デューティ制限値Ｄ_SIjと表記する。）とを取得する。

ステップＳ１１においては、ＬＵＴ作成条件設定モジュール７００がメディアの発色特性評価用のカラーパッチを指定された転用メディアに印刷させる。具体的には、ＬＵＴ作成条件設定モジュール７００がカラーパッチを印刷させるためのパッチ画像データを生成し、該パッチ画像データをプリンター１０のハーフトーン部ＦＷ４に出力する。パッチ画像データは、各画素がインク量Ｉ_j（ｊ＝１〜４）を有する画像データであり、各インクの１次色（Ｃ，Ｍ，Ｙ，Ｋ）と２次色（Ｒ，Ｇ，Ｂ）と全インクを混色したインク量Ｉ_j（ｊ＝１〜８）のグラデーションに基づくカラーパッチを印刷させるものである。例えば、２次色と全インクのグラデーションにおける個々のインク量Ｉ_j（ｊ＝１〜４）は、均等とする。各カラーパッチには、該カラーパッチの印刷につき使用されたインク量Ｉ_j（ｊ＝１〜８）を表す文字が付記されることとする。ステップＳ１２においては、印刷された各カラーパッチを測色し、その測色値（ＣＩＥ−Ｌ^*ａ^*ｂ^*表色系）を取得する。ステップＳ１３においては、ＵＩモジュール７２０が転用メディアの発色特性を指定するためのメディア特性指定ＵＩ画像を表示装置に表示させる。

図２０は、メディア特性指定ＵＩ画像を示す図である。同図において、各グラデーションのカラーパッチ（Ｃ，Ｍ，Ｙ，Ｋのみ）を測色して得られた測色値がグラフによって示されている。有彩色を示すカラーパッチ（Ｃ，Ｍ，Ｙ）についてのグラフでは縦軸が彩度Ｃ^*を表している。無彩色を示すカラーパッチ（Ｋ）についてのグラフでは縦軸が明度Ｌ^*を表している。横軸は、インク量Ｉ_j（ｊ＝１〜８）を表している。各グラフでは、各カラーパッチの測色値に基づく彩度Ｃ^*と明度Ｌ^*とがプロット（白丸）されている。これにより、転用メディアにおけるインク量Ｉ_j（ｊ＝１〜８）に応じた発色特性がグラフに表されたこととなる。また、各グラフにおいては、メディアテーブルＭＴＢから取得した被転用メディアについての発色特性も対比可能にプロット（黒丸）されている。被転用メディア（基本メディア）についても、同様のカラーパッチを印刷・測色されており、その測色結果が予め発色特性データとしてメディアテーブルＭＴＢに格納されている。

また、メディア特性指定ＵＩ画像において、各インクの１次色（Ｃ，Ｍ，Ｙ，Ｋ）と２次色（Ｒ，Ｇ，Ｂ）と全インクのカラーパッチをユーザーが観察した結果、にじみが生じ始めたカラーパッチに付記されたインク量Ｉ_j（ｊ＝１〜８）をそれぞれデューティ制限値Ｄ_Ijとして入力するためのテキストボックスが設けられている。ユーザーの観察によらず、にじみ等により測色値が不安定となり始めたインク量Ｉ_j（ｊ＝１〜８）をＬＵＴ作成条件設定モジュール７００が判定し、デューティ制限値Ｄ_Ijとしてもよい。メディア特性指定ＵＩ画像には、デューティ制限値確定ボタンが設けられており、該ボタンをクリックされたことを受け、設定情報格納モジュール７３０が入力されたデューティ制限値Ｄ_Ijを取得し（ステップＳ１４）、該デューティ制限値Ｄ_Ijを設定テーブルＳＴＢに格納する（ステップＳ１５）。以上により、デューティ制限値Ｄ_Ijが未知の転用メディアについてもデューティ制限値Ｄ_Ijが設定できたこととなる。

デューティ制限値Ｄ_Ijが設定できると、ＬＵＴ作成条件設定モジュール７００は、下記の（３４）式によって、被転用メディアについての各基準デューティ制限値Ｄ_SIjと、転用メディアについての各デューティ制限値Ｄ_Ijとに基づいて正規化比ＲＷを算出する（ステップＳ１６）。

（３４）式が示すように、正規化比ＲＷは、各基準デューティ制限値Ｄ_SIjを各デューティ制限値Ｄ_Ijによって除算した制限値比のうち最も値の小さい成分とされる。正規化比ＲＷが算出できると、インク量コンバーター７１０は、転用メディアに印刷した各カラーパッチの測色値のプロット（白丸）の横軸方向の位置を示すインク量Ｉ_j（ｊ＝１〜４）に正規化比ＲＷを乗算することにより、インク量Ｉ_j（ｊ＝１〜４）を第１仮インク量Ｉ_PPj（ｊ＝１〜４）へと変換する（ステップＳ１７）。すなわち、第１仮インク量Ｉ_PPjとインク量Ｉ_jとの関係は（３５）式で表される。

第１仮インク量Ｉ_PPj（ｊ＝１〜４）の取り得る範囲は、０〜Ｄ_PIjとなる。なお、２次色以上のインク量Ｉ_j（ｊ＝５〜８）は、インク量Ｉ_j（ｊ＝１〜４）の変換に追従して変換されるため、変換の対象としていない。

そして、転用メディアに印刷した各カラーパッチの測色値のプロット（白丸）の横軸方向の位置を第１仮インク量Ｉ_PPj（ｊ＝１〜４）が示す位置に更新する（ステップＳ１８）。このとき、被転用メディアについての発色特性を示すプロット（黒丸）はそのままとしておく。このようにすることにより、転用メディアと被転用メディアの発色特性が見かけ上近似することとなる。本明細書において発色特性とは、インク量Ｉ_ｊのインクを付着させたときに各メディア上において再現される色のＬ^*，Ｃ^*値の推移を意味する。インク量Ｉ_j（ｊ＝１〜４）から第１仮インク量Ｉ_PPj（ｊ＝１〜４）へと変換は、インク量Ｉ_j（ｊ＝１〜４）の取り得る範囲を正規化する意味合いを有し、それにより、転用メディアと被転用メディアの発色特性がある程度近似することが期待できる。次に、インク量コンバーター７１０は、転用メディアの発色特性を、さらに被転用メディアの発色特性に近似させるよう、第１仮インク量Ｉ_PPj（ｊ＝１〜４）を第２仮インク量Ｉ_SPj（ｊ＝１〜４）へと変換する非線形変換関数を設定する（ステップＳ１９）。

図２１（Ａ）〜２１（Ｃ）は、変換関数を設定する様子を示すグラフである。同図においては、Ｋインクの発色特性を示すグラフを示している。図２１（Ａ）は、いずれの変換も行っていない状態の発色特性を示している。図２１（Ｂ）は、転用メディア（実線）のインク量を正規化比ＲＷ（Ｋインクの制限値比が最小であるとする）に基づいて第１仮インク量Ｉ_PPj（ｊ＝１〜４）へ変換した状態を示している。図２１（Ｃ）は、転用メディアのインク量を非線形変換により第２仮インク量Ｉ_PPj（ｊ＝１〜４）へ変換した状態を示している。本実施例では、Ｌ^*＝５０の発色に注目し、転用メディア（実線）においてＬ^*＝５０を再現する第１仮インク量Ｉ_PP4が、被転用メディア（破線）においてＬ^*＝５０を再現するインク量Ｉ₄と等しい値（第２仮インク量Ｉ_SP4）に変換されるようなガンマ関数（γ_j値）を算出する。第２仮インク量Ｉ_SPjは、以下の（３６）式によって与えられる。

γ_j値は、各インクごと（ｊ＝１〜４）に異なる値が求められる。本実施例では、Ｌ^*＝５０の発色に注目したが、最小二乗法等によって全体的に発色特性を近似させるようなγ_j値を算出してもよい。また、グラフに示した発色特性にユーザーがドラッグ＆ドロップ可能な単一または複数の制御点を設け、該ドラッグ＆ドロップされた制御点を通過するようなγ_j値を算出するようにしてもよい。むろん、ガンマ関数以外の非線形関数（スプライン曲線，ベジエ曲線，指数関数等）によって近似してもよい。

図２１（Ａ）〜２１（Ｃ）では、Ｋインクについて例示したが、他の有彩色インクについても、転用メディアにおいて所定のＣ^＊値を再現する第１仮インク量Ｉ_PPjが、被転用メディアにおいて該Ｃ^＊値を再現するインク量Ｉ_jと等しい値（第２仮インク量Ｉ_SPj）に変換されるようなガンマ関数を算出する。なお、前記（３５）式により、Ｉ_PPj＝ＲＷ・Ｉ_j，Ｄ_PIj＝ＲＷ・Ｄ_Ijであるため、前記の（３６）式は、下記の（３７）式に書き換えることができる。

以上により、転用メディアに各インクを付着させた場合の発色特性が、被転用各インクを付着させた場合の発色特性に近似するようにインク量Ｉ_jを変換する変換式が得られたこととなる。ｆ_CVjは、インク量Ｉ_j（ｊ＝１〜４）を第２仮インク量Ｉ_SPjに変換する変換関数を意味する。ここで、前記の（３７）式を満足すれば、第２仮インク量Ｉ_SPjのインクを被転用メディアに付着した場合の発色（Ｌ^*，Ｃ^*値）は、インク量Ｉ_jのインクを転用メディアに付着した場合の発色（Ｌ^*，Ｃ^*値）とほぼ等しくなるということが言える。

ステップＳ２０においては、（３７）式の変換関数ｆ_CVjをインク量コンバーター７１０に設定し、変換関数ｆ_CVj（γ値）を設定テーブルＳＴＢに登録する。これにより、インク量コンバーター７１０は、インク量Ｉ_j（ｊ＝１〜４）を第２仮インク量Ｉ_SPj（ｊ＝１〜４）に変換することが可能となる。また、インク量コンバーター７１０は、変換関数ｆ_CVjの逆変換関数ｆ_CVj ^-1を使用して第２仮インク量Ｉ_SPjをインク量Ｉ_j（ｊ＝１〜４）に逆変換することも可能となる。前記の（３７）式のインク量Ｉ_jに転用メディアのデューティ制限値Ｄ_Ijを代入すると、第２仮インク量Ｉ_SPj＝ＲＷ・Ｄ_Ijとなる。これを前記の（３５）式と比較すると、デューティ制限値Ｄ_Ijについては、Ｉ_PPj＝Ｉ_SPj＝ＲＷ・Ｄ_Ijとなり、第１仮インク量Ｉ_PPjと第２仮インク量Ｉ_SPjとが等しくなる。なお、インク量コンバーター７１０は、複数のインク量Ｉ_jと第２仮インク量Ｉ_SPjとの対応関係を記述したＬＵＴを参照して、ｆ_CVj，ｆ_CVj ^-1に相当する変換・逆変換を行うようにしてもよい。

ステップＳ２１においては、転用メディアに対してＣＭＹインクをそれぞれ単独でデューティ制限値Ｄ_Ij（ｊ＝１〜３）まで付着させたときのａ^*，ｂ^*値（以下、転用メディア色味（ａ_Cj ^*，ｂ_Cj ^*）（ｊ＝１〜３）と表記する。）、および、被転用メディアに対してＣＭＹインクをそれぞれ単独で基準デューティ制限値Ｄ_SIj（ｊ＝１〜３）まで付着させたときのａ^*，ｂ^*値（以下、被転用メディア色味（ａ_Sj ^*，ｂ_Sj ^*）（ｊ＝１〜３）と表記する。）を取得し、これらを解析する。これらの色彩値は、ステップＳ１０，Ｓ１２において、デューティ制限値Ｄ_Ij（ｊ＝１〜３）と基準デューティ制限値Ｄ_SIj（ｊ＝１〜３）に対応するカラーパッチから得られた測色値として取得されている。なお、転用メディア色味（ａ_Cj ^*，ｂ_Cj ^*）と、被転用メディア色味（ａ_Sj ^*，ｂ_Sj ^*）は、所定の明度Ｌ^*や彩度Ｃ^*が再現されるカラーパッチから取得してもよいし、デューティ制限値Ｄ_Ij，Ｄ_SIj以外の所定のインク量Ｉ_j（例えば、デューティ制限値Ｄ_Ij，Ｄ_SIjに対して１５％等。）によるカラーパッチから取得してもよい。さらに、インク量Ｉ_j（ｊ＝１〜３）を等量ずつ各メディアに付着させて形成したコンポジットグレーのカラーパッチから、転用メディア色味（ａ_Cj ^*，ｂ_Cj ^*）と、被転用メディア色味（ａ_Sj ^*，ｂ_Sj ^*）とを取得してもよい。いずれの場合も、各メディアそのものの色彩値を取得するよちも、ある程度インクが付着した状態の各メディアの色彩値が取得されるため、中明度領域の色味を評価することができる。

図２２は、転用メディア色味（ａ_Cj ^*，ｂ_Cj ^*）（ｊ＝１〜３）（白丸）と、被転用メディア色味（ａ_Sj ^*，ｂ_Sj ^*）（ｊ＝１〜３）（黒丸）をａ^*ｂ^*平面にプロットしたグラフである。転用メディアと被転用メディアは、互いに異なる色味を有しているため、転用メディア色味（ａ_Cj ^*，ｂ_Cj ^*）と、被転用メディア色味（ａ_Sj ^*，ｂ_Sj ^*）は完全には一致しない。ステップＳ２０では、（３８）式のように、転用メディア色味（ａ_Cj ^*，ｂ_Cj ^*）（ｊ＝１〜３）を、原点ａ^*，ｂ^*＝０からの位置ベクトルとして、互いに足し合わせることにより転用メディア色味（ａ_C ^*，ｂ_C ^*）（白三角）を算出する。同様に、被転用メディア色味（ａ_Sj ^*，ｂ_Sj ^*）を互いに足し合わせることにより、被転用メディア色味（ａ_S ^*，ｂ_S ^*）（黒三角）を算出する。

さらに、（３９）式のように、転用メディア色味（ａ_C ^*，ｂ_C ^*）（白三角）から被転用メディア色味（ａ_S ^*，ｂ_S ^*）（黒三角）を差し引くことにより、差分色味（ａ_D ^*，ｂ_D ^*）を算出する。

以上のようにして、差分色味ａ_D ^*，ｂ_D ^*を算出できると、差分色味（ａ_D ^*，ｂ_D ^*）の符号を逆にしたベクトル（−ａ_D ^*，−ｂ_D ^*）が示す色相方向を、グレーターゲットの色味（ａ_gt ^*，ｂ_gt ^*）の色相方向として設定する（ステップＳ２２）。すなわち、グレーターゲットの色味（ａ_gt ^*，ｂ_gt ^*）は、ベクトル（−ａ_D ^*，−ｂ_D ^*）に正の係数ｋを乗算したベクトルとなる。係数ｋの大きさは、例えばユーザーによって設定される。ステップＳ２３では、設定情報格納モジュール７３０が、グレーターゲットの色味（ａ_gt ^*，ｂ_gt ^*）を設定テーブルＳＴＢに格納する。さらに、ステップＳ２４において、設定情報格納モジュール７３０は、インク量コンバーター７１０を有効とする旨の有効フラグを設定テーブルＳＴＢに格納させる。以上の処理により、図４に示す設定テーブルＳＴＢには、必要な設定情報が記憶されたこととなる。以降、設定テーブルＳＴＢを参照しながら、転用メディアのためのＬＵＴを作成する処理に移行する。ここでは、基本メディアのためのＬＵＴを作成する処理と異なる部分について順に説明していく。

まず、図６のステップＳ１００において、設定テーブルＳＴＢに記憶された情報に基づいて、各コンバーター等３００，３１０，４１０，１３６を準備（起動）することとなるが、転用メディアの場合には設定テーブルＳＴＢに記憶された被転用メディアについての各コンバーター等３００，３１０，４１０，１３６が準備（起動）されることとなる。すなわち、転用メディアについては各コンバーター等３００，３１０，４１０，１３６が準備されていないため、基本メディアである被転用メディアについての各コンバーター等３００，３１０，４１０，１３６を転用することとする。また、有効フラグが添付されているため、インク量コンバーター７１０が起動し、上述した変換関数ｆ_CVjと逆関数ｆ_CVj ^-1を使用したインク量Ｉ_jと第２仮インク量Ｉ_SPjとの相互変換が有効となる。

図１３のステップＴ１３０におけるスムージング処理（ステップＴ１０７）では、図２のステップＳ２３にて設定テーブルＳＴＢに格納されたグレーターゲットの色味（ａ_gt ^*，ｂ_gt ^*）＝ｋ×（−ａ_D ^*，−ｂ_D ^*）が設定されることとなる。すなわち、基本メディアのＬＵＴを作成する場合は常にＬ^*軸上にグレーターゲットが設定されていたのに対して、転用メディアのＬＵＴを作成する場合は、グレーターゲットが転用メディアの相対的な色味の逆色相方向に遷移することとなる。

図２３は、転用メディアのＬＵＴを作成する場合のグレーターゲットを示す図である。同図に示すように、グレー軸格子点に対応する格子点（白丸）が、Ｌ^*ａ^*ｂ^*色空間のＬ^*軸上の紙黒点と紙白点とを結び、かつ、グレーターゲットの色味（ａ_gt ^*，ｂ_gt ^*）の色相方向に湾曲した線分を１６等分する位置を示すように拘束されている。これにより、スムージング処理後のグレー軸格子点に対応する格子点の位置は、Ｌ^*軸からずれた位置となる。湾曲する量は中明度領域において最大となり、本実施例ではＬ^*＝５０における最大値がグレーターゲットの色味（ａ_gt ^*，ｂ_gt ^*）となる。また、グレーターゲットにユーザーがドラッグ＆ドロップ可能な単一または複数の制御点を設け、該ドラッグ＆ドロップされた制御点を通過するような湾曲する量を規定する係数ｋを決定するようにしてもよい。

さらに、図１３のステップＴ１３０における最適化処理では、設定テーブルＳＴＢに記憶された転用メディアについての各インクのデューティ制限値Ｄ_Ij（ｊ＝１〜８）がインク量コンバーター７１０に設定された変換関数ｆ_CVjによって第２仮インク量Ｉ_SPj（＝仮デューティ制限値Ｄ_PIj）へと変換される。すなわち、設定テーブルＳＴＢに記憶された転用メディアの各インクのデューティ制限値Ｄ_Ij（ｊ＝１〜８）をそのまま使用するのではなく、仮デューティ制限値Ｄ_PIjに変換した上で使用する。従って、仮デューティ制限値Ｄ_PIjの制約条件の下でインク量Ｉ_jの最適化が実行されることとなる。

前記（３５），（３７）式に示すように、仮デューティ制限値Ｄ_PIjは、真のデューティ制限値Ｄ_Ij（ｊ＝１〜８）に対して正規化比ＲＷを一律に乗算したものであり、各インク間の真のデューティ制限値Ｄ_Ijの大きさの相対比は、仮デューティ制限値Ｄ_PIjにおいても維持されることとなる。従って、転用メディアにおける各インクの真のデューティ制限値Ｄ_Ij（ｊ＝１〜８）に応じたインク量Ｉ_jの最適化を実現することができ、転用メディアにおける各インクのデューティ制限値Ｄ_Ijの相対比に応じたガマットを形成することができる。

（３４）式が示すように、正規化比ＲＷは、被転用メディアの基準デューティ制限値Ｄ_SIj（ｊ＝１〜８）を、転用メディアのデューティ制限値Ｄ_Ij（ｊ＝１〜８）で除した制限値比Ｄ_SIj／Ｄ_Ijの最も小さいインク成分である。ここで、制限値比Ｄ_SIj／Ｄ_Ijが最も小さいインク成分ｊ＝ｚとし、それ以外のインク成分ｊ＝ｑとすると、下記の（４０）式が成り立つ。

前記の（４０）式の両辺に、転用メディアに対するインク成分ｊ＝ｚについてのデューティ制限値Ｄ_Izと、転用メディアに対するインク成分ｊ＝ｑについてのデューティ制限値Ｄ_Iqとを乗算すると、それぞれ下記の（４１）（４２）式が成り立つ。

前記の（４１）（４２）式の左辺のＲＷ・Ｄ_Iz，ＲＷ・Ｄ_Iqは、インク成分ｊ＝ｚ，ｑについての仮デューティ制限値Ｄ_Piz，Ｄ_PIqを意味する。インク成分ｊ＝ｚについては、仮デューティ制限値Ｄ_Pizが被転用メディアについての基準デューティ制限値Ｄ_SIzと等しくなる。従って、インク成分ｊ＝ｚについては、被転用メディアの各コンバーター等３００，３１０，４１０，１３６が色彩値等を予測可能なインク量Ｉ_jの範囲の全体にわたってインク量Ｉ_jの最適化を行うことができる。

一方、インク成分ｊ＝ｑについては、仮デューティ制限値Ｄ_Piqが被転用メディアについての基準デューティ制限値Ｄ_SIjよりも小さくなる。従って、被転用メディアの各コンバーター等３００，３１０，４１０，１３６が色彩値等を予測可能なインク量Ｉ_j（ｊ＝１〜８）の範囲内においてインク量Ｉ_jの最適化を行うことができる。いずれの場合も、被転用メディアの各コンバーター等３００，３１０，４１０，１３６が色彩値等を予測可能なインク量Ｉ_jの範囲でインク量Ｉ_jの最適化が行われるため、これらによる予測が破綻を来すことが防止できる。反対に、転用メディアのデューティ制限値Ｄ_Ijが、被転用メディアの基準デューティ制限値Ｄ_Ijを下回る場合にも、被転用メディアの基準デューティ制限値Ｄ_Ijによって制限されたインク量Ｉ_jが取り得る範囲を最大限に利用してインク量Ｉ_j（ｊ＝１〜４）を最適化することができる。

さらに、図６のステップＳ４００におけるベースＬＵＴ５１０を作成においては、スムージング／最適化処理によって最終的に得られた各格子点に対応するインク量Ｉ_j（ｊ＝１〜４）をベースＬＵＴ５１０，５２０の出力値としてそのまま登録するのではなく、該インク量Ｉ_j（ｊ＝１〜４）をインク量コンバーター７１０に設定された逆変換関数ｆ_CVj ^-1によって逆変換した上で、ベースＬＵＴ５１０，５２０の出力値として登録する。上述したように、スムージング／最適化処理によって最終的に得られた各格子点に対応するインク量Ｉ_j（ｊ＝１〜４）は、転用メディアに付着させる真のインク量Ｉ_j（ｊ＝１〜４）の取り得る範囲ではなく、被転用メディアについての各コンバーター等３００，３１０，４１０，１３６が色彩値等を予測可能なインク量Ｉ_jの範囲で最適化されている。この範囲のインク量Ｉ_jを逆変換関数ｆ_CVj ^-1によって逆変換することにより、転用メディアにおいてインク量Ｉ_j（ｊ＝１〜４）が真に取り得る範囲のインク量Ｉ_j（ｊ＝１〜４）へと復元することができる。（３７）式に示す変換関数ｆ_CVjの逆変換関数ｆ_CVj ^-1は、インク量Ｉ_j（ｊ＝１〜４）が取り得る範囲を復元する線形変換（１／ＲＷ）を含むということができる。

また、逆変換関数ｆ_CVj ^-1によって逆変換することにより、転用メディアにおける各インクの発色特性と、被転用メディアにおける各インクの発色特性との差分が補償できる。被転用メディアについての各コンバーター等３００，３１０，４１０，１３６は、被転用メディアにおける各インクの発色特性に対応した予測結果を出力するものであると捉えることができる。例えば、ある格子点について、図２１（Ａ）に発色特性を図示したＫインクのインク量Ｉ₄＝Ｉ_m（黒三角）が最適なインク量Ｉ₄であるとして最適化されたとする。このとき、該格子点については、インク量Ｉ_mが最適であると同時に、インク量Ｉ_mによる発色、すなわちインク量Ｉ_mのＫインクによって再現される明度Ｌ_m ^*が最適であると捉えることができる。

被転用メディアにおいて最適な明度Ｌ_m ^*が再現できるＫインクのインク量はＩ_mであるが、転用メディアにおいて最適な明度Ｌ_m ^*が再現できるＫインクのインク量はＩ_mではない。すなわち、被転用メディアにおいて最適な明度Ｌ_m ^*が再現できるＫインクのインク量Ｉ_mではなく、転用メディアにおいて最適な明度Ｌ_m ^*が再現できるＫインクのインク量Ｉ₄が転用メディアに関して真に最適なインク量であると考えることができる。転用メディアにおいて最適な明度Ｌ_m ^*が再現できるＫインクのインク量Ｉ₄は、被転用メディアにおいて最適な明度Ｌ_m ^*が再現できるＫインクのインク量Ｉ_m（第２仮インク量Ｉ_SPj）を逆変換関数ｆ_CVj ^-1によって逆変換することにより得ることができる。従って、スムージング／最適化処理によって最終的に得られた各格子点のインク量Ｉ_jを第２仮インク量Ｉ_SPjであると捉え、逆変換関数ｆ_CVj ^-1によって逆変換することにより、転用メディアについて最適なインク量Ｉ_jを得ることができる。このようにして得られた最適なインク量Ｉ_jをベースＬＵＴ５１０，５２０の出力値として登録することにより、転用メディアについて最適なインク量Ｉ_jが規定されたベースＬＵＴ５１０，５２０を得ることができる。（３７）式に示す変換関数ｆ_CVjの逆変換関数ｆ_CVj ^-1は、インクの発色特性に基づく非線形変換を含むということができる。

図３に示す重み指定用ＵＩ画像を表示させることにより、インク量Ｉ_jを最適化する際の目的関数Ｅにおける重みｗ_L*，ｗ_a*…は、特にユーザーが設定を変更しない限り、転用メディアの系統に適したデフォルト重みｗ_L*，ｗ_a*…が設定されることとなる。従って、転用メディアの系統に適した画質項目を重視するベースＬＵＴ５１０，５２０を作成することができる。メディアの指定に引き続いて重み指定用ＵＩ画像を表示することにより、メディアの性質や使用目的に応じた重みｗ_L*，ｗ_a*…を設定することができる。

ところで、転用メディアと被転用メディアは、それら自身が有する色味が異なるため、フォワードモデルコンバーター３００（分光プリンティングモデルコンバーター３１０）による色彩値の予測結果にずれが生じることとなる。むろん、Ｌ^*軸まわりについても転用メディアと被転用メディア自身の色味に応じて、色彩値の予測結果にずれが生じる。これに対して、本実施例では、グレーターゲットの色味（ａ_gt ^*，ｂ_gt ^*）を意図的に転用メディアと被転用メディアとの差分色味（ａ_D ^*，ｂ_D ^*）の逆色相方向にずらしているため、グレーターゲットに拘束したグレー軸格子点に対応する格子点の真の色彩値（フォワードモデルコンバーター３００による予測の色彩値ではなく、該格子点に対応するインク量Ｉ_jのインクを転用メディアに付着させた場合に再現される色彩値）をＬ^*軸上に位置させることができる。このようなベースＬＵＴ５１０によれば、真に無彩色（グレー）が再現できるインク量Ｉ_jを使用した補間演算を行うこと可能となる。従って、ベースＬＵＴ５１０，５２０に基づく補間演算を使用して作成した色補正ＬＵＴ６１０，６２０によれば、特に無彩色（グレー）の色再現性や階調性に優れた印刷結果を得ることができる。色補正ＬＵＴ６１０，６２０をプリンター１０のＲＯＭ１５に記憶させることにより、転用メディアについての色変換が可能となり、転用メディアに対して印刷を実行させることができるようになる。

Ｅ．変形例：
なお、この発明は前記の実施例や実施形態に限られるものではなく、その要旨を逸脱しない範囲において種々の態様において実施することが可能であり、例えば次のような変形も可能である。

Ｅ−１．変形例１：
前記実施例においては、発色特性も考慮して（３７）式に示す変換関数ｆ_CVj（ｆ_CVj ^-1）によってベースＬＵＴ５１０，５２０の出力値として登録するインク量Ｉ_jを変換したが、メディア間のデューティ制限値の差異のみを考慮した（３５）式によってインク量Ｉ_jを変換するようにしてもよい。この場合、発色特性の差異に起因する誤差が生じることとなるが、転用メディアにおける各インクのデューティ制限値を反映させたベースＬＵＴ５１０，５２０を作成することができる。

Ｅ−２．変形例２：
図２４は、メディア特性指定ＵＩ画像を示す図である。同図においては、被転用メディアと転用メディアの発色特性（彩度Ｃ^*値）が示されている。一般的に、各メディアにデューティ制限値Ｄ_Ij程度のインクを付着させると、それ以上インクを付着させても発色が変動しなくなる。そのため、Ｃ^*値の傾きの絶対値（｜δＣ^*／δＩ_j｜）が所定の基準値（≒０）と等しくなるインク量Ｉ_jをデューティ制限値Ｄ_Ijとしてもよい。図２４の例では、発色特性のグラフにおいて傾きの絶対値（｜δＣ^*／δＩ_j｜）が基準値と等しくなるインク量Ｉ_jにマーカー（白三角）を表示させるようにしている。これにより、ユーザーが、どのインク量Ｉ_jにデューティ制限値Ｄ_Ijが設定されるかを認識することができる。図２４の例では、ユーザーによるラジオボタンの選択により、デューティ制限値Ｄ_Ijを直接指定することも可能としている。これにより、マーカーの位置が、ユーザーがにじみを感じたカラーパッチのインク量Ｉ_jと明らかに異なる場合等にも対応することができる。

Ｅ−３．変形例３：
前記実施例では、転用メディアのデューティ制限値Ｄ_Ijと被転用メディアの基準デューティ制限値Ｄ_SIjとの比である正規化比ＲＷを使用して、転用メディアのデューティ制限値Ｄ_Ijを仮デューティ制限値Ｄ_PIjに変換することとしたが、正規化比ＲＷは他の指標によって定められてもよい。例えば、被転用メディアと転用メディアの発色特性の傾きの絶対値が一定の基準値（≠０）になったときのインク量Ｉ_jの比を正規化比ＲＷとしてもよい。

Ｅ−４．変形例４：
前記実施例では、機器独立表色系としてＣＩＥ−Ｌａｂ表色系を利用していたが、ＣＩＥ−ＸＹＺ表色系やＣＩＥ−Ｌ^*ｕ^*ｖ^*表色系などの他の任意の機器独立表色系を利用することが可能である。ただし、滑らかな色再現を実現するという意味からは、ＣＩＥ−Ｌａｂ表色系やＣＩＥ−Ｌ^*ｕ^*ｖ^*表色系などの均等色空間である機器独立表色系を用いることが好ましい。

Ｅ−５．変形例５：
前記実施例では、平滑化処理として力学モデルを利用した処理を採用していたが、他の種類の平滑化処理を採用してもよい。例えば、隣接する色彩値同士の間隔を測定し、その平均値になるべく近づくように個々の間隔を調整する平滑化処理を採用することも可能である。

Ｅ−６．変形例６：
本明細書において「インク」とは、インクジェットプリンタやオフセット印刷等に用いられる液体状インクに限らず、レーザプリンタに用いられるトナーも含む広い意味で使用されている。このような「インク」の広い意味を有する他の用語としては、「色材」や「着色材」、「着色剤」を用いることも可能である。

Ｅ−７．変形例７：
前記実施例では、ルックアップテーブルのような色変換プロファイルを作成する方法および装置に関して説明したが、本発明は、こうして得られた色変換プロファイルを印刷装置に組み込む組み込み部を備える印刷装置製造システムにも適用可能である。色変換プロファイルを作成する色変換プロファイル作成装置は、この印刷装置製造システムに含まれるものとしてもよく、他のシステムや装置に含まれるものとしてもよい。なお、この製造システムの組み込み部は、例えば、プリンタドライバのインストーラ（インストールプログラム）として実現することができる。

１０…プリンター、１３…ＣＰＵ、１４…ＲＡＭ、１５…ＲＯＭ、２０…メモリーカードスロット、２１…バス、２２…ＡＳＩＣ、１５ａ…プログラムデータ、１００…ベースＬＵＴ作成モジュール、１２０…スムージング処理初期値設定モジュール、１３０…スムージング処理モジュール、１３２…色点移動モジュール、１３４…インク量最適化モジュール、１３６…画質評価指数コンバーター、１４０…テーブル作成モジュール、２００…色補正ＬＵＴ作成モジュール、３００…フォワードモデルコンバーター、３１０…分光プリンティングモデルコンバーター、３２０…色算出部、４００…ＬＵＴ格納部、４１０…インバースモデル初期ＬＵＴ、５１０…ベース３Ｄ−ＬＵＴ、５１１…逆変換ＬＵＴ、５２０…ベース４Ｄ−ＬＵＴ、５２１…逆変換ＬＵＴ、６１０…色補正３Ｄ−ＬＵＴ、６２０…色補正４Ｄ−ＬＵＴ、７００…ＬＵＴ作成条件設定モジュール、７１０…インク量コンバーター、７２０…ＵＩモジュール、７３０…設定情報格納モジュール。

Claims

プロファイル作成方法であって、
機器独立表色系の複数の格子点の分布を平滑化する平滑化処理と、平滑化された前記複数の格子点が示す色彩値を基本印刷媒体において再現させるインク量を、基本印刷媒体にインクを付着したときの画質を評価する目的関数を用いた最適化によって決定する最適化処理とを実行し、
前記最適化処理における前記インク量の範囲を、前記基本印刷媒体に付着可能な前記インク量の制限値と、前記基本印刷媒体とは異なる転用印刷媒体に付着可能な前記インク量の制限値との比である制限値比に基づいて制限し、
前記平滑化処理と前記最適化処理で決定された前記インク量を前記制限値比に基づいて変換した前記インク量と、該インク量に対応する格子点の色彩値との対応関係を規定することにより前記転用印刷媒体のためのプロファイルを作成することを特徴とするプロファイル作成方法。
前記制限値比は、前記基本印刷媒体に付着可能な前記インク量の制限値を、前記転用印刷媒体に付着可能な前記インク量の制限値で除算したものであり、
前記最適化処理における前記インク量の範囲は、前記転用印刷媒体に付着可能な前記インク量の制限値に前記制限値比を乗算した仮制限値によって制約されることを特徴とする請求項１に記載のプロファイル作成方法。
前記平滑化処理と前記最適化処理で決定された前記インク量を、前記制限値比によって除算した前記インク量と、該インク量に対応する前記機器独立表色系の色彩値との対応関係を規定することにより前記プロファイルを作成することを特徴とする請求項２に記載のプロファイル作成方法。
前記基本印刷媒体に付着可能な前記インク量の制限値と、前記転用印刷媒体に付着可能な前記インク量の制限値は、前記基本印刷媒体と前記転用印刷媒体においてインクによる発色の変動が一定の基準よりも小さくなる前記インク量であることを特徴とする請求項１から請求項３のいずれか一項に記載のプロファイル作成方法。
前記基本印刷媒体に付着可能な前記インク量の制限値と、前記転用印刷媒体に付着可能な前記インク量の制限値は、前記基本印刷媒体と前記転用印刷媒体においてインクによるにじみが観察される前記インク量であることを特徴とする請求項１から請求項３のいずれ一項に記載のプロファイル作成方法。
プロファイル作成プログラムであって、
機器独立表色系の複数の格子点の分布を平滑化する平滑化処理と、平滑化された前記複数の格子点が示す色彩値を基本印刷媒体において再現させるインク量を、基本印刷媒体にインクを付着したときの画質を評価する目的関数を用いた最適化によって決定する最適化処理とをコンピューターに実行させ、
前記最適化処理における前記インク量の範囲を、前記基本印刷媒体に付着可能な前記インク量の制限値と、前記基本印刷媒体とは異なる転用印刷媒体に付着可能な前記インク量の制限値との比である制限値比に基づいて制限する機能と、
前記平滑化処理と前記最適化処理で決定された前記インク量を前記制限値比に基づいて変換した前記インク量と、該インク量に対応する格子点の色彩値との対応関係を規定することにより前記転用印刷媒体のためのプロファイルを作成する機能とをコンピューターに実行させることを特徴とするプロファイル作成プログラム。
プロファイルを参照して色変換を行うことにより得られたインク量のインクを転用メディアに付着させる印刷装置であって、
前記プロファイルは、
機器独立表色系の複数の格子点の分布を平滑化する平滑化処理と、平滑化された前記複数の格子点が示す色彩値を基本印刷媒体において再現させるインク量を、基本印刷媒体にインクを付着したときの画質を評価する目的関数を用いた最適化によって決定する最適化処理とを実行し、
前記最適化処理における前記インク量の範囲を、前記基本印刷媒体に付着可能な前記インク量の制限値と、前記基本印刷媒体とは異なる転用印刷媒体に付着可能な前記インク量の制限値との比である制限値比に基づいて制限し、
前記平滑化処理と前記最適化処理で決定された前記インク量を前記制限値比に基づいて変換した前記インク量と、該インク量に対応する格子点の色彩値との対応関係を規定することにより作成されることを特徴とする印刷装置。
機器独立表色系の複数の格子点の分布を平滑化する平滑化処理と、平滑化された前記複数の格子点が示す色彩値を基本印刷媒体において再現させるインク量を、基本印刷媒体にインクを付着したときの画質を評価する目的関数を用いた最適化によって決定する最適化処理とを実行し、前記最適化処理における前記インク量の範囲を、前記基本印刷媒体に付着可能な前記インク量の制限値と、前記基本印刷媒体とは異なる転用印刷媒体に付着可能な前記インク量の制限値との比である制限値比に基づいて制限し、前記平滑化処理と前記最適化処理で決定された前記インク量を前記制限値比に基づいて変換した前記インク量と、該インク量に対応する格子点の色彩値との対応関係を規定することにより前記転用印刷媒体のためのプロファイルを作成するにあたり、
前記異なる転用印刷媒体に付着可能な前記インク量の制限値の指定を受け付けるための画像を表示装置に表示させることを特徴とする表示制御装置。