(57)【要約】【課題】色変換パラメータによって実現される階調性をより良好にする。【解決手段】本発明の色変換パラメータの作成装置は、入力画像信号値を複数種類の実色材値に変換する色変換処理で用いられる色変換パラメータの作成装置であって、前記実色材値を変換して得られる仮想色材値を用いて前記色変換パラメータの評価値を算出する算出手段と、前記評価値に基づいて前記色変換パラメータを修正する修正手段と、を有し、前記仮想色材値は、前記実色材の分光反射率をｎ個（ｎは３以上の整数）のブロックに分割した場合に、分割された前記各ブロックにおける分光反射率の対数に対応するｎ種類の仮想的な色材の出力値であることを特徴とする。【選択図】図１１
2018191132.tif 000002
【特許請求の範囲】
【請求項１】
入力画像信号値を複数種類の実色材値に変換する色変換処理で用いられる色変換パラメータの作成装置であって、
前記実色材値を変換して得られる仮想色材値を用いて前記色変換パラメータの評価値を算出する算出手段と、
前記評価値に基づいて前記色変換パラメータを修正する修正手段と、を有し、
前記仮想色材値は、前記実色材の分光反射率をｎ個（ｎは３以上の整数）のブロックに分割した場合に、分割された前記各ブロックにおける分光反射率の対数に対応するｎ種類の仮想的な色材の出力値である
ことを特徴とする色変換パラメータの作成装置。
【請求項２】
前記評価値は、前記入力画像信号値の色空間における前記仮想色材値の変化の滑らかさであることを特徴とする請求項１に記載の色変換パラメータの作成装置。
【請求項３】
前記評価値は、前記仮想色材値の色空間における前記実色材値の変化の滑らかさであることを特徴とする請求項１に記載の色変換パラメータの作成装置。
【請求項４】
前記変化の滑らかさは、評価対象位置における前記仮想色材値または前記実色材値の変化の線形性、あるいは、前記評価対象位置の配置の線形性であることを特徴とする請求項２または３に記載の色変換パラメータの作成装置。
【請求項５】
前記変化の滑らかさは、評価対象位置における前記仮想色材値または前記実色材値の曲率であることを特徴とする請求項２または３に記載の色変換パラメータの作成装置。
【請求項６】
前記変化の滑らかさは、評価対象位置と当該評価対象位置の近傍位置とを結ぶ変位ベクトルの和であることを特徴とする請求項２または３に記載の色変換パラメータの作成装置。
【請求項７】
前記色変換パラメータは、前記色変換処理で用いられる数式の係数であり、
前記修正手段は、前記評価値に基づいて前記係数を修正する
ことを特徴とする請求項１から６のいずれか１項に記載の色変換パラメータの作成装置。
【請求項８】
前記色変換パラメータは、前記色変換処理で用いられる３次元色分解ＬＵＴであり、
前記修正手段は、前記評価値に基づいて前記３次元色分解ＬＵＴにおける格子点に保持される前記実色材値を修正する
ことを特徴とする請求項１から７のいずれか１項に記載の色変換パラメータの作成装置。
【請求項９】
前記入力画像信号値を前記３次元色分解ＬＵＴの入力値として用いられる第２の入力画像信号値に非線形変換する非線形変換手段と、
該非線形変換で用いられる非線形変換パラメータを修正する第２の修正手段と、をさらに有し、
前記算出手段は、前記第２の入力画像信号値を可変させた場合の複数の前記実色材値を変換して得られる複数の前記仮想色材値を用いて複数の前記評価値を算出し、
前記第２の修正手段は、前記入力画像信号値を、前記複数の評価値のうち最良の評価値に対応する前記第２の入力画像信号値に変換するように、前記非線形変換パラメータを修正する
ことを特徴とする請求項８に記載の色変換パラメータの作成装置。
【請求項１０】
前記非線形変換パラメータは、前記非線形変換で用いられる１次元非線形変換ＬＵＴであることを特徴とする請求項９に記載の色変換パラメータの作成装置。
【請求項１１】
前記非線形変換パラメータは、前記複数種類の実色材ごとに存在することを特徴とする請求項９または１０に記載の色変換パラメータの作成装置。
【請求項１２】
前記３次元色分解ＬＵＴを用いて変換された前記実色材値を第２の実色材値に非線形変換する非線形変換手段と、
該非線形変換で用いられる第２の非線形変換パラメータを修正する第３の修正手段と、をさらに有し、
前記算出手段は、前記実色材値を可変させた場合の複数の前記第２の実色材値を変換して得られる複数の前記仮想色材値を用いて複数の前記評価値を算出し、
前記第３の修正手段は、前記実色材値を、前記複数の評価値のうち最良の評価値に対応する前記第２の実色材値に変換するように、前記第２の非線形変換パラメータを修正する
ことを特徴とする請求項８から１１のいずれか１項に記載の色変換パラメータの作成装置。
【請求項１３】
前記第２の非線形変換パラメータは、ＯＰＧ処理で用いられる１次元非線形変換ＬＵＴであることを特徴とする請求項１２に記載の色変換パラメータの作成装置。
【請求項１４】
前記第２の非線形変換パラメータは、前記複数種類の実色材ごとに存在することを特徴とする請求項１２または１３に記載の色変換パラメータの作成装置。
【請求項１５】
前記評価値に基づいて前記３次元色分解ＬＵＴにおける前記格子点の数を決定する決定手段をさらに有することを特徴とする請求項８から１４のいずれか１項に記載の色変換パラメータの作成装置。
【請求項１６】
前記算出手段は、前記３次元色分解ＬＵＴにおける前記格子点を減らした場合の前記実色材値を変換して得られる前記仮想色材値を用いて前記評価値を算出し、
前記決定手段は、前記評価値が所定の閾値以上となった場合に、前記３次元色分解ＬＵＴにおける前記格子点の数を、該評価値に対応する前記格子点の数に決定する
ことを特徴とする請求項１５に記載の色変換パラメータの作成装置。
【請求項１７】
前記決定手段は、前記複数種類の実色材ごとに前記３次元色分解ＬＵＴにおける前記格子点の数を決定することを特徴とする請求項１５または１６に記載の色変換パラメータの作成装置。
【請求項１８】
所定の色変換パラメータを用いて入力画像信号を複数種類の実色材の量を示す実色材量信号に変換する色変換手段と、
前記実色材量信号を仮想色材量信号に変換する変換手段と、
前記入力画像信号が表す画像におけるエッジ度と、前記仮想色材量信号が表す画像におけるエッジ度との差に基づいて擬似輪郭画像を生成する生成手段と、
前記擬似輪郭画像を出力する出力手段と、を有し、
前記擬似輪郭画像は、前記エッジ度の差が大きい領域が識別可能な態様で表された画像であり、
前記仮想色材量信号は、前記実色材の分光反射率をｎ個（ｎは３以上の整数）のブロックに分割した場合に、分割された前記各ブロックにおける分光反射率の対数に対応するｎ種類の仮想的な色材の量を示す
ことを特徴とする画像処理装置。
【請求項１９】
入力画像信号値を複数種類の実色材値に変換する色変換処理で用いられる色変換パラメータの修正方法であって、
前記実色材値を変換して得られる仮想色材値を用いて前記色変換パラメータの評価値を算出する算出ステップと、
前記評価値に基づいて前記色変換パラメータを修正する修正ステップと、を有し、
前記仮想色材値は、前記実色材の分光反射率をｎ個（ｎは３以上の整数）のブロックに分割した場合に、分割された前記各ブロックにおける分光反射率の対数に対応するｎ種類の仮想的な色材の出力値である
ことを特徴とする色変換パラメータの修正方法。
【請求項２０】
コンピュータを、請求項１から請求項１７のいずれか１項に記載の色変換パラメータの作成装置の各手段として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、入力画像信号値を複数種類の実色材値に変換する色変換処理で用いられる色変換パラメータの作成装置、色変換パラメータの修正方法などに関する。
【背景技術】
【０００２】
インクジェット方式に代表されるプリンタは、入力画像信号（通常、ＲＧＢ３チャンネルの色信号）を、プリンタが具備する色材（例えば、ＣＭＹＫのインク）の出力量を示す色材量信号に変換することにより印刷データを生成する。この変換は色変換処理または色分解処理と呼ばれ、現在では、入力画像信号と色材量信号との対応を３次元の色分解ルックアップテーブル（以下「ルックアップテーブル」は「ＬＵＴ」と記す）を用いる方式が主流となっている。一般的な色分解ＬＵＴの作成方法によれば、主要な入力画像信号値に対応する色材量値のみがＬＵＴ設計者によって設定され、それ以外の入力画像信号値に対応する色材量値は、補間演算によって算出される。
【０００３】
特許文献１は、補間演算の精度を向上させるために、色分解ＬＵＴにおける格子点の配置を修正する方法を開示している。特許文献１における色分解ＬＵＴの修正方法によれば、色分解ＬＵＴにおける格子点位置がＬ＊ａ＊ｂ＊色空間における座標（Ｌａｂ値）に変換され、Ｌ＊ａ＊ｂ色空間においてこれらＬａｂ値が平滑化される。そして、平滑化されたＬａｂ値に基づいて、色分解ＬＵＴにおける格子点の配置が修正される。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００６−１９７０８０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、特許文献１の色分解ＬＵＴの修正方法では、Ｌａｂ値と該Ｌａｂ値に対応付けられる色材量値とは線形の関係になるとは限らないため、Ｌａｂ値が平滑化されたとしても、色分解ＬＵＴによって実現される階調性が良好とならない場合があった。印刷物にグラデーション領域が含まれている場合、入力画像信号値に対する色材量値の変化が滑らかでないと、このグラデーション領域において擬似輪郭などの不具合が発生してしまうおそれがある。
【０００６】
本発明は上記課題に鑑み、色変換パラメータによって実現される階調性をより良好にすることを目的とする。
【課題を解決するための手段】
【０００７】
本発明の色変換パラメータの作成装置は、入力画像信号値を複数種類の実色材値に変換する色変換処理で用いられる色変換パラメータの作成装置であって、前記実色材値を変換して得られる仮想色材値を用いて前記色変換パラメータの評価値を算出する算出手段と、前記評価値に基づいて前記色変換パラメータを修正する修正手段と、を有し、前記仮想色材値は、前記実色材の分光反射率をｎ個（ｎは３以上の整数）のブロックに分割した場合に、分割された前記各ブロックにおける分光反射率の対数に対応するｎ種類の仮想的な色材の出力値であることを特徴とする。
【発明の効果】
【０００８】
本発明によれば、色変換パラメータによって実現される階調性をより良好にすることができる、という効果を奏する。
【図面の簡単な説明】
【０００９】
【図１】実施形態１における印刷システムの全体構成例を示すブロック図である。
【図２】実施形態１における情報処理装置のハードウェア構成例を示す図である。
【図３】実施形態１における画像形成手順を示すフローチャートである。
【図４】実施形態１における分光反射率の一例を示すグラフである。
【図５】実施形態１におけるブロック濃度の一例を示すグラフである。
【図６】実施形態１において仮想色材量への換算処理を説明する図である。
【図７】実施形態１における３次元色分解ＬＵＴの一例を示す図である。
【図８】実施形態１において色変換特性の一例を示す図である。
【図９】実施形態１における色分解ＬＵＴの修正手順を示すフローチャートである。
【図１０】実施形態１における評価値の算出方法を説明する図である。
【図１１】実施形態１における実色材量が修正される様子を説明する図である。
【図１２】実施形態２における印刷システムの全体構成例を示すブロック図である。
【図１３】実施形態２における非線形変換ＬＵＴの例を示す模式図である。
【図１４】実施形態２における非線形変換ＬＵＴの修正手順を示すフローチャートである。
【図１５】実施形態３における印刷システムの全体構成例を示すブロック図である。
【図１６】実施形態３におけるＯＰＧＬＵＴの例を示す模式図である。
【図１７】実施形態３におけるＯＰＧＬＵＴが修正される様子を説明する図である。
【図１８】実施形態４における色分解ＬＵＴの修正手順を示すフローチャートである。
【図１９】実施形態５における印刷システムの全体構成例を示すブロック図である。
【図２０】実施形態５における擬似輪郭画像の生成手順を示すフローチャートである。
【図２１】実施形態５における擬似輪郭画像の一例を示す図である。
【発明を実施するための形態】
【００１０】
以下、本発明を実施するための形態について図面を参照して説明する。ただし、この実施形態に記載されている構成要素はあくまで例示であり、本発明の範囲をそれらに限定する趣旨のものではない。
【００１１】
［実施形態１］
（印刷システムの全体構成）
図１は、本実施形態における印刷システム１の全体構成例を示すブロック図である。本実施形態の印刷システム１は、画像処理装置１００と、画像形成装置２００と、ＬＵＴ修正装置３００とから構成される。画像処理装置１００と、画像形成装置２００と、ＬＵＴ修正装置３００とは、プリンタインタフェース（以下インタフェースは「Ｉ／Ｆ」とも記す）、回路またはネットワークＩ／Ｆを介して通信可能に接続されている。
【００１２】
本実施形態における画像処理装置１００は、一般的な情報処理装置（パーソナルコンピュータ）にインストールされたプリンタドライバによって実現することができる。あるいは、画像処理装置１００の別構成例では、画像形成装置２００が画像処理装置１００を内包してもよい。画像処理装置１００は、印刷対象の画像を表す画像データの入力を受け付け、カラーマッチング処理部１０１に画像データを送る。この画像データは、例えば８ビットのＲＧＢ画像信号に基づくカラー画像データである。
【００１３】
カラーマッチング処理部１０１は、入力された画像データにカラーマッチング処理を行い、ＲＧＢ画像の色を補正する。カラーマッチング処理により、異なる色再現特性を有するプリンタや記録媒体を用いた場合にも、統一的な色味を再現することができる。カラーマッチング処理が行われる際、カラーマッチングＬＵＴ格納部１０２に格納された３次元のカラーマッチングＬＵＴが参照される。カラーマッチングＬＵＴには、各８ビット（０〜２５５）のＲＧＢ信号値が１７×１７×１７に間引かれた格子点に保持されている。格子点間のＲＧＢ信号値は線形補間により算出される。
【００１４】
色分解処理部１０３は、カラーマッチング処理部１０１で補正された画像データの入力を受け付け、画像形成装置２００が具備する複数種類のインク量（インク値）に対応するインク値画像データを生成する。本実施形態では、画像形成装置２００は６色のインクを具備し、色分解処理部１０３によって、６プレーン８ビットのインク値画像データが生成される。本実施形態において、画像形成装置２００が具備する６色のインクはシアン（Ｃ）、マゼンタ（Ｍ）、イエロー（Ｙ）、ブラック（Ｋ）、淡シアン（Ｌｃ）、淡マゼンタ（Ｌｍ）であるが、上記以外のインクの組み合わせであってもよい。色分解処理部１０３は、色分解ＬＵＴ格納部１０４に格納されている色分解ＬＵＴを参照し、カラーマッチング処理後の画像データに対して色分解処理を行う。色分解ＬＵＴは、９×９×９の格子点を有する３次元ＬＵＴであり、色分解ＬＵＴの格子点には６色のインクの出力量を示すインク値（実色材値）が保持されている。色分解処理部１０３は、この色分解ＬＵＴを参照して、格子点間のインク値を線形補間により算出する。なお、本実施形態における色分解ＬＵＴの設計方法については後述する。
【００１５】
ハーフトーン処理部１０５は、色分解処理部１０３から出力された各色のインク値画像データを２値（または２値以上で入力階調数より少ない階調数の多値）に変換する量子化処理を行う。本実施形態において、ハーフトーン処理を行う手法として、公知のディザマトリクス法が選択される。本実施形態のハーフトーン処理部１０５は、ディザマトリクス格納部１０６に格納されたディザマトリクスを参照して、上記量子化処理を行う。ただし、ハーフトーン処理の方法はこれに限定されず、例えば、公知の誤差拡散法を用いてもよい。ハーフトーン処理部１０５によって生成された２値画像データは、出力Ｉ／Ｆまたは出力端子を介して画像形成装置２００に出力される。
【００１６】
画像形成装置２００は、カット用紙などの記録媒体２０７に対して、記録ヘッド２０５を相対的に縦横に移動することにより、画像処理装置１００で形成された２値画像データが表す画像を記録媒体２０７上に形成する。本実施形態では、画像形成装置２００が、記録ヘッド２０５を備えるインクジェット方式のプリンタである場合を例に説明する。記録ヘッド２０５は、複数の記録素子（ノズル）を有する。本実施形態の記録ヘッド２０５は、Ｃ、Ｍ、Ｙ、Ｋ、Ｌｃ、Ｌｍの６色のインクタンクを搭載する。ヘッド制御部２０４は、駆動部２０８の動作を制御する。駆動部２０８は、ヘッド制御部２０４の制御下で、記録ヘッド２０５を移動する。ヘッド制御部２０４は、搬送部２０６の動作も制御する。搬送部２０６は、ヘッド制御部２０４の制御下で、記録媒体を搬送する。なお、本実施形態では、記録ヘッド２０５が記録媒体２０７に対して複数回の走査を行うことにより画像を形成する。このような画像形成方法は、いわゆるマルチパス記録方式と呼ばれる。パス分解処理部２０１は、画像処理装置１００で生成された各色の２値画像データと、パスマスク格納部２０２とから取得したパスマスクとに基づき、各色の走査データを生成する。なお、走査データとは、複数に分割されたノズル群が各記録走査において記録するパターンを示すデータである。これらパターンは互いに補完関係にあり、全ノズル群のパターンを重ね合わせると全領域における画像形成が完成する。各記録走査が終了するたびに、記録媒体２０７はノズル群の幅分ずつ搬送される。インク色選択部２０３は、生成された走査データに基づき、記録ヘッド２０５に搭載されるインク色の中から、該当インク色を選択し、選択されたインクは記録素子（ノズル）を介して記録媒体２０７に吐出される。
【００１７】
次に、ＬＵＴ修正装置３００について説明する。本実施形態のＬＵＴ修正装置３００は、一般的な情報処理装置（パーソナルコンピュータ）によって実現される。ＬＵＴ作成部３０１は、修正対象となる初期色分解ＬＵＴを作成する。本実施形態の初期色分解ＬＵＴは、例えば、２５６×２５６×２５６の格子点を有し、これら格子点の間隔が等しい３次元色分解ＬＵＴである。初期ＬＵＴ格納部３０２は、修正対象となる初期色分解ＬＵＴが格納される。ＬＵＴ修正部３０３は、初期ＬＵＴ格納部に格納された初期色分解ＬＵＴの修正を行う。また、色分解ＬＵＴ格納部１０４に格納されている既存の色分解ＬＵＴを読み出して、この読み出した色分解ＬＵＴに対して修正を行ってもよい。ＬＵＴ修正方法の詳細は後述する。
【００１８】
（ＬＵＴ修正装置のハードウェア構成）
図２は、本実施形態におけるＬＵＴ修正装置３００のハードウェア構成例を示すブロック図である。ＬＵＴ修正装置３００は、ＣＰＵ３１１、ＲＡＭ３１２、ＲＯＭ３１３、グラフィックコントローラ３１４、表示部３１５、ＨＤＤ（Ｈａｒｄ  Ｄｉｓｋ  Ｄｒｉｖｅ）３１６、を含んで構成される。さらに、ＬＵＴ修正装置３００は、外部接続Ｉ／Ｆ３１７、ネットワークＩ／Ｆ３１８とを含み、各構成部はバス３１９を介して通信可能に接続されている。ＣＰＵ３１１は演算回路からなり、ＬＵＴ修正装置３００を統括制御する。ＣＰＵ３１１はＲＯＭ３１３またはＨＤＤ３１６に記憶されたプログラムをＲＡＭ３１２に読み出し、種々の処理を実行する。ＲＯＭ３１３は、ＬＵＴ修正装置３００の制御に用いられるシステムプログラムなどを記憶する。グラフィックコントローラ３１４は、表示部３１５に表示させるための画面を生成する。ＨＤＤ３１６は記憶領域としての機能を有し、種々の処理を実行するアプリケーションプログラムなどを記憶する。ＨＤＤ３１６は記憶装置の一例であり、ＨＤＤ以外にもＳＳＤ（Ｓｏｌｉｄ  Ｓｔａｔｅ  Ｄｒｉｖｅ）などにより構成することができる。外部接続Ｉ／Ｆ３１７は、ＬＵＴ修正装置３００に種々の機器を接続するためのインタフェースである。例えば、外部接続Ｉ／Ｆ３１７を介して、画像処理装置１００、ディスプレイ、キーボード、マウスなどを接続することができる。ネットワークＩ／Ｆ３１８は、ＣＰＵ３１１の制御に基づいて画像処理装置１００などとネットワークを介して通信を行う。図２に示されるＬＵＴ修正装置３００のハードウェア構成は、情報処理装置として一般的な構成となっている。そのため、画像処理装置１００についても、図２に示されるハードウェア構成を有する情報処理装置によって実現することができる。
【００１９】
（画像形成手順）
図３は、本実施形態において、画像処理装置１００が画像データの入力を受け付けてから画像形成装置２００が画像を出力するまでの一連の処理を示したフローチャートである。図３に示されるフローチャートの処理は、画像処理装置１００および画像形成装置２００に搭載されるＡＳＩＣなどの回路によって実行される。以下の各記号Ｓは、フローチャートにおけるステップであることを意味する。図３以降のフローチャートにおいても同様である。
【００２０】
Ｓ３０１において、画像処理装置１００の外部接続Ｉ／Ｆを介して、ＲＧＢ画像データが入力される。  
Ｓ３０２において、カラーマッチング処理部１０１は、３次元カラーマッチングＬＵＴを参照して、Ｓ３０１で入力されたＲＧＢ画像データのカラーマッチング処理を行う。  
Ｓ３０３において、色分解処理部１０３は、カラーマッチング処理において変換された画像データから、インク値画像データを生成する。色分解処理部１０３は、色分解ＬＵＴ格納部１０４に格納されている３次元の色分解ＬＵＴを参照して上記色分解処理を行う。
Ｓ３０４において、ハーフトーン処理部１０５は、色分解処理後のインク値画像データを２値画像データに変換する。２値画像データは、画像全体または単位記録領域ごとのバンド幅などの任意のサイズで画像形成装置２００に出力される。  
Ｓ３０５において、パス分解処理部２０１は、画像処理装置１００から受信した２値画像データを走査データに変換する。  
Ｓ３０６において、インク色選択部２０３は、走査データに適合するインク色を選択する。記録ヘッド２０５は、記録媒体２０７に対して移動しつつ、選択されたインク色に応じた各ノズルを一定の間隔で駆動する。記録媒体２０７は、記録ヘッド２０５の移動（走査）ごとに所定量だけ搬送される。上述の記録ヘッド２０５の駆動、および、記録媒体２０７の搬送が繰り返されることにより、記録媒体２０７上に画像が形成される。画像形成（Ｓ３０６）が終了すると、本フローチャートの処理は終了する。
【００２１】
（仮想色材量）
ここで、色分解ＬＵＴを修正する手法の説明に先立ち、仮想色材量について説明する。本実施形態では、仮想色材を、減法混色の３原色であるイエロー、マゼンタ、シアンの３色の色材とする。ここで、分光反射率Ｒｅｆ（λ）をｎ個（ｎは３以上の整数）の波長ブロックに区切り、各波長ブロック内の分光反射率を平均化した値をブロック反射率Ｒｅｆ１、Ｒｅｆ２・・・Ｒｅｆｎと定義する。さらに、ブロック反射率Ｒｅｆ１、Ｒｅｆ２・・・Ｒｅｆｎの対数をブロック濃度と定義する。具体的には、ブロック反射率Ｒｅｆ１、Ｒｅｆ２・・・Ｒｅｆｎを、以下の式（１）により変換した値Ｄ１、Ｄ２・・・Ｄｎがブロック濃度と定義される。  
Ｄ＝−ｌｏｇ１０（Ｒｅｆ）  ・・・  式（１）
【００２２】
本実施形態では、分光反射率Ｒｅｆ（λ）を、上記３色の色材が主に吸収する光の波長帯に対応する３つの波長ブロックに分割する。このとき、イエローインクが主に吸収する波長帯（３８０ｎｍ〜４８０ｎｍ）に対応するブロック濃度をＤｙとする。同様に、マゼンタインクが主に吸収する波長帯（４８０ｎｍ〜５８０ｎｍ）に対応するブロック濃度をＤｍとし、シアンインクが主に吸収する波長帯（５８０ｎｍ〜７３０ｎｍ）に対応するブロック濃度をＤｃとする。そして、ブロック濃度Ｄｙ、Ｄｍ、Ｄｃそれぞれに対応する理想的な分光反射率Ｒｅｆｙｉ（λ）、Ｒｅｆｍｉ（λ）、Ｒｅｆｃｉ（λ）を有するｎ種類（３種類）の仮想色材をｙｉ、ｍｉ、ｃｉと定義する。
【００２３】
図４は、仮想色材の分光反射率を示すグラフである。図４（ａ）は、仮想色材ｙｉが、主にイエローインクによって吸収される波長帯（３８０ｎｍ〜４８０ｎｍ）の光のみを吸収し、当該波長帯以外の波長の光を１００％反射することを示している。図４（ｂ）は、仮想色材ｍｉが、主にマゼンタインクによって吸収される波長帯（４８０ｎｍ〜５８０ｎｍ）の光のみを吸収し、当該波長帯以外の波長の光を１００％反射することを示している。図４（ｃ）は、仮想色材ｃｉが、主にシアンインクによって吸収される波長帯（５８０ｎｍ〜７３０ｎｍ）の光のみを吸収し、当該波長帯以外の波長の光を１００％反射することを示している。
【００２４】
上述のとおり、仮想色材ｙｉ、ｍｉ、ｃｉのブロック濃度Ｄｙｉ、Ｄｍｉ、Ｄｃｉは上記式（１）により求められる。例えば、仮想色材ｙｉ、ｍｉ、ｃｉが主に吸収する波長帯の反射率が１０％であったとする。このとき、仮想色材ｙｉのブロック濃度はＤｙｉ＝１．０、Ｄｍｉ＝０．０、Ｄｃｉ＝０．０となる（図５（ａ）を参照）。同様に、仮想色材ｍｉのブロック濃度はＤｙｉ＝０．０、Ｄｍｉ＝１．０、Ｄｃｉ＝０．０、仮想色材ｃｉのブロック濃度ＤｃｉはＤｙｉ＝０．０、Ｄｍｉ＝０．０、Ｄｃｉ＝１．０となる（図５（ｂ）および図５（ｃ）を参照）。
【００２５】
ところで、光散乱を無視できる範囲において、色材量（色材の厚さ）と光学濃度とが比例することがＬａｍｂｅｒｔ法則として既に知られている。この法則が、記録媒体上の仮想色材ｙｉ、ｍｉ、ｃｉに対して常に成り立つとすると、ｙｉ、ｍｉ、ｃｉそれぞれのブロック濃度は記録媒体上の単位面積当たりの色材量に比例する。すなわち、任意のブロック濃度Ｄｙ、Ｄｍ、Ｄｃと仮想色材量Ｖｙｉ、Ｖｍｉ、Ｖｃｉ［％］とは、以下の式（２−１）〜式（２−３）に従い、相互に線形変換が可能である。  
Ｖｙｉ＝（Ｄｙ／Ｄｙｉ）×α  ・・・  式（２−１）
Ｖｍｉ＝（Ｄｍ／Ｄｍｉ）×α  ・・・  式（２−２）
Ｖｃｉ＝（Ｄｃ／Ｄｃｉ）×α  ・・・  式（２−３）
【００２６】
式（２−１）〜式（２−３）において、Ｄｙｉ、Ｄｍｉ、Ｄｃｉはそれぞれ仮想色材のブロック濃度（図５（ａ）〜（ｃ）を参照）であり、αは仮想色材ｙｉ、ｍｉ、ｃｉの濃度に関する比例定数である。また、仮想色材量Ｖ［％］は、平均してＶ％の確率で、仮想色材ｙｉ、ｍｉ、ｃｉが対応画素に打ち込まれてインクドットが形成されることを意味している。ここまでの説明に従えば、任意の分光反射率Ｒｅｆ（λ）は、ブロック反射率を求めた後、式（１）よりブロック濃度Ｄｙ、Ｄｍ、Ｄｃに変換できる。さらに式（２−１）〜式（２−３）により、ブロック濃度Ｄｙ、Ｄｍ、Ｄｃとαとから仮想色材量Ｖｙｉ、Ｖｍｉ、Ｖｃｉ［％］に一意に変換することができる。
【００２７】
ここで、図６（ａ）〜図６（ｃ）を参照して、実色材量を仮想色材量に変換する具体例を説明する。実色材量から仮想色材量に変換するために、まず、画像形成装置２００が具備するＣ、Ｍ、Ｙ、Ｋ、Ｌｃ、Ｌｍの各インクの単色パッチが出力される。次いで、分光測色器などを用いて、出力された単色パッチの分光反射率Ｒｅｆ（λ）が測定される。図６（ａ）は、実シアンインクの単色パッチを測定した結果得られるブロック反射率Ｒｅｆｙ＿ｃ、Ｒｅｆｍ＿ｃ、Ｒｅｆｃ＿ｃを模式的に示したグラフである。本実施形態において、実シアンインクの単色パッチは、画像形成装置２００によって白紙の記録媒体に印刷され、そのインク打ち込み量Ｒ［％］は例えば２５％である。
【００２８】
図６（ｂ）は、図６（ａ）のブロック反射率Ｒｅｆｙ＿ｃ、Ｒｅｆｍ＿ｃ、Ｒｅｆｃ＿ｃに上記の式（１）の演算をそれぞれ施した結果得られるブロック濃度Ｄｙ、Ｄｍ、Ｄｃを模式的に示したグラフである。上述の通り、ブロック濃度Ｄｙ、Ｄｍ、Ｄｃと仮想色材量Ｖｙｉ、Ｖｍｉ、Ｖｃｉとは相互に線形変換が可能であるから、式（２−１）〜式（２−３）に従って、ブロック濃度Ｄｙ、Ｄｍ、Ｄｃから仮想色材量Ｖｙｉ、Ｖｍｉ、Ｖｃｉを算出することができる。図６（ｃ）は、実シアンインクについて算出される仮想色材量Ｖｙｉ、Ｖｍｉ、Ｖｃｉの一例を示している。本実施形態において、実シアンインクについて算出される仮想色材量は、Ｖｙｉ＝７、Ｖｍｉ＝９．５、Ｖｃｉ＝２０［％］である。なお、図６（ｃ）の例では、説明の便宜上α＝１００［％］として算出された仮想色材量が示されている。
【００２９】
さらに、上述の式（２−１）〜式（２−３）に従って求めた仮想色材量Ｖｙｉ、Ｖｍｉ、Ｖｃｉが、単位打ち込み量あたりの仮想色材量Ｖｙｉ´、Ｖｍｉ´、Ｖｃｉ´に換算される。この換算は、仮想色材量Ｖｙｉ、Ｖｍｉ、Ｖｃｉを、単色パッチのインク打ち込み量Ｒ［％］で除算することにより求められる。具体的には、以下の式（３−１）〜式（３−３）に従って単位打ち込み量あたりの仮想色材量Ｖｙｉ´、Ｖｍｉ´、Ｖｃｉ´を算出することができる。
【００３０】
Ｖｙｉ´＝Ｖｙｉ／Ｒ  ・・・  式（３−１）
Ｖｍｉ´＝Ｖｍｉ／Ｒ  ・・・  式（３−２）
Ｖｃｉ´＝Ｖｃｉ／Ｒ  ・・・  式（３−３）
【００３１】
以上説明した通り、実色材の単色パッチを測定した結果得られるブロック反射率に基づいて、その実色材量を仮想色材量に変換することができる。本実施形態では、画像形成装置２００が具備する実色材Ｃ、Ｍ、Ｙ、Ｋ、Ｌｃ、Ｌｍそれぞれについて、実色材量を仮想色材量に換算する処理が予め行われているものとする。このとき、分光反射率Ｒｅｆ（λ）の測定用の単色パッチは、インク打ち込み量Ｒ＝２５％で白紙に印刷されるものとし、単位打ち込み量あたりの仮想色材量Ｖｙｉ´、Ｖｍｉ´、Ｖｃｉ´は、実色材の種別ごとにＨＤＤ３１６などの記憶領域に保存される。
【００３２】
（初期色分解ＬＵＴの作成）
次に、本実施形態において、初期色分解ＬＵＴの作成手順について、図７および図８を参照して説明する。本実施形態における初期色分解ＬＵＴは、ユーザ指示に基づいてＬＵＴ修正装置３００のＬＵＴ作成部３０１によって作成される。
【００３３】
まず、作成対象となる初期色分解ＬＵＴの詳細について、図７を参照して説明する。図７（ａ）は、本実施形態における初期色分解ＬＵＴの構成を示す図である。図７（ａ）に示されるように、初期色分解ＬＵＴは、入力ＲＧＢ信号値について、ＲＧＢ３次元色空間上の立方体（以下「色立方体」とも記す）を格子状に分割した格子点を有する。本実施形態の初期色分解ＬＵＴは、２５６×２５６×２５６の格子点を有する、いわゆる全点ＬＵＴである。色立方体における各格子点について、画像形成装置２００が具備するＣ、Ｍ、Ｙ、Ｋ、Ｌｃ、Ｌｍの色材量（色材値）が格納されている。なお、図７（ａ）のＬＵＴ例では、Ｒ軸Ｇ軸Ｂ軸それぞれについて、ｐ０：０、ｐ１：３２、ｐ２：６４、ｐ３：９６、ｐ４：１２８、ｐ５：１６０、ｐ６：１９２，ｐ７：２２４、ｐ８：２５５の主要格子点が均等に配列されている例が示されている。通常、色分解処理で参照される３次元ＬＵＴは、その記憶容量を節約するために、主要格子点のみに色材値を格納し、これら格子点間の色材値は補間演算により算出される。初期色分解ＬＵＴの格子点を主要格子点のみに間引く処理については後述する。
【００３４】
図７（ｂ）は、ＲＧＢ色空間の概略図であり、図７（ａ）における色立方体の８頂点（Ｗ、Ｃ、Ｍ、Ｙ、Ｒ、Ｇ、Ｂ、Ｋ）と、各頂点同士を結んだラインとを示した図である。ここで、ＲＧＢ値それぞれについて８ビットの入力画像データが入力される場合、図７（ｂ）の色立方体において、Ｗ、Ｃ、Ｍ、Ｙ、Ｒ、Ｇ、Ｂ、Ｋの各頂点における座標は以下の通りに表される。  
Ｗ＝（２５５，２５５，２５５）
Ｃ＝（０，２５５，２５５）
Ｍ＝（２５５，０，２５５）
Ｙ＝（２５５，２５５，０）
Ｒ＝（２５５，０，０）
Ｇ＝（０，２５５，０）
Ｂ＝（０，０，２５５）
Ｋ＝（０，０，０）
【００３５】
まず、図７（ａ）（ｂ）に示されるＲＧＢ色空間において、再現色域が最大となるように、各頂点Ｗ，Ｃ，Ｍ，Ｙ，Ｒ，Ｇ，Ｂ，Ｋに対応する色材量が決定される。例えば、Ｗは紙白であるので、Ｗに対応する色材量は全て０に設定される。１次色の頂点Ｃ，Ｍ，Ｙ，Ｋについては、混色による色の濁りを防ぐため、それぞれ対応する１種類の色材のみの色材量が設定され、当該１種類の色材以外の色材量には全て０に設定される。このとき、Ｃ，Ｍ，Ｙ，Ｋに対応する色材量は、記録媒体の種類や印刷モードなどに応じて印刷可能な最大載り量に対応する。２種類の色材の混合により実現される２次色の頂点Ｒ，Ｇ，Ｂについては、例えばＲの場合、ＭおよびＹの２種の色材を用い、その合計が最大色材量となる複数の組合せのカラーパッチデータを生成する。これらカラーパッチデータは画像形成装置２００で印刷され、その印刷結果をユーザが目視で確認し、Ｒに対応する色として好ましい色となる色材量の組合せが選択される。あるいはまた、測色器を用いてカラーパッチデータの印刷結果を測色して、ＭおよびＹの色相の中間の色相となる色材量の組合せが決定されるようにしてもよい。Ｇ，ＢについてもＲの場合と同様に、ＧはＹとＣとの組み合わせ、ＢはＣとＭとの組み合わせに基づいて、最適な色材量の組み合わせを決定することができる。以下、本明細書において、頂点Ｃ，Ｍ，Ｙ，Ｒ，Ｇ，Ｂをカラープライマリと記す。
【００３６】
次に、頂点Ｗとカラープライマリとを結ぶＷ−Ｃ，Ｗ−Ｍ，Ｗ−Ｙ，Ｗ−Ｒ，Ｗ−Ｇ，Ｗ−Ｂのラインにおける色材量が決定される。図８（ａ）は、Ｗ−Ｃラインにおける色変換特性の例を示す図であり、横軸は入力ＲＧＢ信号値を示し、縦軸は実色材量と最大色材量とを示している。また、図８（ａ）のグラフにおける符号（◇および△）は、色立方体（図７（ａ））における主要格子点（ｐ０〜ｐ８）に対応する色材量であることを意味する。図８（ａ）において、頂点Ｗ（ｐ０）に対応する全色材の色材量は０である。一方、頂点Ｃ（ｐ８）では、色材Ｃの色材量は最大となり、その他の色材の色材量は０である。
【００３７】
ここで、印刷出力における低濃度領域の粒状性を向上させるために、本実施形態の画像形成装置２００が具備するインクの種類には、淡シアン（Ｌｃ）が含まれる。図８（ａ）に示されるグラフでは、頂点Ｗ（ｐ０）を起点として、色材Ｌｃの色材量が最大色材量となるまでが徐々に増加する様子が示されている。そして、色材Ｌｃの色材量が最大色材量に到達すると、色材Ｌｃの色材量は、頂点Ｃ（ｐ８）にかけて徐々に減少する。色材ＬＣの色材量が減少することに応じて、今度はより発色のよい色材である色材Ｃの色材量が頂点Ｃ（ｐ８）にかけて徐々に増加する。このとき、色材Ｌｃの増減が切り替わる点（色材Ｃの増加が開始する点）を、「変曲点」と記す。本実施形態では、ＬＵＴ修正装置３００の表示部３１５に、例えば図８（ａ）に示されるグラフをＵＩとして表示することができる。ユーザは、このＵＩを介して、格子点に対応する色材の実色材量を調整することにより、例えば、例えばＷ−Ｃラインにおける入力ＲＧＢ値と色材量との対応関係（色変換特性）を決定することができる。また、ユーザが決定した色変換特性に基づき画像形成装置２００にカラーパッチを印刷させ、このカラーパッチの印刷結果における粒状性や色味をユーザが目視評価できるようにしてもよい。本実施形態では、Ｗ−Ｃラインにおける色材量を決定する処理について説明したが、Ｗ−Ｍ，Ｗ−Ｙ，Ｗ−Ｒ，Ｗ−Ｇ，Ｗ−Ｂのラインにおける色材量も同様に決定することができる。
【００３８】
次に、カラープライマリと頂点Ｋとを結ぶ各ラインにおける色材量、および、カラープライマリとカラープライマリとを結ぶ各ラインにおける色材量が決定される。カラープライマリと頂点Ｋとを結ぶライン、カラープライマリとカラープライマリとを結ぶラインにおいては、印刷物による再現色域をより広くするため、０〜２５５の各格子点について、対応する全ての色材の総色材量が最大色材量となるように設計される。カラープライマリと頂点Ｋとを結ぶラインにおける色変換特性の例と、カラープライマリとカラープライマリとを結ぶラインにおける色変換特性の例とを、図８（ｂ）と、図８（ｃ）とにそれぞれ示す。
【００３９】
以上説明した通り、初期色分解ＬＵＴの作成においては、２５６×２５６×２５６の格子点を有する初期色分解ＬＵＴ（色立方体）において、先ず各頂点同士を結ぶ各ラインにおける色材量が決定される。次いで、初期色分解ＬＵＴ（図７（ａ））の内部における色材量が、補間処理により順次決定される。補間処理の具体的手法としては、例えば公知の有限要素法を用いた非線形補間法等を適用することができる。
【００４０】
（色分解ＬＵＴの修正手順）
次に、本実施形態における色分解ＬＵＴの修正手順について、図９のフローチャートを参照して説明する。図９に示されるフローチャートの処理は、ＬＵＴ修正装置３００のＬＵＴ修正部３０３によって実行される。
【００４１】
Ｓ９０１において、初期色分解ＬＵＴが取得される。本実施形態では、ＬＵＴ作成部３０１によって作成され、初期ＬＵＴ格納部３０２に格納されている初期色分解ＬＵＴが読み出される。なお、初期色分解ＬＵＴは、図７〜図８で説明した手法によって作成された初期色分解ＬＵＴが読み出されるが、色分解ＬＵＴ格納部１０４に格納されている既存の色分解ＬＵＴが取得されてもよい。
【００４２】
Ｓ９０２において、初期色分解ＬＵＴの主要格子点以外の格子点が間引かれる。Ｓ９０２における格子点の間引き処理は、色分解処理で参照される３次元ＬＵＴの記憶容量の節約のために行われる。本実施形態では、２５６×２５６×２５６の格子点を有する初期色分解ＬＵＴから、Ｒ軸Ｇ軸Ｂ軸それぞれ９つの主要格子点（すなわち９×９×９の主要格子点）以外の格子点が間引かれる。すなわち、上記主要格子点以外の格子点には、入力ＲＧＢ信号値に対応する色材量が格納されない状態となる。そして、Ｓ９０２で間引き処理された色分解ＬＵＴは、その後色分解ＬＵＴ格納部１０４に格納される。
【００４３】
Ｓ９０３において、修正対象となる着目格子点が初期化される。本実施形態では、初期色分解ＬＵＴにおける頂点Ｗ（ＲＧＢ＝（０，０，０））が最初の着目格子点として設定される。
【００４４】
Ｓ９０４において、評価対象位置となる着目格子点および着目格子点の近傍位置に配置される格子点（以下、「近傍格子点」と記す）に保持されている実色材量（実色材値）が取得される。本実施形態では、着目格子点を中心として、Ｒ軸Ｇ軸Ｂ軸方向の３次元における格子点群（３×３×３＝２７点）にそれぞれ保持されている実色材量が取得される。なお、近傍格子点が初期の色分解ＬＵＴの定義域外に位置する場合は、該当近傍格子点に対応する実色材量は取得されない。
【００４５】
Ｓ９０５において、着目格子点における実色材量および近傍格子点における実色材量が、それぞれ仮想色材量を示す仮想色材値に変換される。図６（ａ）〜図６（ｃ）を参照して説明した通り、本実施形態では、画像形成装置２００が具備する実色材Ｃ、Ｍ、Ｙ、Ｋ、Ｌｃ、Ｌｍそれぞれについて、単位打ち込み量あたりの仮想色材量Ｖｙｉ´、Ｖｍｉ´、Ｖｃｉ´が予め算出されている。着目格子点における実色材量および近傍格子点における実色材量は、予め算出されている単位打ち込み量あたりの仮想色材量Ｖｙｉ´、Ｖｍｉ´、Ｖｃｉ´に基づいて、それぞれ仮想色材値（仮想色材量）に変換される。例えば、格子点ｐにおける色材ｉの仮想色材量をＶｙｉ、Ｖｍｉ、Ｖｃｉ、格子点ｐにおける色材ｉの実色材量をＱ（ｐ，ｉ）［％］とすると、Ｖｙｉ、Ｖｍｉ、Ｖｃｉは以下の式（４−１）〜式（４−３）に従って算出することができる。  
Ｖｙｉ（ｐ，ｉ）＝Ｖｙｉ´（ｉ）＊Ｑ（ｐ，ｉ）  ・・・  式（４−１）
Ｖｍｉ（ｐ，ｉ）＝Ｖｍｉ´（ｉ）＊Ｑ（ｐ，ｉ）  ・・・  式（４−２）
Ｖｃｉ（ｐ，ｉ）＝Ｖｃｉ´（ｉ）＊Ｑ（ｐ，ｉ）  ・・・  式（４−３）
【００４６】
実色材量から仮想色材量への変換処理について、図８（ａ）を参照して説明する。上述の通り、Ｗ−Ｃラインにおける頂点Ｃ（格子点ｐ８）では、Ｃインクの色材量が２５５（１００％）、Ｃインク以外の色材量はすべて０（０％）であり、（Ｃ、Ｍ、Ｙ、Ｋ、Ｌｃ、Ｌｍ）＝（１００，０，０，０，０，０，）［％］で表される。このとき、格子点ｐ８（頂点Ｃ）における仮想色材量は、上記式（４−１）〜式（４−３）に従って以下の通りに算出される。  
Ｖｙｉ（ｐ８）＝０．２８＊１００＝２８
Ｖｍｉ（ｐ８）＝０．３８＊１００＝３８
Ｖｃｉ（ｐ８）＝０．８０＊１００＝８０
【００４７】
Ｗ−Ｃラインにおける格子点ｐ６では、Ｃインクの色材量が１０２（４０％）、Ｌｃインクの色材量が１５３（６０％）、これら以外の色材量はすべて０（０％）であり、（Ｃ、Ｍ、Ｙ、Ｋ、Ｌｃ、Ｌｍ）＝（４０，０，０，０，６０，０，）［％］で表される。このとき、格子点ｐ６における仮想色材量は、Ｃインクインクについて算出される仮想色材量と、Ｌｃインクについて算出される仮想色材量とが合算された値となる。
【００４８】
Ｖｙｉ（ｐ６）＝０．２８＊４０＋０．１６＊６０＝２０．８
Ｖｍｉ（ｐ６）＝０．３８＊４０＋０．２２＊６０＝２８．４
Ｖｃｉ（ｐ６）＝０．８０＊４０＋０．４８＊６０＝６０．８
なお、上記計算例において、実Ｌｃインクについて予め算出されている単位打ち込み量あたりの仮想色材量Ｖｙｉ´、Ｖｍｉ´、Ｖｃｉ´は、Ｖｙｉ´＝０．１６、Ｖｍｉ´＝０．２２、Ｖｃｉ´＝０．４８である。単位打ち込み量あたりの仮想色材量Ｖｙｉ´、Ｖｍｉ´、Ｖｃｉ´は、図６（ａ）〜図６（ｃ）および式（３−１）〜式（３−３）で説明した手法により予め算出されているものとする。
【００４９】
次に、Ｓ９０６において、着目格子点における評価値が算出される。本実施形態では、ＲＧＢ色空間における仮想色材量の均等性が評価値として用いられる。具体的には、着目格子点における仮想色材量と、隣接格子点における仮想色材量との間の変化を変位ベクトルで表し、当該ベクトルの和の大きさが評価値として算出される。この評価値は、ＲＧＢ色空間における着目格子点の場所ごとに算出方法が異なる。
【００５０】
図１０は、本実施形態において評価値を算出する方法を説明する図である。図１０（ａ）は、Ｗ−Ｃラインにおける着目格子点と隣接格子点とを示した模式図である。このとき、２つの隣接格子点の間に着目格子点が位置する位置関係は１通りしか存在しないため、評価値Ｅは以下の式（７）に従って算出される。  
Ｅ＝|（Ｖｒ１−Ｖｐ)＋（Ｖｒ０−Ｖｐ）|／|Ｖｒ１−Ｖｒ０|  ・・・  式（５）
【００５１】
式（５）においてＶｐは着目格子点における仮想色材量であり、Ｖｒ０、Ｖｒ１はそれぞれ隣接格子点における仮想色材量である。
【００５２】
図１０（ｂ）は、Ｇ軸とＲ軸とによって区画される初期色分解ＬＵＴ（色立方体）の表面における着目格子点と隣接格子点とを示した模式図である。このとき、２つの隣接格子点の間に着目格子点が位置する位置関係は２通り存在するので、（Ｖｒ０，Ｖｐ，Ｖｒ１）、（Ｖｇ０，Ｖｐ，Ｖｇ１）それぞれの組み合わせについて式（５）を用いた演算が行われる。そして、図１０（ｂ）に示されるように、上記演算結果の和が評価値Ｅとして得られる。
【００５３】
図１０（ｃ）は、初期色分解ＬＵＴ（色立方体）の内部における着目格子点と隣接格子点とを示した図である。このとき、着目格子点が２つの隣接格子点の間に位置する位置関係は３通り存在するので、（Ｖｒ０，Ｖｐ，Ｖｒ１）、（Ｖｇ０，Ｖｐ，Ｖｇ１）、（Ｖｂ０，Ｖｐ，Ｖｂ１）それぞれの組み合わせについて式（５）を用いた演算が行われる。そして、図１０（ｃ）に示されるように、上記演算結果の和が評価値Ｅとして得られる。
【００５４】
Ｓ９０７において、図１０（ａ）〜（ｃ）の手法により算出された評価値Ｅを用いて、着目格子点における実色材量を最適化する。具体的には、まず、評価値Ｅが小さくなるように、着目格子点における仮想色材量が増減される。次いで、増減後の仮想色材量から実色材量に再び変換される。
【００５５】
ここで、Ｓ９０７における実色材量の最適化処理を、図１１（ａ）および図１１（ｂ）を参照して説明する。図１１（ａ）は、Ｓ９０２において主要格子点以外の格子点が間引かれた後の、初期色分解ＬＵＴにおける色変換特性を示す図である。なお、図１１（ａ）の色変換特性において、格子点間の色材量は線形補間より算出されているものとする。一方、図１１（ｂ）は、図１１（ａ）のグラフにおける実Ｃインクの色材量と、実Ｌｃインクの色材量とから変換された（Ｓ９０５）仮想色材量の一例を示す。なお、説明の便宜のため、図１１（ｂ）のグラフでは、仮想色材量Ｖｙｉ、Ｖｍｉ、Ｖｃｉのうち、仮想色材量Ｖｃｉの推移のみが示されている。図１１（ｂ）に示されるグラフにおいて、増減される前の仮想色材量の一例（Ｓ９０７）が破線によって、増減された後の仮想色材量の一例（Ｓ９０７）が実線によってそれぞれ示されている。
【００５６】
図１１（ｂ）のグラフにおいて、破線部は、図１１（ａ）の色変換特性において、実Ｃインクの色材量および実Ｌｃインクの色材量の線形性が崩れてしまっていることを示している。すなわち、図１１（ａ）のグラフにおいて、変曲点近傍の実Ｌｃインクの色材量は、図８（ａ）のグラフと比較して、補間演算により減少してしまっている。同様に、図１１（ａ）のグラフにおいて、変曲点近傍の実Ｃインクの色材量は、図８（ａ）のグラフと比較して、補間演算により増加してしまっている。このように、仮に全点ＬＵＴを参照していれば取得することができた色材量と、補間演算により算出された色材量との間に発生する誤差を、本実施形態では「補間誤差」と記す。本実施形態では、実色材量と仮想色材量とは相互に線形の関係が成立するため、補間誤差などによって実インクの色材量の線形性が崩れてしまった場合であっても、色材量の急峻な変化を仮想色材量の推移によって表すことができる。そして、本実施形態では仮想色材量の特性に着目し、評価値Ｅに基づいて着目格子点における仮想色材量を調整し、当該調整後の仮想色材量から実色材量に再度変換処理を行うことにより、実色材量を最適化する処理を実行するのである。上述の通り、図１１（ｂ）のグラフにおいて、評価値Ｅに基づいて増減された仮想色材量の推移が実線で示されている。Ｓ９０７において、増減された仮想色材量から、実色材量への変換処理が実行される。最適化後の実Ｌｃインクの色材量のグラフ線と、最適化後の実Ｃインクの色材量のグラフ線とが、図１１（ａ）に示されている。
【００５７】
Ｓ９０８において、着目格子点における実色材量を、Ｓ９０７で最適化された実色材量に更新する。このように、本実施形態では、評価値Ｅが小さくなるように着目格子点における実色材量が修正されることにより、初期色分解ＬＵＴが、入力ＲＧＢ信号値に対する実色材量がより良好な階調性となるように修正される。
【００５８】
Ｓ９０９において、全格子点について修正が完了したか否かが判定される。修正が完了している場合（Ｓ９０９：ＹＥＳ）、本フローチャートの処理を終了する。修正が完了していない場合（Ｓ９０９：ＮＯ）、着目格子点を次の格子点に移動させ、再びＳ９０４以降の処理を繰り返す。
【００５９】
以上説明した通り、本実施形態のＬＵＴ修正装置３００によれば、実色材量との関係が線形となる仮想色材量を用いて色分解ＬＵＴの変換特性を評価することができる。さらに、この評価結果に基づいて、入力画像信号値に対する実色材量の変化が滑らかとなるように、色分解ＬＵＴを修正することができる。つまり、本実施形態のＬＵＴ修正装置は、色変換パラメータによって実現される階調性をより良好にすることができる、という効果を奏する。
【００６０】
（変形例）
なお、上述の実施形態は例示であり、本実施形態を以下の変形例に適用することができる。本実施形態では、記録ヘッド２０５に搭載されるインクの種類をＣ、Ｍ、Ｙ、Ｋ、Ｌｃ、Ｌｍの６色としたが、記録ヘッド２０５に搭載されるインクの種類はこれらに限定されない。具体的には、記録ヘッド２０５は、レッド（Ｒ）、グリーン（Ｇ）などの特色インクや、白色インクを搭載してもよい。あるいはまた、記録ヘッド２０５は、無色透明のクリアインクや、金属調のメタリックインクなどを搭載してもよい。
【００６１】
また、本実施形態では、８ビットのＲＧＢカラー画像データが画像処理装置１００に入力される例を説明したが、モノクロ画像データやＣＭＹＫ画像データが画像処理装置１００に入力されてもよい。また、入力される画像データは、その画素ごとに色情報以外の情報を含んでもよく、例えば画素ごとに光沢情報を含んでもよい。この場合、インク付与量に基づく色味の階調性と同様に、発色の階調性についても変換特性の修正処理を実行することができる。
【００６２】
また、本実施形態では、色分解ＬＵＴにおけるグリッド数（Ｒ軸Ｇ軸Ｂ軸それぞれの格子点数）は、間引き前の色分解ＬＵＴについては２５６個、間引き後の色分解ＬＵＴについては９個であった。これら色分解ＬＵＴにおけるグリッド数は上記に限定されない。
【００６３】
また、入力画像信号値から実色材量への色変換特性を示す色変換パラメータは、上記色分解ＬＵＴ以外の形態であってもよく、例えば、行列演算や数式によって表現されてもよい。この場合、色変換特性（色変換パラメータ）の修正処理において、上記行列における係数や、数式における係数が修正される。
【００６４】
また、本実施形態では、ハーフトーン処理の一例として、ディザマトリクス法が適用される例が説明されたが、ディザマトリクス法以外にも、公知の誤差拡散法などの手法を適用してもよい。また、パス分解の手法についても、マルチパス方式だけではなく、シングルパス方式によってパス分解処理が行われてもよい。
【００６５】
また、本実施形態では、ＬＵＴ作成部３０１によって生成された初期色分解ＬＵＴに対して、ＬＵＴ修正部３０３が修正を行う実施例を説明したが、既存の色分解ＬＵＴ（例えば、色分解ＬＵＴ格納部１０４に格納されているＬＵＴなど）を用いてもよい。この場合、Ｓ９０１において、色分解ＬＵＴ格納部１０４に格納されている色分解ＬＵＴが読み出され、読み出された色分解ＬＵＴの格子点が例えば９×９×９の格子点を有する間引き後の色分解ＬＵＴであった場合、Ｓ９０２が省略される。
【００６６】
また、本実施形態では、仮想色材量の評価値を算出する手法として（Ｓ９０６）、ＲＧＢ色空間における変位ベクトルに基づいて評価値を算出する手法を説明したが、評価方法は上記に限定されない。例えば、仮想色材量空間における実色材量の変位ベクトルに基づいて評価値を算出してもよい。また、例えば、変位ベクトルを用いずに曲率に基づいて評価値を算出してもよい。また、変位ベクトルに基づく評価だけではなく、他の評価手法と組み合わせてもよい。例えば、初期色分解ＬＵＴにおける実色材量と、間引きＬＵＴにおける実色材量との差分を導出し、変位ベクトルに基づく評価値と、実色材量との差分値の線形和に基づいて評価処理を実行してもよい。本実施形態では、実色材量を最適化処理の対象としたため（Ｓ９０７）、最適化により色再現域が小さくなることがあるが、実色材量差を用いて評価を行うことにより、色再現域が小さくなることを低減することができる。また、本実施形態では、実色材量の最適化（Ｓ９０７）において、準ニュートン法を適用する実施例を説明したが、例えば公知のレーベンバーグ・マーカート法や遺伝的アルゴリズムなど任意の最適化手法を適用することができる。また、最適化の対象は実色材量に限定されず、例えば、色分解ＬＵＴにおける格子点のＲＧＢ値を最適化対象にした場合にも効果が得られる。この場合、ＲＧＢ空間における色立方体の６表面においては、最適化後のＲＧＢ値を色立方体の表面内部に拘束することにより、色再現域の縮小を抑制することができる。
【００６７】
［実施形態２］
実施形態１では、９×９×９の格子点に間引かれた色分解ＬＵＴの色変換特性を修正する例について説明した。しかしながら、色分解ＬＵＴにおける格子点の間隔は、実施形態１の色分解ＬＵＴのように常に均等な間隔であるとは限らない。例えば、色分解処理の前に、１次元の非線形変換ＬＵＴを参照して入力ＲＧＢ信号値を非線形変換することにより、３次元色分解ＬＵＴの格子点位置を不均等化することができる。本実施形態では、仮想色材量を用いて、このような１次元の非線形変換ＬＵＴを修正する例について説明する。以下、実施形態１と共通する部分については説明を簡略化ないし省略し、以下では本実施形態に特有な点を中心に説明する。
【００６８】
（印刷システムの全体構成）
図１２は、本実施形態における印刷システム１の全体構成例を示すブロック図である。実施形態１の印刷システム１と異なる部分は、画像処理装置１００に、非線形変換部１２０１および非線形変換ＬＵＴ格納部１２０２が追加されている点である。非線形変換部１２０１は、カラーマッチング処理後のＲＧＢ画像データにおけるＲＧＢ信号値のそれぞれに非線形変換を行い、非線形変換されたＲＧＢ画像データを色分解処理部１０３に送信する。非線形変換処理は、画像形成装置２００が具備するインク色ごとに行われる。非線形変換処理において、非線形変換ＬＵＴ格納部１２０２に格納された非線形特性を備える１次元の非線形変換ＬＵＴが用いられる。本実施形態において、非線形変換ＬＵＴにおける格子点には、入力ＲＧＢ信号値に対応する固定値（出力値）が保持されており、これら固定値は色分解ＬＵＴにおける入力値として用いられる。
【００６９】
（非線形変換ＬＵＴ）
図１３は、本実施形態における非線形変換ＬＵＴの例を示す模式図である。図１３に示される各グラフは、横軸を入力ＲＧＢ信号値とし、縦軸を非線形変換後の出力ＲＧＢ信号値とした場合のＲ軸Ｇ軸Ｂ軸の非線形特性を示しており、それぞれＣ、Ｍ、Ｙ、Ｋ、Ｌｃ、Ｌｍのインク種類に対応している。図１３に示される各非線形変換ＬＵＴは、各軸それぞれ１７個の格子点を有し、それぞれの格子点は、８ビット（２５６階調）の入力ＲＧＢ画像データに対応する変換後の出力値が格納されている。格子点間の出力値は線形補間により算出される。
【００７０】
（非線形変換ＬＵＴの修正手順）
次に、本実施形態における非線形変換ＬＵＴの修正手順について、図１４のフローチャートを参照して説明する。図１４に示されるフローチャートの処理は、ＬＵＴ修正装置３００のＬＵＴ修正部３０３によって実行される。
【００７１】
Ｓ１４０１において、修正対象となる非線形変換ＬＵＴが選択される。本実施形態では、ＬＵＴ作成部３０１によって作成され、初期ＬＵＴ格納部３０２に格納されている非線形変換ＬＵＴのうち、まずＣインクのＲ軸非線形変換ＬＵＴが選択される。本実施形態では、初期ＬＵＴ格納部３０２から読み出される初期非線形変換ＬＵＴは、入力ＲＧＢ信号値と出力ＲＧＢ信号値とが等しい、線形の変換特性を有している。
【００７２】
Ｓ１４０２において、Ｓ１４０１で選択した非線形変換ＬＵＴにおける格子点の中から着目格子点ｐが選択される。非線形変換ＬＵＴにおける始点（入力ＲＧＢ信号値＝０の格子点）および終点（入力ＲＧＢ信号値＝２５５の格子点）の位置は固定されており、本実施形態では始点に隣接する格子点が着目格子点ｐとして最初に選択される。
【００７３】
Ｓ１４０３において、パラメータｓが初期化される。本実施形態において、パラメータｓは変動幅を表し、初期化後のパラメータｓの値は「８」となる。
【００７４】
Ｓ１４０４において、初期化されたパラメータｓに基づいて、着目格子点ｐにおける出力ＲＧＢ信号値が可変される。具体的には、着目格子点ｐにおける出力ＲＧＢ信号値に対して、初期化されたパラメータｓを加算する処理と、減算する処理とを実行する。そして、着目格子点ｐについて、出力ＲＧＢ信号値が変更されていないＬ０と、出力ＲＧＢ信号値にパラメータｓが加算されたＬ１と、出力ＲＧＢ信号値からパラメータｓが減算されたＬ２と、の３種類の中間ＬＵＴが作成される。
【００７５】
Ｓ１４０５において、Ｓ１４０４で作成された３種類の中間ＬＵＴと、色分解ＬＵＴとを用いて、評価対象ラインにおける実色材量がそれぞれ算出される。ここで、本実施形態における評価対象ラインとは、図７（ｂ）に示される各頂点Ｗ，Ｃ，Ｍ，Ｙ，Ｒ，Ｇ，Ｂ，Ｋ同士を結ぶラインをいい、評価対象ラインは２８通り存在する。Ｓ１４０５において、まず、評価対象ラインにおけるＲＧＢ入力信号値が３種類の中間ＬＵＴによってそれぞれ変換された後、さらに色分解ＬＵＴによって実色材量に変換される。この結果、１つの評価対象ラインについて、３種類の中間ＬＵＴ（Ｌ０〜Ｌ２）に対応する３種類の実色材量が取得される。なお、Ｓ１４０５において参照される色分解ＬＵＴは、色分解ＬＵＴ格納部１０４に格納されている色分解ＬＵＴが好適に用いられるが、実施形態１で説明した修正後の色分解ＬＵＴが参照されてもよい。色分解ＬＵＴ格納部１０４に格納されている色分解ＬＵＴ、および、修正後の色分解ＬＵＴは、それぞれ主要格子点以外の格子点が間引かれた９×９×９の３次元色分解ＬＵＴである。
【００７６】
Ｓ１４０６において、評価対象ラインにおける実色材量が仮想色材量に変換される。実色材量から仮想色材量への変換処理は、図６（ａ）〜図６（ｃ）の手法および式（４−１）〜式（４−３）に従って実行することができる。Ｓ１４０６における変換処理の結果、３種類の中間ＬＵＴ（Ｌ０〜Ｌ２）に対応する３種類の仮想色材量が取得される。
【００７７】
Ｓ１４０７において、Ｓ１４０６で算出された３種類の仮想色材量それぞれについて評価値Ｅ０、Ｅ１、Ｅ２が算出される。仮想色材量に基づいて評価値Ｅ０〜Ｅ２を算出する手法では、例えば、評価対象ラインにおける仮想色材量の曲率を２乗した値の最大値が用いられる。
【００７８】
Ｓ１４０８において、評価値Ｅ０〜Ｅ２がそれぞれ比較され、３種類の中間ＬＵＴ（Ｌ０〜Ｌ２）のうち、評価値が最良となる中間ＬＵＴが選択される。そして、初期非線形変換ＬＵＴが、この評価値が最良となる中間ＬＵＴに更新される。
【００７９】
Ｓ１４０９において、評価値が収束したか否かが判定される。本実施形態では、パラメータｓによって出力ＲＧＢ信号値が可変されていないＬ０に対応する評価値Ｅ０が算出された場合（Ｓ１４０７）、評価値が収束したと判定される。評価値が収束した場合（Ｓ１４０９：ＹＥＳ）、Ｓ１４１０に移行する。評価値が収束していない場合（Ｓ１４０９：ＮＯ）、再びＳ１４０４に戻る。
【００８０】
Ｓ１４１０において、パラメータｓの値が更新される。本実施形態では、現在のｓの値（例えば「８」）を２で割った値（例えば「４」）に更新される。
【００８１】
Ｓ１４１１において、ｓの値が１．０未満であるか否かが判定される。ｓの値が１．０以上の場合（Ｓ１４１１：ＮＯ）、再び１４０４に戻り、ｓの値が１．０未満の場合（Ｓ１４１１：ＹＥＳ）、Ｓ１４１２に移行する。
【００８２】
Ｓ１４１２において、全ての格子点について最適化処理が完了したか否かが判定される。全ての格子点について最適化処理が完了していない場合（Ｓ１４１２：ＮＯ）、再びＳ１４０２に戻り、全ての格子点について最適化処理が完了した場合（Ｓ１４１２：ＹＥＳ）、Ｓ１４１３に移行する。
【００８３】
Ｓ１４１３において、全ての種類の非線形変換ＬＵＴについて修正処理が完了したか否かが判定される。全ての種類の非線形変換ＬＵＴについて修正処理が完了していない場合（Ｓ１４１３：ＮＯ）、再びＳ１４０１に戻り、全ての種類の非線形変換ＬＵＴについて修正処理が完了している場合（Ｓ１４１３：ＹＥＳ）、本フローチャートの処理を終了する。
【００８４】
以上説明した通り、本実施形態のＬＵＴ修正装置３００によれば、仮想色材量を用いて、色分解処理の前に参照される１次元非線形変換ＬＵＴの変換特性を評価することができる。さらに、この評価結果に基づいて、色分解ＬＵＴにおける格子点の位置が適切となるように、１次元非線形変換ＬＵＴを修正することができる。
【００８５】
（変形例）
なお、上述の実施形態は例示であり、本実施形態を以下の変形例に適用することができる。本実施形態では、入力ＲＧＢ信号値の非線形変換処理をインク色ごとに行う例を説明したが、インク色に共通の処理としてもよい。また、一部のインク色のみについて、共通の非線形変換ＬＵＴを参照してもよい。また、本実施形態では、非線形変換ＬＵＴにおけるグリッド数は１７であったが、非線形変換ＬＵＴにおけるグリッド数は上記に限定されない。特に、非線形ＬＵＴのグリッド数が色分解ＬＵＴのグリッド数よりも多い場合、色分解ＬＵＴの格子点間でインク曲線を曲げることが可能になり、より好適な補間を実現することができる。また、非線形変換ＬＵＴにおいて、格子点間における出力ＲＧＢ信号値は線形補間で算出される例を説明したが、補間手法は上記に限られず、スプライン補間等の補間手法を適用することができる。
【００８６】
また、入力ＲＧＢ信号値から出力ＲＧＢ信号値への変換特性を示す非線形変換パラメータは、上記１次元非線形変換ＬＵＴ以外の形態であってもよく、例えば、行列演算や数式によって表現されてもよい。この場合、非線形変換パラメータの修正処理において、上記行列における係数や、数式における係数が修正される。
【００８７】
また、色分解処理部１０３が、非線形変換部１２０１を含む構成としてもよい。また、本実施形態では、評価ラインにおける仮想色材量の曲率を評価値として算出する例を説明したが、曲率の最大値ではなく、曲率の平均値を評価値として算出してもよいし、これら複数の評価値の線形和を適用することもできる。
【００８８】
［実施形態３］
実施形態１では、仮想色材量を用いて色分解ＬＵＴの色変換特性を評価し、この評価結果に基づいて、入力画像信号値に対する実色材量の変化が滑らかとなるように、色分解ＬＵＴを修正した。しかしながら、画像処理装置１００において最終的な実色材量が決定されるまでに、色分解処理以外の画像処理が行われる場合がある。このような処理の一例として、ＯＰＧ（Ｏｕｔｐｕｔ  Ｇａｍｍａ）処理が知られている。ＯＰＧ処理は、色材量と明度との線形性を改善するために行われる処理であり、ＯＰＧ処理において、１次元のＯＰＧＬＵＴが参照されて実色材量が補正される。つまり、良好な階調特性を得るためには、ＯＰＧＬＵＴにおける変換特性についても適切に設定される必要がある。本実施形態では、仮想色材量を用いてＯＰＧＬＵＴの変換特性を評価し、この評価結果に基づいて、入力画像信号値に対する実色材量の変化が滑らかとなるように、ＯＰＧＬＵＴを修正する例について説明する。以下、実施形態１〜２と共通する部分については説明を簡略化ないし省略し、以下では本実施形態に特有な点を中心に説明する。
【００８９】
（印刷システムの全体構成）
図１５は、本実施形態における印刷システム１の全体構成例を示すブロック図である。実施形態２の印刷システム１と異なる部分は、画像処理装置１００に、ＯＰＧ処理部１５０１およびＯＰＧＬＵＴ格納部１５０２が追加されている点である。ＯＰＧ処理部１５０１は、色分解処理部１０３で生成されたインク値画像データに対し、ＯＰＧ処理を行う。ＯＰＧ処理に際して、ＯＰＧ処理部１５０１は、ＯＰＧＬＵＴ格納部１５０２に格納された１次元のＯＰＧＬＵＴを参照する。
【００９０】
（ＯＰＧＬＵＴ）
図１６は、本実施形態におけるＯＰＧＬＵＴの例を示す模式図である。図１６に示される各グラフは、横軸を色分解処理後の入力色材量とし、縦軸をＯＰＧ処理後の出力色材量としたグラフであり、各グラフはＣ、Ｍ、Ｙ、Ｋ、Ｌｃ、Ｌｍのインク種類に対応している。本実施形態では、入力色材量は８ビット（２５６階調）で表され、図１６の各グラフにおけるグラフ線は、これら２５６階調の入力色材量に対応する出力色材量の推移を示している。
【００９１】
（ＯＰＧＬＵＴの修正手順）
ＯＰＧＬＵＴの修正手順は、実施形態２で説明した非線形変換ＬＵＴの修正手順と共通する部分が多い。そのため、本実施形態におけるＯＰＧＬＵＴの修正手順についても、図１４のフローチャートを参照して説明する。本実施形態において、図１４に示されるフローチャートの処理は、ＬＵＴ修正装置３００のＬＵＴ修正部３０３によって実行される。
【００９２】
Ｓ１４０１において、修正対象となるＯＰＧＬＵＴが選択される。本実施形態では、ＬＵＴ作成部３０１によって作成され、初期ＬＵＴ格納部３０２に格納されているＯＰＧＬＵＴのうち、まずＣインクのＯＰＧＬＵＴが選択される。本実施形態では、初期ＬＵＴ格納部３０２から読み出される初期ＯＰＧＬＵＴは、入力色材量と出力色材量とが等しい、線形の変換特性を有している。なお、ＯＰＧＬＵＴ格納部１５０２に既存のＯＰＧＬＵＴが格納されている場合、ＯＰＧＬＵＴ格納部１５０２から読み出したＯＰＧＬＵＴが選択されてもよい。
【００９３】
Ｓ１４０２において、Ｓ１４０１で選択されたＯＰＧＬＵＴにおける格子点の中から着目格子点ｐが選択される。ＯＰＧＬＵＴにおける始点（入力色材量＝０の格子点）および終点（入力色材量＝２５５の格子点）の位置は固定されており、本実施形態では始点に隣接する格子点が着目格子点ｐとして最初に選択される。Ｓ１４０２において、ＯＰＧＬＵＴにおける全ての格子点が順番に選択されてもよいし、一部の格子点のみを最適化処理の対象としてもよい。一部の格子点のみを最適化処理の対象とする場合、格子点間の実色材量は公知のスプライン補間法などによって算出することができる。
【００９４】
Ｓ１４０３において、パラメータｓが初期化される。本実施形態において、パラメータｓは変動幅を表し、初期化後のパラメータｓの値は「８」となる。
【００９５】
Ｓ１４０４において、初期化されたパラメータｓに基づいて、着目格子点ｐにおける出力色材量が可変される。具体的には、着目格子点ｐにおける出力色材量に対して、初期化されたパラメータｓを加算する処理と、減算する処理とを実行する。そして、着目格子点ｐについて、出力色材量が変更されていないＬ０と、出力色材量にパラメータｓが加算されたＬ１と、出力色材量からパラメータｓが減算されたＬ２と、の３種類の中間ＬＵＴが作成される。
【００９６】
Ｓ１４０５において、色分解ＬＵＴと、Ｓ１４０４で作成された３種類の中間ＬＵＴとを用いて、評価対象ラインにおける実色材量がそれぞれ算出される。Ｓ１４０５において、まず、入力ＲＧＢ信号値が色分解ＬＵＴによって実色材量に変換された後、評価対象ラインにおける実色材量が３種類の中間ＬＵＴによってそれぞれ変換される。この結果、１つの評価対象ラインについて、３種類の中間ＬＵＴ（Ｌ０〜Ｌ２）に対応する３種類の実色材量が取得される。なお、Ｓ１４０５において参照される色分解ＬＵＴは、色分解ＬＵＴ格納部１０４に格納されている色分解ＬＵＴが好適に用いられるが、実施形態１で説明した修正後の色分解ＬＵＴが参照されてもよい。なお、本実施形態では実施形態１とは異なり、色分解ＬＵＴ格納部１０４に格納されている色分解ＬＵＴ、および、修正後の色分解ＬＵＴは、１７×１７×１７の格子点を有する３次元色分解ＬＵＴである。
【００９７】
Ｓ１４０６において、評価対象ラインにおける実色材量が仮想色材量に変換される。実色材量から仮想色材量への変換処理は、図６（ａ）〜図６（ｃ）の手法および式（４−１）〜式（４−３）に従って実行することができる。Ｓ１４０６における変換処理の結果、３種類の中間ＬＵＴ（Ｌ０〜Ｌ２）に対応する３種類の仮想色材量が取得される。
【００９８】
Ｓ１４０７において、Ｓ１４０６で算出された３種類の仮想色材量それぞれについて評価値Ｅ０、Ｅ１、Ｅ２が算出される。仮想色材量に基づいて評価値Ｅ０〜Ｅ２を算出する手法では、例えば、評価対象ラインにおける仮想色材量の曲率を２乗した値の最大値が用いられる。
【００９９】
Ｓ１４０８において、評価値Ｅ０〜Ｅ２がそれぞれ比較され、３種類の中間ＬＵＴ（Ｌ０〜Ｌ２）のうち、評価値が最良となる中間ＬＵＴが選択される。そして、初期ＯＰＧＬＵＴが、この評価値が最良となる中間ＬＵＴに更新される。
【０１００】
Ｓ１４０９において、評価値が収束したか否かが判定される。本実施形態では、パラメータｓによって実色材量が可変されていないＬ０に対応する評価値Ｅ０が算出された場合（Ｓ１４０７）、評価値が収束したと判定される。評価値が収束した場合（Ｓ１４０９：ＹＥＳ）、Ｓ１４１０に移行する。評価値が収束していない場合（Ｓ１４０９：ＮＯ）、再びＳ１４０４に戻る。
【０１０１】
Ｓ１４１０において、パラメータｓの値が更新される。本実施形態では、現在のｓの値（例えば「８」）を２で割った値（例えば「４」）に更新される。
【０１０２】
Ｓ１４１１において、ｓの値が１．０未満であるか否かが判定される。ｓの値が１．０以上の場合（Ｓ１４１１：ＮＯ）、再び１４０４に戻り、ｓの値が１．０未満の場合（Ｓ１４１１：ＹＥＳ）、Ｓ１４１２に移行する。
【０１０３】
Ｓ１４１２において、全ての格子点について最適化処理が完了したか否かが判定される。全ての格子点について最適化処理が完了していない場合（Ｓ１４１２：ＮＯ）、再びＳ１４０２に戻り、全ての格子点について最適化処理が完了した場合（Ｓ１４１２：ＹＥＳ）、Ｓ１４１３に移行する。
【０１０４】
Ｓ１４１３において、全ての種類のＯＰＧＬＵＴについて修正処理が完了したか否かが判定される。全ての種類のＯＰＧＬＵＴについて修正処理が完了していない場合（Ｓ１４１３：ＮＯ）、再びＳ１４０１に戻り、全ての種類のＯＰＧＬＵＴについて修正処理が完了している場合（Ｓ１４１３：ＹＥＳ）、本フローチャートの処理を終了する。
【０１０５】
本実施形態において、仮想色材量に基づいてＯＰＧＬＵＴが修正される具体例を図１７（ａ）〜図１７（ｄ）に示す。図１７（ａ）は、図８（ａ）と同様に、色分解ＬＵＴにおけるＷ−Ｃラインの色変換特性を示す図である。図８（ａ）の色分解ＬＵＴとの相違点は、色分解ＬＵＴにおける格子点の数が１７×１７×１７に間引かれている点である。これにより、図１７（ａ）におけるＷ−Ｃラインでは、変曲点の位置と格子点の位置とが一致している。図１７（ｂ）は、図１６に示されるＯＰＧＬＵＴのうち、ＣインクおよびＬｃインクのそれぞれのＯＰＧＬＵＴを重畳して示した図である。そして、図１７（ｃ）は、ＯＰＧＬＵＴ（図１７（ｂ））によって変換された、色変換後のＣインクの実色材量およびＬｃインクの実色材量（図１７（ａ））の推移を示した図である。
図１７（ｄ）は、図１７（ｃ）における実色材量から変換された仮想色材量の一例を示す図である。図１７（ｄ）のグラフにおいて、破線部は、図１７（ｂ）に示される変換特性によって、実Ｃインクの色材量および実Ｌｃインクの色材量の線形性が崩れてしまっていることを示している。実施形態１と同様に、実色材量と仮想色材量とは相互に線形の関係が成立するため、ＯＰＧ処理によって実インクの色材量の線形性が崩れてしまった場合であっても、色材量の急峻な変化を仮想色材量の推移によって表すことができる。そして、本実施形態では仮想色材量の特性に着目し、評価値に基づいてＯＰＧＬＵＴにおける出力色材量を調整することにより、ＯＰＧＬＵＴを最適化する処理を実行するのである。図１７（ｂ）のグラフにおいて、評価値に基づいて最適化された実ＣインクのＯＰＧＬＵＴおよび実ＬｃインクのＯＰＧＬＵＴの一例が示されている。
【０１０６】
以上説明した通り、本実施形態のＬＵＴ修正装置３００によれば、仮想色材量を用いて、色分解処理以外の画像処理において参照される１次元非線形変換ＬＵＴの変換特性を評価することができる。さらに、この評価結果に基づいて、１次元非線形変換ＬＵＴの変換特性を修正することができる。
【０１０７】
（変形例）
なお、上述の実施形態は例示であり、本実施形態を以下の変形例に適用することができる。本実施形態では、入力ＲＧＢ信号値の非線形変換処理をインク色ごとに行う例を説明したが、インク色に共通の処理としてもよい。また、入力色材量から出力色材量への変換特性を示す非線形変換パラメータは、上記１次元非線形変換ＬＵＴ以外の形態であってもよく、例えば、行列演算や数式によって表現されてもよい。この場合、非線形変換パラメータの修正処理において、上記行列における係数や、数式における係数が修正される。
【０１０８】
［実施形態４］
通常、色分解処理で参照される３次元色分解ＬＵＴは、その記憶容量を節約するために、主要格子点のみに色材量を格納し、これら格子点間の色材量は補間演算により算出される。そのため、色分解ＬＵＴの記憶容量を節約するためには、色変換特性が許容される範囲内において格子点の数が削減されることが望ましい。本実施形態では、仮想色材量を用いて色分解ＬＵＴの変換特性を評価し、この評価結果に基づいて、初期色分解ＬＵＴを、適切な格子点数（グリッド数）を有する色分解ＬＵＴに修正する例について説明する。以下、実施形態１〜３と共通する部分については説明を簡略化ないし省略し、以下では本実施形態に特有な点を中心に説明する。
【０１０９】
（色分解ＬＵＴの修正手順）
本実施形態における色分解ＬＵＴの修正手順について、図１８のフローチャートを参照して説明する。図１８に示されるフローチャートの処理は、ＬＵＴ修正装置３００のＬＵＴ修正部３０３によって実行される。
【０１１０】
Ｓ１８０１において、初期色分解ＬＵＴが取得される。本実施形態では、ＬＵＴ作成部３０１によって作成され、初期ＬＵＴ格納部３０２に格納されている初期色分解ＬＵＴが読み出される。なお、初期色分解ＬＵＴは、図７〜図８で説明した手法によって作成された初期色分解ＬＵＴが読み出されるが、色分解ＬＵＴ格納部１０４に格納されている既存の色分解ＬＵＴが取得されてもよい。本実施形態では、２５６×２５６×２５６の格子点を有する、２５６グリッドの３次元色分解ＬＵＴが取得される。
【０１１１】
Ｓ１８０２において、グリッド数を減らした場合の色分解ＬＵＴにおける格子点位置が取得される。本実施形態において、Ｓ１８０２〜Ｓ１８０７のループが実行されるごとにグリッド数「１」が段階的に削減される。
【０１１２】
Ｓ１８０３において、Ｓ１８０２で取得された格子点位置における実色材量が取得される。このとき、ＬＵＴ修正部３０３は、初期色分解ＬＵＴを参照することにより、入力ＲＧＢ信号値に対応する実色材量を取得することができる。また、Ｓ１８０３において、実施形態２と同様に、評価対象ラインにおける複数の実色材量が取得される。
【０１１３】
Ｓ１８０４において、Ｓ１８０３で取得された実色材量が仮想色材量に変換される。実色材量から仮想色材量への変換処理は、図６（ａ）〜図６（ｃ）の手法および式（４−１）〜式（４−３）に従って実行することができる。
【０１１４】
Ｓ１８０５において、Ｓ１８０４で変換された仮想色材量について評価値が算出される。仮想色材量を用いて評価値を算出する手法は、実施形態２におけるＳ１４０７と同様に、評価対象ラインにおける仮想色材量の曲率を２乗した値の最大値を用いることができる。
【０１１５】
Ｓ１８０６において、Ｓ１８０５で算出された評価値を用いて、各格子点における実色材量が最適化される。具体的には、まず、評価値が小さくなるように、各格子点における仮想色材量が増減される。次いで、増減後の仮想色材量から実色材量に再び変換される。
【０１１６】
Ｓ１８０７において、Ｓ１８０５で算出された評価値が閾値以上であるか否かが判定される。評価値が閾値以上であった場合（Ｓ１８０７：ＹＥＳ）、Ｓ１８０８に移行する。評価値が閾値未満であった場合（Ｓ１８０７：ＮＯ）、再びＳ１８０２に戻る。この場合、続くＳ１８０２において、さらにグリッド数「１」を減らした場合の色分解ＬＵＴにおける格子点位置が取得される。
【０１１７】
本実施形態において、色分解ＬＵＴにおけるグリッド数（格子点数）が多いほど、色変換特性がより線形になり（例えば、図８（ａ）を参照）、Ｓ１８０５で算出される評価値が小さくなる。Ｓ１８０２において色分解ＬＵＴにおけるグリッド数（格子点数）が削減されると、実色材量の線形性が崩れ（例えば、図１１（ａ）における「格子点間引き後の実色材量」を参照）、Ｓ１８０５で算出される評価値が大きくなる。本実施形態では、許容される最大の評価値が閾値に設定され、評価値が閾値に到達した時点でＳ１８０２〜Ｓ１８０７のループを抜ける。
【０１１８】
Ｓ１８０８において、最適化された色分解ＬＵＴが出力される。最適化された色分解ＬＵＴは、初期色分解ＬＵＴから、Ｓ１８０２〜Ｓ１８０７のループが実行されるたびに削減されたグリッド数分の格子点が間引かれている。さらに、最適化された色分解ＬＵＴには、間引き後の格子点それぞれにＳ１８０６で最適化（増加または削減）された実色材量が格納されている。このように最適化された色分解ＬＵＴは、色分解ＬＵＴ格納部１０４に出力され、色分解処理部１０３によって参照される。
【０１１９】
以上説明した通り、本実施形態のＬＵＴ修正装置３００によれば、仮想色材量を用いて、グリッド数（格子点数）を減らした場合における色分解ＬＵＴの変換特性を評価することができる。さらに、この評価結果に基づいて、初期色分解ＬＵＴを、適切なグリッド数（格子点数）を有する色分解ＬＵＴに修正することができる。
【０１２０】
（変形例）
なお、上述の実施形態は例示であり、本実施形態を以下の変形例に適用することができる。本実施形態では、実施形態１と同様に色分解ＬＵＴにおける実色材量を最適化する例について説明したが、実施形態２および実施形態３と同様に、非線形変換ＬＵＴにおける格子点位置や、ＯＰＧＬＵＴを最適化処理の対象としてもよい。さらに、これらの最適化処理は組み合わされてもよい。また、色分解ＬＵＴにおけるグリッド数（格子点数）は、インク種別ごとに異ならせてもよい。
【０１２１】
［実施形態５］
実施形態１〜４では、仮想色材量を用いて色分解ＬＵＴの色変換特性を評価し、当該評価結果に基づいて色分解ＬＵＴなどを修正する例について説明した。本実施形態では、仮想色材量を用いて、入力ＲＧＢ値に対して実色材量が滑らかではない画素を検出し、この検出結果をユーザに提示する例について説明する。以下、実施形態１〜４と共通する部分については説明を簡略化ないし省略し、以下では本実施形態に特有な点を中心に説明する。
【０１２２】
（印刷システムの全体構成）
図１９は、本実施形態における印刷システム１の全体構成例を示すブロック図である。実施形態３の印刷システム１と異なる部分は、画像処理装置１００に、仮想色材画像生成部１９０１と、擬似輪郭画像生成部１９０２と、表示部１９０３と、が追加されている点である。仮想色材画像生成部１９０１は、ＯＰＧ処理部１５０１から出力される、実色材量を示す信号である実色材量信号を取得し、取得した実色材量信号を、仮想的な色材の量を示す仮想色材量信号に変換する。上記「仮想色材量」は、実施形態１〜４と同様のため詳細な説明を省略する。ＯＰＧ処理後の実色材量の一例は図１７（ｃ）に示される通りであり、実色材量信号から仮想色材量信号への変換処理は、図６（ａ）〜図６（ｃ）の手法および式（４−１）〜式（４−３）に従って実行することができる。擬似輪郭画像生成部１９０２は、カラーマッチング処理後のＲＧＢ画像信号と、仮想色材量信号とに基づいて、擬似輪郭画像信号を生成する。表示部１９０３は、ＬＵＴ修正装置３００における表示部３１５と同様に、液晶ディスプレイ（ＬＣＤ）などによって実現され、グラフィックコントローラで生成された画面を表示する。
【０１２３】
（擬似輪郭画像の生成手順）
次に、本実施形態における擬似輪郭画像の生成手順について、図２０のフローチャートを参照して説明する。図２０に示されるフローチャートの処理は、画像処理装置１００の擬似輪郭画像生成部１９０２などによって実行される。
【０１２４】
Ｓ２００１において、カラーマッチング処理部１０１からカラーマッチング処理後のＲＧＢ画像データが取得される。  
Ｓ２００２において、仮想色材画像生成部１９０１から仮想色材画像データが取得される。この仮想色材画像は仮想色材量信号によって表される。  
Ｓ２００３において、Ｓ２００１で取得されたＲＧＢ画像データと、Ｓ２００２で取得された仮想色材画像データそれぞれについて、画素毎のエッジ度が算出される。画像データにおけるエッジ度を算出する処理は、公知のラプラシアンフィルタを用いたフィルタ処理を行うことにより行われる。  
Ｓ２００４において、Ｓ２００３で算出されたエッジ度の差に基づいて、擬似輪郭画像データが生成される。本実施形態では、画素毎のエッジ度の平均と分散とが正規化された後、ＲＧＢ画像データと仮想色材画像データとにおける対応する画素におけるエッジ度の差分が算出されることにより、擬似輪郭画像データが生成される。そして、擬似輪郭画像データが表す擬似輪郭画像が表示部１９０３に表示される。
【０１２５】
図２１は、表示部１９０３に表示された擬似輪郭画像を模式的に示す図である。図２１に示される擬似輪郭画像において、斜線が付与されている領域は、ＲＧＢ画像と仮想色材画像との間でエッジ度の差分が大きい領域を示している。ＲＧＢ画像と仮想色材画像との間でエッジ度の差分が大きい領域では、色材量の急峻な変化が発生する可能性が高くなり、印刷画像において擬似輪郭などの不具合が発生しやすい。本実施形態では、擬似輪郭画像を生成し、生成された擬似輪郭画像において擬似輪郭などの不具合が発生するおそれのある領域を、ユーザが識別可能な態様で（例えば、不具合が発生するおそれのある領域の画素を着色して）表示部１９０３に表示する。ユーザは、図２１に示される擬似輪郭画像を確認することにより、例えば色変換処理において、階調性に問題が発生する可能性の高い色を認識することができる。
【０１２６】
以上説明した通り、本実施形態では、仮想色材量を用いて、入力ＲＧＢ値に対して実色材量が滑らかではない画素を検出し、この検出結果をユーザに提示する。かかる構成により、ユーザは、色変換処理において、階調性に問題が発生する可能性の高い色を認識することができる。
【０１２７】
（変形例）
なお、上述の実施形態は例示であり、本実施形態を以下の変形例に適用することができる。本実施形態では、擬似輪郭画像が表示部に表示される例を説明したが、表示される画像は擬似輪郭画像以外であってもよい。例えば、カラーマッチング後のＲＧＢ画像データが表す画像と、仮想色材画像データが表す画像とを、表示部に並置表示してもよい。
【０１２８】
［その他の実施例］
本発明は、上述の実施形態の１以上の機能を実現するプログラムを、ネットワーク又は記憶媒体を介してシステム又は装置に供給し、そのシステム又は装置のコンピュータにおける１つ以上のプロセッサーがプログラムを読出し実行する処理でも実現可能である。また、１以上の機能を実現する回路（例えば、ＡＳＩＣ）によっても実現可能である。
【符号の説明】
【０１２９】
１・・・・印刷システム
１００・・画像処理装置
１０１・・カラーマッチング処理部
１０２・・カラーマッチングＬＵＴ格納部
１０３・・色分解処理部
１０４・・色分解ＬＵＴ格納部
１０５・・ハーフトーン処理部
３００・・ＬＵＴ修正装置
３０３・・ＬＵＴ修正部
【図１】
2018191132.tif 000003
【図２】
2018191132.tif 000004
【図３】
2018191132.tif 000005
【図４】
2018191132.tif 000006
【図５】
2018191132.tif 000007
【図６】
2018191132.tif 000008
【図７】
2018191132.tif 000009
【図８】
2018191132.tif 000010
【図９】
2018191132.tif 000011
【図１０】
2018191132.tif 000012
【図１１】
2018191132.tif 000013
【図１２】
2018191132.tif 000014
【図１３】
2018191132.tif 000015
【図１４】
2018191132.tif 000016
【図１５】
2018191132.tif 000017
【図１６】
2018191132.tif 000018
【図１７】
2018191132.tif 000019
【図１８】
2018191132.tif 000020
【図１９】
2018191132.tif 000021
【図２０】
2018191132.tif 000022
【図２１】
2018191132.tif 000023