JP2015115738A - ルックアップテーブル生成方法、ルックアップテーブル生成装置、及び、色変換装置 - Google Patents

Abstract

【課題】印刷機の色材に合わせた色材使用量と画像形成装置の色材に合わせた色材使用量との対応関係を規定した一つのルックアップテーブルで特定画像の色に対する高精度の色再現を実現する。【解決手段】印刷機（３００）で印刷物３５０の形成に使用される第一の色材ＣＬ１の使用量と、画像形成装置（１００）で色再現画像１６０の形成に使用される第二の色材ＣＬ２の使用量と、の対応関係を規定したルックアップテーブル２００を生成する際、第一の色材ＣＬ１の色で表される入力色空間のうち特定画像４００に含まれる色ＵＣｉを含む第一の領域Ｒ１について第一の出力特性を示すように対応関係を規定し、入力色空間のうち第一の領域Ｒ１とは異なる第二の領域Ｒ２について第一の出力特性とは異なる第二の出力特性を示すように対応関係を規定して、ルックアップテーブル２００を生成する。【選択図】図３

Description

本発明は、ルックアップテーブル生成方法、ルックアップテーブル生成装置、及び、色変換装置に関する。

印刷本機を使用する前に印刷本機で形成される印刷物の色味等を確認するため、印刷プルーフ機でプルーフを形成して確認することが行われている。印刷本機にオフセット印刷機やグラビア印刷機を用いる場合、ＣＭＹＫ（シアン、マゼンタ、イエロー、ブラック）のみといった４種類の色材で印刷物が形成されている。ＤＤＣＰ（Direct Digital Color Proof）のため、印刷本機と同じＣＭＹＫの色材を印刷プルーフ機に用いることも行われている。一方、インクジェットプリンターには、ＣＭＹＫのインクの他、例えば、ｌｃ（ライトシアン）、ｌｍ（ライトマゼンタ）、ｌｋ（ライトブラック）、ｌｌｋ（ライトライトブラック）、Ｒ（レッド）、Ｏｒ（オレンジ）、Ｇｒ（グリーン）、等のインクの一部又は全部をも用いるものもある。このようなインクジェットプリンターは、高画質の写真画像を印刷可能である。

なお、特許文献１には、色再現性といった画質要素の最適化を目標に、ＬＵＴ（ルックアップテーブル）の格子点のインク量を最適化する技術が開示されている。
特許文献２には、ＣＭＹｌｃｌｍ（シアン、マゼンタ、イエロー、ライトシアン、ライトマゼンタ）の有彩色インク及び３段階の濃度のＫ１，Ｋ２，Ｋ３インクのインク量を規定したＬＵＴと、ＣＭＹｌｃｌｍの有彩色インク及び２段階の濃度のＫ１，Ｋ２インクのインク量を規定したＬＵＴとを使い分けることが開示されている。

特表２００７−５１６６６３号公報 特開２００７−１３７０２１号公報

印刷プルーフ機の色材に印刷本機と同じＣＭＹＫのみの色材を用いる場合、上述したインクジェットプリンターで得られるような高画質の写真印刷を行うことができない。

特許文献１に記載の技術は、色彩値への対応が明確なｓＲＧＢ色空間やＬ_*ａ_*ｂ_*色空間といった３次元色空間で定義された入力画像を想定しており、印刷本機で印刷物の形成に使用される色材のＣＭＹＫの４次元色空間といったデバイス依存の色空間の入力画像を対象としていない。従って、この技術は、ＣＭＹＫ画像といった入力画像の特徴を把握したうえで画質を確立する仕組みになっておらず、また、特定の画像に限定して色再現の精度を最適化する技術思想が含まれていない。
特許文献２に記載の技術は、ＬＵＴを複数保有する必要があり、ＬＵＴ作成に時間がかかるデメリットがある。この技術は、上記同様、色彩値への対応が明確なｓＲＧＢ色空間といった３次元色空間で定義された入力画像を想定しており、印刷本機で印刷物の形成に使用される色材のＣＭＹＫの４次元色空間といったデバイス依存の色空間の入力画像を対象としていない。従って、ＣＭＹＫ画像といった入力画像の特徴を把握したうえで画質を確立する仕組みになっておらず、また、特定の画像に限定して色再現の精度を最適化する技術思想が含まれていない。

印刷プルーフ機の色材を印刷本機の少ない色材に合わせると、高精度の色再現性を得ることができるものの、多色インクを使い分けることで得られる画質の多様性を実現することができない。従って、写真画質とプルーフ画質の両方を一つの印刷プルーフ機で得ることができない。

なお、上記のような問題は、オフセット印刷機及びグラビア印刷機のためのＬＵＴを生成する技術に限らず、種々の技術についても同様に存在する。

以上を鑑み、本発明の目的の一つは、印刷機の色材に合わせた色材使用量と画像形成装置の色材に合わせた色材使用量との対応関係を規定した一つのルックアップテーブルで特定画像の色に対する高精度の色再現を実現することにある。

上記目的の一つを達成するため、本発明は、印刷機で印刷物の形成に使用される第一の色材の使用量と、画像形成装置で色再現画像の形成に使用される第二の色材の使用量と、の対応関係を規定したルックアップテーブルを生成するルックアップテーブル生成方法であって、
前記第一の色材の色で表される入力色空間のうち特定画像に含まれる色を含む第一の領域について第一の出力特性を示すように前記対応関係を規定し、前記入力色空間のうち前記第一の領域とは異なる第二の領域について前記第一の出力特性とは異なる第二の出力特性を示すように前記対応関係を規定して、前記ルックアップテーブルを生成する、態様を有する。

すなわち、生成されるルックアップテーブルは、印刷機で使用される第一の色材の使用量と画像形成装置で使用される第二の色材の使用量との対応関係が規定される。この一つのルックアップテーブルに、入力色空間のうち特定画像に含まれる色を含む第一の領域については第一の出力特性を示すように対応関係が規定され、入力色空間のうち第一の領域とは異なる第二の領域については第一の出力特性とは異なる第二の出力特性を示すように対応関係が規定される。従って、上記態様は、印刷機の色材に合わせた色材使用量と画像形成装置の色材に合わせた色材使用量との対応関係を規定した一つのルックアップテーブルで特定画像の色に対する高精度の色再現を実現することが可能となる。

また、本発明の色変換装置は、印刷機で印刷物の形成に使用される第一の色材の使用量と、画像形成装置で色再現画像の形成に使用される第二の色材の使用量と、の対応関係を規定したルックアップテーブルであって、前記第一の色材の色で表される入力色空間のうち特定画像に含まれる色を含む第一の領域について第一の出力特性を示すように前記対応関係が規定され、前記入力色空間のうち前記第一の領域とは異なる第二の領域について前記第一の出力特性とは異なる第二の出力特性を示すように前記対応関係が規定されたルックアップテーブルを記憶する記憶部と、
前記ルックアップテーブルを参照して色変換を行う色変換部と、を備えた態様を有する。

すなわち、記憶部に記憶されるルックアップテーブルは、印刷機で使用される第一の色材の使用量と画像形成装置で使用される第二の色材の使用量との対応関係が規定されている。この一つのルックアップテーブルに、入力色空間のうち特定画像に含まれる色を含む第一の領域については第一の出力特性を示すように対応関係が規定され、入力色空間のうち第一の領域とは異なる第二の領域については第一の出力特性とは異なる第二の出力特性を示すように対応関係が規定される。従って、上記態様は、印刷機の色材に合わせた色材使用量と画像形成装置の色材に合わせた色材使用量との対応関係を規定した一つのルックアップテーブルで特定画像の色に対する高精度の色再現を実現することが可能となる。

ここで、上記第一の色材は、印刷機で印刷物の形成に使用される色材であればよく、ＣＭＹＫといった４色を組み合わせた色材の他、５色以上を組み合わせた色材、及び、３色以下を組み合わせた色材が含まれる。上記第二の色材は、画像形成装置で色再現画像の形成に使用される色材であればよく、５色以上を組み合わせた色材の他、４色以下を組み合わせた色材が含まれる。
上記特定画像に含まれる色を含む第一の領域の例としては、肌色領域、グレイ領域、カラーチャート領域、等を挙げることができる。第一の領域は、離散的な部分の集まりでもよい。

出力特性は、領域内の全体的な出力の傾向を意味する。例えば、前記入力色空間に含まれる第一及び第二の領域について示される前記第一及び第二の出力特性は、それぞれ、所定の観察光源における前記第一の色材の色彩値と前記第二の色材の色彩値との色差、前記色再現画像の粒状性を表す指標値、前記色再現画像の階調性を表す指標値、の少なくとも一つで表されてもよい。この態様は、特定画像の色に対する高精度の色再現を実現好適な例を提供することができる。なお、観察光源は、複数の光源で構成されてもよい。

ところで、前記特定画像の色分布に基づいて前記入力色空間内の複数の格子点のうち前記特定画像に含まれる色の変換に用いる格子点を特定し、前記第一の領域について前記第一の出力特性を示すように前記特定した格子点について前記第一の出力特性を示すように前記対応関係を規定してもよい。この態様は、ルックアップテーブルの生成を容易にすることができる。

また、前記特定画像から該特定画像に含まれる色を抽出し、前記入力色空間内の複数の格子点のうち前記抽出した色の変換に用いる格子点を特定し、前記第一の領域について前記第一の出力特性を示すように前記特定した格子点について前記対応関係を規定してもよい。この態様は、ルックアップテーブルの生成をさらに容易にすることができる。

さらに、前記特定画像に含まれる色のそれぞれについて前記入力色空間内の複数の格子点のうち前記色の変換に用いる一以上の格子点を抽出し、該抽出した格子点の重複を解消し、前記第一の領域について前記第一の出力特性を示すように重複解消後の格子点について前記対応関係を規定してもよい。この態様は、ルックアップテーブルの生成をさらに容易にすることができる。

さらに、前記ルックアップテーブルに設けられる第一格子点の数Ｎ１よりも少ない数Ｎ２の第二格子点について前記対応関係を規定し、該第二格子点についての対応関係に基づいて前記数Ｎ１の第一格子点について前記対応関係を規定してもよい。ルックアップテーブルに設けられる第一格子点の数Ｎ１が多い場合、いきなり全ての数Ｎ１の第一格子点について対応関係を規定するのは時間がかかり過ぎることがある。本態様は、第一格子点の数Ｎ１よりも少ない数Ｎ２の第二格子点について対応関係を規定すればよいので、ルックアップテーブルの生成を容易にすることができる。

さらに、前記印刷物を観察する複数の観察光源のそれぞれについて、前記第二格子点に、前記印刷物に形成される前記第一の色材の前記観察光源におけるターゲット色彩値を対応付けてもよい。そのうえで、前記複数の観察光源について前記第二格子点において対応付けられたターゲット色彩値への近似性を評価する評価値に基づいて、前記色再現画像に形成される前記第二の色材の色彩値が前記観察光源毎に前記ターゲット色彩値に近似するように前記第二の色材の使用量を前記第二格子点について予測してもよい。そのうえで、前記第二格子点に、予測された前記第二の色材の使用量を対応付け、前記第二格子点についての前記第一及び第二の色材の使用量の対応関係に基づいて前記数Ｎ１の第一格子点について前記ルックアップテーブルの対応関係を規定してもよい。本態様は、複数の観察光源下で印刷機による印刷物に見える色味変化をより忠実に再現可能なルックアップテーブルを生成することができる。

さらに、前記入力色空間において前記第一格子点を囲む複数の前記第二格子点についての前記第一及び第二の色材の使用量の対応関係を用いた補間演算を行うことにより前記第一格子点について前記ルックアップテーブルの対応関係を規定してもよい。この態様は、ルックアップテーブルの生成をさらに容易にすることができる。

さらに、前記第一の色材は、第三の色材と、第四の色材と、を含んでもよい。そのうえで、前記第一格子点に対応する前記第三の色材の使用量のそれぞれについて、前記第四の色材で表される色空間内に前記第二格子点を設けてもよい。そのうえで、前記第三の色材の使用量のそれぞれについて、前記第四の色材で表される色空間において前記第一格子点を囲む複数の前記第二格子点についての前記第一及び第二の色材の使用量の対応関係を用いた補間演算を行うことにより前記第一格子点について前記ルックアップテーブルの対応関係を規定してもよい。本態様は、第一の色材の色で表される入力色空間よりも次元数の少ない色空間において補間演算が行われるので、ルックアップテーブルの生成をさらに容易にすることができる。なお、第三の色材の例としては、Ｋインクといった無彩色インク等を挙げることができる。第四の色材の例としては、ＣＭＹのインクといった有彩色インク等を挙げることができる。
なお、第三の色材の使用量がＮｇ段階あるとき、このＮｇ段階よりも少ない段階についてのみ第四の色材で表される色空間内に設けた第二格子点についての対応関係を用いた補間演算を行うことも、上記態様に含まれる。

さらに、前記第一格子点に対応する前記第三の色材の使用量のそれぞれについて、前記第二格子点に含まれ前記第三の色材の使用量に応じた数の第三格子点を前記入力色空間内に設けてもよい。そのうえで、前記第三格子点のそれぞれについて前記対応関係を規定してもよい。本態様は、入力色空間内に設けられる第三格子点の総数が第三の色材の使用量に応じた数とされるので、ルックアップテーブルの生成をさらに容易にすることができる。むろん、第三の色材の使用量のＮｇ段階よりも少ない段階についてのみ第四の色材で表される色空間内に設けた第二格子点についての対応関係を用いた補間演算を行うことも、前述の態様に含まれる。

さらに、前記第一格子点に対応する前記第三の色材の使用量に第一の使用量と該第一の使用量よりも多い第二の使用量が含まれるとき、前記第三の色材の前記第一の使用量に対応する前記第三格子点の数は、前記第三の色材の前記第二の使用量に対応する前記第三格子点の数よりも多くてもよい。第三の色材の使用量が比較的少ない第一の使用量である場合、第四の色材の色で表される色空間のガマットは広くなる。一方、第三の色材の使用量が比較的多い第二の使用量である場合、第四の色材の色で表される色空間のガマットは狭くなる。そこで、第三の色材の第一の使用量に対応する第三格子点の数が第三の色材の第二の使用量に対応する第三格子点の数よりも多いと、入力色空間内に設けられる第三格子点の総数がより適正な数とされるので、ルックアップテーブルの生成をさらに容易にすることができる。

さらに、本発明は、ルックアップテーブル、ルックアップテーブルを記録したコンピューター読み取り可能な媒体、ルックアップテーブルを搭載した画像形成装置、上記ルックアップテーブル生成方法に対応した手段を備えるルックアップテーブル生成装置、画像形成装置を備えた画像形成システム、前述の装置又はシステムの各部に対応した機能をコンピューターに実現させるプログラム、該プログラムを記録したコンピューター読み取り可能な媒体、等に適用可能である。前述の装置及びシステムは、分散した複数の部分で構成されてもよい。

プルーフシステムＳＹ１の構成を模式的に例示する図。 ルックアップテーブル２００の構造を模式的に例示する図。 ルックアップテーブル２００のデータ構造を模式的に例示する図。 特定画像４００を模式的に例示する図。 コンピューター１０の構成を模式的に例示するブロック図。 （ａ）は少格子点ルックアップテーブル２１０のデータ構造を模式的に例示する図、（ｂ）は部分ルックアップテーブル２３０のデータ構造を模式的に例示する図。 少格子点ルックアップテーブル２１０の色材量予測格子点Ｇ５の配置例を模式的に示す図。 ルックアップテーブル生成処理を例示するフローチャート。 特有色参照格子点抽出処理を例示するフローチャート。 最終的なルックアップテーブル２００を生成する処理を例示するフローチャート。 設定画面５００を例示する図。 特有色参照格子点Ｇ２，Ｇ３を抽出する様子を模式的に例示する図。 第三の色材ＣＬ３の使用量毎の色材量予測格子点Ｇ５の配置例を模式的に示す図。 印刷色プロファイルＰＲ０の構造を模式的に例示する図。 印刷色プロファイルＰＲ０のデータ構造を模式的に例示する図。 印刷色プロファイルＰＲ０の生成方法例を模式的に説明するための図。 分光反射率に基づいて色彩値を算出するための計算例を説明する図。 インク量セットを最適化する処理の流れを模式的に例示する図。 インク量セットが最適化されていく様子を模式的に例示する図。 画像形成装置で行われる色再現画像出力制御処理を例示するフローチャート。 分光反射率データベースの構造を模式的に例示する図。 分光ノイゲバウアモデルを模式的に例示する図。 セル分割ユール・ニールセン分光ノイゲバウアモデルを模式的に例示する図。 元ルックアップテーブル２０５の格子点Ｇ６からルックアップテーブル２００の格子点Ｇ１を補間する様子を模式的に例示する図。 スムージング処理を模式的に例示する図。 色再現画像の粒状性を変える様子を模式的に例示する図。 粒状性指標値ＩＮ１を算出する様子を模式的に例示する図。 色再現画像の階調性を変える様子を模式的に例示する図。

以下、本発明の実施形態を説明する。むろん、以下の実施形態は本発明を例示するものに過ぎず、実施形態に示す特徴の全てが発明の解決手段に必須になるとは限らない。

（１）プルーフシステムの概要：
図１は、印刷本機（印刷機）３００で形成される印刷物３５０のプルーフ１５０をプルーフ機（画像形成装置）１００で形成するプルーフシステムＳＹ１を模式的に例示している。印刷本機で直接、印刷物を形成して色味等を確認すると、高いコストがかかる。そこで、印刷物３５０を形成しなくても印刷物３５０上の印刷画像３６０の色味等を確認可能とするため、プルーフ機１００は、印刷物３５０の色再現画像１６０を形成する。

印刷本機（印刷機）３００には、オフセット印刷機、グラビア印刷機、等が含まれる。図１に示す印刷本機３００は、ＣＭＹＫ（シアン、マゼンタ、イエロー、ブラック）の第一の色材ＣＬ１を使用して被印刷物Ｍ１に印刷画像３６０を形成する。なお、第一の色材ＣＬ１は、Ｋの色材（第三の色材ＣＬ３）と、この色材ＣＬ３の色（Ｋ）に合わせて置き換え可能なＣＭＹの色材（第四の色材ＣＬ４）と、を含む。

プルーフ機（画像形成装置）１００には、インクジェットプリンター、ワイヤドットプリンター、レーザープリンター、ラインプリンター、複写機、ファクシミリ、これらの一部を組み合わせた複合機、等が含まれる。図１に示すプルーフ機１００は、ＣＭＹＫｌｃｌｍの第二の色材ＣＬ２を使用して被印刷物Ｍ２に色再現画像１６０を形成する。ライトシアン（ｌｃ）は、シアンと同系統の色相でシアンよりも明度が高い色である。ライトマゼンタ（ｌｍ）は、マゼンタと同系統の色相でマゼンタよりも明度が高い色である。むろん、第二の色材ＣＬ２は、ＣＭＹＫＲＯｒＧｒの色材等でもよい。レッド（Ｒ）やオレンジ（Ｏｒ）やグリーン（Ｇｒ）は、ＣＭＹに置き換え可能な色である。図１に示すプルーフ機１００は、ＬＵＴ（ルックアップテーブル）２００に従って、ＣＭＹＫの第一の色材ＣＬ１の各使用量Ｄｃ，Ｄｍ，Ｄｙ，ＤｋをＣＭＹＫｌｃｌｍの第二の色材ＣＬ２の各使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmに変換し、得られる使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmに従って色再現画像１６０を形成する。インク使用量を被印刷物Ｍ２上のドットに変換する場合、各使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmを表す階調データに対してディザ法や誤差拡散法や濃度パターン法といった所定のハーフトーン処理を行って前記階調データの階調数を減らし、得られる多値データに従ってインク滴を吐出して被印刷物Ｍ２上にインクドットを形成する。多値データは、ドットの形成状況を表すデータであり、ドットの形成有無を表す２値データでもよいし、大中小の各ドットといった異なるサイズのドットに対応可能な３階調以上の多値データでもよい。２値データは、例えば、ドット形成に１、ドット無しに０、を対応させるデータとすることができる。４値データとしては、例えば、大ドット形成に３、中ドット形成に２、小ドット形成に１、ドット無しに０、を対応させるデータとすることができる。得られる色再現画像は、被印刷物Ｍ２上のドットの形成状況で表現される。

印刷画像３６０や色再現画像１６０は、観察光源Ｌ０の種類によって色味が変わってくる。ここで、符号Ｌ０は、個々の光源Ｌ１〜Ｌ３を総称するときに使用する。規格化された観察光源Ｌ０には、Ｄ５０光源、Ｄ５５光源、Ｄ６５光源、Ｄ７５光源、Ａ光源、Ｆ２光源、Ｆ７光源、Ｆ１０光源、Ｆ１１光源、等が含まれる。色味の変化は、色材の種類にも依存する。
例えば、印刷業界では、実際の照明には無い分光分布を示すＤ５０光源が標準光源である。Ｄ５０光源下で印刷性能を見るため、標準のＤ５０光源で見たときの色精度を保証することは印刷機のプルーフ機として重要なことである。一方、印刷本機による印刷物やプルーフ機による色再現画像を実際に見る環境は、Ｄ５０光源とは異なる光源であることが想定され、Ｄ６５光源とＡ光源の混合環境といった複数の光源が混ざった環境も想定される。従って、Ｄ５０光源Ｌ１下で印刷画像３６０と色再現画像１６０との色味が合っていても、別の光源Ｌ２，Ｌ３下では印刷画像３６０と色再現画像１６０との色味が合わないことがある。また、Ｄ５０光源が実際の照明には無い分光分布を示すため、両画像３６０，１６０の色差を小さくすること、すなわち、定量的な高い色精度の実現を追求しても、目視観察における色精度の向上、すなわち、主観的な高い色精度の実現に繋がらないことがある。本技術は、「定量的な高い色精度の実現」と「主観的な高い色精度の実現」を両立するのに有効である。

なお、ＤＤＣＰ機（Direct Digital Color Proofer；ダイレクト・デジタル・カラー・プルーファー）には、印刷本機と同じＣＭＹＫのみの色材を用いるものもある。印刷プルーフ機に印刷本機と同じＣＭＹＫの色材を用いると、高精度の色再現性を実現することができるものの、より多くの色のインクを用いるインクジェットプリンターの特徴である写真画質といった高画質の色再現画像等、ユーザーの意図を反映した色再現画像を得ることができない。また、このようなＤＤＣＰ機は、高価である。
また、ＤＤＣＰ機は、印刷本機で形成される印刷画像を忠実に再現することができるものの、ＣＭＹＫの色材で形成される印刷画像の特徴を把握したうえでユーザーが出力したい画像コンテンツに応じてメタメリックマッチングの目標を切り替えたり画質を確立したりする仕組みになっていない。

本技術は、上述したＤＤＣＰ機にある課題を解決するものであり、図２，４等に示す特定画像４００の色分布を解析してＬＵＴ２００を生成することによりユーザーが出力したい画像コンテンツに応じた高画質の色再現画像１６０を得るものである。例えば、写真画像に合わせて肌色領域についてはＦ１０，Ｆ２光源下のメタメリックマッチングが非常に良好で残りの領域についてはＤ５０光源のメタメリックマッチングが非常に良好であるといった、インクジェットプルーフ機等のプルーフ機としての良好なメタメリックマッチング機能を得ながらユーザーが出力したい写真画像の見た目を向上させる色再現画像１６０を一つのＬＵＴ２００で得ることができる。
なお、第一の色材ＣＬ１の各使用量Ｄｃ，Ｄｍ，Ｄｙ，Ｄｋと第二の色材ＣＬ２の各使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmとの対応関係を規定したＬＵＴ２００は、印刷本機依存のＣＭＹＫ４次元色空間のデータをプルーフ機依存のＣＭＹＫｌｃｌｍ色空間のデータに変換する点から色変換ＬＵＴと言える。また、ＬＵＴ２００は、ＣＭＹとＫの使用割合を変換し、Ｃとｌｃの使用割合を変換し、Ｍとｌｍの使用割合を変換する点から色分版ＬＵＴとも言える。

（２）ルックアップテーブルの構造の概要：
図２は、ＬＵＴ２００の構造を模式的に例示している。図３は、各格子点（入力点）Ｇ１に対応付けられる指標値とともにＬＵＴ２００のデータ構造を模式的に例示している。なお、ＬＵＴ２００に示した指標値は、模式的な数値であり、出力特性を限定するものではない。また、指標値は、ＬＵＴ２００のデータに含まれなくてもよい。
本技術は、特定画像４００の色（ユニークカラーＵＣｉ）を高精度で再現する比較的高密度の特有色参照格子点Ｇ２，Ｇ３を共存させることにより、特定画像４００に特化したＬＵＴ２００の生成を可能にし、その結果、より印刷本機の色の見えに近いインクジェットプルーフ等のプルーフを実現するものである。

ＬＵＴ２００の入力はＣＭＹＫの４次元であるため、図示できない。そこで、ある１次元を１点に固定した３次元状にＬＵＴ２００を模式的に示すこととし、図２では、Ｋ値をある値に固定してＣＭＹの３次元空間にＬＵＴ２００の格子点（第一格子点）Ｇ１を示している。例えば、ＣＭＹＫそれぞれの使用量についてＮｇ段階（Ｎｇは２以上の整数）の格子点を設ける場合、格子点数はＮｇ⁴個となる。Ｎｇ＝１７である場合、格子点数はＮ１＝１７⁴＝８３５２１個となる。Ｋの使用量がＮｇ段階あるので、図２に示すようなＣＭＹ３次元色空間がＮｇ個存在することになる。
なお、格子点は、ＬＵＴに規定される入力点を意味する一般的用語であり、入力色空間の座標に対応していれば配置に特別な限定は無い。従って、複数の格子点は、入力色空間内で均等に配置されるのみならず、入力色空間内でＣＭＹＫのそれぞれについて異なる位置にあれば不均等に配置されてもよい。

図２に示すように、ＣＭＹＫ色空間内の格子点Ｇ１の座標（位置）は、第一の色材ＣＬ１の各使用量Ｄｃ，Ｄｍ，Ｄｙ，Ｄｋを示している。格子点インク量（第二の色材ＣＬ２の各使用量）ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmは、各格子点Ｇ１に格納される。このようにして、ＬＵＴ２００には、各格子点Ｇ１について、ＣＭＹＫの色材ＣＬ１の使用量とＣＭＹＫｌｃｌｍの色材ＣＬ２の使用量との対応関係が規定されている。

本ＬＵＴ２００は、第一の色材ＣＬ１の色で表されるＣＭＹＫ色空間（入力色空間）のうち特定画像４００に含まれるユニークカラーＵＣｉを含む第一の領域Ｒ１について第一の出力特性を示すように対応関係が規定され、ＣＭＹＫ色空間（入力色空間）のうち第二の領域Ｒ２について第一の出力特性とは異なる第二の出力特性を示すように対応関係が規定されている。ＬＵＴ２００における第一及び第二の領域Ｒ１，Ｒ２は、互いに異なる領域（重複しない領域）である。第一の領域Ｒ１は、ユーザーが出力した画像コンテンツに含まれる色の範囲を意味し、主観的な高い色精度を実現する領域である。第一の領域Ｒ１はＬＵＴ２００に複数設けられてもよく、各領域Ｒ１は異なる出力特性を示してもよい。図２等に示す第二の領域Ｒ２は、ＣＭＹＫ色空間を包括する包括領域Ｒ９のうち第一の領域Ｒ１を除いた領域であり、印刷本機３００で形成される印刷画像３６０の色を忠実に再現する領域、すなわち、定量的な高い色精度を実現する領域である。特性は、全体的な傾向を意味する。第一及び第二の出力特性は、互いに異なる出力特性である。図３に例示する第一の領域Ｒ１の色は、Ｆ１０光源とＦ２光源とを均等に混合した観察光源において印刷本機による印刷画像３６０とプルーフ機による色再現画像１６０との色差ΔＥ_F10-F2が所定の範囲内（例えば１．０以下）となる使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmに変換される。また、図３に例示する第二の領域Ｒ２の色は、Ｄ５０光源において印刷画像３６０と色再現画像１６０との色差ΔＥ_D50が所定の範囲内（例えば１．０以下）となる使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmに変換される。図３の例では、出力特性が所定の観察光源における第一の色材の色彩値と第二の色材の色彩値との色差で表されている。色彩値は、例えば、ＣＩＥ（国際照明委員会）Ｌ^*ａ^*ｂ^*色空間で定義されるＬ^*，ａ^*，ｂ^*で表される。第一の色材の色彩値をＬ_D，ａ_D，ｂ_D、第二の色材の色彩値をＬ_d，ａ_d，ｂ_dで表す場合、色差ΔＥは、ΔＥ＝｛（Ｌ_d−Ｌ_D）²＋（ａ_d−ａ_D）²＋（ｂ_d−ｂ_D）²｝^1/2により算出してもよいし、ＣＩＥ ＤＥ２０００の色差式（ΔＥ₂₀₀₀）に基づいて算出してもよい。

以上のことから、ＣＭＹＫの第一の色材ＣＬ１の使用量をプルーフ機１００で使用される第二の色材ＣＬ２の使用量に変換する一つのＬＵＴ２００に、特定画像４００に含まれる色に対する最適な色再現特性と、特定画像４００に含まれない色に対する色再現特性との、異なる特性が持たされている。

なお、格子点Ｇ１は入力色空間内に離散的に配置されるため、図２に示すように、領域Ｒ１の境界部分Ｒｂは格子点Ｇ１の位置に一致するとは限らない。図２では、補間演算により領域Ｒ１内の色を変換するために用いられる格子点を符号Ｇ２，Ｇ３で示し、領域Ｒ１内の格子点を符号Ｇ２で示している。格子点Ｇ３は、領域Ｒ１外ではあるが出力値ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmの補間のために必要な隣接格子点である。符号Ｒ１ｇは、第一の領域Ｒ１の色の変換に用いる格子点Ｇ２，Ｇ３の範囲を示している。

図４は、特定画像４００を模式的に例示している。特定画像４００は、標準化された画像に基づいて予め用意された色分布解析用の画像であり、高い色精度を得たい画像としてユーザーが選択可能な画像である。特定画像４００には、ポートレート画像４０１、グレイ画像４０２、肌色画像４０３、等が含まれる。ポートレート画像４０１には、ＳＣＩＤ（Standard Color Image Data）といった標準化された画像が含まれる。ＳＣＩＤは、印刷等における作業の標準化や符号化特性の評価等を行うことを目的に国際標準化（ISO12640、JIS X 9201）された標準画像データセットのことであり、ＣＭＹＫの色材量で記述された写真データ等から構成される。グレイ画像４０２は、例えば、ＳＣＩＤからグレイ領域を抽出した画像とすることができる。肌色画像４０３は、例えば、ＳＣＩＤから肌色領域を抽出した画像とすることができる。

特定画像４００の色に基づいて設定される第一の領域Ｒ１は、入力色空間内の一部の領域であればよく、肌色領域、グレイ領域、カラーチャート領域、インクセーブ領域、その他特色領域、等を例示することができ、分割されていないひとまとまりの領域のみならず、離散的な領域も含む。肌色領域やグレイ領域は、人の目に敏感な色の領域の代表例である。カラーチャート領域は、カラーチャートを形成するための色の領域であり、離散的な領域でもよい。インクセーブ領域は、第二の色材ＣＬ２の使用量を節約する領域である。特色領域は、人の目に敏感な色の領域を総称している。

領域Ｒ１，Ｒ２について示される出力特性は、所定の観察光源における第一の色材ＣＬ１の色彩値と第二の色材ＣＬ２の色彩値との色差ΔＥ、色再現画像１６０の粒状性を表す粒状性指標値（ＩＮ１とする。）、色再現画像１６０の階調性を表す階調性指標値（ＩＮ２とする。）、の少なくとも一つ等で定量化される。これらの指標値ΔＥ，ＩＮ１，ＩＮ２等を画質定量指数ＩＮと総称することにする。

図１１に示す設定画面５００の例では、出力特性として、「Ｆ１０−Ｆ２」、「Ｆ８−Ｆ２」、「Ｃ−Ｄ５０」、「Ｄ５０」、「粒状性小」、「階調性」、等が示されている。「Ｆ１０−Ｆ２」は、例えば、見た目を重視することを意図し、Ｆ１０光源とＦ２光源とを１：１で混合した蛍光灯を意味する観察光源の下で印刷画像３６０の色彩値と色再現画像１６０の色彩値との色差ΔＥ_F10-F2が所定の範囲内（例えば１．０など正の閾値以下）となる出力特性とする。むろん、見た目を重視することを意図する観察光源は、Ｄ５０光源とＦ１０光源とＦ２光源の１：１：１の混合光源等の観察光源の下で両画像３６０，１６０の色彩値との色差が所定の範囲内となる出力特性としてもよい。「Ｆ８−Ｆ２」は、例えば、自然画像の色の精度を重視することを意図し、Ｆ８，Ｆ２光源を１：１で混合した蛍光灯を意味する観察光源の下で両画像３６０，１６０の色彩値の色差ΔＥ_F8-F2が所定の範囲内（例えば１．０など正の閾値以下）となる出力特性とする。「Ｃ−Ｄ５０」は、例えば、自然画像以外の色の精度を重視することを意図し、Ｃ，Ｄ５０光源を１：１で混合した観察光源の下で両画像３６０，１６０の色彩値の色差ΔＥ_C-D50が所定の範囲内（例えば１．０など正の閾値以下）となる出力特性とする。「Ｄ５０」は、例えば、Ｄ５０光源の下で両画像３６０，１６０の色彩値の色差ΔＥ_D50が所定の範囲内（例えば１．０など正の閾値以下）となる出力特性とする。「粒状性小」は例えば、色再現画像１６０の粒状性指標値ＩＮ１が所定の範囲内（例えば正の閾値以下）となる出力特性とする。「階調性」は、例えば、色再現画像１６０の階調性指標値ＩＮ２が所定の範囲内（例えば正の閾値以下）となる出力特性とする。

（３）ルックアップテーブル生成装置、及び、色変換装置の概要：
図５は、ＬＵＴ（ルックアップテーブル）生成装置Ｕ０を含むコンピューター１０の構成を例示している。コンピューター１０では、ＣＰＵ（Central Processing Unit）１１、ＲＡＭ（Random Access Memory）１２、ＲＯＭ（Read Only Memory）１３、ハードディスクドライブ（ＨＤＤ）１４、汎用インターフェイス（ＧＩＦ）１５、ビデオインターフェイス（ＶＩＦ）１６、入力インターフェイス（ＩＩＦ）１７、等がバス１８に接続されて互いに情報を入出力可能とされている。ＨＤＤ１４には、オペレーティングシステム（ＯＳ）及びＬＵＴ生成プログラムを含む各種プログラムを実行するためのプログラムデータ１４ａ等が記憶されている。ＨＤＤ１４には、ＬＵＴ生成処理により生成される少格子点ＬＵＴ２１０（図６（ａ）参照）、部分ＬＵＴ２３０（図６（ｂ）参照）、最終ＬＵＴ（最終的なＬＵＴ）２００、等も記憶される。ＣＰＵ１１は、プログラムデータ１４ａを適宜ＲＡＭ１２へ読み出し、プログラムデータ１４ａに従ってコンピューター１０全体を制御する。ＬＵＴ生成プログラムは、ＬＵＴ生成装置Ｕ０の各部Ｕ１〜Ｕ５，Ｕ１０に対応した機能をコンピューター１０に実現させ、コンピューター１０をＬＵＴ生成装置Ｕ０として機能させる。

ＧＩＦ １５には、画像出力装置であるプリンター２０、測色機やスキャナーといった画像入力装置３０、等が接続される。ＧＩＦ １５には、ＵＳＢ（Universal Serial Bus）等を採用することができる。プリンター２０は、プルーフ機１００でもよいし、プルーフ機１００でなくてもよい。ＶＩＦ １６には、画像出力装置であるディスプレイ４０等が接続される。ＩＩＦ １７には、操作入力装置５０であるキーボード５０ａ、同じく操作入力装置５０であるポインティングデバイス５０ｂ、等が接続される。ポインティングデバイス５０ｂには、マウス等を用いることができる。

図５に示すＬＵＴ生成装置Ｕ０は、条件設定部Ｕ１、特有色参照格子点抽出部Ｕ２、部分ＬＵＴ生成部Ｕ３、少格子点ＬＵＴ生成部Ｕ４、最終ＬＵＴ生成部Ｕ５、及び、予測部Ｕ１０を備え、上記ＬＵＴ２００を生成する。条件設定部Ｕ１は、色分布を解析する特定画像４００の選択入力を受け付けたり、特定画像４００に含まれる色を含む領域Ｒ１に対応付ける出力特性の入力を受け付けたりする。特有色参照格子点抽出部Ｕ２は、特定画像４００の色分布（図１２に示すユニークカラーリスト４２０）に基づいてＣＭＹＫ色空間内の複数の格子点Ｇ１のうち特定画像４００に含まれるユニークカラーＵＣｉの変換に用いる特有色参照格子点Ｇ２，Ｇ３を特定する。部分ＬＵＴ生成部Ｕ３は、特有色参照格子点であり色材量予測格子点である格子点Ｇ２，Ｇ３について、予測部Ｕ１０により予測された第二の色材ＣＬ２の使用量を第一の色材ＣＬ１の使用量に対応付けて部分ＬＵＴ２３０を生成する。少格子点ＬＵＴ生成部Ｕ４は、最終ＬＵＴ２００よりも格子点の少ない少格子点ＬＵＴ２１０における色材量予測格子点Ｇ５（図６（ａ）参照）について、予測部Ｕ１０により予測された第二の色材ＣＬ２の使用量を第一の色材ＣＬ１の使用量に対応付けて少格子点ＬＵＴ２１０を生成する。少格子点ＬＵＴ２１０と部分ＬＵＴ２３０は、最終ＬＵＴ２００生成用の元ＬＵＴ２０５となる。最終ＬＵＴ生成部Ｕ５は、元ＬＵＴ２０５の色材量予測格子点（第二格子点）Ｇ２，Ｇ３，Ｇ５に対応付けられた第二の色材ＣＬ２の使用量に基づいて、第一及び第二の色材ＣＬ１，ＣＬ２の使用量の対応関係を規定した最終ＬＵＴ２００を生成する。ＬＵＴ２００は、ＨＤＤ１４に記憶される。予測部Ｕ１０は、観察光源Ｌ０について第一の色材ＣＬ１の使用量（例えばＤｃ，Ｄｍ，Ｄｙ，Ｄｋ）に対応付けられたターゲット色彩値（例えばＬ^*ａ^*ｂ^*値）への近似性を評価する評価値Ｉ（後述）に基づいて、色再現画像１６０に形成される第二の色材ＣＬ２の色彩値（例えばＬ^*ａ^*ｂ^*値）が観察光源Ｌ０毎にターゲット色彩値に近似するように第二の色材ＣＬ２の使用量を予測する。

図６（ａ）は、少格子点ＬＵＴ生成部Ｕ４で生成される少格子点ＬＵＴ２１０のデータ構造を模式的に例示している。図６（ｂ）は、部分ＬＵＴ生成部Ｕ３で生成される部分ＬＵＴ２３０のデータ構造を模式的に例示している。図７は、少格子点ＬＵＴ２１０の色材量予測格子点Ｇ５の配置例を模式的に示している。なお、色材量予測格子点Ｇ５の位置は、最終ＬＵＴ２００の格子点Ｇ１の位置と一致してもよいが、格子点Ｇ１の色材量を予測するために補間演算を行うため、最終ＬＵＴ２００の格子点Ｇ１の位置と一致していなくてもよい。また、図６（ｂ）に示した色材使用量及び指標値は、模式的な数値であり、領域Ｒ１及び出力特性を限定するものではない。むろん、指標値は、部分ＬＵＴ２３０のデータに含まれなくてもよい。

図６（ａ）及び図７に示す少格子点ＬＵＴ２１０は、ＣＭＹＫ色空間の全体にわたって色材ＣＬ１，ＣＬ２の使用量の対応関係を規定しており、最終ＬＵＴ２００の格子点Ｇ１の数Ｎ１＝Ｎｇ⁴よりも少ないＮ３個の色材量予測格子点（第三格子点）Ｇ５を有する。例えば、Ｎｇ＝１７である場合、格子点数はＮ１＝１７⁴＝８３５２１個にもなる。１７⁴個のカラーパッチを用意し測色して１７⁴個の格子点Ｇ１について印刷画像３６０と色再現画像１６０との色味が合うようにｃｍｙｋｌｃｌｍの色材量を求めるのは、非常に労力がかかる。そこで、測色するカラーパッチの数Ｎ２を最終的な格子点数Ｎ１よりも少なくしてＮ２個の色材量予測格子点（第二格子点）Ｇ２，Ｇ３，Ｇ５についてｃｍｙｋｌｃｌｍの色材量を予測し、補間演算により最終的なＮ１個の格子点Ｇ１についてｃｍｙｋｌｃｌｍの色材量を求めてもよい。以下、この思想に従って、説明する。

元ＬＵＴ２０５を構成する色材量予測格子点Ｇ２，Ｇ３，Ｇ５には、少格子点ＬＵＴ２１０の格子点Ｇ５と、部分ＬＵＴ２３０の格子点Ｇ２，Ｇ３と、が含まれている。少格子点ＬＵＴ２１０の格子点Ｇ５は、図７に示すように、包括領域Ｒ９にわたって分散した配置とされている。部分ＬＵＴ２３０の格子点Ｇ２，Ｇ３は、第一の領域Ｒ１及びその近傍において密に配置されている。これは、ユーザーに選択された特定画像４００に含まれるユニークカラーＵＣｉを特に精度良く色変換するためである。

少格子点ＬＵＴ２１０は、部分ＬＵＴ２３０の格子点Ｇ２，Ｇ３及び最終ＬＵＴ２００の格子点Ｇ１よりも低い密度で配置された色材量予測格子点（第三格子点）Ｇ５について、所定の観察光源下で印刷画像３６０と色再現画像１６０との色味が合うように対応関係が規定されている。例えば、Ｄ５０光源下で印刷画像３６０と色再現画像１６０との色味を極力忠実に合わせる場合、印刷画像３６０の色彩値と色再現画像１６０の色彩値との色差ΔＥ_D50が所定の範囲内（例えば１．０など正の閾値以下）となるように対応関係を規定した少格子点ＬＵＴ２１０を用意すればよい。
むろん、色味を合わせる観察光源は、Ｄ５０光源以外の光源でもよいし、複数の光源の組合せでもよい。少格子点ＬＵＴの生成方法の詳細は、後述する。

図６（ｂ）に示す部分ＬＵＴ２３０は、最終ＬＵＴ２００の格子点数Ｎ１よりも少ない数Ｎｐの特有色参照格子点Ｇ２，Ｇ３のそれぞれについて色材ＣＬ１，ＣＬ２の使用量の対応関係を規定している。格子点Ｇ２，Ｇ３は、少格子点ＬＵＴ２１０の格子点Ｇ５よりも高い密度で配置されている。部分ＬＵＴ２３０は領域Ｒ１外の隣接格子点Ｇ３についても色材ＣＬ１，ＣＬ２の使用量の対応関係を規定しているが、これは補間演算により領域Ｒ１内全ての色を変換するためである。従って、格子点Ｇ２，Ｇ３を有する部分ＬＵＴ２３０は、入力色空間に含まれる第一の領域Ｒ１について色材ＣＬ１，ＣＬ２の使用量の対応関係を規定していることになる。

部分ＬＵＴ２３０は、ＣＭＹＫ色空間のうち特定画像４００に含まれるユニークカラーＵＣｉを含む第一の領域Ｒ１について第一の出力特性を示すように対応関係が規定されている。少格子点ＬＵＴ２１０は、ＣＭＹＫ色空間のうち第一の領域Ｒ１とは異なる第二の領域Ｒ２について第一の出力特性とは異なる第二の出力特性を示すように対応関係が規定されている。ＬＵＴ２００を生成するときには、部分ＬＵＴ２３０に基づいてＣＭＹＫ色空間内の第一の領域Ｒ１について第一の出力特性を示すように対応関係を規定し、少格子点ＬＵＴ２１０に基づいてＣＭＹＫ色空間内の第二の領域Ｒ２について第二の出力特性を示すように対応関係を規定する。

図５に示すコンピューター１０は、ＨＤＤ１４に記憶されたＬＵＴ２００をＲＡＭ１２に読み出し、この読み出したＬＵＴ２００を参照して第一の色材ＣＬ１の使用量Ｄｃ，Ｄｍ，Ｄｙ，Ｄｋから第二の色材ＣＬ２の使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmへの色変換を実行可能である。この場合、コンピューター１０に色変換装置が含まれ、ＨＤＤ１４及びＲＡＭ１２が記憶部Ｕ２１を構成し、プログラムデータ１４ａに含まれる色変換プログラムを実行するコンピューター１０が色変換部Ｕ２２を構成する。むろん、色変換装置は、ＬＵＴ生成装置Ｕ０とは別のコンピューターに含まれてもよい。

（４）ルックアップテーブルの生成例：
次に、本技術のＬＵＴ生成方法に従ったＬＵＴ２００の生成例を説明する。

（４−１）ルックアップテーブル生成処理の全体フロー：
図８は、コンピューター１０（ＬＵＴ生成装置Ｕ０）で行われるＬＵＴ生成処理のフローチャートを例示している。図９は、ステップＳ１４で行われる特有色参照格子点抽出処理のフローチャートを例示している。図１０は、ステップＳ２０で行われる最終ＬＵＴ生成処理のフローチャートを例示している。以下、「ステップ」の記載を省略する。
ＬＵＴ生成処理が開始されると、コンピューター１０は、ＬＵＴ生成の条件設定を受け付ける（Ｓ１０）。

図１１は、コンピューター１０がディスプレイ４０に表示させる設定画面５００を例示している。設定画面５００には、色分布を解析する対象の特定画像４００のファイル表示欄５１１，５１２，…、画像選択ボタン５２１，５２２，…、出力特性の選択欄５３１，５３２，…、印刷本機３００の種類選択欄５４１、作成ボタン５５０、等が設けられている。色分布解析対象の特定画像４００を一つに限定する場合には、ファイル表示欄５１２，…と画像選択ボタン５３１，…と出力特性選択端５３２，…を省略してもよい。最終的に生成されるＬＵＴ２００は、被印刷物Ｍ１の種類に合わせられるので、被印刷物Ｍ１の種類選択欄５４２等を設定画面５００に設けてもよい。
印刷本機の選択欄５４１には、オフセット印刷機、グラビア印刷機、といった複数の項目の中から印刷本機の種類を選択入力可能である。被印刷物の選択欄５４２には、コート紙、ケント紙、といった複数の項目の中から被印刷物の種類を選択入力可能である。

ファイル表示欄５１１，５１２，…には、予め用意された色分布解析用の特定画像４００のファイル名が表示される。画像選択ボタン５２１，５２２，…は、複数の特定画像ファイルの中から色分布解析対象の特定画像４００のファイルを選択するための操作ボタンである。画像選択ボタン５２１，５２２，…が選択操作されると、複数の特定画像ファイルが表示され、これらの特定画像ファイルの中から色分布解析対象の特定画像のファイルを選択する操作が行われると、選択された特定画像４００のファイル名がファイル表示欄５１１，５１２，…に表示される。特定画像ファイルの選択入力を受け付けて特定画像４００を読み出す処理は、特定画像４００を選択するＳ１１の処理である。

出力特性選択欄５３１，５３２，…には、「Ｆ１０−Ｆ２」、「Ｆ８−Ｆ２」、「Ｃ−Ｄ５０」、「Ｄ５０」、「粒状性小」、「階調性」、といった項目の中から第一の領域Ｒ１に対応付ける出力特性を選択入力可能である。図１１の例では、「ＡＢＣ．ｔｉｆ」というファイルに格納された特定画像４００に「Ｆ１０−Ｆ２」という第一の出力特性が対応付けられていることが示されている。出力特性選択欄への操作入力を受け付けて選択された項目を記憶する処理は、第一の領域Ｒ１に対応付ける出力特性を選択する入力を受け付けるＳ１２の処理である。このＳ１２の処理が行われることにより、特定画像４００に含まれる色の再現性について主観的な高い色精度が最終ＬＵＴ２００に付与される。

作成ボタン５５０が操作されると、コンピューター１０は、特定画像４００の色分布に基づいてＣＭＹＫ色空間内の複数の格子点Ｇ１のうち特定画像４００に含まれるユニークカラーＵＣｉの変換に用いる特有色参照格子点Ｇ２，Ｇ３を特定する特有色参照格子点抽出処理を行う（Ｓ１４）。Ｓ１６では、図１５に示す印刷色プロファイルＰＲ０（後述）に基づいて、第一の領域Ｒ１について第一の出力特性を示すように特有色参照格子点Ｇ２，Ｇ３について対応関係を規定した部分ＬＵＴ２３０を生成する。Ｓ１８では、印刷色プロファイルＰＲ０に基づいて、第二の領域Ｒ２について第二の出力特性を示すように対応関係を規定した少格子点ＬＵＴ２１０を生成する。その後、コンピューター１０は、部分ＬＵＴ２３０と少格子点ＬＵＴ２１０に基づいて、Ｎ１個の格子点Ｇ１について対応関係を規定したＬＵＴ２００を生成してＨＤＤ１４に保存する最終ＬＵＴ生成処理を行い（Ｓ２０）、ＬＵＴ生成処理を終了させる。

（４−２）特有色参照格子点抽出処理：
図１２は、図９に示す特有色参照格子点抽出処理に従って特有色参照格子点Ｇ２，Ｇ３を抽出する例を模式的に示している。図９に示す特有色参照格子点抽出処理が開始されると、コンピューター１０は、必要に応じて、図８のＳ１１で選択された特定画像４００を縮小して縮小画像４１０を生成する（Ｓ１０２）。この縮小処理は、例えば、特定画像４００の複数の画素ＰＸ１の画素値（Ｃ，Ｍ，Ｙ，Ｋ）を平均する処理とすることができる。図１２では、２×２個の画素ＰＸ１の画素値を例えば相加平均して１個の画素ＰＸ２の画素値（ΣＣ／４，ΣＭ／４，ΣＹ／４，ΣＫ／４）にする様子が示されている。むろん、平均は加重平均等でもよく、縮小率も縦横各５０％に限定されず、縦横の縮小率が異なってもよい。Ｓ１０２の処理は、任意の処理であり、省略可能であるため、例えば、後の処理に時間がかかりすぎなければを省略してもよい。

Ｓ１０４では、特定画像４００から特定画像４００に含まれるユニークカラーＵＣｉを抽出する。例えば、図１２の右上に示すように、縮小画像４１０の各画素ＰＸ２に画素値（Ｃ_i，Ｍ_i，Ｙ_i，Ｋ_i）のユニークカラーＵＣｉ（ｉは１から７の整数）が含まれ、ユニークカラーＵＣ１，ＵＣ７の画素が１個だけであり、ユニークカラーＵＣ２〜ＵＣ６の画素が複数あるものとする。なお、二つのユニークカラーについて画素値のＣ，Ｍ，Ｙ，Ｋのいずれかが異なる場合、これらのユニークカラーは異なる。図１２の例の場合、重複したユニークカラーＵＣ２〜ＵＣ６をそれぞれ一つと数えてユニークカラー数Ｎｃ＝７のユニークカラーＵＣ１〜ＵＣ７が得られる。これらのユニークカラーＵＣ１〜ＵＣ７の一覧は、ユニークカラーリスト（特定画像の色分布）４２０であり、ＲＡＭ１２に格納される。Ｓ１０６では、抽出されたユニークカラーＵＣｉの数Ｎｃを取得する。

Ｓ１０８では、ＣＭＹＫ色空間内における最終ＬＵＴ２００の格子点Ｇ１の座標値を取得する。Ｓ１１０では、ＣＭＹＫ色空間内の複数の格子点Ｇ１のうち、Ｓ１０４で抽出したユニークカラーＵＣｉの変換に用いる特有色参照格子点Ｇ２，Ｇ３を特定する。ＣＭＹＫ色空間内でユニークカラーＵＣｉの位置を内挿する格子点Ｇ１は、ＣＭＹＫそれぞれ２箇所の計２⁴＝１６点となる。図１２の中段には、ＣＭＹＫ色空間の一部を模式的に２次元平面で表し、格子点ＵＧ１〜ＵＧ９を格子点Ｇ１として示し、ユニークカラーＵＣｉの位置を示している。例えば、図１２に示されるユニークカラーＵＣ１を内挿するのは、格子点ＵＧ１，ＵＧ２，ＵＧ４，ＵＧ５を含む１６点である。ユニークカラーＵＣ７を内挿するのは、格子点ＵＧ５，ＵＧ６，ＵＧ８，ＵＧ９を含む１６点である。従って、格子点ＵＧ５は、両ユニークカラーＵＣ１，ＵＣ７を内挿する格子点である。格子点ＵＧ２，ＵＧ３，ＵＧ５，ＵＧ６を含む１６点は、互いに異なるユニークカラーＵＣ２，ＵＣ３を内挿する。Ｓ１１２では、ユニークカラーＵＣｉについて特定した格子点Ｇ１の座標値（Ｃ，Ｍ，Ｙ，Ｋ）をＲＡＭ１２にスタックする。図１２に示す特有色参照格子点スタック４３０は、座標値がスタックされた格子点Ｇ１の一部を格子点ＵＧ１〜ＵＧ９の符号で示している。Ｓ１１４では、ユニークカラーＵＣｉの全てについてＳ１１０〜Ｓ１１２の処理を行ったか否かを判断する。未処理のユニークカラーＵＣｉが残っていれば、コンピューター１０は、処理をＳ１１０に戻す。
以上のようにして、特定画像４００に含まれるユニークカラーＵＣｉのそれぞれについてＣＭＹＫ色空間内の複数の格子点Ｇ１のうちユニークカラーＵＣｉの変換に用いる一以上の特有色参照格子点Ｇ２，Ｇ３が抽出される。

ユニークカラーＵＣｉの全てについてＳ１１０〜Ｓ１１２の処理が行われた場合、コンピューター１０は、特有色参照格子点スタック４３０にスタックされた格子点座標値の中で重複した座標値を一つにまとめる（Ｓ１１６）。図１２の下には、重複が解消された特有色参照格子点Ｇ２，Ｇ３の座標値を格納した特有色参照格子点リスト４４０が示されている。このようにして、Ｓ１１０〜Ｓ１１２で抽出した特有色参照格子点Ｇ２，Ｇ３の重複を解消する。Ｓ１０２〜Ｓ１１６の処理で特定画像４００の色分布（４２０）に基づいてＮｐ個の特有色参照格子点Ｇ２，Ｇ３を特定すると、コンピューター１０は、これらＮｐ個の特有色参照格子点Ｇ２，Ｇ３をＨＤＤ１４に記録して（Ｓ１１８）、特有色参照格子点抽出処理を終了させる。その後行われる図８のＳ１６では、第一の領域Ｒ１について第一の出力特性を示すようにＮｐ個の特有色参照格子点Ｇ２，Ｇ３について対応関係を規定して部分ＬＵＴ２３０を生成する処理を行う。すなわち、特有色参照格子点Ｇ２，Ｇ３は、予測部Ｕ１０でｃｍｙｋｌｃｌｍの色材量を予測する色材量予測格子点である。部分ＬＵＴ生成処理の詳細は、後述する。

（４−３）少格子点ルックアップテーブルの構造：
図８のＳ１８で行われる少格子点ＬＵＴ生成処理では、図６（ａ）及び図７で示したように、最終ＬＵＴ２００よりも格子点の少ない少格子点ＬＵＴ２１０におけるＮ３個の格子点（第三格子点）Ｇ５について、第一及び第二の色材ＣＬ１，ＣＬ２の使用量を対応付けて少格子点ＬＵＴ２１０を生成する処理を行う。格子点Ｇ５は、部分ＬＵＴ２３０の特有色参照格子点Ｇ２，Ｇ３とともに、予測部Ｕ１０でｃｍｙｋｌｃｌｍの色材量を予測する色材量予測格子点である。なお、特有色参照格子点Ｇ２，Ｇ３の範囲Ｒ１ｇに格子点Ｇ５がある場合、該格子点Ｇ５については色材量の予測を行わなくてもよく、前記範囲Ｒ１ｇに特有色参照格子点Ｇ２，Ｇ３の予測色材量を用いて最終ＬＵＴ２００を生成することにしている。格子点Ｇ５の数Ｎ３は、前記範囲Ｒ１ｇにある格子点を除いた数とする。
上述したように、ＣＭＹＫ色空間は４次元であるため、ＣＭＹＫそれぞれの使用量についてＮｇ段階の格子点を設ける場合、最終ＬＵＴ２００の格子点数はＮ１＝Ｎｇ⁴個となり、Ｎｇ＝１７であればＮ１＝８３５２１個にもなる。そこで、少格子点ＬＵＴ２１０の格子点Ｇ５は、最終ＬＵＴ２００の格子点数Ｎ１よりも少ないＮ３個にしている。

少格子点ＬＵＴ２１０の格子点Ｇ５の数Ｎ３は、ＣＭＹＫそれぞれの色材使用量についてＮｓ段階（Ｎｓは２≦Ｎｓ＜Ｎｇを満たす整数）としてＮｓ⁴個にする考えもある。この思想も本技術に含まれるが、ＣＭＹＫのいずれか１色（選択色とする。）の色材使用量を固定して残りの３色で形成される３次元色空間のガマットは選択色の色材使用量に応じて変わるので、Ｎ３＜Ｎｓ⁴個にしてもよい。ガマットとは、色再現領域の大きさ若しくは体積のことであり、例えば、Ｌ^*ａ^*ｂ^*色空間内の大きさ若しくは体積で表される。例えば、Ｋの色材ＣＬ３は、明度及び彩度を下げる効果をもたらす色材であり、（Ｃ，Ｍ，Ｙ，Ｋ）＝（５０％，３０％，１０％，０％）のＣＭＹの色に対して、（Ｃ，Ｍ，Ｙ，Ｋ）＝（５０％，３０％，１０％，１００％）のＣＭＹＫの色は、暗く、彩度が低く再現されることになる。このようなことから、Ｋの色材使用量をある一つの有限値（例えばＫ＝１０％）に固定したＣＭＹＫ色空間のガマットは、Ｋを全く使わないＣＭＹ色空間のガマットよりも小さくなる。Ｋの色材を一定量加えることで、加えないときよりもガマットが小さくなるとも言える。
以上のことから、Ｋの色材使用量が多くて小さいガマットしか持たないＣＭＹ色空間に対しては、多くの格子点を持たせる必要が無い。逆に、Ｋを使用しないＣＭＹ色空間等、Ｋの色材使用量が少なくてガマットが大きいＣＭＹ色空間に対しては、離散的な格子点で色空間内の色を正確に記述するために多くの格子点を持たせる必要がある。そこで、図１３の例のように、最終ＬＵＴ２００の格子点Ｇ１に対応するＫの色材ＣＬ３の使用量のそれぞれ、すなわち、格子点Ｇ１のＫ軸の各位置について、Ｋの色材ＣＬ３の使用量に応じた数の格子点Ｇ５をＣＭＹの色材ＣＬ４で表されるＣＭＹ色空間内に設けることにしている。具体的には、Ｋの色材使用量が少ないときにはガマットが大きいのでＣＭＹ色空間内の格子点の数を多くし、Ｋの色材使用量が多いときにはガマットが小さいのでＣＭＹ色空間内の格子点の数を少なくすることにしている。

図１３は、Ｋの色材ＣＬ３の使用量毎の色材量予測格子点（第三格子点）Ｇ５の配置例を模式的に示している。最終ＬＵＴ２００の格子点Ｇ１のＫの色材使用量がＮｇ段階あるとき、ＣＭＹ３次元色空間はＮｇ個となる。図１３の例では、Ｋの色材使用量が０％である場合、ＣＭＹ色空間のガマットは、比較的広い範囲である。Ｋの色材使用量が２５％である場合、ＣＭＹ色空間のガマットは、Ｋ＝０％の場合よりも狭い範囲である。そこで、Ｋ＝０％（第一の使用量の例）に対応する第三格子点Ｇ５の数はＫ＝２５％（第二の使用量の例）に対応する第三格子点Ｇ５の数よりも多くしている。また、Ｋの色材使用量が５０％である場合、ＣＭＹ色空間のガマットは、Ｋ＝２５％の場合よりも狭い範囲である。従って、Ｋ＝０％及びＫ＝２５％（第一の使用量の例）に対応する第三格子点Ｇ５の数はＫ＝５０％（第二の使用量の例）に対応する第三格子点Ｇ５の数よりも多くしている。さらに、Ｋの色材使用量が７５％、１００％の順に、ＣＭＹ色空間のガマットがさらに狭くなる。従って、Ｋの色材使用量が多くなるほどＣＭＹ色空間の第三格子点Ｇ５の数を減らしている。なお、補間演算のためにＣＭＹ色空間に最低限２³＝８個の格子点Ｇ５を設けることにしている。

なお、Ｋにおける第一及び第二の使用量は、相対的なものであるため、例えば、第一の使用量としてＫ＝７５％とし第二の使用量としてＫ＝１００％としてもよい。むろん、各ＣＫＹ色空間に設けられる格子点Ｇ５は、均等に配置されてもよいが、不均等に配置されてもよい。

最終的な格子点Ｇ１のＫの各位置についてＣＭＹの格子点Ｇ５の各位置を決めると、図７の例のように少格子点ＬＵＴ２１０となる色材量予測格子点Ｇ５の配置が決まったことになる。色材量予測格子点Ｇ５が特有色参照格子点（部分ＬＵＴ２３０の色材量予測格子点）Ｇ２，Ｇ３の範囲Ｒ１ｇにある場合、該範囲Ｒ１ｇの色材量予測格子点Ｇ５については色材量の予測を行わなくてもよい。図８のＳ１８では、第二の領域Ｒ２について第二の出力特性を示すようにＮ３個の色材量予測格子点Ｇ５について対応関係を規定して少格子点ＬＵＴ２１０を生成する処理を行う。少格子点ＬＵＴ生成処理の詳細は、後述する。

（４−４）印刷色プロファイルの生成例：
図１４，１５は、第一の色材ＣＬ１の使用量とターゲット色彩値との対応関係を規定した印刷色プロファイルＰＲ０を示している。図１４では、部分ＬＵＴ２３０生成用の印刷色プロファイルと少格子点ＬＵＴ２１０生成用の印刷色プロファイルをまとめて示している。すなわち、図１４に示す印刷色プロファイルＰＲ０は、部分ＬＵＴ２３０となる特有色参照格子点（色材量予測格子点）Ｇ２，Ｇ３と、少格子点ＬＵＴ２１０となる色材量予測格子点Ｇ５と、を有している。色材量予測格子点Ｇ２，Ｇ３，Ｇ５は、最終ＬＵＴ２００に設けられる第一格子点Ｇ１の数Ｎ１よりも少ない数Ｎ２の第二格子点Ｇ６であり、最終ＬＵＴ２００生成用の元ＬＵＴ２０５の格子点Ｇ６である。

印刷色プロファイルＰＲ０は、第一の色材ＣＬ１の使用量Ｄｃ，Ｄｍ，Ｄｙ，Ｄｋと、印刷物３５０に形成される前記使用量Ｄｃ，Ｄｍ，Ｄｙ，Ｄｋの第一の色材ＣＬ１の観察光源Ｌ０下におけるターゲット色彩値Ｌ_D，ａ_D，ｂ_Dと、の対応関係が規定されている。複数の観察光源Ｌ０のそれぞれについてターゲット色彩値に近似するように第二の色材ＣＬ２の使用量を予測する場合、各観察光源Ｌ０について印刷色プロファイルＰＲ１，ＰＲ２，ＰＲ３，…を設けるとよい。

ここで、印刷色プロファイルＰＲ０に規定されるターゲット色彩値、及び、色再現画像１６０に形成される第二の色材ＣＬ２の色彩値は、デバイス非依存の色空間（機器独立色空間）、且つ、均等色空間の色彩値が好ましいものの、デバイス依存の色空間（機器従属色空間）の色彩値にすることも可能であるし、均等色空間でない色空間の色彩値にすることも可能である。デバイス非依存の均等色空間は、ＣＩＥ Ｌ^*ａ^*ｂ^*色空間の他、ＣＩＥ Ｌ^*ｕ^*ｖ^*色空間等でもよい。Ｌ^*ａ^*ｂ^*色空間のＬ^*は明度を表し、ａ^*及びｂ^*は色相及び彩度を示す色度を表す。

図１６は、少格子点ＬＵＴ２１０及び部分ＬＵＴ２３０を生成するための印刷色プロファイルＰＲ０を生成する方法を模式的に例示している。
印刷本機３００は、図１６に例示するように、印刷色プロファイルＰＲ０を作成するためのチャート印刷物３５１を形成する。印刷物３５１は、被印刷物Ｍ１にカラーチャート３６１を形成したものである。チャート３６１には、複数のパッチ３６２が２次元に並べられている。パッチ３６２は、色票とも呼ばれ、測色機８００による測色の単位となる領域を意味する。各パッチ３６２は、印刷色プロファイルＰＲ０の色材量予測格子点Ｇ２，Ｇ３，Ｇ５に対応している。従って、チャート３６１には、Ｎ２個のパッチ３６２が配置される。観察光源Ｌ０毎に印刷色プロファイルＰＲ０を生成する場合、同じ印刷物３５１を用い観察光源を切り替えてチャート３６１を測色すればよい。例えば、Ｄ５０，Ｄ５５，Ｄ６５，Ａ，Ｆ１１光源と多くの観察光源下でＣＭＹＫ色空間の印刷画像３６０とＣＭＹＫｌｃｌｍ色空間の色再現画像１６０との色味を極力合わせる少格子点ＬＵＴ２１０を生成する場合、Ｄ５０，Ｄ５５，Ｄ６５，Ａ，Ｆ１１光源（擬似的な光源を含む。以下、同様。）を順次切り替えてチャート３６１を測色し、Ｄ５０，Ｄ５５，Ｄ６５，Ａ，Ｆ１１光源のそれぞれについて第一の色材の使用量Ｄｃ，Ｄｍ，Ｄｙ，Ｄｋとターゲット色彩値Ｌ_D，ａ_D，ｂ_Dとを対応付けた色プロファイルＰＲ０を生成すればよい。

例えば、観察光源Ｌ０にＤ５０光源Ｌ１を用い各パッチ３６２を測色して複数組の色彩値Ｌ_D-D50，ａ_D-D50，ｂ_D-D50をそれぞれ第一の色材の使用量Ｄｃ，Ｄｍ，Ｄｙ，Ｄｋに対応付けると、図１５に示す印刷色プロファイルＰＲ１を生成することができる。この処理は、印刷色プロファイルＰＲ１のＮ４個の格子点全てについて行ってもよいが、格子点Ｇ１の一部について複数組の色彩値Ｌ_D-D50，ａ_D-D50，ｂ_D-D50をそれぞれ第一の色材の使用量Ｄｃ，Ｄｍ，Ｄｙ，Ｄｋに対応付けた後、残りの格子点については補間演算により第一の色材の使用量Ｄｃ，Ｄｍ，Ｄｙ，Ｄｋを補間してもよい。印刷色プロファイルＰＲ１は、例えば、出力特性「Ｄ５０」の部分ＬＵＴ２３０又は少格子点ＬＵＴ２１０を生成するために使用される。

また、観察光源Ｌ０に例えばＦ１０光源を用い各パッチ３６２を測色して複数組の色彩値Ｌ_D-F10，ａ_D-F10，ｂ_D-F10をそれぞれ第一の色材の使用量Ｄｃ，Ｄｍ，Ｄｙ，Ｄｋに対応付けると、図１５に示す印刷色プロファイルＰＲ２を生成することができる。観察光源Ｌ０に例えばＦ２光源を用い各パッチ３６２を測色して複数組の色彩値Ｌ_D-F2，ａ_D-F2，ｂ_D-F2をそれぞれ第一の色材の使用量Ｄｃ，Ｄｍ，Ｄｙ，Ｄｋに対応付けると、図１５に示す印刷色プロファイルＰＲ３を生成することができる。これらの印刷色プロファイルＰＲ２，ＰＲ３の生成にも補間演算を用いてもよい。印刷色プロファイルＰＲ２，ＰＲ３は、例えば、出力特性「Ｆ１０−Ｆ２」の部分ＬＵＴ２３０又は少格子点ＬＵＴ２１０を生成するために使用される。

（４−５）第二の色材の使用量を予測する例：
ＬＵＴ生成装置Ｕ０は、決定した光源数分のターゲット色彩値を同時に再現するインク量をインク量の最適量探索手法（最適化アルゴリズム）で算出することができる。予測部Ｕ１０は、複数の観察光源Ｌ０について印刷色プロファイルＰＲ１に規定されたターゲット色彩値Ｌ_D，ａ_D，ｂ_Dへの近似性を評価する評価値Ｉに基づいて、色再現画像１６０に形成されるＣＭＹＫｌｃｌｍの色材ＣＬ２の色彩値が観察光源毎にターゲット色彩値Ｌ_D，ａ_D，ｂ_Dに近似するようにＣＭＹＫｌｃｌｍの色材ＣＬ２の使用量を予測する。

図１７は、ある分光反射率を有するターゲット（パッチ）から複数の観察光源下でそれぞれ色彩値を得る様子を模式的に説明している。ターゲットの分光反射率Ｒ_t（λ）は、通常、全可視波長領域において均一でない分布を有している。各光源は、それぞれ異なる分光エネルギーＰ（λ）の分布を有している。ターゲットに光源を照射したときの各波長の反射光の分光エネルギーは、ターゲット分光反射率Ｒ_t（λ）と分光エネルギーＰ（λ）を各波長について掛け合わせた値となる。さらに、反射光の分光エネルギーのスペクトルに対して人間の分光感度特性に応じた等色関数ｘ（λ），ｙ（λ），ｚ（λ）をそれぞれ畳み込み積分し、係数ｋによって正規化することにより、３刺激値Ｘ，Ｙ，Ｚを得ることができる。

上記３刺激値Ｘ，Ｙ，Ｚを所定の変換式によって変換することにより、色彩値Ｌ^*ａ^*ｂ^*を得ることができる。

図１７に示すように、光源毎に分光エネルギーＰ（λ）のスペクトルが異なるので、最終的に得られるターゲット色彩値も光源に応じて異なることとなる。

図１８は、ターゲット色彩値Ｌ_D，ａ_D，ｂ_Dと同様の色を再現可能なインク量セットφを算出する最適インク量算出モジュール群の処理の流れを模式的に例示している。第二の色材ＣＬ２がＣＭＹＫｌｃｌｍの色材である場合、インク量セットφは、吐出するＣＭＹＫｌｃｌｍインクの使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmの組合せを意味する。

最適インク量算出モジュール群（予測部Ｕ１０）は、インク量セット算出モジュール（ＩＣＭ）Ｐ３ａ１、分光反射率予測モジュール（ＲＰＭ）Ｐ３ａ２、色算出モジュール（ＣＣＭ）Ｐ３ａ３、評価値算出モジュール（ＥＣＭ）Ｐ３ａ４、を備える。
インク量セット算出モジュール（ＩＣＭ）Ｐ３ａ１は、入力がＣＭＹＫの色材使用量である４次元の印刷色プロファイルＰＲ０から一つの格子点Ｇ１を選択し、該格子点Ｇ１に対応付けられているターゲット色彩値Ｌ_D，ａ_D，ｂ_Dを取得する。この点は、入力がＲＧＢの画像出力である特開2009-200820号公報記載の印刷システムと大きく異なる。
分光反射率予測モジュール（ＲＰＭ）Ｐ３ａ２は、ＩＣＭ Ｐ３ａ１からのインク量セットφ、具体的にはインク使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmの入力に応じて、当該インク量セットφに基づいてプルーフ機１００が印刷用紙といった被印刷物Ｍ２にインクを吐出させたときの分光反射率Ｒ（λ）を予測分光反射率Ｒ_s（λ）として予測する。インク量セットφを指定すれば被印刷物Ｍ２上における各インクドットの形成状態が予測できるため、ＲＰＭ Ｐ３ａ２は一意に予測分光反射率Ｒ_s（λ）を算出することができる。

ここで、図２０〜２３を参照して、ＲＰＭ Ｐ３ａ２に用いられる予測モデル（分光プリンティングモデル）を説明する。図２０には、プルーフ機１００の記録ヘッド２１を模式的に例示している。記録ヘッド２１は、ＣＭＹＫｌｃｌｍのインク毎に複数のノズル２１ａを備えている。プルーフ機１００は、ＣＭＹＫｌｃｌｍのインク毎の使用量をインク量セットφ（ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lm）にする制御を行う。各ノズル２１ａから吐出されたインク滴は被印刷物Ｍ２上において多数のドットの集まりとなり、これによってインク量セットφ（ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lm）に応じたインク被覆率の色再現画像１６０が被印刷物Ｍ２上に形成される。

ＲＰＭ Ｐ３ａ２に用いられる予測モデル（分光プリンティングモデル）は、任意のインク量セットφ（ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lm）で印刷を行った場合の分光反射率Ｒ（λ）を予測分光反射率Ｒ_s（λ）として予測するものである。分光プリンティングモデルにおいては、インク量空間における複数の代表点についてカラーパッチを印刷し、その分光反射率Ｒ（λ）を分光反射率計によって測定することにより得られた分光反射率データベースＲＤＢを用意する。この分光反射率データベースＲＤＢを使用したセル分割ユール・ニールセン分光ノイゲバウアモデル(Cellular Yule-Nielsen Spectral Neugebauer Model)による予測を行うことにより、任意のインク量セットφで印刷を行った場合の予測分光反射率Ｒ_s（λ）を正確に予測する。

図２１は、分光反射率データベースＲＤＢの構造を模式的に例示している。本実施形態のインク量空間は６次元であるが、図の簡略化のためＣＭ面のみ図示している。分光反射率データベースＲＤＢは、インク量空間における複数の格子点のインク量セット（ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lm）について実際に印刷及び測定をして得られた分光反射率Ｒ（λ）が記述されたＬＵＴとなっている。このＬＵＴは、各インク量軸を複数に分割する格子点を有している。なお、一部の格子点のみ実際に印刷及び測定をし、他の格子点については実際に印刷及び測定を行った格子点の分光反射率Ｒ（λ）に基づいて分光反射率Ｒ（λ）を予測してもよい。これにより、実際に印刷及び測定を行うカラーパッチの個数を低減させることができる。

分光反射率データベースＲＤＢは、被印刷物の種類毎に用意される。分光反射率Ｒ（λ）は被印刷物上に形成されたインク膜（ドット）による分光透過率と被印刷物の反射率によって決まるものであり、被印刷物の表面物性（ドット形状が依存）や反射率の影響を大きく受けるからである。

ＲＰＭ Ｐ３ａ２は、ＩＣＭ Ｐ３ａ１の要請に応じて分光反射率データベースＲＤＢを使用したセル分割ユール・ニールセン分光ノイゲバウアモデルによる予測を実行する。この予測にあたっては、ＩＣＭ Ｐ３ａ１から予測条件を取得し、この予測条件を設定する。例えば、被印刷物やインク量セットφを印刷条件として設定する。光沢紙を印刷用紙として予測を行う場合には、光沢紙にカラーパッチを印刷することにより作成した分光反射率データベースＲＤＢが設定される。

分光反射率データベースＲＤＢの設定ができると、ＩＣＭ Ｐ３ａ１から入力されたインク量セットφ（ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lm）を分光プリンティングモデルに適用する。セル分割ユール・ニールセン分光ノイゲバウアモデルは、よく知られた分光ノイゲバウアモデルとユール・ニールセンモデルとに基づいている。簡略化のためＣＭＹの３種類のインクを用いた場合のモデルについて説明するが、同様のモデルを本実施形態のＣＭＹＫｌｃｌｍのインクセットを用いたモデルに拡張することができる。セル分割ユール・ニールセン分光ノイゲバウアモデルについては、Color Ｒes Appl 25, 4-19, 2000及びＲ Balasubramanian, Optimization of the spectral Neugebauer model for printer characterization, J. Electronic Imaging 8(2), 156-166 (1999)を参照。

図２２は、分光ノイゲバウアモデルを模式的に例示している。分光ノイゲバウアモデルでは、任意のインク量セット（ｄ_c，ｄ_m，ｄ_y）で印刷したときの印刷物の予測分光反射率Ｒ_s（λ）が以下の式で与えられる。

ここで、ａ_iはｉ番目の領域の面積率であり、Ｒ_i（λ）はｉ番目の領域の分光反射率である。添え字ｉは、インクの無い領域（ｗ）と、シアンインクのみの領域（ｃ）と、マゼンタインクのみの領域（ｍ）と、イエローインクのみの領域（ｙ）と、マゼンタインクとイエローインクが吐出される領域（ｒ）と、イエローインクとシアンインクが吐出される領域（ｇ）と、シアンインクとマゼンタインクが吐出される領域（ｂ）と、ＣＭＹの３つのインクが吐出される領域（ｋ）をそれぞれ意味している。また、ｆ_c，ｆ_m，ｆ_yは、ＣＭＹ各インクを１種類のみ吐出したときにそのインクで覆われる面積の割合（「インク被覆率(Ink area coverage)」と呼ぶ。）である。

インク被覆率ｆ_c，ｆ_m，ｆ_yは、図２２（Ｂ）に示すマーレイ・デービスモデルで与えられる。マーレイ・デービスモデルでは、例えばシアンインクのインク被覆率ｆ_cは、シアンのインク量ｄ_cの非線形関数であり、例えば１次元ルックアップテーブルによってインク量ｄ_cをインク被覆率ｆ_cに換算することができる。インク被覆率ｆ_c，ｆ_m，ｆ_yがインク量ｄ_c，ｄ_m，ｄ_yの非線形関数となる理由は、単位面積に少量のインクが吐出された場合にはインクが十分に広がるが、多量のインクが吐出された場合にはインクが重なり合うためにインクで覆われる面積があまり増加しないためである。他の種類のＭＹインクについても同様である。

分光反射率に関するユール・ニールセンモデルを適用すると、前記（２）式は以下の（３ａ）式又は（３ｂ）式に書き換えられる。

ここで、ｎは１以上の所定の係数であり、例えばｎ＝１０に設定することができる。（３ａ）式及び（３ｂ）式は、ユール・ニールセン分光ノイゲバウアモデル(Yule-Nielsen Spectral Neugebauer Model)を表す式である。

本実施形態で採用するセル分割ユール・ニールセン分光ノイゲバウアモデル(Cellular Yule-Nielsen Spectral Neugebauer Model)は、上述したユール・ニールセン分光ノイゲバウアモデルのインク量空間を複数のセルに分割したものである。

図２３（Ａ）は、セル分割ユール・ニールセン分光ノイゲバウアモデルにおけるセル分割の例を示している。ここでは、説明の簡略化のため、ＣＭインクのインク量ｄ_c，ｄ_mの２つの軸を含む２次元インク量空間でのセル分割を描いている。インク被覆率ｆ_c，ｆ_mは、上述したマーレイ・デービスモデルにてインク量ｄ_c，ｄ_mと一意の関係にあるため、インク被覆率ｆ_c，ｆ_mを示す軸と考えることもできる。白丸は、セル分割のグリッド点（「格子点」と呼ぶ。）であり、２次元のインク量（被覆率）空間が９つのセルＣ１〜Ｃ９に分割されている。各格子点に対応するインク量セット（ｄ_c，ｄ_m）は、分光反射率データベースＲＤＢに規定された格子点に対応するインク量セットとされている。すなわち、上述した分光反射率データベースＲＤＢを参照することにより、各格子点の分光反射率Ｒ（λ）を得ることができる。従って、各格子点の分光反射率Ｒ（λ）₀₀，Ｒ（λ）₁₀，Ｒ（λ）₂₀・・・Ｒ（λ）₃₃は、分光反射率データベースＲＤＢから取得することができる。

本実施形態では、セル分割がＣＭＹＫｌｃｌｍの６次元インク量空間で行われ、各格子点の座標も６次元のインク量セットφ（ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lm）によって表される。各格子点のインク量セットφに対応する格子点の分光反射率Ｒ（λ）は、分光反射率データベースＲＤＢ（例えばコート紙のもの）から取得されることとなる。

図２３（Ｂ）は、セル分割モデルにて使用するインク被覆率ｆ_cとインク量ｄ_cとの関係を示している。ここでは、１種類のインクのインク量の範囲０〜ｄ_cmaxも３つの区間に分割され、区間毎に０から１まで単調に増加する非線形の曲線によってセル分割モデルにて使用する仮想的なインク被覆率ｆ_cが求められる。他のインクについても同様にインク被覆率ｆ_m，ｆ_yが求められる。

図２３（Ｃ）は、図２３（Ａ）の中央のセルＣ５内にある任意のインク量セット（ｄ_c，ｄ_m）にて印刷を行った場合の予測分光反射率Ｒ_s（λ）の算出方法を示している。インク量セット（ｄ_c，ｄ_m）にて印刷を行った場合の予測分光反射率Ｒ_s（λ）は、以下の式で与えられる。

ここで、（４）式におけるインク被覆率ｆ_c，ｆ_mは図２３（Ｂ）のグラフで与えられる値である。また、セルＣ５を囲む４つの格子点に対応する分光反射率Ｒ（λ）₁₁，（λ）₁₂，（λ）₂₁，（λ）₂₂は分光反射率データベースＲＤＢを参照することにより取得することができる。これにより、（４）式の右辺を構成するすべての値を確定することができ、その計算結果として任意のインク量セットφ（ｄ_c，ｄ_m）にて印刷を行った場合の予測分光反射率Ｒ_s（λ）を算出することができる。波長λを可視波長域にて順次シフトさせていくことにより、可視光領域における予測分光反射率Ｒ_s（λ）を得ることができる。インク量空間を複数のセルに分割すれば、分割しない場合に比べて予測分光反射率Ｒ_s（λ）をより精度良く算出することができる。
以上のようにして、ＲＰＭ Ｐ３ａ２は、ＩＣＭ Ｐ３ａ１の要請に応じて予測分光反射率Ｒ_s（λ）を予測する。

予測分光反射率Ｒ_s（λ）が得られると、色算出モジュール（ＣＣＭ）Ｐ３ａ３は当該予測分光反射率Ｒ_s（λ）の物体に複数の観察光源Ｌ０を照射したときの予測色彩値を算出する。この予測色彩値には、例えば、ＣＩＥ Ｌ^*ａ^*ｂ^*色空間のＬ^*ａ^*ｂ^*値が用いられる。予測色彩値を算出する流れは、図１７及び上記（１）式に示したものと同様である。

（５）式に示すように、予測分光反射率Ｒ_s（λ）に各光源の分光エネルギーのスペクトルを乗算し、等色関数による畳み込み積分をし、さらに３刺激値をＬ^*ａ^*ｂ^*値に変換することにより、予測色彩値Ｌ_d，ａ_d，ｂ_dが求まる。予測色彩値は観察光源Ｌ０毎に算出される。

評価値算出モジュール（ＥＣＭ）Ｐ３ａ４は、各観察光源についてターゲット色彩値Ｌ_D，ａ_D，ｂ_Dと予測色彩値Ｌ_d，ａ_d，ｂ_dとの色差ΔＥを算出する。色差は、ΔＥ＝｛（Ｌ_d−Ｌ_D）²＋（ａ_d−ａ_D）²＋（ｂ_d−ｂ_D）²｝^1/2により算出してもよいし、ＣＩＥ ＤＥ２０００の色差式（ΔＥ₂₀₀₀）に基づいて算出してもよい。観察光源Ｌ０としてＤ５０，Ｄ６５，Ａ光源が選択された場合、各光源の色差をΔＥ_D50，ΔＥ_D65，ΔＥ_Aと表記することにする。ターゲット色彩値Ｌ_D，ａ_D，ｂ_Dへの近似性を評価する評価値Ｉ（φ）は、インク使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmに依存する評価関数であり、例えば、以下の式によって算出することができる。

ここで、ｉは、観察光源を表している。上述した例では、ｉ＝１はＤ５０光源を表し、ｉ＝２はＤ６５光源を表し、ｉ＝３はＡ光源を表す。Ｎは、選択された観察光源の数を表している。ΔＥ_iは、観察光源ｉの下でターゲット色彩値Ｌ_D，ａ_D，ｂ_Dと予測色彩値Ｌ_d，ａ_d，ｂ_dとの色差を表している。ｗ_iは、各観察光源下での色差ΔＥ_iに対する重みを表している。本実施形態では重みｗ_iが均等であるものとして説明するが、均等でなくてもよい。

評価値Ｉ（φ）は、各色差ΔＥ_iが小さくなると小さくなり、ターゲット色彩値と予測色彩値とが各観察光源において総合的に近いほど小さくなる性質を有している。ＩＣＭＰ３ａ１がインク量セットφをＲＰＭ Ｐ３ａ２とＣＣＭ Ｐ３ａ３とＥＣＭ Ｐ３ａ４に出力することにより、最終的に評価値Ｉ（φ）がＩＣＭ Ｐ３ａ１に返される。ＩＣＭ Ｐ３ａ１は、インク量セットφに対応する評価値Ｉ（φ）を算出することを繰り返すことにより、目的関数としての評価値Ｉ（φ）が極小化するようなインク量セットφの最適解を算出する。この最適解を算出する手法としては、例えば勾配法といった非線形最適化手法を用いることができる。

図１９は、ＣＩＥ Ｌ^*ａ^*ｂ^*色空間において、各観察光源下のターゲット色彩値Ｌ_D，ａ_D，ｂ_Dと、インク量セットφが最適化されていく際の各観察光源下の予測色彩値Ｌ_d，ａ_d，ｂ_dの推移を示している。インク量セットφ（ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lm）は、各色差ΔＥ_iが次第に小さくなるように最適化されていく。このようにして、観察光源毎にターゲット色彩値Ｌ_D，ａ_D，ｂ_Dに近似する色彩値の色を色再現画像１６０に再現させることが可能なインク量セットφが算出される。

最適化処理の終了条件は、例えば、色差ΔＥ_iが閾値ＴＥ_i以下となったときとすることができる。閾値ＴＥ_iは、各観察光源ｉに設けられ、同じ値のみならず、互いに異なる値でもよい。全ての観察光源ｉについてΔＥ_i≦ＴＥ_iとなると、最適化処理が終了する。このようにして、予測部Ｕ１０は、第二の色材ＣＬ２の使用量を予測する。部分ＬＵＴ生成部Ｕ３は、ＣＭＹＫ色空間のうち特定画像４００に含まれるユニークカラーＵＣｉを含む第一の領域Ｒ１において第一の出力特性を示すような第二の色材ＣＬ２の予測使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmを第一の色材ＣＬ１の使用量Ｄｃ，Ｄｍ，Ｄｙ，Ｄｋに対応付けて部分ＬＵＴ２３０を生成する。少格子点ＬＵＴ生成部Ｕ４は、ＣＭＹＫ色空間のうち第一の領域Ｒ１とは異なる第二の領域Ｒ２において第二の出力特性を示すような第二の色材ＣＬ２の予測使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmを第一の色材ＣＬ１の使用量Ｄｃ，Ｄｍ，Ｄｙ，Ｄｋに対応付けて少格子点ＬＵＴ２１０を生成する。

以下の具体例を説明するとき、第一の領域Ｒ１がＦ１０，Ｆ２光源下で両画像３６０，１６０の色味が合う「Ｆ１０−Ｆ２」の第一の出力特性を有し、第二の領域Ｒ２がＤ５０光源下で両画像３６０，１６０の色味が合う「Ｄ５０」の第二の出力特性を有するものとする。この場合、Ｆ１０，Ｆ２光源が第一の観察光源となり、Ｄ５０光源が第二の観察光源となる。Ｆ１０，Ｆ２光源下の色差ΔＥ₁，ΔＥ₂に対する閾値ＴＥ₁，ＴＥ₂は、混合比率１：１を示すＴＥ₁＝ＴＥ₂とする。第二の出力特性は、Ｄ５０光源のメタメリックマッチングを優先しながら多数の観察光源（例えばＤ５０，Ｄ５５，Ｄ６５，Ａ，Ｆ１１光源）の下で両画像３６０，画像１６０の色味が合う出力特性でもよい。この場合、Ｄ５０光源下の色差ΔＥ_D50については他の光源下の色差よりも小さくなるようにすればよい。

（４−６）最終ルックアップテーブル生成処理：
図１０に示す最終ＬＵＴ生成処理では、色材量予測格子点Ｇ２，Ｇ３，Ｇ５である元ＬＵＴ２０５のＮ２個の第二格子点Ｇ６（図２，７，１４参照）についての対応関係に基づいてＮ１個（Ｎ１＞Ｎ２）の第一格子点Ｇ１について対応関係を規定して最終ＬＵＴ２００を生成する。なお、第一の領域Ｒ１については第一の出力特性が示されるように、特有色参照格子点Ｇ２，Ｇ３の範囲Ｒ１ｇの格子点Ｇ５は元ＬＵＴ２０５に含まれていないものとする。
図１０に示す最終ＬＵＴ生成処理が開始されると、コンピューター１０は、まず、最終ＬＵＴ２００の色材使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmの格納領域をＲＡＭ１２に確保し、最終ＬＵＴ２００の格子点Ｇ１の座標値（Ｃ_i，Ｍ_i，Ｙ_i，Ｋ_i）を取得する（Ｓ２０２）。Ｓ２０４では、Ｋの色材ＣＬ３の使用量のそれぞれについて、ＣＭＹの色材ＣＬ４で表されるＣＭＹ３次元色空間をドローネー三角分割法（Delaunay Triangulation）により元ＬＵＴ２０５の格子点Ｇ６を頂点とする三角形に分割する。Ｋの使用量がＮｇ段階あると、図２，７に示すようなＣＭＹ３次元色空間のＮｇ個についてＣＭＹ３次元色空間において四面体補間を行うための三角形分割を行うことになる。基本的には、分割後の三角形の最小角が最も大きくなるように三角形分割を行う。Ｓ２０６では、最終ＬＵＴ２００の格子点Ｇ１の数Ｎ１を取得する。ここで、最終ＬＵＴ２００の各格子点Ｇ１を識別する変数をｉとする。ｉは、１からＮ１までの整数となる。

図２４は、Ｓ２０８〜Ｓ２１０の処理で元ＬＵＴ２０５の格子点Ｇ６から最終ＬＵＴ２００の格子点Ｇ１を補間する例を模式的に示している。各格子点Ｇ６には、第二の色材ＣＬ２の使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmが格納されている。Ｓ２０８では、最終ＬＵＴ２００の格子点Ｇ１ｉを内挿する元ＬＵＴ２０５の四面体ＴＥを特定する。図２４中、元ＬＵＴ２０５中の８個の格子点Ｇ６ｈはＣＭＹ色空間において格子点Ｇ１ｉを囲む最小の六面体の頂点を示し、元ＬＵＴ２０５中の４個の格子点Ｇ６ｔはＣＭＹ色空間を三角形分割したときに８個の格子点Ｇ６ｈの中で格子点Ｇ１ｉを囲む最小の四面体ＴＥの頂点を示している。この四面体ＴＥが格子点Ｇ１ｉを内挿する四面体である。Ｓ２１０では、４個の格子点Ｇ６ｔにおける第一の色材ＣＬ１の使用量Ｄｃ，Ｄｍ，Ｄｙ，Ｄｋと第二の色材ＣＬ２の使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmとの対応関係を用いて公知の四面体補間（補間演算）を行い、格子点Ｇ１ｉについて第二の色材ＣＬ２の使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmを算出する。このようにして、Ｋの色材使用量のそれぞれについて、ＣＭＹ色空間において第一格子点Ｇ１を囲む複数の第二格子点Ｇ６についてのＫとＣＭＹの色材使用量の対応関係を用いた補間演算を行うことにより第一格子点Ｇ１について最終ＬＵＴ２００の対応関係を規定していく。

Ｓ２１２では、最終ＬＵＴ２００の格子点Ｇ１ｉに色材使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmを格納する。むろん、格子点Ｇ１ｉの位置が元ＬＵＴの格子点Ｇ６の位置と一致する場合には、Ｓ２０８〜Ｓ２１０の処理を行わず格子点Ｇ６の色材使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmをそのまま格子点Ｇ１ｉに格納してもよい。Ｓ２１４では、格子点Ｇ１ｉの全てについてＳ２０８〜Ｓ２１２の処理を行ったか否かを判断する。未処理の格子点Ｇ１ｉが残っていれば、コンピューター１０は、処理をＳ２０８に戻す。
以上のようにして、ＣＭＹＫ色空間において最終ＬＵＴの第一格子点Ｇ１を囲む元ＬＵＴの複数の第二格子点Ｇ６についての色材ＣＬ１，ＣＬ２の使用量の対応関係を用いた補間演算を行うことにより第一格子点Ｇ１について最終ＬＵＴ２００の対応関係を規定する。

格子点Ｇ１ｉの全てについてＳ２０８〜Ｓ２１２の処理が行われた場合、コンピューター１０は、必要に応じて、ＣＭＹＫ色空間の中で第一の領域Ｒ１の境界部分Ｒｂ（図２参照）に規定される色材使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmをスムージングする（Ｓ２１６）。Ｓ２１６のスムージング処理は、任意の処理であり、省略可能である。

図２５は、上記スムージング処理を模式的に例示している。分かり易く示すため、横軸にＣＭＹＫ色空間のＣの入力使用量Ｄｃ、すなわち、格子点のＣの座標を示し、縦軸にｃｍｙｋｌｃｌｍ色空間のｃの出力使用量ｄ_c、すなわち、格子点に格納されるｃの第二の色材ＣＬ２の使用量を示している。むろん、横軸にはＣＭＹＫ色空間のＭ，Ｙ，Ｋの入力使用量Ｄｍ，Ｄｙ，Ｄｋを適用することが可能であり、縦軸にはｃｍｙｋｌｃｌｍ色空間のｍ，ｙ，ｋ，ｌｃ，ｌｍの出力使用量ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmを適用することが可能である。図２５の例では、境界部分Ｒｂにおいて、第一の領域Ｒ１内の格子点ＧＲ１における出力使用量ｄ_cと、第二の領域Ｒ２内の格子点ＧＲ２における出力使用量ｄ_cと、が滑らかに繋がっていないことが示されている。このような場合、格子点ＧＲ１，ＧＲ２の出力使用量Ｄｃ，ｄ_cをスムージングして滑らかに繋げると、色再現画像１６０の画質を向上させることができる。

スムージング処理は、例えば、第一の色材ＣＬ１の使用量から対応する第二の色材ＣＬ２の使用量を求める近似の関数ｆ（例えばｄ_c＝ｆ（Ｄｃ））を最小自乗法等により求め、第二の色材ＣＬ２の使用量を関数ｆの値に置き換える処理とすることができる。スムージング処理は、入力使用量Ｄｃに対する出力使用量ｄ_c、入力使用量Ｄｍに対する出力使用量ｄ_m、入力使用量Ｄｙに対する出力使用量ｄ_y、入力使用量Ｄｋに対する出力使用量ｄ_k、等と色別に行ってもよい。また、関数ｆを使用せず、第一の色材ＣＬ１の使用量（例えばＤｃ）の変化に対する第二の色材ＣＬ２の使用量（例えばｄ_c）の変化が極力変動しないように境界部分Ｒｂの第二の色材ＣＬ２の使用量（例えばｄ_c）を変える処理を行ってもよい。
第一の領域Ｒ１の境界部分Ｒｂに規定される第二の色材ＣＬ２の使用量がスムージングされると、ユーザーが意図する色再現目標が好適に実現される。

ＬＵＴ格納領域の全格子点に第二の色材ＣＬ２の使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmが格納されると、Ｎ１個の格子点Ｇ１について対応関係を規定した最終ＬＵＴ２００が生成されたことになる。コンピューター１０は、ＲＡＭ１２に格納されている最終ＬＵＴ２００をＨＤＤ１４に保存し（Ｓ２１８）、最終ＬＵＴ生成処理を終了させる。

上述した具体例に従って生成した最終ＬＵＴ２００は、図３に示すように、ＣＭＹＫ色空間のうち特定画像４００に含まれるユニークカラーＵＣｉを含む第一の領域Ｒ１について例えば「Ｆ１０−Ｆ２」の第一の出力特性を示すように対応関係が規定され、ＣＭＹＫ色空間のうち第一の領域Ｒ１とは異なる第二の領域Ｒ２について第一の出力特性とは異なる例えば「Ｄ５０」の第二の出力特性を示すように対応関係が規定されている。第一の領域Ｒ１の色を変換する格子点の出力使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmは、部分ＬＵＴ２３０に規定された出力使用量に基づく。第二の領域Ｒ２の色を変換する格子点の出力使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmは、少格子点ＬＵＴ２１０に規定された出力使用量に基づく。

以上のことから、図３に示すように、最終ＬＵＴ２００に含まれる第一の領域Ｒ１は、１：１のＦ１０，Ｆ２光源下で印刷画像３６０の色彩値と色再現画像１６０の色彩値との色差ΔＥ_F10-F2が正の閾値ＴＨ１以下（第一の範囲内）となる第一の出力特性を有する。最終ＬＵＴ２００に含まれる第二の領域Ｒ２は、Ｄ５０光源下で両画像３６０，１６０の色彩値の色差ΔＥ_D50が正の閾値ＴＨ２以下（第二の範囲内）となる第二の出力特性を有する。閾値ＴＨ１，ＴＨ２は、特に限定されないが、例えば１．０〜２．０程度とすることができ、異なる値でもよい。なお、第一の範囲内にはΔＥ_F10-F2＜ＴＨ１とすることも含まれ、第二の範囲内にはΔＥ_D50＜ＴＨ２とすることも含まれる。
第一の領域Ｒ１では、各格子点Ｇ２，Ｇ３についてΔＥ_F10-F2≦ＴＨ１であることを重視した結果、ΔＥ_D50＞ＴＨ２（第二の範囲外）となる格子点が存在する。ただし、ΔＥ_D50≦ＴＨ２（第二の範囲内）となる格子点が一部に存在してもよい。第二の領域Ｒ２では、各格子点Ｇ２，Ｇ３についてΔＥ_D50≦ＴＨ２であることを重視した結果、ΔＥ_F10-F2＞ＴＨ１（第一の範囲内）となる格子点が存在する。ただし、ΔＥ_F10-F2≦ＴＨ１（第一の範囲内）となる格子点が一部に存在してもよい。

すなわち、図３に示すように、第一の領域Ｒ１の出力使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmは、第一の出力特性を表すΔＥ_F10-F2≦ＴＨ１（第一の範囲内）となるように設定され、少なくとも一部がΔＥ_D50＞ＴＨ２（第二の範囲外）となるように設定される。第二の領域Ｒ２の出力使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmは、第二の出力特性を表すΔＥ_D50≦ＴＨ２（第二の範囲内）となるように設定され、少なくとも一部がΔＥ_F10-F2＞ＴＨ１（第一の範囲外）となるように設定される。

以上説明したように、一つのＬＵＴ２００に、ＣＭＹＫ色空間内のうち、特定画像４００に含まれる色（ＵＣｉ）を含む第一の領域Ｒ１と、該領域Ｒ１とは異なる第二の領域Ｒ２と、で異なる出力特性を示すように対応関係が規定される。従って、印刷本機３００のＣＭＹＫの色材ＣＬ１に合わせた色材使用量Ｄｃ，Ｄｍ，Ｄｙ，Ｄｋとプルーフ機１００のｃｍｙｋｌｃｌｍの色材ＣＬ２に合わせた色材使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmとの対応関係を規定した一つのＬＵＴ２００で特定画像４００の色に対する高精度の色再現が実現される。その結果、例えば、特定画像に含まれる色の再現性については主観的な高い色精度が実現され、特定画像に含まれない色の再現性については客観的な高い色精度が実現される。

また、元ＬＵＴ２０５の格子点Ｇ６のうち、特定画像に含まれる色を変換する格子点Ｇ２，Ｇ３については格子点Ｇ５よりも高い密度で配置されているので、この点でも特定画像の色再現を高精度で行うことができる。
一方、元ＬＵＴ２０５の格子点Ｇ５については最終ＬＵＴ２００の格子点Ｇ１よりも低い密度で配置されているので、測色するカラーパッチの数が少なくて済み、ｃｍｙｋｌｃｌｍの色材量を予測する格子点の数が少なくて済む。従って、本技術は、最終ＬＵＴの生成を容易にすることができる。

なお、出力特性は、領域Ｒ１，Ｒ２内の全体的な出力の傾向を意味する。このため、最終ＬＵＴ２００に含まれる第一の領域Ｒ１内の複数箇所の色における印刷画像３６０と色再現画像１６０の色差ΔＥ_F10-F2の相加平均といった平均が正の閾値ＴＨ１’以下（第一の範囲内）となるように設定され、両画像３６０，１６０の色差ΔＥ_D50の平均が正の閾値ＴＨ２’よりも大（第二の範囲外）となるように設定されてもよい。また、最終ＬＵＴ２００に含まれる第二の領域Ｒ２内の複数箇所の色における両画像３６０，１６０の色差ΔＥ_D50の相加平均といった平均が閾値ＴＨ２’以下（第二の範囲内）となるように設定され、両画像３６０，１６０の色差ΔＥ_F10-F2の平均が閾値ＴＨ１’よりも大（第一の範囲外）となるように設定されてもよい。

（４−７）色再現画像を形成する例：
図２０は、上記ＬＵＴ２００を記憶したプルーフ機（画像形成装置）１００で行われる色再現画像出力制御処理を例示している。この処理は、例えば、色再現画像１６０の形成要求をプルーフ機１００が受け付けたときに開始する。
処理が開始すると、プルーフ機１００は、印刷本機３００で印刷物３５０の形成に使用されるＣＭＹＫの色材ＣＬ１の使用量Ｄｃ，Ｄｍ，Ｄｙ，Ｄｋに対応した入力データを取得する（Ｓ３０２）。必要に応じて、Ｓ３０４では、色変換処理に使用するＬＵＴ２００を選択する。

Ｓ３０６では、ＬＵＴ２００を参照して、ＣＭＹＫの色材使用量Ｄｃ，Ｄｍ，Ｄｙ，Ｄｋのデータを、プルーフ機１００で色再現画像１６０の形成に使用されるＣＭＹＬｌｃｌｍの色材ＣＬ２の使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmのデータに変換する。色材使用量Ｄｃ，Ｄｍ，Ｄｙ，Ｄｋに一致するＬＵＴ２００の入力点が無い場合には、ＣＭＹＫ色空間において色材使用量Ｄｃ，Ｄｍ，Ｄｙ，Ｄｋに近隣する複数の入力点のそれぞれに対応したＣＭＹＬｌｃｌｍの色材使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmを用いてＤｃ，Ｄｍ，Ｄｙ，Ｄｋに対応するｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmを補間すればよい。このようにして、ＬＵＴ２００に従ってＣＭＹＫの色材使用量Ｄｃ，Ｄｍ，Ｄｙ，ＤｋをＣＭＹＬｌｃｌｍの色材使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmに変換することができる。

Ｓ３０８では、色材使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmのデータに対してハーフトーン処理を行い、使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmに対応したドットの形成状況を表す多値データを生成する。ハーフトーン処理には、ディザ法、誤差拡散法、濃度パターン法、等を用いることができる。Ｓ３１０では、得られたドットの形成状況を表す多値データを記録ヘッド２１の各走査パス及び各ノズル２１ａに割り振って出力制御データを生成する。記録ヘッド２１は、前記出力制御データに従って被印刷物Ｍ２にインクドットを形成し、色再現画像１６０を有するプルーフ１５０を形成する。

ここで、印刷画像３６０を１：１のＦ１０，Ｆ２光源下で測色機により測色して得られる色彩値をＬ_D1，ａ_D1，ｂ_D1とし、色再現画像１６０を１：１のＦ１０，Ｆ２光源下で測色機により測色して得られる色彩値をＬ_d1，ａ_d1，ｂ_d1とする。両画像３６０，１６０の第一の領域Ｒ１（特定画像のユニークカラーＵＣｉを含む領域）については、Ｌ_D1，ａ_D1，ｂ_D1とＬ_d1，ａ_d1，ｂ_d1から求まる色差ΔＥ_F10-F2がＴＨ１以下（第一の範囲内）となる。両画像３６０，１６０の第二の領域Ｒ２については、色差ΔＥ_F10-F2の少なくとも一部がＴＨ１よりも大（第一の範囲外）となることがある。
また、印刷画像３６０をＤ５０光源下で測色機により測色して得られる色彩値をＬ_D2，ａ_D2，ｂ_D2とし、色再現画像１６０をＤ５０光源下で測色機により測色して得られる色彩値をＬ_d2，ａ_d2，ｂ_d2とする。両画像３６０，１６０の第二の領域Ｒ２については、Ｌ_D2，ａ_D2，ｂ_D2とＬ_d2，ａ_d2，ｂ_d2から求まる色差ΔＥ_D50がＴＨ２以下（第二の範囲内）となる。両画像３６０，１６０の第一の領域Ｒ１（特定画像のユニークカラーＵＣｉを含む領域）については、色差ΔＥ_D50の少なくとも一部がＴＨ２よりも大（第二の範囲外）となることがある。

以上説明したように、一つのＬＵＴ２００に、入力色空間内の第一の領域Ｒ１については第一の出力特性を示すように対応関係が規定され、入力色空間内の第二の領域Ｒ２については第一の出力特性とは異なる第二の出力特性を示すように対応関係が規定される。従って、本技術は、印刷本機３００の色材ＣＬ１に合わせた色材使用量とプルーフ機１００の色材ＣＬ２に合わせた色材使用量との対応関係を規定した一つのＬＵＴ２００で特定画像４００の色に対する高精度の色再現を実現することができる。

（４−８）出力特性が色再現画像の粒状性を表す粒状性指標値ＩＮ１で表される例：
入力色空間に含まれる領域について示される出力特性は、上述してきた色差ΔＥで表される以外にも、粒状性指標値ＩＮ１や階調性指標値ＩＮ２等でも表される。
粒状性指標値ＩＮ１は、小さくなるほど色再現画像１６０の粒状感が弱くなる指標値である。写真画像等では一般に粒状感が弱くなるように色変換をすると画質を向上させることができる。一方、グラフ画像等では鮮やかな色を出すため色変換としては粒状感が強くなる方向にすると画質が向上する可能性がある。

図２６は、色再現画像１６０の粒状性を変える様子を模式的に例示している。
色再現画像１６０の粒状性を変えた最終ＬＵＴは、例えば、上述したＬＵＴ２００の出力使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmを変えることにより生成することができる。図２６に示す例は、分かり易く示すため、ｋの１ドットがｃの１ドットとｍの１ドットとｙの１ドットの組合せに等価である場合にｋの使用量ｄ_kの少なくとも一部とｃｍｙの使用量の少なくとも一部とを置き換える様子を示している。第一の領域Ｒ１又は第二の領域Ｒ２のｋの使用量ｄ_kの少なくとも一部Δｄをｃｍｙの使用量に置き換えると、第一の領域Ｒ１又は第二の領域Ｒ２について、粒状性指標値ＩＮ１の小さい、すなわち、比較的明度の低いｋのドットが少ないか無くて比較的明度の高いｃｍｙのドットが多い粒状感の小さい色再現画像１６０が形成される最終ＬＵＴ２５１を生成することができる。また、第一の領域Ｒ１又は第二の領域Ｒ２のｃｍｙの使用量ｄ_c，ｄ_m，ｄ_yの少なくとも一部Δｄをｋの使用量に置き換えると、第一の領域Ｒ１又は第二の領域Ｒ２について、粒状性指標値ＩＮ１の大きい、すなわち、比較的明度の高いｃｍｙのドットが少ないか無くて比較的明度の低いｋのドットが多い粒状感の大きい色再現画像１６０が形成される最終ＬＵＴ２５２を生成することができる。

色再現画像１６０の粒状感を変えるためには、ｃの１ドットとｌｃの２以上のドットとを置き換えたり、ｍの１ドットとｌｍの２以上のドットとを置き換えたりしてもよい。第一の領域Ｒ１又は第二の領域Ｒ２のｃの使用量ｄ_cの少なくとも一部Δｄをｌｃの使用量に置き換えたりｍの使用量ｄ_mの少なくとも一部Δｄをｌｍの使用量に置き換えたりすると、第一の領域Ｒ１又は第二の領域Ｒ２について、粒状性指標値ＩＮ１の小さい、すなわち、比較的明度の低いｃｍのドットが少ないか無くて比較的明度の高いｌｃｌｍのドットが多い粒状感の小さい色再現画像１６０が形成されるＬＵＴ２５１を生成することができる。
また、粒状性指標値ＩＮ１に基づいてＬＵＴ２５１，２５２を生成してもよい。

色再現画像１６０の粒状性を定量化する粒状性指標値ＩＮ１は、例えば、以下の式で表すことができる。

ここで、ａ_Lは明度補正項、ＷＳ（ｕ）は画像のウイナースペクトラム、ＶＴＦは視覚の空間周波数特性、ｕは空間周波数、である。粒状性指標値ＩＮ１の思想は、例えば、Makoto Fujino,Image Quality Evaluation of Inkjet Prints, Japan Hardcopy '99, p.291-294に示されている。なお、（７）式は一次元で示しているが、空間周波数ｕ，ｖの関数として二次元画像の空間周波数を算出することは容易である。

粒状性指標値ＩＮ１は、ある色再現画像を観察者が視認したときに、その観察者が感じる粒状感（あるいはノイズの程度）である。ＩＮ１が小さい程、観察者が感じる粒状感は小さくなる。むろん、ＩＮ１は、色再現画像を形成したときの粒状性を評価する指標値であればよく、他の式を用いることも可能である。粒状性指標値ＩＮ１を求める詳細は、例えば、特表2007-516663号公報に示されている。

粒状性指標値ＩＮ１は、例えば、図２０のＳ３０８のハーフトーン処理により得られる色再現画像１６０形成用の多値データ、すなわち、色材使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmに対応したドットの形成状況を表す多値データに基づいて算出することができる。この多値データがあれば、実際に色再現画像１６０を被印刷物Ｍ２に形成しなくてもドットの形成状況をシミュレートすることによってＩＮ１を求めることができる。従って、ＬＵＴ２００を構成する各格子点Ｇ１に対応付けられた出力使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmのベタの画像（巨視的に一様な画像）に対して前記のハーフトーン処理を行って得られる多値データからＩＮ１を求めると、各格子点Ｇ２，Ｇ３に粒状性指標値ＩＮ１を対応付けることができる。むろん、被印刷物Ｍ２に形成された色再現画像１６０をスキャナーといった画像読み取り装置等で読み取り、該読み取った情報に基づいて色再現画像１６０形成用の多値データを取得し、この多値データに基づいてＩＮ１を算出してもよい。

図２７は、粒状性指標値ＩＮ１を算出する様子を模式的に例示している。まず、色再現画像１６０形成用の多値データに基づいて、色再現画像１６０の各画素について明度Ｌ（ｘ，ｙ）を求める。座標（ｘ，ｙ）にドットが形成されない場合には、被印刷物Ｍ２自体の明度がＬ（ｘ，ｙ）となる。座標（ｘ，ｙ）にドットが形成される場合、画素に対するドットの面積比と、ドットの部分の明度と、被印刷物Ｍ２に自体の明度と、に基づいて（例えば明度の比例配分）、明度Ｌ（ｘ，ｙ）を算出する。図２７の左端に示すような明度Ｌ（ｘ，ｙ）が求まると、明度Ｌ（ｘ，ｙ）に対してＦＦＴ（Fast Fourier Transformation）を実施し、空間周波数のスペクトルＳ（ｕ，ｖ）を得る。スペクトルＳ（ｕ，ｖ）は実部Ｒｅ（ｕ，ｖ）と虚部Ｉｍ（ｕ，ｖ）とからなり、Ｓ（ｕ，ｖ）＝Ｒｅ（ｕ，ｖ）＋ｊＩｍ（ｕ，ｖ）である。このスペクトルＳ（ｕ，ｖ）は、上述のウイナースペクトラムに相当する。

次に、座標（ｕ，ｖ）に対応した空間周波数（cycle/mm）の大きさｆ（ｕ，ｖ）を算出するため、シミュレート対象となっている色再現画像１６０の最低周波数を算出する。この最低周波数は、シミュレート対象の多値データで形成される色再現画像１６０において１回振動する周波数であり、主走査方向（Ｘ方向）と副走査方向（Ｙ方向）とのそれぞれについて定義される。例えば、１インチを２５．４mmとして、主走査方向の最低周波数ｅ_uはＸ解像度／（Ｘ方向の画素数×２５．４）、副走査方向の最低周波数ｅ_vはＹ解像度／（Ｙ方向の画素数×２５．４）と定義される。任意の座標（ｕ，ｖ）における空間周波数の大きさｆ（ｕ，ｖ）は、｛（ｅ_u・ｕ）²＋（ｅ_v・ｖ）²｝^1/2として算出可能である。

一方、人間の目は、空間周波数の大きさｆ（ｕ，ｖ）に応じて明度に対する感度が異なる。当該視覚の空間周波数特性は、例えば、図２７の中央下部に示すＶＴＦ（ｆ）のような特性を示す。このＶＴＦ（ｆ）は、ＶＴＦ（ｆ）＝5.05×exp（−0.138・ｄ・π・ｆ／１８０）×｛１−exp（−0.1・ｄ・π・ｆ／１８０）｝である。ここで、ｄは色再現画像と目の距離であり、ｆは空間周波数の大きさｆである。このｆは（ｕ，ｖ）の関数として表現されているので、視覚の空間周波数特性ＶＴＦは（ｕ，ｖ）の関数ＶＴＦ（ｕ，ｖ）とすることができる。ＶＴＦが得られたら、当該ＶＴＦの２乗とスペクトルＳ（ｕ，ｖ）の２乗とを乗じ、全画素（ｕ，ｖ）について積分を行う。得られる積分値（Ｐｏｗとする。）の平方根を全画素数で除して規格化すると、元の多値データの画素数に依存しない客観的な指数（Ｉｎｔとする。）が算出される。
Ｉｎｔ＝Ｐｏｗ^1/2／全画素数
なお、多値データの画素数を常に同じ画素数として粒状性を評価するのであれば、必ずしも規格化しなくてもよい。

粒状性指標値ＩＮ１は、色再現画像１６０全体の明度による影響を考慮してもよい。空間周波数のスペクトルが同じであっても色再現画像全体が明るい場合と暗い場合とでは人間の目に異なった印象を与えるので、全体が明るい方が粒状性を感じやすいものとして補正を行ってもよい。このため、例えば、色再現画像１６０の全画素の明度の平均（Ａｖｅとする。）を算出し、色再現画像全体の明るさによる補正係数ａ（Ｌ）＝｛（Ａｖｅ＋１６）／１１６｝^0.8を上記Ｉｎｔに乗じて粒状性指標値ＩＮ１としてもよい。
ＩＮ１＝ａ（Ｌ）・Ｉｎｔ
補正係数ａ（Ｌ）は、上述の明度補正項ａ_Lに相当する。

以上のことから、図１１で示した設定画面５００において、第一の領域Ｒ１について「粒状性小」を選択してもよい。例えば、ＬＵＴ２５１の第一の領域Ｒ１については、色再現画像１６０の粒状性指標値ＩＮ１が正の閾値ＴＨ３以下（第三の範囲内）となる「粒状性小」の第一の出力特性を有するものとする。ＬＵＴ２５１の第二の領域Ｒ２については、例えばＤ５０光源下の両画像３６０，１６０の色彩値の色差を所定範囲内にすることを重視し粒状性を重視してない結果、色再現画像１６０の粒状性指標値ＩＮ１が閾値ＴＨ３よりも大（第三の範囲外）となることがある第二の出力特性を有するものとする。なお、第三の範囲内にはＩＮ１＜ＴＨ３とすることも含まれる。

図３を参照して説明すると、ＬＵＴ２５１に含まれる第一の領域Ｒ１の出力使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmは、ＩＮ１≦ＴＨ３（第三の範囲内）となる「粒状性小」の第一の出力特性を示すように設定される。ＬＵＴ２５１に含まれる第二の領域Ｒ２の出力使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmは、少なくとも一部がＩＮ１＞ＴＨ３（第三の範囲外）となる第二の出力特性を示すように設定される。

ここで、特定画像４００に含まれる色（ＵＣｉ）からＬＵＴ２５１に従って変換された色再現画像１６０の粒状性指標値ＩＮ１は、ＴＨ３以下（第三の範囲内）となる。また、特定画像４００に含まれない色からＬＵＴ２５１に従って変換された色再現画像１６０の粒状性指標値ＩＮ１は、ＴＨ３よりも大（第三の範囲外）となることがある。
なお、第一の色材ＣＬ１の使用量（例えばＤｃ，Ｄｍ，Ｄｙ，Ｄｋ）を徐々に増加又は減少させたグラデーションを有する色再現画像１６０を形成したとき、色再現画像１６０の中である濃度範囲内の領域の粒状性指標値ＩＮ１がＴＨ３以下となり、別の濃度範囲内の領域で粒状性指標値ＩＮ１がＴＨ３よりも大きくなれば、「粒状性小」の出力特性を示す領域と「粒状性小」以外の出力特性を示す別の領域とが一つのＬＵＴ２５１に存在することになる。

以上説明したように、一つのＬＵＴ２５１に、入力色空間内のある領域については粒状性指標値ＩＮ１が第三の範囲内となるように対応関係が規定され、入力色空間内の別の領域については「粒状性小」とは異なる出力特性を示すように対応関係が規定される。従って、本技術は、印刷本機３００の色材ＣＬ１に合わせた色材使用量とプルーフ機１００の色材ＣＬ２に合わせた色材使用量との対応関係を規定した一つのＬＵＴ２５１で入力色空間の領域毎にユーザーが意図する粒状性の色再現目標を実現することができる。

（４−９）出力特性が色再現画像の階調性を表す階調性指標値ＩＮ２で表される例：
階調性指標値ＩＮ２は、小さくなるほど色再現画像１６０の階調変化が滑らかになる指標値である。写真画像等では一般に階調変化が滑らかになるように色変換をすると画質を向上させることができる。一方、グラフ画像等では鮮やかな色を出すため色変換としては階調変化が粗くなる方向にすると画質が向上する可能性がある。

図２８は、色再現画像１６０の階調性を変える様子を模式的に例示している。
色再現画像１６０の階調性を変えた最終ＬＵＴは、例えば、上述したＬＵＴ２００の出力使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmを変えることにより生成することができる。図２８に示す例は、分かり易く示すため、横軸にＣＭＹＫ色空間のＣの入力使用量Ｄｃ、すなわち、格子点のＣの座標を示し、縦軸にｃｍｙｋｌｃｌｍ色空間のｃの出力使用量ｄ_c、すなわち、格子点に格納されるｃの第二の色材ＣＬ２の使用量を示している。むろん、横軸にはＣＭＹＫ色空間のＭ，Ｙ，Ｋの入力使用量Ｄｍ，Ｄｙ，Ｄｋを適用することが可能であり、縦軸にはｃｍｙｋｌｃｌｍ色空間のｍ，ｙ，ｋ，ｌｃ，ｌｍの出力使用量ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmを適用することが可能である。図２８中、最終ＬＵＴ２５３，２５４に示される白丸は、元のＬＵＴ２００における出力使用量ｄ_cを示している。図２８の例では、Ｄｃ軸に等間隔に入力使用量変数Ｄ_i（ｉは１からｎまでの整数、ｎは３以上の整数）を設け、Ｄ_iに対応する出力使用量関数ｄ_c（Ｄ_i）を示している。

階調変化が滑らかであることは、出力使用量関数ｄ_c（Ｄ_i）の微小変化Δｄ_c（Ｄ_i）＝｛ｄ_c（Ｄ_i+1）−ｄ_c（Ｄ_i）｝／ΔＤの変動が小さいことを意味する。また、階調変化が粗いことは、微小変化Δｄ_c（Ｄ_i）の変動が大きいことを意味する。そこで、階調性指標値ＩＮ２は、微小変化Δｄ_c（Ｄ_i）の微小変化ΔΔｄ_c（Ｄ_i）＝｛Δｄ_c（Ｄ_i+1）−Δｄ_c（Ｄ_i）｝／ΔＤの絶対値の平均と定義してもよい。
ＩＮ２＝Σ｛｜ΔΔｄ_c（Ｄ_i）｜｝／（ｎ−２）
ただし、Σはｉ＝１からｉ＝ｎ−２までの｛｜ΔΔｄ_c（Ｄ_i）｜｝の総和を意味する。

得られる階調性指標値ＩＮ２は、小さくなるほど階調変化が滑らかになることを意味し、大きくなるほど階調変化が粗くなることを意味する。従って、第一の色材ＣＬ１の使用量（例えばＤｃ）の変化に対する第二の色材ＣＬ２の使用量（例えばｄ_c）の変化が極力変動しないように第二の色材ＣＬ２の使用量（例えばｄ_c）を変えると、階調性指標値ＩＮ２の小さい、すなわち、階調変化が滑らかな色再現画像１６０が形成される最終ＬＵＴ２５３を生成することができる。また、第一の色材ＣＬ１の使用量（例えばＤｃ）の変化に対する第二の色材ＣＬ２の使用量（例えばｄ_c）の変化がより変動するように第二の色材ＣＬ２の使用量（例えばｄ_c）を変えると、階調性指標値ＩＮ２の大きい、すなわち、階調変化が粗い色再現画像１６０が形成される最終ＬＵＴ２５４を生成することができる。
また、階調性指標値ＩＮ２に基づいてＬＵＴ２５３，２５４を生成してもよい。

以上のことから、図１１で示した設定画面５００において、第一の領域Ｒ１について「階調性」を選択してもよい。例えば、ＬＵＴ２５３の第一の領域Ｒ１については、色再現画像１６０の階調性指標値ＩＮ２が正の閾値ＴＨ４以下（第四の範囲内）となる「階調性」の第一の出力特性を有するものとする。ＬＵＴ２５３の第二の領域Ｒ２については、例えばＤ５０光源下の両画像３６０，１６０の色彩値の色差を所定範囲内にすることを重視し階調性を重視してない結果、色再現画像１６０の階調性指標値ＩＮ２が閾値ＴＨ４よりも大（第四の範囲外）となることがある第二の出力特性を有するものとする。なお、第四の範囲内にはＩＮ２＜ＴＨ４とすることも含まれる。

図３を参照して説明すると、ＬＵＴ２５３に含まれる第一の領域Ｒ１の出力使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmは、ＩＮ２≦ＴＨ４（第四の範囲内）となる「階調性」の第一の出力特性を示すように設定される。ＬＵＴ２５３に含まれる第二の領域Ｒ２の出力使用量ｄ_c，ｄ_m，ｄ_y，ｄ_k，ｄ_lc，ｄ_lmは、少なくとも一部がＩＮ２＞ＴＨ４（第四の範囲外）となる第二の出力特性を示すように設定される。

ここで、特定画像４００に含まれる色（ＵＣｉ）からＬＵＴ２５３に従って変換された色再現画像１６０の階調性指標値ＩＮ２は、ＴＨ４以下（第四の範囲内）となる。また、特定画像４００に含まれない色からＬＵＴ２５３に従って変換された色再現画像１６０の階調性指標値ＩＮ２は、ＴＨ４よりも大（第四の範囲外）となることがある。
なお、第一の色材ＣＬ１の使用量（例えばＤｃ，Ｄｍ，Ｄｙ，Ｄｋ）を徐々に増加又は減少させたグラデーションを有する色再現画像１６０を形成したとき、色再現画像１６０の中である濃度範囲内の領域の階調性指標値ＩＮ２がＴＨ４以下となり、別の濃度範囲内の領域で階調性指標値ＩＮ２がＴＨ４よりも大きくなれば、「階調性」の出力特性を示す領域と「階調性」以外の出力特性を示す別の領域とが一つのＬＵＴ２５３に存在することになる。

以上説明したように、一つのＬＵＴ２５３に、入力色空間内のある領域については階調性指標値ＩＮ２が第四の範囲内となるように対応関係が規定され、入力色空間内の別の領域については「階調性」とは異なる第二の出力特性を示すように対応関係が規定される。従って、本技術は、印刷本機３００の色材ＣＬ１に合わせた色材使用量とプルーフ機１００の色材ＣＬ２に合わせた色材使用量との対応関係を規定した一つのＬＵＴ２５３で入力色空間の領域毎にユーザーが意図する階調性の色再現目標を実現することができる。

（５）変形例：
本発明は、種々の変形例が考えられる。
例えば、ＬＵＴ２００を記憶しているコンピューターにプリンターといった記録装置が接続され、前記ＬＵＴ２００に従って第一の色材の使用量から第二の色材の使用量を求め、求めた第二の色材の使用量となるように記録装置で色再現画像を形成する場合、コンピューターと記録装置が画像形成装置を構成する。画像形成装置は、記録装置以外にも、ディスプレイといった画像出力装置等が含まれる。色再現画像には、画像出力装置の画面上の画像が含まれる。

印刷機で印刷物の形成に使用される第一の色材は、ＣＭＹＫの組合せ以外にも、５色の組合せ、３色の組合せ、等でもよい。画像形成装置で色再現画像の形成に使用される第二の色材は、ＣＭＹＫｌｃｌｍの組合せ以外にも、７色以上の組合せでもよいし、５色以下の組合せでもよい。第二の色材として使用可能な色材の色には、ｃｍｙｋｌｃｌｍ以外にも、ｏｒ（オレンジ）、ｇｒ（グリーン）、ｂ（ブルー）、ｖ（バイオレット）、ｄｙ（ダークイエロー）、ｌｋ（ライトブラック）、ｌｌｋ（ライトライトブラック）、無着色、等が含まれる。無着色の色材には、被印刷物に光沢を付与する色材、有色の色材を保護する色材、等が含まれる。
第一の色材に含まれる第三の色材と第四の色材の組合せは、ＫとＣＭＹの組合せ以外にも、Ｋの濃淡の組合せ、Ｃの濃淡の組合せ、Ｍの濃淡の組合せ、Ｙの濃淡の組合せ、等でもよい。
上述した処理は、順番を入れ替える等、適宜、変更可能である。例えば、図８のＬＵＴ生成処理において、Ｓ１６の部分ＬＵＴ生成処理の前にＳ１８の少格子点ＬＵＴ生成処理を行ってもよい。

（６）結び：
以上説明したように、本発明によると、種々の態様により、印刷機の色材に合わせた色材使用量と画像形成装置の色材に合わせた色材使用量との対応関係を規定した一つのルックアップテーブルで特定画像の色に対する高精度の色再現を実現することが可能な技術等を提供することができる。むろん、従属請求項に係る構成要件を有しておらず独立請求項に係る構成要件のみからなる技術等でも、上述した基本的な作用、効果が得られる。
また、上述した実施形態及び変形例の中で開示した各構成を相互に置換したり組み合わせを変更したりした構成、公知技術並びに上述した実施形態及び変形例の中で開示した各構成を相互に置換したり組み合わせを変更したりした構成、等も実施可能である。本発明は、これらの構成等も含まれる。

１０…コンピューター（色変換装置）、１００…プルーフ機（画像形成装置）、１５０…プルーフ、１６０…色再現画像、２００，２５１〜２５４…ルックアップテーブル、２０５…元ルックアップテーブル、２１０…少格子点ルックアップテーブル、２３０…部分ルックアップテーブル、３００…印刷本機（印刷機）、３５０，３５１…印刷物、３６０…印刷画像、３６１…チャート、３６２…パッチ、４００…特定画像、４０１…ポートレート画像、４０２…グレイ画像、４０３…肌色画像、４１０…縮小画像、４２０…ユニークカラーリスト（特定画像の色分布）、４３０…特有色参照格子点スタック、４４０…特有色参照格子点リスト、５００…設定画面、８００…測色機、ＣＬ１…第一の色材、ＣＬ２…第二の色材、ＣＬ３…第三の色材、ＣＬ４…第四の色材、Ｇ１…第一格子点、Ｇ２，Ｇ３…特有色参照格子点、Ｇ５…第三格子点、Ｇ６…第二格子点、Ｌ０，Ｌ１，Ｌ２，Ｌ３…観察光源、Ｍ１，Ｍ２…被印刷物、ＰＲ０〜ＰＲ３…印刷色プロファイル、Ｒ１…第一の領域、Ｒ２…第二の領域、Ｒ９…包括領域、Ｒ１ｇ…領域の色の変換に用いる格子点の範囲、Ｒｂ…境界部分、ＳＹ１…プルーフシステム、ＵＣｉ…ユニークカラー（特定画像に含まれる色）。

Claims (14)

1. 印刷機で印刷物の形成に使用される第一の色材の使用量と、画像形成装置で色再現画像の形成に使用される第二の色材の使用量と、の対応関係を規定したルックアップテーブルを生成するルックアップテーブル生成方法であって、
   前記第一の色材の色で表される入力色空間のうち特定画像に含まれる色を含む第一の領域について第一の出力特性を示すように前記対応関係を規定し、前記入力色空間のうち前記第一の領域とは異なる第二の領域について前記第一の出力特性とは異なる第二の出力特性を示すように前記対応関係を規定して、前記ルックアップテーブルを生成する、ルックアップテーブル生成方法。
2. 前記特定画像の色分布に基づいて前記入力色空間内の複数の格子点のうち前記特定画像に含まれる色の変換に用いる格子点を特定し、前記第一の領域について前記第一の出力特性を示すように前記特定した格子点について前記対応関係を規定する、請求項１に記載のルックアップテーブル生成方法。
3. 前記特定画像から該特定画像に含まれる色を抽出し、前記入力色空間内の複数の格子点のうち前記抽出した色の変換に用いる格子点を特定し、前記第一の領域について前記第一の出力特性を示すように前記特定した格子点について前記対応関係を規定する、請求項２に記載のルックアップテーブル生成方法。
4. 前記特定画像に含まれる色のそれぞれについて前記入力色空間内の複数の格子点のうち前記色の変換に用いる一以上の格子点を抽出し、該抽出した格子点の重複を解消し、前記第一の領域について前記第一の出力特性を示すように重複解消後の格子点について前記対応関係を規定する、請求項２又は請求項３に記載のルックアップテーブル生成方法。
5. 前記ルックアップテーブルに設けられる第一格子点の数Ｎ１よりも少ない数Ｎ２の第二格子点について前記対応関係を規定し、該第二格子点についての対応関係に基づいて前記数Ｎ１の第一格子点について前記対応関係を規定する、請求項１〜請求項４のいずれか一項に記載のルックアップテーブル生成方法。
6. 前記印刷物を観察する複数の観察光源のそれぞれについて、前記第二格子点に、前記印刷物に形成される前記第一の色材の前記観察光源におけるターゲット色彩値を対応付け、
   前記複数の観察光源について前記第二格子点において対応付けられたターゲット色彩値への近似性を評価する評価値に基づいて、前記色再現画像に形成される前記第二の色材の色彩値が前記観察光源毎に前記ターゲット色彩値に近似するように前記第二の色材の使用量を前記第二格子点について予測し、
   前記第二格子点に、予測された前記第二の色材の使用量を対応付け、前記第二格子点についての前記第一及び第二の色材の使用量の対応関係に基づいて前記数Ｎ１の第一格子点について前記ルックアップテーブルの対応関係を規定する、請求項５に記載のルックアップテーブル生成方法。
7. 前記入力色空間において前記第一格子点を囲む複数の前記第二格子点についての前記第一及び第二の色材の使用量の対応関係を用いた補間演算を行うことにより前記第一格子点について前記ルックアップテーブルの対応関係を規定する、請求項５又は請求項６に記載のルックアップテーブル生成方法。
8. 前記第一の色材は、第三の色材と、第四の色材と、を含み、
   前記第一格子点に対応する前記第三の色材の使用量のそれぞれについて、前記第四の色材で表される色空間内に前記第二格子点を設け、
   前記第三の色材の使用量のそれぞれについて、前記第四の色材で表される色空間において前記第一格子点を囲む複数の前記第二格子点についての前記第一及び第二の色材の使用量の対応関係を用いた補間演算を行うことにより前記第一格子点について前記ルックアップテーブルの対応関係を規定する、請求項５〜請求項７のいずれか一項に記載のルックアップテーブル生成方法。
9. 前記第一の色材は、第三の色材と、第四の色材と、を含み、
   前記第一格子点に対応する前記第三の色材の使用量のそれぞれについて、前記第二格子点に含まれ前記第三の色材の使用量に応じた数の第三格子点を前記入力色空間内に設け、
   該設けた第三格子点のそれぞれについて前記対応関係を規定する、請求項５〜請求項８のいずれか一項に記載のルックアップテーブル生成方法。
10. 前記第一格子点に対応する前記第三の色材の使用量に第一の使用量と該第一の使用量よりも多い第二の使用量が含まれ、
    前記第三の色材の前記第一の使用量に対応する前記第三格子点の数は、前記第三の色材の前記第二の使用量に対応する前記第三格子点の数よりも多い、請求項９に記載のルックアップテーブル生成方法。
11. 前記入力色空間に含まれる第一及び第二の領域について示される前記第一及び第二の出力特性は、それぞれ、所定の観察光源における前記第一の色材の色彩値と前記第二の色材の色彩値との色差、前記色再現画像の粒状性を表す指標値、前記色再現画像の階調性を表す指標値、の少なくとも一つで表される、請求項１〜請求項１０のいずれか一項に記載のルックアップテーブル生成方法。
12. 印刷機で印刷物の形成に使用される第一の色材の使用量と、画像形成装置で色再現画像の形成に使用される第二の色材の使用量と、の対応関係を規定したルックアップテーブルを生成するルックアップテーブル生成装置であって、
    前記第一の色材の色で表される入力色空間のうち特定画像に含まれる色を含む第一の領域について第一の出力特性を示すように前記対応関係を規定し、前記入力色空間のうち前記第一の領域とは異なる第二の領域について前記第一の出力特性とは異なる第二の出力特性を示すように前記対応関係を規定して、前記ルックアップテーブルを生成する、ルックアップテーブル生成装置。
13. 印刷機で印刷物の形成に使用される第一の色材の使用量と、画像形成装置で色再現画像の形成に使用される第二の色材の使用量と、の対応関係を規定したルックアップテーブルであって、前記第一の色材の色で表される入力色空間のうち特定画像に含まれる色を含む第一の領域について第一の出力特性を示すように前記対応関係が規定され、前記入力色空間のうち前記第一の領域とは異なる第二の領域について前記第一の出力特性とは異なる第二の出力特性を示すように前記対応関係が規定されたルックアップテーブルを記憶する記憶部と、
    前記ルックアップテーブルを参照して色変換を行う色変換部と、を備えた色変換装置。
14. 前記入力色空間に含まれる第一及び第二の領域について示される前記第一及び第二の出力特性は、それぞれ、所定の観察光源における前記第一の色材の色彩値と前記第二の色材の色彩値との色差、前記色再現画像の粒状性を表す指標値、前記色再現画像の階調性を表す指標値、の少なくとも一つで表される、請求項１３に記載の色変換装置。

Images