JP4645666B2

JP4645666B2 - 色処理装置およびプログラム

Info

Publication number: JP4645666B2
Application number: JP2008072534A
Authority: JP
Inventors: 信佐々木
Original assignee: Fuji Xerox Co Ltd; Fujifilm Business Innovation Corp
Current assignee: Fujifilm Business Innovation Corp
Priority date: 2008-03-19
Filing date: 2008-03-19
Publication date: 2011-03-09
Anticipated expiration: 2028-03-19
Also published as: US20090237416A1; US8395639B2; JP2009231989A

Description

本発明は、色処理装置およびプログラムに関する。

プリンタ等の画像形成装置や液晶ディスプレイ等の表示装置、スキャナ等の画像読取装置等といったカラー画像を入出力する装置では、入力した例えばＲＧＢ色空間の色信号を装置に依存しない例えばＬ^＊ａ^＊ｂ^＊色空間の色信号に変換し、さらに各装置での出力色空間である例えばＣＭＹＫ色空間の色信号に変換して、カラー画像の出力処理を行う。その際に、各装置では、例えばＣＭＹＫ色空間およびＬ^＊ａ^＊ｂ^＊色空間にて実際に測色された色データの対(実データ対)に基づきＣＭＹＫ色空間とＬ^＊ａ^＊ｂ^＊色空間との対応関係を規定する色変換特性モデルを予測して、実データには含まれない色の生成が行われる。
例えば特許文献１には、複数の入力信号の実データとそれに対応する出力信号の実データとからなるデータ対から、内挿によってデータ対の数を増やした後、入出力信号が３次元の場合は入出力の局所的な３角錘空間を線形行列で結び付け、任意の出力信号に対する入力信号を予測する方法が記載されている。

特開平２−２２６８７０号公報

ここで一般に、異なる色空間の色の対応関係を規定する色変換特性モデルは、変換元の色空間における色標本(カラーパッチ)の色データと、カラーパッチ各々について変換先の色空間にて実際に測色して得られた色データとを一対とする実データ対が複数用いられる。そのため、使用される色データ対の色空間での分散状態等といった色データ対の特性によって色変換特性モデルの色変換精度が左右されるという問題がある。
本発明は、異なる色空間の色の対応関係を規定する色変換特性モデルについて色域全域での色変換精度を向上させることを目的とする。

請求項１に記載の発明は、第１色空間の変換対象色データを取得する変換対象色データ取得部と、前記第１色空間の色と変換先である第２色空間の色との対応関係を規定する色変換特性モデルの算出に際し用いる前記第１色空間の色データと、当該色変換特性モデルの算出に際し用いる当該第２色空間の色データとを対とする色データ対を取得する色データ対取得部と、前記色データ対取得部にて取得された前記色データ対の前記第１色空間の色データと前記変換対象色データ取得部にて取得した前記変換対象色データとの色差距離に応じて当該第１色空間の色データに重み付けを行い、当該第１色空間の色データが重み付けされた当該色データ対を用いて前記色変換特性モデルを算出する色変換特性モデル算出部と、前記色データ対を構成する前記第１色空間の色データに関する当該第１色空間にて所定間隔毎に配列された基準色座標の周辺領域での分散と、当該周辺領域内の当該第１色空間の色データと対をなす前記第２色空間の色データに関する分散とを算出し、算出した当該分散の相対比により当該色データ対に関する前記色変換特性モデルの色変換精度への影響度を評価する評価値を算出する評価値算出部とを備え、前記色変換特性モデル算出部は、前記評価値算出部にて算出された前記評価値に基づいて、前記色変換特性モデルの算出に際し前記第１色空間の色データに施す前記重み付けを調整することを特徴とする色処理装置である。

請求項２に記載の発明は、前記色変換特性モデル算出部は、前記第１色空間の色データに施す前記重み付けの調整を、当該第１色空間の色データと前記変換対象色データとの色差距離と当該重み付けとの対応関係を前記評価値に基づいて変更することにより行うことを特徴とする請求項１記載の色処理装置である。

請求項３に記載の発明は、前記評価値算出部は、前記第１色空間にて所定間隔毎に配列された基準色座標の周辺領域毎に前記評価値を算出し、算出された当該評価値を当該基準色座標の一の座標成分と対応付けて前記色変換特性モデル算出部に送り、前記色変換特性モデル算出部は、当該一の座標成分と対応付けられた前記評価値に基づき前記重み付けを調整することを特徴とする請求項１記載の色処理装置である。
請求項４に記載の発明は、前記評価値算出部は、前記基準色座標の周辺領域毎に算出された複数の前記評価値に関する前記一の座標成分の座標値毎の最大値、最小値および平均値の何れかを、当該一の座標成分と対応付けられた当該評価値として規定することを特徴とする請求項３記載の色処理装置である。
請求項５に記載の発明は、前記色変換特性モデル算出部は、前記評価値算出部にて算出された前記評価値に対して前記色変換特性モデルの色変換精度が所定値を満たすように統計的に対応付けた前記重み付けを設定する前記調整を行うことを特徴とする請求項１記載の色処理装置である。

請求項６に記載の発明は、コンピュータに、第１色空間の変換対象色データを取得する機能と、前記第１色空間の色と変換先である第２色空間の色との対応関係を規定する色変換特性モデルの算出に際し用いる当該第１色空間の色データと、当該色変換特性モデルの算出に際し用いる当該第２色空間の色データとを対とする色データ対を取得する機能と、前記色データ対の前記第１色空間の色データと前記変換対象色データとの色差距離に応じて当該第１色空間の色データに重み付けを行い、当該第１色空間の色データが重み付けされた当該色データ対を用いて前記色変換特性モデルを算出する機能と、前記色データ対を構成する前記第１色空間の色データに関する当該第１色空間にて所定間隔毎に配列された基準色座標の周辺領域での分散と、当該周辺領域内の当該第１色空間の色データと対をなす前記第２色空間の色データに関する分散とを算出し、算出した当該分散の相対比により当該色データ対に関する前記色変換特性モデルの色変換精度への影響度を評価する評価値を算出する機能と、前記評価値に基づいて、前記色変換特性モデルの算出に際し前記第１色空間の色データに施す前記重み付けを調整する機能とを実現させることを特徴とするプログラムである。

請求項７に記載の発明は、前記重み付けを調整する機能は、前記第１色空間の色データに施す前記重み付けの調整を、当該第１色空間の色データと前記変換対象色データとの色差距離と当該重み付けとの対応関係を前記評価値に基づいて変更することにより行う請求項６記載のプログラムである。

請求項８に記載の発明は、前記評価値を算出する機能は、前記第１色空間にて所定間隔毎に配列された基準色座標の周辺領域毎に前記評価値を算出し、算出された当該評価値を当該基準色座標の一の座標成分と対応付けるとともに、前記重み付けを調整する機能は、当該一の座標成分と対応付けられた当該評価値に基づき前記重み付けを調整することを特徴とする請求項６記載のプログラムである。
請求項９に記載の発明は、前記評価値を算出する機能は、前記基準色座標の周辺領域毎に算出された複数の前記評価値に関する前記一の座標成分の座標値毎の最大値、最小値および平均値の何れかを、当該一の座標成分と対応付けられた当該評価値として規定することを特徴とする請求項８記載のプログラムである。
請求項１０に記載の発明は、前記重み付けを調整する機能は、算出された前記評価値に対して前記色変換特性モデルの色変換精度が所定値を満たすように統計的に対応付けた前記重み付けを設定する前記調整を行うことを特徴とする請求項６記載のプログラムである。

本発明の請求項１によれば、本発明を採用しない場合に比べて、異なる色空間の色の対応関係を規定する色変換特性モデルについて色域全域での色変換精度を向上させることができる。
本発明の請求項１によれば、色データ対に関する色変換特性モデルの色変換精度への影響度を客観的に評価でき、本発明を採用しない場合に比べて、色変換精度を向上するような色データ対に施す重み付けの設定を行うことができる。
本発明の請求項２によれば、色差距離に応じて設定される重み付けを色差距離に応じて変更することができ、本発明を採用しない場合に比べて、色データ対に評価値に対応した重み付けを設定できる。

本発明の請求項３によれば、評価値が色変換特性モデルの色変換精度に影響を及ぼし易い座標成分と対応付けられるので、本発明を採用しない場合に比べて、色変換精度の向上を図る上での色データ対に施す重み付けの設定をより効果的に行うことができる。
本発明の請求項４によれば、色空間内の領域毎に算出された複数の評価値を領域毎に対応付けることができるので、本発明を採用しない場合に比べて、評価値の色空間内の領域に関する特性を色空間内で一律に表現することができる。
本発明の請求項５によれば、評価値と重み付け値との対応関係を統計的に算出しておくことで、本発明を採用しない場合に比べて、評価値や色変換特性モデルの算出の際に生じる不確定な要素に対応して安定的な色変換特性モデルを生成できる。

本発明の請求項６によれば、本発明を採用しない場合に比べて、異なる色空間の色の対応関係を規定する色変換特性モデルについて色域全域での色変換精度を向上させることができる。
本発明の請求項６によれば、色データ対に関する色変換特性モデルの色変換精度への影響度を客観的に評価でき、本発明を採用しない場合に比べて、色変換精度を向上するような色データ対に施す重み付けの設定を行うことができる。
本発明の請求項７によれば、色差距離に応じて設定される重み付けを色差距離に応じて変更することができ、本発明を採用しない場合に比べて、色データ対に評価値に対応した重み付けを設定できる。

本発明の請求項８によれば、評価値が色変換特性モデルの色変換精度に影響を及ぼし易い座標成分と対応付けられるので、本発明を採用しない場合に比べて、色変換精度の向上を図る上での色データ対に施す重み付けの設定をより効果的に行うことができる。
本発明の請求項９によれば、色空間内の領域毎に算出された複数の評価値を領域毎に対応付けることができるので、本発明を採用しない場合に比べて、評価値の色空間内の領域に関する特性を色空間内で一律に表現することができる。
本発明の請求項１０によれば、評価値と重み付け値との対応関係を統計的に算出しておくことで、本発明を採用しない場合に比べて、評価値や色変換特性モデルの算出の際に生じる不確定な要素に対応して安定的な色変換特性モデルを生成できる。

以下、添付図面を参照して、本発明の実施の形態について詳細に説明する。
図１は、本実施の形態が適用される色処理装置１の全体構成を示すブロック図である。図１に示す色処理装置１は、格子点データの入力を受け付ける格子点データ取得部１０、格子点データを色変換処理する際に用いる色変換特性モデルを生成する色変換特性モデル生成部２０、色変換特性モデル生成部２０にて生成された色変換特性モデルを用いて格子点データを色変換処理する色変換処理部３０、格子点データ取得部１０にて受け付けた格子点データと色変換処理部３０にて色変換処理された色データとを対応付けて出力する色変換データ出力部４０、色変換特性モデル生成部２０にて色変換特性モデルを生成する際に使用する実データを評価する評価値を算出する評価値算出部の一例としての実データ評価部５０を備えている。

格子点データ取得部１０は、色変換される対象となる変換対象色データを取得する変換対象色データ取得部の一例であって、例えばカラープリンタ等のカラー画像形成装置等にて処理される色空間(例えば、ＣＭＹＫ色空間やＬ^＊ａ^＊ｂ^＊色空間)内にて規則的に配列された色データ(変換対象色データ)を取得する。この規則的に配列された色データは一般に「格子点」と呼ばれ、この格子点の色座標を本明細書では「格子点データ」と称することとする。
例えば４次元のＣＭＹＫ色空間では、一般に、色座標値Ｃ,Ｍ,Ｙ,Ｋは、それぞれ０〜１００の範囲で設定される。そこで、各色座標軸における色座標値をそれぞれ任意のステップ(例えば、１０ステップ)で区切り、区切られた色座標値のすべての組み合わせを作成して、これをＣＭＹＫ色空間での格子点データとする。例えば１０ステップで色座標値を区切った場合には、(０,０,０,０)、(０,１０,０,０)、(０,０,１０,０)、(０,０,０,１０)、・・・が格子点データとなる。
また、例えば３次元のＬ^＊ａ^＊ｂ^＊色空間では、一般に、明度Ｌ^＊は０〜１００、色度ａ^＊, ｂ^＊はそれぞれ−３００〜３００の範囲で設定される。そこで、同様に、各色座標軸における色座標値をそれぞれ任意のステップで区切った色座標値のすべての組み合わせが格子点データとなる。
格子点データ取得部１０では、例えばこのようなＣＭＹＫ色空間やＬ^＊ａ^＊ｂ^＊色空間内の格子点の色座標が格子点データとして入力される。そして、取得された格子点データは色変換特性モデル生成部２０と色変換処理部３０とに出力される。

色変換特性モデル生成部２０は、格子点データ取得部１０から入力される格子点データの色空間(例えばＣＭＹＫ色空間)内で種々の色座標値が組み合わされて生成された色標本(カラーパッチ)の色データと、生成されたカラーパッチを変換先の色空間(例えばＬ^＊ａ^＊ｂ^＊色空間)において実際に測定(測色)して得られた色データとの対である実データ対(色データ対)を用いて、色変換処理部３０にて色変換処理を行う際に用いられる色変換特性モデルを生成する。
ここで、「色変換特性モデル」とは、例えば第１の色空間の色信号(第１色信号)を第２の色空間の色信号(第２色信号)に変換する際に、第１色信号と第２色信号との対応関係を定式化して規定するものをいう。本実施の形態の色変換特性モデル生成部２０では、実データ評価部５０での実データ対に関する評価結果に基づいて、色変換特性モデルを生成する際に実データ対に対して施すパラメータを調整する。この点に関しては、後段で詳述する。

色変換処理部３０は、色変換特性モデル生成部２０にて生成された色変換特性モデルに基づいて、格子点データ取得部１０から取得した格子点データに対する色変換処理を行う。すなわち、格子点データ取得部１０から取得した例えばＣＭＹＫ色空間の格子点データを、色変換特性モデル生成部２０にて生成された色変換特性モデル(正変換)に基づいてデバイスに依存しない色空間(デバイス非依存色空間)である例えばＬ^＊ａ^＊ｂ^＊色空間の色データに変換する。また、格子点データ取得部１０から取得した例えばＬ^＊ａ^＊ｂ^＊色空間の格子点データを、色変換特性モデル生成部２０にて生成された色変換特性モデル(逆変換)に基づいてデバイスに依存する色空間(デバイス依存色空間)である例えばＣＭＹＫ色空間の色データに変換する。そして、色変換処理された色データを色変換データ出力部４０に出力する。

色変換データ出力部４０は、格子点データ取得部１０にて受け付けた格子点データ(第１色信号)と、格子点データ取得部１０にて受け付けた格子点データ各々について色変換処理部３０にて色変換処理された色データ(第２色信号)とを対応付ける。そして、対応付けられた格子点データと色変換処理された色データとを一対とするデータ群を出力する。このデータ群は、例えばカラープリンタ等のカラー画像形成装置等にてパーソナルコンピュータ(ＰＣ)等の外部機器から入力された画像データを色変換処理する際に用いられるＤＬＵＴ(Direct Look-Up Table)と呼ばれる色変換テーブルの生成に使用される。

次に、本実施の形態の実データ評価部５０について説明する。
上記したように、色変換特性モデル生成部２０では、ＣＭＹＫ色空間で生成されたカラーパッチの色データと、カラーパッチをＬ^＊ａ^＊ｂ^＊色空間において測色して得られた色データとの対である実データ対を用いて色変換特性モデルを生成する。その際に、色変換特性モデルの色変換精度は、生成の際に使用される実データ対のＣＭＹＫ色空間やＬ^＊ａ^＊ｂ^＊色空間での分散状態に応じて定まるという特性を有している。そこで、実データ評価部５０では、使用される実データ対のＣＭＹＫ色空間やＬ^＊ａ^＊ｂ^＊色空間での統計学上の分散(以下、単に「分散」という)を求める。そして、求められた実データ対の分散に基づいて、実データ対により生成される色変換特性モデルの色変換精度を評価する評価値を算出する。
実データ評価部５０にて算出された評価値は、色変換特性モデル生成部２０に送られる。そして、色変換特性モデル生成部２０にて評価値に基づき、色変換特性モデルを生成する際に実データ対に対して作用させるパラメータの設定値(例えば加重)の調整が行われる。

図２は、実データ評価部５０の機能構成を示すブロック図である。図２に示したように、実データ評価部５０は、実データ対取得部５１、実データ対記憶部５２、評価用格子点生成部５３、実データ対評価値算出部５４、評価値出力部５５を備えている。
実データ対取得部５１は、色変換特性モデル生成部２０にて色変換特性モデルを生成する際に使用する実データ対(色データの対)を取得する。
実データ対記憶部５２は、実データ対取得部５１にて取得された実データ対を記憶する。すなわち、色標本(カラーパッチ)に関するＣＭＹＫ色空間(第１色空間)で表された実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)と、それに対応するＬ^＊ａ^＊ｂ^＊色空間(第２色空間)にて測色して得られた実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)とが一組となったｎ組の実データ対を記憶する。ここで、ｎ,ｉは整数であってｉ＝１〜ｎとする。

評価用格子点生成部５３は、実データ対を色変換特性モデルの色変換精度への影響度の観点から評価する際の色空間内での評価基準位置となる基準色座標を生成する。すなわち、評価用格子点生成部５３は、色変換特性モデルの生成時に使用する実データ対を評価する際に用いるＣＭＹＫ色空間(第１色空間)での評価用格子点データ(Ｃ_Ｌ,Ｍ_Ｌ,Ｙ_Ｌ,Ｋ_Ｌ)を生成する。評価用格子点生成部５３にて生成する評価用格子点データは、実データ対を評価する際に用いるものであるため、通常、格子点データ取得部１０から入力される格子点データと同等数を用意する必要はない。そこで、評価用格子点生成部５３は、例えば格子点データ取得部１０から入力される格子点データよりも粗いステップ(例えば、２０ステップ)で色座標値を区切って評価用格子点データを生成する。なお、例えば色空間内の特定の領域において実データ対の詳細な評価を行う場合には、その特定の領域に関しては格子点データ取得部１０から入力される格子点データよりも細かいステップ(例えば、５ステップ)で色座標値を区切って評価用格子点データを生成してもよい。また、格子点データ取得部１０から入力される格子点データと同等のステップ(例えば、１０ステップ)で評価用格子点データを生成してもよい。
すなわち、評価用格子点生成部５３にて生成される評価用格子点データ(Ｃ_Ｌ,Ｍ_Ｌ,Ｙ_Ｌ,Ｋ_Ｌ)は、色変換特性モデルを生成する際に用いる実データ対が色変換特性モデルに対して如何なる誤差を生じさせるかの目安(評価値)を算出するための色空間内での評価位置となるものである。

実データ対評価値算出部５４は、評価用格子点生成部５３からＣＭＹＫ色空間(第１色空間)の評価用格子点データ(Ｃ_Ｌ,Ｍ_Ｌ,Ｙ_Ｌ,Ｋ_Ｌ)を取得する。さらに、評価用格子点生成部５３にて生成されたＣＭＹＫ色空間での評価用格子点の周囲に位置するＣＭＹＫ色空間の実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)と、それに対応するＬ^＊ａ^＊ｂ^＊色空間(第２色空間)の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)との対である実データ対を実データ対記憶部５２から読み出す。そして、読み出した実データ対に関して、ＣＭＹＫ色空間における実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)の分散と、Ｌ^＊ａ^＊ｂ^＊色空間における実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)の分散とを算出する。

ここで、各色空間での実データの分散と、実データ対を用いて生成される色変換特性モデルの色変換精度との関係について述べる。
まず図３は、色変換特性モデルの機能を説明する図である。図３に示したように、色変換特性モデルはＣＭＹＫ色空間(第１色空間)とＬ^＊ａ^＊ｂ^＊色空間(第２色空間)との間の写像を規定する関数であって、ＣＭＹＫ色空間の色データ(変換対象色)は、色変換特性モデルによりＬ^＊ａ^＊ｂ^＊色空間の色データ(変換色)に変換される。なお、第１色空間であるＣＭＹＫ色空間から第２色空間であるＬ^＊ａ^＊ｂ^＊色空間への色変換を正方向として、この場合の色変換特性モデルを「正変換の色変換特性モデル」とする。一方、Ｌ^＊ａ^＊ｂ^＊色空間の色データ(変換色)は、第２色空間から第１色空間への「逆変換の色変換特性モデル」によりＣＭＹＫ色空間の色データ(変換対象色)に変換される。
ここで、正変換の色変換特性モデルの生成には、主として変換対象色の周囲に位置するＣＭＹＫ色空間の実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)と、これらに対応するＬ^＊ａ^＊ｂ^＊色空間の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)とが用いられることとなる。一方、逆変換の色変換特性モデルの生成には、主として変換色の周囲に位置するＬ^＊ａ^＊ｂ^＊色空間の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)と、これらに対応するＣＭＹＫ色空間の実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)とが用いられる。

続いて図４は、正変換により変換された変換色が、逆変換によって元のＣＭＹＫ色空間の色座標に変換される際に生じる誤差を説明する図である。図４に示したように、ＣＭＹＫ色空間の変換対象色が正変換の色変換特性モデルによりＬ^＊ａ^＊ｂ^＊色空間の変換色に変換された場合に、その逆変換の色変換特性モデルによりＬ^＊ａ^＊ｂ^＊色空間の変換色が元のＣＭＹＫ色空間の色座標に変換されると、逆変換された色座標は、元の変換対象色と完全に一致せず、元の変換対象色との間に誤差が生じる場合がある。
この誤差は、正変換の色変換特性モデルと逆変換の色変換特性モデルとが生成される際に使用される実データがそれぞれ異なることに起因するものであって、モデル間の相違により生じる「モデル誤差」と呼ばれるものである。すなわち、上記したように、正変換の色変換特性モデルでは、主として変換対象色の周囲に位置するＣＭＹＫ色空間の実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)と、これらに対応するＬ^＊ａ^＊ｂ^＊色空間の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)とが用いられる。それに対して、逆変換の色変換特性モデルでは、主として変換色の周囲に位置するＬ^＊ａ^＊ｂ^＊色空間の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)と、これらに対応するＣＭＹＫ色空間の実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)とが用いられる。そのため、以下に述べるように、ＣＭＹＫ色空間の実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)に対応するＬ^＊ａ^＊ｂ^＊色空間の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)がＬ^＊ａ^＊ｂ^＊色空間内で偏在する場合には、逆変換の色変換特性モデルの色変換精度が低下し、モデル誤差が大きくなる。

図５および図６は、ＣＭＹＫ色空間の実データおよびこれらに対応するＬ^＊ａ^＊ｂ^＊色空間の実データに関するぞれぞれの色空間内での分散状態を表す図である。ＣＭＹＫ色空間にて種々の色座標値が組み合わされて生成されたカラーパッチは、個々の色座標(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)がＣＭＹＫ色空間内で適度に分散するように設定される。それに対応して、図５に示したように、カラーパッチをＬ^＊ａ^＊ｂ^＊色空間において測色して得られた色座標(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)もＬ^＊ａ^＊ｂ^＊色空間にて適度に分散した状態であれば、個々の色座標(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)はＬ^＊ａ^＊ｂ^＊色空間内でそれぞれ適度な色差距離を持った位置に存在することとなる。
ところが、カラーパッチのＬ^＊ａ^＊ｂ^＊色空間での測色値である色座標(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)が、図６(ａ)に示したような小さな領域に集中して分布する場合や、図６(ｂ)に示したような偏った領域に分布する場合等がある。このような分散状態では、個々の色座標(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)はＬ^＊ａ^＊ｂ^＊色空間内でそれぞれ色差距離の小さい近接した位置に存在することとなる。このような状態は、Ｋ値の大きなカラーパッチにて生じ易くなる。

ここで図７は、Ｌ^＊ａ^＊ｂ^＊色空間での実データの分散状態と、実データ対を用いて生成される色変換特性モデルの色変換精度との関係を説明する図である。図７に示した領域Ａは、図５に相当する領域であって、ＣＭＹＫ色空間の実データの分散ｘ１に対応するＬ^＊ａ^＊ｂ^＊色空間での実データの分散ｙ１は、比較的大きい。それにより、この領域Ａで生成される色変換特性モデルは適度な傾き(変化率)を有することとなり、ＣＭＹＫ色空間の色信号からＬ^＊ａ^＊ｂ^＊色空間の色信号への変換誤差は小さなものとなる。
これに対して、図７に示した領域Ｂは、図６(ａ)や(ｂ)に相当する領域であって、ＣＭＹＫ色空間の実データの分散ｘ２に対応するＬ^＊ａ^＊ｂ^＊色空間での実データの分散ｙ２は、比較的小さい。それにより、この領域Ｂで生成される色変換特性モデルは傾き(変化率)が小さくなり、ＣＭＹＫ色空間の色信号からＬ^＊ａ^＊ｂ^＊色空間の色信号への変換誤差は大きなものとなる。また、Ｌ^＊ａ^＊ｂ^＊色空間での色データの測定誤差が色変換特性モデルの色変換精度に与える影響も大きくなる。
そのため、例えば図７の領域Ｂにおいて色変換特性モデルを生成するに際して、領域Ａと同様の手法で(例えば、領域Ａと同様のパラメータ値を用いて)色変換特性モデルを生成するとした場合に、領域Ｂの色変換特性モデルの色変換精度は低くなることが予測される。

そこで、実データ対評価値算出部５４では、ＣＭＹＫ色空間(第１色空間)の実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)の分散と、Ｌ^＊ａ^＊ｂ^＊色空間(第２色空間)の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)の分散とを求める。さらに、ＣＭＹＫ色空間の実データの分散に対するＬ^＊ａ^＊ｂ^＊色空間の実データの分散の相対比(相対分散値)を算出する。そして、算出された相対分散値を実データに関する評価値とする。それにより、算出した相対分散値(評価値)を参照することで、実データ対により生成される色変換特性モデルの色変換精度が予測される。
すなわち、図７に示したように、例えば領域Ａのような相対分散値(ｙ１/ｘ１)が高い領域では、生成される色変換特性モデルは適度な傾き(変化率)を有する。すなわち、ＣＭＹＫ色空間の色データとＬ^＊ａ^＊ｂ^＊色空間の色データとの相関性は高い。そのため、色変換特性モデルにおける領域Ａでは、ＣＭＹＫ色空間の色データとＬ^＊ａ^＊ｂ^＊色空間の色データとの間の変換誤差は小さく、この領域Ａでの色変換特性モデルの色変換精度は高いものとなる。その一方で、例えば領域Ｂのような相対分散値(ｙ２/ｘ２)が低い領域では、生成される色変換特性モデルの傾き(変化率)は小さい。すなわち、ＣＭＹＫ色空間の色データとＬ^＊ａ^＊ｂ^＊色空間の色データとの相関性は低い。そのため、色変換特性モデルにおける領域Ｂでは、ＣＭＹＫ色空間の色データとＬ^＊ａ^＊ｂ^＊色空間の色データとの間に変換誤差が生じ易く、この領域Ｂでの色変換特性モデルの色変換精度は低いものとなる。
このように、ＣＭＹＫ色空間(第１色空間)の実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)の分散と、Ｌ^＊ａ^＊ｂ^＊色空間(第２色空間)の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)の分散との相対分散値(評価値)を算出することで、生成される色変換特性モデルの色変換精度が、色空間内での領域毎に予測される。

具体的には、実データ対評価値算出部５４は、次のような処理を行って相対分散値(評価値)を算出する。
図８は、実データ対評価値算出部５４にて相対分散値を算出する際の処理の手順の一例を示したフローチャートである。実データ対評価値算出部５４は、まず、評価用格子点生成部５３から評価用格子点データ(Ｃ_Ｌ,Ｍ_Ｌ,Ｙ_Ｌ,Ｋ_Ｌ)を取得する(ステップ１０１)。また、実データ対記憶部５２から色変換特性モデルを生成する際に用いる実データ対を取得する(ステップ１０２)。
そして実データ対評価値算出部５４は、実データ対記憶部５２から取得したＣＭＹＫ色空間の実データの中から、評価用格子点生成部５３から取得した一の評価用格子点の周囲の所定の領域Ｆに存在する実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)を抽出する(ステップ１０３)。例えば、評価用格子点を基点とした所定の色差距離の領域Ｆの内側に位置する実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)を抽出する。なお、評価用格子点の周囲の領域Ｆを規定する場合において、評価用格子点を基点として設定する色差距離は、各座標成分Ｃ,Ｍ,Ｙ,Ｋ毎に異なる距離を設定してもよいし、同じ距離を設定してもよい。
またステップ１０３においては、評価用格子点と実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)との色差距離を算出し、算出された色差距離の小さなものから順に所定個数を抽出してもよい。

ここで図９は、ＣＭＹＫ色空間にて評価用格子点の周囲の領域Ｆ内に存在する実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)と、それらに対応するＬ^＊ａ^＊ｂ^＊色空間の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)とを説明する図である。図９に示したように、ステップ１０３では、評価用格子点を基点とした領域Ｆの内側に位置する実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)が抽出される。

次に、実データ対評価値算出部５４は、ステップ１０３にて抽出されたＣＭＹＫ色空間の領域Ｆ内の実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)に関する分散σ^２ _ＣＭＹＫを算出する(ステップ１０４)。すなわち、領域Ｆの内側に位置する実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)について、実データ対評価値算出部５４は次の(１)式により分散σ^２ _ＣＭＹＫを算出する。ただし、(１)式において、ｘ_ｉ＝ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ、ｘ＝ｘ_１,ｘ_２,ｘ_３,ｘ_４である。また、ｍｘは、次の(２)式に示した領域Ｆの内側に位置する実データの平均値、Ｎは領域Ｆの内側に位置する実データの個数である。(２)式により算出された実データの平均値により構成される色座標(ｍｘ_１,ｍｘ_２,ｍｘ_３,ｍｘ_４)は、図９に示した平均値座標となる。
このステップ１０４にて算出される分散σ^２ _ＣＭＹＫは、ＣＭＹＫ色空間の評価用格子点の周辺領域Ｆにおける平均値座標からの実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)の散らばり度合いを示している。

続いて、実データ対評価値算出部５４は、ステップ１０３にて抽出されたＣＭＹＫ色空間の領域Ｆ内の実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)各々に対応する、Ｌ^＊ａ^＊ｂ^＊色空間の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)を抽出する(ステップ１０５)。ステップ１０５にて抽出される実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)は、図９に示したＬ^＊ａ^＊ｂ^＊色空間の領域Ｇ内に位置するものであるとする。

そして実データ対評価値算出部５４は、Ｌ^＊ａ^＊ｂ^＊色空間の領域Ｇ内の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)に関する分散σ^２ _Ｌ*ａ*ｂ*を算出する(ステップ１０６)。すなわち、Ｌ^＊ａ^＊ｂ^＊色空間の領域Ｇ内に位置する実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)について、実データ対評価値算出部５４は次の(３)式により分散σ^２ _Ｌ*ａ*ｂ*を算出する。ただし、(３)式において、ｙ_ｉ＝ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ、ｙ＝ｙ_１,ｙ_２,ｙ_３である。また、ｍｙは、次の(４)式に示した領域Ｇの内側に位置する実データの平均値、Ｎは領域Ｇの内側に位置する実データの個数であり、領域Ｆの内側に位置する実データの個数と同数となる。(４)式により算出された実データの平均値により構成される色座標(ｍｙ_１ｉ,ｍｙ_２ｉ,ｍｙ_３ｉ)は、図９に示した平均値座標となる。
このステップ１０６にて算出される分散σ^２ _Ｌ*ａ*ｂ*は、ＣＭＹＫ色空間の評価用格子点の周辺領域Ｆに対応するＬ^＊ａ^＊ｂ^＊色空間の領域Ｇにおける、平均値座標からの実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)の散らばり度合いを示している。

次いで実データ対評価値算出部５４は、ステップ１０４にて算出された分散σ^２ _ＣＭＹＫとステップ１０６にて算出された分散σ^２ _Ｌ*ａ*ｂ*とから、ＣＭＹＫ色空間の実データの分散σ^２ _ＣＭＹＫに対するＬ^＊ａ^＊ｂ^＊色空間の実データの分散σ^２ _Ｌ*ａ*ｂ*の相対比(相対分散値)を算出する。そして、算出された相対分散値を、色変換特性モデルを生成する際に用いる実データ対に関する評価値Ｓとする。すなわち、実データ対評価値算出部５４は、次の(５)式により実データ対に関する評価値Ｓを算出する(ステップ１０７)。
この評価値Ｓは、上記したように、色変換特性モデルを生成する際に用いる実データ対が色変換特性モデルに対して如何なる誤差を生じさせるかの目安となる。

そして、実データ対評価値算出部５４は、ステップ１０７にて算出した評価値Ｓを、評価値Ｓが算出された評価用格子点データ(Ｃ_Ｌ,Ｍ_Ｌ,Ｙ_Ｌ,Ｋ_Ｌ)と対応付けて評価値出力部５５に送る。また、ステップ１０７にて評価値Ｓを算出した後は、ステップ１０１にて取得した評価用格子点の中から次の評価用格子点を選択して、次の評価用格子点の周辺の領域Ｆにて同様に評価値Ｓを算出する。そして、すべての評価用格子点の周辺の領域Ｆの評価値Ｓを算出する。
その後、評価値出力部５５は、評価用格子点データ(Ｃ_Ｌ,Ｍ_Ｌ,Ｙ_Ｌ,Ｋ_Ｌ)と対応付けられた評価値Ｓを色変換特性モデル生成部２０に出力する。評価値Ｓを受け取った色変換特性モデル生成部２０では、評価値Ｓに基づき、色変換特性モデルを生成する際に実データ対に対して施すパラメータの設定値(加重)を、変換対象色である格子点データが位置する色空間内での領域(評価用格子点の周辺領域Ｆ)毎に調整する。
すなわち、実データ評価部５０にて算出された評価値Ｓを用いれば、生成された色変換特性モデルにより生じるモデル誤差を色空間内での領域毎に予測できる。このことから、色変換特性モデルの生成の際に使用される実データ対に関して言えば、この実データ対により生成される色変換特性モデルの色変換精度への影響度の観点からの当該実データ対の評価が色空間内での領域毎に行われる。

ここで図１０および図１１は、評価値出力部５５から色変換特性モデル生成部２０に出力される評価用格子点データ(Ｃ_Ｌ,Ｍ_Ｌ,Ｙ_Ｌ,Ｋ_Ｌ)と対応付けられた評価値Ｓの一例を示した図である。図１０および図１１では、評価用格子点データ(Ｃ_Ｌ,Ｍ_Ｌ,Ｙ_Ｌ,Ｋ_Ｌ)のＫ座標値(Ｋ_Ｌ)によりＣＭＹＫ色空間での領域を区分けし、同一のＫ座標値(Ｋ_Ｌ)を有するＣＭＹＫ色空間での領域と評価値Ｓとを関係を示している。また、図１０は画像形成装置(プリンタ)Ａでの評価値Ｓ、図１１は画像形成装置Ａとは異なる色変換特性を有する画像形成装置(印刷機)Ｂでの評価値Ｓをそれぞれ表している。
実データ評価部５０では、図１０および図１１に示したように、それぞれの画像形成装置に備えられた色処理装置１において色変換特性モデルの生成時に用いられる実データを例えばＫ座標値(Ｋ_Ｌ)により区分けされたＣＭＹＫ色空間での領域毎に生成する。そして、ＣＭＹＫ色空間でのＫ座標値(Ｋ_Ｌ)であるカバレッジ値(Ｋ色の画像面積率％)と評価値Ｓとの関係を評価値特性として表し、色変換特性モデル生成部２０に出力する。

なお、本実施の形態の実データ評価部５０においては、各Ｋ座標値(Ｋ_Ｌ)における評価値Ｓの最小値を用いて評価値特性を表している。その他として、評価値特性を各Ｋ座標値(Ｋ_Ｌ)における評価値Ｓの最大値を用いて表してもよい。さらには、各Ｋ座標値(Ｋ_Ｌ)における評価値Ｓの平均値を用いて表してもよい。すなわち、各Ｋ座標値(Ｋ_Ｌ)における評価値Ｓを客観的に表現する方法であれば、何れの方法によって各Ｋ座標値(Ｋ_Ｌ)での評価値Ｓを表してもよい。

次の図１２は、各種の画像形成装置の評価値Ｓに関する評価値特性を示した図である。画像形成装置においては、使用される実データ対が異なる。また、実データを測色する測定器も異なる。さらには、装置構成や動作設定条件、使用される現像剤の特性等もそれぞれ異なっている。そのため、図１２に示したように、画像形成装置毎に評価値特性は異なる。そこで、各画像形成装置に備えられた色処理装置１にて色変換特性モデルを生成する際には、それぞれ使用する実データ対の評価値特性に応じて色変換特性モデルを生成することとなる。

なお、本実施の形態の実データ評価部５０では、ＣＭＹＫ色空間の実データの分散σ^２ _ＣＭＹＫとＬ^＊ａ^＊ｂ^＊色空間の実データの分散σ^２ _Ｌ*ａ*ｂ*とを用いて評価値Ｓを算出したが、これに限定されず、実データ対により生成される色変換特性モデルの色変換精度を評価することができれば、他のファクターを用いて評価値Ｓしてもよい。

次に、本実施の形態の色変換特性モデル生成部２０について詳細に説明する。
図１３は、色変換特性モデル生成部２０の機能構成を示すブロック図である。色変換特性モデル生成部２０は、実データ対入力部２１、実データ対記憶部２２、加重生成部２３、算出部２４を備えている。
実データ対入力部２１は、色データ対取得部の一例であって、格子点データ取得部１０から入力される格子点データの色空間(第１色空間、ＣＭＹＫ色空間)において種々の色座標値が組み合わされて生成された色標本(カラーパッチ)の色データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)と、カラーパッチを変換先の色空間(第２色空間、Ｌ^＊ａ^＊ｂ^＊色空間)において実際に測定(測色)して得られた色データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)との対である実データ対(色データ対)を取得する。
実データ対記憶部２２は、実データ対入力部２１から入力された実データ対を記憶する。すなわち、カラーパッチに関するＣＭＹＫ色空間で表された実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)と、それに対応するＬ^＊ａ^＊ｂ^＊色空間の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)とが対となったｎ組の実データ対を記憶する。ここで、ｎ,ｉは整数であってｉ＝１〜ｎとする。上記した実データ評価部５０の実データ対取得部５１は、この実データ対記憶部２２から実データ対を取得する。

加重生成部２３は、色変換特性モデル算出部の一部を構成する機能部の一例であって、実データを用いて例えば回帰分析等の統計処理により色変換特性モデルを算出する際に、実データに対して施す加重(重み)を生成する。加重生成部２３では、後段で説明するように、逆変換の色変換特性モデルを生成する場合に、実データ評価部５０から取得した評価値Ｓに基づき、実データ対への重み付けを実データ対に対して施す加重によってＬ^＊ａ^＊ｂ^＊色空間内での領域毎に調整する。
算出部２４は、色変換特性モデル算出部の一部を構成する機能部の一例であって、例えば回帰分析等の統計処理によりＣＭＹＫ色空間の色データとＬ^＊ａ^＊ｂ^＊色空間の色データとの対応関係を規定する正変換および逆変換両方向の色変換特性モデルを算出する。すなわち、算出部２４は、格子点データ取得部１０から格子点データと、実データ対記憶部２２から実データ対とを取得する。また、加重生成部２３から格子点データとの色差距離に応じた加重を取得する。そして算出部２４は、実データ対記憶部２２から取得した実データに対して、加重生成部２３にて生成された加重を施し、加重が施された実データを用いて線形回帰分析を行う。それにより、第１色空間の色(格子点データ)と、取得した第１色空間の色に対応する第２色空間の色(変換色データ)との対応関係を規定する正変換の色変換特性モデルを算出する。また、第２色空間の色(格子点データ)と、取得した第２色空間の色に対応する第１色空間の色(変換色データ)との対応関係を規定する逆変換の色変換特性モデルを算出する。そして、算出された正逆両方向の色変換特性モデルを色変換処理部３０に送信する。

次に、色変換特性モデル生成部２０にて行われる色変換特性モデルの生成について説明する。まず、第１色信号として４次元であるＣＭＹＫ色空間での色データ(Ｃ,Ｍ,Ｙ,Ｋ)から、第２色信号として３次元であるＬ^＊ａ^＊ｂ^＊色空間での色データ(Ｌ^＊,ａ^＊,ｂ^＊)への正変換の色変換特性モデルを生成する場合について述べる。

格子点データ取得部１０にて入力される格子点データがＣＭＹＫ色空間(第１色空間)の４次元データであり、色変換処理部３０にて色変換処理されて出力される変換色データがＬ^＊ａ^＊ｂ^＊色空間(第２色空間)の３次元データである場合には、次の(６)式に示したように、ｎ組のＣＭＹＫ色空間の実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)と、それに対応するＬ^＊ａ^＊ｂ^＊色空間での出力変換値(ｙ´_１ｉ,ｙ´_２ｉ,ｙ´_３ｉ)とは、それぞれが定数項を含む線形の関係を構成する、ｍ_１１,ｍ_１２,ｍ_１３,ｍ_１４,…を成分とした(７)式の行列Ｍ_ｉにより結び付けられる。
そこで、色変換特性モデル生成部２０では、(６)式により算出される出力変換値(ｙ´_１ｉ,ｙ´_２ｉ,ｙ´_３ｉ)が、ＣＭＹＫ色空間の実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)と対をなすＬ^＊ａ^＊ｂ^＊色空間の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)に極力近似するような(７)式の行列Ｍ_ｉを正変換の色変換特性モデルとして算出する。

図１４は、色変換特性モデル生成部２０にて色変換特性モデルを生成する際の処理の手順の一例を示したフローチャートである。図１４に示したように、色変換特性モデル生成部２０における色変換特性モデルの生成処理では、初期加重生成工程(ステップ２０１)、行列成分算出工程(ステップ２０２)、変換値算出工程(ステップ２０３)、加重再生成工程(ステップ２０４)、行列成分再算出工程(ステップ２０５)、変換値算出工程(ステップ２０６)が実行される。

[初期加重生成工程]
色変換特性モデル生成部２０においては、まず、加重生成部２３にて初期加重Ｗ´_ｉｊを生成する初期加重生成工程(ステップ２０１)が行われる。加重生成部２３は、格子点データ取得部１０から変換対象色である格子点データ(ｘ_１ｊ,ｘ_２ｊ,ｘ_３ｊ,ｘ_４ｊ)を取得し、また、実データ対記憶部２２からＣＭＹＫ色空間の実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)を取得する。そして、格子点データ(ｘ_１ｊ,ｘ_２ｊ,ｘ_３ｊ,ｘ_４ｊ)と実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)との信号成分毎の差分を規格化し、この規格化された差分成分からなる規格化信号に関するユークリッド距離に応じて単調減少する、後段の(１０)式に示した第１の関数Ｆ_ｉｊで設定される初期加重Ｗ´_ｉｊを生成する。なお、(１０)式のｘ_１０,ｘ_２０,ｘ_３０,ｘ_４０は規格化定数であって、ユークリッド距離に応じて初期加重Ｗ´_ｉｊの幅を規定する。したがって、規格化定数ｘ_１０,ｘ_２０,ｘ_３０,ｘ_４０は、実データに施す加重を調整するパラメータの一例である。

[行列成分算出工程]
初期加重生成工程(ステップ２０１)に続く行列成分算出工程(ステップ２０２)では、算出部２４により、(６)式の行列Ｍ_ｉの各成分ｍ_１１,ｍ_１２,ｍ_１３,ｍ_１４,…を求める処理を行う。
算出部２４は、格子点データ取得部１０から格子点データ(ｘ_１ｊ,ｘ_２ｊ,ｘ_３ｊ,ｘ_４ｊ) を取得する。また、実データ対記憶部２２からＣＭＹＫ色空間の実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)およびそれに対応するＬ^＊ａ^＊ｂ^＊色空間の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)を取得する。さらに、加重生成部２３から初期加重Ｗ´_ｉｊを取得する。
そして、算出部２４は、(６)式にて算出される出力変換値(ｙ´_１ｉ,ｙ´_２ｉ,ｙ´_３ｉ)と、(６)式にて出力変換値(ｙ´_１ｉ,ｙ´_２ｉ,ｙ´_３ｉ)を算出する際に用いたＣＭＹＫ色空間の実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)と対をなすＬ^＊ａ^＊ｂ^＊色空間の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)との間のユークリッド距離に、初期加重生成工程(ステップ２０１)にて生成された初期加重Ｗ´_ｉｊを乗じた２乗和Ｅ_ｊ(次の(８)式)を最小にするという条件の下で、線形回帰分析の一例である最小２乗法により、(６)式の行列Ｍ_ｉの各成分ｍ_１１,ｍ_１２,ｍ_１３,ｍ_１４,…を求める。なお、行列成分算出工程(ステップ２０２)においては、(８)式の加重Ｗ_ｉｊが初期加重Ｗ´_ｉｊに置き換えられて算出される。

[変換値算出工程]
次の変換値算出工程(ステップ２０３)では、算出部２４が、ステップ２０２の行列成分算出工程にて算出された行列Ｍ_ｉの各成分ｍ_１１,ｍ_１２,ｍ_１３,ｍ_１４,…を用いて、次の(９)式により、格子点データ取得部１０にて入力された格子点データ(ｘ_１ｊ,ｘ_２ｊ,ｘ_３ｊ,ｘ_４ｊ)から変換色データ(ｙ_１ｊ,ｙ_２ｊ,ｙ_３ｊ)を算出する。

[加重再生成工程]
次の加重再生成工程(ステップ２０４)では、加重生成部２３により加重Ｗ_ｉｊが再度生成される。ここでは、ステップ２０１の初期加重生成工程にて初期加重Ｗ´_ｉｊを生成する際に用いた第１の関数Ｆ_ｉｊ((１０)式)に加えて、第２の関数Ｇ_ｉｊ((１１)式)を用いる。
第１の関数Ｆ_ｉｊでは、初期加重生成工程(ステップ２０１)と同様に、変換対象である格子点データ(ｘ_１ｊ,ｘ_２ｊ,ｘ_３ｊ,ｘ_４ｊ)と実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)との信号成分毎の差分が規格化される。そして、第１の関数Ｆ_ｉｊにより、この規格化された差分成分からなる規格化信号に関するユークリッド距離に応じて単調減少する第１の加重Ｗ_１ｉｊを設定する。なお、加重再生成工程(ステップ２０４) にて用いる(１０)式のｘ_１０,ｘ_２０,ｘ_３０,ｘ_４０は規格化定数であって、ユークリッド距離に応じて第１の加重Ｗ_１ｉｊの幅を規定する。すなわち同様に、ここでの規格化定数ｘ_１０,ｘ_２０,ｘ_３０,ｘ_４０も、実データに施す加重を調整するパラメータの一例である。
第２の関数Ｇ_ｉｊでは、格子点データ(ｘ_１ｊ,ｘ_２ｊ,ｘ_３ｊ,ｘ_４ｊ)と実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)との信号成分毎の差分が(９)式の行列Ｍ_ｉの各成分を用いて感度が加味された差分成分に変換され、さらに、これらの差分成分は規格化される。そして、第２の関数Ｇ_ｉｊにより、行列Ｍ_ｉの各成分を用いて感度が加味された差分成分からなる規格化信号に関するユークリッド距離に応じて単調減少する第２の加重Ｗ_２ｉｊを設定する。なお、(１１)式のｙ_１０,ｙ_２０,ｙ_３０は規格化定数であって、ユークリッド距離に応じて第２の加重Ｗ_２ｉｊの幅を規定する。したがって、規格化定数ｙ_１０,ｙ_２０,ｙ_３０は、実データに施す加重を調整するパラメータの一例である。
そして、加重生成部２３は、(１２)式に示したように、第１の加重Ｗ_１ｉｊと第２の加重Ｗ_２ｉｊとを結合した加重Ｗ_ｉｊを生成する。

このように加重再生成工程(ステップ２０４)では、色データ(色座標)間のユークリッド距離に関する２つの単調減少関数である第１の関数Ｆ_ｉｊおよび第２の関数Ｇ_ｉｊによって加重Ｗ_ｉｊを設定する。それによって、色空間での距離の差(色差距離)が大きい実データについては加重Ｗ_ｉｊを小さくして影響を少なくし、色差が小さい実データについては加重Ｗ_ｉｊを大きくして重要なデータとして扱う。すなわち、加重Ｗ_ｉｊは実データ対への重み付けを色差距離に応じて設定している。

ここで、第１の関数Ｆ_ｉｊおよび第２の関数Ｇ_ｉｊとして用いられる色データ(色座標)間のユークリッド距離に関する具体的な単調減少関数の例として、次の(１３)式が用いられる。また図１５は、例えば(１３)式により設定される加重の一例を示した図である。加重生成部２３では、初期加重生成工程(ステップ２０１)や加重再生成工程(ステップ２０４)において、図１５に示したような変換対象色(格子点データ)からの色差距離ｄに応じた加重が設定される。なお、上記した規格化定数により規定される「加重の幅」は、例えば図１５に示した加重曲線の半値幅×２として定義できる。

[行列成分再算出工程]
続く行列成分再算出工程(ステップ２０５)において算出部２４は、(６)式にて算出される出力変換値(ｙ´_１ｉ,ｙ´_２ｉ,ｙ´_３ｉ)と、(６)式にて出力変換値(ｙ´_１ｉ,ｙ´_２ｉ,ｙ´_３ｉ)を算出する際に用いたＣＭＹＫ色空間の実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)と対をなすＬ^＊ａ^＊ｂ^＊色空間の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)との間のユークリッド距離に、加重再生成工程(ステップ２０４)にて生成された加重Ｗ_ｉｊを乗じた２乗和Ｅ_ｊ((８)式)を最小にするという条件の下で、線形回帰分析の一例である最小２乗法により、(６)式の行列Ｍ_ｉの各成分ｍ_１１,ｍ_１２,ｍ_１３,ｍ_１４,…を再度求める。
ただし、加重Ｗ_ｉｊが行列Ｍ_ｉの成分または変換色であるＬ^＊ａ^＊ｂ^＊色空間の色データに依存する場合は、行列Ｍ_ｉの各成分ｍ_１１,ｍ_１２,ｍ_１３,ｍ_１４,…を一意的に定めることができない。そこで、その場合には、(８)式のＥ_ｊを最小にするという条件の下で、加重Ｗ_ｉｊと行列Ｍ_ｉの各成分ｍ_１１,ｍ_１２,ｍ_１３,ｍ_１４,…と変換色の最適値とを逐次近似の手法を用いて決定する。

[変換値算出工程]
行列成分再算出工程(ステップ２０５)に続く変換値算出工程(ステップ２０６)では、算出部２４により、行列成分再算出工程(ステップ２０５)により求められた(７)式の行列Ｍ_ｉを適用した(９)式に基づき、変換対象色である格子点データ(ｘ_１ｊ,ｘ_２ｊ,ｘ_３ｊ,ｘ_４ｊ)から変換色データ(ｙ_１ｊ,ｙ_２ｊ,ｙ_３ｊ)を算出する。算出部２４は、さらに、この計算された変換色データ(ｙ_１ｊ,ｙ_２ｊ,ｙ_３ｊ)を加重生成部２３に送る。加重生成部２３は、算出部２４にて算出された変換色データ(ｙ_１ｊ,ｙ_２ｊ,ｙ_３ｊ)を参照して、再び加重Ｗ_ｉｊを算出する。そして、算出部２４は、(８)式のＥ_ｊを最小とする成分ｍ_１１,ｍ_１２,ｍ_１３,ｍ_１４,…からなる行列Ｍ_ｉを求める。算出部２４は、この行列Ｍ_ｉを用いて変換色データ(ｙ_１ｊ,ｙ_２ｊ,ｙ_３ｊ)を再度算出する。色変換特性モデル生成部２０では、ステップ２０４の加重再生成工程からステップ２０６の変換値算出工程までの処理を繰り返すことによって変換色データ(ｙ_１ｊ,ｙ_２ｊ,ｙ_３ｊ)を収束させ、精度の高い行列Ｍ_ｉを得る。
そして算出部２４は、求めた(７)式の行列Ｍ_ｉを正変換の色変換特性モデルとして色変換処理部３０に設定する。それにより、色変換処理部３０は、(９)式に基づき、格子点データ取得部１０にて入力された変換対象色である格子点データ(ｘ_１ｊ,ｘ_２ｊ,ｘ_３ｊ,ｘ_４ｊ)から変換色である変換色データ(ｙ_１ｊ,ｙ_２ｊ,ｙ_３ｊ)を算出する。

引き続いて、上記の図１４のフローチャートに示した処理により、変換対象色としての３次元のＬ^＊ａ^＊ｂ^＊色空間での色データ(Ｌ^＊,ａ^＊,ｂ^＊)から、変換色としての４次元のＣＭＹＫ色空間での色データ(Ｃ,Ｍ,Ｙ,Ｋ)への逆変換の色変換特性モデルを生成する場合について述べる。

Ｌ^＊ａ^＊ｂ^＊色空間の格子点データからＣＭＹＫ色空間の変換色データへの逆変換の色変換特性モデルを生成する場合は、上記の(９)式を逆に解く(逆演算する)こととなる。ところが、(９)式においては、変換対象色が３次元の格子点データ(ｙ_１ｊ,ｙ_２ｊ,ｙ_３ｊ)であって、変換色が４次元の変換色データ(ｘ_１ｊ,ｘ_２ｊ,ｘ_３ｊ,ｘ_４ｊ)であることから、既知数不足となり、(９)式の逆演算を一意に解くことはできない。そこで、この場合は、変換色データの一部である例えば(ｘ_４ｊ)を変換対象色として設定し、(９)式を逆演算して、残りの変換色データである例えば(ｘ_１ｊ,ｘ_２ｊ,ｘ_３ｊ)を求める。具体的には、始めにＣＭＹＫ色空間の“Ｋ”を算出しておく手法が一般的である。例えば、Ｌ^＊ａ^＊ｂ^＊色空間の色データ(Ｌ^＊,ａ^＊,ｂ^＊)に対応するＣＭＹＫ色空間のＫ座標値を予め設計しておく方法や、色データ(Ｌ^＊,ａ^＊,ｂ^＊)とＫ座標値とが入力信号(変換対象色)として直接与えられる方法が用いられる。

ここで、Ｌ^＊ａ^＊ｂ^＊色空間の色データ(Ｌ^＊,ａ^＊,ｂ^＊)からＣＭＹＫ色空間の“Ｋ”を算出する処理を含んでいない場合には、変換対象色が３次元のＬ^＊ａ^＊ｂ^＊色空間の色データ(Ｌ^＊,ａ^＊,ｂ^＊)であって、変換色が３次元のＣＭＹ色空間の色データ(Ｃ,Ｍ,Ｙ)であるので、(９)式の逆演算によって第２色信号が生成される。この場合には、図１における格子点データ取得部１０にて格子点データ(Ｌ^＊,ａ^＊,ｂ^＊,Ｋ)が入力され、色変換処理部３０から第２色空間(ＣＭＹＫ色空間)の変換色データ(Ｃ,Ｍ,Ｙ,Ｋ)が出力される。

そこで、この場合には、色変換特性モデル生成部２０の加重生成部２３は、Ｌ^＊ａ^＊ｂ^＊色空間において加重ｗ_ｉｊを生成する。例えば、図１４の加重再生成工程(ステップ２０４)では、次の(１４)式に示したように、変換対象色である格子点データ(ｙ_１ｊ,ｙ_２ｊ,ｙ_３ｊ)とＬ^＊ａ^＊ｂ^＊色空間の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)との信号成分毎の差分を規格化し、求められた規格化信号のユークリッド距離に関する単調減少関数である第２の関数Ｇ_ｉｊにより、第２の加重ｗ_２ｉｊが設定される。また、次の(１５)式に示したように、変換色データ(ｘ_１ｊ,ｘ_２ｊ,ｘ_３ｊ,ｘ_４ｊ)の中の１つの信号成分である例えば(ｘ_４ｊ)(＝“Ｋ”)と、その実データ(ｘ_４ｉ)との差分を規格化し、変換色データ(“Ｋ”)の差分成分から求められる規格化信号のユークリッド距離に関する単調減少関数である第１の関数Ｆ_ｉｊにより、第１の加重ｗ_１ｉｊが設定される。そして、第１の加重ｗ_１ｉｊと第２の加重ｗ_２ｉｊとを(１６)式で示した関数Ｈで合成し、加重ｗ_ｉｊを生成する。なお、初期加重生成工程(ステップ２０１)でも、同様にして初期加重ｗ´_ｉｊが生成される。

そして、ステップ２０５の行列成分再算出工程では、(１６)式の加重ｗ_ｉｊが施された実データを用いた線形回帰分析が実施される。それにより、(９)式を逆演算する逆行列Ｍ_ｉ ^−１(逆変換の色変換特性モデル)を算出する。同様に、図１４のステップ２０２の行列成分算出工程でも、(１６)式の加重ｗ_ｉｊが施された実データを用いた線形回帰分析が実施される。
なお、(１４)式のｙ_１０,ｙ_２０,ｙ_３０は規格化定数であって、ユークリッド距離に応じて第２の加重ｗ_２ｉｊの幅を規定する。また、(１５)式のｘ_４０も規格化定数であって、ユークリッド距離に応じて第１の加重ｗ_１ｉｊの幅を規定する。したがって、規格化定数ｙ_１０,ｙ_２０,ｙ_３０,ｘ_４０は、実データに施す加重を調整するパラメータの一例である。

続いて、色変換特性モデル生成部２０にて逆変換の色変換特性モデルを生成する際に、初期加重生成工程(ステップ２０１)や加重再生成工程(ステップ２０４)にて実データに加重を施す場合に設定する加重の幅について述べる。
上記したように、色変換特性モデル生成部２０は、実データ評価部５０から評価用格子点データと対応付けられた評価値Ｓを評価値特性(図１０〜図１２参照)として取得する。この実データ評価部５０にて生成された評価値Ｓに関する評価値特性は、色変換特性モデルの色変換精度への影響度の観点から実データ対を色空間内の領域毎に評価するファクターである。そこで、色変換特性モデル生成部２０では、逆変換の色変換特性モデルを生成する際に、格子点データ(ｙ_１ｊ,ｙ_２ｊ,ｙ_３ｊ)と実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)とのユークリッド距離に応じて施す加重ｗ_ｉｊの幅を規定するパラメータの設定値(規格化定数の値)を、評価値Ｓに関する評価値特性に応じて調整する。それにより、逆変換の色変換特性モデルを生成する際に実データに施す加重ｗ_ｉｊを調整し、逆変換の色変換特性モデルの色変換精度の向上を図っている。

すなわち、評価値Ｓが小さい領域は、図６に示したようにＬ^＊ａ^＊ｂ^＊色空間での測色値である実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)がＬ^＊ａ^＊ｂ^＊色空間で偏在して分布する領域である。そのため、上記の図７を用いて説明したように、評価値Ｓが小さい領域においては、実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)に加重ｗ_ｉｊを施して生成される色変換特性モデルの傾き(変化率)は小さく、Ｌ^＊ａ^＊ｂ^＊色空間の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)とＣＭＹＫ色空間の実データ(ｘ_１ｉ,ｘ_２ｉ,ｘ_３ｉ,ｘ_４ｉ)との相関性は低くなる。そのため、評価値Ｓが小さい領域にて生成された色変換特性モデルにおいては、変換誤差が生じ易くなる。
そこで、色変換特性モデル生成部２０は、評価値Ｓが小さい領域において色変換特性モデルを生成する際には、実データに施す加重ｗ_ｉｊの幅を広げるように設定している。それにより、評価値Ｓが小さく実データ相互の相関性の低い領域では、色変換特性モデルの生成に使用するＬ^＊ａ^＊ｂ^＊色空間の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)が位置する領域は広げられる。それによって、かかる領域内での実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)の偏在が緩和され、評価値Ｓが小さい領域にておいても、生成される色変換特性モデルの色変換精度が向上される。

図１６は、評価値Ｓに応じて設定するＬ^＊ａ^＊ｂ^＊色空間の実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)に施す加重の幅を説明する図である。図１６に示したように、評価値Ｓが大きい領域Ａでは、Ｌ^＊ａ^＊ｂ^＊色空間での実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)の分散性は高いので、加重の幅を相対的に小さく設定しても、生成される色変換特性モデルの色変換精度は確保される。
ところが、評価値Ｓが小さい領域Ｂにおいては、Ｌ^＊ａ^＊ｂ^＊色空間での実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)の分散性は低いので、評価値Ｓが大きい領域Ａと同程度の加重の幅を設定すると、生成される色変換特性モデルの色変換精度は低下する。そこで、評価値Ｓが小さい領域Ｂでは、加重の幅を評価値Ｓが大きい領域Ａよりも大きく設定する。それにより、色変換特性モデルの生成に使用する実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ)に関するＬ^＊ａ^＊ｂ^＊色空間の領域は広がるため、この領域全体としての分散性は高まり、生成される色変換特性モデルの色変換精度は向上される。

次の図１７は、色変換特性モデル生成部２０にて用いられる評価値Ｓと加重の幅を規定するパラメータ値との対応関係を設定する関係曲線の一例を示した図である。色変換特性モデル生成部２０の加重生成部２３は、実データ評価部５０から取得した評価値Ｓに関する評価値特性(図１０および図１１)と、例えば図１７に示した関係曲線とを用いて、加重の幅を規定するパラメータの設定値(パラメータ値)を評価値Ｓに応じて設定する。
この場合の関係曲線としては、評価値Ｓが小さいほど加重の幅が大きく設定されるように、評価値Ｓに応じてパラメータ値(加重の幅)が単調減少する関係を規定する関数が用いられる。またその際に、図１７の関係曲線により評価値Ｓに応じて対応付けられるパラメータ値(加重の幅)は、パラメータ値(加重の幅)により生成される色変換特性モデルの色変換精度が所定値を満たすものであるように統計的に算出されて設定される。それにより、評価値や色変換特性モデルの算出の際に生じる不確定な要素に対応させて、安定的な色変換特性モデルの生成を図っている。

続いて、色変換特性モデル生成部２０において実行されるパラメータ値(規格化定数の値)の設定処理について説明する。
図１８は、色変換特性モデル生成部２０がパラメータ値を設定する際の処理の手順の一例を示したフローチャートである。図１８に示したように、色変換特性モデル生成部２０は、実データ評価部５０から、Ｋ座標値(カバレッジ値Ｋ(％))と対応付けられた評価値Ｓに関する評価値特性(図１０および図１１参照)を取得する(ステップ３０１)。
また、色処理装置１の記憶部(不図示)には、評価値Ｓと加重の幅を規定するパラメータ値との対応関係を設定する例えば図１７に示した関係曲線が記憶されている。そして、色変換特性モデル生成部２０は、記憶部から図１７の関係曲線を取得する(ステップ３０２)。
そして、取得した評価値Ｓに関する評価値特性と関係曲線とを加重生成部２３に設定する(ステップ３０３)。それにより、色変換特性モデル生成部２０の加重生成部２３は、カバレッジ値Ｋ(％)と対応付けられた評価値Ｓに関する評価値特性と、評価値Ｓとパラメータ値との対応関係を設定する関係曲線とを用いて、Ｋ座標値(カバレッジ値Ｋ(％))に応じたパラメータ値(規格化定数値)を算出して設定する(ステップ３０４)。

図１９は、ステップ３０４にて加重生成部２３が行うＫ座標値(カバレッジ値Ｋ(％))に応じたパラメータ値(規格化定数値)を算出する処理を説明する図である。図１９は、上記の図１０に示した画像形成装置(プリンタ)Ａでの評価値Ｓに関する評価値特性と図１７に示した関係曲線とから、画像形成装置(プリンタ)Ａに設定されるカバレッジ値Ｋ(％)に対応するパラメータ値を算出する場合を示している。
図１９に示したように、加重生成部２３は、図１９の第３象限位置に表した評価値特性(図１０)から各カバレッジ値Ｋ(％)に対応する評価値Ｓを読み出す。そして、図１９の第２象限位置に表した関係曲線(図１７)を用いて、読み出した評価値Ｓに対応するパラメータ値を算出する。それにより、各カバレッジ値Ｋ(％)に対応するパラメータ値を求める。さらには、求めた各カバレッジ値Ｋ(％)に対応するパラメータ値から、図１９の第１象限位置に表したカバレッジ値Ｋ(％)の全域に亘るパラメータ値Ｐを設定するパラメータ値設定関数を算出する。カバレッジ値Ｋ(％)の全域に亘るパラメータ値設定関数の算出に際しては、次の(１７)式を用いる。

ここで、(１７)式には、パラメータ値設定関数の形状を調整する３つの係数α、β、Ｋ_ｉｎｆが設定されている。αは、パラメータ値設定関数の最大値を調整するための係数である。βは、パラメータ値設定関数の勾配を調整するための係数である。また、Ｋ_ｉｎｆは、パラメータ値設定関数の変曲点を調整するための係数である。加重生成部２３では、求めた各カバレッジ値Ｋ(％)に対応するパラメータ値を用いて３つの係数α、β、Ｋ_ｉｎｆを算出し、パラメータ値設定関数を設定する。
そして、画像形成装置(プリンタ)Ａでは、図１９に示した画像形成装置(プリンタ)Ａでのパラメータ値設定関数により、カバレッジ値Ｋ(％)に対応したパラメータ値、すなわち実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ,ｘ_４０)に施される加重ｗ_ｉｊの幅を規定する規格化定数ｙ_１０,ｙ_２０,ｙ_３０,ｘ_４０を設定する。

同様に、図２０は、ステップ２０４にて加重生成部２３が行うＫ座標値(カバレッジ値Ｋ(％))に応じたパラメータ値(規格化定数)を算出する処理を説明する図である。図２０は、上記の図１１に示した画像形成装置(印刷機)Ｂでの評価値Ｓに関する評価値特性と図１７に示した関係曲線とから、画像形成装置(印刷機)Ｂに設定されるカバレッジ値Ｋ(％)に対応するパラメータ値を算出する場合を示している。
この場合にも、図１９の場合と同様にして、画像形成装置(印刷機)Ｂでのパラメータ値設定関数を設定する。そして、画像形成装置(印刷機)Ｂでは、図２０に示した画像形成装置(印刷機)Ｂでのパラメータ値設定関数により、カバレッジ値Ｋ(％)に対応したパラメータ値、すなわち実データ(ｙ_１ｉ,ｙ_２ｉ,ｙ_３ｉ,ｘ_４０)に施される加重ｗ_ｉｊの幅を規定する規格化定数ｙ_１０,ｙ_２０,ｙ_３０,ｘ_４０を設定する。

図２１は、画像形成装置(プリンタ)Ａでのパラメータ値設定関数と、画像形成装置(印刷機)Ｂでのパラメータ値設定関数とを比較した図である。図２１に示したように、画像形成装置(プリンタ)Ａおよび画像形成装置(印刷機)Ｂにて使用される実データ対の相違により、設定されるパラメータ値設定関数は異なるものとなる。
さらに、図２２は、各種の画像形成装置にて算出されるパラメータ値設定関数を示した図である。図２２に示したように、各画像形成装置では使用される実データ対によりそれぞれ異なるパラメータ値設定関数が設定され、設定されたパラメータ値設定関数により、カバレッジ値Ｋ(％)に対応した規格化定数ｙ_１０,ｙ_２０,ｙ_３０,ｘ_４０が設定される。それにより、評価値Ｓが小さい領域においても、生成される逆変換の色変換特性モデルの色変換精度が向上する。

このように、本実施の形態の色処理装置１では、色変換特性モデルの生成の際に使用される実データ対を、この実データ対により生成される色変換特性モデルの色変換精度への影響度の観点から評価し、色空間内での領域毎の実データ対の評価値Ｓを算出する。さらに、算出された色空間内での領域(カバレッジ値Ｋ(％))毎の評価値Ｓに応じたパラメータ値Ｐを設定するパラメータ値設定関数を設定する。そして、設定されたパラメータ値設定関数に基づき、逆変換の色変換特性モデルの生成の際に使用される実データ対に施される重み付けを加重の幅を規定するパラメータ値(規格化定数値)によって設定し、逆変換の色変換特性モデルを生成する際に使用する実データ対に施す重み付けを調整する。それにより、評価値Ｓが小さい領域においても逆変換の色変換特性モデルの色変換精度が向上され、正変換および逆変換双方の色変換特性モデルの色変換精度が色域全域で高められる。

図２３は、パラメータ値設定関数に基づき設定されたパラメータ値(規格化定数)を用いて逆変換の色変換特性モデルを生成した場合のモデル誤差(図４参照)を、従来の方法に基づき設定された所定のパラメータ値を用いた場合と比較した図である。図２３(ａ)は、各種の画像形成装置Ｃ〜Ｌでのモデル誤差を二乗平均平方根(ＲＭＳ)を用いて比較したものであり、(ｂ)は、各種の画像形成装置Ｃ〜Ｌでのモデル誤差をモデル誤差の最大値(最大誤差)を用いて比較したものである。また、図２３は、ＣＭＹＫ色空間の９×９×９×９＝６５６１個の格子点データにおいてモデル誤差を算出したものである。
図２３に示したように、パラメータ値設定関数に基づきパラメータ値(規格化定数)を設定することで、モデル誤差は改善される。

次に、図２４は、本実施の形態の色処理装置１にて生成されたＤＬＵＴ(Direct Look-Up Table)を搭載する画像形成装置２００の構成を示したブロック図である。
図２４に示した画像形成装置２００は、例えばデジタルカラープリンタであって、外部機器から入力された画像データに対して所定の画像処理を施す画像処理装置２１０、画像形成装置２００全体の動作を制御する制御部２３０、各色成分からなる画像データに基づき記録媒体上に画像形成を行う例えば電子写真方式のプリンタエンジンにて実現される画像形成部２４０を備えている。

画像処理装置２１０は、例えばパーソナルコンピュータ(ＰＣ)３００等といった外部機器からの画像データの入力を受け付ける入力インターフェース２１１、入力インターフェース２１１にて受け付けた画像データを一時記憶する入力バッファ２１２、ＰＤＬ（Page Description Language：ページ記述言語）形式の画像データを解析して中間データを生成するＰＤＬ解析部２１３、ＰＤＬ解析部２１３にて生成された中間データを画素の並びで表現された印刷用の画像データ(ラスタ画像データ)に展開(レンダリング)するレンダリング処理部２１４、レンダリング処理部２１４でのレンダリング処理に際して作業領域として使用される中間バッファ２１５、レンダリングされた画像データを色変換処理する色処理手段である色変換処理部２１６、色変換された画像データに対してスクリーン処理を行うスクリーン処理部２１７を備えている。
そして、画像処理装置２１０の色変換処理部２１６には、本実施の形態の色処理装置１にて生成されたＤＬＵＴが備えられる。そして、レンダリング処理部２１４にてレンダリングされた画像データをＤＬＵＴを用いて色変換処理する。色変換処理部２１６にて色変換処理された画像データは、スクリーン処理部２１７にて所定のスクリーン処理が行われた後、画像形成部２４０にて用紙上に画像形成される。

また、図２５は、本実施の形態の色処理装置１のハードウェア構成を示した図である。図２５に示したように、色処理装置１は、色変換処理を実行するに際して、予め定められた処理プログラムに従ってデジタル演算処理を実行する演算手段の一例としてのＣＰＵ１０１、ＣＰＵ１０１により実行される処理プログラム等が格納されるＲＡＭ１０２、ＣＰＵ１０１により実行される処理プログラム等にて用いられる設定値等のデータが格納されるＲＯＭ１０３、書き換え可能で電源供給が途絶えた場合にもデータを保持できるＥＥＰＲＯＭやフラッシュメモリ等の不揮発性メモリ１０４、色処理装置１に接続される各部との信号の入出力を制御するインターフェース部１０５を備えている。
また、色処理装置１には外部記憶装置６０が接続されている。そして、外部記憶装置６０には、色処理装置１により実行される処理プログラムが格納されており、色処理装置１がこの処理プログラムを読み込むことによって、本実施の形態の色処理装置１での色変換処理が実行される。

すなわち、上記した格子点データ取得部１０、色変換特性モデル生成部２０、色変換処理部３０、色変換データ出力部４０、実データ評価部５０の各機能を実現するプログラムを外部記憶装置６０から色処理装置１内のＲＡＭ１０２に読み込む。そして、ＲＡＭ１０２に読み込まれた処理プログラムに基づいて、ＣＰＵ１０１が各種処理を行う。この処理プログラムは、例えば外部記憶装置６０としてのハードディスクやＤＶＤ−ＲＯＭ等に格納された処理プログラムが、ＲＡＭ１０２にロードされて提供される。また、その他の提供形態として、予めＲＯＭ１０３に格納された状態にてＲＯＭ１０３からＲＡＭ１０２にロードされる形態がある。さらに、ＥＥＰＲＯＭ等の書き換え可能なＲＯＭ１０３を備えている場合には、色処理装置１がアッセンブリされた後に処理プログラムがＲＯＭ１０３に提供され、ＲＯＭ１０３からＲＡＭ１０２にロードされる形態がある。また、インターネット等のネットワークを介して色処理装置１に処理プログラムが伝送され、色処理装置１のＲＡＭ１０２にインストールされる形態がある。

以上説明したように、本実施の形態の色処理装置１においては、色変換特性モデルの生成の際に使用される実データ対を、この実データ対により生成される色変換特性モデルの色変換精度への影響度の観点から評価し、色空間内での領域毎の実データ対の評価値Ｓを算出する。さらに、算出された色空間内での領域毎の評価値Ｓに応じたパラメータ値Ｐを設定するパラメータ値設定関数を設定する。そして、設定されたパラメータ値設定関数に基づき、逆変換の色変換特性モデルの生成の際に使用される実データ対に施される重み付けを加重の幅を規定するパラメータ値(規格化定数値)によって設定し、逆変換の色変換特性モデルを生成する際に使用する実データ対に施す重み付けを調整する。それにより、評価値Ｓが小さい領域においても逆変換の色変換特性モデルの色変換精度が向上され、正変換および逆変換双方の色変換特性モデルの色変換精度が色域全域で高められる。

本実施の形態が適用される色処理装置の全体構成を示すブロック図である。実データ評価部の機能構成を示すブロック図である。色変換特性モデルの機能を説明する図である。正変換により変換された変換色が、逆変換によって元のＣＭＹＫ色空間の色座標に変換される際に生じる誤差を説明する図である。ＣＭＹＫ色空間の実データおよびこれらに対応するＬ^＊ａ^＊ｂ^＊色空間の実データに関するぞれぞれの色空間内での分散状態を表す図である。ＣＭＹＫ色空間の実データおよびこれらに対応するＬ^＊ａ^＊ｂ^＊色空間の実データに関するぞれぞれの色空間内での分散状態を表す図である。Ｌ^＊ａ^＊ｂ^＊色空間での実データの分散状態と、実データ対を用いて生成される色変換特性モデルの色変換精度との関係を説明する図である。実データ対評価値算出部にて相対分散値を算出する際の処理の手順の一例を示したフローチャートである。ＣＭＹＫ色空間にて評価用格子点の周囲の領域内に存在する実データと、それらに対応するＬ^＊ａ^＊ｂ^＊色空間の実データとを説明する図である。評価値出力部から色変換特性モデル生成部に出力される評価用格子点データと対応付けられた評価値の一例を示した図である。評価値出力部から色変換特性モデル生成部に出力される評価用格子点データと対応付けられた評価値の一例を示した図である。各種の画像形成装置の評価値特性を示した図である。色変換特性モデル生成部の機能構成を示すブロック図である。色変換特性モデル生成部にて色変換特性モデルを生成する際の処理の手順の一例を示したフローチャートである。加重の一例を示した図である。評価値に応じて設定するＬ^＊ａ^＊ｂ^＊色空間の実データに施す加重の幅を説明する図である。色変換特性モデル生成部にて用いられる評価値と加重の幅を規定するパラメータ値との対応関係を設定する関係曲線の一例を示した図である。色変換特性モデル生成部がパラメータ値を設定する際の処理の手順の一例を示したフローチャートである。加重生成部が行うＫ座標値に応じたパラメータ値を算出する処理を説明する図である。加重生成部が行うＫ座標値に応じたパラメータ値を算出する処理を説明する図である。パラメータ値設定関数を比較した図である。各種の画像形成装置にて算出されるパラメータ値設定関数を示した図である。パラメータ値設定関数に基づき設定されたパラメータ値を用いて逆変換の色変換特性モデルを生成した場合のモデル誤差を、所定のパラメータ値による従来の生成方法を用いた場合と比較した図である。色処理装置にて生成されたＤＬＵＴを搭載する画像形成装置の構成を示したブロック図である。色処理装置のハードウェア構成を示した図である。

符号の説明

１…色処理装置、１０…格子点データ取得部、２０…色変換特性モデル生成部、２１…実データ対入力部、２２…実データ対記憶部、２３…加重生成部、２４…算出部、３０…色変換処理部、４０…色変換データ出力部、５０…実データ評価部、５１…実データ対取得部、５２…実データ対記憶部、５３…評価用格子点生成部、５４…実データ対評価値算出部、５５…評価値出力部、２００…画像形成装置、２１６…色変換処理部

Claims

第１色空間の変換対象色データを取得する変換対象色データ取得部と、
前記第１色空間の色と変換先である第２色空間の色との対応関係を規定する色変換特性モデルの算出に際し用いる前記第１色空間の色データと、当該色変換特性モデルの算出に際し用いる当該第２色空間の色データとを対とする色データ対を取得する色データ対取得部と、
前記色データ対取得部にて取得された前記色データ対の前記第１色空間の色データと前記変換対象色データ取得部にて取得した前記変換対象色データとの色差距離に応じて当該第１色空間の色データに重み付けを行い、当該第１色空間の色データが重み付けされた当該色データ対を用いて前記色変換特性モデルを算出する色変換特性モデル算出部と、
前記色データ対を構成する前記第１色空間の色データに関する当該第１色空間にて所定間隔毎に配列された基準色座標の周辺領域での分散と、当該周辺領域内の当該第１色空間の色データと対をなす前記第２色空間の色データに関する分散とを算出し、算出した当該分散の相対比により当該色データ対に関する前記色変換特性モデルの色変換精度への影響度を評価する評価値を算出する評価値算出部とを備え、
前記色変換特性モデル算出部は、前記評価値算出部にて算出された前記評価値に基づいて、前記色変換特性モデルの算出に際し前記第１色空間の色データに施す前記重み付けを調整することを特徴とする色処理装置。
前記色変換特性モデル算出部は、前記第１色空間の色データに施す前記重み付けの調整を、当該第１色空間の色データと前記変換対象色データとの色差距離と当該重み付けとの対応関係を前記評価値に基づいて変更することにより行うことを特徴とする請求項１記載の色処理装置。
前記評価値算出部は、前記第１色空間にて所定間隔毎に配列された基準色座標の周辺領域毎に前記評価値を算出し、算出された当該評価値を当該基準色座標の一の座標成分と対応付けて前記色変換特性モデル算出部に送り、前記色変換特性モデル算出部は、当該一の座標成分と対応付けられた前記評価値に基づき前記重み付けを調整することを特徴とする請求項１記載の色処理装置。
前記評価値算出部は、前記基準色座標の周辺領域毎に算出された複数の前記評価値に関する前記一の座標成分の座標値毎の最大値、最小値および平均値の何れかを、当該一の座標成分と対応付けられた当該評価値として規定することを特徴とする請求項３記載の色処理装置。
前記色変換特性モデル算出部は、前記評価値算出部にて算出された前記評価値に対して前記色変換特性モデルの色変換精度が所定値を満たすように統計的に対応付けた前記重み付けを設定する前記調整を行うことを特徴とする請求項１記載の色処理装置。
コンピュータに、
第１色空間の変換対象色データを取得する機能と、
前記第１色空間の色と変換先である第２色空間の色との対応関係を規定する色変換特性モデルの算出に際し用いる当該第１色空間の色データと、当該色変換特性モデルの算出に際し用いる当該第２色空間の色データとを対とする色データ対を取得する機能と、
前記色データ対の前記第１色空間の色データと前記変換対象色データとの色差距離に応じて当該第１色空間の色データに重み付けを行い、当該第１色空間の色データが重み付けされた当該色データ対を用いて前記色変換特性モデルを算出する機能と、
前記色データ対を構成する前記第１色空間の色データに関する当該第１色空間にて所定間隔毎に配列された基準色座標の周辺領域での分散と、当該周辺領域内の当該第１色空間の色データと対をなす前記第２色空間の色データに関する分散とを算出し、算出した当該分散の相対比により当該色データ対に関する前記色変換特性モデルの色変換精度への影響度を評価する評価値を算出する機能と、
前記評価値に基づいて、前記色変換特性モデルの算出に際し前記第１色空間の色データに施す前記重み付けを調整する機能と
を実現させることを特徴とするプログラム。
前記重み付けを調整する機能は、前記第１色空間の色データに施す前記重み付けの調整を、当該第１色空間の色データと前記変換対象色データとの色差距離と当該重み付けとの対応関係を前記評価値に基づいて変更することにより行う請求項６記載のプログラム。
前記評価値を算出する機能は、前記第１色空間にて所定間隔毎に配列された基準色座標の周辺領域毎に前記評価値を算出し、算出された当該評価値を当該基準色座標の一の座標成分と対応付けるとともに、前記重み付けを調整する機能は、当該一の座標成分と対応付けられた当該評価値に基づき前記重み付けを調整することを特徴とする請求項６記載のプログラム。
前記評価値を算出する機能は、前記基準色座標の周辺領域毎に算出された複数の前記評価値に関する前記一の座標成分の座標値毎の最大値、最小値および平均値の何れかを、当該一の座標成分と対応付けられた当該評価値として規定することを特徴とする請求項８記載のプログラム。
前記重み付けを調整する機能は、算出された前記評価値に対して前記色変換特性モデルの色変換精度が所定値を満たすように統計的に対応付けた前記重み付けを設定する前記調整を行うことを特徴とする請求項６記載のプログラム。