JPH0244436B2 - - Google Patents

Info

Publication number
JPH0244436B2
JPH0244436B2 JP60122684A JP12268485A JPH0244436B2 JP H0244436 B2 JPH0244436 B2 JP H0244436B2 JP 60122684 A JP60122684 A JP 60122684A JP 12268485 A JP12268485 A JP 12268485A JP H0244436 B2 JPH0244436 B2 JP H0244436B2
Authority
JP
Japan
Prior art keywords
color
masking
patch
weight
values
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
JP60122684A
Other languages
Japanese (ja)
Other versions
JPS61281768A (en
Inventor
Joji Tajima
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NEC Corp
Original Assignee
Nippon Electric Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nippon Electric Co Ltd filed Critical Nippon Electric Co Ltd
Priority to JP60122684A priority Critical patent/JPS61281768A/en
Publication of JPS61281768A publication Critical patent/JPS61281768A/en
Publication of JPH0244436B2 publication Critical patent/JPH0244436B2/ja
Granted legal-status Critical Current

Links

Landscapes

  • Ink Jet (AREA)
  • Facsimile Image Signal Circuits (AREA)
  • Color Image Communication Systems (AREA)

Description

【発明の詳細な説明】[Detailed description of the invention]

〔産業上の利用分野〕 本発明はカラーマスキングパラメータ決定装置
に関し、特にカラープリンタにおいてスキヤナか
ら得られた赤,緑,青の輝度信号を、シアン、マ
ゼンタ,イエローのインクの主濃度信号に変換す
るカラーマスキングのパラメータ決定装置に関す
るものである。 〔従来技術とその問題点〕 従来カラープリンタにおいて、色分解カラース
キヤナから得られた原稿の三原色(赤,緑,青)
の輝度信号(R,G,B)から、これを再現する
ための三原色(シアン,マゼンタ,イエロー)の
インクの主濃度信号(C,M,Y)を得る場合、
まず式(1)により三原色輝度濃度信号(R,G,
B)を三原色濃度信号(Dr,Dg,Db)に式(1)に
より変換し、 Dr=−logR Dg=−logG Db=−logB (1) 次に式(2)のように行列{aij}を用いて三原色濃
度信号(Dr,Dg,Db)をインクの主濃度信号の
近似値C′,M′,Y′に変換することが行われる。
この{aij)をカラーマスキングパラメータとよ
ぶ。 C′ M′ Y′=a11 a12 a13 a21 a22 a23 a31 a32 a33Dr Dg Db (2) この変換をできるだけ正確に行うことによつて
主濃度信号(C′,M′,Y′)のインクで印刷され
た色が元の(R,G,B)の色と等しい色に見え
るようにするために、従来では主濃度に関する最
小自乗法を実行してカラーマスキングパラメータ
{aTj}を求めることが行われていた。 具体的にはN個(Nは自然数)の既知のCk
Mk,Yk(k=1,……,N)を持つカラーパツ
チを印刷し、このカラーパツチの三原色濃度
(Drk,Dgk,Dbk)をカラースキヤナより測定し、
例えばシアンインクに関しては、 e2NK=1 (Ck―a11Drk―a12Dgk―a13Dbk2 (3) を最小にするa11,a12,a13を求める。これは連立
方程式(4)を解くことにより求められる。 マゼンタインク,イエローインクに関しても同
様にして{aij}の9つのカラーマスキングパラメ
ータを得ることができる。 しかしながら、上記の方法によつて得られたカ
ラーマスキングパラメータは、用いられるカラー
パツチの色のセツトによる影響を受け、実際に変
換される画像に最適なものとはならない。また、
最終的に人の眼で見て色差が最小であるような評
価にはL*u*V*系又はL*a*b*系のような均等色空
間が用いられるにも拘らず式(3)のe2の最小化は三
原色インクの主濃度の差に関する最小化であると
いう意味で最適なカラーマスキングパラメータの
決定法としては不十分なものであつた。 〔発明の目的〕 本発明の目的は、多数のカラーパツチを実際に
印刷することなく、仮想的なカラーパツチを用い
シミユレータにより人の眼で見て色差を最小とし
画像の色分布に対して再現色を最適に近似できる
カラーマスキングパラメータを決定するカラーマ
スキングパラメータ決定装置を提供することにあ
る。 〔発明の構成〕 本発明は、均等色空間中で三原色輝度値(R,
G,B)及びインク三原色主濃度(C,M,Y)
の両者で表現可能な色領域から一定間隔で複数の
仮想的カラーパツチを選択するカラーパツチ選択
手段と、該仮想的カラーパツチの均等色空間にお
ける座標値を記憶するカラーパツチ記憶手段と、
該座標値から各カラーパツチの三原色輝度値
(R,G,B)を計算する輝度値計算手段と、対
象とする画像の三原色輝度値(R,G,B)を均
等色空間における座標値に変換する色変換手段
と、該座標値と前記仮想的カラーパツチの座標値
とを比較して、評価値を計算するための各カラー
パツチの重みを計算する重み計算手段と、該重み
計算手段により得られた重みを各カラーパツチに
ついて記憶する重み記憶手段と、カラーマスキン
グパラメータを記憶するパラメータ記憶手段と、
前記輝度値計算手段により得られた三原色輝度値
と前記カラーマスキングパラメータからマスキン
グ計算によつてインク三原色主濃度の近似値
(C′,M′,Y′)を計算するマスキング計算手段
と、該主濃度の近似値(C′,M′,Y′)から均等
色空間中での座標値の近似値を計算し、前記カラ
ーパツチ記憶手段に記憶されている座標値との色
差を得、更に前記重み記憶手段に記憶されている
各カラーパツチに対する重みを用いて評価値を計
算する評価値計算手段と、該評価値を最小化する
ように前記カラーマスキングパラメータを更新す
る制御手段とから成り、収束計算にり対象とする
画像に最適なカラーマスキングパラメータを決定
することを特徴する。 〔発明の原理〕 本発明に係るカラーマスキングパラメータ決定
装置を実現する原理について説明する。カラー印
刷はシアン,マゼンタ,イエローの三原色インク
の網点印刷によつて行われるものとする。このと
き、各インクの実効面積率をそれぞれc,m,y
とすると、各インクの三原色主濃度信号C,M,
Yとは式(5)で関係づけられる。 C=−log(l−c) M=−log(l−m) Y=−log(l−y) (5) 網点印刷による再現色(X,Y,Z)は、上記
実効面積率c,m,yを用いるときClE―
1931XYZ系によつて式(6)のノイゲバウアー方程
式によつて予測される。
[Industrial Application Field] The present invention relates to a color masking parameter determination device, and in particular, to a color masking parameter determination device, which converts red, green, and blue luminance signals obtained from a scanner into main density signals of cyan, magenta, and yellow inks in a color printer. The present invention relates to a parameter determining device for color masking. [Prior art and its problems] In conventional color printers, the three primary colors (red, green, blue) of the original obtained from a color separation color scanner
When obtaining the main density signals (C, M, Y) of ink of the three primary colors (cyan, magenta, yellow) to reproduce this from the luminance signals (R, G, B) of
First, by formula (1), the three primary color luminance density signals (R, G,
B) into three primary color density signals (Dr, Dg, Db) using equation (1), Dr=-logR Dg=-logG Db=-logB (1) Next, the matrix {a ij } is used to convert the three primary color density signals (Dr, Dg, Db) into approximate values C', M', Y' of the main density signals of the ink.
This {a ij ) is called a color masking parameter. C′ M′ Y′=a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33 D r D g D b (2) By performing this conversion as accurately as possible, the main concentration signal (C In order to make the color printed with ink (', M', Y') appear to be the same as the original (R, G, B) color, conventionally, a least squares method regarding the principal density is performed. The color masking parameter {a Tj } was calculated. Specifically, N (N is a natural number) known C k ,
Print a color patch with M k , Y k (k=1,...,N), measure the three primary color densities (D rk , D gk , D bk ) of this color patch using a color scanner,
For example, for cyan ink, a 11 , a 12 , a 13 that minimizes e 2 = NK=1 (C k − a 11 D rk − a 12 D gk − a 13 D bk ) 2 (3) demand. This is obtained by solving simultaneous equations (4). Similarly, nine color masking parameters {a ij } can be obtained for magenta ink and yellow ink. However, the color masking parameters obtained by the above method are affected by the color set of the color patch used and are not optimal for the image actually converted. Also,
Although a uniform color space such as the L * u * V * system or the L * a * b * system is used to evaluate the minimum color difference when seen by the human eye, the formula (3 ) is an insufficient method for determining the optimal color masking parameter in the sense that the minimization of e 2 is a minimization of the differences in the main densities of the three primary color inks. [Object of the Invention] The object of the present invention is to use a simulator to minimize the color difference as seen by the human eye and to adjust the reproduced colors to the color distribution of an image using virtual color patches without actually printing a large number of color patches. An object of the present invention is to provide a color masking parameter determination device that determines a color masking parameter that can be optimally approximated. [Configuration of the Invention] The present invention provides three primary color luminance values (R,
G, B) and ink primary color main density (C, M, Y)
color patch selection means for selecting a plurality of virtual color patches at regular intervals from a color region that can be expressed by both; color patch storage means for storing coordinate values of the virtual color patches in a uniform color space;
A brightness value calculation means that calculates the three primary color brightness values (R, G, B) of each color patch from the coordinate values, and converts the three primary color brightness values (R, G, B) of the target image into coordinate values in a uniform color space. color conversion means for calculating the weight of each color patch for calculating an evaluation value by comparing the coordinate values with the coordinate values of the virtual color patch; weight storage means for storing weights for each color patch; parameter storage means for storing color masking parameters;
a masking calculation means for calculating approximate values (C', M', Y') of three primary color main densities by masking calculation from the three primary color brightness values obtained by the brightness value calculation means and the color masking parameters; An approximate coordinate value in a uniform color space is calculated from the approximate density value (C', M', Y'), a color difference with the coordinate value stored in the color patch storage means is obtained, and the weight is calculated. It consists of an evaluation value calculation means that calculates an evaluation value using the weight for each color patch stored in the storage means, and a control means that updates the color masking parameter so as to minimize the evaluation value, and performs convergence calculation. It is characterized by determining the optimal color masking parameters for the target image. [Principle of the Invention] The principle for realizing the color masking parameter determination device according to the present invention will be explained. Color printing is performed by dot printing using three primary color inks: cyan, magenta, and yellow. At this time, the effective area ratio of each ink is c, m, y, respectively.
Then, the three primary color main density signals C, M,
It is related to Y by equation (5). C=-log(l-c) M=-log(l-m) Y=-log(ly) (5) The reproduced color (X, Y, Z) by halftone printing is based on the above effective area ratio c , m, y, ClE-
It is predicted by the Neugebauer equation of equation (6) using the 1931XYZ system.

【表】 但し、(Xi,Yi,Zi)(i=1,……,8)は、
各インクのベタ印刷のすべての組み合わせに対し
て測定された(X,Y,Z)値であり、白に対し
てY=1に正規化されている。 iと印刷されるインクの対応は表1の通りであ
る。
[Table] However, (X i , Y i , Z i ) (i=1,...,8) is
(X, Y, Z) values measured for all combinations of solid prints for each ink, normalized to Y=1 for white. Table 1 shows the correspondence between i and the printed ink.

〔実施例〕〔Example〕

第1図は前述のカラーマスキングパラメータ決
定方法の原理に基づいて構成されたカラーマスキ
ングパラメータ決定装置のブロツク構成図であ
る。第1図において、まずカラーパツチ選択手段
1は仮想的カラーパツチセツトをL*u*v*空間中
で等間隔にサンプルし、(R,C,B)及び(c,
m,y)で表現できる色のみを選択し、これによ
つて選択された(L* k,u* k,v* k)をカラーパツ
チ記憶手段2に記憶する。上記カラーパツチ選択
手段1は、具体的にはマイクロコンピユータ等に
より第2図に示されるフローを実行することによ
つて実現される。 第2図に基づき上記(L* k,u* k,v* k)を選択
する方法を具体的に説明する。まずL*u*v*空間
を所要のΔdlに間隔で三次元的に走査し、最初カ
ラーパツチ候補色(L* P,u* P,v* P)を選ぶ(ス
テツプS1)。その走査範囲は、通常、 0≦L*≦100,−200≦u*≦300, −200≦v*≦200 (13) 程度をとれば十分である。次に上記の選ばれた各
候補についてR,G,B及びc,m,yで表現可
能か否かが検査される。まず各候補のL*,u*
v*は前記の式(8)の逆変換である式(14)によつ
てX,Y,Zに変換される(ステツプS2)。こ
の変換で求められた値を(XP,YP,ZP)とする。 上記(XP,YP,ZP)に対応する(R,G,B)
を式(7)により求め(ステツプS3)、求められた
R,G,Bが式(10)の条件を満たしているかを調べ
ることによりR,G,Bによる表現可能性が検査
できる(ステツプS4)。 次に上記式(14)により求められた(VP,YP
ZP)に対応する(cP,mP,yP)を解析的に求め
ることはできないので、第3図に示す如き逐次計
算で求める。まず(c,m,y)において表現可
能な立方体空間1を設定し、この立方体空間11
を等間隔d0から成る立方体領域に粗く分割し、各
領域を代表する(c,m,y)値について式(6)に
より(X,Y,Z)を求める。これらの(X,
Y,Z)の中で前記(XP,YP,ZP)に最も近い
(X,Y,Z)に有する領域12を更に細かく間
隔d1で分割する。かかる逐次計算をd1≦dnioにな
るまで繰返し、(XP,YP,ZP)に最も近い(X,
Y,Z)を有する領域の(c,m,y)を(cP
mP,yP)とする。ここでdnioは通常0.001程度で
ある。以上の逐次計算をステツプS5において行
う。このようにして得たcP,mP,yPのいずれか
が0または1の場合、すなわち立方体空間11の
面上にあるときは(XP,YP,ZP)はc,m,y
により表現不可能であつたとみなし、そうでない
場合には(XP,YP,ZP)はc,m,yにより表
現可能であるとみなす(ステツプS6)。上記の
各ステツプを繰返すことによりカラーパツチ候補
色(L* P,u* P,v* P)のうち(R,G,B)及び
(c,m,y)によつて表現可能な色のみが(L*
,u* k,v* k)としてカラーパツチ記憶手段2に
出力される(ステツプS7)。 第1図に戻り、色変換手段3は画像を走査し、
各画素の(R,G,B)を式(7)の逆変換及び式(8)
によつて(L*,u*,v*)に変換する。重み計算
手段4は、全画素の(L*,u*v*)についてカラ
ーパツチ記憶手段2に記憶されている各仮想的カ
ラーパツチ(L* k,u* k,v* k)のうち最も近いも
ののkに対する度数分布wkを計算し、重み記憶
手段5に記憶する。 輝度値計算手段6は、仮想的カラーパツチの
(L* k,u* k,v* k)を式(14),(7)を用いて三原色
輝度値(Rk,Gk,Bk)に変換する。 パラメータ記憶手段7には初期値として適当な
カラーマスキングパラメータ{aij}が格納されて
いる。マスキング計算手段8は輝度値計算手段6
から得られる三原色輝度値に対して式(15)によ
つてまず三原色濃度を求め、 Drk=−logRk Dgk=−logGk Dbk=−logBk (15) 次に式(16)によつてインク三原色主濃度の近
似値(Ck′,Mk′,Yk′)を求める。 Ck′ Mk′ Yk′=a11 a12 a13 a21 a22 a23 a31 a23 a33Drk Dgk Dbk (16) 評価値計算手段9は仮想的カラーパツチの色
(L* k,u* k,v* k)と、インク三原色主濃度の近似
値(Ck′,Mk′,Yk′)が示す色との色差の自乗
に重みwkを掛けて足し合わせた評価値を計算す
るもので、第4図示す演算ブロツクから成る。第
4図において演算ブロツク91は(Ck,Mk′,
Yk′)に対する各インクの実効面積率(Ck′,
mk′,yk′)式(17)に従つて求める。 ck′=1−10-Ck′ mk′=1−10-Mk′ Yk′=1−10-Yk′ (17) 演算ブロツク92は式6のノイゲバウアー方程
式を(ck′,mk′,yk′)について計算し(Xk′,
Yk′,Zk′)を得るものである。 α1k′=(1−ck′)(1−mk′)(1−yk′) α2k′=ck′(1−mk′)(1−yk′) α3k′=(1−ck′)mk′(1−yk′) α4k′=ck′mk′(1−yk′) α5k′=1−ck′)(1−mk′)yk′ α6k′=ck′(1−mk′)yk′ α7k′=(1−ck′)mk′yk′ α8k′=ck′mk′yk′ (18) 演算ブロツク93は式(8)によつて(Xk′,Yk′,
Zk′)から均等色空間における座標値(L* k′,u*
′,v* k′)を求める。但し、ここでは、式(8)の
X,Y,Z及びL*,u*,v*がXk′,Yk′,Zk′及
びL* k′,u* k′,v* k′に置きかられる。 演算ブロツク94は、上記計算で得られた
(L* k′,u* k′,v* k′)とカラーパツチ記憶手段2
からの仮想的カラーパツチの色(L* k,u* k,v* k
との均等色空間における距離を計算し、更に重み
記憶手段5から得られたこの仮想的カラーパツチ
の重みwkをかけ合わせて加算することによつて
評価値E2を得る。 制御手段10は、以上で得られた評価値を収束
計算によつて最小化するものであり、例えば非線
型数理計画法によりパラメータ記憶手段7中のカ
ラーマスキングパラメータ{aij}を更新する。制
御手段10は周知の電子計算機などで実現するこ
とができる。 以上のマスキング計算以降の処理は更新された
カラーマスキングパラメータ{aij}によつて繰り
返され、これ以上の改善が望めないことを制御手
段10が適当な基準により判断した時点で、最適
なカラーマスキングパラメータはパラメータ記憶
手段7中に格納されているものとして決定され
る。 上記のカラーパツチ記憶手段2、重み記憶手段
5、パラメータ記憶手段7は公知の半導体メモリ
等で実現できる。また色変換手段3、重み計算手
段4、輝度値計算手段6、マスキング計算手段
8、評価値計算手段9はすべて、公知の演算ユニ
ツトの組み合わせによつて容易に実現可能であ
る。 尚、上記の実施例では均等色空間としてL*u*
v*系を用いたが、L*a*b*系を用いても全く同様
に本発明を構成することができる。L*a*b*系の
XYZ系との関係は式(19)で与えられる。 L*=116(Y/Y01/3−16 (Y/Y0>0.008856) 903.29 (Y/Y0) (Y/Y0≦0.008856) a*=500〔(X/X01/3−(Y/Y01/3〕 b*=200〔(Y/Y01/3−(Z/Z01/3〕 (19) ここでX0,Y0,Z0は標準光源のX,Y,Zで
ある。 また、本実施例では、重みwkの決定に、対象
画素の全画素について最も近い仮想的カラーパツ
チを求めているが、通常の画像では近傍の画素
は、似た色を持つているので、n:1(nは任意
の数)に画素をまびいて同様の計算を行なつても
目的は達せられる。重みは式(12)の形で導入するの
が最も簡単であるが、応用によつて種々の変形が
可能である。例えば式(20)(21)のようなもの
である。 E2NK=1 wk 2{(L* k−L* k′)2+(u* k −u* k′)2+(v* k−v* k′)2} (20) E2NK=1 (logwk){(L* k−L* k′)2+ (u* k−u* k′)2+(v* k′)2} (21) 〔発明の効果〕 本発明によれば、輝度信号からインクの主濃度
信号に変換するカラーマスキングパラメータ決定
装置において実際にカラーパツチを印刷すること
なく、仮想的なカラーパツチを用いてカラーマス
キングパラメータを得るようにし、更に対象画像
の色分布に応じた重みを導入することによつてカ
ラーマスキングパラメータを得るようにしたた
め、人の眼で見て色差が最小となるような対象画
像の色分布に対して最適なカラーマスキングパラ
メータを決定することができる効果がある。
FIG. 1 is a block diagram of a color masking parameter determining apparatus constructed based on the principle of the color masking parameter determining method described above. In FIG. 1, the color patch selection means 1 first samples a virtual color patch set at equal intervals in L * u * v * space, and samples (R, C, B) and (c,
Only the colors that can be expressed as (L * k , u * k , v * k ) are selected and stored in the color patch storage means 2. Specifically, the color patch selection means 1 is realized by executing the flow shown in FIG. 2 using a microcomputer or the like. A method for selecting the above (L * k , u * k , v * k ) will be specifically explained based on FIG. First, the L * u * v * space is three-dimensionally scanned at intervals of a required Δdl, and first color patch candidate colors (L * P , u * P , v * P ) are selected (step S1). It is usually sufficient that the scanning range is approximately 0≦L * ≦100, −200≦u * ≦300, −200≦v * ≦200 (13). Next, it is checked whether each of the above selected candidates can be expressed using R, G, B, c, m, and y. First, each candidate's L * , u * ,
v * is converted into X, Y, and Z by equation (14), which is the inverse transformation of equation (8) above (step S2). Let the values obtained by this conversion be (X P , Y P , Z P ). (R, G, B) corresponding to the above (X P , Y P , Z P )
is obtained by equation (7) (step S3), and the expressibility of R, G, and B can be checked by checking whether the obtained R, G, and B satisfy the condition of equation (10) (step S4). ). Next, (V P , Y P ,
Since it is not possible to analytically obtain (c P , m P , y P ) corresponding to Z P ), it is obtained by sequential calculations as shown in FIG. First, we set a cubic space 1 that can be expressed in (c, m, y), and this cubic space 11
is roughly divided into cubic regions having equal intervals d 0 , and (X, Y, Z) is determined using equation (6) for the (c, m, y) values representative of each region. These (X,
The region 12 located at (X, Y, Z) closest to (X P , Y P , Z P ) among the regions Y, Z) is further divided into smaller intervals d 1 . This sequential calculation is repeated until d 1 ≦d nio , and the closest ( X , Y P , Z P ) is found.
(c, m, y) of the area with (Y, Z) as (c P ,
m P , y P ). Here, d nio is usually about 0.001. The above sequential calculations are performed in step S5. If any of c P , m P , and y P obtained in this way is 0 or 1, that is, when it is on the surface of the cubic space 11, (X P , Y P , Z P ) is c, m, y
If not, it is assumed that (X P , Y P , Z P ) can be expressed by c, m, and y (step S6). By repeating each of the above steps, only the colors that can be expressed by (R, G, B) and (c, m, y) among the color patch candidate colors (L * P , u * P , v * P ) are selected. (L *
k , u * k , v * k ) to the color patch storage means 2 (step S7). Returning to FIG. 1, the color conversion means 3 scans the image,
(R, G, B) of each pixel is inversely transformed by equation (7) and equation (8)
Convert to (L * , u * , v * ) by. The weight calculation means 4 calculates the closest one of each virtual color patch (L * k , u * k , v * k ) stored in the color patch storage means 2 for all pixels (L * , u * v * ) . A frequency distribution w k for k is calculated and stored in the weight storage means 5. The brightness value calculation means 6 converts (L * k , u * k , v * k ) of the virtual color patch into three primary color brightness values (R k , G k , B k ) using equations (14) and (7). Convert. The parameter storage means 7 stores appropriate color masking parameters {a ij } as initial values. The masking calculation means 8 is the brightness value calculation means 6
First, calculate the three primary color densities using equation (15) for the three primary color luminance values obtained from D rk = -logR k D gk = -logG k D bk = -logB k (15) Next, use equation (16). Therefore, approximate values (C k ′, M k ′, Y k ′) of the main densities of the three primary colors of ink are determined. C k ′ M k ′ Y k ′=a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 23 a 33 D rk D gk D bk (16) The evaluation value calculation means 9 calculates the color of the virtual color patch (L * k , u * k , v * k ) and the color indicated by the approximate values of the three primary color principal densities (C k ′, M k ′, Y k ′), multiplied by the weight w k and added together. It calculates the evaluated value and consists of the calculation blocks shown in Fig. 4. In FIG. 4, the calculation block 91 is (C k , M k ′,
Effective area ratio of each ink (C k ′, Y k ′)
m k ′, y k ′) is calculated according to equation (17). c k ′=1-10 -Ck ′ m k ′=1-10 -Mk ′ Y k ′=1-10 -Yk ′ (17) The calculation block 92 converts the Neugebauer equation of equation 6 into (c k ′, m k ′, y k ′) and (X k ′,
Y k ′, Z k ′). α 1k ′=(1−c k ′)(1−m k ′)(1−y k ′) α 2k ′=c k ′(1−m k ′)(1−y k ′) α 3k ′= (1−c k ′)m k ′(1−y k ′) α 4k ′=c k ′m k ′(1−y k ′) α 5k ′=1−c k ′)(1−m k ′ )y k ′ α 6k ′=c k ′ (1−m k ′)y k ′ α 7k ′=(1−c k ′)m k ′y k ′α 8k ′=c k ′m k ′y k ′ (18) The calculation block 93 calculates (X k ′, Y k ′,
Z k ′) to coordinate values (L * k ′, u *
k ′, v * k ′). However, here, X, Y, Z and L * , u * , v * in equation (8) are X k ′, Y k ′, Z k ′ and L * k ′, u * k ′, v * k ’. The calculation block 94 stores (L * k ', u * k ', v * k ') obtained in the above calculation and the color patch storage means 2.
The color of the virtual color patch from (L * k , u * k , v * k )
The evaluation value E 2 is obtained by calculating the distance in the uniform color space between the two and further multiplying by the weight w k of this virtual color patch obtained from the weight storage means 5 and adding the result. The control means 10 minimizes the evaluation value obtained above by convergence calculation, and updates the color masking parameter {a ij } in the parameter storage means 7 by, for example, nonlinear mathematical programming. The control means 10 can be realized by a well-known electronic computer or the like. The processing after the masking calculation described above is repeated using the updated color masking parameter {a ij }, and when the control means 10 determines based on an appropriate criterion that no further improvement can be expected, the optimum color masking is performed. The parameters are determined as stored in the parameter storage means 7. The color patch storage means 2, weight storage means 5, and parameter storage means 7 described above can be realized by a known semiconductor memory or the like. Further, the color conversion means 3, the weight calculation means 4, the brightness value calculation means 6, the masking calculation means 8, and the evaluation value calculation means 9 can all be easily realized by a combination of known calculation units. In the above embodiment, L * u * is used as a uniform color space.
Although the v * system was used, the present invention can be constructed in exactly the same way using the L * a * b * system. L * a * b * system
The relationship with the XYZ system is given by equation (19). L * = 116 (Y/Y 0 ) 1/3 -16 (Y/Y 0 >0.008856) 903.29 (Y/Y 0 ) (Y/Y 0 ≦0.008856) a * = 500 [(X/X 0 ) 1 /3 - (Y/Y 0 ) 1/3 ] b * = 200 [(Y/Y 0 ) 1/3 - (Z/Z 0 ) 1/3 ] (19) Here, X 0 , Y 0 , Z 0 is the standard light source's X, Y, and Z. Furthermore, in this embodiment, to determine the weight w k , the closest virtual color patch is found for all pixels of the target pixel, but in a normal image, neighboring pixels have similar colors, so n The purpose can be achieved even if the same calculation is performed by multiplying pixels by :1 (n is an arbitrary number). It is easiest to introduce the weights in the form of equation (12), but various modifications are possible depending on the application. For example, equations (20) and (21). E 2 = NK=1 w k 2 {(L * k −L * k ′) 2 + (u * k −u * k ′) 2 + (v * k −v * k ′) 2 } (20 ) E 2 = NK=1 (logw k ) {(L * k −L * k ′) 2 + (u * k −u * k ′) 2 + (v * k ′) 2 } (21) [ [Effects of the Invention] According to the present invention, in a color masking parameter determination device that converts a luminance signal into an ink main density signal, a color masking parameter is obtained using a virtual color patch without actually printing a color patch. In addition, the color masking parameters are obtained by introducing weights according to the color distribution of the target image, so that the color masking parameters are optimal for the color distribution of the target image so that the color difference as seen by the human eye is minimized. The effect is that color masking parameters can be determined.

【図面の簡単な説明】[Brief explanation of drawings]

第1図は本発明の一実施例を示すカラーマスキ
ングパラメータの決定装置のブロツク構成図、第
2図はカラーパツチ選択手段をマイクロコンピユ
ータで実現した場合の処理の流れを示すフローチ
ヤート、第3図はc,m,y立方体空間を分割す
る処理の説明図、第4図は評価値計算手段の詳細
なブロツク図である。 1……カラーパツチ選択手段、2……カラーパ
ツチ記憶手段、3……色変換手段、4……重み計
算手段、5……重み記憶手段、6……輝度値計算
手段、7……パラメータ記憶手段、8……マスキ
ング計算手段、9……評価値計算手段、10……
制御手段。
FIG. 1 is a block diagram of a color masking parameter determining device showing an embodiment of the present invention, FIG. 2 is a flowchart showing the process flow when the color patch selection means is realized by a microcomputer, and FIG. An explanatory diagram of the process of dividing the c, m, y cubic space, FIG. 4 is a detailed block diagram of the evaluation value calculation means. 1... Color patch selection means, 2... Color patch storage means, 3... Color conversion means, 4... Weight calculation means, 5... Weight storage means, 6... Luminance value calculation means, 7... Parameter storage means, 8... Masking calculation means, 9... Evaluation value calculation means, 10...
control means.

Claims (1)

【特許請求の範囲】[Claims] 1 均等色空間中で三原色輝度値(R,G,B)
及びインク三原色主濃度(C,M,Y)の両者で
表現可能な色領域から一定間隔で複数の仮想的カ
ラーパツチを選択するカラーパツチ選択手段と、
該仮想的カラーパツチの均等色空間における座標
値を記憶するカラーパツチ記憶手段と、該座標値
から各カラーパツチの三原色輝度値(R,G,
B)を計算する輝度値計算手段と、対象とする画
像の三原色輝度値(R,G,B)を均等色空間に
おける座標値に変換する色変換手段と、該座標値
と前記仮想的カラーパツチの座標値とを比較し
て、評価値を計算するための各カラーパツチの重
みを計算する重み計算手段と、該重み計算手段に
より得られた重みを各カラーパツチについて記憶
する重み記憶手段と、カラーマスキングパラメー
タを記憶するパラメータ記憶手段と、前記輝度値
計算手段により得られた三原色輝度値と前記カラ
ーマスキングパラメータからマスキング計算によ
つてインク三原色主濃度の近似値(C′,M′,Y′)
を計算するマスキング計算手段と、該主濃度の近
似値(C′,M′,Y′)から均等色空間中での座標
値の近似値を計算し、前記カラーパツチ記憶手段
に記憶されている座標値との色差を得、更に前記
重み記憶手段に記憶されている各カラーパツチに
対する重みを用いて評価値を計算する評価値計算
手段と、該評価値を最小化するように前記カラー
マスキングパラメータを更新する制御手段とから
成り、収束計算により対象とする画像に最適なカ
ラーマスキングパラメータを決定するカラーマス
キングパラメータ決定装置。
1 Luminance values of three primary colors (R, G, B) in uniform color space
and color patch selection means for selecting a plurality of virtual color patches at regular intervals from a color region that can be expressed by both the three primary color main densities (C, M, Y) of the ink;
a color patch storage means for storing the coordinate values of the virtual color patch in the uniform color space; and a color patch storage means for storing the coordinate values of the virtual color patch in the uniform color space;
B); a color conversion means that converts the three primary color brightness values (R, G, B) of the target image into coordinate values in a uniform color space; a weight calculation means for calculating the weight of each color patch for calculating an evaluation value by comparing the coordinate values with each other; a weight storage means for storing the weight obtained by the weight calculation means for each color patch; and a color masking parameter. and approximate values (C', M', Y') of the three primary color main densities of the ink by performing masking calculations from the three primary color luminance values obtained by the luminance value calculation means and the color masking parameters.
masking calculation means for calculating the coordinate values in the uniform color space from the approximate values (C', M', Y') of the main densities, and calculating the coordinate values stored in the color patch storage means; an evaluation value calculation means for calculating an evaluation value using a weight for each color patch stored in the weight storage means; and updating the color masking parameter so as to minimize the evaluation value. A color masking parameter determination device that determines the optimal color masking parameter for a target image by convergence calculation.
JP60122684A 1985-06-07 1985-06-07 Color masking parameter deciding device Granted JPS61281768A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP60122684A JPS61281768A (en) 1985-06-07 1985-06-07 Color masking parameter deciding device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP60122684A JPS61281768A (en) 1985-06-07 1985-06-07 Color masking parameter deciding device

Publications (2)

Publication Number Publication Date
JPS61281768A JPS61281768A (en) 1986-12-12
JPH0244436B2 true JPH0244436B2 (en) 1990-10-03

Family

ID=14842065

Family Applications (1)

Application Number Title Priority Date Filing Date
JP60122684A Granted JPS61281768A (en) 1985-06-07 1985-06-07 Color masking parameter deciding device

Country Status (1)

Country Link
JP (1) JPS61281768A (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2575134B2 (en) * 1987-04-10 1997-01-22 コニカ株式会社 Color estimation method by color patch
JPS63254887A (en) * 1987-04-10 1988-10-21 Konica Corp Color patch forming method
JP2529404B2 (en) * 1989-08-02 1996-08-28 松下電器産業株式会社 Masking coefficient determination device
EP0454275B1 (en) * 1990-04-26 1996-07-03 Bayer Corporation Color separation method
JP2543249B2 (en) * 1990-10-22 1996-10-16 松下電器産業株式会社 Color image forming method and apparatus and correction coefficient determining method
US5473439A (en) * 1992-10-23 1995-12-05 At&T Corp. Model-based halftoning of color images
JP2560972B2 (en) * 1993-05-31 1996-12-04 日本電気株式会社 Automatic color gamut conversion method, automatic color gamut conversion device and automatic color conversion device
JPH08116460A (en) * 1994-10-19 1996-05-07 Fuji Xerox Co Ltd Color masking parameter determining device

Also Published As

Publication number Publication date
JPS61281768A (en) 1986-12-12

Similar Documents

Publication Publication Date Title
US6435654B1 (en) Color calibration for digital halftoning
US7190487B2 (en) Color conversion with hue straightening using multiple look-up tables and interpolation
EP0495914B2 (en) Color digital halftoning with vector error diffusion
US5296947A (en) System for softproofing a color reproduction
EP0574905A2 (en) Method and apparatus for forming color images
US7889403B2 (en) Methods and apparatus for calibrating a digital color imaging device that uses multi-hue colorants
JPH1155536A (en) Color transformation table, color transformation, color signal generator, method for generating color transformation table, color transformation method and method for deciding color area
US7411701B2 (en) N-colorant gamut construction
US5781709A (en) Method of and apparatus for generating proof
JPH05183742A (en) Device for deciding color compensation parameter
JPH0244436B2 (en)
US7103216B2 (en) Color chart and color conversion definition creating method
JP2006340012A (en) Color adjustment method
JP4189188B2 (en) Color conversion definition correction method, color conversion definition correction apparatus, and color conversion definition correction program
JPH0518305B2 (en)
JP2554366B2 (en) Method of estimating color reproduction characteristics using color patches
JPH08116460A (en) Color masking parameter determining device
RU2190306C2 (en) Method and device for processing half-tone color picture using vector diffusion of chroma error
JPH10164352A (en) Image signal processing method
JP3430178B2 (en) Color correction method and apparatus for image processing system
JPH0517740B2 (en)
KR100503517B1 (en) Color sample formation method of printer
JPH0517739B2 (en)
Gonzalez et al. Color Management—how accurate need it be?
JPH0517742B2 (en)

Legal Events

Date Code Title Description
EXPY Cancellation because of completion of term