(57)【要約】 【課題】  【解決手段】  カラーインクジェットプリンタ、カラーレーザプリンタ等に用いられる色変換プロファイルの生成方法および装置。一実施形態では、分光プリンティングモデルコンバータを用いて、サンプルインク量データを分光反射率Rsmp(λ)に変換する。そして、この分光反射率Rsmp(λ)を用いて２つの異なる観察条件における測色値ＣＶ１，ＣＶ２を計算し、測色値ＣＶ１，ＣＶ２の間の色差指数を含む評価指数ＥＩを算出する。他の実施形態では、サンプルインク量データで表されるサンプル色と、参照カラーパッチから得られる参照色との間の色差指数を含む評価指数ＥＩが算出される。そして、最良の評価指数ＥＩを有するサンプルインク量データが、代表サンプルとして選択される。こうして選択された複数の代表サンプルからプロファイルが作成される。【選択図】  図１
2007511175.tif 000002
【特許請求の範囲】
【請求項１】
プリンタで使用可能な複数のインクのインク量セットを表すインク量データと、測色値との対応関係を規定するプロファイルを作成する方法であって、
（ａ）インク量データを、前記インク量データに従って印刷されるカラーパッチの分光反射率に変換するように構成された分光プリンティングモデルコンバータを準備する工程と、
（ｂ）複数のインクのインク量セットをそれぞれ表す複数のサンプルインク量データを準備する工程と、
（ｃ）前記分光プリンティングモデルコンバータを用いて、各サンプルインク量データを、前記サンプルインク量データによって表されるインク量で印刷されるべき仮想サンプルパッチの分光反射率に変換する工程と、
（ｄ）各サンプルインク量データについて、前記分光反射率から算出されるサンプル色と、比較の基礎として選択された比較色との間の色差を表す色差指数を含む評価指数を算出する工程と、
（ｅ）前記評価指数に基づいて、複数のサンプルインク量データを選択する工程と、
（ｆ）前記選択された複数のサンプルインク量データに基づいて、測色値とインク量データとの対応関係を規定するプロファイルを作成する工程と、
を備える方法。
【請求項２】
請求項１記載の方法であって、
前記工程（ｄ）は、
第１の観察条件下における前記仮想サンプルカラーパッチの第１の測色値を、前記サンプル色の測色値として算出する工程と、
第２の観察条件下における前記仮想サンプルカラーパッチの第２の測色値を、前記比較色の測色値として算出する工程と、
前記第１と第２の測色値から前記評価指数を算出する工程と、
を含む、方法。
【請求項３】
請求項２記載の方法であって、
前記評価指数を算出する工程は、
前記第１と第２の測色値を、共通する観察条件下における前記サンプル色と前記比較色のそれぞれの対応色の第１と第２の対応測色値に変換する工程と、
前記第１と第２の対応測色値の間の色差を前記評価指数として算出する工程と、
を含む、方法。
【請求項４】
請求項１記載の方法であって、
前記工程（ｄ）は、
第１の観察条件下における前記仮想サンプルカラーパッチの第１の測色値を、前記サンプル色の測色値として算出する工程と、
前記第１の観察条件下における所定の参照カラーパッチの第２の測色値を、前記比較色の測色値として算出する工程と、
前記第１と第２の測色値から前記評価指数を算出する工程と、
を含む、方法。
【請求項５】
請求項１記載の方法であって、
前記工程（ｄ）は、各サンプルインク量データについて、特定の色空間における前記仮想サンプルカラーパッチの測色値を算出する工程を含み、
前記工程（ｅ）は、
前記特定の色空間を複数のセルに分割する工程と、
少なくとも１つのサンプルインク量データに対する少なくとも１つの測色値をそれぞれ含む複数のセルを選択する工程と、
前記評価指数に基づいて、前記特定の色空間の前記選択されたセルのそれぞれに関して１つのサンプルインク量データを選択する工程と、
を含む、方法。
【請求項６】
請求項１記載の方法であって、
前記工程（ｅ）は、
１つのサンプルインク量データの評価指数が、所定の基準を満足するか否かを判断する工程と、
前記サンプルインク量データの評価指数が前記所定の基準を満足しないときには、
（ｉ）前記サンプルインク量データを修正する工程と、
（ｉｉ）前記修正されたサンプルインク量データの評価指数を算出する工程と、
を、前記修正されたサンプルインク量データの評価指数が前記基準を満足するまで繰り返す工程と、
前記所定の基準を満たす評価指数を有するサンプルインク量データを選択する工程と、
を含む、方法。
【請求項７】
プリンタで使用可能な複数のインクのインク量セットを表すインク量データと、測色値との対応関係を規定するプロファイルを作成する装置であって、
インク量データを、前記インク量データに従って印刷されるカラーパッチの分光反射率に変換する分光プリンティングモデルコンバータであって、複数のサンプルインク量データのそれぞれを、前記サンプルインク量データによって表されるインク量で印刷されるべき仮想サンプルパッチの分光反射率に変換する分光プリンティングモデルコンバータと、
各サンプルインク量データについて、前記分光反射率から算出されるサンプル色と、比較の基礎として選択された比較色との間の色差を表す色差指数を含む評価指数を算出する評価指数生成部と、
前記評価指数に基づいて、複数のサンプルインク量データを選択する選択部と、
前記選択された複数のサンプルインク量データに基づいて、測色値とインク量データとの対応関係を規定するプロファイルを作成するプロファイル作成部と、
を備える装置。
【請求項８】
印刷方法であって、
請求項１ないし７のいずれかに従って作成されたプロファイルに基づいて、カラー画像データをインク量データに変換するためのルックアップテーブルを作成する工程と、
前記プロファイルを用いて、与えられたカラー画像データをインク量データに変換する工程と、
前記インク量データをハーフトーン処理して印刷データを作成する工程と、
前記印刷データに従って画像を印刷する工程と、
を備える印刷方法。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、カラーインクジェットプリンタ、カラーレーザプリンタ等に用いられる色変換プロファイルの生成方法および装置に関し、また、その色変換プロファイルを用いるプリンタに関する。
【背景技術】
【０００２】
近年、カラープリンタとしてインクジェットプリンタやレーザプリンタなどが普及してきている。カラープリンタでは、入力カラー画像データを、対応する複数色のインクの量に変換する色変換ルックアップテーブルが用いられる。従来の色変換ルックアップテーブルの作成方法では、(i) 多数のサンプルカラーパッチを印刷し、(ii) 各カラーパッチを測色して測色値を求め、(iii) 各カラーパッチの測色値と各カラーパッチの印刷に使用されたインク量との対応関係を表すルックアップテーブルを作成していた。
【発明の開示】
【発明が解決しようとする課題】
【０００３】
しかし、測色値は、測色時の観察条件（viewing condition）に依存する。従って、ある観察条件で好ましい色再現を有する印刷物が得られるように色変換ルックアップテーブルを作成したとしても、その色変換ルックアップテーブルを用いて作成された印刷物が他の観察条件では好ましい色再現を達成できない場合がある。そこで、従来から、種々の観察条件で好ましい色再現を達成できる色変換ルックアップテーブルを作成することが望まれていた。また、例えば絵画の複製を作るために、原画の色の見えをなるべく忠実に再現することのできる色変換ルックアップテーブルを作成したいという要望もあった。これらの要望は、色変換ルックアップテーブルに限らず、一般に、色変換に使用される各種のプロファイルに共通する要望であった。
【課題を解決するための手段】
【０００４】
本発明は、種々の観察条件で好ましい色再現を達成できる色変換プロファイルを作成する方法及び装置を提供する。また、本発明は、原画の色の見えなるべく忠実に再現することのできる色変換プロファイルを作成する方法及び装置を提供する。
【０００５】
より詳細に言えば、本発明は、プリンタで使用可能な複数のインクのインク量セットを表すインク量データと、測色値との対応関係を規定するプロファイルを作成する方法を提供する。この方法は、（ａ）インク量データを、前記インク量データに従って印刷されるカラーパッチの分光反射率に変換するように構成された分光プリンティングモデルコンバータを準備する工程と、（ｂ）複数のインクのインク量セットをそれぞれ表す複数のサンプルインク量データを準備する工程と、（ｃ）前記分光プリンティングモデルコンバータを用いて、各サンプルインク量データを、前記サンプルインク量データによって表されるインク量で印刷されるべき仮想サンプルパッチの分光反射率に変換する工程と、（ｄ）各サンプルインク量データについて、前記分光反射率から算出されるサンプル色と、比較の基礎として選択された比較色との間の色差を表す色差指数を含む評価指数を算出する工程と、（ｅ）前記評価指数に基づいて、複数のサンプルインク量データを選択する工程と、（ｆ）前記選択された複数のサンプルインク量データに基づいて、測色値とインク量データとの対応関係を規定するプロファイルを作成する工程と、を備える。複数のインクは、１つ以上の特色インクを含むことが好ましい。
【０００６】
この方法によれば、適切な評価指数を設定することによって、種々の観察条件で好ましい色再現を達成できる色変換プロファイルを生成することができ、また、原画の色の見えをなるべく忠実に再現することのできる色変換プロファイルを作成することができる。
【０００７】
なお、評価指数としては、非色恒常性指数(Color Inconstancy Index)を含むものや、メタメリズム指数(Metamerism Index)を含むものなどを利用することができる。非色恒常性指数を含む評価指数を用いたときには、種々の観察条件で好ましい色再現を達成できる色変換プロファイルを作成することができる。一方、メタメリズム指数を用いたときには、原画の色の見えをなるべく忠実に再現することのできる色変換プロファイルを作成することができる。メタメリズム指数は、２つの色の色差を表す色差評価指数ＣＤＩである。非色恒常性指数は、１つのサンプルを２つの観察条件で観察したときの色差を表す色差評価指数ＣＤＩである。
【０００８】
なお、本発明は、種々の形態で実現することが可能であり、例えば、プロファイル作成方法およびプロファイル作成装置、プロファイルを用いた印刷方法および印刷装置、そのプロファイルを用いた印刷方法および印刷装置、それらの方法または装置の機能を実現するためのコンピュータプログラム、そのコンピュータプログラムを記録した記録媒体、そのコンピュータプログラムを含み搬送波内に具現化されたデータ信号、等の形態で実現することができる。
【発明を実施するための最良の形態】
【０００９】
次に、本発明の実施の形態を以下の順序で説明する。
Ａ．第１実施形態
Ｂ．第２実施形態
Ｃ．第３実施形態
Ｄ．分光プリンティングモデルの一例
Ｅ．変形例
【００１０】
Ａ．第１実施形態：
図１は、本発明の好ましい実施形態におけるプロファイル作成システムの構成を示すブロック図である。このシステムは、分光プリンティングモデルコンバータ１００と、評価指数生成部１２０と、選択部１３０と、プロファイル生成部１４０と、ガマットマッピング処理部１６０とを備えているが、このシステムはこれ以外の数や種類の要素を異なる構成で備えていてもよい。分光プリンティングモデルコンバータ１００は、インク量データを、そのインク量データに応じて印刷されるカラーパッチの分光反射率Rsmp(λ)に変換する。なお、本明細書において「カラーパッチ」とは、有彩色のパッチに限らず、無彩色のパッチも含む広い意味を有している。この実施形態では、シアン（Ｃ）、マゼンタ（Ｍ）、イエロー（Ｙ）、ブラック（Ｋ）、オレンジ（Ｏ）、および、グリーン（Ｇ）の６色のインクを利用可能なカラープリンタを想定しており、分光プリンティングモデルコンバータ１００もこの６種類のインクの吐出量を入力としている。分光プリンティングモデルについては、Ｄ節で詳述する。以下では、「分光プリンティングモデル」を「フォワードモデル」とも呼ぶ。
【００１１】
評価指数生成部１２０は、サンプル色算出部１２２と、比較色算出部１２４と、評価指数算出部１２６とを含んでいるが、評価指数生成部１２０はこれ以外の数や種類の要素を異なる組合せで備えていてもよい。サンプル色算出部１２２は、サンプルインク量データに対する分光反射率Rsmp(λ)を用いて、第１の観察条件における測色値ＣＶ１を算出する。本実施形態では、第１の観察条件として標準の光Ｄ５０を用いる。この第１の観察条件で得られた測色値ＣＶ１で表される色を「サンプル色」とも呼ぶ。比較色算出部１２４は、サンプルインク量データに対する分光反射率Rsmp(λ)を用いて、第２の観察条件における測色値ＣＶ２を算出する。本実施形態では、第２の観察条件として標準の光Ｆ１１を用いる。以下では、第２の観察条件で得られた測色値ＣＶ２で表される色を「比較色」とも呼ぶ。
【００１２】
上述の説明から理解できるように、第１実施形態において、サンプル色算出部１２２と比較色算出部１２４は、同じ分光反射率Rsmp(λ)を用いて異なる観察条件における測色値ＣＶ１，ＣＶ２をそれぞれ算出する。評価指数算出部１２６は、これらの測色値ＣＶ１，ＣＶ２を用いて、そのサンプルインク量データの良否を判定するための評価指数ＥＩ1を算出する。評価指数ＥＩ1の具体例については後述する。
【００１３】
選択部１３０は、良好な評価指数ＥＩ1を有するサンプルインク量データを選択する。プロファイル生成部１４０は、選択されたサンプルインク量データと、それらのサンプルインク量データを用いて印刷されるカラーパッチの測色値（Ｌ*ａ*ｂ*値）とを用いて、インクプロファイル１４２を作成する。このインクプロファイル１４２は、測色値（Ｌ*ａ*ｂ*値）とＣＭＹＫＯＧのインク量との対応関係を示すルックアップテーブルである。なお、「インクプロファイル」を「出力デバイスプロファイル」とも呼ぶ。本明細書において、「プロファイル」とは、色空間の変換を行うための変換規則を具現化したものを意味しており、各種のデバイスプロファイルとルックアップテーブルとを含む広い意味を有している。
【００１４】
ガマットマッピング処理部１６０は、このインクプロファイル１４２と、予め準備された入力デバイスプロファイル１６２とを用いて、プリンタルックアップテーブル１８０を作成する。ここで、入力デバイスプロファイル１６２としては、例えばｓＲＧＢ色空間をＬ*ａ*ｂ*色空間に変換するプロファイルを使用することができる。プリンタルックアップテーブル１８０は、入力カラー画像データ（例えばｓＲＧＢデータ）をインク量データに変換するためのものである。
【００１５】
図２は、第１実施形態の処理手順を示すフローチャートである。ステップＳ１では、分光プリンティングモデルを決定して、コンバータ１００を作成する。一実施例では、分光プリンティングモデルとしてセル分割ユール・ニールセン分光ノイゲバウアモデル(Cellular Yule-Nielsen Spectral Neugebauer Model)を使用する。その詳細についはＤ節で説明する。
【００１６】
ステップＳ２では、多数の仮想サンプルを設定する。ここで、「仮想サンプル」とは、プロファイルの作成処理において使用される仮のインク量データ、および、そのインク量データに応じて印刷される仮想的なカラーパッチを意味している。なお、以下では、仮想サンプルを単に「サンプル」とも呼ぶ。一実施例では、ＣＭＹＫＯＧの各インクに関して、０〜１００％の範囲で１０％おきに１１点のインク量をそれぞれ設定し、６種類のインクの量のすべての組合せを仮想サンプル（サンプルインク量データ）として準備した。この結果、１１6＝１，７７１，５１６個の仮想サンプルが準備された。なお、「インク量１００％」とは、１種類のインクでベタ打ちとなるインク量を意味する。
【００１７】
ステップＳ３では、分光プリンティングモデルコンバータ１００を用いて各仮想サンプルのサンプルインク量データを分光反射率Rsmp(λ)に変換し、この分光反射率Rsmp(λ)からＣＩＥＬＡＢ表色系の測色値Ｌ*ａ*ｂ*を算出する。一実施例では、ＣＩＥ標準の光Ｄ５０、および、ＣＩＥ１９３１  ２°観測者の観察条件で測色値を算出した。以下では、仮想サンプルを特定の観察条件で観察したときの色を「サンプル色」と呼ぶ。図３（Ａ）〜３（Ｃ）は、一実施例で算出されたサンプル色の分布を示している。図３（Ａ）の横軸はＣＩＥＬＡＢ表色系のａ*軸であり、縦軸はｂ*軸である。図３（Ｂ），（Ｃ）の横軸はａ*軸およびｂ*軸であり、縦軸はＬ*軸である。これから理解できるように、１１6個のサンプル色は、明度Ｌ*の低いところにより集中しており、明度Ｌ*の高いところではまばらに分布している。サンプル色をより均一に分布させるためには、サンプルインク量を、比較的低いインク量の範囲ではより細かい間隔で設定し、比較的高いインク量の範囲ではより粗い間隔で設定するようにしてもよい。
【００１８】
ステップＳ４では、測色値の色空間（ここではＣＩＥＬＡＢ空間）を複数のセルに分割し、複数のサンプル色をセルに関してソート（分類）する。一実施例では、ＣＩＥＬＡＢ空間を１６×１６×１６個のセルに均等に分割した。
【００１９】
ステップＳ５では、好ましいサンプルの選択に使用する評価指数ＥＩ1を設定する。第１実施形態で使用される評価指数ＥＩ1は、以下の（１）式で表される。
【数１】
2007511175.tif 000003
ここで、ＣＩＩは非色恒常性指数(Color Inconstancy Index)であり、ｋは係数、αはその他の指数である。
【００２０】
ＣＩＩは、例えば以下の式で表現される。
【数２】
2007511175.tif 000004
ここで、ΔＬ*はＣＶ１とＣＶ２の条件下におけるサンプルの明度差、ΔＣ*abは彩度差、ΔＨ*abは色相差を示す。ＣＩＩの計算時には、ＣＶ１とＣＶ２のＣＩＥＬＡＢ値は、図４に即して後で説明するように、色順応変換（ＣＡＴ）を用いてＤ６５のような標準観察条件に変換される。ＣＩＩについては、Billmeyer and Saltzman's Principles of Color Technology, 3rd edition, John Wiley & Sons, Inc, 2000, p.129, pp. 213-215を参照。
【００２１】
なお、（２）式の右辺は、ＣＩＥ１９９４年色差式において、明度と彩度の係数ｋL，ｋCの値を２に設定し、色相の係数ｋHの値を１に設定した色差ΔＥ*94(2:2)に相当する。ＣＩＥ１９９４年色差式では、（２）式の右辺の分母の係数ＳL，Ｓc，ＳHは以下の（３）式で与えられる。
【数３】
2007511175.tif 000005
【００２２】
なお、ＣＩＩの算出に使用する色差式としては、他の式を用いることも可能である。
【００２３】
ＣＩＩは、あるカラーパッチを第１と第２の異なる観察条件下で観察したときの色の見えの差として定義されている。従って、ＣＩＩが小さいサンプルは、異なる観察条件での色の見えの差が小さいという点で好ましい。
【００２４】
一実施例においては、以下の（４）式で与えられる評価指数ＥＩ1が使用された。
【数４】
2007511175.tif 000006
ここで、max(CII)は、全サンプル中のＣＩＩの最大値である。また、Ｔinkは、そのサンプルで使用されるインク量の合計値である。例えば、６種類のインクのインク量がすべて２０％のときには、Ｔink／６の値は１２０％／６＝０．２になる。合計インク量Ｔinkは、合計インク量を最小化してコストを節約するために使用されている。従って、（４）式は、ＣＩＩがほぼ同じ場合には、合計インク量Ｔinkが少ないほど良好な（すなわち小さな）評価指数ＥＩ1が得られる式となっている。なお、第１実施形態の評価指数ＥＩ1としては、（１）式や（４）式に示したものに限らず、ＣＩＩを含む任意の関数を利用することが可能である。
【００２５】
図２のステップＳ６では、評価指数生成部１２０が、各サンプルに対する評価指数ＥＩ1を算出し、選択部１３０がこの評価指数ＥＩ1に応じてＣＩＥＬＡＢ色空間の各セル内で最良のサンプルを選択する。
【００２６】
図４は、ステップＳ６の詳細手順を示すフローチャートである。ステップＳ１０では、あるセル内の１つのサンプルを選択する。ステップＳ１１では、分光プリンティングモデルコンバータ１００を用いて、そのサンプルの分光反射率Rsmp(λ)を算出する。次のステップＳ１２〜Ｓ１４は評価指数生成部１２０のサンプル色算出部１２２（図１）によって実行され、ステップＳ１５〜Ｓ１７は比較色算出部１２４によって実行される。
【００２７】
サンプル色算出部１２２は、ステップＳ１２において、分光反射率Rsmp(λ)を用いて第１の観察条件下で三刺激値ＸＹＺを算出する。一実施例では、ＣＩＥ標準の光Ｄ５０、および、ＣＩＥ１９３１  ２°観測者の観察条件で三刺激値ＸＹＺを算出した。本明細書において、「観察条件」とは照明光と観測者の組合せを意味しているが、特に言及しない限り観測者としてＣＩＥ１９３１  ２°観測者を使用する。ステップＳ１３では、この三刺激値ＸＹＺに色順応変換を適用して、標準観察条件での対応色を算出する。一実施例では、標準観察条件の光源として標準の光Ｄ６５を用い、色順応変換としてCIECAT02を利用した。CIECAT02については、例えば"The CIECAM02 Color Appearance Model", Nathan Moroney et al., IS&T/SID Tenth Color Imaging Conference, pp.23-27, および、"The performance of CIECAM02", Changjun Li et al., IS&T/SID Tenth Color Imaging Conference, pp.28-31に記載されている。但し、色順応変換としては、フォン・クリースの色順応予測式などの他の任意の色順応変換を用いることも可能である。ステップＳ１４では、この対応色のＣＩＥＬＡＢ表色系の測色値CV1=（L*a*b*)D50→D65を算出する。この下付き文字「Ｄ５０→Ｄ６５」は、標準の光Ｄ５０の下での色の見えを、標準の光Ｄ６５の対応色で表現した測色値であることを意味している。
【００２８】
比較色算出部１２４も、サンプル色算出部１２２と同様の演算を実行するが、第２の観察条件の下で行う。すなわち、ステップＳ１５において、分光反射率Rsmp(λ)を用いて第２の観察条件下で三刺激値ＸＹＺを算出する。一実施例では、ＣＩＥ標準の光Ｆ１１、および、ＣＩＥ１９３１  ２°観測者の観察条件で三刺激値ＸＹＺを算出した。ステップＳ１６では、この三刺激値ＸＹＺに色順応変換を適用して、標準観察条件での対応色を算出する。そして、ステップＳ１７では、この対応色のＣＩＥＬＡＢ表色系の測色値CV2=（L*a*b*)F11→D65を算出する。
【００２９】
サンプル色の測色値CV1=（L*a*b*)D50→D65と、比較色の測色値CV2=（L*a*b*)F11→D65は、同一の標準観察条件におけるそれぞれの対応色の測色値なので、それらの色差ΔＥであるＣＩＩ（上記（２）式参照）は、サンプル色と比較色の色の見えの違いをかなり正確に表現する値となる。
【００３０】
なお、標準観察条件は、標準の光Ｄ６５に限らず、任意の照明光下での観察条件を採用することができる。例えば、標準の光Ｄ５０を標準観察条件として採用した場合には、図４のステップＳ１３は不要であり、また、ステップＳ１６では標準の光Ｄ５０に対する色順応変換が実行される。但し、ＣＩＥＬＡＢ表色系を用いて算出される色差ΔＥは、標準の光Ｄ６５を用いたときに最も信頼性の高い値を示す。この点からは、標準観察条件として標準の光Ｄ６５を用いることが好ましい。
【００３１】
ステップＳ１８では、評価指数算出部１２６（図１）が、これらの２つの測色値ＣＶ１，ＣＶ２を用い、上記（４）式に従って評価指数ＥＩ1を算出する。
【００３２】
ステップＳ１９では、そのセルに含まれるすべてのサンプル色に関して評価指数ＥＩ1の算出が終了したか否かが判断される。こうして、ステップＳ１０〜Ｓ１９が繰り返し実行されて、そのセル内のすべてのサンプル色に関して評価指数ＥＩ1が算出される。ステップＳ２０では、選択部１３０が、そのセル内のサンプル色のうちで、評価指数ＥＩ1が最良であるサンプルを、そのセルに関する代表サンプルとして選択する。この結果、少なくとも１つのサンプルを含む各セルに関して、１つの代表サンプルがそれぞれ選択される。以下では、代表サンプルを「高評価サンプル」とも呼ぶ。
【００３３】
なお、ステップＳ４で分割された複数のセルの中には、サンプル色を全く含まないセルも存在する。従って、図４の処理は、少なくとも１つのサンプル色を含むようなセルを対象として実行され、サンプル色を１つも含まないセルは処理対象から除外される。
【００３４】
図２のステップＳ７では、プロファイル生成部１４０が、代表サンプルに基づいて不等間隔補間を行うことによって、予備インクプロファイルを作成する。この予備インクプロファイルは、ＣＩＥＬＡＢ測色値をインク量に変換する色変換ルックアップテーブルである。「予備的」という接頭語は、ステップＳ４で分割された比較的粗いセルに関するプロファイルであることを意味している。
【００３５】
図５（Ａ）は、ステップＳ７における不等間隔補間を示す。この図は、ＣＩＥＬＡＢ空間を示しており、図中の丸は代表サンプル色の測色値の位置を示し、網模様は細かなセルのグリッドを示している。ステップＳ７では、グリッド点（網模様の交点）におけるインク量が、複数の代表サンプル色のインク量を不等間隔補間することによって算出される。図５（Ｂ）および図５（Ｃ）は、ａ*ｂ*色平面上のL*=23．8での補間前と補間後の代表サンプル点の例をそれぞれ示している。この不等間隔補間は、例えばMATLAB（MathWorks, Inc.の商標）のgriddata関数を用いて実行することができる。一実施例では、予備プロファイルとして、ＣＩＥＬＡＢ空間の６４×６４×６４個のグリッド点を入力とするプロファイルを作成した。なお、不等間隔補間としては、非線形補間を用いても良く、あるいは線形補間を用いても良い。通常は、非線形補間の方が線形補間よりも精度が高く、処理速度は遅い傾向にある。
【００３６】
ステップＳ８では、プロファイル生成部１４０が、予備プロファイルを線形補間することによって、最終インクプロファイル１４２（図１）を作成する。この最終インクプロファイル１４２は、予備プロファイルよりもさらに細かなセルのグリッド点を入力とするプロファイルである。一実施例では、ＣＩＥＬＡＢ空間の２５６×２５６×２５６個のグリッド点を入力とするプロファイルを作成した。上述したように、予備プロファイルはＣＩＥＬＡＢ空間の６４×６４×６４個のグリッド点を入力としているので、予備プロファイルを線形補間することによって最終インクプロファイル１４２を容易に得ることができる。最終インクプロファイル１４２として、２５６×２５６×２５６個のグリッド点を入力とするプロファイルを作成すれば、ＣＩＥＬＡＢのすべての入力値に対するインク量を直ちに得ることができる。従って、その後のルックアップテーブルの作成に要する処理時間を短縮することができる。但し、最終インクプロファイル１４２の代わりに予備インクプロファイルを用いて、後述するプリンタルックアップテーブル１８０を作成することも可能である。
【００３７】
ステップＳ９では、ガマットマッピング処理部１６０（図１）が、インクプロファイル１４２とｓＲＧＢプロファイル１６２とに基づいてガマットマッピングを行い、プリンタルックアップテーブル１８０を作成する。ガマットマッピングを行う理由は、プリンタで実現可能な色空間（「インク色空間」とも呼ぶ）の色域と、入力色空間（この実施形態ではｓＲＧＢ空間）で実現可能な色域とに差があるためである。インク色空間の色域はインクプロファイル１４２で規定されており、入力色空間の色域はｓＲＧＢプロファイル１６２で規定されている。一般に、入力色空間とインク色空間には食い違いがあるので、入力色空間の色域をインク色空間の色域にマッピングする必要がある。
【００３８】
図６（Ａ）および６（Ｂ）は、ガマットマッピングの一例を示している。ここでは、いわゆるガマットクリッピングと呼ばれる方法が採用されている。具体的には、図６（Ａ）に示すように、インク色空間の色域の外側にあるｓＲＧＢ色空間の色は、色相を保った状態で彩度を低下させるようにマッピングされる。明度Ｌ*に関しては、インク色空間の明度範囲内にある色は、明度がそのまま維持される。インク色空間の明度の最大値Ｌmaxよりも大きな明度を有する色は、最大値Ｌmaxにマッピングされる。一方、インク色空間の明度の最小値Ｌminよりも大きな明度を有する色は、最小値Ｌminにマッピングされる。なお、ガマットマッピングの方法としては、従来から種々の方法が知られており、そのいずれの方法を採用してもよい。
【００３９】
こうしてガマットマッピングが行われると、最終的なプリンタルックアップテーブル１８０が完成する。このルックアップテーブル１８０は、ｓＲＧＢデータを入力とし、６種類のインクのインク量を出力とするルックアップテーブルである。このプリンタルックアップテーブル１８０をプリンタ又はコンピュータ内に存在する関連するプリンタドライバに実装すれば、非色恒常性の小さい（すなわち、異なる観察条件における色の見えの変化が小さな）印刷物を作成することが可能である。実際の印刷の際には、プリンタルックアップテーブル１８０が入力カラー画像データ（ここではｓＲＧＢデータ）をインク量データに変換し、ハーフトーン処理部がこのインク量データをハーフトーン処理することによってインクドットの記録状態を示す印刷データを生成し、印刷ヘッドが印刷データに応じた量のインク滴を吐出することによってカラー印刷物を作成する。
【００４０】
図７は、一実施例と比較例における色差の算出結果を示すグラフである。実線は、オリジナルカラーパッチと、オリジナルカラーパッチを一実施例のＣＭＹＫＧ０プリンタで印刷したカラーパッチとの色差ΔＥを示している。破線は、オリジナルパッチと、オリジナルカラーパッチを比較例としてのＣＭＹＫｃｍプリンタで印刷したカラーパッチとの色差ΔＥを示している。ＣＭＹＫｃｍプリンタは、ＣＭＹＫの４種類のインクの他に、ライトシアンインクｃとライトマゼンタインクｍを使用するセイコーエプソン社製の６色プリンタである。図７の横軸は、オリジナルパッチの番号である。ここでは、４４個のオリジナルカラーパッチを使用した。最初の２０個のカラーパッチはコダック社のグレースケールであり、残りの２４個のカラーパッチはMacbeth社の GretagMacbeth Color Checker Color Rendition Chart（Macbeth社の商標）である。プリンタでオリジナルカラーパッチを印刷する際には、まず、オリジナルカラーパッチの分光反射率を測定し、その測色値（ＣＩＥＬＡＢ値）を算出した。そして、その測色値をプリンタドライバに入力してパッチを印刷した。なお、この実験の際には、プリンタルックアップテーブル１８０（図１）の代わりにインクプロファイル１４２を用いて、測色値（ＣＩＥＬＡＢ値）からインク量を求めた。
【００４１】
図７の結果では、ほぼすべてのカラーパッチに関して、比較例のＣＭＹＫｃｍプリンタよりも実施例のＣＭＹＫＧＯプリンタの方がオリジナルと再現パッチとの色差が小さい。ＣＭＹＫｃｍプリンタで色差が大きいのは、ＣＭＹＫｃｍプリンタにおいて通常のカラープロファイルが使用されているからである。
【００４２】
図８は、一実施例と比較例における非色恒常性指数（Color Inconstancy Index）の算出結果を示すグラフである。実線は、一実施例のＣＭＹＫＧ０プリンタで印刷したカラーパッチのＣＩＩを示し、破線は、比較例のＣＭＹＫｃｍプリンタで印刷したカラーパッチのＣＩＩを示している。また、一点鎖線はオリジナルカラーパッチのＣＩＩを示している。なお、図８におけるオリジナルカラーパッチと、実施例および比較例のカラーパッチは、図７のものと同じである。
【００４３】
図８の結果では、ほぼすべてのカラーパッチに関して、比較例のＣＭＹＫｃｍプリンタよりも実施例のＣＭＹＫＧＯプリンタの方がＣＩＩが小さい。従って、実施例のＣＭＹＫＧＯプリンタの方が、比較例のＣＭＹＫｃｍプリンタよりも、異なる観察条件における色の見えの変化が小さな印刷物を作成することができる。また、図８の例では、実施例のＣＭＹＫＧＯプリンタで作成されたカラーパッチの方が、オリジナルカラーパッチよりもＣＩＩが小さい場合がある。この点から見ても、実施例のＣＭＹＫＧＯプリンタが非色恒常性を小さくするという目的を十分良く達成できていることが理解できる。
【００４４】
以上のように、第１実施形態では、ＣＩＩ（非色恒常性指数）を含む評価指数ＥＩ1を用いて、サンプルインク量データの中から高い評価値を有するものを選択し、選択されたサンプルインク量データに基づいてインクプロファイルを作成する。この結果、非色恒常性が小さく、種々の異なる観察条件下における色の見えの変化が小さな印刷物を作成可能な色変換プロファイルを作成することができる。
【００４５】
Ｂ．第２実施形態：
図９は、本発明の第２実施形態としてのシステム構成を示すブロック図である。図１に示した第１実施形態のシステムとの大きな違いは、参照カラーパッチ１０２を用いる点と、評価指数ＥＩとして非色恒常性指数（ＣＩＩ）の代わりにメタメリズム指数（後述する）を含む指数を使用する点である。
【００４６】
図１０は、第２実施形態の処理手順を示すフローチャートである。ステップＳ２１では、分光プリンティングモデルを決定してコンバータ１００を作成する。このステップは、図２のステップＳ１と同じである。ステップＳ２２では、複数の比較色を準備する。図９の例では、複数の比較色を与えるための原画として参照カラーパッチ１０２を準備し、これらの参照カラーパッチ１０２の分光反射率Rref(λ)を図示しない分光反射率計で測定する。なお、参照カラーパッチ１０２は、複数のカラーパッチを含んでいるオリジナルカラーパッチである。参照カラーパッチ１０２の代わりに、例えば絵画から複数の比較色を取得するようにしてもよい。なお、参照カラーパッチ１０２や絵画から得られる色を「参照色」とも呼ぶ。ステップＳ２３では、比較色算出部１２４ａが、参照色（比較色）の分光反射率Rref(λ)を用いて、複数の観察条件における測色値ＣＶ２をそれぞれ算出する。本実施形態では、複数の観察条件として３種類の標準の光Ｄ５０，Ｆ１１，Ａを用いる。
【００４７】
ステップＳ２４では、各比較色に対して、その比較色と条件等色度が高いと期待される複数の仮想サンプル（サンプルインク量データ）を準備する。ここで条件等色度とは、一対の物体色（サンプル色と比較色）を複数の照明光下で観察したときの条件等色の安定性（いわゆる照明光条件等色度）を意味している。ステップＳ２５では、これらの仮想サンプルに対する測色値ＣＶ１を算出する。具体的には、第１実施形態と同様に、各仮想サンプルのインク量データを分光プリンティングモデルコンバータ１００で分光反射率Rsmp(λ)に変換し、この分光反射率Rsmp(λ)を用いて、比較色に対して用いたのと同じ複数の観察条件における測色値ＣＶ１をそれぞれ算出する。
【００４８】
ステップＳ２６では、次の（５）式で表される評価指数ＥＩ2を設定する。
【数５】
2007511175.tif 000007
ここで、演算子aveは括弧内の平均を取る演算を意味しており、ＭＩiはｉ番目の観察条件下（すなわちｉ番目の照明光下）におけるメタメリズム指数、ｋは係数、αはその他の指数である。
【００４９】
ＭＩiは、例えば以下の式で表現される。
【数６】
2007511175.tif 000008
【００５０】
このメタメリズム指数ＭＩiは、ｉ番目の照明光下におけるサンプル色と比較色の色差である。この色差は、サンプル色と比較色とが標準観察条件でマッチするようにパラメリック補正がなされた後に算出される。上記（５）式の右辺第１項は、照明光条件等色度を表す指数であることが理解できる。メタメリズム指数については、Billmeyer and Saltzman's Principles of Color Technology, 3rd edition, John Wiley & Sons, Inc, 2000, p.127, pp. 213-215を参照。
【００５１】
上述した（２）式と比較すれば理解できるように、ＭＩを与える式としては、ＣＩＩを与える式と同じ色差式を用いることができる。ＭＩとＣＩＩとの違いは、前者は２つの物体色を同一の観察条件下で観察したときの色の差であり、後者は１つの物体色を異なる観察条件下で観察したときの色の差である点である。なお、ＭＩの算出に使用する色差式として（６）式以外の式を用いることも可能である。（５）式の右辺の第２項αとしては任意の指数を用いることができ、また、αを省略してもよい。
【００５２】
ステップＳ２７では、評価指数生成部１２０ａが、各サンプルに対する評価指数ＥＩ2を算出し、選択部１３０がこの評価指数ＥＩ2に応じて各参照色に対する最良のサンプルを選択する。ステップＳ２８では、各参照色に対して選択されたサンプルを用いて、インクプロファイル１４２およびプリンタルックアップテーブル１８０が生成される。このステップＳ２８は、図２のステップＳ７〜Ｓ９とほぼ同じ処理である。なお、絵画を忠実に再現するためのルックアップテーブとして、例えば絵画の色の色番号を入力とし、インク量を出力とするテーブルを作成してもよい。この場合には、ステップＳ７，Ｓ８の補間演算やステップＳ９のガマットマッピングは省略可能である。
【００５３】
このように、第２実施形態では、メタメリズム指数ＭＩを含む評価指数ＥＩ2を用いて、サンプルインク量データの中から高い評価値を有するものを選択し、選択されたサンプルインク量データに基づいてインクプロファイルを作成する。この結果、プリンタで使用された場合に、種々の異なる観察条件下において参照色との色差が小さな色を有する印刷物を作成可能な色変換プロファイルを作成することができる。この第２実施形態は、特に、原画を忠実に再現するための色変換プロファイルの作成に有用である。
【００５４】
こうして作成されたプロファイルを用いて再現された印刷物の品質は、分光反射率の差異を表す指数、例えばＲＭＳによって評価できる。ＲＭＳ（Root Mean Square）は、２つの分光反射率の差の根２乗平均である。なお、２つの分光反射率の差は、所定の波長幅毎に（例えば５ｎｍまたは１０ｎｍ毎に）計算される。
【００５５】
Ｃ．第３実施形態：
図１１は、本発明の第３実施形態としてのシステム構成を示すブロック図である。図１に示した第１実施形態のシステムとの大きな違いは、選択部１３０ａが基準判定部２００を含む点と、サンプルデータ修正部２１０が追加されている点である。この第３実施形態では、あるサンプルが所定の評価基準を満たさない場合には、サンプルデータ修正部２１０がそのサンプルインク量データを修正し、修正されたサンプルインク量データに対する評価指数ＥＩの算出が再度実行される。そして、評価基準を満足したサンプルを利用して、色変換プロファイルが作成される。
【００５６】
図１２は、第３実施形態の処理手順を示すフローチャートである。ステップＳ３１では、分光プリンティングモデルを作成する。このステップＳ３１は、図２のステップＳ１と同じものである。ステップＳ３２では、ＣＩＥＬＡＢ空間を複数のセルに分割する。ここでは、第１実施形態のステップＳ４で使用されたものと同じ１６×１６×１６個のセル分割を利用できる。
【００５７】
ステップＳ３３では、サンプルの良否を決定するための評価指数ＥＩを設定する。第３実施形態では、第１実施形態の評価指数ＥＩ1＝ｆ（ＣＩＩ）を利用する。この代わりに、第２実施形態の評価指数ＥＩ2＝ｆ（ＭＩ）を用いることも可能である。
【００５８】
ステップＳ３４〜Ｓ３７は、各セル毎に１つの代表サンプルを選択するための再帰的ルーチンである。ステップＳ３４では、ＣＩＥＬＡＢ空間内の１つのセルを処理対象（対象セル）として選択し、その対象セルにおける初期サンプルインク量データを設定する。この初期サンプルインク量データは、そのインク量データに応じて印刷されるサンプル色の測色値（Ｌ*ａ*ｂ*値）が、対象セル内に存在するものである。サンプル色の測色値は、第１の観察条件（例えば標準の光Ｄ５０およびＣＩＥ１９３１  ２°観測者）で算出される。設定された初期サンプルインク量データの測色値が対象セル内に入っていない場合には、対象セル内の測色値が得られるまで初期サンプルインク量データが修正される。
【００５９】
なお、セルによっては、セル内の測色値を与えるようなインク量データが存在しない場合がある。例えば、高明度または低明度で彩度が高いセルの色は再現できない。このような場合には、そのセルは、処理の対象外としてその後の処理から除外される。
【００６０】
ステップＳ３５では、評価指数算出部１２０が、初期サンプルインク量データの評価指数ＥＩを算出する。この評価指数ＥＩは、例えば上記（４）式に従って算出してもよい。ステップＳ３６では、基準判断部２００が、評価指数ＥＩが予め設定された評価基準を満足するか否かを判断する。評価基準としては、例えば、以下の（７）式を用いることができる。
【数７】
2007511175.tif 000009
ここで、δは評価指数ＥＩの許容上限値である。
【００６１】
上記（７）式を用いる場合には、評価指数ＥＩが許容上限値δ以下の場合に基準を満足するものと判定される。なお、１つの評価指数ＥＩを使用する代わりに、１つのサンプルインク量データに関して複数種類の異なる評価指数を算出し、各評価指数がそれぞれの評価基準をすべて満足する場合に、そのサンプルインク量データが基準を満足すると判断するようにしてもよい。
【００６２】
初期サンプルインク量データが評価基準を満足しない場合には、ステップＳ３７においてサンプルデータ修正部２１０が初期サンプルインク量データを修正する。このとき、修正後のサンプルインク量データに関しては、以下のようないくつかの制限条件を付けることが好ましい。
（制限条件１）：修正後のサンプルインク量データで与えられる測色値が、対象セル内にある。
（制限条件２）：修正後のサンプルインク量データで表されるインク量が、インクデューティ制限値を満足する。
【００６３】
制限条件１は、対象セルに対する代表サンプルを求めるために必要な条件である。制限条件２は、修正後のサンプルインク量データが、現実の印刷に使用できるインク量を表すことを保証するための条件である。なお、インクデューディ制限値とは、印刷媒体上の単位面積当たりに吐出可能なインク量であり、印刷媒体の種類毎に予め設定され、また、インクの滲みを考慮して設定される。通常のインクデューティ制限値は、各インクのインク量の最大値と、全インクの合計インク量の最大値とを含んでいる。なお、上記制限条件１，２以外の制限条件を追加するようにしてもよい。
【００６４】
こうしてサンプルインク量データが修正されると、修正されたサンプルインク量データを用いて上述したステップＳ３５，Ｓ３６の処理が再度実行される。こうして、ステップＳ３５〜Ｓ３７の処理が再帰的に実行され、評価基準を満たすサンプルが、その対象セルに対する代表サンプルとして選択される。なお、ある対象セルに関する再帰的処理を所定回行っても評価基準を満足するサンプルが得られない場合も考えられる。この場合には、その対象セルに関してそれまで調べられた複数のサンプルの中で、最も評価基準に近いサンプル（最も評価指数の良いサンプル）を代表サンプルとして選択してもよい。あるいは、その対象セルに関しては代表サンプルを選択しないものとしてもよい。
【００６５】
ステップＳ３８では、すべてのセルに関する処理が完了したか否かが判断され、完了していなければステップＳ３４に戻り、次のセルに関する処理を開始する。こうして、すべてのセルに関して処理が終了すると、ステップＳ３９において、選択された代表サンプルを用いてインクプロファイル１４２およびプリンタルックアップテーブル１８０が作成される。このステップＳ３９の処理は、図２のステップＳ７〜Ｓ９と同じものである。
【００６６】
このように、第３実施形態では、所定の測色値の色空間（上記の例ではＣＩＥＬＡＢ色空間）を複数のセルに分割し、所定の評価基準を満足する代表サンプルを各セル毎に再帰的プロセスを用いて追求し、これらの代表サンプルを用いてプロファイルやルックアップテーブルを作成する。従って、第１実施形態に比べて、全く実サンプルを含まないセルの数を低減することができる。この結果、より広い色域を有するプリンタルックアップテーブル１８０を得ることができる。また、色再現特性の点でもより優れたプリンタルックアップテーブル１８０を得ることが可能である。
【００６７】
Ｄ．分光プリンティングモデルの一例：
以下では、分光プリンティングモデルの一例としてセル分割ユール・ニールセン分光ノイゲバウアモデル(Cellular Yule-Nielsen Spectral Neugebauer Model)を説明する。このモデルは、よく知られた分光ノイゲバウアモデルとユール・ニールセンモデルとに基づいている。なお、以下の説明では、ＣＭＹの３種類のインクを用いた場合のモデルについて説明するが、これを任意の数のインクを用いたモデルに拡張することは容易である。セル分割ユール・ニールセン分光ノイゲバウアモデルについてははWyble and Berns Color Res. Appl. 25, 4-19, 2000に記載されており、その開示内容はすべての目的のためにここに組み込まれる。
【００６８】
図１３は、分光ノイゲバウアモデルを示す図である。分光ノイゲバウアモデルでは、任意の印刷物の分光反射率Ｒ（λ）は、以下の（８）式で与えられる。
【数８】
2007511175.tif 000010
ここで、aiはｉ番目の領域の面積率であり、Ri(λ)はｉ番目の領域の分光反射率である。添え字ｗ，ｃ，ｍ，ｙ，ｒ，ｇ，ｂ，ｋは、インクの無い領域（ｗ)と、シアンインクのみの領域（ｃ）と、マゼンタインクのみの領域（ｍ）と、イエローインクのみの領域（ｙ）と、マゼンタインクとイエローインクが吐出される領域（ｒ）と、イエローインクとシアンインクが吐出される領域（ｇ）と、シアンインクとマゼンタインクが吐出される領域（ｂ）と、ＣＭＹの３つのインクが吐出される領域（ｋ）をそれぞれ意味している。また、ｆc，ｆm，ｆyは、ＣＭＹ各インクを１種類のみ吐出したときにそのインクで覆われる面積の割合（「インク被覆率(Ink area coverage)」と呼ぶ）である。分光反射率Ri(λ)は、カラーパッチを分光反射率計で測定することによって取得される。
【００６９】
インク被覆率ｆc，ｆm，ｆyは、図１３（Ｂ）に示すマーレイ・デービスモデルで与えられる。マーレイ・デービスモデルでは、例えばシアンインクの面積率ｆcは、シアンのインク吐出量ｄcの非線形関数であり、１次元ルックアップテーブルの形で与えられる。インク被覆率がインク吐出量の非線形関数となる理由は、単位面積に少量のインクが吐出された場合にはインクが十分に広がるが、多量のインクが吐出された場合にはインクの広がりによってインクが重なり合うためにインクで覆われる面積があまり増加しないためである。
【００７０】
分光反射率に関するユール・ニールセンモデルを適用すると、上記（８）式は以下の（９ａ）式または（９ｂ）式に書き換えられる。
【数９】
2007511175.tif 000011
ここで、ｎは１以上の所定の係数であり、例えばｎ＝１０に設定することができる。（９ａ）式および（９ｂ）式は、ユール・ニールセン分光ノイゲバウアモデル(Yule-Nielsen Spectral Neugebauer Model)を表す式である。
【００７１】
セル分割ユール・ニールセン分光ノイゲバウアモデル(Cellular Yule-Nielsen Spectral Neugebauer Model)は、上述したユール・ニールセン分光ノイゲバウアモデルのインク色空間を複数のセルに分割したものである。
【００７２】
図１４（Ａ）は、セル分割ユール・ニールセン分光ノイゲバウアモデルにおけるセル分割の例を示している。ここでは、簡単のために、シアンのインク被覆率ｆcとマゼンタのインク被覆率ｆmの２つの軸を含む２次元空間でのセル分割を描いている。なお、これらの軸ｆc，ｆmは、インク吐出量ｄc，ｄmを示す軸と考えることもできる。白丸は、セル分割のグリッド点（「ノード」と呼ぶ）であり、２次元空間が９つのセルＣ１〜Ｃ９に分割されている。１６個のノードにおける印刷物（カラーパッチ）に対しては、分光反射率R00, R10, R20, R30, R01, R11…R33がそれぞれ予め決定される。
【００７３】
図１４（Ｂ）は、このセル分割に対応するインク被覆率ｆc（ｄ）の形状を示している。ここでは、１種類のインクのインク量の範囲０〜ｄmaxが３つの区間に分割されており、インク被覆率ｆc（ｄ）は、各区間毎に０から１まで単調に増加する曲線によって表されている。
【００７４】
図１４（Ｃ）は、図１４（Ａ）の中央のセルＣ５内にあるサンプルの分光反射率Rsmp(λ)の算出方法を示している。分光反射率Rsmp(λ)は、以下の（１０）式で与えられる。
【数１０】
2007511175.tif 000012
【００７５】
ここで、インク被覆率ｆc，ｆmは図１４（Ｂ）のグラフで与えられる値であり、このセルＣ５内で定義された値である。また、セルＣ５の４つの頂点における分光反射率R11(λ)，R12(λ)，R21(λ)，R22(λ)の値は、上記（１０）式に従ってサンプル分光反射率Rsmp(λ)を正しく与えるように調整されている。
【００７６】
このように、インク色空間を複数のセルに分割すれば、分割しない場合に比べてサンプルの分光反射率Rsmp(λ)をより精度良く算出することができる。図１５は、一実施例で採用されたセル分割のノード値を示している。この例に示されているように、セル分割のノード値は、各インク毎に独立に設定することが好ましい。
【００７７】
ところで、図１４（Ａ）に示すモデルにおいて、すべてのノードにおける分光反射率をカラーパッチの測定で得ることはできないのが普通である。この理由は、多量のインクを吐出すると滲みが発生してしまい、均一な色のカラーパッチを印刷できないからである。図１６は、測定できない分光反射率を求める方法を示している。これは、シアンとマゼンタの２種類のインクのみを使用する場合の例である。シアンとマゼンタの２種類のインクで印刷される任意のカラーパッチの分光反射率R(λ)は、以下の（１１）式で与えられる。
【数１１】
2007511175.tif 000013
【００７８】
上記（１１）式に含まれる複数のパラメータのうちで、シアンインクとマゼンタインクの両方が１００％吐出量であるときの分光反射率Rb(λ)のみが未知であり、他のパラメータの値は既知であると仮定する。このとき、（１１）式を変形すれば、以下の（１２）式が得られる。
【数１２】
2007511175.tif 000014
【００７９】
上述したように右辺の各項はすべて既知である。従って、（１２）式を解くことによって、未知の分光反射率Rb(λ)を算出することができる。印刷できない色の分光反射率の算定に関しては、R Balasubramanian, Optimization of the spectral Neugebauer model for printer characterization, J. Electronic Imaging 8(2), 156-166 (1999)に記載されており、その開示内容はすべての目的のためにここに組み込まれる。
【００８０】
シアンとマゼンタの２次色以外の他の２次色の分光反射率も同様にして求めることが可能である。また、複数の２次色の分光反射率が求まれば、複数の３次色の分光反射率も同様にして求めることができる。こうして、より高次色の分光反射率を順次求めてゆくことによって、セル分割されたインク色空間の各ノードにおける分光反射率をすべて求めることが可能である。
【００８１】
図１に示す分光プリンティングコンバータ１００は、図１４（Ａ）に示すようにセル分割されたインク色空間の各ノードにおける分光反射率の値と、図１４（Ｃ）に示すインク被覆率を示す１次元ルックアップテーブルとを有しており、これらを用いて任意のサンプルインク量データに対する分光反射率Rsmp(λ)を算出するように構成されている。
【００８２】
なお、一般に、印刷されたカラーパッチの分光反射率は、インクセットと印刷媒体とに依存する。従って、図１に示す分光プリンティングモデルコンバータ１００は、インクセットと印刷媒体との組合せ毎に作成される。また、インクプロファイル１４２やプリンタルックアップテーブル１８０も、インクセットと印刷媒体との組合せ毎に作成される。
【００８３】
Ｅ．変形例：
Ｅ１．変形例１：
上記実施形態では、インクとしてＣＭＹＫＯＧの６種類のインクを用いていたが、インクの種類としてはこれに限らず、任意の複数種類のインクを使用することができる。但し、オレンジインクやグリーンインク、ブルーインクなどの任意の特色インクを利用すれば、再現可能な分光反射率の形状の自由度がより大きくなる点で好ましい。
【００８４】
Ｅ２．変形例２：
上述した第１，第２実施形態では、測色値の色空間を複数のセルに分割し、各セルにおいて評価指数ＥＩが最良であるサンプルを代表サンプルとして選択し、一方、第３実施形態では、各セルに対して評価指数ＥＩが所定の基準を満足するようにサンプルを最適化することによって代表サンプルを決定していた。しかし、色変換プロファイルを作成する際に用いる複数の代表サンプルを選択する方法は、これらに限らず、一般には、評価指数ＥＩに基づいて複数の代表サンプルを選択することが可能である。例えば、測色値の色空間を複数のセルに分割せずに複数の代表サンプルを選択することも可能である。具体的には、測色値の色空間内で複数のグリッド点（ノード）を設定し、各ノードの近傍において所定の評価基準を満足するサンプルをそのノードの代表サンプルとして選択することができる。
【００８５】
本願発明が詳細に記述され図示されてきたが、これらは単なる説明のための例示であり、限定として解釈されるべきではなく、本願発明の精神とスコープはクレームの用語のみによって限定されることが理解される。
【図面の簡単な説明】
【００８６】
【図１】本発明の第１実施形態としての全体のシステム構成を示すブロック図である。
【図２】第１実施形態の処理手順を示すフローチャートである。
【図３】一実施例におけるＣＩＥＬＡＢ空間内のサンプル色を示すグラフである。
【図４】図２のステップＳ６の詳細手順を示すフローチャートである。
【図５】ステップＳ７における不等間隔補間を示す図である。
【図６】ステップＳ８におけるガマットマッピングを示す図である。
【図７】一実施例と比較例における色差の算出結果を示すグラフである。
【図８】一実施例と比較例における非色恒常性指数（Color Inconstancy Index）の算出結果を示すグラフである。
【図９】本発明の第２実施形態としてのシステム構成を示すブロック図である。
【図１０】第２実施形態の処理手順を示すフローチャートである。
【図１１】本発明の第４実施形態としての全体のシステム構成を示すブロック図である。
【図１２】第４実施形態の処理手順を示すフローチャートである。
【図１３】分光ノイゲバウアモデルを示す図である。
【図１４】セル分割ユール・ニールセン分光ノイゲバウアモデル(Cellular Yule-Nielsen Spectral Neugebauer Model)を示す図である。
【図１５】セル分割ユール・ニールセン分光ノイゲバウアモデルにおけるセル分割の格子点座標を示す図である。
【図１６】セル分割モデルを示す図である。
【図１】
2007511175.tif 000015
【図２】
2007511175.tif 000016
【図３】
2007511175.tif 000017
【図４】
2007511175.tif 000018
【図５】
2007511175.tif 000019
【図６】
2007511175.tif 000020
【図７】
2007511175.tif 000021
【図８】
2007511175.tif 000022
【図９】
2007511175.tif 000023
【図１０】
2007511175.tif 000024
【図１１】
2007511175.tif 000025
【図１２】
2007511175.tif 000026
【図１３】
2007511175.tif 000027
【図１４】
2007511175.tif 000028
【図１５】
2007511175.tif 000029
【図１６】
2007511175.tif 000030
【国際調査報告】
2007511175.tif 000031
2007511175.tif 000032
2007511175.tif 000033