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
Links
- 230000000873 masking effect Effects 0.000 claims description 40
- 238000011156 evaluation Methods 0.000 claims description 14
- 239000003086 colorant Substances 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 239000000976 ink Substances 0.000 description 16
- 238000000034 method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Landscapes
- Ink Jet (AREA)
- Facsimile Image Signal Circuits (AREA)
- Color Image Communication Systems (AREA)
Description
〔産業上の利用分野〕
本発明はカラーマスキングパラメータ決定装置
に関し、特にカラープリンタにおいてスキヤナか
ら得られた赤,緑,青の輝度信号を、シアン、マ
ゼンタ,イエローのインクの主濃度信号に変換す
るカラーマスキングのパラメータ決定装置に関す
るものである。
〔従来技術とその問題点〕
従来カラープリンタにおいて、色分解カラース
キヤナから得られた原稿の三原色(赤,緑,青)
の輝度信号(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)をカラースキヤナより測定し、
例えばシアンインクに関しては、
e2=N
〓K=1
(Ck―a11Drk―a12Dgk―a13Dbk)2 (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 = N 〓 K=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.
第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*
k,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*
k′,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/Y0)1/3−16
(Y/Y0>0.008856)
903.29 (Y/Y0)
(Y/Y0≦0.008856)
a*=500〔(X/X0)1/3−(Y/Y0)1/3〕
b*=200〔(Y/Y0)1/3−(Z/Z0)1/3〕
(19)
ここでX0,Y0,Z0は標準光源のX,Y,Zで
ある。
また、本実施例では、重みwkの決定に、対象
画素の全画素について最も近い仮想的カラーパツ
チを求めているが、通常の画像では近傍の画素
は、似た色を持つているので、n:1(nは任意
の数)に画素をまびいて同様の計算を行なつても
目的は達せられる。重みは式(12)の形で導入するの
が最も簡単であるが、応用によつて種々の変形が
可能である。例えば式(20)(21)のようなもの
である。
E2=N
〓K=1
wk 2{(L* k−L* k′)2+(u* k
−u* k′)2+(v* k−v* k′)2} (20)
E2=N
〓K=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 = N 〓 K=1 w k 2 {(L * k −L * k ′) 2 + (u * k −u * k ′) 2 + (v * k −v * k ′) 2 } (20 ) E 2 = N 〓 K=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.
第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)
及びインク三原色主濃度(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.
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)
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 |
-
1985
- 1985-06-07 JP JP60122684A patent/JPS61281768A/en active Granted
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 |