JP6155766B2 - Print reproduction color prediction method - Google Patents
Print reproduction color prediction method Download PDFInfo
- Publication number
- JP6155766B2 JP6155766B2 JP2013072945A JP2013072945A JP6155766B2 JP 6155766 B2 JP6155766 B2 JP 6155766B2 JP 2013072945 A JP2013072945 A JP 2013072945A JP 2013072945 A JP2013072945 A JP 2013072945A JP 6155766 B2 JP6155766 B2 JP 6155766B2
- Authority
- JP
- Japan
- Prior art keywords
- color
- spectral reflectance
- spectral
- effective area
- area ratio
- 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.)
- Active
Links
Images
Description
本発明は面積変調により色を再現する印刷物において再現される色すなわち再現色を予測する印刷再現色予測方法に関するものである。 The present invention relates to a print reproduction color prediction method for predicting a color to be reproduced in a printed material that reproduces a color by area modulation, that is, a reproduction color .
従来、印刷の再現色を予測する方法として、例えばシアン（Ｃ）、マゼンタ（Ｍ）、イエロー（Ｙ）、ブラック（Ｋ）の各単色のデバイス制御値と実際にその掛け合わせを出力した際の色を測色した測色値との関係からルックアップテーブル（以下、ＬＵＴと称する）を作成する方法が用いられている。例えば、非特許文献１に、ＬＵＴを用いた再現色予測方法が記載されている。
Conventionally, as a method for predicting the reproduction color of printing, for example, when a device control value of each single color of cyan (C), magenta (M), yellow (Y), and black (K) is actually output by multiplying them. A method of creating a look-up table (hereinafter referred to as LUT) from the relationship with colorimetric values obtained by measuring colors is used. For example, Non-Patent
この方法においては、例えば各原色を適当な間隔で区切ったすべての掛け合わせを出力し、その結果を測色することで、図２に示すようなＬＵＴを作成する。図２において各格子点は適当な間隔に区切った入力値に対応し、格子点上には実際に出力して測定した測色値が格納されている。格子点が存在しない掛け合わせ、すなわち実際に測定を行っていない掛け合わせに対しては、周辺の格子点が持つ値から、補間などの手法を用いて、任意の掛け合わせに対する再現色の測色値を求めることができる。 In this method, for example, all multiplications obtained by dividing each primary color at an appropriate interval are output, and the result is measured to create an LUT as shown in FIG. In FIG. 2, each grid point corresponds to an input value divided at an appropriate interval, and a colorimetric value actually output and measured is stored on the grid point. For multiplications where there are no grid points, that is, multiplications that are not actually measured, color measurement of the reproduced color for any multiplication is performed using a method such as interpolation from the values of surrounding grid points. The value can be determined.
もう１つの方法として、ノイゲバウア方程式を利用する方法が知られている。例えば、特許文献１に、ノイゲバウア方程式を用いた色予測の方法が記載されている。この方法では、各原色のベタおよびベタ同士のすべての掛け合わせと用紙の測色値をあらかじめ測定し、以下の式で再現色を予測する。
As another method, a method using the Neugebauer equation is known. For example,
ＸＹＺ＝ｐｃＸＹＺｃｉ＋ｐｍＸＹＺｍｉ＋ｐｙＸＹＺｙｉ＋ｐｋＸＹＺｋｉ＋ｐｃｍＸＹＺｃｍｉ＋ｐｍｙＸＹＺｍｙｉ＋ｐｃｙＸＹＺｃｙｉ＋ｐｃｋＸＹＺｃｋｉ＋ｐｍｋＸＹＺｍｋｉ＋ｐｙｋＸＹＺｙｋｉ＋ｐｃｍｙＸＹＺｃｍｙｉ＋ｐｃｍｋＸＹＺｃｍｋｉ＋ｐｍｙｋＸＹＺｍｙｋｉ＋ｐｃｙｋＸＹＺｃｙｋｉ＋ｐｃｙｍｋＸＹＺｃｙｍｋｉ＋ｐｐＸＹＺｐ
ｐｃ＝ｃ（１−ｍ）（１−ｙ）（１−ｋ）
ｐｍ＝（１−ｃ）ｍ（１−ｙ）（１−ｋ）
ｐｙ＝（１−ｃ）（１−ｍ）ｙ（１−ｋ）
ｐｋ＝（１−ｃ）（１−ｍ）（１−ｙ）ｋ
ｐｃｍ＝ｃｍ（１−ｙ）（１−ｋ）
ｐｍｙ＝（１−ｃ）ｍ（１−ｂ）（１−ｋ）
ｐｃｙ（λ）＝ｃ（１−ｍ）ｙ（１−ｋ）
ｐｃｋ（λ）＝（１−ｃ）ｍｙ（１−ｋ）
ｐｍｋ（λ）＝（１−ｃ）ｍ（１−ｙ）ｋ
ｐｙｋ（λ）＝（１−ｃ）（１−ｍ）ｙｋ
ｐｃｍｙ（λ）＝ｃｍｙ（１−ｋ）
ｐｃｍｋ（λ）＝ｃｍ（１−ｙ）ｋ
ｐｍｙｋ（λ）＝（１−ｃ）ｍｙｋ
ｐｃｙｋ（λ）＝ｃ（１−ｍ）ｙｋ
ｐｃｙｍｋ（λ）＝ｃｍｙｋ
ｐｐ（λ）＝（１−ｃ）（１−ｍ）（１−ｙ）（１−ｋ）
（ただし、ｃ，ｍ，ｙ，ｋは各々原色の実効面積率、
Ｐｊは４色重ね刷り時の各原色ベタ，二次色ベタ，三次色ベタ，四次色ベタ，用紙の実効
面積率、
ＸＹＺｊｉは各原色ベタ，二次色ベタ，三次色ベタ，四次色ベタ，用紙の測色値ベクトルを示す）
以下、この式を（式１）と称する。
XYZ = p c XYZ ci + p m XYZ mi + p y XYZ yi + p k XYZ ki + p cm XYZ cmi + p my XYZ myi + p cy XYZ cyi + p ck XYZ cki + p mk XYZ mki + p yk XYZ yki + p cmy XYZ cmyi + p cmk XYZ cmki + P myk XYZ myki + p cyk XYZ cyki + p cymk XYZ cymki + pp XYZ p
p c = c (1-m) (1-y) (1-k)
p m = (1-c) m (1-y) (1-k)
p y = (1-c) (1-m) y (1-k)
p k = (1-c) (1-m) (1-y) k
p cm = cm (1-y) (1-k)
p my = (1-c) m (1-b) (1-k)
p cy (λ) = c (1-m) y (1-k)
p ck (λ) = (1-c) my (1-k)
p mk (λ) = (1-c) m (1-y) k
p yk (λ) = (1-c) (1-m) yk
p cmy (λ) = cmy (1-k)
p cmk (λ) = cm (1-y) k
p myk (λ) = (1-c) myk
p cyk (λ) = c (1−m) yk
p cymk (λ) = cmyk
p p (λ) = (1-c) (1-m) (1-y) (1-k)
(Where c, m, y, and k are the primary area effective area ratios,
P j each primary color solid at the time I 4 color superimposition printing, secondary color solid, tertiary color solid, fourth color solid, the effective area ratio of the paper,
XYZ ji indicates each primary color solid, secondary color solid, tertiary color solid, quaternary color solid, and a colorimetric value vector of the paper)
Hereinafter, this equation is referred to as (Equation 1).
また、この（式１）において各原色の測色値の代わりに分光反射率ベクトルを用いて再現色の分光反射率を推定し、観察光源の分光分布と等色関数を用いた計算により、再現色の測色値を求める方法をとることもある。このとき、（式１）は以下に示す式となり、分光ノイゲバウア方程式と呼ばれる。 Also, in this (Equation 1), the spectral reflectance of the reproduced color is estimated using the spectral reflectance vector instead of the colorimetric value of each primary color, and reproduced by calculation using the spectral distribution and color matching function of the observation light source. A method for obtaining a colorimetric value of a color may be used. At this time, (Expression 1) becomes an expression shown below and is called a spectral Neugebauer equation.
Ｒ＾（λ）＝ｐｃ（λ）Ｒｃｉ（λ）＋ｐｍ（λ）Ｒｍｉ（λ）＋ｐｙ（λ）Ｒｙｉ（λ）＋ｐｋ（λ）Ｒｋｉ（λ）＋ｐｃｍ（λ）Ｒｃｍｉ（λ）＋ｐｍｙ（λ）Ｒｍｙｉ（λ）＋ｐｃｙ（λ）Ｒｃｙｉ（λ）＋ｐｃｋ（λ）Ｒｃｋｉ（λ）＋ｐｍｋ（λ）Ｒｍｋｉ（λ）＋ｐｙｋ（λ）Ｒｙｋｉ（λ）＋ｐｃｍｙ（λ）Ｒｃｍｙｉ（λ）＋ｐｃｍｋ（λ）Ｒｃｍｋｉ（λ）＋ｐｍｙｋ（λ）Ｒｍｙｋｉ（λ）＋ｐｃｙｋ（λ）Ｒｃｙｋｉ（λ）＋ｐｃｙｍｋ（λ）Ｒｃｙｍｋｉ（λ）＋ｐｐ（λ）Ｒｐ（λ）（ただし、Ｒｉ（λ）は各原色ベタ，二次色ベタ，三次色ベタ，四次色ベタ，用紙の分光反射率ベクトルを示す）
以下、この式を（式２）と称する。
R ^ (λ) = p c (λ) R ci (λ) + p m (λ) R mi (λ) + p y (λ) R yi (λ) + p k (λ) R ki (λ) + p cm (λ ) R cmi (λ) + p my (λ) R myi (λ) + p cy (λ) R cyi (λ) + p ck (λ) R cki (λ) + p mk (λ) R mki (λ) + p yk (λ ) R yki (λ) + p cmy (λ) R cmyi (λ) + p cmk (λ) R cmki (λ) + p myk (λ) R myki (λ) + p cyk (λ) R cyki (λ) + p cymk (λ ) R cymki (λ) + p p (λ) R p (λ) (where R i (λ) is the primary solid, secondary color solid, tertiary color solid, quaternary color solid, and the spectral reflectance vector of the paper. Show)
Hereinafter, this equation is referred to as (Equation 2).
分光反射率と測色値の関係は以下の式で与えられる。 The relationship between the spectral reflectance and the colorimetric value is given by the following equation.
Ｘ＝ｋ∫Ｒ（λ）Ｅ（λ）ｘ￣（λ）ｄλ
Ｙ＝ｋ∫Ｒ（λ）Ｅ（λ）ｙ￣（λ）ｄλ
Ｚ＝ｋ∫Ｒ（λ）Ｅ（λ）ｚ￣（λ）ｄλ
ｋ＝１００／∫Ｅ（λ）ｙ￣（λ）ｄλ
（ただし、Ｒ（λ）は分光反射率、Ｅ（λ）は観察光源の分光エネルギー分布、ｘ￣（λ）、ｙ￣（λ）、ｚ￣（λ）は各々の等色関数を示す）
以下、この式を（式３）と称する。
X = k∫R (λ) E (λ) x￣ (λ) dλ
Y = k∫R (λ) E (λ) y￣ (λ) dλ
Z = k∫R (λ) E (λ) z￣ (λ) dλ
k = 100 / ∫E (λ) y￣ (λ) dλ
(Where R (λ) is the spectral reflectance, E (λ) is the spectral energy distribution of the observation light source, and x￣ (λ), y￣ (λ), and z￣ (λ) are the respective color matching functions)
Hereinafter, this equation is referred to as (Equation 3).
しかし、上述したＬＵＴを用いる方法では、測定するサンプルの数を増加させることにより、補間による誤差が減少し、比較的高い精度での再現色予測が可能となる反面、測定数の増加により、作業時間や手間が増大するという問題があった。また、数多くの掛け合わせを出力するには非常に大きな面積が必要となり、印刷物面内のばらつきが大きい印刷機においては、その面内ばらつきの誤差がＬＵＴに含まれ精度が低下するという問題があった。 However, in the method using the LUT described above, by increasing the number of samples to be measured, errors due to interpolation can be reduced and reproduction color prediction can be performed with relatively high accuracy. There was a problem of increasing time and labor. In addition, a very large area is required to output a large number of multiplications, and in a printing press with large variations in the printed surface, errors in the variations in the surface are included in the LUT, resulting in a decrease in accuracy. It was.
一方、ノイゲバウア方程式を用いる方法では、測色が必要な掛け合わせは各色ベタとベタ同士の掛け合わせとなるため非常に少なくなるものの、実際の印刷では光学的ドットゲインの影響により、この方程式が成立することはほとんどなく、再現色予測の精度に問題があった。 On the other hand, in the method using the Neugebauer equation, the multiplication that requires color measurement is very little because it is a multiplication between each color and solid, but this equation is established in actual printing due to the effect of optical dot gain. There was little to do, and there was a problem in the accuracy of reproduction color prediction.
光学的ドッドゲインの影響により生じるノイゲバウア方程式と実測の差異を補正するために、特許文献２のような方法が提案されている。
In order to correct the difference between the Neugebauer equation and the actual measurement caused by the influence of the optical Dodd gain, a method as disclosed in
この方法では、１つの原色ベタに対し、光学的ドットゲインにより生じるベタとは異なる仮想的な原色が１つ存在するものと仮定し、（式１）を修正することにより、より高い精度で再現色を予測することを可能としている。 In this method, it is assumed that there is one virtual primary color different from the solid color generated by optical dot gain for one primary color solid, and it is reproduced with higher accuracy by correcting (Equation 1). It is possible to predict the color.
しかしながら、上述した方法は、商業印刷で主に用いられている１７５線程度のＡＭ印
刷を対象としており、網点形状が大きく異なるＦＭ印刷や網点の大きさが小さい高細線のＡＭ印刷では、光学的ドットゲインの影響が補償しきれず、再現色予測の精度が低下するという問題があった。
However, the above-described method is intended for AM printing of about 175 lines, which is mainly used in commercial printing, and in FM printing with greatly different halftone dot shapes and high thin line AM printing with small halftone dot sizes, There was a problem that the effect of optical dot gain could not be fully compensated and the accuracy of reproduction color prediction was reduced.
本発明は、上述の従来の技術の問題点を解決すべくなされたものであって、印刷の面内ばらつきの影響を受けない比較的狭い面積に印刷可能な少数の掛け合わせの分光測定値から、高い精度で任意の掛け合わせのデバイス制御値に対する分光反射率を推定し、印刷物の再現色を容易に推定する方法を提供することを目的としている。 The present invention has been made to solve the above-mentioned problems of the prior art, and is based on a small number of multiplying spectroscopic measurement values that can be printed in a relatively small area that is not affected by in-plane variation of printing. Another object of the present invention is to provide a method for estimating a spectral reflectance for an arbitrarily multiplied device control value with high accuracy and easily estimating a reproduced color of a printed matter .
前記課題を解決するために、第１の発明は、面積変調で色を再現する印刷物の任意の掛け合わせに対する分光反射率を推定し、再現色を予測する方法において、
用紙、各色ベタ、ベタの掛け合わせについて分光反射率を測定するステップと、
一次色の任意のデバイス制御値に対する分光反射率を測定するステップと、
一次色の任意のデバイス制御値に対する分光反射率を、用紙の分光反射率、ベタの分光反射率、ドットゲインを考慮し、インク層を光が１回のみ通過する部分に対応する第１の仮想原色と、インク層を一部のみ光が透過する部分に対応する第２の仮想原色からなる２つの仮想的原色の分光反射率の４つの分光反射率をそれぞれの実効面積率で重み付けした
和として表現し、測定した一次色の実測値の存在するデバイス制御値に対する分光反射率との誤差が最小となるように、ベタおよび仮想的原色の実効面積率を求めるステップと、
前段のステップを複数の実測値の存在するデバイス制御値について繰り返し行い、一次色のデバイス制御値とベタ、第１の仮想原色、第２の仮想原色の実効面積率の関係を近似式で表し、デバイス制御値と実効面積率の関係をもとめるステップと、
一次色の任意のデバイス制御値に対する分光反射率を、
上記表現とは異なる表現として用紙の分光反射率とベタの分光反射率の線形和で表現し、
ベタと用紙それぞれの係数となる実効面積率が波長の関数である分光実効面積率であるとして、この分光実効面積率を求めるステップと、
任意の設定値の掛け合わせに対する分光反射率を、
上記の分光実効面積率を分光ノイゲバウア方程式に代入することにより求めるステップと、
観察光源、等色関数を用いて、前段のステップで求めた分光反射率から再現色の測色値を計算するステップと、
からなることを特徴とした印刷再現色予測方法である。
In order to solve the above-mentioned problem, the first invention is a method for estimating a spectral reflectance for an arbitrary multiplication of a printed material that reproduces a color by area modulation, and predicting a reproduced color.
Measuring spectral reflectance for the paper, each color solid, solid multiplication,
Measuring the spectral reflectance for any device control value of the primary color;
The spectral reflectance for an arbitrary device control value of the primary color is considered in consideration of the spectral reflectance of the paper, the solid spectral reflectance, and the dot gain, and the first virtual corresponding to the portion where the light passes only once through the ink layer As a sum of four spectral reflectances of the primary colors and the spectral reflectances of the two virtual primary colors corresponding to the part of the ink layer through which light passes only partially , weighted by the respective effective area ratios. Expressing and measuring the effective area ratio of the solid and virtual primary colors so as to minimize the error from the spectral reflectance with respect to the device control value in which the measured value of the measured primary color is present ;
The previous step is repeated for a device control value having a plurality of actual measurement values, and the relationship between the primary color device control value and the solid, the first virtual primary color, and the effective area ratio of the second virtual primary color is expressed by an approximate expression, A step for determining the relationship between the device control value and the effective area ratio;
Spectral reflectance for any device control value of the primary color,
Expressed as a linear sum of the spectral reflectance of the paper and the spectral reflectance of the solid as an expression different from the above expression ,
Assuming that the effective area ratio that is a coefficient of each of the solid and the paper is a spectral effective area ratio that is a function of the wavelength, obtaining the spectral effective area ratio ;
Spectral reflectance for multiplying any set value
Obtaining the spectral effective area ratio by substituting it into the spectral Neugebauer equation;
A step of calculating a colorimetric value of the reproduced color from the spectral reflectance obtained in the previous step using an observation light source and a color matching function;
A print reproduction color prediction method characterized by comprising:
また、前記課題を解決するために第２の発明は、第１の発明において、一次色の中間調の分光反射率を、Ｒ＾（ｓ，λ）＝Ｒｉ（λ）ａ１（ｓ）＋Ｒo１（λ）ａ２（ｓ）＋Ｒo２（λ）ａ３（ｓ）＋Ｒｐ（λ）（１−ａ１（ｓ）−ａ２（ｓ）−ａ３（ｓ））（但し、ａ１（ｓ）は単色ベタの実効面積率、ａ２（ｓ）は第１の仮想的原色の実効面積率、ａ３（ｓ）は第２の仮想的原色の実効面積率、Ｒｉ（λ）は単色ベタの分光反射率、Ｒｐ（λ）は用紙の分光反射率、Ｒo１（λ）は第１の仮想原色の分光反射率、Ｒo２（λ）は第２の仮想原色の分光反射率、Ｒ＾（ｓ，λ）はデバイス制御値ｓに対応する単色の推定分光反射率を示す）で示し、
実効面積率をｂ（ｓ，λ）＝（Ｒ＾（ｓ，λ）−Ｒｐ（λ））／（Ｒｉ（λ）−Ｒｐ（λ））としたとき、
分光ノイゲバウア方程式として以下の式を用いることを特徴とした印刷再現色予測方法である。
In order to solve the above-mentioned problem, the second invention is the first invention, wherein the spectral reflectance of the halftone of the primary color is expressed as R ^ (s, λ) = R i (λ) a 1 (s). + R o1 (λ) a 2 (s) + R o2 (λ) a 3 (s) + R p (λ) (1-a 1 (s) -a 2 (s) -a 3 (s)) (where a 1 (s) is the effective area ratio of the single solid color, a 2 (s) is the effective area ratio of the first virtual primary color, a 3 (s) is the effective area ratio of the second virtual primary color, and R i (λ ) Is the spectral reflectance of the solid solid color, R p (λ) is the spectral reflectance of the paper, R o1 (λ) is the spectral reflectance of the first virtual primary color, and R o2 (λ) is the spectral spectrum of the second virtual primary color. The reflectance, R ^ (s, λ) indicates the monochromatic estimated spectral reflectance corresponding to the device control value s),
When the effective area ratio is b (s, λ) = (R ^ (s, λ) −R p (λ)) / (R i (λ) −R p (λ)),
A print reproduction color prediction method using the following equation as a spectral Neugebauer equation.
Ｒ＾（λ）＝ｐｃ（λ）Ｒｃｉ（λ）＋ｐｍ（λ）Ｒｍｉ（λ）＋ｐｙ（λ）Ｒｙｉ（λ）＋ｐｋ（λ）Ｒｋｉ（λ）＋ｐｃｍ（λ）Ｒｃｍｉ（λ）＋ｐｍｙ（λ）Ｒｍｙｉ（λ）＋ｐｃｙ（λ）Ｒｃｙｉ（λ）＋ｐｃｋ（λ）Ｒｃｋｉ（λ）＋ｐｍｋ（λ）Ｒｍｋｉ（λ）＋ｐｙｋ（λ）Ｒｙｋｉ（λ）＋ｐｃｍｙ（λ）Ｒｃｍｙｉ（λ）＋ｐｃｍｋ（λ）Ｒｃｍｋｉ（λ）＋ｐｍｙｋ（λ）Ｒｍｙｋｉ（λ）＋ｐｃｙｋ（λ）Ｒｃｙｋｉ（λ）＋ｐｃｙｍｋ（λ）Ｒｃｙｍｋｉ（λ）＋ｐｐ（λ）Ｒｐ（λ）
ｐｃ（λ）＝ｂｃ（ｓｃ，λ）（１−ｂｍ（ｓｍ，λ））（１−ｂｙ（ｓｙ，λ））（１−ｂｋ（ｓｋ，λ））
ｐｍ（λ）＝（１−ｂｃ（ｓｃ，λ））ｂｍ（ｓｍ，λ）（１−ｂｙ（ｓｙ，λ））（１−ｂｋ（ｓｋ，λ））
ｐｙ（λ）＝（１−ｂｃ（ｓｃ，λ））（１−ｂｍ（ｓｍ，λ））ｂｙ（ｓｙ，λ）（１−ｂｋ（ｓｋ，λ））
ｐｋ（λ）＝（１−ｂｃ（ｓｃ，λ））（１−ｂｍ（ｓｍ，λ））（１−ｂｙ（ｓｙ，λ））ｂｋ（ｓｋ，λ）
ｐｃｍ（λ）＝ｂｃ（ｓｃ，λ）ｂｍ（ｓｍ，λ）（１−ｂｙ（ｓｙ，λ））（１−ｂｋ（ｓｋ，λ））
ｐｍｙ（λ）＝（１−ｂｃ（ｓｃ，λ））ｂｍ（ｓｍ，λ）ｂｙ（ｓｙ，λ）（１−ｂｋ（ｓｋ，λ））
ｐｃｙ（λ）＝ｂｃ（ｓｃ，λ）（１−ｂｍ（ｓｍ，λ））ｂｙ（ｓｙ，λ）（１−ｂｋ（ｓｋ，λ））
ｐｃｋ（λ）＝ｂｃ（ｓｃ，λ）（１−ｂｍ（ｓｍ，λ））（１−ｂｙ（ｓｙ，λ））ｂｋ（ｓｋ，λ）
ｐｍｋ（λ）＝（１−ｂｃ（ｓｃ，λ））ｂｍ（ｓｍ，λ）（１−ｂｙ（ｓｙ，λ））ｂｋ（ｓｋ，λ）
ｐｙｋ（λ）＝（１−ｂｃ（ｓｃ，λ））（１−ｂｍ（ｓｍ，λ））ｂｙ（ｓｙ，λ）ｂｋ（ｓｋ，λ）
ｐｃｍｙ（λ）＝ｂｃ（ｓｃ，λ）ｂｍ（ｓｍ，λ）ｂｙ（ｓｙ，λ）（１−ｂｋ（ｓｋ，λ））
ｐｃｍｋ（λ）＝ｂｃ（ｓｃ，λ）ｂｍ（ｓｍ，λ）（１−ｂｙ（ｓｙ，λ））ｂｋ（ｓｋ，λ）
ｐｍｙｋ（λ）＝（１−ｂｃ（ｓｃ，λ））ｂｍ（ｓｍ，λ）ｂｙ（ｓｙ，λ）ｂｋ（ｓｋ，λ）
ｐｃｙｋ（λ）＝ｂｃ（ｓｃ，λ）（１−ｂｍ（ｓｍ，λ））ｂｙ（ｓｙ，λ）ｂｋ（ｓｋ，λ）
ｐｃｙｍｋ（λ）＝ｂｃ（ｓｃ，λ）ｂｍ（ｓｍ，λ）ｂｙ（ｓｙ，λ）ｂｋ（ｓｋ，λ）
ｐｐ（λ）＝（１−ｂｃ（ｓｃ，λ））（１−ｂｍ（ｓｍ，λ））（１−ｂｙ（ｓｙ，λ））（１−ｂｋ（ｓｋ，λ））
R ^ (λ) = p c (λ) R ci (λ) + p m (λ) R mi (λ) + p y (λ) R yi (λ) + p k (λ) R ki (λ) + p cm (λ ) R cmi (λ) + p my (λ) R myi (λ) + p cy (λ) R cyi (λ) + p ck (λ) R cki (λ) + p mk (λ) R mki (λ) + p yk (λ ) R yki (λ) + p cmy (λ) R cmyi (λ) + p cmk (λ) R cmki (λ) + p myk (λ) R myki (λ) + p cyk (λ) R cyki (λ) + p cymk (λ ) R cymki (λ) + p p (λ) R p (λ)
p c (λ) = b c (s c , λ) (1-b m (s m , λ)) (1-b y (s y , λ)) (1-b k (s k , λ))
p m (λ) = (1−b c (s c , λ)) b m (s m , λ) (1−b y (s y , λ)) (1−b k (s k , λ))
p y (λ) = (1−b c (s c , λ)) (1−b m (s m , λ)) b y (s y , λ) (1−b k (s k , λ))
p k (λ) = (1−b c (s c , λ)) (1−b m (s m , λ)) (1−b y (s y , λ)) b k (s k , λ)
p cm (λ) = b c (s c , λ) b m (s m , λ) (1-b y (s y , λ)) (1-b k (s k , λ))
p my (λ) = (1−b c (s c , λ)) b m (s m , λ) b y (s y , λ) (1−b k (s k , λ))
p cy (λ) = b c (s c , λ) (1−b m (s m , λ)) b y (s y , λ) (1−b k (s k , λ))
p ck (λ) = b c (s c , λ) (1−b m (s m , λ)) (1−b y (s y , λ)) b k (s k , λ)
p mk (λ) = (1−b c (s c , λ)) b m (s m , λ) (1−b y (s y , λ)) b k (s k , λ)
p yk (λ) = (1−b c (s c , λ)) (1−b m (s m , λ)) b y (s y , λ) b k (s k , λ)
p cmy (λ) = b c (s c , λ) b m (s m , λ) b y (s y , λ) (1−b k (s k , λ))
p cmk (λ) = b c (s c , λ) b m (s m , λ) (1−b y (s y , λ)) b k (s k , λ)
p myk (λ) = (1−b c (s c , λ)) b m (s m , λ) b y (s y , λ) b k (s k , λ)
p cyk (λ) = b c (s c , λ) (1−b m (s m , λ)) b y (s y , λ) b k (s k , λ)
p cymk (λ) = b c (s c , λ) b m (s m , λ) b y (s y , λ) b k (s k , λ)
p p (λ) = (1-b c (s c , λ)) (1-b m (s m , λ)) (1-b y (s y , λ)) (1-b k (s k , Λ) )
以上、説明したように、本発明によれば、少数の一次色ステップとベタの掛け合わせおよび用紙の測定値から任意の掛け合わせのデバイス制御値に対する分光反射率を精度良く推定することができ、従って任意の掛け合わせのデバイス制御値に対する印刷物の再現色を容易に推定することが可能となる。また、必要とする測定が少ないことから、少ない面積で測定のためのパッチを出力することができ、印刷の面内ばらつき等の影響を受けないという特徴もある。 As described above, according to the present invention, it is possible to accurately estimate the spectral reflectance with respect to a device control value of an arbitrary multiplication from a small number of primary color steps and a solid multiplication and a measured value of the paper, Accordingly, it is possible to easily estimate the reproduced color of the printed material for an arbitrarily multiplied device control value. In addition, since less measurement is required, it is possible to output a patch for measurement in a small area, and there is a feature that it is not affected by in-plane variation of printing.
また、本発明に記載の設定値算出方法によれば、上記印刷再現色予測方法を使用することで、印刷の面内ばらつきによる影響を受けにくいといった特徴を引き継ぎつつ、任意の所望色の測色値、もしくは分光反射率を再現するためのデバイス制御値を容易に算出することができる。 Further, according to the set value calculation method described in the present invention, by using the print reproduction color prediction method, it is possible to perform color measurement of an arbitrary desired color while taking over the feature that it is not easily affected by in-plane variation of printing. The device control value for reproducing the value or the spectral reflectance can be easily calculated.
以下、本発明の一実施形態による印刷再現色予測方法を、図面を参照して説明する。 Hereinafter, a print reproduction color prediction method according to an embodiment of the present invention will be described with reference to the drawings.
図１は本発明の一実施形態による印刷再現色予測方法の処理の流れを表すフローチャートである。このフローチャートにそって各ステップでの処理の詳細を説明する。なお、この実施例ではＣＭＹＫの４原色の面積変調により色を再現する印刷を対象とする。 FIG. 1 is a flowchart showing a processing flow of a print reproduction color prediction method according to an embodiment of the present invention. Details of processing in each step will be described with reference to this flowchart. Note that this embodiment is intended for printing that reproduces colors by area modulation of the four primary colors of CMYK.
（ステップ１）用紙、単色ベタ、ベタ掛け合わせの測定（フローチャートの１０１）
用紙、単色ベタ（Ｃ，Ｍ，Ｙ，Ｋ）、ベタの掛け合わせ（ＣＭ，ＭＹ，ＣＹ，ＣＫ，ＭＫ，ＹＫ，ＣＭＹ，ＣＭＫ，ＭＹＫ，ＣＹＫ，ＣＭＹＫ）の計１６色について、分光反射率を測定する。
(Step 1) Measurement of paper, single color solid, solid multiplication (101 in the flowchart)
Spectral reflectance for a total of 16 colors of paper, single color solid (C, M, Y, K), and solid multiplication (CM, MY, CY, CK, MK, YK, CMY, CMK, MYK, CYK, CMYK) Measure.
（ステップ２） 単色ステップ測定（フローチャートの１０２）
ＣＭＹＫの各単色のデバイス制御値を任意のステップで変化させたパッチの分光反射率を測定する。本実施例では各単色のステップは１０％刻みとする。また、面積率の小さい用紙付近の色については、後述するように他の部分と特性が異なるため、その特性をより正確に把握するために、１０％まで２％刻みの一次色を同様に測定する。
(Step 2) Monochromatic step measurement (102 in the flowchart)
The spectral reflectance of a patch obtained by changing the device control value of each single color of CMYK in an arbitrary step is measured. In the present embodiment, each single color step is set in increments of 10%. In addition, the color near the paper with a small area ratio is different from the other parts as described later. Therefore, in order to grasp the characteristic more accurately, the primary color is measured in increments of 2% up to 10%. To do.
（ステップ３） 設定値と実効面積率の関係を求める（フローチャートの１０３）
ステップ１およびステップ２の各ステップで得られた分光反射率データから、任意の一次色のデバイス制御値に対する分光反射率を求めるため、任意のデバイス制御値ｓに対応する分光反射率を次式のように表現する。
(Step 3) Obtaining the relationship between the set value and the effective area ratio (103 in the flowchart)
In order to obtain the spectral reflectance for the device control value of an arbitrary primary color from the spectral reflectance data obtained in each step of
Ｒ＾（ｓ，λ）＝Ｒｉ（λ）ａ１（ｓ）＋Ｒo１（λ）ａ２（ｓ）＋Ｒo２（λ）ａ３（
ｓ）＋Ｒｐ（λ）（１−ａ１（ｓ）−ａ２（ｓ）−ａ３（ｓ））とし（但し、ａ１（ｓ）は単色ベタの実効面積率、ａ２（ｓ）は第１の仮想的原色の実効面積率、ａ３（ｓ）は第２の仮想的原色の実効面積率、Ｒｉ（λ）は単色ベタの分光反射率、Ｒｐ（λ）は用紙の分光反射率、Ｒo１（λ）は第１の仮想原色の分光反射率、Ｒo２（λ）は第２の仮想原色の分光反射率、Ｒ＾（ｓ，λ）はデバイス制御値ｓに対応する単色の推定分光反射率を示す）
以下、この式を（式４）と称する。
R ^ (s, λ) = R i (λ) a 1 (s) + R o1 (λ) a 2 (s) + R o2 (λ) a 3 (
s) + R p (λ) (1-a 1 (s) -a 2 (s) -a 3 (s)) (where a 1 (s) is the effective area ratio of a solid monochrome image and a 2 (s) Is the effective area ratio of the first virtual primary color, a 3 (s) is the effective area ratio of the second virtual primary color, R i (λ) is the spectral reflectance of the monochromatic solid, and R p (λ) is the paper Spectral reflectance, R o1 (λ) corresponds to the spectral reflectance of the first virtual primary color, R o2 (λ) corresponds to the spectral reflectance of the second virtual primary color, and R ^ (s, λ) corresponds to the device control value s. Shows the estimated spectral reflectance of a single color)
Hereinafter, this equation is referred to as (Equation 4).
上記（式４）のＲo１（λ）、Ｒo２（λ）について説明する。このパラメータは、面積変調により色を再現する印刷で発生する光学的ドットゲインと呼ばれる現象を考慮するために加えられている。図３は光学的ドットゲインが発生する原理を示したものである。図３に示したように、通常インクが乗った部分では、インク部分から入射した光が再びインク部分から射出してくるが、インク部と用紙部の境界部分では、インク部から入射した光が用紙部から射出したり、用紙部から入射した光がインク部から射出することにより、観測されるインクの面積率は実際のインク面積率と異なることとなる。 R o1 (λ) and R o2 (λ) in the above (Formula 4) will be described. This parameter is added in order to consider a phenomenon called optical dot gain that occurs in printing that reproduces colors by area modulation. FIG. 3 shows the principle of generation of optical dot gain. As shown in FIG. 3, the light incident from the ink portion is emitted from the ink portion again at the portion where the normal ink is placed, but the light incident from the ink portion is emitted at the boundary portion between the ink portion and the paper portion. When the light emitted from the paper portion or light incident from the paper portion is emitted from the ink portion, the observed ink area ratio is different from the actual ink area ratio.
本実施例では図４に示すように光学的ドットゲインによる発色をベタ部とは異なる分光反射率を持った仮想的な原色が存在するものと想定し、その分光反射率をＲo１（λ）Ｒo２（λ）としている。 In this embodiment, as shown in FIG. 4, it is assumed that there is a virtual primary color having a spectral reflectance different from that of the solid portion for color development by optical dot gain, and the spectral reflectance is R o1 (λ). R o2 (λ).
次に、Ｒo１（λ）Ｒo２（λ）を既知のデータで表現するために以下の式を考える。 Next, in order to express R o1 (λ) R o2 (λ) with known data, the following equation is considered.
Ｋｉ（λ）＝（−１／２）ｌｎ（Ｒｉ（λ）／Ｒｐ（λ））
以下、この式を（式５）と称する。
K i (λ) = (− 1/2) ln (R i (λ) / R p (λ))
Hereinafter, this equation is referred to as (Equation 5).
この式で表されるＫｉ（λ）はインクの吸光率と考えることができる。先に示した仮想的な原色についても同様な式が定義できる。 K i (λ) represented by this equation can be considered as the absorbance of the ink. Similar expressions can be defined for the virtual primary colors shown above.
Ｋo１（λ）＝（−１／２）ｌｎ（Ｒo１（λ）／Ｒｐ（λ））
以下、この式を（式６）と称する。
K o1 (λ) = (− 1/2) ln (R o1 (λ) / R p (λ))
Hereinafter, this equation is referred to as (Equation 6).
Ｋo２（λ）＝（−１／２）ｌｎ（Ｒo２（λ）／Ｒｐ（λ））
以下、この式を（式７）と称する。
K o2 (λ) = (− 1/2) ln (R o2 (λ) / R p (λ))
Hereinafter, this equation is referred to as (Equation 7).
光学的ドットゲインにより発生する仮想的な原色について、ベタ部ではインク層を２回通過するのに対し（図３のＡ）、光学的ドットゲインが発生する部分ではインク層を１回のみ通過（図３のＢ）もしくは、インク層を一部のみ通過（図３のＣ）することから、上記吸光率が半分および、その半分の１／４になるとすると、Ｋo１（λ）＝−１／２Ｒｉ（λ）
（以下、この式を（式８）と称する）、および、Ｋo２（λ）＝−１／４Ｋｉ（λ）（以下、この式を（式９）と称する）、とすることができ、この（式８）および（式９）に（式５）、（式６）、（式７）を代入することによりＲo１（λ）およびＲo２（λ）は既知のデータから次の式で表される。
For the virtual primary color generated by the optical dot gain, the solid portion passes through the ink layer twice (A in FIG. 3), whereas the portion where the optical dot gain occurs passes only once through the ink layer ( B) in FIG. 3 or only a part of the ink layer (C in FIG. 3) passes, so if the absorbance is half and 1/4 of the half, K o1 (λ) = − 1 / 2R i (λ)
(Hereinafter, this equation is referred to as (Equation 8)) and K o2 (λ) = − 1 / 4K i (λ) (hereinafter, this equation is referred to as (Equation 9)), By substituting (Equation 5), (Equation 6), and (Equation 7) into (Equation 8) and (Equation 9), R o1 (λ) and R o2 (λ) can be obtained from the known data by the following equation: expressed.
Ｒo１（λ）＝（Ｒｉ（λ）／Ｒｐ（λ））１／２Ｒｐ（λ）
以下、この式を（式１０）と称する。
R o1 (λ) = (R i (λ) / R p (λ)) 1/2 R p (λ)
Hereinafter, this equation is referred to as (Equation 10).
Ｒo２（λ）＝（Ｒｉ（λ）／Ｒｐ（λ））１／４Ｒｐ（λ）
以下、この式を（式１１）と称する。
R o2 (λ) = (R i (λ) / R p (λ)) 1/4 R p (λ)
Hereinafter, this equation is referred to as (Equation 11).
この式を（式４）に代入することにより、任意のデバイス制御値に対する分光反射率の推定値は次式で表される。 By substituting this equation into (Equation 4), an estimated value of the spectral reflectance with respect to an arbitrary device control value is expressed by the following equation.
Ｒ＾（ｓ，λ）＝Ｒｉ（λ）ａ１（ｓ）＋（Ｒｉ（λ）／Ｒｐ（λ））１／２Ｒｐ（λ）ａ２（ｓ）＋Ｒｐ（λ）ａ２（ｓ）＋（Ｒｉ（λ）／Ｒｐ（λ））１／４Ｒｐ（λ）ａ３（ｓ）＋Ｒｐ（λ）（１−ａ１（ｓ）−ａ２（ｓ）−ａ３（ｓ））
以下、この式を（式１２）と称する。
R ^ (s, λ) = R i (λ) a 1 (s) + (R i (λ) / R p (λ)) 1/2 R p (λ) a 2 (s) + R p (λ) a 2 (s) + (R i (λ) / R p (λ)) 1/4 R p (λ) a 3 (s) + R p (λ) (1-a 1 (s) −a 2 (s -A 3 (s))
Hereinafter, this equation is referred to as (Equation 12).
ここで、（ステップ２）で実際に測定した一次色ステップの実測値Ｒ（ｓ，λ）と推定値Ｒ＾（ｓ，λ）の誤差が最小となるように（式１３）のｒｍｓが最小となるようなａ１，ａ２，ａ３を最適化手法を用いて計算することにより、実測値の存在するデバイス制御値に対する実効面積率ａ１，ａ２，ａ３が求められ、図５に示すような傾向を示す。 Here, the rms in (Equation 13) is minimized so that the error between the measured value R (s, λ) and the estimated value R ^ (s, λ) of the primary color step actually measured in (Step 2) is minimized. become such a 1, a 2, by calculating with a 3 optimization techniques, the effective area ratio a 1, a 2, a 3 are determined for the device control values existing in the measured value, FIG. 5 The tendency as shown in is shown.
ｒｍｓ＝Σ√（（Ｒ（ｓ，λ）−Ｒ＾（ｓ，λ））２／ｎ）
以下、この式を（式１３）と称する。
rms = Σ√ ((R (s, λ) −R ^ (s, λ)) 2 / n)
Hereinafter, this equation is referred to as (Equation 13).
ａ１，ａ２，ａ３をそれぞれデバイス制御値ｓの関数として近似することにより、任意のデバイス制御値ｓに対する実効網点面積率を計算することができるようになり、（式１２）から任意のデバイス制御値に対する一次色の分光反射率を推定することができるようになる。なお、本実施例ではａ１，ａ２，ａ３ともにｓの２次関数で近似するものとする。ただし、図５に示すように、デバイス制御値の小さい部分において、ａ１は０となる傾向がある。これは図６に示すように、デバイス制御値の小さい部分では描画される網点の大きさが小さく、インク層を２回通過する光が非常に少ないためであると考えられる。この特性を把握し、近似式に反映することが肝要であるため、本実施例では（ステップ２）の単色ステップ測定のステップにおいて、１０％までの２％刻みのステップを測定した。これらの結果から、ａ１，ａ２，ａ３を２次関数で近似すると、近似値に負となる部分が生じる。しかし、ａ１，ａ２，ａ３は網点面積率であり、負となることはありえない。これにより、ａ１，ａ２，ａ３の近似式は非負の条件付きの２次式となる。 By approximating a 1 , a 2 , and a 3 as functions of the device control value s, the effective halftone dot area ratio for an arbitrary device control value s can be calculated. The spectral reflectance of the primary color with respect to the device control value can be estimated. In the present embodiment, all of a 1 , a 2 , and a 3 are approximated by a quadratic function of s. However, as shown in FIG. 5, a 1 tends to be 0 in a portion where the device control value is small. As shown in FIG. 6, this is considered to be because the halftone dot to be drawn is small in the portion where the device control value is small, and very little light passes through the ink layer twice. Since it is important to grasp this characteristic and reflect it in the approximate expression, in this example, steps of 2% increments up to 10% were measured in the step of monochromatic step measurement of (Step 2). From these results, when a 1 , a 2 , and a 3 are approximated by a quadratic function, a negative portion is generated in the approximate value. However, a 1 , a 2 , and a 3 are halftone dot area ratios and cannot be negative. Thereby, the approximate expression of a 1 , a 2 , a 3 becomes a quadratic expression with a non-negative condition.
これらの計算を各原色について実行することにより、それぞれの原色でのデバイス制御値と実効面積率の関係を表す関数ａ１ｊ（ｓ），ａ２ｊ（ｓ），ａ３ｊ（ｓ）（ｊ＝Ｃ，Ｍ，Ｙ，Ｋ）を求めることが可能となり、これらを用いて各原色の任意のデバイス制御値に対する分光反射率を推定することができる。 By executing these calculations for each primary color, functions a 1j (s), a 2j (s), a 3j (s) (j = C) representing the relationship between the device control value and the effective area ratio for each primary color. , M, Y, K), which can be used to estimate the spectral reflectance of each primary color with respect to an arbitrary device control value.
（ステップ４）分光実効面積率の計算（フローチャートの１０４）
先に示した光学的ドットゲインによる影響を、ここでは実効面積率が波長の関数として変化するとして以下の式で表すこととする。
(Step 4) Calculation of spectral effective area ratio (104 in the flowchart)
The influence of the optical dot gain shown above is expressed by the following equation assuming that the effective area ratio changes as a function of wavelength.
Ｒ＾（ｓ，λ）＝Ｒｉ（λ）ｂ（ｓ，λ）＋Ｒｐ（λ）（１−ｂ（ｓ，λ））
以下、この式を（式１４）と称する。
R ^ (s, λ) = R i (λ) b (s, λ) + R p (λ) (1-b (s, λ))
Hereinafter, this equation is referred to as (Equation 14).
この（式１４）においてｂ（ｓ，λ）は、各原色の分光実効面積率と考えることができる。この式をｂ（ｓ，λ）について解くと次式を得る。 In this (Expression 14), b (s, λ) can be considered as the spectral effective area ratio of each primary color. When this equation is solved for b (s, λ), the following equation is obtained.
ｂ（ｓ，λ）＝（Ｒ＾（ｓ，λ）−Ｒｐ（λ））／（Ｒｉ（λ）−Ｒｐ（λ））
以下、この式を（式１５）と称する。
b (s, λ) = (R ^ (s, λ) −R p (λ)) / (R i (λ) −R p (λ))
Hereinafter, this equation is referred to as (Equation 15).
ステップ３では任意のデバイス制御値ｓに対して求めることができることから、（式１５）を用いて任意のデバイス制御値ｓに対する分光実効面積率を計算することができる。 Since it is possible to obtain an arbitrary device control value s in step 3, the spectral effective area ratio for the arbitrary device control value s can be calculated using (Equation 15).
（ステップ５）ノイゲバウア方程式による分光反射率の推定（フローチャートの１０５）
（ステップ４）で求めた分光実効面積率を分光ノイゲバウア方程式に代入することにより、以下の式を得る。
(Step 5) Estimation of spectral reflectance by Neugebauer equation (105 in the flowchart)
By substituting the spectral effective area ratio obtained in (Step 4) into the spectral Neugebauer equation, the following formula is obtained.
Ｒ＾（λ）＝ｐｃ（λ）Ｒｃｉ（λ）＋ｐｍ（λ）Ｒｍｉ（λ）＋ｐｙ（λ）Ｒｙｉ（λ）＋ｐｋ（λ）Ｒｋｉ（λ）＋ｐｃｍ（λ）Ｒｃｍｉ（λ）＋ｐｍｙ（λ）Ｒｍｙｉ（λ）＋ｐｃｙ（λ）Ｒｃｙｉ（λ）＋ｐｃｋ（λ）Ｒｃｋｉ（λ）＋ｐｍｋ（λ）Ｒｍｋｉ（λ）＋ｐｙｋ（λ）Ｒｙｋｉ（λ）＋ｐｃｍｙ（λ）Ｒｃｍｙｉ（λ）＋ｐｃｍｋ（λ）Ｒｃｍｋｉ（λ）＋ｐｍｙｋ（λ）Ｒｍｙｋｉ（λ）＋ｐｃｙｋ（λ）Ｒｃｙｋｉ（λ）＋ｐｃｙｍｋ（λ）Ｒｃｙｍｋｉ（λ）＋ｐｐ（λ）Ｒｐ（λ）
ｐｃ（λ）＝ｂｃ（ｓｃ，λ）（１−ｂｍ（ｓｍ，λ））（１−ｂｙ（ｓｙ，λ））（１−ｂｋ（ｓｋ，λ））
ｐｍ（λ）＝（１−ｂｃ（ｓｃ，λ））ｂｍ（ｓｍ，λ）（１−ｂｙ（ｓｙ，λ））（１−ｂｋ（ｓｋ，λ））
ｐｙ（λ）＝（１−ｂｃ（ｓｃ，λ））（１−ｂｍ（ｓｍ，λ））ｂｙ（ｓｙ，λ）（１−ｂｋ（ｓｋ，λ））
ｐｋ（λ）＝（１−ｂｃ（ｓｃ，λ））（１−ｂｍ（ｓｍ，λ））（１−ｂｙ（ｓｙ，λ））ｂｋ（ｓｋ，λ）
ｐｃｍ（λ）＝ｂｃ（ｓｃ，λ）ｂｍ（ｓｍ，λ）（１−ｂｙ（ｓｙ，λ））（１−ｂｋ（ｓｋ，λ））
ｐｍｙ（λ）＝（１−ｂｃ（ｓｃ，λ））ｂｍ（ｓｍ，λ）ｂｙ（ｓｙ，λ）（１−ｂｋ（ｓｋ，λ））
ｐｃｙ（λ）＝ｂｃ（ｓｃ，λ）（１−ｂｍ（ｓｍ，λ））ｂｙ（ｓｙ，λ）（１−ｂｋ（ｓｋ，λ））
ｐｃｋ（λ）＝ｂｃ（ｓｃ，λ）（１−ｂｍ（ｓｍ，λ））（１−ｂｙ（ｓｙ，λ））ｂｋ（ｓｋ，λ）
ｐｍｋ（λ）＝（１−ｂｃ（ｓｃ，λ））ｂｍ（ｓｍ，λ）（１−ｂｙ（ｓｙ，λ））ｂｋ（ｓｋ，λ）
ｐｙｋ（λ）＝（１−ｂｃ（ｓｃ，λ））（１−ｂｍ（ｓｍ，λ））ｂｙ（ｓｙ，λ）ｂｋ（ｓｋ，λ）
ｐｃｍｙ（λ）＝ｂｃ（ｓｃ，λ）ｂｍ（ｓｍ，λ）ｂｙ（ｓｙ，λ）（１−ｂｋ（ｓｋ，λ））
ｐｃｍｋ（λ）＝ｂｃ（ｓｃ，λ）ｂｍ（ｓｍ，λ）（１−ｂｙ（ｓｙ，λ））ｂｋ（ｓｋ，λ）
ｐｍｙｋ（λ）＝（１−ｂｃ（ｓｃ，λ））ｂｍ（ｓｍ，λ）ｂｙ（ｓｙ，λ）ｂｋ（ｓｋ，λ）
ｐｃｙｋ（λ）＝ｂｃ（ｓｃ，λ）（１−ｂｍ（ｓｍ，λ））ｂｙ（ｓｙ，λ）ｂｋ（ｓｋ，λ）
ｐｃｙｍｋ（λ）＝ｂｃ（ｓｃ，λ）ｂｍ（ｓｍ，λ）ｂｙ（ｓｙ，λ）ｂｋ（ｓｋ，λ）
ｐｐ（λ）＝（１−ｂｃ（ｓｃ，λ））（１−ｂｍ（ｓｍ，λ））（１−ｂｙ（ｓｙ，λ））（１−ｂｋ（ｓｋ，λ））（ただし、ｂｊ（ｓ，λ）は（１５）式で求めた各原色の実効面積率である）
以下、この式を（式１６）と称する。
R ^ (λ) = p c (λ) R ci (λ) + p m (λ) R mi (λ) + p y (λ) R yi (λ) + p k (λ) R ki (λ) + p cm (λ ) R cmi (λ) + p my (λ) R myi (λ) + p cy (λ) R cyi (λ) + p ck (λ) R cki (λ) + p mk (λ) R mki (λ) + p yk (λ ) R yki (λ) + p cmy (λ) R cmyi (λ) + p cmk (λ) R cmki (λ) + p myk (λ) R myki (λ) + p cyk (λ) R cyki (λ) + p cymk (λ ) R cymki (λ) + p p (λ) R p (λ)
p c (λ) = b c (s c , λ) (1-b m (s m , λ)) (1-b y (s y , λ)) (1-b k (s k , λ))
p m (λ) = (1−b c (s c , λ)) b m (s m , λ) (1−b y (s y , λ)) (1−b k (s k , λ))
p y (λ) = (1−b c (s c , λ)) (1−b m (s m , λ)) b y (s y , λ) (1−b k (s k , λ))
p k (λ) = (1−b c (s c , λ)) (1−b m (s m , λ)) (1−b y (s y , λ)) b k (s k , λ)
p cm (λ) = b c (s c , λ) b m (s m , λ) (1-b y (s y , λ)) (1-b k (s k , λ))
p my (λ) = (1−b c (s c , λ)) b m (s m , λ) b y (s y , λ) (1−b k (s k , λ))
p cy (λ) = b c (s c , λ) (1−b m (s m , λ)) b y (s y , λ) (1−b k (s k , λ))
p ck (λ) = b c (s c , λ) (1−b m (s m , λ)) (1−b y (s y , λ)) b k (s k , λ)
p mk (λ) = (1−b c (s c , λ)) b m (s m , λ) (1−b y (s y , λ)) b k (s k , λ)
p yk (λ) = (1−b c (s c , λ)) (1−b m (s m , λ)) b y (s y , λ) b k (s k , λ)
p cmy (λ) = b c (s c , λ) b m (s m , λ) b y (s y , λ) (1−b k (s k , λ))
p cmk (λ) = b c (s c , λ) b m (s m , λ) (1−b y (s y , λ)) b k (s k , λ)
p myk (λ) = (1−b c (s c , λ)) b m (s m , λ) b y (s y , λ) b k (s k , λ)
p cyk (λ) = b c (s c , λ) (1−b m (s m , λ)) b y (s y , λ) b k (s k , λ)
p cymk (λ) = b c (s c , λ) b m (s m , λ) b y (s y , λ) b k (s k , λ)
p p (λ) = (1-b c (s c , λ)) (1-b m (s m , λ)) (1-b y (s y , λ)) (1-b k (s k , Λ)) (where b j (s, λ) is the effective area ratio of each primary color obtained by equation (15)).
Hereinafter, this equation is referred to as (Equation 16).
この式を用いることにより、任意のデバイス制御値の掛け合わせに対する分光反射率を推定することができる。 By using this equation, it is possible to estimate the spectral reflectance for multiplication of an arbitrary device control value.
（ステップ６） 色値の計算（フローチャートの１０６）
前段のステップで求めた分光反射率から、観察光源の分光分布と、等色関数をもちいて以下の式により測色値Ｘ，Ｙ，Ｚを求めることができる。
(Step 6) Color value calculation (106 in the flowchart)
From the spectral reflectance obtained in the previous step, the colorimetric values X, Y, and Z can be obtained by the following formula using the spectral distribution of the observation light source and the color matching function.
Ｘ＝ｋ∫Ｒ（λ）Ｅ（λ）ｘ￣（λ）ｄλ
Ｙ＝ｋ∫Ｒ（λ）Ｅ（λ）ｙ￣（λ）ｄλ
Ｚ＝ｋ∫Ｒ（λ）Ｅ（λ）ｚ￣（λ）ｄλ
ｋ＝１００／∫Ｅ（λ）ｙ￣（λ）ｄλ
（ただし、Ｅ（λ）は観察光源の分光エネルギー分布、ｘ￣（λ）、ｙ￣（λ）、ｚ￣（λ）は各々の等色関数を示す）
以下、この式を（式１７）と称する。
X = k∫R (λ) E (λ) x￣ (λ) dλ
Y = k∫R (λ) E (λ) y￣ (λ) dλ
Z = k∫R (λ) E (λ) z￣ (λ) dλ
k = 100 / ∫E (λ) y￣ (λ) dλ
(Where E (λ) is the spectral energy distribution of the observation light source, and x￣ (λ), y￣ (λ), and z￣ (λ) indicate the respective color matching functions)
Hereinafter, this equation is referred to as (Equation 17).
以上、本発明の再現色予測方法の一実施例について説明したが、以下のような変形例も考えることができる。 As mentioned above, although one Example of the reproduction color prediction method of this invention was described, the following modifications can also be considered.
本実施例では、ベタ、ベタの掛け合わせ、用紙の測定を実際に行っているが、これらがデータとして取得可能な場合には、測定を行わずこれらのデータを使用しても良い。 In this embodiment, solid, solid multiplication, and paper measurement are actually performed. However, when these can be acquired as data, these data may be used without performing measurement.
本実施例ではデバイス制御値と実効網点面積率の関係を、デバイス制御値の二次関数で近似しているが、より次数の高い多項式や、その他の関数で近似してもよい。 In this embodiment, the relationship between the device control value and the effective halftone dot area ratio is approximated by a quadratic function of the device control value, but may be approximated by a higher order polynomial or other functions.
本実施例では、網点面積率の小さい部分での特性を把握するために０％から１０％のデバイス制御値に対して２％刻みの一次色を測定しているが、この範囲をより広く（例えば０％から２０％）とったり、測定の間隔を小さく（例えば１％間隔）するなどしても良い。 In this embodiment, the primary color is measured in increments of 2% with respect to the device control value of 0% to 10% in order to grasp the characteristics in the portion where the halftone dot area ratio is small. (For example, 0% to 20%) or the measurement interval may be reduced (for example, 1% interval).
また、測色値としてＸＹＺ値を計算しているが、さらにそれらの値から、ＣＩＥＬＡＢやＣＩＥＬＵＶなどの異なる測色値を計算するようにしてもよい。 Further, although XYZ values are calculated as colorimetric values, different colorimetric values such as CIELAB and CIELV may be calculated from these values.
本実施例では、ＣＭＹＫを用いて説明を行っているが、この限りではなく、特色を含む４色印刷や４色以上の印刷にこの手法を用いてもよい。 In the present embodiment, the description is given using CMYK, but the present invention is not limited to this, and this method may be used for four-color printing including special colors or printing of four or more colors.
本実施例では、仮想的原色を２つとしているが、この限りではなく、３つ以上の仮想的原色を用いてもよい。その場合、ベタの吸光度と第ｎの仮想的原色の吸収度の関係は、１／２ｎとする。 In this embodiment, two virtual primary colors are used. However, the present invention is not limited to this, and three or more virtual primary colors may be used. In this case, the relationship between the absorbance of the solid and the absorbance of the nth virtual primary color is 1 / 2n .
１０１〜１０６ ステップ 101-106 steps
Claims (2)
用紙、各色ベタ、ベタの掛け合わせについて分光反射率を測定するステップと、
一次色の任意のデバイス制御値に対する分光反射率を測定するステップと、
一次色の任意のデバイス制御値に対する分光反射率を、用紙の分光反射率、ベタの分光反射率、ドットゲインを考慮し、インク層を光が１回のみ通過する部分に対応する第１の仮想原色と、インク層を一部のみ光が透過する部分に対応する第２の仮想原色からなる２つの仮想的原色の分光反射率の４つの分光反射率をそれぞれの実効面積率で重み付けした和として表現し、測定した一次色の実測値の存在するデバイス制御値に対する分光反射率との誤差が最小となるように、ベタおよび仮想的原色の実効面積率を求めるステップと、
前段のステップを複数の実測値の存在するデバイス制御値について繰り返し行い、一次色のデバイス制御値とベタ、第１の仮想原色、第２の仮想原色の実効面積率の関係を近似式で表し、デバイス制御値と実効面積率の関係をもとめるステップと、
一次色の任意のデバイス制御値に対する分光反射率を、
上記表現とは異なる表現として用紙の分光反射率とベタの分光反射率の線形和で表現し、
ベタと用紙それぞれの係数となる実効面積率が波長の関数である分光実効面積率であるとして、この分光実効面積率を求めるステップと、
任意の設定値の掛け合わせに対する分光反射率を、
上記の分光実効面積率を分光ノイゲバウア方程式に代入することにより求めるステップと、
観察光源、等色関数を用いて、前段のステップで求めた分光反射率から再現色の測色値を計算するステップと、
からなることを特徴とした印刷再現色予測方法。 In the method of estimating the spectral reflectance for an arbitrary multiplication of a printed material that reproduces the color by area modulation, and predicting the reproduced color,
Measuring spectral reflectance for the paper, each color solid, solid multiplication,
Measuring the spectral reflectance for any device control value of the primary color;
The spectral reflectance for an arbitrary device control value of the primary color is considered in consideration of the spectral reflectance of the paper, the solid spectral reflectance, and the dot gain, and the first virtual corresponding to the portion where the light passes only once through the ink layer As a sum of four spectral reflectances of the primary colors and the spectral reflectances of the two virtual primary colors corresponding to the part of the ink layer through which light passes only partially , weighted by the respective effective area ratios. Expressing and measuring the effective area ratio of the solid and virtual primary colors so as to minimize the error from the spectral reflectance with respect to the device control value in which the measured value of the measured primary color is present ;
The previous step is repeated for a device control value having a plurality of actual measurement values, and the relationship between the primary color device control value and the solid, the first virtual primary color, and the effective area ratio of the second virtual primary color is expressed by an approximate expression, A step for determining the relationship between the device control value and the effective area ratio;
Spectral reflectance for any device control value of the primary color,
Expressed as a linear sum of the spectral reflectance of the paper and the spectral reflectance of the solid as an expression different from the above expression ,
Assuming that the effective area ratio that is a coefficient of each of the solid and the paper is a spectral effective area ratio that is a function of the wavelength, obtaining the spectral effective area ratio ;
Spectral reflectance for multiplying any set value
Obtaining the spectral effective area ratio by substituting it into the spectral Neugebauer equation;
A step of calculating a colorimetric value of the reproduced color from the spectral reflectance obtained in the previous step using an observation light source and a color matching function;
A print reproduction color prediction method characterized by comprising:
一次色の中間調の分光反射率を、Ｒ＾（ｓ，λ）＝Ｒｉ（λ）ａ１（ｓ）＋Ｒo１（λ）ａ２（ｓ）＋Ｒo２（λ）ａ３（ｓ）＋Ｒｐ（λ）（１−ａ１（ｓ）−ａ２（ｓ）−ａ３（ｓ））（但し、ａ１（ｓ）は単色ベタの実効面積率、ａ２（ｓ）は第１の仮想的原色
の実効面積率、ａ３（ｓ）は第２の仮想的原色の実効面積率、Ｒｉ（λ）は単色ベタの分光反射率、Ｒｐ（λ）は用紙の分光反射率、Ｒo１（λ）は第１の仮想原色の分光反射率、Ｒo２（λ）は第２の仮想原色の分光反射率、Ｒ＾（ｓ，λ）はデバイス制御値ｓに対応する単色の推定分光反射率を示す）で示し、
実効面積率をｂ（ｓ，λ）＝（Ｒ＾（ｓ，λ）−Ｒｐ（λ））／（Ｒｉ（λ）−Ｒｐ（λ））としたとき、
分光ノイゲバウア方程式として以下の式を用いることを特徴とした印刷再現色予測方法。
Ｒ＾（λ）＝ｐｃ（λ）Ｒｃｉ（λ）＋ｐｍ（λ）Ｒｍｉ（λ）＋ｐｙ（λ）Ｒｙｉ（λ）＋ｐｋ（λ）Ｒｋｉ（λ）＋ｐｃｍ（λ）Ｒｃｍｉ（λ）＋ｐｍｙ（λ）Ｒｍｙｉ（λ）＋ｐｃｙ（λ）Ｒｃｙｉ（λ）＋ｐｃｋ（λ）Ｒｃｋｉ（λ）＋ｐｍｋ（λ）Ｒｍｋｉ（λ）＋ｐｙｋ（λ）Ｒｙｋｉ（λ）＋ｐｃｍｙ（λ）Ｒｃｍｙｉ（λ）＋ｐｃｍｋ（λ）Ｒｃｍｋｉ（λ）＋ｐｍｙｋ（λ）Ｒｍｙｋｉ（λ）＋ｐｃｙｋ（λ）Ｒｃｙｋｉ（λ）＋ｐｃｙｍｋ（λ）Ｒｃｙｍｋｉ（λ）＋ｐｐ（λ）Ｒｐ（λ）
ｐｃ（λ）＝ｂｃ（ｓｃ，λ）（１−ｂｍ（ｓｍ，λ））（１−ｂｙ（ｓｙ，λ））（１−ｂｋ（ｓｋ，λ））
ｐｍ（λ）＝（１−ｂｃ（ｓｃ，λ））ｂｍ（ｓｍ，λ）（１−ｂｙ（ｓｙ，λ））（１−ｂｋ（ｓｋ，λ））
ｐｙ（λ）＝（１−ｂｃ（ｓｃ，λ））（１−ｂｍ（ｓｍ，λ））ｂｙ（ｓｙ，λ）（１−ｂｋ（ｓｋ，λ））
ｐｋ（λ）＝（１−ｂｃ（ｓｃ，λ））（１−ｂｍ（ｓｍ，λ））（１−ｂｙ（ｓｙ，λ））ｂｋ（ｓｋ，λ）
ｐｃｍ（λ）＝ｂｃ（ｓｃ，λ）ｂｍ（ｓｍ，λ）（１−ｂｙ（ｓｙ，λ））（１−ｂｋ（ｓｋ，λ））
ｐｍｙ（λ）＝（１−ｂｃ（ｓｃ，λ））ｂｍ（ｓｍ，λ）ｂｙ（ｓｙ，λ）（１−ｂｋ（ｓｋ，λ））
ｐｃｙ（λ）＝ｂｃ（ｓｃ，λ）（１−ｂｍ（ｓｍ，λ））ｂｙ（ｓｙ，λ）（１−ｂｋ（ｓｋ，λ））
ｐｃｋ（λ）＝ｂｃ（ｓｃ，λ）（１−ｂｍ（ｓｍ，λ））（１−ｂｙ（ｓｙ，λ））ｂｋ（ｓｋ，λ）
ｐｍｋ（λ）＝（１−ｂｃ（ｓｃ，λ））ｂｍ（ｓｍ，λ）（１−ｂｙ（ｓｙ，λ））ｂｋ（ｓｋ，λ）
ｐｙｋ（λ）＝（１−ｂｃ（ｓｃ，λ））（１−ｂｍ（ｓｍ，λ））ｂｙ（ｓｙ，λ）ｂｋ（ｓｋ，λ）
ｐｃｍｙ（λ）＝ｂｃ（ｓｃ，λ）ｂｍ（ｓｍ，λ）ｂｙ（ｓｙ，λ）（１−ｂｋ（ｓｋ，λ））
ｐｃｍｋ（λ）＝ｂｃ（ｓｃ，λ）ｂｍ（ｓｍ，λ）（１−ｂｙ（ｓｙ，λ））ｂｋ（ｓｋ，λ）
ｐｍｙｋ（λ）＝（１−ｂｃ（ｓｃ，λ））ｂｍ（ｓｍ，λ）ｂｙ（ｓｙ，λ）ｂｋ（ｓｋ，λ）
ｐｃｙｋ（λ）＝ｂｃ（ｓｃ，λ）（１−ｂｍ（ｓｍ，λ））ｂｙ（ｓｙ，λ）ｂｋ（ｓｋ，λ）
ｐｃｙｍｋ（λ）＝ｂｃ（ｓｃ，λ）ｂｍ（ｓｍ，λ）ｂｙ（ｓｙ，λ）ｂｋ（ｓｋ，λ）
ｐｐ（λ）＝（１−ｂｃ（ｓｃ，λ））（１−ｂｍ（ｓｍ，λ））（１−ｂｙ（ｓｙ，λ））（１−ｂｋ（ｓｋ，λ）） The print reproduction color prediction method according to claim 1,
The spectral reflectance of the primary color halftone is represented by R ^ (s, λ) = R i (λ) a 1 (s) + R o1 (λ) a 2 (s) + R o2 (λ) a 3 (s) + R p (λ) (1-a 1 (s) -a 2 (s) -a 3 (s)) (where a 1 (s) is the effective area ratio of the solid monochrome, and a 2 (s) is the first The effective area ratio of the virtual primary color, a 3 (s) is the effective area ratio of the second virtual primary color, R i (λ) is the spectral reflectance of the solid monochrome, R p (λ) is the spectral reflectance of the paper, R o1 (λ) is the spectral reflectance of the first virtual primary color, R o2 (λ) is the spectral reflectance of the second virtual primary color, and R ^ (s, λ) is a monochromatic estimate corresponding to the device control value s. Spectral reflectivity)
When the effective area ratio is b (s, λ) = (R ^ (s, λ) −R p (λ)) / (R i (λ) −R p (λ)),
A printing reproduction color prediction method using the following equation as a spectral Neugebauer equation:
R ^ (λ) = p c (λ) R ci (λ) + p m (λ) R mi (λ) + p y (λ) R yi (λ) + p k (λ) R ki (λ) + p cm (λ ) R cmi (λ) + p my (λ) R myi (λ) + p cy (λ) R cyi (λ) + p ck (λ) R cki (λ) + p mk (λ) R mki (λ) + p yk (λ ) R yki (λ) + p cmy (λ) R cmyi (λ) + p cmk (λ) R cmki (λ) + p myk (λ) R myki (λ) + p cyk (λ) R cyki (λ) + p cymk (λ ) R cymki (λ) + p p (λ) R p (λ)
p c (λ) = b c (s c , λ) (1-b m (s m , λ)) (1-b y (s y , λ)) (1-b k (s k , λ))
p m (λ) = (1−b c (s c , λ)) b m (s m , λ) (1−b y (s y , λ)) (1−b k (s k , λ))
p y (λ) = (1−b c (s c , λ)) (1−b m (s m , λ)) b y (s y , λ) (1−b k (s k , λ))
p k (λ) = (1−b c (s c , λ)) (1−b m (s m , λ)) (1−b y (s y , λ)) b k (s k , λ)
p cm (λ) = b c (s c , λ) b m (s m , λ) (1-b y (s y , λ)) (1-b k (s k , λ))
p my (λ) = (1−b c (s c , λ)) b m (s m , λ) b y (s y , λ) (1−b k (s k , λ))
p cy (λ) = b c (s c , λ) (1−b m (s m , λ)) b y (s y , λ) (1−b k (s k , λ))
p ck (λ) = b c (s c , λ) (1−b m (s m , λ)) (1−b y (s y , λ)) b k (s k , λ)
p mk (λ) = (1−b c (s c , λ)) b m (s m , λ) (1−b y (s y , λ)) b k (s k , λ)
p yk (λ) = (1−b c (s c , λ)) (1−b m (s m , λ)) b y (s y , λ) b k (s k , λ)
p cmy (λ) = b c (s c , λ) b m (s m , λ) b y (s y , λ) (1−b k (s k , λ))
p cmk (λ) = b c (s c , λ) b m (s m , λ) (1−b y (s y , λ)) b k (s k , λ)
p myk (λ) = (1−b c (s c , λ)) b m (s m , λ) b y (s y , λ) b k (s k , λ)
p cyk (λ) = b c (s c , λ) (1−b m (s m , λ)) b y (s y , λ) b k (s k , λ)
p cymk (λ) = b c (s c , λ) b m (s m , λ) b y (s y , λ) b k (s k , λ)
p p (λ) = (1-b c (s c , λ)) (1-b m (s m , λ)) (1-b y (s y , λ)) (1-b k (s k , Λ))
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2013072945A JP6155766B2 (en) | 2013-03-29 | 2013-03-29 | Print reproduction color prediction method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2013072945A JP6155766B2 (en) | 2013-03-29 | 2013-03-29 | Print reproduction color prediction method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2014197790A JP2014197790A (en) | 2014-10-16 |
| JP6155766B2 true JP6155766B2 (en) | 2017-07-05 |
Family
ID=52358312
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2013072945A Active JP6155766B2 (en) | 2013-03-29 | 2013-03-29 | Print reproduction color prediction method |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6155766B2 (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017178329A1 (en) | 2016-04-12 | 2017-10-19 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Audio encoder for encoding an audio signal, method for encoding an audio signal and computer program under consideration of a detected peak spectral region in an upper frequency band |
| CN106483080A (en) * | 2016-10-12 | 2017-03-08 | 东莞市美盈森环保科技有限公司 | A kind of printed matter color method for evaluating quality |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002016806A (en) * | 2000-06-30 | 2002-01-18 | Fuji Photo Film Co Ltd | Method and device for creating gradation correction curve |
| JP4725190B2 (en) * | 2005-05-25 | 2011-07-13 | 凸版印刷株式会社 | Color information processing method, color information processing apparatus, and program thereof |
| JP4586668B2 (en) * | 2005-08-01 | 2010-11-24 | 凸版印刷株式会社 | Color information processing method, color information processing apparatus, and program thereof |
-
2013
- 2013-03-29 JP JP2013072945A patent/JP6155766B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2014197790A (en) | 2014-10-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR101148381B1 (en) | Method for generating a color chart | |
| JP5381380B2 (en) | Image processing apparatus, image forming system, and program | |
| JP6156401B2 (en) | Color conversion method, program, and image processing apparatus | |
| JP2017195522A (en) | Image formation apparatus and program for color conversion control, and method for color conversion control | |
| JP5644345B2 (en) | Color conversion apparatus, image forming apparatus, and color conversion program | |
| US9122978B2 (en) | Image processing device and method, and non-transitory computer readable medium | |
| JP6155766B2 (en) | Print reproduction color prediction method | |
| JP4725190B2 (en) | Color information processing method, color information processing apparatus, and program thereof | |
| JP5152515B2 (en) | Color correction coefficient generation apparatus and color correction coefficient generation program | |
| JP4600673B2 (en) | Color gamut calculation method and color gamut calculation device | |
| JP5699608B2 (en) | Image processing apparatus, image processing method, and image processing program | |
| US7944595B2 (en) | Color processing method and apparatus | |
| JP4123354B2 (en) | Calibration apparatus and calibration method | |
| JP5324405B2 (en) | Color processing apparatus and method, and image forming apparatus | |
| JP4586668B2 (en) | Color information processing method, color information processing apparatus, and program thereof | |
| JP2007251373A (en) | Method for making color estimation model | |
| JP6717058B2 (en) | Reproduction color profile correction system, reproduction color profile correction method and program | |
| US10836181B2 (en) | Color prediction method | |
| JP2005047186A (en) | Image processor and its method | |
| JP2014143522A (en) | Printed matter spectral reflection factor prediction method and printed matter spectral reflection factor predictor | |
| JP2010130485A (en) | Color processing device and program | |
| US8437541B2 (en) | Color correction coefficient generating apparatus, color correction coefficient generation method, and computer readable medium storing program | |
| JP2008017293A (en) | Image processor and image processing program | |
| JP2007074521A (en) | Calibration method and image processor | |
| JP5942896B2 (en) | Color processing apparatus, color processing method, and program |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20160219 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20161220 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20170131 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20170328 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20170509 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20170522 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6155766 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |