(57)【要約】【課題】入力表色系の座標値に対応する複数種類のインクで構成されるインク表色系のインク量セットを探索する最適化処理において、落ち着かせるミニマムポイントをコントロール可能にする。【解決手段】入力表色系の座標値を複数種類のインクで構成されるインク表色系のインク量セットに変換するために色変換テーブルを作成するにあたり、上記入力表色系の座標値に対応する上記インク表色系のインク量を、各々に重み付けした複数の画質評価指数の組み合わせにより表現される目的関数Ｅを用いた最適化を複数回行うことにより決定する最適化工程と、上記複数の画質評価指数に対する重みを上記最適化の各回ごとに変化させる重み管理工程と、を行う。【選択図】図１４
2011223298.tif 000002
【特許請求の範囲】
【請求項１】
入力表色系の座標値を複数種類のインクで構成されるインク表色系のインク量セットに変換するための色変換テーブルの作成方法であって、
上記入力表色系の座標値に対応する上記インク表色系のインク量を、各々に重み付けした複数の画質評価指数の組み合わせにより表現される所定の目的関数を用いた最適化を複数回行うことにより決定する最適化工程と、
上記複数の画質評価指数に対する重みを上記最適化の各回ごとに変化させる重み管理工程と、
上記最適化工程によって決定されたインク量に基づいて、上記入力表色系の座標値を上記インク表色系のインク量に変換するための色変換テーブルを作成する色変換テーブル作成工程と、
を備えることを特徴とする色変換テーブルの作成方法。
【請求項２】
上記重み管理工程では、上記最適化工程の最適化の回数が進むにつれて重要度の高い上記画質評価指数ほど大きくなるように重みを徐々に変化させる請求項１に記載の色変換テーブルの作成方法。
【請求項３】
上記重み管理工程では、上記最適化工程の初期において重要度の高い上記画質評価指数ほど大きくなるように上記重みを設定し、上記最適化工程の終了時において重要度の低い上記画質評価指数ほど、上記重みが上記最適化工程の初期における上記重みに対して相対的に大きくなるように徐々に変化させる請求項１または請求項２に記載の色変換テーブルの作成方法。
【請求項４】
上記重み管理工程では、上記最適化工程の初期において複数の上記画質評価指数で上記重みを平均的に設定し、上記最適化工程の終了時において重要度の高い上記画質評価指数ほど、上記重みが上記最適化工程の初期における上記重みに対して相対的に大きくなるように徐々に変化させる請求項１または請求項２に記載の色変換テーブルの作成方法。
【請求項５】
入力表色系の座標値を複数種類のインクで構成されるインク表色系のインク量セットに変換するための色変換テーブルの作成プログラムであって、
上記入力表色系の座標値に対応する上記インク表色系のインク量を、各々に重み付けした複数の画質評価指数の組み合わせにより表現される所定の目的関数を用いた最適化を複数回行うことにより決定する最適化機能と、
上記複数の画質評価指数に対する重みを上記最適化の各回ごとに変化させる重み管理機能と、
上記最適化機能によって決定されたインク量に基づいて、上記入力表色系の座標値を上記インク表色系のインク量に変換するための色変換テーブルを作成する色変換テーブル作成機能と、
を備えることを特徴とする色変換テーブルの作成プログラム。
【請求項６】
請求項１〜請求項４のいずれか１項に記載の色変換テーブルの作成方法にて作成された色変換テーブルを記憶した記憶部を備え、上記記憶部に記憶された色変換テーブルに基づいて印刷データにおける上記入力表色系の座標値を上記インク表色系のインク量セットに変換し、該インク量セットに基づいて印刷する印刷装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、色変換テーブルの作成方法、色変換テーブルの作成プログラムおよび印刷装置に関し、特に、入力表色系の座標値を複数種類のインクで構成されるインク表色系のインク量セットに変換するための色変換テーブルの作成方法、色変換テーブルの作成プログラム、および上記作成方法により作成された色変換テーブルを備えた印刷装置に関する。
【背景技術】
【０００２】
色変換対応情報は、入力表色系と出力表色系との対応関係を示す情報であり、例えば、色変換テーブルや色変換関数などの形式で利用される。色変換テーブルの入力表色系の座標値は、入力表色系の色空間内の点の位置を表しており、出力表色系の座標値は出力表色系の色空間内の点の位置を表している。この出力表色系の色空間内の点は、上記インク表色系のインク量セットに対応する。なお、本明細書において、任意の色空間内の点を「色点」または「格子点」とも呼ぶ。また、色変換テーブルに登録されている入力値で表される色点および出力値で表される色点を、それぞれ「入力格子点」および「出力格子点」とも呼ぶ。
【０００３】
色変換テーブルの入力格子点に対応する出力格子点を最適化する技術として、例えば本願の出願人により開示された特許文献１に記載されたものがある。この文献に記載された最適化では、Ｌａｂ表色系の格子点を移動させた後、複数の画質評価指数の組み合わせにより構成される目的関数を用いた最適化処理を利用して、入力表色系の格子点に対応する最適なインク量を決定している。この最適なインク量は、目的関数を最小化するインク量として決定される。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００８−２６３５７９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、上述した特許文献１に記載の最適化処理では、最適なインク量までたどり着かずにローカルミニマムに落ち着いてしまうことがあった。また、逆に、特定のローカルミニマムに落ち着かせたいにも関わらず、最小値に落ち着いてしまったり他のローカルミニマムに陥ったりすることがあった。
本発明は、入力表色系の座標値に対応する複数種類のインクで構成されるインク表色系のインク量セットを探索する最適化処理において、最適化されるミニマムポイントをコントロール可能にした色変換テーブルの作成方法、色変換テーブルの作成プログラム、および作成された色変換テーブルを実装した印刷装置の提供を目的とする。
【課題を解決するための手段】
【０００６】
本発明は、上述の課題の少なくとも一部を解決するためになされたものであり、以下の形態として実現することが可能である。
【０００７】
本発明の一態様として、入力表色系の座標値を複数種類のインクで構成されるインク表色系のインク量セットに変換するための色変換テーブルの作成方法であって、上記入力表色系の座標値に対応する上記インク表色系のインク量を、各々に重み付けした複数の画質評価指数の組み合わせにより表現される所定の目的関数を用いた最適化を複数回行うことにより決定する最適化工程と、上記複数の画質評価指数に対する重みを上記最適化の各回ごとに変化させる重み管理工程と、上記最適化工程によって決定されたインク量に基づいて、上記入力表色系の座標値を上記インク表色系のインク量に変換するための色変換テーブルを作成する色変換テーブル作成工程と、を備える構成とすることができる。
この構成によれば、入力表色系の座標値に対応する複数種類のインクで構成されるインク表色系のインク量セットを探索する最適化処理において、最適化の結果として得られる最適解をコントロールすることが可能となる。
【０００８】
また、本発明の選択的な一態様として、上記重み管理工程では、上記最適化工程の最適化の回数が進むにつれて、重要度の高い上記画質評価指数ほど大きくなるように重みを徐々に変化させる構成としてもよい。
この構成によれば、初期の重み係数を与えられた目的関数により得られる最適解から、最終的な重み係数を与えられた目的関数により得られる最適解へと、徐々に最適解を誘導することができる。例えば、固定された重みを複数の画質評価指数のそれぞれに加重して表現された目的関数を用いて最適化処理を行ったときにローカルミニマムに陥りやすい場合には、初期においてこのローカルミニマムの深さを浅くするような重み付けを行ったり、このローカルミニマム以外の箇所により深いミニマムポイントが生じるような重み付けを行うことにより、特定のローカルミニマムに陥りにくくコントロールすることができる。
【０００９】
また、本発明の選択的な一態様として、上記重み管理工程では、上記最適化工程の初期において重要度の高い上記画質評価指数ほど大きくなるように上記重みを設定し、上記最適化工程の終了時に重要度の低い上記画質評価指数ほど、上記重みが上記最適化工程の初期における上記重みに対して相対的に大きくなるように徐々に変化させる構成としてもよい。
この構成によれば、最終的な重み係数を与えられた目的関数において重要度の高い画質評価指数に起因して形成されるローカルミニマムを、初期の重み係数を与えられた目的関数において強調する（より深くする）ことにより、初期の最適化処理により得られる最適解を該ローカルミニマムに誘導することができる。その後、徐々に、初期の重み係数から最終的な重み係数へと、重みが上記最適化工程の初期における上記重みに対して相対的に大きくなるように変化させていくことにより、最適解が上記ローカルミニマムから移動しにくく、最終的な最適化の結果として得られる最適解が上記ローカルミニマムの近辺に収束しやすくなるようにコントロールすることができる。
【００１０】
また、本発明の選択的な一態様として、上記重み管理工程では、上記最適化工程の初期において複数の上記画質評価指数で上記重みを平均的に設定し、上記最適化工程の終了時において重要度の高い上記画質評価指数ほど、上記重みが上記最適化工程の初期における上記重みに対して相対的に大きくなるように徐々に変化させる構成としてもよい。
この構成によれば、初期の最適化処理により得られる最適解を、各画質評価指数をバランスよく考慮した平均的な最適解へと誘導することができる。その後、徐々に、初期の重み係数から最終的な重み係数へと、重みが上記最適化工程の初期における上記重みに対して相対的に大きくなるように変化させていくことにより、最終的に得られる最適解を各画質評価指数をバランスよく考慮した平均的な最適解の近辺にあるミニマムポイントへ収束しやすくなるようにコントロールすることができる。
【００１１】
上述した色変換テーブルの作成方法は、他の方法の一環として実施されたり各工程に対応する手段を備えた色変換テーブルの作成装置として実現されたりする等の各種の態様を含む。また、本発明は上記色変換テーブルの作成装置を備える色変換テーブルの作成システム、前述した方法の構成に対応した機能をコンピュータに実現させる色変換テーブルの作成プログラム、該プログラムを記録したコンピュータ読み取り可能な記録媒体、等としても実現可能である。さらに、上記作成方法等によって作成された色変換テーブルを実装された印刷装置が実行する印刷によって作成される印刷物において、上記画質評価指数の重みのコントロール結果を反映した色を再現することができるため、該印刷装置も本発明の一態様であると言える。すなわち、これら色変換テーブルの作成装置、色変換テーブルの作成システム、色変換テーブルの作成プログラム、該プログラムを記録した媒体、印刷装置、の発明も、前述した作用、効果を奏する。
【図面の簡単な説明】
【００１２】
【図１】本発明の一実施例におけるルックアップテーブル作成装置の構成を示すブロック図である。
【図２】実施例の全体処理手順を示すフローチャートである。
【図３】ベース３Ｄ−ＬＵＴを作成する場合の処理内容を示す説明図である。
【図４】入力表色系であるＲＧＢ表色系の色点とＬａｂ表色系の色点との対応関係を示す説明図である。
【図５】図２のステップＳ１００〜Ｓ３００によってベース４Ｄ−ＬＵＴを作成する場合の処理内容を示す説明図である。
【図６】ベースＬＵＴを用いた色補正ＬＵＴの作成方法を示す説明図である。
【図７】実施例のスムージング処理に利用される力学モデルを示す説明図である。
【図８】スムージング処理の典型的な処理手順を示すフローチャートである。
【図９】図８のステップＴ１００の詳細手順を示すフローチャートである。
【図１０】図８のステップＴ１２０〜Ｔ１５０の処理内容を示す説明図である。
【図１１】ＬＵＴ作成条件の設定処理の流れを示すフローチャートである。
【図１２】ステップＴ９１にて表示されるＵＩ画面を示している。
【図１３】総合的な最適解へ最適化する場合の重みバランスを説明する図である。
【図１４】最も重要な画質要素に関する最適化を優先する場合の重みバランスを説明する図である。
【図１５】最適化処理（図８のステップＴ１３０）の詳細手順を示すフローチャートである。
【発明を実施するための形態】
【００１３】
次に、本発明の実施の形態を以下の順序で説明する。
（１）装置構成と全体処理手順：
（２）力学モデル：
（３）スムージング処理（平滑化／最適化処理）：
（４）ＬＵＴ作成条件の設定：
（５）最適化処理：
（６）変形例：
【００１４】
（１）装置構成と全体処理手順：
図１は、本発明の一実施例におけるルックアップテーブル作成装置の構成を示すブロック図である。この装置は、ベースＬＵＴ作成モジュール１００と、色補正ＬＵＴ作成モジュール２００と、フォワードモデルコンバーター３００と、ＬＵＴ格納部４００と、を備えている。「ＬＵＴ」は、色変換プロファイルとしてのルックアップテーブルの略語である。本実施例においては、この色変換プロファイルが色変換テーブルを構成する。これらのモジュール１００，２００やコンバーター３００の機能は、メモリに格納されたコンピュータプログラムであるＬＵＴ作成プログラムをコンピュータが実行することによってそれぞれ実現される。また、ＬＵＴ格納部４００は、ハードディスク装置などの記録媒体によって実現される。
【００１５】
ベースＬＵＴ作成モジュール１００は、スムージング処理初期値設定モジュール１２０と、スムージング処理モジュール１３０と、テーブル作成モジュール１４０とを有している。スムージング処理モジュール１３０は、色点移動モジュール１３２と、最適化条件設定モジュール１３３と、インク量最適化モジュール１３４と、画質評価指数算出モジュール１３６とを有している。フォワードモデルコンバーター３００は、分光プリンティングモデルコンバーター３１０と、色算出部３２０とを備えている。これらの各部の機能については後述する。
【００１６】
ＬＵＴ格納部４００は、インバースモデル初期ＬＵＴ４１０や、ベース３Ｄ−ＬＵＴ５１０，ベース４Ｄ−ＬＵＴ５２０，色補正３Ｄ−ＬＵＴ６１０，色補正４Ｄ−ＬＵＴ６２０などを格納するためのものである。但し、インバースモデル初期ＬＵＴ４１０以外のＬＵＴは、ベースＬＵＴ作成モジュール１００や色補正ＬＵＴ作成モジュール２００によって作成されるものである。ベース３Ｄ−ＬＵＴ５１０は、ＲＧＢ表色系を入力とし、インク量を出力とする色変換ルックアップテーブルである。一方、ベース４Ｄ−ＬＵＴ５２０は、ＣＭＹＫ表色系を入力とし、インク量を出力とする色変換ルックアップテーブルである。なお、「３Ｄ」や「４Ｄ」は、入力値の数を意味している。これらのベース３Ｄ−ＬＵＴ５１０，ベース４Ｄ−ＬＵＴ５２０は、例えば色補正３Ｄ−ＬＵＴ６１０，色補正４Ｄ−ＬＵＴ６２０を作成する際に使用される。「ベースＬＵＴ」という名前は、色補正ＬＵＴを作成する基礎として用いられるからである。色補正３Ｄ−ＬＵＴ６１０，色補正４Ｄ−ＬＵＴ６２０は、標準的な機器依存表色系（例えばｓＲＧＢ表色系やＪＡＰＡＮ  ＣＯＬＯＲ  ２００１表色系）を、特定のプリンターのインク量に変換するためのルックアップテーブルである。インバースモデル初期ＬＵＴ４１０については後述する。
【００１７】
図２は、実施例の全体処理手順を示すフローチャートである。図３（Ａ）〜（Ｃ）は、図２のステップＳ１００〜Ｓ３００によってベース３Ｄ−ＬＵＴを作成する場合の処理内容を示す説明図である。ステップＳ１００では、フォワードモデルコンバーター３００とインバースモデル初期ＬＵＴ４１０とが準備される。ここで、「フォワードモデル」とは、インク量を機器非依存表色系の色彩値（測色値）に変換する変換モデルを意味し、「インバースモデル」とは、逆に、機器非依存表色系の色彩値をインク量に変換する変換モデルを意味している。実施例では、機器非依存表色系としてＣＩＥ−Ｌａｂ表色系を使用する。なお、以下では、ＣＩＥ−Ｌａｂ表色系の色彩値を、単に「Ｌ*ａ*ｂ*値」または「Ｌａｂ値」とも呼ぶ。
【００１８】
図３（Ａ）に示すように、フォワードモデルコンバーター３００は、前段に分光プリンティングモデルコンバーター３１０、後段に色算出部３２０を備える。
分光プリンティングモデルコンバーター３１０は、複数種類のインクのインク量を、そのインク量に応じて印刷されるカラーパッチの分光反射率Ｒ（λ）に変換する。なお、本明細書において「カラーパッチ」という用語は、有彩色のパッチに限らず、無彩色のパッチも含む広い意味で使用される。本実施例では、シアン（Ｃ）、マゼンタ（Ｍ）、イエロー（Ｙ）、ブラック（Ｋ）、淡シアン（Ｌｃ）、淡マゼンタ（Ｌｍ）、淡ブラック（Ｌｋ）、極淡ブラック（ＬＬｋ）の８種類のインクを利用可能なカラープリンターを想定しており、分光プリンティングモデルコンバーター３１０もこの８種類のインクのインク量を入力としている。但し、プリンターで使用する複数種類のインクとしては、任意のインクセットを利用することが可能である。なお、分光プリンティングモデルコンバーター３１０を作成する方法としては、例えば特表２００７−５１１１７５号公報に記載された方法を採用することが可能である。
色算出部３２０は、分光反射率Ｒ（λ）からＬａｂ表色系の色彩値を算出する。この色彩値の算出には、予め選択された光源（例えば標準の光Ｄ５０）がカラーパッチの観察条件として使用される。
【００１９】
インバースモデル初期ＬＵＴ４１０は、Ｌ*ａ*ｂ*値を入力とし、インク量を出力とするルックアップテーブルである。このインバースモデル初期ＬＵＴ４１０は、例えば、Ｌ*ａ*ｂ*空間を複数の小セルに区分し、各小セル毎に最適なインク量を選択して登録したものである。この選択は、例えば、そのインク量で印刷されるカラーパッチの画質を考慮して行われる。一般に、或る１つのＬ*ａ*ｂ*値を再現するインク量の組み合わせは多数存在する。そこで、インバースモデル初期ＬＵＴ４１０では、ほぼ同じＬ*ａ*ｂ*値を再現する多数のインク量の組み合わせの中から、画質等の所望の観点から最適なインク量を選択したものが登録されている。このインバースモデル初期ＬＵＴ４１０の入力値であるＬ*ａ*ｂ*値は各小セルの代表値である。一方、出力値であるインク量はそのセル内のいずれかのＬ*ａ*ｂ*値を再現するものである。
【００２０】
従って、このインバースモデル初期ＬＵＴ４１０では、入力値であるＬ*ａ*ｂ*値と出力値であるインク量とが厳密に対応したものとなっておらず、出力値のインク量をフォワードモデルコンバーター３００でＬ*ａ*ｂ*値に変換すると、インバースモデル初期ＬＵＴ４１０の入力値とは多少異なる値が得られる。むろん、インバースモデル初期ＬＵＴ４１０として、入力値と出力値とが完全に対応するものを利用してもよい。また、変換式等でＬ*ａ*ｂ*値をインク量に変換できれば、インバースモデル初期ＬＵＴ４１０を用いずにベースＬＵＴを作成することも可能である。なお、小セル毎に最適なインク量を選択してインバースモデル初期ＬＵＴ４１０を作成する方法としては、例えば前記特表２００７−５１１１７５号公報に記載された方法を採用することが可能である。
【００２１】
図２のステップＳ２００では、ベースＬＵＴ作成のための初期入力値がユーザーによって設定される。図３（Ｂ）は、ベース３Ｄ−ＬＵＴ５１０の構成とその初期入力値設定の例を示している。ベース３Ｄ−ＬＵＴ５１０の入力値としては、ＲＧＢの各値として予め定められたほぼ等間隔の値が設定される。１組のＲＧＢ値はＲＧＢ色空間内の点を表していると考えられるので、１組のＲＧＢ値を「入力格子点」とも呼ぶ。
【００２２】
ステップＳ２００においては、複数の入力格子点のうちから予め選択されたいくつかの少数の入力格子点に対するインク量の初期値がユーザーによって入力される。この初期入力値が設定される入力格子点としては、ＲＧＢ色空間における３次元色立体の頂点に相当する入力格子点を少なくとも選択することが好ましい。この３次元色立体の頂点では、ＲＧＢの各値がその定義範囲の最小値または最大値を取る。
【００２３】
具体的には、ＲＧＢの各値を８ビットで表現した場合には、（Ｒ，Ｇ，Ｂ）＝（０，０，０）、（０，０，２５５）、（０，２５５，０）、（２５５，０，０）、（０，２５５，２５５）、（２５５，０，２５５）、（２５５，２５５，０）、（２５５，２５５，２５５）である８つの入力格子点に関してインク量の初期入力値が設定される。なお、（Ｒ，Ｇ，Ｂ）＝（２５５，２５５，２５５）の入力格子点に対するインク量は、すべて０に設定される。他の入力格子点に対するインク量の初期入力値は任意であり、例えば０に設定される。図３（Ｂ）の例では、（Ｒ，Ｇ，Ｂ）＝（０，０，３２）の入力格子点に対するインク量が０以外の値になっているが、これはベース３Ｄ−ＬＵＴ５１０が完成したときの値である。
【００２４】
図２のステップＳ３００では、スムージング処理モジュール１３０（図１）が、ステップＳ２００で設定された初期入力値に基づいてスムージング処理（平滑化／最適化処理）を実行する。図３（Ｃ）は、ステップＳ３００の処理内容を示している。図３（Ｃ）の左側には、スムージング処理前の状態における複数の色点の分布が２重丸と白丸とで示されている。これらの色点は、Ｌ*ａ*ｂ*空間における３次元色立体ＣＳを構成している。各色点のＬ*ａ*ｂ*座標値は、ベース３Ｄ−ＬＵＴ５１０の複数の入力格子点におけるインク量を、フォワードモデルコンバーター３００（図３（Ａ））を用いてＬ*ａ*ｂ*値に変換した値である。上述したように、ステップＳ２００では一部の少数の入力格子点についてのみインク量の初期入力値が設定される。そこで、他の入力格子点に対するインク量の初期値は、初期入力値に基づいてスムージング処理初期値設定モジュール１２０（図１）によって設定される。この初期値設定方法については後述する。
【００２５】
Ｌａｂ表色系の３次元色立体ＣＳは、以下の８つの頂点（図３（Ｃ）の２重丸の点）を有している。
・点ＰK：（Ｒ，Ｇ，Ｂ）＝（０，０，０）に対応する紙黒点。
・点ＰW：（Ｒ，Ｇ，Ｂ）＝（２５５，２５５，２５５）に対応する紙白点。
・点ＰC：（Ｒ，Ｇ，Ｂ）＝（０，２５５，２５５）に対応するシアン点。
・点ＰM：（Ｒ，Ｇ，Ｂ）＝（２５５，０，２５５）に対応するマゼンタ点。
・点ＰY：（Ｒ，Ｇ，Ｂ）＝（２５５，２５５，０）に対応するイエロー点。
・点ＰR：（Ｒ，Ｇ，Ｂ）＝（２５５，０，０）に対応するレッド点。
・点ＰG：（Ｒ，Ｇ，Ｂ）＝（０，２５５，０）に対応するグリーン点。
・点ＰB：（Ｒ，Ｇ，Ｂ）＝（０，０，２５５）に対応するブルー点。
【００２６】
図３（Ｃ）の右側は、スムージング処理後の色点の分布を示している。スムージング処理は、Ｌ*ａ*ｂ*空間における複数の色点を移動させて、それらの色点の分布を等間隔に近い平滑なものにする処理である。スムージング処理では、さらに、移動後の各色点のＬ*ａ*ｂ*値を再現するために最適なインク量も決定される。この最適なインク量がベース３Ｄ−ＬＵＴ５１０の出力値として登録されると、ベース３Ｄ−ＬＵＴ５１０が完成する。
【００２７】
図４（Ａ）〜（Ｃ）は、入力表色系の色点（すなわち入力格子点）とＬａｂ表色系の色点との対応関係を示している。Ｌａｂ表色系の３次元色立体ＣＳの頂点は、ベース３Ｄ−ＬＵＴ５１０の入力表色系の３次元色立体の頂点と一対一に対応している。また、各頂点を結ぶ辺（稜線）も、両方の色立体で互いに対応しているものと考えることができる。スムージング処理前のＬａｂ表色系の各色点は、ベース３Ｄ−ＬＵＴ５１０の入力格子点にそれぞれ対応付けられており、従って、スムージング処理後のＬａｂ表色系の各色点もベース３Ｄ−ＬＵＴ５１０の入力格子点にそれぞれ対応付けられる。なお、ベース３Ｄ−ＬＵＴ５１０の入力格子点はスムージング処理によって変化しない。スムージング処理後のＬａｂ表色系の３次元色立体ＣＳは、ベース３Ｄ−ＬＵＴ５１０の出力表色系を構成するインクセットで再現可能な色域（カラーガマット）の全体に対応している。従って、ベース３Ｄ−ＬＵＴ５１０の入力表色系は、このインクセットで再現可能な色域の全体を表す表色系としての意義を有している。
【００２８】
ベース３Ｄ−ＬＵＴ５１０を作成する際に、Ｌ*ａ*ｂ*空間においてスムージング処理を行う理由は以下の通りである。ベース３Ｄ−ＬＵＴ５１０では、なるべく大きな色域を再現できるように出力表色系のインク量を設定したいという要望がある。特定のインクセットで再現可能な色域は、インクデューティ制限（一定面積に吐出可能なインク量の制限）などの所定の制限条件を考慮して決定される。一方、上述したフォワードモデルコンバーター３００は、これらの制限条件が考慮されておらず、再現可能な色域とは無関係に作成されている。そこで、スムージング処理の際にインクデューティ制限等の制限条件を考慮してＬ*ａ*ｂ*空間内の色点の取り得る範囲を決定すれば、特定のインクセットで再現可能な色域を決定することが可能となる。なお、色点の移動を行うアルゴリズムとしては、例えば、後述する力学モデルを使用したものが利用される。
【００２９】
図２のステップＳ４００では、スムージング処理の結果を用いて、テーブル作成モジュール１４０がベース３Ｄ−ＬＵＴ５１０を作成する。すなわち、テーブル作成モジュール１４０は、各入力格子点に対応付けられたＬａｂ表色系の色点を再現するための最適なインク量をベース３Ｄ−ＬＵＴ５１０（図３（Ｃ））の出力値として登録する。なお、スムージング処理では、その計算負荷を軽減するために、ベース３Ｄ−ＬＵＴ５１０の入力格子点の一部のみに対応する色点のみを処理対象として選択することも可能である。例えば、ベース３Ｄ−ＬＵＴ５１０の入力格子点におけるＲＧＢ値の間隔が１６である場合に、スムージング処理の対象となる入力格子点におけるＲＧＢ値の間隔を３２に設定すれば、スムージング処理の負荷を半減することができる。この場合には、テーブル作成モジュール１４０は、スムージング処理結果を補間することによってベース３Ｄ−ＬＵＴ５１０のすべての入力格子点に対するインク量を決定して登録する。
【００３０】
図５（Ａ）〜（Ｃ）は、図２のステップＳ１００〜Ｓ３００によってベース４Ｄ−ＬＵＴ５２０を作成する場合の処理内容を示す説明図である。図５（Ａ）は、図３（Ａ）と同じである。図５（Ｂ）に示すベース４Ｄ−ＬＵＴ５２０は、入力がＣＭＹＫ表色系である点が図３（Ｂ）に示したベース３Ｄ−ＬＵＴ５１０と異なっている。このベース４Ｄ−ＬＵＴ５２０の初期入力値としては、（Ｃ，Ｍ，Ｙ，Ｋ）＝（０，０，０，０），（０，０，２５５，０），（０，２５５，０，０），（０，２５５，２５５，０），（２５５，０，０，０），（２５５，０，２５５，０），（２５５，２５５，０，０），（２５５，２５５，２５５，０），（０，０，０，２５５），（０，０，２５５，２５５），（０，２５５，０，２５５），（０，２５５，２５５，２５５），（２５５，０，０，２５５），（２５５，０，２５５，２５５），（２５５，２５５，０，２５５），（２５５，２５５，２５５，２５５）である１６個の入力格子点に関してインク量の初期値が設定される。他の入力格子点に対するインク量の初期入力値は任意であり、例えば０に設定される。
【００３１】
図５（Ｃ）は、スムージング処理の様子を示している。なお、Ｌ*ａ*ｂ*空間においてベース４Ｄ−ＬＵＴ５２０に対応する色立体としては、図５（Ｃ）の右端に示すように、入力値のうちのＫ値のそれぞれの値に対して１つの３次元色立体ＣＳが存在する。この例では、Ｋ＝０の色立体とＫ＝３２の色立体とを含む複数の色立体ＣＳが図示されている。本明細書では、これらの個々の色立体ＣＳを「Ｋレイヤ」とも呼ぶ。この理由は、各色立体ＣＳが、ＣＭＹＫ値のうちのＫ値が一定でＣ，Ｍ，Ｙ値が可変である入力層に対応するものと考えることができるからである。複数の色立体ＣＳは、Ｋ値が大きいほど暗い色域を表現するものとなっている。これらの複数の色立体ＣＳは、入力表色系のＫ値が大きいほど濃ブラックインクＫのインク量が多くなるように濃ブラックインクＫのインク量を決定することによって実現できる。上述したように、再現可能な色域はインクデューティ制限値によって制限される。このインクデューティ制限値としては、個々のインクのインク量と、全インクの合計インク量と、の２つの制限値が課されるのが普通である。一方、暗い色を再現する方法としては、濃ブラックインクＫなどの無彩色インクを用いる方法と、コンポジットブラックを用いる方法とがある。しかし、コンポジットブラックは合計インク量が多くなるので、濃ブラックインクＫに比べてインクデューティ制限値に抵触する可能性が高く、暗い色を再現するのには不利である。従って、入力表色系のＫ値が大きく濃ブラックインクＫのインク量が多い色立体の方が、入力表色系のＫ値が小さく濃ブラックインクＫのインク量が少ない色立体に比べてより暗い色を再現することが可能となる。
【００３２】
図６（Ａ），（Ｂ）は、ベースＬＵＴをベースとして行う色補正ＬＵＴの作成方法の説明図である。図６（Ａ）に示すように、ベース３Ｄ−ＬＵＴ５１０は、ＲＧＢ値とインク量Ｉjとの対応関係を規定している。このインク量Ｉjは、図３（Ｂ）に示した８種類のインクのインク量を表しているため、インク量Ｉjの添え字ｊは１〜８である。ベース３Ｄ−ＬＵＴ５１０を参照してＲＧＢ値から変換されたインク量Ｉjは、フォワードモデルコンバーター３００によってＬ*ａ*ｂ*値に変換される。
一方、ｓＲＧＢ値は、既知の変換式に従ってＬ*ａ*ｂ*値に変換される。この変換後のＬ*ａ*ｂ*値は、その色域が、フォワードモデルコンバーター３００で変換されたＬ*ａ*ｂ*値の色域と一致するようにガマットマッピングされる。
【００３３】
一方、ベース３Ｄ−ＬＵＴ５１０とフォワードモデルコンバーター３００を通じて、ＲＧＢ値から変換したＬ*ａ*ｂ*値を、逆方向ルックアップテーブルとして、逆変換ＬＵＴ５１１を作成する。ガマットマッピングされたＬ*ａ*ｂ*値は、この逆変換ＬＵＴ５１１によってＲＧＢ値に変換される。このＲＧＢ値は、さらに、ベース３Ｄ−ＬＵＴ５１０によってインク量Ｉjに再度変換される。この最後のインク量Ｉjと最初のｓＲＧＢ値の対応関係をルックアップテーブルに登録することによって、色補正３Ｄ−ＬＵＴ６１０を作成することができる。この色補正３Ｄ−ＬＵＴ６１０は、ｓＲＧＢ表色系をインク表色系に変換する色変換テーブルである。
【００３４】
図６（Ｂ）は、色補正４Ｄ−ＬＵＴ６２０を作成する方法を示している。図６（Ａ）との違いは、ベース３Ｄ−ＬＵＴ５１０およびその逆変換ＬＵＴ５１１の代わりに、ベース４Ｄ−ＬＵＴ５２０およびその逆変換ＬＵＴ５２１を利用している点と、ｓＲＧＢ表色系をＬ*ａ*ｂ*値に変換する既知変換式の代わりにＪＡＰＡＮ  ＣＯＬＯＲ表色系（図中では「ｊＣＭＹＫ」と記したもの）をＬ*ａ*ｂ*値に変換する既知変換式を使用している点である。良く知られているように、ＪＡＰＡＮ  ＣＯＬＯＲは、ＣＭＹＫの４色で構成される表色系である。なお、図６（Ｂ）の方法では、逆変換ＬＵＴ５２１において、Ｌ*ａ*ｂ*値からＣＭＹＫ値に変換する際に、既知変換前の最初のｊＣＭＹＫ値のＫ値から、逆変換ＬＵＴ５２１のＫレイヤ（Ｋ値が一定を取る部分）が選択される。従って、色補正４Ｄ−ＬＵＴ６２０として、ベース４Ｄ−ＬＵＴ５２０のうちのＫレイヤにおける特性を反映したものを作成することが可能である。以上のようにテーブル作成モジュール１４０が実行するステップＳ４００が、本実施例において色変換テーブル作成工程を構成する。
【００３５】
なお、通常は、ベース３Ｄ−ＬＵＴ５１０，ベース４Ｄ−ＬＵＴ５２０がプリンタドライバに実装されており、色補正ＬＵＴの作成処理以外の処理にも活用されているが、ここでは他の活用例の説明は省略する。以下では、実施例のスムージング処理（平滑化／最適化処理）に利用される力学モデルについて簡単に説明した後に、スムージング処理の処理手順、および、最適化処理の内容について順次説明する。
【００３６】
（２）力学モデル：
図７は、本実施例のスムージング処理（平滑化／最適化処理）に利用される力学モデルを示す説明図である。ここでは、Ｌ*ａ*ｂ*色空間内に複数の色点（白丸および２重丸）が配列されている様子を示している。ただし、ここでは説明の便宜上、色点の配置を２次元的に描いている。この力学モデルでは、着目色点ｇに対して次式の仮想的な力Ｆpgがかかるものと仮定する。
【数１】
2011223298.tif 000003
ここで、Ｆgは着目色点ｇが隣接色点ｇｎ（ｎは１〜Ｎ）から受ける引力の合計値、Ｖgは着目色点ｇの速度ベクトル、−ｋvＶgは速度に応じた抵抗力、Ｘgは着目色点ｇの位置ベクトル、Ｘgnは隣接色点ｇｎの位置ベクトル、ｋp，ｋvは係数である。係数ｋp，ｋvは予め一定の値に設定される。なお、文中では、ベクトルを示す矢印は省略される。
【００３７】
この力学モデルは、バネで互いに結ばれた質点の減衰振動モデルである。すなわち、着目色点ｇに係る仮想合力Ｆpgは、着目色点ｇと隣接色点ｇｎとの距離が大きいほど大きくなるバネ力Ｆgと、着目色点ｇの速度が大きいほど大きくなる抵抗力−ｋvＶgとの合計値である。この力学モデルでは、各色点について、位置ベクトルＸgと速度ベクトルＦgの初期値を設定した後に、微小時間経過後の位置ベクトルＸgと速度ベクトルＦgを順次算出してゆく。なお、複数の色点の速度ベクトルＶgの初期値は、例えば０に設定される。このような力学モデルを利用すれば、各色点が徐々に移動して、平滑な色点分布を得ることが可能である。
【００３８】
なお、各色点に係る力としては、バネ力Ｆgと抵抗力−ｋvＶg以外の力を用いても良い。例えば、本出願人により開示された特開２００６−１９７０８０号公報で説明されている他の種々の力をこの力学モデルで利用してもよい。また、力学モデルを適用して各色点を移動させる際に、特定の色点は、力学モデルによって移動しない拘束点として取り扱うことも可能である。
【００３９】
（３）スムージング処理：
図８は、図２のステップＳ３００で行うスムージング処理（平滑化／最適化処理）の典型的な処理手順を示すフローチャートである。ステップＴ９０では、最適化条件設定モジュール１３３が、スムージング処理におけるＬＵＴ作成条件を受付ける。このＬＵＴ作成条件は最適化処理の指針であり、後述のステップＴ１３０の最適化処理において用いる目的関数の作成条件に相当する。ＬＵＴ作成条件については、後述することにする。次に、ステップＴ１００では、スムージング処理初期値設定モジュール１２０（図１）が、スムージング処理の対象とする複数の色点を初期設定する。
【００４０】
図９は、図８のステップＴ１００の詳細手順を示すフローチャートである。ステップＴ１０２では、インク量の初期入力値（図３（Ｂ），図５（Ｂ））から、スムージング処理の対象となる各色点の仮インク量が決定される。例えば、３Ｄ−ＬＵＴ用のスムージング処理では、次の（２）式、（３）式に従って、各入力格子点に対する仮インク量Ｉ(R,G,B)が決定される。
【数２】
2011223298.tif 000004
【数３】
2011223298.tif 000005
ここで、Ｉ(R,G,B)は、入力格子点のＲＧＢ値に対するインクセット全体のインク量（図３の例では８種類のインクのインク量）を表している。ＲＧＢ値が０または２５５を取る入力格子点に対するインク量は、図２のステップＳ２００においてユーザーによって予め入力された値である。前記（２）式および（３）式によれば、任意のＲＧＢ値における仮インク量Ｉ(R,G,B)を求めることが可能である。
【００４１】
４Ｄ−ＬＵＴ用のスムージング処理では、次の（４）式、（５）式に従って、各入力格子点に対する仮インク量Ｉ(C,M,Y,K)が決定される。
【数４】
2011223298.tif 000006
【数５】
2011223298.tif 000007
【００４２】
なお、（４）式からも理解できるように、４Ｄ−ＬＵＴ用のインク量の初期入力値は１６個存在するので、初期入力値の設定が煩雑である。そこで、例えば、インク量の初期入力値を設定する入力格子点を、Ｋ＝０の８個の頂点、すなわち、（Ｃ，Ｍ，Ｙ，Ｋ）＝（０，０，０，０），（０，０，２５５，０），（０，２５５，０，０），（０，２５５，２５５，０），（２５５，０，０，０），（２５５，０，２５５，０），（２５５，２５５，０，０），（２５５，２５５，２５５，０）の８個の頂点と、Ｋ＝２５５の１個の頂点、例えば、（Ｃ，Ｍ，Ｙ，Ｋ）＝（０，０，０，２５５）の頂点のみとし、Ｋ＝２５５の色点のインク量を次の（６）式または（７）式で決定してもよい。
【数６】
2011223298.tif 000008
【数７】
2011223298.tif 000009
ここで、Ｉ(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を各行要素として有する行列（ベクトル）を意味することとする。
【００４３】
図９のステップＴ１０４では、フォワードモデルコンバーター３００を用いて、仮インク量に対応する色彩値Ｌ*ａ*ｂ*を求める。この演算は、以下の（８）式または（９）式で表すことができる。
【数８】
2011223298.tif 000010
【数９】
2011223298.tif 000011
ここで、Ｌ* (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に近いＬ*ａ*ｂ*値を再現するインク量Ｉjの中で、複数のパラメータがなるべく小さいインク量が最適なインク量として決定される。また、最適なインク量の探索は、ステップＴ１００で設定された各入力格子点の初期インク値から開始され、移動後の格子点について最適なインク量Ｉjを決定すると、次に力学モデルによって格子点を移動したときはこのインク量Ｉjから最適なインク量を探索する。従って、探索で得られるインク量は、この初期インク量を修正した値となる。なお、ステップＴ１３０の詳細手順や目的関数Ｅの内容については後述する。
【００４９】
図８のステップＴ１４０では、ステップＴ１３０で探索されたインク量Ｉjに対応するＬ*ａ*ｂ*値が、フォワードモデルコンバーター３００で再計算される（図１０（Ｄ）参照）。ステップＴ１３０で探索されたインク量Ｉjは目的関数Ｅを最小とするインク量なので、そのインク量Ｉjで再現されるＬ*ａ*ｂ*値は、最適化処理のターゲット値ＬＡＢtから多少ずれている。インク量Ｉｊに対応するＬ*ａ*ｂ*値を再計算することにより、探索されたインク量Ｉjの再現するＬ*ａ*ｂ*値が正確に算出される。こうして再計算されたＬ*ａ*ｂ*値が、各色点の移動後の座標値として使用される。
【００５０】
ステップＴ１５０では、各色点の座標値の移動量の平均値（ΔＬａｂ） aveが、予め設定された閾値ε以下であるか否かが判定される。移動量の平均値（ΔＬａｂ） aveが閾値εよりも大きい場合には、ステップＴ１２０に戻りステップＴ１２０〜Ｔ１５０のスムージング処理が継続される。一方、移動量の平均値（ΔＬａｂ） aveが閾値ε以下の場合には、色点の分布が十分に平滑になっているので、スムージング処理が終了する。なお、閾値εは、予め適切な値が実験的に決定される。
【００５１】
このように、本実施例の典型的なスムージング処理（平滑化／最適化処理）では、力学モデルによって各色点を微小時間間隔毎に移動させつつ、移動後の色点に対応する最適なインク量を最適化手法で探索する。そして、色点の移動量が十分に小さくなるまでそれらの処理が継続される。この結果、図３（Ｃ）に示したように、スムージング処理によって、平滑な色点分布を得ることが可能である。
【００５２】
（４）ＬＵＴ作成条件の設定：
図１１は、ＬＵＴ作成条件の設定処理（ステップＴ９０）の流れを示すフローチャートである。同図に示す処理は最適化条件設定モジュール１３３により実行される。ステップＴ９１では、最適化条件設定モジュール１３３がＵＩ画面をディスプレイに表示させ、印刷に使用する印刷用紙の種類と、印刷物を観察する観察光源の指定と、を受付ける。
【００５３】
図１２は、ステップＴ９１にて表示されるＵＩ画面を示している。同図に示すＵＩ画面では、印刷用紙の種類と、印刷物を観察する観察光源をと指定することができる。このＬＵＴ作成条件に基づいて作成されるベースＬＵＴは、ここで指定された印刷用紙の種類と観察光源とに最適化されることになる。また、同図に示すＵＩ画面では、複数の画質要素の優先順位と、総合的な最適化と最も重要な画質要素に関する最適化とのいずれを優先するかの設定、とが可能になっている。同図に示すＵＩ画面では、複数の画質要素として、粒状性、色恒常性、階調性、色再現ガマット、の優先順位を設定できるようになっている。なお、本実施例では、ここで設定する優先順位が、各画質要素の重要度を示すことになる。
【００５４】
ユーザーがマウスやキーボードを操作して、印刷用紙の種類と観察光源とを選択し、画質要素の優先順位の設定と、総合的な最適化と最も重要な画質要素に関する最適化とのいずれを優先するかの設定とを行い、決定ボタンをクリックすると、最適化条件設定モジュール１３３がこれを検知し、これらの設定を受付ける（Ｔ９２）。なお、以下では、印刷用紙として光沢紙が選択され、観察光源としてＤ６５が選択され、粒状性、色恒常性、階調性、色再現ガマットの順に優先するように設定されたものとして説明する。
【００５５】
ステップＴ９３では、最適化条件設定モジュール１３３が、ユーザーが設定した印刷用紙の種類と観察光源とに基づいて、下記の式（１０）における画質評価関数を設定し、画質要素の優先順位と総合的な最適化と最も重要な画質要素に関する最適化とのいずれを優先するかの設定とに基づいて、下記の式（１０）における重み係数ｗgr、ｗCII，ｗS，ｗgaを設定する。
【数１０】
2011223298.tif 000012
【００５６】
上記の（１０）式においては、Ｅは目的関数を示し、目的関数Ｅが小さければ小さいほど、ステップＳ２００でユーザーが設定したＬＵＴ作成条件の下で、総合的な印刷パフォーマンスが高くなるという性質を有している。（１０）式において、ψは、設定されたインクセット（本実施形態ではＣＭＹＫｌｃｌｍｌｋｌｌｋ）の各インクのインク量の組み合わせを意味するインク量セット（ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlm，ｄlk，ｄllk）を表している。上記の（１０）式の第１項は印刷物の粒状性の性能を要求する項である。第２項は印刷色の光源変動に対する色恒常性の性能を要求する項である。第３項は印刷物の階調性の性能を要求する項である。第４項は色再現ガマットの性能を要求する項である。
【００５７】
いずれの項も同一の大きさで正規化されたスカラー量であり、値が小さいほど印刷パフォーマンスが高い。また、各画質要素に対応する第１項〜第４項を個別の重み係数ｗgr、ｗCII，ｗS，ｗgaによって重みを調整しつつ線形結合することにより、総合的な印刷パフォーマンスが評価可能な目的関数Ｅを定義している。すなわち、重み係数ｗgr、ｗCII，ｗS，ｗgaは、どの画質要素を重視するかを調整する値としての意味を有している。
【００５８】
また、上記（１０）式の重み係数ｗgr、ｗCII，ｗS，ｗgaは、下記式（１１）で表される。
【数１１】
2011223298.tif 000013
【００５９】
上記（１１）式においてｎは最適化処理の進行度合いを示すカウンターであり、本実施例では、最適化処理を行った回数とする。また、１つの格子点に対応するインク量セットを最適化する最適化処理を実行する回数は、予め規定されているものとする（後述の最適化処理では１０回としてある）。ｗ gr0、ｗ CII0，ｗ S0，ｗ ga0は初期値であり、最適化処理の初期に利用される重み係数である。ａgr、ａCII，ａS，ａgaは重み補正係数であり、最適化処理の進行度合いに応じて重み係数を徐々に変化させるための係数である。例えば、重み係数は、上記（１１）式に示すように、最適化処理を繰り返すたびに、直前の最適化処理において利用された重み係数に対して、重み補正係数ａgr、ａCII，ａS，ａgaを乗じた値となっていく。すなわち、重み係数は、最適化処理が進むにつれて段階的に変更されていく。
【００６０】
また、重み係数ｗgr，ｗCII，ｗS，ｗgaは、最適化処理を実行した回数が規定回数に近づくにつれて、目標重み係数ｗgr_t，ｗCII_t，ｗS_t，ｗga_tに段階的に収束していく。すなわち、最適化処理の進行度合いに応じて、重み係数ｗgr，ｗCII，ｗS，ｗgaは、初期値ｗ gr0、ｗ CII0，ｗ S0，ｗ ga0における重みバランスから、目標重み係数ｗgr_t，ｗCII_t，ｗS_t，ｗga_tの重みバランスへと段階的に収束するようになっている。本実施例では、上記（１１）式に示すように、最適化回数が規定回数に到達したときに目標重み係数に到達するようにしてあるが、むろん、規定回数に到達する前に目標重み係数ｗgr_t，ｗCII_t，ｗS_t，ｗga_tに収束するようにして目標重み係数に到達した後の最適化処理では重み係数を変化させないようにしてもよいし、毎回で所定回数の最適化処理が終了するごとに重み係数ｗgr，ｗCII，ｗS，ｗgaを重み補正係数ａgr、ａCII，ａS，ａgaによって変化させるようにすることもできる。
また、本実施例では、重み補正係数を幾何級数的に変化させているが、むろん、その他の手法で重み補正係数を変化させてもよく、例えば、重み係数を算術級数的に変化させてもよい。この場合は、１回の最適化処理が終了するごとに初期値に対して所定の重み補正量を加算していくことになる。
【００６１】
また、各画質要素に対応する初期値ｗ gr0、ｗ CII0，ｗ S0，ｗ ga0と重み補正係数ａgr、ａCII，ａS，ａgaは、総合的な最適化と最も重要な画質要素に関する最適化とのいずれを優先するかの設定に応じて、大小関係が設定される。図１３は、総合的な最適解へ最適化する場合の重みバランスを説明する図であり、図１４は、最も重要な画質要素に関する最適化を優先する場合の重みバランスを説明する図である。
【００６２】
まず、総合的な最適解へ最適化する場合は、図１３に示すように、優先順位の高い画質要素にかかる初期値ｗ gr0、ｗ CII0，ｗ S0，ｗ ga0が、各重み係数に対して最終的に設定したい目標重み係数ｗgr_t，ｗCII_t，ｗS_t，ｗga_tにおける重みバランスに比べて、相対的に低い重みを有するような重みバランスとなるように設定される。
【００６３】
例えば、目標重み係数の重みバランスが、ｗgr_t：ｗCII_t：ｗS_t：ｗga_t＝５：４：３：２であるときには、初期値の重みバランスをｗ gr0：ｗ CII0：ｗ S0：ｗ ga0＝１：１：１：１などのように設定する。この例では、優先順位のもっとも高い画質要素である粒状性の重み係数ｗｇｒは、最終的に設定したい重みバランスでは他の重み係数ｗＣＩＩ，ｗＳ，ｗｇａに対してそれぞれ５／４倍、５／３倍、５／２倍の重みを有するのに対し、初期値では他の重み係数に対して１倍と平均的な重みバランスになっている。すなわち、粒状性の重み係数ｗｇｒは、初期値では目標重み係数よりも相対的に低い重みバランスとなっている。
【００６４】
このような重みバランスで初期値を設定したときには、重み補正係数は、画質要素の優先順位に従った大小関係とされる。すなわち、上述のように、最終的に設定したい重みバランスがｗgr＞ｗCII＞ｗS＞ｗgaのとき、重み補正係数の大小関係は、画質要素の優先順位と同様にａgr＞ａCII＞ａS＞ａgaとされる。このように重み係数の初期値の重みバランスと重み補正係数の大小関係とを設定することにより、初期の最適化処理においては、重視したい画質要素を十分に最適化できない状態を許容しつつ、その後、段階的に目標重み係数の重みバランスとなるように重み係数が調整されていく。よって、目標重み係数を設定された目的関数において総合的に最適なインク量セットへと徐々に導かれるため、ローカルミニマムに陥りにくく、所望の画質要素が良好なインク量セットを探索することが可能となる。
【００６５】
一方、最も重要な画質要素に関する最適化を優先する場合は、図１４に示すように、優先順位の高い重み係数の初期値ｗ gr0、ｗ CII0，ｗ S0，ｗ ga0が、各重み係数に対して最終的に設定したい目標重み係数ｗgr_t，ｗCII_t，ｗS_t，ｗga_tにおける重みバランスに比べて、相対的に高い重みを有するような重みバランスとなるように設定される。
【００６６】
例えば、目標重み係数の重みバランスが、ｗgr_t：ｗCII_t：ｗS_t：ｗga_t＝５：４：３：２であるときには、初期値の重みバランスをｗ gr0：ｗ CII0：ｗ S0：ｗ ga0＝９：７：５：３などのように設定する。この例では、優先順位のもっとも高い画質要素である粒状性の重み係数ｗｇｒは、最終的に設定したい重みバランスでは他の重み係数ｗＣＩＩ，ｗＳ，ｗｇａに対してそれぞれ５／４倍、５／３倍、５／２倍の重みを有するが、初期値では他の重み係数に対してそれぞれ９／７倍、９／５倍、９／３倍の重みとなっている。すなわち、粒状性の重み係数ｗｇｒは、初期値では目標重み係数よりも相対的に高い重みバランスとなっている。
【００６７】
このような重みバランスで初期値を設定したときには、重み補正係数は、画質要素の優先順位とは逆の大小関係とされる。すなわち、上述のように、最終的に設定したい重みバランスがｗgr＞ｗCII＞ｗS＞ｗgaのときは、重み補正係数の大小関係は、画質要素の優先順位とは逆のａgr＜ａCII＜ａS＜ａgaとされる。このように重み係数の初期値の重みバランスと重み補正係数の大小関係とを設定することにより、初期の最適化処理においては、重視したい画質要素について優先的に最適化される状態とし、その後、段階的に全画質要素のトータルで本来設定したい重みバランスとなるように重み係数が調整されていく。よって、仮に初期の最適化処理においてローカルミニマムに陥ると、最適化処理の終了時においてもそのローカルミニマムに落ち着いている可能性が高く、優先順位の最も高い画質要素について特に良好な画質を再現するＬＵＴを作成することができる。
【００６８】
また、目的関数Ｅの粒状性指数ＧＲ，色恒常性指数ＣＩＩ、階調性指数Ｓ、色再現ガマットの評価指数ΔＥは、それぞれ下記のようにして算出することができる。
【００６９】
粒状性指数ＧＲは、各種の粒状性予測モデルを用いて算出可能であり、例えば以下の（１２）式で算出することができる。
【数１２】
2011223298.tif 000014
ここで、ａＬは明度補正係数、ＷＳ（ｕ）はカラーパッチの印刷に利用されるハーフトーンデータが示す画像のウイナースペクトラム、ＶＴＦ（ｕ）は視覚の空間周波数特性、ｕは空間周波数である。ハーフトーンデータは、カラーパッチのインク量Ｉjからハーフトーン処理（本実施例で作成するＬＵＴが搭載されるプリンターが実行するハーフトーン処理と同一のものとする）によって決定される。前記（１２）式は一次元で表現しているが、空間周波数の関数として二次元画像の空間周波数を算出することは容易である。粒状性指数ＧＩの計算方法としては、例えば、本出願人により開示された特表２００７−５１１１６１号公報に記載された方法や、Makoto Fujino， Ｉmage Quality Evaluation of Inkjet Prints， Japan Hardcopy '99， p.291-294に記載された方法を利用することができる。
【００７０】
色非恒常性指数ＣＩＩは、例えば以下の（１３）式で与えられる。
【数１３】
2011223298.tif 000015
ここで、ΔＬ*は２つの異なる観察条件下（異なる光源下）におけるカラーパッチの明度差、ΔＣ*abは彩度差、ΔＨ*abは色相差を示す。色非恒常性指数ＣＩＩの計算時には、２つの異なる観察条件下でのＬ*ａ*ｂ*値は、色順応変換（ＣＡＴ）を用いて標準観察条件（例えば標準の光Ｄ６５の観察下）に変換される。ＣＩＩについては、Billmeyer and Saltzman's Principles of Color Technology， 3rd edition， John Wiley & Sons， Ｉnc， 2000， p.129， pp. 213-215を参照。
【００７１】
階調性指数Ｓは、例えば以下の（１４）式で与えられる。
【数１４】
2011223298.tif 000016
ここで、Ｌｐは注目する格子点の位置ベクトル、Ｌ a1〜Ｌ a6は当該格子点に隣接する６個の格子点の位置ベクトルを示す。ここで、隣接する格子点を結ぶ線（ベクトルＬ a1〜ベクトルＬp〜ベクトルＬ a2が示す格子点を通る線等）が直線に近く、また格子点が均等に配置されるほどＣＩＥＬＡＢ色空間における格子点の配置が平滑化される傾向にある。すなわち階調性指数Ｓが小さくなればなるほど、平滑程度が高く、良好な階調性が期待できる。
【００７２】
色再現ガマットの評価指数ΔＥは、例えば、ＬＵＴの適用対象であるプリンターの色再現ガマットデータＧＤから色再現ガマットを取得し、当該色再現ガマットの外面上や稜線上や頂点上のＬ*ａ*ｂ*値と、分光プリンティングモデルコンバーター３１０と色算出部３２０とにより算出した一部の格子点のＬ*ａ*ｂ*値との色差ΔＥとして算出される。なお、色再現ガマットデータＧＤは、ハードディスク装置などに予め格納されているものとする。色差ΔＥの算出の対象は、格子点のうちインク量空間の外縁に存在するものとされ、内側の格子点については色差ΔＥ＝０とする。インク量空間の外縁に存在する格子点を分光プリンティングモデルコンバーター３１０と色算出部３２０とによって変換したＬ*ａ*ｂ*値も同様にＣＩＥＬＡＢ空間において外縁に存在させるためである。これにより、格子点のうちインク量空間の外縁に存在するものを分光プリンティングモデルコンバーターと色算出部３２０とによって順次変換したＬ*ａ*ｂ*値について、上記プリンターが再現可能な色再現ガマットの外面上や稜線上や頂点上までの色ずれを定量化することができる。
【００７３】
以上で説明したように、ユーザーからＬＵＴ作成条件の設定を受付けると、目的関数における画質評価指数に対する重み係数ｗgr、ｗCII，ｗS，ｗgaの初期値と重み補正係数とが設定できるようになり、目的関数Ｅが確定されるため、以下に説明する最適化処理の準備が整ったことになる。
【００７４】
（５）最適化処理：
図１５は、最適化処理の流れを示すフローチャートである。同図に示す処理は、図８のステップＴ１３０に相当する処理であり、力学モデルにしたがってターゲットＬＡＢｔが変更される度にインク量最適化モジュール１３４によって実行される。処理が開始されると、ステップＴ１３１では、ステップＴ９０において設定されたＬＵＴ作成条件、すなわち目的関数Ｅおよび重み係数の初期値ｗ gr0、ｗ CII0，ｗ S0，ｗ ga0と重み補正係数ａgr、ａCII，ａS，ａgaを取得する。
【００７５】
そして、格子点を１つ選択し（Ｔ１３２）、重み係数を管理するためのカウンターｎを初期化し（Ｔ１３３）、カウンターｎに基づいて目的関数Ｅの重み係数を設定する（Ｔ１３４）。すなわち、目的関数の重み係数は、上述した式（１１）に示すように、最適化処理の初期においては初期値ｗ gr0、ｗ CII0，ｗ S0，ｗ ga0となり、最適化が進むにつれて重み補正係数ａgr、ａCII，ａS，ａgaによって段階的に変更されていく。本実施例では、カウンターｎが上述した式（１１）のｎに相当するものとし、各格子点の最適化処理を１回実行するたびにカウンターｎを増加し（Ｔ１３６）、このカウンターｎを重み補正係数の乗数として利用している。
【００７６】
ステップＴ１３５では、ステップＴ１３２で選択された格子点について初期のインク量セットを最適化していく。具体的には、インク量空間における初期のインク量セットの位置から局所的にインク量セットを移動させ、その際に目的関数Ｅを極小化させるインク量セットを各格子点について算出し、ステップＴ１３２で選択された格子点について目的関数Ｅを極小化させるインク量セットを探索する。このとき、最適化のアルゴリズム次第ではあるが、例えば、インク量空間において、その時点で目的関数を最小化するインク量セットと隣接するインク量セットの中から最小化するインク量セットを逐次選択し、さらに目的関数を低下させるインク量セットが無くなるまでインク量セットの探索を継続する手法を用いる場合には、本実施例のように重みバランスを調整することにより、インク量空間における最適化の方向性を、所定のＬＵＴ作成条件に合致するようにコントロールすることができる。
【００７７】
このようにしてステップＴ１３２で選択された格子点について１回の最適化処理が終了すると（Ｔ１３５）、上述したようにカウンターｎをインクリメントし（Ｔ１３６）最適化処理の回数が所定の回数ｎｍａｘに到達したか否かを判断し（Ｔ１３７）、到達していれば当該格子点についての最適化処理を終了し（Ｔ１３７：Ｙｅｓ）、到達していなければ重み係数を更新してさらに最適なインク量セットを探索する（Ｔ１３７：Ｎｏ）。この繰り返し回数ｎｍａｘは、例えば、１０回等に設定されている。
【００７８】
各回の最適化処理では、その回において設定された重み係数を適用した目的関数を最小化するインク量セットを探索し、その次の回の最適化処理では、直前の最適化処理で目的関数を最小化したインク量セットを基点として、さらに目的関数を最小化するインク量セットを探索することになる。本実施例では目的関数が徐々に変化していくため、各回の目的関数における最適解も徐々に変化していく。
【００７９】
このとき、図１３に示す例のように、各画質要素に対する重みバランスが徐々に強調されるように目的関数を変化させることにより、最適化処理が繰り返されるにつれて最終的な目的関数の最適解へと徐々に誘導することができる。このようにして最適化処理が所定の回数だけ繰り返し実行されると、最終的には各格子点についての目的関数Ｅが極めて小さくなる（総合的な印刷パフォーマンスが高い）インク量セットに格子点を最適化することができる。むろん、目的関数Ｅの値が所定の閾値を下回ることをもって格子点の最適化を完了させてもよい。
【００８０】
また、逆に、図１４に示す例のように、各画質要素に対する重みバランスが徐々に緩和されるように目的関数を変化させることにより、最終的な目的関数におけるローカルミニマムに相当する解を初期の最適化処理の最適解とし、初期の最適化処理により得られる最適解を最終的な最適化処理におけるローカルミニマムと同一の解へと誘導することができる。このようにして最適化処理を所定の回数だけ繰り返し実行されると、最終的な目的関数による最適化処理においても初期に陥ったローカルミニマムから脱出しにくくなる。よって、最も重要視する画質要素を優先的に最適化しつつ、この状態で他の画質要素を最適化したインク量セットに格子点を最適化することができる。
【００８１】
ステップＴ１３８では、全格子点について最適化処理を完了したか否かを判断する。ステップＴ１２０の力学モデルにしたがって決定されたＬＡＢｔ点の全てについて最適化処理を行って最適なインク量セットの探索を完了していれば最適化処理を終了し（Ｔ１３８：Ｙｅｓ）、完了していなければステップＴ１３２に戻って未処理の格子点について最適化処理を行う（Ｔ１３８：Ｎｏ）。以上、本実施例においては、インク量最適化モジュール１３４が実行する最適化処理が最適化工程を構成し、ステップＴ１３１，Ｔ１３４やＬＵＴ作成条件の設定処理において最適化条件設定モジュール１３３がＬＵＴ作成条件を受付ける工程が重み管理工程を構成する。
【００８２】
以上のようにして各格子点が最適化されると、スムージング処理のステップＴ１５０おいて、最適化された格子点のインク量セットに対応した色（Ｌ*ａ*ｂ*）を分光プリンティングモデルコンバーター３１０および色算出部３２０によって算出し、各色点の移動量の平均値（ΔＬａｂ）を評価することになる。
【００８３】
以上のように、各画質評価指数の重視度合いを重み係数ｗgr、ｗCII，ｗS，ｗgaによって設定された総合的な印刷パフォーマンスの指標となる目的関数Ｅを使用して各格子点を最適化することにより、複雑な分版規則を規定することなく好適なインク量セットを探し出すことができる。さらに、ユーザーの設定した画質評価指数の優先順位に従って、目的関数における各画質評価指数の重み係数を最適化処理が進むにつれて段階的に変化させていくことにより、より尊重したい画質要素のローカルミニマムへ陥りにくくしたり、逆に尊重したい画質要素のローカルミニマムへ陥りやすくしたり、と最適化処理をコントロールできるようになる。
（６）変形例：
なお、この発明は前記の実施例や実施形態に限られるものではなく、その要旨を逸脱しない範囲において種々の態様において実施することが可能であり、例えば次のような変形も可能である。
【００８４】
（６−１）変形例１：
本明細書において「インク」とは、インクジェットプリンタやオフセット印刷等に用いられる液体状インクに限らず、レーザプリンタに用いられるトナーも含む広い意味で使用されている。このような「インク」の広い意味を有する他の用語としては、「色材」や「着色材」、「着色剤」を用いることも可能である。
【００８５】
（６−２）変形例２：
前記実施例では、ルックアップテーブルのような色変換プロファイルを作成する方法および装置に関して説明したが、本発明は、こうして得られた色変換プロファイルを印刷装置に組み込む組み込み部を備える印刷装置製造システムにも適用可能である。色変換プロファイルを作成する色変換プロファイル作成装置は、この印刷装置製造システムに含まれるものとしてもよく、他のシステムや装置に含まれるものとしてもよい。なお、この製造システムの組み込み部は、例えば、プリンタドライバのインストーラ（インストールプログラム）として実現することができる。
【００８６】
（６−３）変形例３：
上記実施例では、ルックアップテーブルのような色変換プロファイルを作成する方法および装置に関して説明したが、本発明は、こうして得られた色変換プロファイルそのものであったり、色変換プロファイルを記憶する記憶部を備えており、入力された印刷データを該色変換プロファイルに基づいて変換して印刷を実行する印刷装置によっても実現することができる。
【符号の説明】
【００８７】
１００…ベースＬＵＴ作成モジュール、１２０…スムージング処理初期値設定モジュール、１３０…スムージング処理モジュール、１３２…色点移動モジュール、１３３…最適化条件設定モジュール、１３４…インク量最適化モジュール、１３６…画質評価指数算出モジュール、１４０…テーブル作成モジュール、２００…色補正ＬＵＴ作成モジュール、３００…フォワードモデルコンバーター、３１０…分光プリンティングモデルコンバーター、３２０…色算出部、４００…ＬＵＴ格納部、４１０…インバースモデル初期ＬＵＴ、５１０…ベース３Ｄ−ＬＵＴ、５２０…ベース４Ｄ−ＬＵＴ、６１０…色補正３Ｄ−ＬＵＴ、６２０…色補正４Ｄ−ＬＵＴ、５１１…逆変換ＬＵＴ、５２１…逆変換ＬＵＴ
【図１】
2011223298.tif 000017
【図２】
2011223298.tif 000018
【図３】
2011223298.tif 000019
【図４】
2011223298.tif 000020
【図５】
2011223298.tif 000021
【図６】
2011223298.tif 000022
【図７】
2011223298.tif 000023
【図８】
2011223298.tif 000024
【図９】
2011223298.tif 000025
【図１０】
2011223298.tif 000026
【図１１】
2011223298.tif 000027
【図１２】
2011223298.tif 000028
【図１３】
2011223298.tif 000029
【図１４】
2011223298.tif 000030
【図１５】
2011223298.tif 000031