(57)【要約】【課題】複数の観察光源において色の見えが一致する印刷結果を提供する。【解決手段】ＬＵＴを参照しつつ入力画像データＩＤを変換して生成した記録材量データに基づいてプリンター３０に印刷を実行させるコンピューター２０において、入力画像データの表色系とプリンター３０のインク量データの表色系とを直接に対応付けるとともにそれぞれ特定の観察光源に対応付けて作成された複数の色変換テーブルを記憶した記憶媒体２４を備えており、複数の色変換テーブルは、色変換テーブルを参照しつつ印刷された印刷結果を該色変換テーブルに対応付けられた特定の観察光源下で観察したときに同一の入力画像データに基づく印刷結果同士で色の見えが一致するように作成されている。【選択図】図３
2011139213.tif 000002
【特許請求の範囲】
【請求項１】
色変換テーブルを参照しつつ入力画像データを変換して生成した記録材量データに基づいて印刷装置に印刷を実行させる印刷制御装置であって、
前記入力画像データの表色系と前記記録材量データの表色系とを直接に対応付けるとともにそれぞれ特定の観察光源に対応付けて作成された複数の色変換テーブルを記憶した記憶媒体を備えており、
前記複数の色変換テーブルは、色変換テーブルを参照しつつ印刷された印刷結果を該色変換テーブルに対応付けられた特定の観察光源下で観察したときに同一の入力画像データに基づく印刷結果同士で色の見えが一致するように作成されていることを特徴とする印刷制御装置。
【請求項２】
前記観察光源の種類を指定する操作入力を受付ける操作入力部を備え、
前記印刷制御装置は、前記操作入力部で指定された観察光源に対応する色変換テーブルを参照しつつ前記入力画像データを変換して生成した記録材量データに基づいて前記印刷装置に印刷を実行させる請求項１に記載の印刷制御装置。
【請求項３】
前記複数の色変換テーブルは、前記入力画像データの色域と前記印刷装置の色域とを同一のガマットマッピングによって対応付けて作成されている請求項１または請求項２に記載の印刷制御装置。
【請求項４】
前記複数の色変換テーブルにそれぞれ対応付けられた複数の観察光源の中から２以上の観察光源の選択を受付けるとともにその混合比率を受付けるＬＵＴ合成条件受付部と、
前記ＬＵＴ合成条件受付部の受付けた観察光源に対応付けられた色変換テーブルを前記混合比率で合成して新たな色変換テーブルを作成するＬＵＴ合成部とを備える請求項１〜請求項３の何れか１項に記載の印刷制御装置。
【請求項５】
前記ＬＵＴ合成条件受付部は、前記新たな色変換テーブルに対応する観察光源と受付けた各観察光源との間の色温度差の比を前記混合比率とする請求項４に記載の印刷制御装置。
【請求項６】
前記複数の色変換テーブルにそれぞれ対応付けられている複数の観察光源は、色温度が約３０００Ｋ〜約７０００Ｋの範囲の全域に離散的に分布している請求項１〜請求項５の何れか１項に記載の印刷制御装置。
【請求項７】
色変換テーブルを参照しつつ入力画像データを変換して生成した記録材量データに基づいて印刷装置に印刷を実行させる機能をコンピューターに実現させるための印刷制御プログラムであって、
前記入力画像データの表色系と前記記録材量データの表色系とを直接に対応付けるとともにそれぞれ特定の観察光源に対応付けて作成された複数の色変換テーブルを備えており、
前記複数の色変換テーブルは、色変換テーブルを参照しつつ印刷された印刷結果を該色変換テーブルに対応付けられた特定の観察光源下で観察したときに同一の入力画像データに基づく印刷結果同士で色の見えが一致するように、作成されていることを特徴とする印刷制御プログラム。
【請求項８】
色変換テーブルを参照しつつ入力画像データを変換して生成した記録材量データに基づいて印刷を実行する印刷装置であって、
前記入力画像データの表色系と前記記録材量データの表色系とを直接に対応付けるとともにそれぞれ特定の観察光源に対応付けて作成された複数の色変換テーブルを記憶した記憶媒体を備えており、
前記複数の色変換テーブルは、色変換テーブルを参照しつつ印刷された印刷結果を該色変換テーブルに対応付けられた特定の観察光源下で観察したときに同一の入力画像データに基づく印刷結果同士で色の見えが一致するように、作成されていることを特徴とする印刷装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、印刷制御装置、印刷制御プログラムおよび印刷装置に関し、特に、色変換テーブルを参照しつつ印刷データに基づく印刷を印刷装置に実行させる印刷制御装置、印刷制御プログラムおよび印刷装置に関する。
【背景技術】
【０００２】
色の見えが観察光源しだいで変化するという光源依存性が知られている。特に印刷物においては、作成者の調整した色を正しく観察できることが重要である。引用文献１には、撮影された被写体を所望光源下で観察する場合の色を、プリンター上で容易に再現する技術が開示されている。引用文献２には、メタメリズム改善用インクを用いずに、光源が変わっても出力画像の無彩色が同じ色に見えるようにする技術が開示されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開平９−２８４５８１号公報
【特許文献２】特開２００５−５７５７０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
上述した特許文献１の技術は、メタメリズムを改善するものであって、撮影環境ごとに色の見えを正しく再現するためのものではなかった。また、上述した特許文献２の技術は、無彩色部分について各光源間での色相のズレが最小限になるようなインク構成を用いるため、撮影環境毎に色の見えを正しく再現するものではなかった。
本発明は上記課題に鑑みてなされたものであり、複数の観察光源において色の見えが一致する印刷結果を提供可能な印刷制御装置、印刷制御プログラムおよび印刷装置を提供する事を目的とする。
【課題を解決するための手段】
【０００５】
前記課題を解決するために、本発明にかかる印刷制御装置は、色変換テーブルを参照しつつ入力画像データを変換して生成した記録材量データに基づいて印刷装置に印刷を実行させるために、複数の色変換テーブルを記憶した記憶媒体を備えている。該記憶媒体に記憶された複数の色変換テーブルは、入力画像データの表色系と記録材料データの表色系とを直接に対応付けたデータとして作成されており、それぞれ特定の観察光源に対応させて作成されたものである。すなわち、一の入力画像データについてそれぞれ複数の色変換テーブルのひとつを参照しつつ印刷して複数の印刷結果を作成したとき、一の色変換テーブルを参照しつつ印刷された印刷結果を該一の色変換テーブルに対応付けられた特定の観察光源下で観察したときの色の見えと、他の一の色変換テーブルを参照しつつ印刷された印刷結果を該他の一の色変換テーブルに対応付けられた特定の観察光源下で観察したときの色の見えとが、一致するようになっている。したがって、複数の光源下でそれぞれ観察したときに、同じ色の見えが再現される印刷結果を作成できるようになる。また、印刷に利用された色変換テーブルに対応する観察光源下で印刷結果を観察すれば、入力画像データの作成者の意図を反映した色の見えが再現されることになる。
【０００６】
また、本発明の選択的な一態様として、前記観察光源の種類を指定する操作入力を受付ける操作入力部を備え、該操作入力部で観察光源の種類が指定されると、前記印刷制御装置は、前記操作入力部で指定された観察光源に対応する色変換テーブルを参照しつつ前記入力画像データを変換して生成した記録材量データに基づいて前記印刷装置に印刷を実行させるようにしてもよい。このように構成すると、ユーザの希望する観察光源に対応した印刷物を作成し、ユーザが想定していない不要な観察光源に対応した印刷物は生成されない。むろん、複数の観察光源の種類を指定することにより、複数の観察光源下のそれぞれで同じ見えになる印刷物を作成することもできる。
【０００７】
また、本発明の選択的な一態様として、前記複数の色変換テーブルは、前記入力画像データの色域と前記印刷装置の色域とを同一のガマットマッピングによって対応付けて作成されていてもよい。色変換テーブルの入力と出力とが複数の色変換テーブルの間で一致しているので、これに加えてガマットマッピングの種類も一致していることにより、印刷結果に再現される色の見えもより一致しやすくなる。
【０００８】
また、本発明の選択的な一態様として、前記複数の色変換テーブルに対応付けられた複数の観察光源の中から２以上の観察光源の選択を受付けるとともにその混合比率を受付けるＬＵＴ合成条件受付部と、前記ＬＵＴ合成条件受付部の受付けた観察光源に対応付けられた色変換テーブルを前記混合比率で合成して新たな色変換テーブルを作成するＬＵＴ合成部とを備えるようにしてもよい。このように予め前記記憶媒体に記憶されている複数の色変換テーブルから新たな観察光源に対応した色変換テーブルを合成できるようにすることにより、記憶媒体のリソースを節約することができる。また、多様な観察光源に対応した色変換テーブルに対応することができるようになる。
【０００９】
また本発明の選択的な一態様として、前記複数の色変換テーブルにそれぞれ対応付けられている複数の観察光源は、色温度が約３０００Ｋ〜約７０００Ｋの範囲の全域に離散的に分布するように構成してもよい。このように構成すると、利用者が希望する大方の観察光源に対応する色変換テーブルを精度よく合成することができるようになる。
【００１０】
上述した印刷制御装置は、他の機器に組み込まれた状態で実施されたり他の方法とともに実施されたりする等の各種の態様を含む。また、本発明は前記印刷制御装置を備える印刷制御システム、上述した装置の構成に対応した工程を有する印刷制御方法、上述した装置の構成に対応した機能をコンピューターに実現させる印刷制御プログラム、該プログラムを記録したコンピューター読み取り可能な記録媒体、上述した印刷制御装置と同等の機能を備えつつ印刷機能をも備えている印刷装置等としても実現可能である。これら印刷システム、印刷方法、印刷制御方法、印刷プログラム、該プログラムを記録した媒体、の発明も、上述した作用、効果を奏する。むろん、請求項２〜６に記載した構成も、前記システムや前記方法や前記プログラムや前記記録媒体に適用可能である。
【図面の簡単な説明】
【００１１】
【図１】本実施形態の概略構成を示す説明図である。
【図２】コンピューター２０のハードウェア構成を示すブロック図である。
【図３】ＬＵＴと観察光源の対応関係を説明する図である。
【図４】印刷制御処理の流れを示すフローチャートである。
【図５】光源選択ＵＩの一例を示す図である。
【図６】ＬＵＴ合成処理の流れを示すフローチャートである。
【図７】合成条件受付ＵＩの例を示す図である。
【図８】ミックスＬＵＴ作成を説明する図である。
【図９】コンピュータ１０の構成を示すブロック図である。
【図１０】ＬＵＴ作成処理の手順を示すフローチャートである。
【図１１】ベース３Ｄ−ＬＵＴを作成する場合の処理内容を示す説明図である。
【図１２】ＬＵＴ作成条件ＵＩの一例を示す図である。
【図１３】入力格子点とＬａｂ表色系の色点との対応関係を示す図である。
【図１４】ベース４Ｄ−ＬＵＴを作成する場合の処理内容を示す説明図である。
【図１５】ベースＬＵＴを用いた色補正ＬＵＴの作成方法を示す説明図である。
【図１６】力学モデルの説明図である。
【図１７】スムージング処理の典型的な処理手順を示すフローチャートである。
【図１８】ステップＴ１００の詳細手順を示すフローチャートである。
【図１９】ステップＴ１２０〜Ｔ１５０の処理内容を示す説明図である。
【図２０】調整係数の一例を示す模式図である。
【図２１】最適化処理の詳細手順を示すフローチャートである。
【図２２】プリンター３０の構成を示す図である。
【図２３】プリンター３０の印刷方式を示す模式図である。
【図２４】ファームウェアのソフトウェア構成を示す図である。
【発明を実施するための形態】
【００１２】
以下、下記の順序に従って本発明の実施形態を説明する。
（１）本実施形態の構成：
（２）印刷制御処理：
（３）ミックスＬＵＴの作成：
（４）ＬＵＴ作成処理：
（５）印刷装置の構成
（６）その他変形例：
（７）まとめ：
【００１３】
（１）本実施形態の構成： 
図１は本実施形態にかかる概略構成を示す図である。同図には、印刷制御装置としてのコンピューター２０と、印刷装置としてのプリンター３０と、コンピューター２０に格納されている色変換テーブル（以下、ＬＵＴと略す）を作成するＬＵＴ作成装置としてのコンピューター１０とを示してある。なお、コンピューター２０とプリンター３０の組合せが印刷装置を構成するとも言える。コンピューター２０にはＬＵＴが記憶されており、このＬＵＴを参照しつつコンピューター２０がプリンター３０の印刷を制御する。
【００１４】
図２は、コンピューター２０のハードウェア構成を示すブロック図である。同図において、コンピューター２０は、ハードディスク（ＨＤ）等で構成された不揮発性の記憶媒体２４と、ユーザーの操作入力を受付けるためのキーボードやマウス等の操作入力部２２と、液晶ディスプレイやＣＲＴ等で構成される表示部２３と、プリンター３０と通信する際のインターフェース２５（Ｉ／Ｆ２５）と、プリンター３０の印刷を制御する印刷制御部２１とを備えている。記憶媒体２４には、複数のＬＵＴが記憶されている。
【００１５】
記憶媒体２４に記憶されているＬＵＴは、標準的な機器依存表色系のデータ（以下、「機器依存入力データ」と記載する。）を入力とし、特定のプリンターのインク量データを出力とする対応関係を規定している。標準的な機器依存表色系とは、例えばｓＲＧＢ表色系やＪａｐａｎＣｏｌｏｒ２００１表色系などである。インク量データとは、特定のプリンターに搭載されているインク種類ごとのインク量の組合せで構成されるデータであり、本実施形態のプリンター２０はＣ（シアン），Ｍ（マゼンタ），Ｙ（イエロー），Ｋ（ブラック），Ｌｃ（淡シアン），Ｌｍ（淡マゼンタ），Ｌｋ（淡ブラック），ＬＬｋ（極淡ブラック）の８種類のインクを利用可能なカラープリンターを想定しているため、インク量データはこれら８種類のインクのインク量の組合せで構成されているものとして以下の説明を行うことにする。
【００１６】
本実施形態においては記憶媒体２４に、ＬＵＴ（Ｄ５０） sRGB、ＬＵＴ（Ａ） sRGB、ＬＵＴ（Ｃ） sRGB、ＬＵＴ（Ｄ５０） JP2001、ＬＵＴ（Ａ） JP2001、ＬＵＴ（Ｃ） JP2001、が記憶されているものとして説明を行う。ＬＵＴ（Ｄ５０） sRGB、ＬＵＴ（Ａ） sRGB、ＬＵＴ（Ｃ） sRGBは、ｓＲＧＢ表色系を入力とし、ＣＭＹＫＬｃＬｍＬｋＬＬｋを出力としている。ＬＵＴ（Ｄ５０） JP2001、ＬＵＴ（Ａ） JP2001、ＬＵＴ（Ｃ） JP2001は、ＪａｐａｎＣｏｌｏｒ２００１表色系を入力とし、ＣＭＹＫＬｃＬｍＬｋＬＬｋを出力としている。かっこ内の文字はＬＵＴの対応する観察光源を示し、下付きの添字は入力表色系を示している。すなわち、記憶媒体２４には、入力表色系と出力表色系とが同一の複数のＬＵＴが記憶されており、これら複数のＬＵＴはたがいに異なる観察光源に対応している。
【００１７】
なお説明の簡略のため、ここでは印刷用紙の種類やＬＵＴ作成時のガマットマッピングの種類を１種類のみとして説明を行うが、印刷用紙やガマットマッピングの種類を選択するようにしてもよい。異なる種類の印刷用紙ではインクの発色特性が異なるため、印刷用紙ごとに異なるＬＵＴが用意される。すなわち、印刷用紙の種類ごとに上述したＬＵＴ（Ｄ５０） sRGB、ＬＵＴ（Ａ） sRGB、ＬＵＴ（Ｃ） sRGB、ＬＵＴ（Ｄ５０） JP2001、ＬＵＴ（Ａ） JP2001、ＬＵＴ（Ｃ） JP2001と同じ数のＬＵＴが作成されて記憶媒体２４に記憶されている。従って、印刷用紙の種類（普通紙、マット紙、光沢紙など）を選択すると、その印刷用紙に最適化されたＬＵＴを参照して印刷が行われる。また、一般的に入力表色系と出力表色系の色域は異なるため、ＬＵＴを作成する際にいわゆるガマットマッピングが行われるが、このガマットマッピングの種類ごとにも異なるＬＵＴが用意される。すなわち、ガマットマッピングの種類ごとに上述したＬＵＴ（Ｄ５０） sRGB、ＬＵＴ（Ａ） sRGB、ＬＵＴ（Ｃ） sRGB、ＬＵＴ（Ｄ５０） JP2001、ＬＵＴ（Ａ） JP2001、ＬＵＴ（Ｃ） JP2001と同じ数のＬＵＴが作成されて記憶媒体２４に記憶されている。従って、ガマットマッピングの種類（明度維持、彩度維持など）を選択すると、そのガマットマッピングで作成されたＬＵＴを参照して印刷が行われる。
【００１８】
また、本実施形態では、コンピューター１０で作成された複数のＬＵＴが予め記憶媒体２４に保存されているが、コンピューター２０に後述のＬＵＴ作成処理を実行するためのプログラムを備えさせ、後述の光源選択部Ｍ２に対して観察光源が指定されてからコンピューター２０においてＬＵＴを作成するようにしてもよい。このときのＬＵＴ作成手法は、後述のＬＵＴ作成処理と同様の手順で行う。この場合も、色変換の精度の低下は防止されるし、多様な観察光源に対応したＬＵＴを作成可能である。
【００１９】
図３は、ＬＵＴと観察光源の対応関係を説明する図である。同図では、同一の入力画像データＩＤをそれぞれＬＵＴ（Ｄ５０） sRGB，ＬＵＴ（Ａ） sRGB，ＬＵＴ（Ｃ） sRGBで色変換してプリンター３０で印刷した印刷結果を、対応する観察光源下で観察した場合について対比してある。
【００２０】
同図において、１つめの印刷結果ＰＲ１は、ＬＵＴ（Ｄ５０） sRGBを参照しつつ入力画像データＩＤを後述の色変換部Ｍ３で色変換して印刷したものであり、印刷結果ＰＲ１の所定部位における分光反射率分布をＲ t1（λ）としてある。なお、所定部位は、印刷結果ＰＲ１〜ＰＲ３のそれぞれにおいて対応して印刷された部位であり、印刷結果ＰＲ１〜ＰＲ３のそれぞれで入力画像データＩＤの共通するデータで印刷された部位である。ここで、印刷結果ＰＲ１に観察光源Ｄ５０を照射して観察すると、人の目には下記（１）式の３刺激値Ｘ１，Ｙ１，Ｚ１に相当する色が見えることになる。下記式（１）において、Ｐ D50（λ）は観察光源Ｄ５０の分光強度分布であり、ｋは正規化のための係数である。
【数１】
2011139213.tif 000003
【００２１】
同様に、２つめの印刷結果ＰＲ２は、ＬＵＴ（Ａ） sRGBを参照しつつ入力画像データＩＤを後述の色変換部Ｍ３で色変換して印刷したものであり、印刷結果ＰＲ２の所定部位における分光反射率分布をＲ t2（λ）としてある。印刷結果ＰＲ２に観察光源Ａを照射して観察すると、人の目には下記（２）式の３刺激値Ｘ２，Ｙ２，Ｚ２に相当する色が見えることになる。下記式（２）において、ＰA（λ）は観察光源Ａの分光強度分布であり、ｋは正規化のための係数である。
【数２】
2011139213.tif 000004
【００２２】
同様に、３つめの印刷結果ＰＲ３は、ＬＵＴ（Ｃ） sRGBを参照しつつ入力画像データＩＤを色変換部Ｍ３で色変換して印刷したものであり、印刷結果ＰＲ３の所定部位における分光反射率分布をＲ t3（λ）としてある。印刷結果ＰＲ３に観察光源Ａを照射して観察すると、人の目には下記（３）式の３刺激値Ｘ３，Ｙ３，Ｚ３に相当する色が見えることになる。下記式（３）において、ＰC（λ）は観察光源Ａの分光強度分布であり、ｋは正規化のための係数である。
【数３】
2011139213.tif 000005
【００２３】
以上のようにして各印刷結果ＰＲ１〜ＰＲ３の測色結果に基づいて算出される三刺激値は、ほぼＸ１＝Ｘ２＝Ｘ３、Ｙ１＝Ｙ２＝Ｙ３、Ｚ１＝Ｚ２＝Ｚ３となっており、互いに等しい色を示すようになっている。より正確には、各三刺激値を所定の変換式で変換することにより得られる色彩値が、互いに等しい色を示すようになっている。すなわち、各印刷結果ＰＲ１〜ＰＲ３を、各印刷結果の印刷に利用したＬＵＴの対応する光源下で人間が観察すると、色の見えが等しくなることになる。このように互いに異なる観察光源に対応しつつ入力表色系と出力表色系とが同一になっている複数のＬＵＴを備えたときの印刷制御処理は、以下のようにして行う。
【００２４】
（２）印刷制御処理：
図４は、コンピューター２０にて実行される印刷制御処理の流れを示すフローチャートである。この印刷制御処理は、コンピューター２０で実行されるアプリケーションプログラムなどからプリンター３０を利用した印刷が指示されたときに、印刷制御部２１が実行する。印刷制御部２１は、例えば、コンピューター２０でプリンター３０のプリンタードライバープログラムを実行することにより実現される。むろん、印刷制御処理は、印刷制御部２１と同等の機能を備えたＡＳＩＣ等のカスタムＩＣのようなハードウェアによって実行されても構わない。
【００２５】
印刷制御部２１は、印刷制御処理を実行するために、アプリケーションプログラムなどから画像や文書などの印刷データを印刷対象として取得して該印刷データに基づいて印刷に使用する入力画像データＩＤを生成するデータ取得部Ｍ１と、いずれの観察光源に対応したＬＵＴを参照して色変換を実行するかを選択するための光源選択部Ｍ２と、選択されたＬＵＴを参照しつつ入力画像データＩＤからインク量データを生成する色変換部Ｍ３と、ディザ法や誤差拡散法等によってインク量データからドットのオンオフを示す二値データを生成するハーフトーン処理部Ｍ４と、二値データに基づいてラスタライズ処理等を行って印刷制御データを生成してプリンター３０に出力するコマンド作成部Ｍ５と、に相当する機能を備えている。本実施形態においては、インク量データが記録材量データを構成する。
【００２６】
処理が開始されると、データ取得部Ｍ１はアプリケーションプログラムから印刷データを取得するとともに入力画像データＩＤを生成する（ステップＳ１０）。次に、光源選択部Ｍ２は図５に示すような光源選択ＵＩを表示部２３に表示する（ステップＳ２０）。光源選択ＵＩには、Ｄ５０，Ａ，Ｃなどの観察光源の種類が表示されており、ユーザーは操作入力部２２を操作することにより所望の観察光源の種類を選択することができる。この光源選択ＵＩにおいては、複数の観察光源を選択してもよく、複数の観察光源を選択すると各観察光源に対応した印刷結果がそれぞれに作成されることになる。
【００２７】
観察光源が指定されると、色変換部Ｍ３はステップＳ２０で指定された印刷条件に適合するＬＵＴを記憶媒体２４から取得する（ステップＳ３０）。入力画像データＩＤにはどのような表色系で作成されたデータであるかを示す情報がヘッダ等に記載されているので、色変換部Ｍ３はこの情報に基づいて、入力表色系と観察光源の一致するＬＵＴを記憶媒体２４から取得する。そして、色変換部Ｍ３は、ステップＳ３０で取得したＬＵＴを参照しつつ入力画像データＩＤをプリンター３０のインク量データに変換する（ステップＳ４０）。
【００２８】
その後、ハーフトーン処理部Ｍ４がステップＳ４０で作成されたインク量データにハーフトーン処理を行って２値データを作成し（ステップＳ５０）、２値データに基づいてプリンター３０の印刷制御データを作成してプリンター３０に出力する（ステップＳ６０）。このようにして印刷された印刷物は、ＬＵＴの対応する観察光源下で観察すると、入力画像データＩＤの作成者の意図を良好に反映した色の見えとなる。また、ステップＳ２０において複数の観察光源を選択して印刷した場合は、各印刷結果をそれぞれの対応する観察光源下で観察すると、その両者で色の見えが一致することになる。
【００２９】
以上の実施形態によれば、プロファイル変換処理が不要である。一般的なプリンターでは、例えばＩＣＣプロファイルのように、標準的な機器依存表色系とＣＩＲＬＡＢなどの機器独立表色系との対応関係を規定したＬＵＴを記憶しており、色変換を行う際に、例えばソースプロファイルとデスティネーションプロファイルの双方を取得するとともにガマットマッピングを行い、標準的な機器依存表色系のデータをプリンター２００のインク量データに変換するためのプロファイルを作る処理が必要になる。これに対し、本実施形態では、このような処理が不要であるためリソースが節約できる。また、一般的なプリンターでは、例えば入力データがｓＲＧＢであれば、ＲＧＢ値からＬａｂ値への変換とＬａｂ値からＣＭＹＫＬｃＬｍＬｋＬＬｋへの変換の２回の変換が発生するため変換精度が低下する可能性があったが、本実施形態では複数の色変換が発生しないため色変換の精度の低下を防止することができる。
【００３０】
（３）ミックスＬＵＴの作成：
以上説明した実施形態では、予め作成したＬＵＴを参照して印刷を行っていた。しかしながら、多くの観察光源に対応したＬＵＴを予め作成して保存しておくと、リソースを圧迫する。また、観察光源を指定されてから後述のＬＵＴ作成処理をおこなってＬＵＴを作成する場合も、ＬＵＴ作成プログラムを記憶するためのリソースが必要になるため、やはりリソースを圧迫することになる。そこで、本変形例では、予め記憶媒体２４に記憶しておくＬＵＴを必要最低限の種類数の観察光源に対応したＬＵＴとしておいて、他の種類の観察光源に対応するＬＵＴについては作成済みのＬＵＴを利用して合成できるようにしてある。
【００３１】
図６は、ＬＵＴ合成処理の流れを示すフローチャートである。ＬＵＴ合成処理は、記憶媒体２４に記憶されているＬＵＴの対応する観察光源以外の観察光源に対応するＬＵＴを、記憶媒体２４に記憶されているＬＵＴを利用して作成する処理である。この処理は、コンピューター２０において実行される処理であり、例えば、上述した印刷制御処理のステップＳ２０において呼び出されて実行される。ステップＳ２０でＬＵＴ合成処理を呼び出して実行するためには、例えば、光源選択ＵＩに「他の光源」などのボタンを設けておき、このボタンを操作入力部２２で操作したときに実行するなどすればよい。むろん、印刷制御処理とは別個にＬＵＴ合成処理を実行しても構わない。印刷制御部２１は、ＬＵＴ合成処理を実行するためにＬＵＴ合成条件受付部Ｍ６とＬＵＴ合成部Ｍ７とを備えている。
【００３２】
処理が開始されると、ＬＵＴ合成条件受付部Ｍ６は合成条件受付ＵＩを表示部３２に表示する（ステップＳ３００）。図７（ａ）〜（ｃ）は、２以上のＬＵＴの合成比率を指定するための合成条件受付ＵＩの例を示す図である。これらのＵＩで指定されるＬＵＴの合成比率は、合成したいＬＵＴの観察光源（以下、ミックス光源と呼ぶことにする）を、作成済みのＬＵＴに対応付けられた観察光源（以下、既存光源と呼ぶことにする）を合成して再現する際の合成比率に対応している。すなわち、合成条件受付ＵＩでは、ユーザーの想定している観察光源（ミックス光源）を何らかの方法で指定することになる（ステップＳ３１０）。
【００３３】
例えば、図７（ａ）に示す例では、２つの既存光源を一覧から選択するとともに、観察光源の合成比率を数値で指定するようになっている。むろん２以上の既存光源を選択して、合成比率を指定するようにしても構わない。ユーザーは、例えば、所望の観察光源と類似する２つの光源を指定し、どのような比率で２つの光源を混合すれば所望の観察光源を再現できるか判断し、その比率を混合比率として入力する。
【００３４】
また、図７（ｂ）、（ｃ）に示す例では、色温度により所望の観察光源を指定できるようになっている。図７（ｂ）では、光源の色温度軸をスライダーバーとして表示してあり、スライダーＳＬで色温度軸の１点を指定することにより、所望の観察光源の色温度を指定できるようになっている。また、図６（ｃ）では、ｘｙ色度図を２次元で表示してあり、このｘｙ色度図上の１点をカーソルで指定することにより、所望の観察光源の色温度を指定できるようになっている。色温度は比較的ユーザーにとって馴染みのある概念であり、所望の観察光源を指定しやすくなると考えられる。また、色温度は、１つの数値で任意の光源を指定できるため、観察光源の名称などを知らないユーザーであっても指定が容易である。むろん、色温度を数値で入力して指定できるようにしてもよいし、観察光源を特定するための情報として各光源のスペクトル分布を示したり、各光源の名称を示したり、分光光度計や測色計等によって実際の環境光を測定したデータによって観察光源を特定できるようにしても構わない。
【００３５】
以上のようにしてこのようにしてユーザーの想定している観察環境における観察光源が特定されてＬＵＴ合成条件受付部Ｍ６がＬＵＴ合成条件を受付けると（ステップＳ３１０）、ＬＵＴ合成部Ｍ７は、ミックスＬＵＴの合成を行う（ステップＳ３２０）。まず、ＬＵＴ合成部Ｍ７は、予め作成して保持しているＬＵＴ（Ｄ５０） sRGB、ＬＵＴ（Ａ） sRGB、ＬＵＴ（Ｃ） sRGB、ＬＵＴ（Ｄ５０） JP2001、ＬＵＴ（Ａ） JP2001、ＬＵＴ（Ｃ） JP2001で想定されているＤ５０光源，Ａ光源，Ｃ光源によってユーザー所望の観察光源ＬＳｘを再現するための混合比率を算出する。すなわち、下記式（１）における、係数ｃ１，ｃ２，ｃ３を特定する。なお、図６（ａ）のようにユーザーが予め観察光源の混合比率を指定している場合は、混合比率の算出は行わない。
【数４】
2011139213.tif 000006
【００３６】
混合比率が求まると、ＬＵＴ合成部Ｍ７は、特定された混合比率で既存ＬＵＴの各出力値におのおの重み付けして加算した値を出力値とするとともに入力値を既存ＬＵＴと同様にしたＬＵＴを作成し、これをミックス光源に対応するミックスＬＵＴとする。
図８は、既存ＬＵＴの合成によるミックスＬＵＴの作成を説明する図である。同図は、Ｄ５０光源とＡ光源を１：２の比率で混合してミックス光源を合成する場合を示してあり、入力表色系が一致しているＬＵＴ（Ｄ５０） sRGBとＬＵＴ（Ａ） sRGBを１：２の比率で混合してミックスＬＵＴを作成している。図８で作成されるミックスＬＵＴでは入力値に既存のＬＵＴ（Ｄ５０） sRGBやＬＵＴ（Ａ） sRGBと同様の値を利用しており、出力値は同じ入力値に対応している出力値同士を前記混合比率で重み付けした値を線形的に足し合わせた値としてある。このようにして作成されるミックスＬＵＴは、観察光源以外の作成条件が全て同一のＬＵＴ同士を合成して作成するため、出力値同士を単純に線形結合するだけで精度の高いＬＵＴを作成することができる。
【００３７】
ところで、以上のようなミックスＬＵＴを精度よく作成するためには、予め作成されているＬＵＴに対応する既存光源の色温度軸上での分布範囲やｘｙ色度図上での分布範囲が広い方が好ましいといえる。合成される既存ＬＵＴ同士の相違が大きいほど、例えば後述のＬＵＴ作成処理により作成した理想的なＬＵＴとミックスＬＵＴとの間の誤差を小さく出来るからである。また、予め作成されているＬＵＴに対応する既存光源が互いに大きく相違していると、広範な光源に対応したミックスＬＵＴを精度よく作成できるようになるからである。
【００３８】
例えば、色温度軸上での既存光源の分布範囲としては、色温度空間において既存光源が３０００Ｋ〜７０００Ｋの範囲の全域で離散的に分布するように選択することが好ましい。この範囲の色温度を再現可能であれば大方の観察光源が再現できるようになるからである。例えば、既存光源にＡ光源（色温度が２８５６Ｋ（約３０００Ｋ））やＣ光源（６７７４Ｋ（約７０００Ｋ））を含むように構成すると色温度のレンジが３０００Ｋ〜７０００Ｋの範囲の全域に離散的に分布するようになる。その他、Ｄ５０光源（５００３Ｋ）、Ｄ６５光源（６５０４Ｋ）、Ｆ２光源（４２３０Ｋ）、Ｆ１１光源（４０００Ｋ）、Ｆ７光源（６５００Ｋ）、Ｆ８光源（５０００Ｋ）、Ｄ５５光源（５５００Ｋ）、ＨＭＩ Ｌｉｇｈｔ４８００Ｋ（４８００Ｋ）、Ｔｕｎｇｓｔｅｎ ３２００Ｋ（３２００Ｋ）、Ｔｕｎｇｓｔｅｎ ３４００Ｋ（３４００Ｋ）等を適宜含むように構成すると、より精度の高いミックスＬＵＴを合成できるようになる。
【００３９】
以上のようにしてミックスＬＵＴが作成されると、色変換部Ｍ４は、ミックスＬＵＴを参照しつつ色変換を行ってインク量データを作成し、ハーフトーン処理部Ｍ５に出力する事になる。以上のＬＵＴ合成処理によれば、既存ＬＵＴを利用して新たな観察光源に対応するＬＵＴを作成できるため、ＬＵＴを記憶するリソースを増加することなく多様な観察光源に対応した印刷を行えるようになる。
【００４０】
（４）ＬＵＴ作成処理： 
（４−Ａ）．ＬＵＴ作成装置の構成と全体処理手順：
次に、コンピューター２０にあらかじめ記憶されている複数のＬＵＴの作成手法について説明する。以下の説明では、コンピューター２０に記憶されている複数のＬＵＴは、ＬＵＴ作成装置としてのコンピューター１０で作成される場合について説明する。むろん、プリンター３０でパッチを印刷して、このパッチを後述のステップＳ１００で指定された複数の観察光源下において測色機で測色し、同じパッチ同士の色の見えが一致するように複数のＬＵＴを作成してもよい。
【００４１】
図９は、コンピューター１０の構成を示すブロック図である。コンピューター１０は、ベースＬＵＴ作成モジュール１００と、色補正ＬＵＴ作成モジュール２００と、フォワードモデルコンバーター３００と、ＬＵＴ格納部４００と、を備えている。これらのモジュール１００，２００やコンバーター３００の機能は、メモリに格納されたコンピュータプログラムをコンピューターが実行することによってそれぞれ実現される。また、ＬＵＴ格納部４００は、ハードディスク装置などの記録媒体によって実現される。
【００４２】
ベースＬＵＴ作成モジュール１００は、スムージング処理初期値設定モジュール１２０と、テーブル作成モジュール１４０と、ＬＵＴ作成条件受付部１４５とを有している。スムージング処理モジュール１３０は、色点移動モジュール１３２と、インク量最適化モジュール１３４と、画質評価指数算出モジュール１３６とを有している。フォワードモデルコンバーター３００は、分光プリンティングモデルコンバーター３１０と、色算出部３２０とを備えている。これらの各部の機能については後述する。
【００４３】
ＬＵＴ格納部４００は、インバースモデル初期ＬＵＴ４１０や、ベース３Ｄ−ＬＵＴ５１０，ベース４Ｄ−ＬＵＴ５２０，色補正３Ｄ−ＬＵＴ６１０，色補正４Ｄ−ＬＵＴ６２０などを格納するためのものである。但し、インバースモデル初期ＬＵＴ４１０以外のＬＵＴは、ベースＬＵＴ作成モジュール１００や色補正ＬＵＴ作成モジュール２００によって作成されるものである。ベース３Ｄ−ＬＵＴ５１０は、ＲＧＢ表色系を入力とし、インク量を出力とする色変換ルックアップテーブルである。一方、ベース４Ｄ−ＬＵＴ５２０は、ＣＭＹＫ表色系を入力とし、インク量を出力とする色変換ルックアップテーブルである。なお、「３Ｄ」や「４Ｄ」は、入力値の数を意味している。これらのベースＬＵＴ５１０，５２０の入力表色系であるＲＧＢ表色系やＣＭＹＫ表色系は、いわゆる機器依存表色系では無く、特定のデバイスとは無関係に設定された仮想の表色系（あるいは抽象的な表色系）である。これらのベースＬＵＴ５１０，５２０は、例えば色補正ＬＵＴ６１０，６２０を作成する際に使用される。「ベースＬＵＴ」という名前は、色補正ＬＵＴを作成する基礎として用いられるからである。色補正ＬＵＴ６１０，６２０は、標準的な機器依存表色系（例えばｓＲＧＢ表色系やＪＡＰＡＮ  ＣＯＬＯＲ  ２００１表色系）を、特定のプリンターのインク量に変換するためのルックアップテーブルである。インバースモデル初期ＬＵＴ４１０については後述する。
【００４４】
図１０は、ＬＵＴ作成処理の全体処理手順を示すフローチャートである。図１１（Ａ）〜（Ｃ）は、図１０のステップＳ１２０〜Ｓ１４０によってベース３Ｄ−ＬＵＴを作成する場合の処理内容を示す説明図である。
処理が開始されると、ステップＳ１００で、ＬＵＴ作成者がベースＬＵＴ作成モジュール１００にＬＵＴ作成条件の指定を行う。そのために、ベースＬＵＴ作成モジュール１００は、ＬＵＴ作成条件受付部１４５を備えている。ＬＵＴ作成条件受付部１４５は、ＬＵＴ作成条件ＵＩをコンピューター１０の表示部１２に表示して操作入力部１３を介してユーザーからＬＵＴ作成条件の指定を受付ける。
【００４５】
図１２は、ＬＵＴ作成条件受付部１４５により表示部１２に表示されるＬＵＴ作成条件ＵＩの一例である。同図に示すように、ＬＵＴ作成条件ＵＩは、２以上の複数の観察光源を指定する観察光源選択欄Ｂ１と、ガマットマッピングの種類を指定するガマットマッピング指定択欄Ｂ２と、印刷用紙の種類を指定するための印刷用紙指定欄Ｂ３と、印刷品質を指定するための印刷品質指定欄Ｂ４とを備えている。観察光源選択欄Ｂ１は、コンピューター２０で印刷された印刷結果をユーザーが観察するときの観察条件としての観察光源（Ｄ５０，Ａ，Ｃ等）を指定する欄である。ガマットマッピング指定欄Ｂ２は、入力画像データの表色系と標準的な表示部２３の機器依存表色系との間で行うガマットマッピングの手法を選択する欄である。ガマットマッピングは、公知の各種の手法（「明度重視」、「彩度重視」等）を採用可能である。
【００４６】
ＬＵＴ作成条件の指定が終了してＬＵＴ作成実行ボタンＢ５が選択されると、ＬＵＴ作成条件受付部１４５は、受付けたベースＬＵＴ作成モジュール１００に出力する。すると、ベースＬＵＴ作成モジュール１００は、観察光源の種類やガマットマッピングの種類、印刷用紙の種類、印刷品質の種類、の全ての組合せでそれぞれＬＵＴが作成されるようにステップＳ１２０〜Ｓ１５０の処理を繰り返し実行することになる（Ｓ１１０）。全ての組合せについてＬＵＴを作成するとＬＵＴ作成処理を終了し（Ｓ１１０：Ｙｅｓ）、未作成の組合せがあればステップＳ１２０〜Ｓ１５０のＬＵＴ作成処理を繰り返す（Ｓ１１０：Ｎｏ）。
【００４７】
ステップＳ１２０では、フォワードモデルコンバーター３００とインバースモデル初期ＬＵＴ４１０とが準備される。ここで、「フォワードモデル」とは、インク量を機器独立表色系の色彩値（測色値）に変換する変換モデルを意味し、「インバースモデル」とは、逆に、機器独立表色系の色彩値をインク量に変換する変換モデルを意味している。実施例では、機器独立表色系としてＣＩＥ−Ｌａｂ表色系を使用する。なお、以下では、ＣＩＥ−Ｌａｂ表色系の色彩値を、単に「Ｌ*ａ*ｂ*値」または「Ｌａｂ値」とも呼ぶ。
【００４８】
図１１（Ａ）に示すように、フォワードモデルコンバーター３００の前段を構成する分光プリンティングモデルコンバーター３１０は、複数種類のインクのインク量を、そのインク量に応じて印刷されるカラーパッチの分光反射率Ｒ（λ）に変換する。なお、本明細書において「カラーパッチ」という用語は、有彩色のパッチに限らず、無彩色のパッチも含む広い意味で使用される。本実施例では、シアン（Ｃ）、マゼンタ（Ｍ）、イエロー（Ｙ）、ブラック（Ｋ）、淡シアン（Ｌｃ）、淡マゼンタ（Ｌｍ）、淡ブラック（Ｌｋ）、極淡ブラック（ＬＬｋ）の８種類のインクを利用可能なカラープリンターを想定しており、分光プリンティングモデルコンバーター３１０もこの８種類のインクのインク量を入力としている。但し、プリンターで使用する複数種類のインクとしては、任意のインクセットを利用することが可能である。色算出部３２０は、分光反射率Ｒ（λ）からＬａｂ表色系の色彩値を算出する。この色彩値の算出には、ステップＳ１００においてＬＵＴ作成者によって選択された光源がカラーパッチの観察条件として使用される。なお、分光プリンティングモデルコンバーター３１０を作成する方法としては、例えば特表２００７−５１１１７５号公報に記載された方法を採用することが可能である。
【００４９】
インバースモデル初期ＬＵＴ４１０は、Ｌ*ａ*ｂ*値を入力とし、インク量を出力とするルックアップテーブルである。この初期ＬＵＴ４１０は、例えば、Ｌ*ａ*ｂ*空間を複数の小セルに区分し、各小セル毎に最適なインク量を選択して登録したものである。この選択は、例えば、そのインク量で印刷されるカラーパッチの画質を考慮して行われる。一般に、或る１つのＬ*ａ*ｂ*値を再現するインク量の組合せは多数存在する。そこで、初期ＬＵＴ４１０では、ほぼ同じＬ*ａ*ｂ*値を再現する多数のインク量の組合せの中から、画質等の所望の観点から最適なインク量を選択したものが登録されている。この初期ＬＵＴ４１０の入力値であるＬ*ａ*ｂ*値は各小セルの代表値である。一方、出力値であるインク量はそのセル内のいずれかのＬ*ａ*ｂ*値を再現するものである。従って、この初期ＬＵＴ４１０では、入力値であるＬ*ａ*ｂ*値と出力値であるインク量とが厳密に対応したものとなっておらず、出力値のインク量をフォワードモデルコンバーター３００でＬ*ａ*ｂ*値に変換すると、初期ＬＵＴ４１０の入力値とは多少異なる値が得られる。但し、初期ＬＵＴ４１０として、入力値と出力値とが完全に対応するものを利用してもよい。また、初期ＬＵＴ４１０を用いずにベースＬＵＴを作成することも可能である。なお、小セル毎に最適なインク量を選択して初期ＬＵＴ４１０を作成する方法としては、例えば前記特表２００７−５１１１７５号公報に記載された方法を採用することが可能である。
【００５０】
図１０のステップＳ１３０では、ベースＬＵＴ作成のための初期入力値がユーザーによって設定される。図１１（Ｂ）は、ベース３Ｄ−ＬＵＴ５１０の構成とその初期入力値設定の例を示している。ベース３Ｄ−ＬＵＴ５１０の入力値としては、ＲＧＢの各値として予め定められたほぼ等間隔の値が設定される。１組のＲＧＢ値はＲＧＢ色空間内の点を表していると考えられるので、１組のＲＧＢ値を「入力格子点」とも呼ぶ。ステップＳ２００においては、複数の入力格子点のうちから予め選択されたいくつかの少数の入力格子点に対するインク量の初期値がユーザーによって入力される。この初期入力値が設定される入力格子点としては、ＲＧＢ色空間における３次元色立体の頂点に相当する入力格子点を少なくとも選択することが好ましい。この３次元色立体の頂点では、ＲＧＢの各値がその定義範囲の最小値または最大値を取る。具体的には、ＲＧＢの各値を８ビットで表現した場合には、（Ｒ，Ｇ，Ｂ）＝（０，０，０）、（０，０，２５５）、（０，２５５，０）、（２５５，０，０）、（０，２５５，２５５）、（２５５，０，２５５）、（２５５，２５５，０）、（２５５，２５５，２５５）である８つの入力格子点に関してインク量の初期入力値が設定される。なお、（Ｒ，Ｇ，Ｂ）＝（２５５，２５５，２５５）の入力格子点に対するインク量は、すべて０に設定される。他の入力格子点に対するインク量の初期入力値は任意であり、例えば０に設定される。図１１（Ｂ）の例では、（Ｒ，Ｇ，Ｂ）＝（０，０，３２）の入力格子点に対するインク量が０以外の値になっているが、これはこのＬＵＴ５１０が完成したときの値である。
【００５１】
図１０のステップＳ１４０では、スムージング処理モジュール１３０（図９）が、ステップＳ１２０で設定された初期入力値に基づいてスムージング処理（平滑化／最適化処理）を実行する。図１１（Ｃ）は、ステップＳ１４０の処理内容を示している。図１１（Ｃ）の左側には、スムージング処理前の状態における複数の色点の分布が２重丸と白丸とで示されている。これらの色点は、Ｌ*ａ*ｂ*空間における３次元色立体ＣＳを構成している。各色点のＬ*ａ*ｂ*座標値は、ベース３Ｄ−ＬＵＴ５１０の複数の入力格子点におけるインク量を、フォワードモデルコンバーター３００（図１１（Ａ））を用いてＬ*ａ*ｂ*値に変換した値である。上述したように、ステップＳ１３０では一部の少数の入力格子点についてのみインク量の初期入力値が設定される。そこで、他の入力格子点に対するインク量の初期値は、初期入力値からスムージング処理初期値設定モジュール１２０（図９）によって設定される。この初期値設定方法については後述する。
【００５２】
Ｌａｂ表色系の３次元色立体ＣＳは、以下の８つの頂点（図１１（Ｃ）の２重丸の点）を有している。
・点ＰK：（Ｒ，Ｇ，Ｂ）＝（０，０，０）に対応する紙黒点。
・点ＰW：（Ｒ，Ｇ，Ｂ）＝（２５５，２５５，２５５）に対応する紙白点。
・点ＰC：（Ｒ，Ｇ，Ｂ）＝（０，２５５，２５５）に対応するシアン点。
・点ＰM：（Ｒ，Ｇ，Ｂ）＝（２５５，０，２５５）に対応するマゼンタ点。
・点ＰY：（Ｒ，Ｇ，Ｂ）＝（２５５，２５５，０）に対応するイエロー点。
・点ＰR：（Ｒ，Ｇ，Ｂ）＝（２５５，０，０）に対応するレッド点。
・点ＰG：（Ｒ，Ｇ，Ｂ）＝（０，２５５，０）に対応するグリーン点。
・点ＰB：（Ｒ，Ｇ，Ｂ）＝（０，０，２５５）に対応するブルー点。
【００５３】
図１１（Ｃ）の右側は、スムージング処理後の色点の分布を示している。スムージング処理は、Ｌ*ａ*ｂ*空間における複数の色点を移動させて、それらの色点の分布を等間隔に近い平滑なものにする処理である。スムージング処理では、さらに、移動後の各色点のＬ*ａ*ｂ*値を再現するために最適なインク量も決定される。この最適なインク量がベースＬＵＴ５１０の出力値として登録されると、ベースＬＵＴ５１０が完成する。
【００５４】
図１３（Ａ）〜（Ｃ）は、入力表色系の色点（すなわち入力格子点）とＬａｂ表色系の色点との対応関係を示している。Ｌａｂ表色系の３次元色立体ＣＳの頂点は、ベースＬＵＴ５１０の入力表色系の３次元色立体の頂点と一対一に対応している。また、各頂点を結ぶ辺（稜線）も、両方の色立体で互いに対応しているものと考えることができる。スムージング処理前のＬａｂ表色系の各色点は、ベースＬＵＴ５１０の入力格子点にそれぞれ対応付けられており、従って、スムージング処理後のＬａｂ表色系の各色点もベースＬＵＴ５１０の入力格子点にそれぞれ対応付けられる。なお、ベースＬＵＴ５１０の入力格子点はスムージング処理によって変化しない。スムージング処理後のＬａｂ表色系の３次元色立体ＣＳは、ベースＬＵＴ５１０の出力表色系を構成するインクセットで再現可能な色域（カラーガマット）の全体に対応している。従って、ベースＬＵＴ５１０の入力表色系は、このインクセットで再現可能な色域の全体を表す表色系としての意義を有している。
【００５５】
ベースＬＵＴ５１０を作成する際に、Ｌ*ａ*ｂ*空間においてスムージング処理を行う理由は以下の通りである。ベースＬＵＴ５１０では、なるべく大きな色域を再現できるように出力表色系のインク量を設定したいという要望がある。特定のインクセットで再現可能な色域は、インクデューティ制限（一定面積に吐出可能なインク量の制限）などの所定の制限条件を考慮して決定される。一方、上述したフォワードモデルコンバーター３００は、これらの制限条件が考慮されておらず、再現可能な色域とは無関係に作成されている。そこで、スムージング処理の際にインクデューティ制限等の制限条件を考慮してＬ*ａ*ｂ*空間内の色点の取り得る範囲を決定すれば、特定のインクセットで再現可能な色域を決定することが可能となる。なお、色点の移動を行うアルゴリズムとしては、例えば、後述する力学モデルを使用したものが利用される。
【００５６】
図１０のステップＳ４００では、スムージング処理の結果を用いて、テーブル作成モジュール１４０がベースＬＵＴ５１０を作成する。すなわち、テーブル作成モジュール１４０は、各入力格子点に対応付けられたＬａｂ表色系の色点を再現するための最適なインク量をベースＬＵＴ５１０（図１１（Ｂ））の出力値として登録する。なお、スムージング処理では、その計算負荷を軽減するために、ベースＬＵＴ５１０の入力格子点の一部のみに対応する色点のみを処理対象として選択することも可能である。例えば、ベースＬＵＴ５１０の入力格子点におけるＲＧＢ値の間隔が１６である場合に、スムージング処理の対象となる入力格子点におけるＲＧＢ値の間隔を３２に設定すれば、スムージング処理の負荷を半減することができる。この場合には、テーブル作成モジュール１４０は、スムージング処理結果を補間することによってベースＬＵＴ５１０のすべての入力格子点に対するインク量を決定して登録する。
【００５７】
図１４（Ａ）〜（Ｃ）は、図１０のステップＳ１２０〜Ｓ１４０によってベース４Ｄ−ＬＵＴ５２０を作成する場合の処理内容を示す説明図である。図１４（Ａ）は、図１１（Ａ）と同じである。図１４（Ｂ）に示すベース４Ｄ−ＬＵＴ５２０は、入力がＣＭＹＫ表色系である点が図１１（Ｂ）に示したベース３Ｄ−ＬＵＴ５１０と異なっている。このベース４Ｄ−ＬＵＴ５２０の初期入力値としては、（Ｃ，Ｍ，Ｙ，Ｋ）＝（０，０，０，０），（０，０，２５５，０），（０，２５５，０，０），（０，２５５，２５５，０），（２５５，０，０，０），（２５５，０，２５５，０），（２５５，２５５，０，０），（２５５，２５５，２５５，０），（０，０，０，２５５），（０，０，２５５，２５５），（０，２５５，０，２５５），（０，２５５，２５５，２５５），（２５５，０，０，２５５），（２５５，０，２５５，２５５），（２５５，２５５，０，２５５），（２５５，２５５，２５５，２５５）である１６個の入力格子点に関してインク量の初期値が設定される。他の入力格子点に対するインク量の初期入力値は任意であり、例えば０に設定される。
【００５８】
図１４（Ｃ）は、スムージング処理の様子を示している。なお、Ｌ*ａ*ｂ*空間においてベース４Ｄ−ＬＵＴ５２０に対応する色立体としては、図１４（Ｃ）の右端に示すように、入力値のうちのＫ値のそれぞれの値に対して１つの３次元色立体ＣＳが存在する。この例では、Ｋ＝０の色立体とＫ＝３２の色立体とを含む複数の色立体ＣＳが図示されている。本明細書では、これらの個々の色立体ＣＳを「Ｋレイヤ」とも呼ぶ。この理由は、各色立体ＣＳが、ＣＭＹＫ値のうちのＫ値が一定でＣ，Ｍ，Ｙ値が可変である入力層に対応するものと考えることができるからである。複数の色立体ＣＳは、Ｋ値が大きいほど暗い色域を表現するものとなっている。これらの複数の色立体ＣＳは、入力表色系のＫ値が大きいほど濃ブラックインクＫのインク量が多くなるように濃ブラックインクＫのインク量を決定することによって実現できる。上述したように、再現可能な色域はインクデューティ制限値によって制限される。このインクデューティ制限値としては、個々のインクのインク量と、全インクの合計インク量と、の２つの制限値が課されるのが普通である。一方、暗い色を再現する方法としては、濃ブラックインクＫなどの無彩色インクを用いる方法と、コンポジットブラックを用いる方法とがある。しかし、コンポジットブラックは合計インク量が多くなるので、濃ブラックインクＫに比べてインクデューティ制限値に抵触する可能性が高く、暗い色を再現するのには不利である。従って、入力表色系のＫ値が大きく濃ブラックインクＫのインク量が多い色立体の方が、入力表色系のＫ値が小さく濃ブラックインクＫのインク量が少ない色立体に比べてより暗い色を再現することが可能となる。
【００５９】
図１５（Ａ），（Ｂ）は、ベースＬＵＴを用いた色補正ＬＵＴの作成方法を示す説明図である。図１５（Ａ）に示すように、ベース３Ｄ−ＬＵＴ５１０は、ＲＧＢ値をインク量Ｉjに変換する。このインク量Ｉjは、図１１（Ｂ）に示した８種類のインクのインク量を表している。このとき、インク量Ｉjの添え字ｊは１〜８である。変換後のインク量Ｉjは、フォワードモデルコンバーター３００によってＬ*ａ*ｂ*値に変換される。一方、ｓＲＧＢ値は、既知の変換式に従ってＬ*ａ*ｂ*値に変換される。この変換後のＬ*ａ*ｂ*値は、その色域が、フォワードモデルコンバーター３００で変換されたＬ*ａ*ｂ*値の色域と一致するようにガマットマッピングされる。このガマットマッピングの手法は、図１１のＬＵＴ作成条件ＵＩで指定されたガマットマッピングの手法である。一方、ベース３Ｄ−ＬＵＴ５１０とフォワードモデルコンバーター３００を通じて、ＲＧＢ値から変換したＬ*ａ*ｂ*値を、逆方向ルックアップテーブルとして、逆変換ＬＵＴ５１１を作成する。ガマットマッピングされたＬ*ａ*ｂ*値は、この逆変換ＬＵＴ５１１によってＲＧＢ値に変換される。このＲＧＢ値は、さらに、ベース３Ｄ−ＬＵＴ５１０によってインク量Ｉjに再度変換される。この最後のインク量Ｉjと最初のｓＲＧＢ値の対応関係をルックアップテーブルに登録することによって、色補正３Ｄ−ＬＵＴ６１０を作成することができる。この色補正３Ｄ−ＬＵＴ６１０は、ｓＲＧＢ表色系をインク表色系に変換するＬＵＴである。
【００６０】
図１５（Ｂ）は、色補正４Ｄ−ＬＵＴ６２０を作成する方法を示している。図１５（Ａ）との違いは、３Ｄ−ＬＵＴ５１０およびその逆変換ＬＵＴ５１１の代わりに、４Ｄ−ＬＵＴ５２０およびその逆変換ＬＵＴ５２１を利用している点と、ｓＲＧＢ表色系をＬ*ａ*ｂ*値に変換する既知変換式の代わりにＪＡＰＡＮ  ＣＯＬＯＲ表色系（図中では「ｊＣＭＹＫ」と記したもの）をＬ*ａ*ｂ*値に変換する既知変換式を使用している点である。良く知られているように、ＪＡＰＡＮ  ＣＯＬＯＲは、ＣＭＹＫの４色で構成される表色系である。なお、図１５（Ｂ）の方法では、逆変換ＬＵＴ５２１において、Ｌ*ａ*ｂ*値からＣＭＹＫ値に変換する際に、既知変換前の最初のｊＣＭＹＫ値のＫ値から、逆変換ＬＵＴ５２１のＫレイヤ（Ｋ値が一定を取る部分）が選択される。従って、色補正４Ｄ−ＬＵＴ６２０として、ベース４Ｄ−ＬＵＴ５２０のうちのＫレイヤにおける特性を反映したものを作成することが可能である。
【００６１】
なお、通常は、ベースＬＵＴ５１０，５２０がプリンタドライバーに実装されており、色補正ＬＵＴの作成処理以外の処理にも活用されているが、ここでは他の活用例の説明は省略する。以下では、実施例のスムージング処理（平滑化／最適化処理）に利用される力学モデルについて簡単に説明した後に、スムージング処理の処理手順、および、最適化処理の内容について順次説明する。
【００６２】
（４−Ｂ）力学モデル：
図１６は、本実施例のスムージング処理（平滑化／最適化処理）に利用される力学モデルを示す説明図である。ここでは、Ｌ*ａ*ｂ*色空間内に複数の色点（白丸および２重丸）が配列されている様子を示している。ただし、ここでは説明の便宜上、色点の配置を２次元的に描いている。この力学モデルでは、着目色点ｇに対して次式の仮想的な力Ｆpgが係るものと仮定する。
【数５】
2011139213.tif 000007
ここで、Ｆgは着目色点ｇが隣接色点ｇｎ（ｎは１〜Ｎ）から受ける引力の合計値、Ｖgは着目色点ｇの速度ベクトル、−ｋvＶgは速度に応じた抵抗力、Ｘgは着目色点ｇの位置ベクトル、Ｘgnは隣接色点ｇｎの位置ベクトル、ｋp，ｋgは係数である。係数ｋp，ｋgは予め一定の値に設定される。なお、文中では、ベクトルを示す矢印は省略される。
【００６３】
このモデルは、バネで互いに結ばれた質点の減衰振動モデルである。すなわち、着目色点ｇに係る仮想合力Ｆpgは、着目色点ｇと隣接色点ｇｎとの距離が大きいほど大きくなるバネ力Ｆgと、着目色点ｇの速度が大きいほど大きくなる抵抗力−ｋvＶgとの合計値である。この力学モデルでは、各色点について、位置ベクトルＸgと速度ベクトルＦgの初期値を設定した後に、微小時間経過の位置ベクトルＸgと速度ベクトルＦgを順次算出してゆく。なお、複数の色点の速度ベクトルＶgの初期値は、例えば０に設定される。このような力学モデルを利用すれば、各色点が徐々に移動して、平滑な色点分布を得ることが可能である。
【００６４】
なお、各色点に係る力としては、バネ力Ｆgと抵抗力−ｋvＶg以外の力を用いても良い。例えば、本出願人により開示された特開２００６−１９７０８０号公報で説明されている他の種々の力をこの力学モデルで利用してもよい。また、力学モデルを適用して各色点を移動させる際に、特定の色点は、力学モデルによって移動しない拘束点として取り扱うことも可能である。
【００６５】
（４−Ｃ）スムージング処理（平滑化／最適化処理）の処理手順：
図１７は、スムージング処理（図１０のステップＳ１４０）の典型的な処理手順を示すフローチャートである。ステップＴ１００では、スムージング処理初期値設定モジュール１２０（図９）が、スムージング処理の対象とする複数の色点を初期設定する。
【００６６】
図１８は、ステップＴ１００の詳細手順を示すフローチャートである。ステップＴ１０２では、インク量の初期入力値（図１１（Ｂ），図１４（Ｂ））から、スムージング処理の対象となる各色点の仮インク量が決定される。例えば、３Ｄ−ＬＵＴ用のスムージング処理では、次の（６）式、（７）式に従って、各入力格子点に対する仮インク量Ｉ(R,G,B)が決定される。
【数６】
2011139213.tif 000008
【数７】
2011139213.tif 000009
ここで、Ｉ(R,G,B)は、入力格子点のＲＧＢ値に対するインクセット全体のインク量（図１０の例では８種類のインクのインク量）を表している。ＲＧＢ値が０または２５５を取る入力格子点に対するインク量は、図１０のステップＳ１３０においてユーザーによって予め入力された値である。前記（６）式および（７）式によれば、任意のＲＧＢ値における仮インク量Ｉ(R,G,B)を求めることが可能である。
【００６７】
４Ｄ−ＬＵＴ用のスムージング処理では、次の（８）式、（９）式に従って、各入力格子点に対する仮インク量Ｉ(C,M,Y,K)が決定される。
【数８】
2011139213.tif 000010
【数９】
2011139213.tif 000011
【００６８】
なお、（８）式からも理解できるように、４Ｄ−ＬＵＴ用のインク量の初期入力値は１６個存在するので、初期入力値の設定が煩雑である。そこで、例えば、インク量の初期入力値を設定する入力格子点を、Ｋ＝０の８個の頂点、すなわち、（Ｃ，Ｍ，Ｙ，Ｋ）＝（０，０，０，０），（０，０，２５５，０），（０，２５５，０，０），（０，２５５，２５５，０），（２５５，０，０，０），（２５５，０，２５５，０），（２５５，２５５，０，０），（２５５，２５５，２５５，０）の８個の頂点と、Ｋ＝２５５の１個の頂点、例えば、（Ｃ，Ｍ，Ｙ，Ｋ）＝（０，０，０，２５５）の頂点のみとし、Ｋ＝２５５の色点のインク量を次の（１０）式または（１１）式で決定してもよい。
【数１０】
2011139213.tif 000012
【数１１】
2011139213.tif 000013
ここで、Ｉ(C,M,Y,0)は、Ｋ＝０の８個の頂点におけるインク量の初期入力値から、前記（６）式と同様の式で算出されたインク量である。（１０）式の関数ｆ D1は値Ｉ(C,M,Y,0)と値Ｉ(0,0,0,255)の合計値がインクデューティ制限値をオーバーする場合に、値Ｉ(C,M,Y,0)を減じることによって、インク量Ｉ(C,M,Y,255)がインクデューティ制限値内に納まるようにする関数である。また、（１１）式の関数ｆ D2は、値Ｉ(C,M,Y,0)と値Ｉ(0,0,0,255)の合計値がインクデューティ制限値をオーバーする場合に、合計値（Ｉ(C,M,Y,0)＋Ｉ(0,0,0,255)）の全体を減じることによって、インク量Ｉ(C,M,Y,255)がインクデューティ制限値内に納まるようにする関数である。なお、インク量Ｉj（Ｉ j(R,G,B)，ΔＩj，Ｉjr，ｈjも含む）を下付文字ｊを付すことなく示す場合は、各インクのインク量Ｉjを各行要素として有する行列（ベクトル）を意味することとする。
【００６９】
図１８のステップＴ１０４では、フォワードモデルコンバーター３００を用いて、仮インク量に対応する色彩値Ｌ*ａ*ｂ*を求める。この演算は、以下の（１２）式または（１３）式で表すことができる。
【数１２】
2011139213.tif 000014
【数１３】
2011139213.tif 000015
ここで、Ｌ* (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はフォワードモデルコンバーター３００による変換を意味している。なお、これらの式からも理解できるように、この変換後の色彩値Ｌ*ａ*ｂ*は、ベースＬＵＴの入力値であるＲＧＢ値またはＣＭＹＫ値に対応付けられている。
【００７０】
図１８のステップＴ１０６では、ステップＴ１０４で得られた色彩値Ｌ*ａ*ｂ*を、インバースモデル初期ＬＵＴ４１０（図１１（Ａ））を用いてインク量に再度変換する。ここで、インバースモデル初期ＬＵＴ４１０を用いてインク量に再度変換する理由は、インク量の初期入力値や、ステップＴ１０２で決定された仮インク量が、Ｌ*ａ*ｂ*値を再現するインク量として必ずしも好ましいインク量では無いからである。一方、インバースモデル初期ＬＵＴ４１０では、画質等を考慮した好ましいインク量が登録されているので、これを用いてＬ*ａ*ｂ*値をインク量に再度変換すれば、そのＬ*ａ*ｂ*値を実現するための好ましいインク量を初期値として得ることができる。但し、ステップＴ１０６を省略してもよい。
【００７１】
上述のステップＴ１００の処理の結果、スムージング処理の対象となる色点について、以下の初期値が決定される。
（１）ベースＬＵＴの入力格子点の値：（Ｒ，Ｇ，Ｂ）または（Ｃ，Ｍ，Ｙ，Ｋ）
（２）各入力格子点に対応するＬ*ａ*ｂ*空間の色点の初期座標値：（Ｌ* (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に対するインク量の最適値を探索する（図１９（Ｃ）参照）。この目的関数Ｅを用いた最適化では、力学モデルで微小量だけ移動した後の色点の座標値ＬＡＢtに近いＬ*ａ*ｂ*値を再現するインク量Ｉjの中で、複数のパラメータΔＬ*，Δａ*，ΔＧＩ，ΔＣＩＩ，ΔＴＩの２乗誤差の和がなるべく小さいインク量が最適なインク量として決定される。また、最適なインク量の探索は、ステップＴ１００で設定された各入力格子点の初期インク値から開始される。従って、探索で得られるインク量は、この初期インク量を修正した値となる。後で詳述するように、（ＥＱ１）式で与えられる目的関数Ｅは、（ＥＱ２）式のようなインク量ベクトルＩに関する２次形式の関数として書き表すことができる。インク量の最適化は、このような２次形式の目的関数Ｅを用いて、２次計画法に従って実行される。なお、ステップＴ１３０の詳細手順や目的関数Ｅの内容については後述する。
【００７５】
図１７のステップＴ１４０では、ステップＴ１３０で探索されたインク量Ｉjに対応するＬ*ａ*ｂ*値が、フォワードモデルコンバーター３００で再計算される（図１９（Ｄ）参照）。ここで、Ｌ*ａ*ｂ*値を再計算する理由は、探索されたインク量Ｉjが目的関数Ｅを最小とするインク量なので、そのインク量Ｉjで再現されるＬ*ａ*ｂ*値は、最適化処理のターゲット値ＬＡＢtから多少ずれているからである。こうして再計算されたＬ*ａ*ｂ*値が、各色点の移動後の座標値として使用される。
【００７６】
ステップＴ１５０では、各色点の座標値の移動量の平均値（ΔＬａｂ） aveが、予め設定された閾値ε以下であるか否かが判定される。移動量の平均値（ΔＬａｂ） aveが閾値εよりも大きい場合には、ステップＴ１２０に戻りステップＴ１２０〜Ｔ１５０のスムージング処理が継続される。一方、移動量の平均値（ΔＬａｂ） aveが閾値ε以下の場合には、色点の分布が十分に平滑になっているので、スムージング処理が終了する。なお、閾値εは、予め適切な値が実験的に決定される。
【００７７】
このように、本実施例の典型的なスムージング処理（平滑化／最適化処理）では、力学モデルによって各色点を微小時間間隔毎に移動させつつ、移動後の色点に対応する最適なインク量を最適化手法で探索する。そして、色点の移動量が十分に小さくなるまでそれらの処理が継続される。この結果、図１１（Ｃ）に示したように、スムージング処理によって、平滑な色点分布を得ることが可能である。
【００７８】
（４−Ｄ）最適化処理の内容：
最適化処理の目的関数Ｅ（図１９（Ｃ）参照）は、インク量の関数である色彩値（Ｌ*ａ*ｂ*値）および画質評価指数に関するヤコビ行列Ｊを用いて表現することが可能である。ヤコビ行列Ｊは、例えば以下の（１４）式で表される。
【００７９】
【数１４】
2011139213.tif 000016
【００８０】
（１４）式の右辺の第１行〜第３行は、色彩値Ｌ*ａ*ｂ*を個々のインク量Ｉjで偏微分した値を示している。また、第４行以下は、１組のインク量Ｉj（ｊ＝１〜８）で印刷されるカラーパッチの画質を表す画質評価指数（粒状性指数ＧＩ（Graininess Ｉndex）と、色非恒常性指数ＣＩＩ（Color Inconstancy Index）と、合計インク量ＴＩを個々のインク量Ｉjで偏微分した値を示している。なお、画質評価指数ＧＩ，ＣＩＩ，ＴＩは、その値が小さいほど、インク量Ｉjで再現されるカラーパッチの画質が良い傾向にあることを示す指数である。
【００８１】
色彩値Ｌ*ａ*ｂ*は、フォワードモデルコンバーター３００を用いて、以下の（１５）式でインク量Ｉjから変換される。
【数１５】
2011139213.tif 000017
【００８２】
画質評価指数ＧＩ，ＣＩＩ，ＴＩも、一般にインク量Ｉjの関数としてそれぞれ表現できる。
【数１６】
2011139213.tif 000018
【数１７】
2011139213.tif 000019
【数１８】
2011139213.tif 000020
【００８３】
なお、（１７）式の色非恒常性指数ＣＩＩillの下付文字「ill」は、光源の種類を表している。上述した（１４）式では、光源の種類として、標準の光Ａと標準の光Ｆ１２とを用いている。なお、色非恒常性指数ＣＩＩの計算方法の例は後述するが、色非恒常性指数ＣＩＩとしては一種類または複数種類の任意の標準光源に関するものを利用することが可能である。
【００８４】
粒状性指数ＧＩは、各種の粒状性予測モデルを用いて算出可能であり、例えば以下の（１９）式で算出することができる。
【数１９】
2011139213.tif 000021
ここで、ａＬは明度補正係数、ＷＳ（ｕ）はカラーパッチの印刷に利用されるハーフトーンデータが示す画像のウイナースペクトラム、ＶＴＦ（ｕ）は視覚の空間周波数特性、ｕは空間周波数である。ハーフトーンデータは、カラーパッチのインク量Ｉjからハーフトーン処理（コンピューター２０やプリンター３０が実行するハーフトーン処理と同一のものとする）によって決定される。前記（１９）式は一次元で表現しているが、空間周波数の関数として二次元画像の空間周波数を算出することは容易である。粒状性指数ＧＩの計算方法としては、例えば、本出願人により開示された特表２００７−５１１１６１号公報に記載された方法や、Makoto Fujino， Ｉmage Quality Evaluation of Inkjet Prints， Japan Hardcopy '99， p.291-294に記載された方法を利用することができる。
【００８５】
色非恒常性指数ＣＩＩは、例えば以下の（２０）式で与えられる。
【数２０】
2011139213.tif 000022
ここで、ΔＬ*は２つの異なる観察条件下（異なる光源下）におけるカラーパッチの明度差、ΔＣ*abは彩度差、ΔＨ*abは色相差を示す。色非恒常性指数ＣＩＩの計算時には、２つの異なる観察条件下でのＬ*ａ*ｂ*値は、色順応変換（ＣＡＴ）を用いて標準観察条件（例えば標準の光Ｄ６５の観察下）に変換される。ＣＩＩについては、Billmeyer and Saltzman's Principles of Color Technology， 3rd edition， John Wiley & Sons， Ｉnc， 2000， p.129， pp. 213-215を参照。
【００８６】
ヤコビ行列Ｊの複数の成分（「要素」とも呼ぶ）のうち、例えばＬ*値に関する成分は、（２１）式で与えられる。
【数２１】
2011139213.tif 000023
ここで、ｆL*FMは、フォワードモデルによるインク量ＩからＬ*値への変換関数、Ｉrはインク量Ｉの現在値（平滑化／最適化処理前のインク量）、ｈjはｊ番目のインク量Ｉjの微小変動量である。ヤコビ行列Ｊの最下行を除く他の成分も同じ形式で表される。前記（１８）、（２１）式に準じて、ヤコビ行列Ｊの最下行の要素を算出すると、ヤコビ行列Ｊの最下行の要素はすべて１となる。あるインクのインク量Ｉjが微小変動量ｈjだけ変動した場合の合計インク量ＴＩの変動量もｈjとなるからである。しかしながら、本実施例では、（２２）式に示すように、ヤコビ行列Ｊの最下行の要素を１とせず、調整係数Ｕjを使用することとする。
【数２２】
2011139213.tif 000024
前記（２２）式右辺は、合計インク量ＴＩに関する成分を除いた本来のヤコビ行列の最下行に調整係数Ｕjの行を追加したものであると捉えることができる。また、前記（２２）式右辺は、前記（１８）式に示した現実の合計インク量ＴＩではなく、下記（２３）式に示す仮想的な合計インク量ＴＩについての偏微分要素を最下行に有するヤコビ行列Ｊであるとも捉えることができる。本明細書では、前記（２２）式の右辺全体をヤコビ行列Ｊと称呼するものとする。
【数２３】
2011139213.tif 000025
【００８７】
図２０（Ａ），（Ｂ）は、調整係数Ｕjの一例を模式的に示している。図２０（Ａ）はＣＩＥ−Ｌａｂ色空間を示しており、ＣＩＥ−Ｌａｂ色空間においてターゲット値Ｌt*ａt*ｂt*が示す位置を２重丸で示している。なお、ターゲット値Ｌt*ａt*ｂt*は、上述した力学モデルに基づくものである。ＣＩＥ−Ｌａｂ色空間を概ね上下（高明度側と低明度側）に分断する基準面Ｂ１，Ｂ２が設けられており、ターゲット値Ｌt*ａt*ｂt*を通過しＬ*軸と平行な直線Ｌと、基準面Ｂ１，Ｂ２との交点がそれぞれ制御点ＣＰ１，ＣＰ２として白丸で示されている。図２０（Ｂ）はターゲット値Ｌt*ａt*ｂt*（２重丸）と制御点ＣＰ１，ＣＰ２（白丸）を通過し、Ｌ*軸と平行な直線に沿った調整係数Ｕj（ｊ＝１，５：Ｃインク，Ｌｃインク）の変動を示している。
【００８８】
Ｃインクについての調整係数Ｕ1は制御点ＣＰ１よりも低明度の領域において一定の正値Ｕ c1となり、制御点ＣＰ１よりも高明度の領域において明度が大きくなるほど正値Ｕ c1から線形的に減少していく。そして、所定の明度よりも高明度の領域では、調整係数Ｕ1は一定の負値Ｕ c2となる。一方、Ｌｃインクについての調整係数Ｕ5は制御点ＣＰ２よりも高明度の領域において一定の正値Ｕ c1となり、制御点ＣＰ２よりも低明度の領域において明度が小さくなるほど正値Ｕ c1から線形的に減少していく。そして、調整係数Ｕ5が所定の明度よりも低明度の領域では、調整係数Ｕ5は一定の負値Ｕ c2となる。図２０（Ｂ）では、ＣインクとＬｃインクを例に挙げたが、調整係数Ｕjは、濃インクについては低明度域で淡インクよりも大きい値とされ、淡インクについては高明度域で濃インクよりも大きい値とされる。なお、基準面Ｂ１，Ｂ２は、Ｌ*軸に対して垂直ではないため、ターゲット値Ｌt*ａt*ｂt*が示す色相・彩度に応じて制御点ＣＰ１，ＣＰ２の明度が変動する。図２０（Ａ）の例では、基準面Ｂ１，Ｂ２を互いに平行としたが、これらは互いに平行であるものに限られない。基準面Ｂ１，Ｂ２を非平行とした場合、色相・彩度に応じて制御点ＣＰ１，ＣＰ２の明度の大小関係が反転することもある。
【００８９】
最適化の目的関数Ｅは、例えば以下の（２４）式で与えられる。
【数２４】
2011139213.tif 000026
ここで、右辺の各項の最初に記載されているｗ L*，ｗ a*等は、各項の重みである。各項の重みｗ L*、ｗ a*…は、予め設定されている。ヤコビ行列Ｊの最下行の合計インク量ＴＩに対する重みｗTIは、１としておく。
【００９０】
（２４）式の右辺第１項ｗＬ*（ΔＬ*−ΔＬ*t）2は、色彩値Ｌ*の変動量ΔＬ*，ΔＬ*tに関する２乗誤差である。これらの変動量ΔＬ*，ΔＬ*tは、次の式で与えられる。
【数２５】
2011139213.tif 000027
【数２６】
2011139213.tif 000028
【００９１】
前記（２５）式の右辺における偏微分値はヤコビ行列（（１４）式）で与えられる値であり、Ｉjは最適化処理の結果として得られるインク量、Ｉjrは現在のインク量である。第１の変動量ΔＬ*は、最適化処理によるインク量の変動量ΔＩjを、ヤコビ行列の成分である偏微分値で線形変換した量である。一方、第２の変動量ΔＬ*tは、ステップＴ１２０の平滑化処理で得られたターゲット値Ｌ*tと、現在インク量Ｉrで与えられる色彩値Ｌ*（Ｉr）との差分である。なお、第２の変動量ΔＬ*tは、平滑化処理の前後におけるＬ*値の差分と考えることが可能である。
【００９２】
前記（２４）式の右辺の第２項以降の各項も、前記（２５）式および（２６）式と同様の式で与えられる。すなわち、目的関数Ｅは、最適化処理によるインク量の変動量ΔＩjをヤコビ行列の成分で線形変換して得られる第１の変動量ΔＬ*，Δａ*，Δｂ*，ΔＧＩ…と、パラメータＬ*，ａ*，ｂ*，ＧＩ…に関する平滑化処理の前後における第２の変動量ΔＬ*t，Δａ*t，Δｂ*t，ΔＧＩt…と、の２乗誤差の和として与えられている。
【００９３】
ところで、第１の変動量ΔＬ*，Δａ*，Δｂ*，ΔＧＩ…は、行列を用いて以下の（２７）式および（２８）式の形式に書き表すことが可能である。
【数２７】
2011139213.tif 000029
【数２８】
2011139213.tif 000030
【００９４】
前記（２４）式は、行列を用いて（２９）式のように表記できる。
【数２９】
2011139213.tif 000031
【００９５】
ここで、Ｔは行列の転置を表している。行列ＷMはそれぞれ対角要素に重みを配置した対角行列（（３０）式参照）であり、行列ΔＭは各パラメータに関する目標変動量ベクトル（（３１）式参照）である。
【数３０】
2011139213.tif 000032
【数３１】
2011139213.tif 000033
【００９６】
（３１）式の右辺は、各パラメータＬ*，ａ*，ｂ*，ＣＩＩ…（「要素」とも呼ぶ）に関するターゲット値と、現在のインク量Ｉrで与えられる各パラメータ値との差分である。各パラメータのターゲット値のうち、色彩値Ｌ*t，ａ*t，ｂ*tは平滑化処理（ステップＴ１２０）で決定される。前記（３０）式における行列ＷMのｗIT＝１と設定する。
画質評価指数のターゲット値と現在の画質評価指数から求められる目標変動量ΔＧＩt，ΔＣＩＩt，ΔＴＩtについては、いくつかの決定方法がある。第１の方法は、目標変動量ΔＧＩt，ΔＣＩＩt，ΔＴＩtとして所定の定数（例えばΔＧＩt＝−２，ΔＣＩＩt＝−１，ΔＴＩt＝１）を使用する方法である。なお、定数としてマイナスの値を使用する理由は、これらの画質評価指数は、より小さいほど高画質であることを示す指数だからである。また、粒状性指数ＧＩのターゲット値ＧＩtは、ゼロとすることも好ましい。第２の方法は、ターゲット値ＧＩt，ＣＩＩt，ＴＩtを色彩値のターゲット値Ｌ*t，ａ*t，ｂ*tの関数として定義しておく方法である。以上のように、各パラメータのターゲット値は最適化処理前に決められているので、目標変動量ベクトルΔＭの各成分はすべて定数である。
【００９７】
前記（２９）式の右辺の各項のうち、第３項（ＩrTＪT+ΔＭT）ＷM（ＪＩr+ΔＭ）、は、最適化処理の結果として得られるインク量Ｉを含まないので定数である。一般に、最適化の目的関数Ｅとしては定数項は不要である。そこで、前記（２９）式から定数項を削除して全体に１／２を乗じると、次の（３２）式が得られる。
【００９８】
【数３２】
2011139213.tif 000034
【００９９】
ここで、以下の（３３）式および（３４）式のように行列Ａおよびベクトルｇを定義すると、前記（３２）式は（３５）式のように書き表せる。
【数３３】
2011139213.tif 000035
【数３４】
2011139213.tif 000036
【数３５】
2011139213.tif 000037
【０１００】
（３５）式で与えられる目的関数Ｅは、最適化で得られるインク量ベクトルＩに関する２次形式であることが理解できる。図１９（Ｃ）に示した（ＥＱ１）式と（ＥＱ２）式は、（２４）式と（３５）式とそれぞれ同じものである。
【０１０１】
本実施例の最適化処理では、（３５）式のような２次形式の目的関数Ｅを用いるので、最適化手法として２次計画法を使用することが可能である。ここで、「２次計画法」とは、逐次２次計画法を含まない狭義の２次計画法を意味している。２次形式の目的関数を用いた２次計画法を利用すれば、準ニュートン法や逐次２次計画法などの他の非線形計画法を利用する場合に比べて、処理を大幅に高速化することが可能である。
【０１０２】
次に、目的関数Ｅにおいて、調整係数Ｕjを含む項について検討する。前記（２９）式の右辺から、調整係数Ｕjを含んだ項のみを抽出すると、下記式のようになる。
【数３６】
2011139213.tif 000038
ここで、Ｚはインク量調整項を表す。ｗTI＝１，ΔＴＩt＝１であるため、前記（３６）式は下記（３７）式で表すことができる。
【数３７】
2011139213.tif 000039
前記（３７）式が示すようにインク量調整項Ｚは、現在のインク量Ｉjrからインク量Ｉjへの変動量（Ｉj−Ｉjr）を、調整係数Ｕjを重みとして線形結合した値の２乗に応じた値である。前記（３７）式をさらに展開した場合、調整係数Ｕjの２乗が各インクのインク量変動量（Ｉj−Ｉjr）の２次の係数となることが理解できる。
【０１０３】
ここで、あるインクの調整係数Ｕjが正値である場合、最適化における現在のインク量Ｉjrからインク量Ｉjへの増加量が大きければ大きいほど、該インクのインク量Ｉjに起因してインク量調整項Ｚは小さい値になる。インク量調整項Ｚが小さければ目的関数Ｅの値も小さくなるため、最適化において当該インクのインク量Ｉjの増加は促進されることとなる。さらに、正値のインクの調整係数Ｕjが大きければ大きいほど、インク量調整項Ｚおよび目的関数Ｅに対する寄与が大きくなるため、当該インクのインク量Ｉjの増加に対する促進度合いは大きくなると言える。反対に、あるインクの調整係数Ｕjが負値である場合、最適化における現在のインク量Ｉjrからインク量Ｉjへの増加量が大きければ大きいほど、該インクのインク量Ｉjに起因してインク量調整項Ｚは大きい値になる。インク量調整項Ｚが大きければ目的関数Ｅの値も大きくなるため、最適化において当該インクのインク量Ｉjの増加は抑制されることとなる。さらに、負値のインクの調整係数Ｕjの絶対値が大きければ大きいほど、インク量調整項Ｚおよび目的関数Ｅに対する寄与が大きくなるため、当該インクのインク量Ｉjの増加に対する抑制度合いは大きくなると言える。このとき、前記（３７）式は、インク量変動量（Ｉj−Ｉjr）と調整係数Ｕjの線形結合が１となることで、最もＺを小さくするが、Ｕjが正値の場合はインク量Ｉjが平均的に大きくなることで前記線形結合が１に近づき、Ｕjが負値の場合はインク量Ｉjが平均的に小さくなることで前記線形結合が１に近づく。
【０１０４】
以上説明したように、調整係数Ｕjが大きいインクほど、最適化におけるインク量Ｉjの増加が促進される。従って、印刷に多く使用させたいインクについては調整係数Ｕjを大きい値に設定すればよい。本実施例では、比較的濃度の大きいＣインクについて、低明度の領域において調整係数Ｕjを大きくし、高明度の領域において調整係数Ｕjを小さく設定している。そのため、高明度の領域におけるＣインクのインク量Ｉjを抑制することができる。逆の傾向となるように調整係数Ｕjを設定したＬｃインクについては、高明度の領域においてインク量Ｉjを多くさせることができる。また、図２０に示すように調整係数Ｕjをなだらかに変動させるようにしたため、急激にインク量Ｉjの傾向が変動することが防止でき、なめらかな階調性を実現することができる。以上説明したように、本実施例では、各インクについて個別に調整係数Ｕjを設定するようにしたため、各インク間でインク量Ｉjの相対的なバランスを調整することができる。
【０１０５】
前記（３０）式で示した行列ＷMは、各画質評価指数ＧＩ，ＣＩＩ，ＴＩの全体に対する重みを設定することができるが、重みｗTIによって個々のインクについて重みを設定することはできない。各インクのインク量Ｉjに対してターゲットインク量を設定し、これらの２乗誤差を各インクに応じた重みによって重み付けた項を目的関数Ｅに追加することも考えられるが、目的関数Ｅの複雑化や処理負荷の増加を招いてしまう。これに対して、本実施例のように前記（１４）式の本来のヤコビ行列の最下行の要素（インク量合計ＴＩの偏微分値）を前記（２２）式のように調整係数Ｕjに置き換えることにより、容易に各インクのインク量Ｉjの相対的なバランスを調整することができる。
【０１０６】
ところで、本実施例における最適化処理によるインク量の探索は、以下の条件の下で実行される。
（最適化条件）目的関数Ｅを最小とする。
（制約条件）インクデューティ制限を守る。
【０１０７】
インクデューティ制限値としては、例えば、個々のインクのインク量の最大許容値と、合計インク量の最大許容値とが使用される。例えば、各インクのインク量を８ビットで表現したとき、８種類の個々のインクのインク量Ｉjの最大許容値を１８０に設定し、合計インク量ΣＩjの最大許容値を２４０に設定してもよい。
【０１０８】
インクデューティ制限は、次の（３８）式で表すことができる。
【数３８】
2011139213.tif 000040
ここで、ベクトルｂは、デューティ制限の対象となるインク種類を識別するための係数であり、要素に０か１を持つベクトルである。例えば、１種類のインクに関するデューティ制限の場合には、ベクトルｂの１個の要素のみが１となる。一方、全インクの合計インク量に関するデューティ制限の場合には、ベクトルｂのすべての要素が１となる。（３８）式の右辺のlimIは、デューティ制限値である。
【０１０９】
各インク量Ｉjには、負でないという制約も存在する。この非負制限は、以下の（３９）式で表せる。
【数３９】
2011139213.tif 000041
【０１１０】
前記（３８）式と（３９）式とを合体すると、インクデューティ制限は、次の（４０）式で与えられる。
【数４０】
2011139213.tif 000042
【０１１１】
この（４０）式で表される制約は、線形不等号制約である。一般に、２次計画法は線形制約の下で実行することが可能である。すなわち、本実施例における最適化処理では、（４０）式の制約の下で、前記（３５）式で与えられる２次形式の目的関数Ｅを用いた２次計画法を実行することによって、最適なインク量を探索する。この結果、この線形制約を厳密に満足しつつ、インク量探索を高速に実行することが可能である。
【０１１２】
図２１は、最適化処理（図１７のステップＴ１３０）の詳細手順を示すフローチャートである。ステップＴ１３２では、まず、前記（３１）式で与えられる目標変動量ΔＭを求める。この目標変動量ΔＭは、前述したように、ステップＴ１２０（平滑化処理）で得られたターゲット値Ｌ*t，ａ*t，ｂ*tと現在インク量Ｉrに基づいて決定される。
【０１１３】
ステップＴ１３４では、前記（１４）式で与えられるヤコビ行列Ｊを算出する。なお、ヤコビ行列Ｊの各成分は、前記（２１）式で例示されるように、インク量の現在値Ｉr（平滑化／最適化前の値）に関して算出される値である。
【０１１４】
ステップＴ１３６では、ヤコビ行列Ｊによる線形変換の結果ΔＬ*，Δａ*，Δｂ*，ΔＧＩ…と、目標変動量ΔＭ（ΔＬ*t，Δａ*t，Δｂ*t，ΔＧＩt…）との差が最小になるように、インク量の最適化を実行する。この最適化は、前記（３５）式で与えられる２次形式の目的関数Ｅを用いた２次計画法を実行することによって実現される。
【０１１５】
なお、図１７のフローチャートにおいて既に説明したように、ステップＴ１３０の最適化処理の後、収束が不十分と判断される場合には、平滑化処理（ステップＴ１２０）および最適化処理（ステップＴ１３０）が再度実行される。この際、平滑化／最適化処理の初期値としては、その前の平滑化／最適化処理で得られた値が利用される。なお、このような繰り返し処理は必須ではなく、少なくとも１回の平滑化／最適化処理を行えばよい。
【０１１６】
このように、本実施例では、２次形式の目的関数Ｅを用いた２次計画法を実行することによって最適なインク量を探索するので、インク量探索を高速に実行することが可能である。発明者らの実測では、従来の準ニュートン法を用いた場合に比べて約１／１０の時間で処理を完了できることが見いだされた。
【０１１７】
（５）印刷装置の構成：
上述した実施形態では、コンピューター２０に複数のＬＵＴを備えさせていたが、プリンター３０のファームウェアに複数のＬＵＴを備えさせることにより、プリンター３０単体で本発明の印刷装置を構成することも可能である。図２２は、プリンター３０の構成を示している。同図において、プリンター３０はＣＰＵ３３とＲＡＭ３４とＲＯＭ３５とメモリーカードスロット４０とバス４１とＡＳＩＣ４２を備えている。ＲＯＭ３５に記憶されたプログラムデータ３５ａをＲＡＭ３４に展開しつつＣＰＵ３３がプログラムデータ３５ａにしたがった演算を行うことによりプリンター３０を制御するためのファームウェアＦＷが実行される。
【０１１８】
ファームウェアＦＷは、メモリーカードスロット４０に装着されたメモリーカードＭＣに記憶された印刷データＰＤに基づいて駆動データを生成する。ＡＳＩＣ４２は駆動データを取得し、紙送り機構３７やキャリッジモータ３８や印刷ヘッド３９の駆動信号を生成する。ＲＯＭ３５においては、色補正３Ｄ−ＬＵＴ６１０が記憶されている。プリンター３０はキャリッジ４１を備えており、キャリッジ３１はインクカートリッジ３２，３２・・・が取り付け可能なカートリッジホルダー３１ａ，３１ａ・・・を備えている。そして、各インクカートリッジ３２，３２・・・をカートリッジホルダー３１ａ，３１ａ・・・に挿入して固定することができる。
【０１１９】
本実施形態においてカートリッジホルダー３１ａ，３１ａ・・・は４個備えられているものとし、カートリッジホルダー３１ａ，３１ａ・・・にシアン（Ｃ）、マゼンタ（Ｍ）、イエロー（Ｙ）、ブラック（Ｋ）、淡シアン（Ｌｃ）、淡マゼンタ（Ｌｍ）、淡ブラック（Ｌｋ）、極淡ブラック（ＬＬｋ）のインクをそれぞれ充填した８種類のインクのインクカートリッジ３２，３２・・・が取り付けられる。なお、ＣＭＹインクは本発明の有彩色インクに相当し、Ｋインクは本発明の無彩色インクに相当する。キャリッジ３１は、インクカートリッジ３２，３２・・・から供給されるインクを多数のインクノズルＮｚから吐出する印刷ヘッド３９を備えている。
【０１２０】
図２３は、本実施形態のプリンター３０の印刷方式を模式的に示している。同図において、プリンター３０は、ＣＭＹＫのインクごとに複数のノズルＮｚ・・・を備えた印刷ヘッド３９を備えており、ノズルＮｚ・・・が吐出するインク量をＣＭＹＫごとに制御する処理が駆動データに基づいて行われる。各ノズルＮｚ・・・が吐出したインク滴は印刷用紙上において微細なドットとなり、多数のドットの集まりによって各インク量Ｉjに応じたインク被覆率の印刷画像が印刷用紙上に形成されることとなる。なお、各インクのインク量をＩjと表すものとし、下付文字ｊは各インクの種類を識別するための自然数（ｊ＝１〜８）とする。本明細書において、各インク量Ｉjに基づいて印刷することは、各インク量Ｉjに応じたインク被覆率の画像が印刷用紙上に再現されることを意味する。本実施形態におけるプリンター３０はインクジェット方式のプリンターであるが、インクジェット方式の他にも種々のプリンターに対して本発明を適用可能である。なお、本実施形態において、各インク量Ｉjは０〜２５５の８ビットで与えられるものとする。
【０１２１】
図２４は、ファームウェアＦＷのソフトウェア構成を示している。ファームウェアＦＷは、画像データ取得部ＦＷ１とレンダリング部ＦＷ２と色変換部ＦＷ３とハーフトーン部ＦＷ４とラスタライズ部ＦＷ５とから構成されている。画像データ取得部ＦＷ１は、メモリーカードＭＣに記憶された印刷データＰＤを印刷対象として取得する。印刷データＰＤは、文書データやグラッフィックデータであってもよいし、写真画像データであってもよい。レンダリング部ＦＷ２は、印刷データＰＤに基づいて印刷に使用する入力画像データＩＤを生成する。入力画像データＩＤは、印刷解像度（例えば２８８０×２８８０ｄｐｉ）に対応した画素数（印刷解像度×印刷実サイズ）の画素で構成されており、各画素が８ビット（０〜２５５）のＣＩＥ−ｓＲＧＢ色空間に準拠したＲＧＢ値で表現されている。
【０１２２】
色変換部ＦＷ３は、入力画像データＩＤを取得し、該入力画像データＩＤを色変換する。具体的には、色変換部ＦＷ３は、色補正３Ｄ−ＬＵＴ６１０を参照しつつ補間演算を実行することにより、ＲＧＢ値を各インクのインク量Ｉjに変換する。ハーフトーン部ＦＷ４は、色変換部ＦＷ３が出力した各インクのインク量Ｉjに基づくハーフトーン処理を実行する。ラスタライズ部ＦＷ５は、ハーフトーン処理後のハーフトーンデータの各画素（吐出可否）を印刷ヘッド３９の各主走査および各インクノズルＮｚに割り当て、駆動データを生成する。駆動データはＡＳＩＣ４２に出力され、ＡＳＩＣ４２が紙送り機構３７やキャリッジモータ３８や印刷ヘッド３９の駆動信号を生成する。
【０１２３】
（６）その他変形例：
なお、この発明は前記の実施例や実施形態に限られるものではなく、その要旨を逸脱しない範囲において種々の態様において実施することが可能であり、例えば次のような変形も可能である。
【０１２４】
前記実施例では、ＣＩＥ−Ｌａｂ色空間を概ね上下に分断する基準面Ｂ１，Ｂ２が設けることとしたが、色相に応じて調整係数Ｕjを制御するようにしてもよい。例えば、ある色相については、あるインクのインク量Ｉjの使用を他のインクよりも促進することが可能である。
【０１２５】
前記実施例では、機器独立表色系としてＣＩＥ−Ｌａｂ表色系を利用していたが、ＣＩＥ−ＸＹＺ表色系やＣＩＥ−Ｌ*ｕ*ｖ*表色系などの他の任意の機器独立表色系を利用することが可能である。但し、滑らかな色再現を実現するという意味からは、ＣＩＥ−Ｌａｂ表色系やＣＩＥ−Ｌ*ｕ*ｖ*表色系などの均等色空間である機器独立表色系を用いることが好ましい。
【０１２６】
前記実施例では、平滑化処理として力学モデルを利用した処理を採用していたが、他の種類の平滑化処理を採用してもよい。例えば、隣接する色点同士の間隔を測定し、その平均値になるべく近づくように個々の間隔を調整する平滑化処理を採用することも可能である。
【０１２７】
本明細書において「インク」とは、インクジェットプリンタやオフセット印刷等に用いられる液体状インクに限らず、レーザプリンタに用いられるトナーも含む広い意味で使用されている。このような「インク」の広い意味を有する他の用語としては、「色材」や「着色材」、「着色剤」、「記録材」を用いることも可能である。
【０１２８】
（７）まとめ： 
以上説明した実施形態によれば、ＬＵＴを参照しつつ入力画像データＩＤを変換して生成した記録材量データに基づいてプリンター３０に印刷を実行させるコンピューター２０において、入力画像データの表色系とプリンター３０のインク量データの表色系とを直接に対応付けるとともにそれぞれ特定の観察光源に対応付けて作成された複数の色変換テーブルを記憶した記憶媒体２４を備えており、複数の色変換テーブルは、色変換テーブルを参照しつつ印刷された印刷結果を該色変換テーブルに対応付けられた特定の観察光源下で観察したときに同一の入力画像データに基づく印刷結果同士で色の見えが一致するように作成されている。よって、複数の観察光源において色の見えが一致する印刷結果を提供することができるようになる。
【０１２９】
なお、本発明は上述した実施形態や変形例に限られず、上述した実施形態および変形例の中で開示した各構成を相互に置換したり組合せを変更したりした構成、公知技術並びに上述した実施形態および変形例の中で開示した各構成を相互に置換したり組合せを変更したりした構成、等も含まれる。
【符号の説明】
【０１３０】
１０…コンピューター、２０…コンピューター、２１…印刷制御部、２２…操作入力部、２３…表示部、２４…記憶媒体、２５…Ｉ／Ｆ、３０…プリンター、Ｍ１…データ取得部、Ｍ２…光源選択部、Ｍ３…色変換部、Ｍ４…ハーフトーン処理部、Ｍ５…コマンド作成部、１００…ベースＬＵＴ作成モジュール、１２０…スムージング処理初期値設定モジュール、１３０…スムージング処理モジュール、１３２…色点移動モジュール、１３４…インク量最適化モジュール、１３６…画質評価指数算出モジュール、１４０…テーブル作成モジュール、１４５…ＬＵＴ作成条件受付部、２００…色補正ＬＵＴ作成モジュール、３００…フォワードモデルコンバーター、３１０…分光プリンティングモデルコンバーター、３２０…色算出部、４００…ＬＵＴ格納部、４１０…インバースモデル初期ＬＵＴ、５１０…ベース３Ｄ−ＬＵＴ、５２０…ベース４Ｄ−ＬＵＴ、６１０…色補正３Ｄ−ＬＵＴ、６２０…色補正４Ｄ−ＬＵＴ 
【図１】
2011139213.tif 000043
【図２】
2011139213.tif 000044
【図３】
2011139213.tif 000045
【図４】
2011139213.tif 000046
【図５】
2011139213.tif 000047
【図６】
2011139213.tif 000048
【図７】
2011139213.tif 000049
【図８】
2011139213.tif 000050
【図９】
2011139213.tif 000051
【図１０】
2011139213.tif 000052
【図１１】
2011139213.tif 000053
【図１２】
2011139213.tif 000054
【図１３】
2011139213.tif 000055
【図１４】
2011139213.tif 000056
【図１５】
2011139213.tif 000057
【図１６】
2011139213.tif 000058
【図１７】
2011139213.tif 000059
【図１８】
2011139213.tif 000060
【図１９】
2011139213.tif 000061
【図２０】
2011139213.tif 000062
【図２１】
2011139213.tif 000063
【図２２】
2011139213.tif 000064
【図２３】
2011139213.tif 000065
【図２４】
2011139213.tif 000066