(57)【要約】【課題】  色空間に分布する標本点の変化の傾向に従って標本点を滑らかに接続することができ、階調性の良好な色再現が実現可能とする【解決手段】  入力色信号値と出力色信号値の対応関係を示す複数の標本点を入力する手段と、前記標本点の入力信号値を媒介変数に変換するためのパラメータを設定する手段と、前記入力色信号の表色系の次元ごとに基底関数を算出し、前記標本点および該基底関数に基づき係数を設定する手段と、前記パラメータを用いて入力色値を媒介変数に変換する手段と、前記変換された媒介変数に対して、前記係数と前記入力色信号の表色系の次元ごとに算出された基底関数の積とに基づく多重総和演算を行うことにより出力色値を算出する手段とを有することを特徴とする。【選択図】  図５
2008294931.tif 000002
【特許請求の範囲】
【請求項１】
入力色信号値と出力色信号値の対応関係を示す複数の標本点を入力する手段と、
前記標本点の入力信号値を媒介変数に変換するためのパラメータを設定する手段と、
前記入力色信号の表色系の次元ごとに基底関数を算出し、前記標本点および該基底関数に基づき係数を設定する手段と、
前記パラメータを用いて入力色値を媒介変数に変換する手段と、
前記変換された媒介変数に対して、前記係数と前記入力色信号の表色系の次元ごとに算出された基底関数の積とに基づく多重総和演算を行うことにより出力色値を算出する手段とを有することを特徴とする色処理装置。
【請求項２】
前記パラメータは、前記標本点の入力信号値の次元ごとに前記媒介変数が均等な間隔になるように設定されることを特徴とする請求項１記載の色処理装置。
【請求項３】
前記標本点は、入力色域の内部および境界近傍の点であることを特徴とする請求項１または２記載の色処理装置。
【請求項４】
前記係数はスプライン係数であり、前記多重総和演算はスプライン関数に応じた演算であることを特徴とする請求項１乃至３のいずれかに記載の色処理装置。
【請求項５】
前記入力色信号の表色系と前記出力色信号の表色系は異なることを特徴とする請求項１乃至４のいずれかに記載の色処理装置。
【請求項６】
請求項１乃至５のいずれかに記載の色処理装置をコンピュータを用いて実現するために、コンピュータが読み取り可能に記憶媒体に記憶されたプログラム。
【請求項７】
入力色信号値と出力色信号値の対応関係を示す複数の標本点を入力し、
前記標本点の入力信号値を媒介変数に変換するためのパラメータを設定し、
前記入力色信号の表色系の次元ごとに基底関数を算出し、前記標本点および該基底関数に基づき係数を設定し、
前記パラメータを用いて入力色値を媒介変数に変換し、
前記変換された媒介変数に対して、前記係数と前記入力色信号の表色系の次元ごとに算出された基底関数の積とに基づく多重総和演算を行うことにより出力色値を算出することを特徴とする色処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、色を算出するための色処理方法および装置に関する。
【背景技術】
【０００２】
近年、カラー画像を扱うデバイスの普及により、同じカラー画像を異なるデバイスで扱う機会が増加している。例えば、ディジタルカメラで撮影した画像の出力先としてモニタやプリンタが考えられる。しかし、デバイス毎に色再現が可能な範囲（以下、色域という）が異なるため、同じ画像であってもモニタで表示される色とプリンタで出力される色が異なるという問題がある。そこで、異なるデバイス間での色域の差を吸収し、可能な限り同等の色再現を実現するためのＣＭＳ（Ｃｏｌｏｒ  Ｍａｎａｇｅｍｅｎｔ  Ｓｙｓｔｅｍ）が重要になっている。
【０００３】
一般的なＣＭＳは、入出力デバイスの特性を記述したデバイスプロファイルを用いて実現される。このデバイスプロファイルには、デバイス依存の色空間（ＲＧＢ、ＣＭＹＫなど）とデバイス非依存の色空間（ＣＩＥＬＡＢ、ＣＩＥＬＵＶなど）を相互変換する変換式や変換テーブル（ＬＵＴ：ＬｏｏｋＵｐ  Ｔａｂｌｅ）などが記述されている。このデバイス非依存の色空間はＰＣＳ（Ｐｒｏｆｉｌｅ  Ｃｏｎｎｅｃｔｉｏｎ  Ｓｐａｃｅ）と呼ばれ、ある光源の下で一意に定まるため、ＰＣＳ上でカラーマッチングを行うことでデバイス間の色域の差を吸収することが可能となる。より具体的には、カラーモニタのデバイスプロファイルを用いてカラーモニタのデバイス依存色であるＲＧＢ値をデバイス非依存色であるＬ＊ａ＊ｂ＊値に変換する。そして、出力するプリンタのデバイスプロファイルを用いて、変換したＬ＊ａ＊ｂ＊値をプリンタのデバイス依存色であるＣＭＹＫ値に変換する。
【０００４】
ここで、ディジタルカメラ、スキャナ、モニタ等のＲＧＢ値をデバイス色空間とする機器とプリンタ等のＣＭＹＫ値をデバイス色空間とする機器間のカラーマッチングでは、色空間の間で非線形性が強いためＬＵＴを用いられる。但し、例えば一般的なカラーモニタで表現可能な色の数はＲＧＢ各２５６階調、全１６７０万色以上にも及ぶためその全ての色のＲＧＢ値に対応するＬ＊ａ＊ｂ＊値をＬＵＴに記述することは非現実的である。そこでＬＵＴにはいくつかの代表色の対応関係を記述しておき、それ以外の色はＬＵＴに記述された代表色の関係から補間演算により算出する。この補間演算には立方体補間や四面体補間等の線形補間が用いられる事が多い。立方体補間は、変換する色を囲む立方格子上の格子点の８点で補間する方法である。また、四面体補間は図１６に示すように、前述の立方体を六分割して得られる四面体の４頂点（点Ｐ１，Ｐ３，Ｐ４，Ｐ８）を用いて補間する方法であり、図の点Ｐに対応する四面体補間式は次式になる。
Ｐ＝Ｐ１＋ＤＬ＊（Ｐ８−Ｐ１）＋Ｄａ＊（Ｐ３−Ｐ１）＋Ｄｂ＊（Ｐ４−Ｐ１）
ここで、ＤＬ、Ｄａ、Ｄｂは点Ｐを囲む立方体において、点Ｐ１を基準として点Ｐの位置を示す正規化した各軸方向の距離である。この四面体補間は補間演算に用いる点が４点と少なく、演算も容易なため広く利用されている。
【０００５】
しかしながら、ＬＵＴの格子間隔が粗い場合に線形補間を用いると、立方体や四面体で囲まれた局所的な領域から値を算出するため、色算出精度が低下する。また、立方体や四面体の表面を境界として線形補間が切り替わることに起因して、値が急峻に変化して滑らかさが失われ、再現画像に擬似輪郭が現れることがある。
【０００６】
このような問題を解決する手法として、特許文献１において非線型演算による色算出法が提案されている。この手法では色を算出する際に、複数の標本点の三次元的な変化の傾向を自然に反映した滑らかな非線型補間式で色を算出するため、少ない標本点数でも精度良く色を算出することが可能である。
【特許文献１】特開２００５−１７６２０３号公報
【発明の開示】
【発明が解決しようとする課題】
【０００７】
しかしながら、特許文献１では、非線型補間式を算出する際に標本点を均等に配置する必要があった。従って、入力色域の最大明度および最大彩度を再現可能にするためには、標本点を立方体あるいは円筒形に設定して入力色域全体を包含する必要があった。その結果、図１８に示すように、色域外にも多くの標本点が設定されることになり、設定した標本点数から期待される補間精度を十分に得ることができなかった。また、入力色域全体を包含するように標本点を設定すると色域内部と色域境界部を比較した場合、相対的に色域境界部の標本点数が少なくなるため、高彩度の色の再現性が低下するという問題があった。
【０００８】
本発明は、色空間に分布する標本点の変化の傾向に従って標本点を滑らかに接続することができ、階調性の良好な色再現が実現可能とすることを目的とする。
【課題を解決するための手段】
【０００９】
以上のような目的を達成するために、本発明は、入力色信号値と出力色信号値の対応関係を示す複数の標本点を入力する手段と、前記標本点の入力信号値を媒介変数に変換するためのパラメータを設定する手段と、前記入力色信号の表色系の次元ごとに基底関数を算出し、前記標本点および該基底関数に基づき係数を設定する手段と、前記パラメータを用いて入力色値を媒介変数に変換する手段と、前記変換された媒介変数に対して、前記係数と前記入力色信号の表色系の次元ごとに算出された基底関数の積とに基づく多重総和演算を行うことにより出力色値を算出する手段とを有することを特徴とする。
【発明の効果】
【００１０】
本発明によれば、色空間に分布する標本点の変化の傾向に従って標本点を滑らかに接続することができ、階調性の良好な色再現が実現可能である。
【発明を実施するための最良の形態】
【００１１】
（第１実施例）
図１は、本実施例にかかる色処理装置のブロック図である。１０１はＣＰＵ、１０２はメインメモリ、１０３はＳＣＳＩインタフェース、１０４はネットワークインタフェース、１０５はＨＤＤ、１０６はグラフィックアクセラレータである。１０７はカラーモニタ、１０８はＵＳＢコントローラ、１０９はカラープリンタ、１１０はキーボード／マウスコントローラ、１１１はキーボード、１１２はマウス、１１３はローカルエリアネットワーク、１１４はＰＣＩバスである。例えば、コンピュータとその周辺機器で構成することができる。
【００１２】
本実施例では、ＰＣＳとしてＣＩＥＬＡＢ色空間用いる。そして、入力プロファイルおよび出力プロファイルを用いて、カラーモニタの色空間であるＲＧＢ値で表される画像をプリンタの色空間であるＣＭＹＫ値で表される画像に変換する処理について説明する。
【００１３】
本実施例では、色変換アプリケーションを用いて画像の色を変換する。ユーザが色処理装置に対して色変換アプリケーションの起動を指示すると、ユーザの指示を受けたＯＳプログラムは記憶媒体であるＨＤＤ１０５に記憶されているプログラムである色変換アプリケーションを起動し、メインメモリ１０２に展開する。色変換アプリケーションが起動されると図２に示すウィンドウがカラーモニタ１０７上に表示される。
【００１４】
図２に示した色変換アプリケーションのウィンドウについて説明する。２０１は変換する画像ファイルの入力を指示するためのボタンである。２０２は入力デバイスプロファイル（以下、Ｓｒｃプロファイル）の入力を指示するためのボタン、２０３は出力デバイスプロファイル（以下、Ｄｓｔプロファイル）の入力を指示するためのボタンである。２０４は色変換後の画像の出力を指示するためのボタン、２０５は画像の色変換の様子を表示するための表示領域、２０６は画像の色変換を指示するためのボタン、２０７はアプリケーションの終了ボタンである。
【００１５】
色変換アプリケーションの動作について図４を用いて説明する。ステップ４０１では、色変換する画像ファイルを読み込む。色変換アプリケーションは、ユーザによって画像ファイル入力ボタン２０１により画像の読み込みが指示されると、指定された画像がＨＤＤ１０５もしくはネットワークＩ／Ｆ１０４経由で読み込まれる。読み込まれた画像は２０５（ａ）の表示領域に表示される。
【００１６】
ステップ４０２では、カラーモニタデバイス依存の色信号であるＲＧＢ値をＰＣＳの色信号であるＬ＊ａ＊ｂ＊値に変換するＬＵＴを持つＳｒｃプロファイルを読み込む。色変換アプリケーションにおいて、ユーザによってＳｒｃプロファイル入力ボタン２０２によりＳｒｃプロファイルの読み込みが指示されると指定されたプロファイルがＨＤＤ１０５もしくはネットワークＩ／Ｆ１０４経由で読み込まれる。ここで、ＳｒｃプロファイルのＬＵＴが、Ｒ、Ｇ、Ｂ各５スライスであり計１２５個のＲＧＢ値と、それに対応するＬ＊ａ＊ｂ＊値との組を持つとすると、図３（ａ）のようなデータ構造で構成される。
【００１７】
ステップ４０３では、ＰＣＳの色信号であるＬ＊ａ＊ｂ＊値からプリンタデバイス依存の色信号であるＣＭＹＫ値に変換するＬＵＴを持つＤｓｔプロファイルを読み込む。色変換アプリケーションにおいて、ユーザによってＤｓｔプロファイル入力ボタン２０３よりＤｓｔプロファイルの読み込みが指示されると、指定されたプロファイルがＨＤＤ１０５もしくはネットワークＩ／Ｆ１０４経由で読み込まれる。ここで、Ｄｓｔプロファイルも、Ｓｒｃプロファイルと同様に１２５点のＬ＊ａ＊ｂ＊値とＣＭＹＫ値との組を持つとすると、図３（ｂ）のようなデータ構造で構成される。
【００１８】
以上の処理に続き、変換実行ボタン２０６の指示に基づき、読み込んだ画像の色変換を実行する。ステップ４０４では、ステップ４０２で入力したＳｒｃプロファイルとステップ４０３で入力したＤｓｔプロファイルに基づき、ステップ４０１で入力した画像の色を変換する。変換方法については後述する。
【００１９】
ステップ４０５では、変換後の画像を２０５（ｂ）の表示領域に表示する。その後、画像ファイル出力ボタン２０４の指示に基づき変換後の画像を出力する。そして、終了ボタン２０７の指示に基づき全てのデータ、並びにアプリケーションをメインメモリ１０２から消去し、アプリケーションを終了する。
【００２０】
以下では、ステップ４０４における画像の色変換について、図５のフローチャートを用いて説明する。ステップ５０１では、スプライン関数の生成に利用する標本点の入力色信号値と出力色信号値を設定する。設定方法の詳細は後述する。
【００２１】
ステップ５０２では、標本点の入力色信号値に対して媒介変数を設定する。設定方法の詳細は後述する。
【００２２】
ステップ５０３では、ステップ５０２で設定した媒介変数と、ステップ５０１で設定した標本点に基づき基底関数およびスプライン係数を設定する。スプライン係数の算出方法の詳細は後述する。
【００２３】
ステップ５０４では、ステップ５０３で設定したスプライン係数に基づいて生成したスプライン関数を利用してステップ４０１で入力した画像の色を変換する。画像の色の変換方法の詳細については後述する。
【００２４】
以下では、ステップ５０１における標本点の設定について、図６および図７のフローチャートを用いて説明する。
【００２５】
まず、図６について説明する。図６はＣＩＥＬＡＢ色空間をａ＊ｂ＊平面、Ｌ＊ａ＊平面でそれぞれ示した模式図であり、入力色域境界、出力色域境界、標本点の入力色信号、標本点の出力色信号をそれぞれ図示している。
【００２６】
次に、図７のフローチャートについて説明する。ステップ７０１では、ステップ４０２で入力したＳｒｃプロファイルを基にＣｏｎｖｅｘ  Ｈｕｌｌ等の所望の方法で入力色域境界を算出する。
【００２７】
ステップ７０２では、ステップ７０１で算出した入力色域境界に基づいて、標本点の入力色信号を入力色域の内部および入力色域境界近傍に設定する。本実施例では、入力色域に対して明度方向Ｎ１点、彩度方向Ｎ２点、色相方向Ｎ３点の計Ｎ１＊Ｎ２＊Ｎ３点の標本点を設定する。明度方向、彩度方向に関しては標本点を等間隔に設定する。色相方向に関しては入力色域のプライマリカラー（レッド、グリーン、ブルー、シアン、マゼンタ、イエロー）の色相に標本点を設定する。但し、標本点の入力色信号の設定方法や標本点数はこの限りではなく、標本点はａ＊軸、ｂ＊軸に設定してもよいし、標本点数を増減しても構わない。また、標本点の間隔も任意に設定可能である。尚、彩度Ｃおよび色相ＨはＬ＊ａ＊ｂ＊値から次式で求める。
【００２８】
【数１】
2008294931.tif 000003
【００２９】
Ｈ＝ｔａｎ−１（ｂ／ａ）
ステップ７０３では、ステップ４０３で入力したＤｓｔプロファイルを基にＣｏｎｖｅｘ  Ｈｕｌｌ等の所望の方法で出力色域境界を算出する。
【００３０】
ステップ７０４では、ＳｒｃプロファイルとＤｓｔプロファイルの関係およびステップ７０１で算出した入力色域境界とステップ７０３で算出した出力色域境界の関係を用いて所望の方法で標本点の出力色信号値を算出する。例えば、ＳｒｃプロファイルとＤｓｔプロファイルを用いて四面体補間等の補間演算によって算出してもよいし、入力色域境界と出力色域境界を基に色差が最小となるように算出してもよい。算出された標本点の出力色信号値は、図６に示すように、出力色域の内部および出力色域境界近傍の色信号値になる。
【００３１】
ステップ７０５では、ステップ７０２で設定した標本点の入力色信号全てに関して出力色信号が算出されたかどうか判定する。全ての標本点の出力色信号が算出されていれば処理を終了し、全ての標本点の出力色信号が算出されていなければステップ７０４に移行して処理を継続する。
【００３２】
以下では、ステップ５０２における媒介変数の設定について、図８を用いて説明する。
【００３３】
まず、一般的な一次元のスプライン関数について説明する。図８（ａ）に示す一次元のスプライン関数は、入力信号ｘ、出力信号ｙに対して次式になる。
【００３４】
【数２】
2008294931.tif 000004
【００３５】
ここで、Ｎは標本点（入力信号：ｘ＝｛３，４，７，１１，１３｝）の数、αはスプライン係数、Ｂは基底関数である。
【００３６】
次に図８（ｂ）に示すように、入力信号を媒介変数ｔ＝｛０，１，２，３，４｝に置き換えるとスプライン関数は次式になる。
【００３７】
【数３】
2008294931.tif 000005
【００３８】
以上のように、図８（ａ）のような不等間隔な標本点の入力信号に対して、図８（ｂ）のように、整数値等を用いた等間隔の媒介変数を対応させる。つまり、入力色信号を媒介変数に変換するためのパラメータを求める。これにより、ステップ５０１において、標本点の入力色信号を入力色域の形状に合わせて設定しても、標本点を格子状や円筒状に設定していた従来法と同様の演算を適用することが可能となる。
【００３９】
そこで、上述の一次元のスプライン関数を三次元色空間に当てはまるように拡張する。ｓ番目の標本点の入力色信号値（ｌｓ，ｃｓ，ｈｓ）が明度方向ｘ番目の標本点ｌｘ（＝ｌｓ）、彩度方向ｙ番目の標本点ｃｙ（＝ｃｓ）、色相方向ｚ番目の標本点ｈｚ（＝ｈｚ）からなるとすると、出力色信号値（Ｌｓ，ａｓ，ｂｓ）を表すスプライン関数は次式になる。
【００４０】
【数４】
2008294931.tif 000006
【００４１】
【数５】
2008294931.tif 000007
【００４２】
【数６】
2008294931.tif 000008
【００４３】
ここで、Ｎ１、Ｎ２、Ｎ３は各次元の標本点の数、α、β、γは各次元のスプライン係数、Ｂは非周期的な基底関数、Ｂｃは周期的な基底関数である。ここで、入力色信号値を媒介変数（ｕｘ，ｖｙ，ｗｚ）に置き換えるとスプライン関数は次式になる。
【００４４】
【数７】
2008294931.tif 000009
【００４５】
【数８】
2008294931.tif 000010
【００４６】
【数９】
2008294931.tif 000011
【００４７】
以下では、ステップ５０３におけるスプライン係数の設定について、図１２のフローチャートを用いて説明する。
【００４８】
ステップ１２０１では、スプライン関数の節点を適当となるように設定する。節点は媒介変数で表され、シェーンバーグ・ホイットニの条件を満たす限り、その値は任意に決定してかまわない。例えば、標本点と同じ値に設定しても良いし、連続する標本点の中点の値に設定しても良い。但し、周期スプラインの場合は、節点の始点（媒介変数が最小の点）と終点（媒介変数が最大の点）は、周期の始点と終点と一致する必要がある。即ち、図９（ａ）に示すように、標本点の入力信号値の始点と終点が周期の始点と終点と一致している場合は、節点の始点と終点も標本点の始点と終点と一致する。また、図９（ｂ）や図９（ｃ）に示すように、標本点の入力信号値の始点あるいは終点もしくはその両方が周期の始点と終点と異なる場合の節点の始点と終点は、標本点の媒介変数の範囲を超えて周期の始点と終点と一致する媒介変数を設定する必要がある。
【００４９】
ステップ１２０２では、ステップ１２０１で設定した節点に基づいて基底関数を算出して係数行列を生成する。基底関数の算出には一般的なＤｅｂｏｏｒ−Ｃｏｘのアルゴリズムを用いる。算出する基底関数が三位のＢ−Ｓｐｌｉｎｅの基底とすると、４点の節点（ｑ０，ｑ１，ｑ２，ｑ３）で生成される非周期Ｂ−Ｓｐｌｉｎｅ基底は図１０のようになる。同様に、５点の節点（ｑ０，ｑ１，ｑ２，ｑ３，ｑ４）で生成される周期Ｂ−Ｓｐｌｉｎｅ基底は図１１のようになる。算出した基底関数にＮ点の標本点の媒介変数を代入することでＮ行Ｎ列の係数行列Ｍが生成される。係数行列Ｍのｓ行ｔ列成分ｍｓｔは次式になる。
【００５０】
【数１０】
2008294931.tif 000012
【００５１】
従って、ｓ番目の標本点の出力色信号値（Ｌｓ，ａｓ，ｂｓ）を算出するスプライン関数は次式になる。
【００５２】
【数１１】
2008294931.tif 000013
【００５３】
【数１２】
2008294931.tif 000014
【００５４】
【数１３】
2008294931.tif 000015
【００５５】
さらに、上式を行列式に変換すると次式になる。
Ｌ＝Ｍα
Ｃ＝Ｍβ
Ｈ＝Ｍγ
ここで、Ｌ、Ｃ、Ｈは各次元のＮ次出力色信号値ベクトル、α、β、γは各次元のＮ次スプライン係数ベクトルである。
【００５６】
ステップ１２０３では、ステップ１２０２で算出した行列式をＬＵ分解等で解き、スプライン係数ベクトルを算出する。
α＝Ｍ−１Ｌ＝［α０，α１・・・，αＮ１＊Ｎ２＊Ｎ３−１］Ｔ
β＝Ｍ−１Ｃ＝［β０，β１・・・，βＮ１＊Ｎ２＊Ｎ３−１］Ｔ
γ＝Ｍ−１Ｈ＝［γ０，γ１・・・，γＮ１＊Ｎ２＊Ｎ３−１］Ｔ
ステップ１２０４では、ステップ１２０３で算出したスプライン係数値をスプライン関数に設定する。
【００５７】
以下では、ステップ５０４の画像の色変換について、図１３のフローチャートを用いて説明する。
【００５８】
ステップ１３０１では、入力された画像より、あるピクセルのＲＧＢ値を読み込む。
【００５９】
ステップ１３０２では、ステップ１３０１で読み込んだＲＧＢ値をステップ４０２で読み込んだＳｒｃプロファイルを用いてＬ＊ａ＊ｂ＊値に変換すると共に彩度Ｃおよび色相Ｈを算出する。
【００６０】
ステップ１３０３では、ステップ５０２で標本点の入力色信号値に設定した媒介変数に基づいて、ステップ１３０２で入力した色信号値を媒介変数に変換する。一次元のある入力信号ｘ＝｛ｘ０，ｘ１，ｘ２，ｘ３，ｘ４｝に対して、媒介変数ｔ＝｛ｔ０，ｔ１，ｔ２，ｔ３，ｔ４｝が対応しており、
【００６１】
【数１４】
2008294931.tif 000016
【００６２】
が入力されたとすると、その媒介変数ｔａは次式になる。
【００６３】
【数１５】
2008294931.tif 000017
【００６４】
同様に、ステップ１３０２で入力した色信号値が（ｌａ，ａａ，ｂａ）がそれぞれ
【００６５】
【数１６】
2008294931.tif 000018
【００６６】
の区間に存在するものとすると対応する媒介変数（ｕａ，ｖａ，ｗａ）は次式になる。
【００６７】
【数１７】
2008294931.tif 000019
【００６８】
【数１８】
2008294931.tif 000020
【００６９】
【数１９】
2008294931.tif 000021
【００７０】
ステップ１３０４では、ステップ１２０４で設定したスプライン係数に基づくスプライン関数にステップ１３０３で算出した媒介変数値を代入して色を変換する。すなわち、スプライン係数と入力色信号の表色系の次元ごとに算出された基底関数の積とに基づく多重総和演算を用いて色を変換する。ステップ１２０４で算出した明度、彩度、色相それぞれの成分のスプライン係数をα、β、γ、ステップ１３０３で算出した媒介変数値を（ｕ，ｖ，ｗ）とすると、出力色信号値（Ｌ＊，ａ＊，ｂ＊）は次式になる。
【００７１】
【数２０】
2008294931.tif 000022
【００７２】
【数２１】
2008294931.tif 000023
【００７３】
【数２２】
2008294931.tif 000024
【００７４】
ステップ１３０５では、ステップ１３０４で算出されたＬ＊ａ＊ｂ＊値を、ステップ４０３で読み込んだＤｓｔプロファイルを用いてＣＭＹＫ値に変換する。
【００７５】
ステップ１３０６では、入力された画像の全ピクセルの色が変換されたかどうか確認する。全ピクセルの色変換が終了していれば処理を終了し、色変換が終了していなければステップ１３０１に移行して処理を継続する。
【００７６】
以上で説明した方法によってＲＧＢ値で表される画像がＣＭＹＫ値で表される画像に変換される。
【００７７】
以上の構成によれば、非線型の色変換式を算出するための標本点を色域形状に合わせて設定することで、低彩度色から高彩度色までの全色域に渡って高精度に色を算出することが可能となる。また、色空間に分布する標本点の変化の傾向に従って、標本点を滑らかに接続する非線型の色変換式を用いて色を算出するため、階調性の良好な色再現を得ることが可能となる。
【００７８】
（第２実施例）
第１実施例では、標本点の入力色信号値を補間するスプライン関数を用いて色変換する例を示した。しかしながら、より階調性の良好な色再現を得るためには補間関数よりも滑らかな関数を用いたほうが効果的である。補間関数の場合、図１４（ａ）に示すような関数が生成される可能性がある。このような関数を用いて色変換した場合、標本点が密な部分で色の変動が激しくなり画像の階調悪化が懸念される。それに対して、図１４（ｂ）のような関数であれば、階調が良好に再現できると考えられる。そこで、図１４（ｂ）のような関数を生成するためには、上記ステップ５０３においてスプライン係数の設定方法を変更すればよい。本実施例では、第一実施例で示したスプライン関数よりも滑らかな関数を生成するスプライン係数の設定例を示す。尚、ステップ５０３以外の処理については、第一実施例と同様とする。
【００７９】
以下では、第一実施例で示したステップ５０３におけるスプライン係数設定について、図１５のフローチャートを用いて説明する。
【００８０】
ステップ１５０１では、ステップ１２０１と同様にスプライン関数の節点を設定する。
【００８１】
ステップ１５０２では、ステップ１２０２と同様に基底関数を算出する。但し、節点の設定および基底関数の算出は、以降の色変換処理のために設定するものであり、本実施例におけるスプライン係数の設定とは無関係の処理である。
【００８２】
ステップ１５０３では、ステップ５０１で設定した標本点の出力色信号値をスプライン係数として設定する。従って、ｓ番目の標本点の出力色信号値を（ＬＳ，ａｓ，ｂＳ）とするとｓ番目のスプライン係数値（αＳ，βｓ，γＳ）は次式になる。
αＳ＝ＬＳ
βｓ＝ａｓ
γＳ＝ｂＳ
以上の手法で設定したスプライン係数を用いることによって、より滑らかに変化するスプライン関数を得ることが可能であり、階調性の良好な色変換が可能になる。
【００８３】
上記実施例では色変換を実行する色空間としてＣＩＥＬＡＢ色空間を用いた。しかしながら、色変換を実行する色空間には何ら制限はなくＲＧＢ色空間、ＣＭＹＫ色空間、ＪＣｈ色空間、ＸＹＺ色空間、ｘｙ色空間等を用いることが可能である。さらに、インク色分解などにも適用可能であり、ＣＭＹＫに淡シアンならびに淡マゼンタを加えた６次元色空間等も用いることができる。また、式に与える色空間と、式より算出される色空間は実施例のように同一であっても良いし、異なってもよい。
【図面の簡単な説明】
【００８４】
【図１】色処理装置のブロック図である。
【図２】色処理アプリケーションのインタフェースを表す図である。
【図３】プロファイルのＬＵＴのフォーマットを表す図である
【図４】色算出アプリケーション動作を表すフローチャート図である。
【図５】色変換処理のフローチャート図である。
【図６】標本点の設定方法を示す模式図である。
【図７】標本点設定方法のフローチャート図である。
【図８】入力信号とスプライン関数の関係を示す模式図である。
【図９】媒介変数とスプライン関数の関係を示す模式図である。
【図１０】非周期的なスプライン基底関数を表す図である。
【図１１】周期的なスプライン基底関数を表す図である。
【図１２】スプライン係数設定のフローチャート図である。
【図１３】画像の色変換処理のフローチャート図である。
【図１４】スプライン係数の異なるスプライン関数の模式図である。
【図１５】スプライン係数設定のフローチャート図である。
【図１６】四面体補間における四面体の生成方法を示す模式図である。
【図１７】ＲＧＢ色空間とＣＩＥＬＡＢ色空間の対応を示す模式図である。
【図１８】従来法の標本点の設定方法を示す模式図である。
【図１】
2008294931.tif 000025
【図２】
2008294931.tif 000026
【図３】
2008294931.tif 000027
【図４】
2008294931.tif 000028
【図５】
2008294931.tif 000029
【図６】
2008294931.tif 000030
【図７】
2008294931.tif 000031
【図８】
2008294931.tif 000032
【図９】
2008294931.tif 000033
【図１０】
2008294931.tif 000034
【図１１】
2008294931.tif 000035
【図１２】
2008294931.tif 000036
【図１３】
2008294931.tif 000037
【図１４】
2008294931.tif 000038
【図１５】
2008294931.tif 000039
【図１６】
2008294931.tif 000040
【図１７】
2008294931.tif 000041
【図１８】
2008294931.tif 000042