JP3954663B2

JP3954663B2 - Color signal conversion method

Info

Publication number: JP3954663B2
Application number: JP28489295A
Authority: JP
Inventors: 伸青木
Original assignee: Ricoh Co Ltd
Current assignee: Ricoh Co Ltd
Priority date: 1995-11-01
Filing date: 1995-11-01
Publication date: 2007-08-08
Anticipated expiration: 2015-11-01
Also published as: JPH09130626A

Description

【０００１】
【発明の属する技術分野】
本発明は、カラー画像信号を、属性の異なるカラー画像信号に変換する色信号変換方法に関する。
【０００２】
【従来の技術】
カラー画像入出力システムにおいては色変換などの技術が不可欠である。ここで色変換とは、カラー画像信号が色を表す表色空間（ＲＧＢ、ＣＭＹＫ、Ｌ＊ａ＊ｂ＊など）を変換する処理のことであり、例えば、ＮＴＳＣ−ＲＧＢで表現されたカラー画像信号をＣＩＥ−Ｌ＊ａ＊ｂ＊信号などに変換する処理をいう。
【０００３】
従来の色変換法としては、多数のサンプル出力を測色し、その入出力関係を多項式で近似する方法が知られている。例えば、ｃｍｙｋの多項式により、出力ＸＹＺを予測するもので、多項式の係数は実測値との最小２乗近似によって決定する再現色予測方法がある（特開平４−３３７９６５号公報を参照）。
【０００４】
この多項式近似による色変換法では、多項式の次数を増やすことで近似式の自由度が増し、カラープリンタなど非線形性の強い変換に対してもより精度よく近似を行うことが可能となる。しかし、一般に多項式近似は次数が高くなると不要な振動を起こしやすくなるという性質がある。図４は、その例を示す。図では簡単のために横軸が入力信号、縦軸が出力信号の１次元の信号変換を示している。図中（ａ）は、サンプル出力の測定値であり、（ｂ）、（ｃ）はそれぞれ測定値を最小２乗近似する３次、４次多項式曲線である。３次の曲線は滑らかであるが、４次の曲線は波打っていることが分かる。
【０００５】
【発明が解決しようとする課題】
近似多項式が振動すると、例えば、なだらかな濃度変化のある画像を変換すると、縞状の濃淡模様が生じてしまうことになり、これはたとえ平均的な色差が小さくても画質上、大きな問題となる。
【０００６】
また、別の方法として格子点補間法がある。これは入力空間を格子状に分割し、格子の頂点での出力値を予め実測しておく。そして格子内部の入力に対しては、頂点の値を使って、補間によって出力値を計算する方法である（例えば、特開平２−２２６８６７号公報を参照）。この方法では、各格子の内部は多項式により補間することになるが、空間を分割しているため、各領域内では低次の多項式を使用しても精度が保てる。
【０００７】
しかし、各格子点での出力値を補間する多項式はその形式が限られていて、任意の次数の項を含む多項式を使用することができなかった。そのため、実現できる変換特性の自由度が限られ、色変換の精度向上に限界があった。
【０００８】
本発明の目的は、高次の多項式を使用することによる振動を回避しながら、非線形性の強いプリンタに対しても精度良く色変換特性を近似するための多くの自由度を持つ色信号変換方法を提供することにある。
【０００９】
【課題を解決するための手段】
前記目的を達成するために、請求項１記載の発明では、２次多項式を用いて第１のカラー画像信号を、属性の異なる第２のカラー画像信号に変換する色信号変換方法において、前記第１のカラー画像信号が黒信号を含むとき、前記第１のカラー画像信号から黒信号を減算した信号からなる３次元空間を、前記黒信号を減算した各信号の値が０となる平面を境界とする複数の領域に分割し、該複数の領域毎に異なる係数を用い、かつ領域の境界で１次微分値が連続となる係数を用いた前記多項式によって、前記第１のカラー画像信号を第２のカラー画像信号に変換することを特徴としている。
【００１０】
【発明の実施の形態】
上記課題で説明したように区分的な多項式を使用した色変換方法では、領域の境界で連続となる多項式係数であっても、その微分値が連続になるとは限らない。ここで、区分的な多項式とは、例えばｘとｙの関係が、ｘ＞０の領域ではｙ＝ｘの多項式で表され、ｘ≦０の領域ではｙ＝ｘ²の多項式で表されるように、複数の多項式で近似したものをいう。
【００１１】
微分値が不連続ということは、境界で変換特性が急激に変化する、つまり「角」を持つことを意味する。一般にプリンタなどの色変換特性は、なめらかな変化をするものが多いので、このような急激な変化は好ましくない。
【００１２】
本発明では、領域の境界で多項式の微分値を連続にすることで、なめらかな色変換特性を持つようにする。
【００１３】
ところで、カラー画像はｃ，ｍ，ｙ，ｋのインクやトナーを重ねることによって再現されるが、この再現方法として、各色毎に異なるスクリーン角を持たせて重ねる方法と、各色を同心円的に重ねる方法（同心円モデル）がある。
【００１４】
理想的な同心円モデルのプリンタでは、各色材の面積の大小関係によって場合分けすれば、ｃｍｙｋの面積率からＸＹＺへの色変換が、それぞれの場合で線形になることが知られている。例えばｃ＞ｍ＞ｙ＞ｋのときｃとＸの関係がリニアになる。しかし、実際のプリンタは、各版の位置ずれ、ドットの形状のゆがみ、紙内部の光散乱などのため、理想的な同心円モデルとは異なる特性を持つ。そこで、
本発明では、理想的な同心円モデルでは線形となるべき各領域（つまり、上記したｃ＞ｍ＞ｙ＞ｋの領域など）を、それぞれ異なる多項式で近似することにより、より現実のプリンタにあった色信号の変換特性となる。
【００１５】
多項式の係数を最小２乗法などで設定するためには、各領域について多数のサンプルデータが必要となる。サンプルデータの少ない領域については、正確な係数を設定することが難しい。そこで色変換式としては、できるだけ分割する領域の数が少ないことが望ましい。
【００１６】
ｃｍｙｋ４種の色材を使用するプリンタでは、全ての大小関係を考慮して分類すると、４つの要素の順列は４×３×２×１＝２４通りあるので、２４個の領域に分割される。
【００１７】
しかし、全ての組み合わせが同じ重要性を持つわけではない。例えば、シアンと黒、およびシアンとマゼンタの関係を比べてみる。図５は、理想的な同心円モデルで、一定のマゼンタ（ｍ）に対し、シアン（ｃ）が小さい場合（ｃ＜ｍ）、等しい場合（ｃ＝ｍ）、シアンが大きい場合（ｃ＞ｍ）の様子を示している。マゼンタよりもシアンが小さい場合は、マゼンタと（マゼンタとシアンの重なりである）ブルーの２つの領域からなり、マゼンタとシアンが同じ大きさになると全てブルーになり、さらにシアンが増えるとシアンとブルーの領域ができる。
【００１８】
このようにシアンが増えるに従い、ｃ＜ｍの領域ではマゼンタがブルーに置き換わることで、また、ｃ＞ｍの領域では白がシアンに置き換わることで、全体として色が変化する。そのため、一定の大きさのマゼンタの上にシアンを重ねた時の全体としての色度（Ｘ）変化は、図５の下のグラフのように（ｃ＝ｍ）を境に異なる傾きを持つ。
【００１９】
一方、図６は、図５と同様に一定の黒に、３通りの大きさのシアンが重なった様子を示している。この場合、シアンが黒より小さいときは、黒は全波長域を吸収するため黒の上にシアンがのってもほぼ黒のままなので、シアンの量によらずほぼ一定の色を表し、シアンが黒より大きくなって初めてシアン色材の量が全体として色に大きく影響を及ぼすようになる。
【００２０】
そのため一定の大きさの黒の上にシアンを重ねた時の全体としての色度（Ｘ）変化は、図６の下のグラフのように、ｃ＜ｋでは一定（黒）のままであり、ｃ＞ｋになって初めてｃの量が全体の色に影響を与えるようになる。
【００２１】
ここではシアンが、マゼンタと黒に重なった例を挙げたが、シアン、マゼンタを他の色材に変えてもほぼ同じことが起こる。これは黒色材が全ての波長の光を吸収するために、他のどの色材とも吸収帯域が重なるためである。
【００２２】
このように、各カラー信号の間の大小関係の中で、特にその境界で色再現特性の振舞が変わるのは黒とその他の色の関係である。そのため、区分的多項式を構成するために領域を分割するならば、黒と他の色との大小関係が特に重要である。つまり、図５のグラフでは、境界（ｃ＝ｍ）の左右では色変換特性の違いが少ないので、全体を一つの多項式で近似することができるが、図６のグラフのようにｃ＝ｋを境界に二つの領域で大きく異なる色変換特性を持つ場合は、異なる多項式を使うべきである。
【００２３】
本発明では、各信号の関係の内、その大小関係により、色再現特性の大きく変化する部分だけを分割することにより、不必要な領域数の増加を抑え、少数のサンプルで正確に係数を設定することができる。
【００２４】
上記したように、黒信号（ｋ）と他のカラー信号（ｃ，ｍ，ｙ）との大小関係で領域を分割する場合、境界は、
ｃ＝ｋ
ｍ＝ｋ
ｙ＝ｋ
となる。ｃ，ｍ，ｙ，ｋそのままの４次元空間では、これはどの軸とも直交しない、いわば斜めの超平面であり、各境界での連続性の扱いが難しい。図７は、例えばｃとｋの関係を示し、ｃ＝ｋが斜めの超平面となる。
【００２５】
本発明では、座標系を変換して各領域の境界が各軸と直交する超平面とすることで、連続性の保証を簡単にしている。すなわち、座標系（ｃ，ｋ）を、座標系（ｃ−ｋ，ｋ）に変換する。図８は、ｃ−ｋとｋの関係を示し、ｃ−ｋ＝０（ｃ＝ｋ）の超平面は、（ｃ−ｋ）軸に直交する。以下、本発明の一実施例を図面を用いて具体的に説明する。
本発明は、色変換式として区分多項式を使用し、各領域毎の多項式係数を少数のサンプルの実測値に対する最小２乗近似によって決定するものである。以下の実施例ではＣＭＹＫをＬ＊ａ＊ｂ＊に変換する変換法について説明するが、ｒｇｂをＬｕｖに変換するなど、その他の属性、次元数の変換についても同様の方法で実施することができる。
【００２６】
〈実施例〉
（色変換式）
図１は、本発明の実施例の構成を示す。本実施例では、ｃ，ｍ，ｙ，ｋそれぞれの面積率（０〜１）を入力とし、実際にプリンタでそれを出力し、測色器で測定したときのＬ＊ａ＊ｂ＊値に変換するものとする。図において、１はｃ，ｍ，ｙからｋを引く減算器、２はｃ−ｋ，ｍ−ｋ，ｙ−ｋの大小関係を基に領域を判別する領域判別部、３は領域毎に設定された係数を記憶する係数記憶部、４は判定された領域の係数を使った多項式によってＬ＊ａ＊ｂ＊を出力する多項式演算部である。
【００２７】
まず、減算器１は、座標系を変換するために、ｃ，ｍ，ｙそれぞれからｋを引き、ｃ−ｋ（＝ｃ’）、ｍ−ｋ（＝ｍ’）、ｙ−ｋ（＝ｙ’）、ｋの４次元信号を作り、その内ｃ’，ｍ’，ｙ’の３次元空間をｃ’＝０、ｍ’＝０、ｙ’＝０を境界として分割する。その結果として、それぞれの入力信号は、図２に示すように８つの領域に判定される。図３は、図２で判定される各領域を示す。
【００２８】
このような分割は、ｃ’＝ｃ−ｋ＝０より、元々のｃ，ｍ，ｙ，ｋの４次元空間では、ｃ＝ｋ，ｍ＝ｋ，ｙ＝ｋの３つの超平面を境界に分割する、つまりｋとｃ，ｍ，ｙとの大小関係により分割することに相当する。
【００２９】
次に、ｃ’，ｍ’，ｙ’，ｋの４つの信号から、後述するように各領域で予め決められた係数を使った多項式によって、Ｌ＊，ａ＊，ｂ＊の３つを計算する。多項式としては、ｃ’，ｍ’，ｙ’，ｋの双２次形式、つまり、
【００３０】
【数１】

【００３１】
のように各項が、
【００３２】
【数２】

【００３３】
の各行から１つづつ要素を取り出して掛け合わせたものとする。各変数について３通りあるので、全ての組み合わせは３⁴＝８１通りあり、次のような８１項からなる多項式になる。
【００３４】
【数３】

【００３５】
ただし、ａ_iは、後述する方法で係数記憶部３に設定される各領域毎の係数である。使用する多項式としては、例えば上記した８１項の内からｍ項だけを選択して、一部分だけを使うなど、異なる形式とすることもできる。
【００３６】
（多項式係数の設定）
測定データから最小２乗法により各領域での係数を設定する。ここで各領域の多項式係数には以下の性質が要求される。すなわち、
（１）サンプルデータ（ｃ，ｍ，ｙ，ｋ−Ｌ＊ａ＊ｂ＊）を良く近似する。
近似式の精度が当然求められる。
（２）領域の境界で連続である。
本来プリンタの色変換特性は、連続的な関数と考えられるので、近似の都合で不連続点が生じてはならない。
【００３７】
そこで、各領域での変換の連続性を保証するため、次のような方法をとる。まず、簡単のために１次元の信号について説明し、次いでｃｍｙｋ４次元の場合を説明する。
【００３８】
ｃ’の１次元信号について、
【００３９】
【数４】

【００４０】
の４つの項を考える。ただし、
【００４１】
【数５】

【００４２】
とする（なお、（式３）は切断ベキ関数という）。
【００４３】
この４つから、任意の係数（ａ０，ａ１，ａ２，ａ３）を使った多項式を作ると、
【００４４】
【数６】

【００４５】
となる（１次元の区分多項式）。これは、ｃ’が０以上の場合と０未満の場合に分けて考えれば、それぞれが
【００４６】
【数７】

【００４７】
なる２次多項式である。また、その境界（ｃ’＝０）の前後で異なるのは、２次の項の係数（ａ２，ａ３）だけであり、ｃ’²は境界点で０となるので、これら２つの２次多項式は境界で滑らかにつながる。正確にはその値および１次微分値が連続となる。
【００４８】
従って、式（４）は、境界で滑らかに接続する区分多項式であるので、係数（ａ０からａ３）を、測定値との平均２乗誤差が最小となるように決めればよい。そして、領域（ｃ’＜０）には係数ａ０，ａ１，ａ２が設定され、領域（ｃ’≧０）には係数ａ０，ａ１，ａ３が設定される。
【００４９】
ｃ’ｍ’ｙ’ｋ４次元の場合も式（２）と同様に、ｍ’，ｙ’，ｋについても４つの項をとり、それぞれを掛け合わせる。つまり、
【００５０】
【数８】

【００５１】
の各行から１つづつを掛け合わせたものとする（ただし、（ｋ＞０）なのでｋ²＿は必要がなく、ｋについて３項となる）。そうすると、４³×３＝１９２項の式ができる。
【００５２】
この多項式の第ｉ項をｐ_i（ｃ，ｍ，ｙ，ｋ）とし、係数ｑ_iを使って、
【００５３】
【数９】

【００５４】
を考える。これは１次元の場合の式（４）に対応するものであり、任意の係数（ｑ_i）に対して、図２に示す各領域で双２次式であり、また境界（ｃ＝ｋ，ｍ＝ｋ，ｙ＝ｋ）で関数値および１次微分値が連続な関数となる。
【００５５】
上記した１９２項は、各領域の多項式として、式（１）を使用する場合に対応するもの、つまり式（１）のうち、ｘ²を含む項（ｘはｃ，ｍ，ｙ，ｋのどれか)を
【００５６】
【数１０】

【００５７】
と
【００５８】
【数１１】

【００５９】
で置き換えたものである。各領域の多項式として、式（１）以外の形式を使用する場合にも、同様な置き換えをすればよい。
【００６０】
この係数（ｑ_i）は、多数のパッチの測定値から最小２乗法で決める。つまり、まず多数（Ｎ組）の（ｃ，ｍ，ｙ，ｋ）の組み合わせでプリンタ出力を行い、出力を測定し、測色値（Ｌ＊_n，ａ＊_n，ｂ＊_n）を得る。また、出力したそれぞれのパッチについて、（ｃ，ｍ，ｙ，ｋ）からｐ_i（ｃ，ｍ，ｙ，ｋ）を求める。第ｎパッチについての第ｉ項をＰ_i，_nとすると、
【００６１】
【数１２】

【００６２】
となる近似式が得られる。式（５）の第１の行列の各行がｐ_i（ｃ，ｍ，ｙ，ｋ)であり、第２の行列の各列がｑ_iであり、式（５）の右辺の行列の各行が測定値である。これは係数（ｑ_i）についての線形式であり、周知の線形最小２乗法で係数を求めることができる。
【００６３】
係数を求めた後、ｃｍｙｋからＬ＊ａ＊ｂ＊を計算するときには、各領域毎に３⁴＝８１項の双２次式を使う。例えば、図２、３で、領域０に属するｃ，ｍ，ｙ，ｋの入力に対しては、ｃ’≧０，ｍ’≧０，ｙ’≧０’より、
【００６４】
【数１３】

【００６５】
は必ず０となるので、これらを含む項は０となる。従って、（ｃ’，ｍ’，ｙ’，ｋ）がそれぞれ、
【００６６】
【数１４】

【００６７】
３種類の組み合わせからなる項の係数だけを選んで使用すればよい。各領域についても同様に、
【００６８】
【数１５】

【００６９】
のうち０とならない方だけを選べば、８１項の多項式（つまり、８１個のｑ_iを領域毎に選んだ多項式）となる。
【００７０】
このように予め係数（ｑ_i）から、領域毎に必要な係数を選択し、各領域毎に８１個の係数（ａｉ）を係数記憶部３に記憶する。そして、入力される信号（ｃ，ｍ，ｙ，ｋ）の領域を領域判別部２で判別し、その領域の係数ａｉを係数記憶部３から読み出し、多項式演算部４に与える。多項式演算部４は読み出された係数ａｉを使って、式（１）の多項式を計算し、Ｌ＊ａ＊ｂ＊を出力する。
【００７１】
【発明の効果】
以上、説明したように、本発明によれば、入力空間を複数の領域に分割し、各領域毎に異なる係数を使った多項式により色変換を行うので、各領域毎の多項式としては低次の式を用い、高次多項式に起こりがちな不要な振動を抑えることができる。また、領域の境界で多項式の微分値が連続となるので、滑らかな特性の色変換方法を実現することができる。
【００７２】
本発明によれば、任意の次数の項を含む多項式を用いるので、各領域内での自由度を利用して、高精度な色変換を実現することができる。
【００７３】
本発明によれば、各カラー信号の大小関係を基に領域を分割し、同心円モデルに近いプリンタに対しては、各領域での色変換特性を多項式近似することによって、そのプリンタの理想的な特性とのずれを吸収し、精度よい近似を行うことができる。
【００７４】
本発明によれば、黒信号と他のカラー信号との大小関係を基に領域を分割するので、近似多項式の自由度を保ちながら領域数を抑えることができる。
【００７５】
本発明によれば、黒以外のカラー信号から黒信号を減算した信号を作り、それら黒減算信号と黒信号を、多項式によって色変換するので、黒と他の色との大小関係による領域分割の境界で、近似多項式が連続となる係数を簡単に設定することができる。
【図面の簡単な説明】
【図１】本発明の実施例の構成を示す。
【図２】入力信号の判定結果を示す。
【図３】図２で判定される各領域を示す。
【図４】多項式近似の問題点を説明する図である。
【図５】同心円モデルにおけるシアンとマゼンタの関係を説明する図である。
【図６】同心円モデルにおけるシアンと黒の関係を説明する図である。
【図７】ｃとｋの関係を示す図である。
【図８】ｃ−ｋとｋの関係を示す図である。
【符号の説明】
１減算器
２領域判別部
３係数記憶部
４多項式演算部[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a color signal conversion method for converting a color image signal into a color image signal having different attributes.
[0002]
[Prior art]
In color image input / output systems, techniques such as color conversion are indispensable. Here, the color conversion is a process of converting a color space (RGB, CMYK, L * a * b *, etc.) in which a color image signal represents a color. For example, a color image expressed in NTSC-RGB. A process for converting a signal into a CIE-L * a * b * signal or the like.
[0003]
As a conventional color conversion method, a method of measuring a number of sample outputs and approximating the input / output relationship with a polynomial is known. For example, there is a reproduction color prediction method in which the output XYZ is predicted by a cmyk polynomial, and the coefficient of the polynomial is determined by least square approximation with an actual measurement value (see Japanese Patent Laid-Open No. 4-337965).
[0004]
In this color conversion method by polynomial approximation, the degree of freedom of the approximate expression increases by increasing the degree of the polynomial, and it is possible to perform approximation with higher accuracy even for conversion with strong nonlinearity such as a color printer. However, in general, polynomial approximation has the property that unnecessary vibrations are likely to occur as the order increases. FIG. 4 shows an example. In the figure, for the sake of simplicity, the horizontal axis represents the one-dimensional signal conversion of the input signal and the vertical axis represents the output signal. In the figure, (a) is a measured value of the sample output, and (b) and (c) are third-order and fourth-order polynomial curves that respectively approximate the measured value to the least square. It can be seen that the cubic curve is smooth, but the quartic curve is wavy.
[0005]
[Problems to be solved by the invention]
When the approximate polynomial vibrates, for example, when an image having a gentle density change is converted, a striped shading pattern is generated. This is a serious problem in image quality even if the average color difference is small. .
[0006]
Another method is a lattice point interpolation method. In this method, the input space is divided into grids, and the output values at the vertices of the grid are measured in advance. For the input inside the grid, the output value is calculated by interpolation using the value of the vertex (see, for example, Japanese Patent Laid-Open No. 2-226867). In this method, the interior of each lattice is interpolated by a polynomial, but since the space is divided, the accuracy can be maintained even if a low-order polynomial is used in each region.
[0007]
However, the form of the polynomial that interpolates the output value at each grid point is limited, and a polynomial including an arbitrary-order term cannot be used. For this reason, the degree of freedom of conversion characteristics that can be realized is limited, and there is a limit to improving the accuracy of color conversion.
[0008]
An object of the present invention is to provide a color signal conversion method having many degrees of freedom for accurately approximating a color conversion characteristic even for a printer with strong nonlinearity while avoiding vibration caused by using a higher-order polynomial. Is to provide.
[0009]
[Means for Solving the Problems]
In order to achieve the above object, according to the first aspect of the present invention, in the color signal conversion method for converting the first color image signal to the second color image signal having a different attribute using a second-order polynomial, When one color image signal includes a black signal, a three-dimensional space formed by subtracting the black signal from the first color image signal is bounded by a plane where the value of each signal obtained by subtracting the black signal is zero into a plurality of regions and, by the polynomial with coefficients different for each region of the plurality of, and the primary differential value at the boundary of the region with a continuous coefficient, the first color image signal the It is characterized in that it is converted into two color image signals.
[0010]
DETAILED DESCRIPTION OF THE INVENTION
As described in the above problem, in the color conversion method using the piecewise polynomial, even if the polynomial coefficient is continuous at the boundary of the region, the differential value is not always continuous. Here, the piecewise polynomial is such that, for example, the relationship between x and y is represented by a polynomial of y = x in the region of x> 0, and represented by a polynomial of y = x ^{2 in} the region of x ≦ 0. The one approximated by a plurality of polynomials.
[0011]
The fact that the differential value is discontinuous means that the conversion characteristic changes abruptly at the boundary, that is, has an “angle”. In general, color conversion characteristics of a printer or the like often change smoothly, and such a rapid change is not preferable.
[0012]
In the present invention, a smooth color conversion characteristic is obtained by making the differential value of the polynomial continuous at the boundary of the region.
[0013]
By the way, a color image is reproduced by superimposing c, m, y, and k inks and toners. As this reproduction method, a method of superimposing each color with different screen angles and a method of superimposing each color concentrically. There is a method (concentric circle model).
[0014]
In an ideal concentric circle model printer, it is known that the color conversion from the area ratio of cmyk to XYZ is linear in each case if classified according to the size relationship of the area of each color material. For example, when c>m>y> k, the relationship between c and X is linear. However, an actual printer has characteristics different from an ideal concentric model due to misalignment of each plate, distortion of dot shape, light scattering inside the paper, and the like. Therefore,
In the present invention, each region that should be linear in the ideal concentric circle model (that is, the region of c>m>y> k described above) is approximated by different polynomials, so that the printer can be more realistic. This is a color signal conversion characteristic.
[0015]
In order to set the coefficients of the polynomial by the least square method or the like, a large number of sample data is required for each region. It is difficult to set an accurate coefficient for an area with a small amount of sample data. Therefore, as a color conversion formula, it is desirable that the number of regions to be divided is as small as possible.
[0016]
In a printer using cmyk4 color materials, classification is performed in consideration of all the magnitude relations, and there are 4 × 3 × 2 × 1 = 24 permutations of the four elements, so that the area is divided into 24 areas.
[0017]
However, not all combinations have the same importance. For example, compare the relationship between cyan and black, and cyan and magenta. FIG. 5 shows an ideal concentric circle model, where cyan (c) is small (c <m), equal (c = m), and cyan is large (c> m) with respect to constant magenta (m). The state of is shown. If cyan is smaller than magenta, it consists of two areas, magenta and blue (which is the overlap of magenta and cyan). The area is made.
[0018]
As cyan increases, the color changes as a whole by replacing magenta with blue in the region of c <m, and replacing white with cyan in the region of c> m. Therefore, the chromaticity (X) change as a whole when cyan is superimposed on magenta having a certain size has different slopes with respect to (c = m) as shown in the lower graph of FIG.
[0019]
On the other hand, FIG. 6 shows a state in which three different sizes of cyan are superimposed on a certain black as in FIG. In this case, when cyan is smaller than black, black absorbs the entire wavelength range, so even if cyan is placed on black, it remains almost black. The amount of the cyan color material has a large influence on the color as a whole only when the color becomes larger than black.
[0020]
Therefore, the overall chromaticity (X) change when cyan is superimposed on black of a certain size remains constant (black) when c <k, as shown in the lower graph of FIG. Only when c> k, the amount of c affects the overall color.
[0021]
Here, an example is given in which cyan overlaps magenta and black, but the same thing happens if cyan and magenta are changed to other color materials. This is because the black material absorbs light of all wavelengths, so that the absorption band overlaps with any other color material.
[0022]
As described above, among the magnitude relationships between the color signals, it is the relationship between black and other colors that changes the behavior of the color reproduction characteristics particularly at the boundary. Therefore, the size relationship between black and other colors is particularly important if the region is divided to construct a piecewise polynomial. That is, in the graph of FIG. 5, there is little difference in color conversion characteristics on the left and right of the boundary (c = m), so the whole can be approximated by one polynomial, but c = k can be expressed as in the graph of FIG. Different polynomials should be used if the boundary has very different color conversion characteristics in the two regions.
[0023]
In the present invention, by dividing only the portion where the color reproduction characteristics greatly change due to the magnitude relationship among the relationships of each signal, the increase in the number of unnecessary areas is suppressed, and the coefficient is set accurately with a small number of samples. can do.
[0024]
As described above, when the region is divided by the magnitude relationship between the black signal (k) and the other color signals (c, m, y), the boundary is
c = k
m = k
y = k
It becomes. In a four-dimensional space with c, m, y, and k as they are, this is an oblique hyperplane that is not orthogonal to any axis, so it is difficult to handle continuity at each boundary. FIG. 7 shows the relationship between c and k, for example, where c = k is an oblique hyperplane.
[0025]
In the present invention, the guarantee of continuity is simplified by converting the coordinate system into a hyperplane where the boundary of each region is orthogonal to each axis. That is, the coordinate system (c, k) is converted to the coordinate system (c−k, k). FIG. 8 shows the relationship between ck and k, and the hyperplane of ck = 0 (c = k) is orthogonal to the (c−k) axis. Hereinafter, an embodiment of the present invention will be described in detail with reference to the drawings.
According to the present invention, a piecewise polynomial is used as a color conversion formula, and a polynomial coefficient for each region is determined by least square approximation with respect to actually measured values of a small number of samples. In the following embodiments, a conversion method for converting CMYK to L * a * b * will be described. However, conversion of other attributes and dimensions such as conversion of rgb to Luv can be performed in the same manner. .
[0026]
<Example>
(Color conversion formula)
FIG. 1 shows the configuration of an embodiment of the present invention. In this embodiment, the respective area ratios (0 to 1) of c, m, y, and k are input, and are actually output by a printer and converted into L * a * b * values when measured by a colorimeter. Shall be converted. In the figure, 1 is a subtracter that subtracts k from c, m, and y, 2 is a region discriminating unit that discriminates a region based on the magnitude relationship of ck, mk, and yk, and 3 is set for each region. A coefficient storage unit 4 for storing the determined coefficients is a polynomial calculation unit for outputting L * a * b * by a polynomial using the coefficients of the determined area.
[0027]
First, in order to convert the coordinate system, the subtracter 1 subtracts k from each of c, m, and y to obtain ck (= c ′), mk (= m ′), and yk (= y). '), A four-dimensional signal of k is created, and a three-dimensional space of c', m ', and y' among them is divided with c '= 0, m' = 0, and y '= 0 as boundaries. As a result, each input signal is determined as eight regions as shown in FIG. FIG. 3 shows each area determined in FIG.
[0028]
Such division is based on the three hyperplanes c = k, m = k, y = k in the original four-dimensional space of c, m, y, k from c ′ = c−k = 0. This corresponds to dividing, that is, dividing according to the magnitude relationship between k and c, m, and y.
[0029]
Next, three of L *, a *, and b * are calculated from the four signals c ′, m ′, y ′, and k using a polynomial that uses coefficients predetermined in each region as will be described later. To do. The polynomial is a biquadratic form of c ′, m ′, y ′, k, that is,
[0030]
[Expression 1]

[0031]
Each term like
[0032]
[Expression 2]

[0033]
It is assumed that one element is taken out from each row and multiplied. Since triplicate for each variable, all combinations There are 3 ⁴ = 81, the polynomial comprising a following 81 Section.
[0034]
[Equation 3]

[0035]
Here, a _i is a coefficient for each region set in the coefficient storage unit 3 by a method described later. As a polynomial to be used, for example, only m terms can be selected from the above 81 terms and only a part can be used.
[0036]
(Polynomial coefficient setting)
Coefficients in each region are set from the measurement data by the least square method. Here, the following properties are required for the polynomial coefficient of each region. That is,
(1) Approximate sample data (c, m, y, k-L * a * b *) well.
Naturally, the accuracy of the approximate expression is required.
(2) It is continuous at the boundary of the region.
Originally, the color conversion characteristics of a printer are considered to be a continuous function, so that discontinuities should not occur due to approximation.
[0037]
In order to guarantee the continuity of conversion in each region, the following method is taken. First, for the sake of simplicity, a one-dimensional signal will be described, and then a cmyk four-dimensional case will be described.
[0038]
For the one-dimensional signal of c ′
[0039]
[Expression 4]

[0040]
Consider the following four terms. However,
[0041]
[Equation 5]

[0042]
(Equation 3 is called a cutting power function).
[0043]
From these four, if a polynomial using arbitrary coefficients (a0, a1, a2, a3) is made,
[0044]
[Formula 6]

[0045]
(One-dimensional piecewise polynomial). This can be divided into a case where c ′ is 0 or more and a case where it is less than 0.
[Expression 7]

[0047]
Is a second order polynomial. Further, only the coefficients (a2, a3) of the second-order terms are different before and after the boundary (c ′ = 0), and c ′ ² is 0 at the boundary point. Is smoothly connected at the boundary. To be precise, the value and the first derivative value are continuous.
[0048]
Therefore, since Equation (4) is a piecewise polynomial that is smoothly connected at the boundary, the coefficients (a0 to a3) may be determined so that the mean square error with the measured value is minimized. Then, coefficients a0, a1, a2 are set in the area (c ′ <0), and coefficients a0, a1, a3 are set in the area (c ′ ≧ 0).
[0049]
In the case of c′m′y′k four-dimensional, similarly to the equation (2), m ′, y ′, and k take four terms and multiply them. That means
[0050]
[Equation 8]

[0051]
Are multiplied by one from each row (however, since (k> 0), k ² _ is not necessary and there are three terms for k). Then, the equation of 4 ³ × 3 = 192 terms is obtained.
[0052]
The i-th term of the polynomial _{p i (c, m, y} , k) and, using the coefficients q _i,
[0053]
[Equation 9]

[0054]
think of. This corresponds to equation (4) in the one-dimensional case, and for an arbitrary coefficient (q _i ), it is a biquadratic equation in each region shown in FIG. 2, and a boundary (c = k, m = k, y = k), the function value and the first derivative are continuous functions.
[0055]
The term 192 described above corresponds to the case where the expression (1) is used as the polynomial in each region, that is, the expression (1) including the term x ² (where x is c, m, y, k). )] [0056]
[Expression 10]

[0057]
And [0058]
[Expression 11]

[0059]
It has been replaced with. The same replacement may be made when a form other than the expression (1) is used as the polynomial in each region.
[0060]
This coefficient (q _i ) is determined by the least square method from the measured values of a large number of patches. That is, first, printer output is performed with a large number (N sets) of (c, m, y, k), the output is measured, and colorimetric values (L * _n , a * _n , b * _n ) are obtained. Also, for each of the patches output, determine the (c, m, y, k ) from _{p i (c, m, y} , k). Assuming that the i-th term for the n-th patch is _Pi , _n ,
[0061]
[Expression 12]

[0062]
An approximate expression is obtained. Each row of the first matrix of Equation (5) is p _i (c, m, y, k), each column of the second matrix is q _i , and each row of the right-hand side matrix of Equation (5) is It is a measured value. This is a linear format for the coefficient (q _i ), and the coefficient can be obtained by a well-known linear least square method.
[0063]
After calculating the coefficients, when calculating L * a * b * from cmyk, a biquadratic equation of 3 ⁴ = 81 terms is used for each region. For example, in FIGS. 2 and 3, for c, m, y, and k belonging to region 0, from c ′ ≧ 0, m ′ ≧ 0, y ′ ≧ 0 ′,
[0064]
[Formula 13]

[0065]
Is always 0, so terms containing these are 0. Therefore, (c ′, m ′, y ′, k)
[0066]
[Expression 14]

[0067]
It is only necessary to select and use the coefficient of the term consisting of three types of combinations. Similarly for each area,
[0068]
[Expression 15]

[0069]
If you choose the only person who does not become 0 out of, and of the 81 Section polynomial (in other words, chose polynomial 81 of q _i for each region).
[0070]
In this way, necessary coefficients are selected for each area from the coefficients (q _i ) in advance, and 81 coefficients (ai) are stored in the coefficient storage unit 3 for each area. Then, the region of the input signal (c, m, y, k) is discriminated by the region discriminating unit 2, and the coefficient ai of the region is read from the coefficient storage unit 3 and given to the polynomial arithmetic unit 4. The polynomial calculation unit 4 calculates the polynomial of the expression (1) using the read coefficient ai, and outputs L * a * b *.
[0071]
【The invention's effect】
As described above, according to the present invention, the input space is divided into a plurality of regions, and color conversion is performed by a polynomial using a different coefficient for each region. Therefore, a low-order polynomial is used for each region. Using the formula, it is possible to suppress unnecessary vibrations that tend to occur in higher-order polynomials. In addition, since the differential value of the polynomial is continuous at the boundary of the region, a color conversion method with smooth characteristics can be realized.
[0072]
According to the present invention, since a polynomial including a term of an arbitrary order is used, high-precision color conversion can be realized using the degree of freedom in each region.
[0073]
According to the present invention, an area is divided based on the magnitude relationship of each color signal, and for a printer close to a concentric model, the color conversion characteristics in each area are approximated by a polynomial, so that the ideal of the printer is obtained. The deviation from the characteristic can be absorbed and accurate approximation can be performed.
[0074]
According to the present invention, since the regions are divided based on the magnitude relationship between the black signal and the other color signals, the number of regions can be suppressed while maintaining the degree of freedom of the approximate polynomial.
[0075]
According to the present invention, a signal obtained by subtracting a black signal from a color signal other than black is generated, and the black subtraction signal and the black signal are color-converted by a polynomial. Therefore, the region division based on the magnitude relationship between black and other colors is performed. It is possible to easily set a coefficient at which the approximate polynomial is continuous at the boundary.
[Brief description of the drawings]
FIG. 1 shows a configuration of an embodiment of the present invention.
FIG. 2 shows a determination result of an input signal.
FIG. 3 shows each area determined in FIG.
FIG. 4 is a diagram illustrating a problem of polynomial approximation.
FIG. 5 is a diagram illustrating the relationship between cyan and magenta in a concentric circle model.
FIG. 6 is a diagram illustrating the relationship between cyan and black in a concentric circle model.
FIG. 7 is a diagram illustrating a relationship between c and k.
FIG. 8 is a diagram illustrating a relationship between ck and k.
[Explanation of symbols]
DESCRIPTION OF SYMBOLS 1 Subtractor 2 Area | region discrimination | determination part 3 Coefficient memory | storage part 4 Polynomial operation part

Claims

In a color signal conversion method for converting a first color image signal into a second color image signal having different attributes using a quadratic polynomial, when the first color image signal includes a black signal, A three-dimensional space composed of a signal obtained by subtracting a black signal from a color image signal is divided into a plurality of regions whose boundaries are planes where the value of each signal obtained by subtracting the black signal is 0, and each region is different. using the coefficients, and by the polynomial with coefficients primary differential value at the boundary region becomes continuously, color signal conversion and converting the first color image signal into a second color image signal Method.