JP3362281B2 - Color reproduction device and color reproduction method - Google Patents

Color reproduction device and color reproduction method

Info

Publication number
JP3362281B2
JP3362281B2 JP26671892A JP26671892A JP3362281B2 JP 3362281 B2 JP3362281 B2 JP 3362281B2 JP 26671892 A JP26671892 A JP 26671892A JP 26671892 A JP26671892 A JP 26671892A JP 3362281 B2 JP3362281 B2 JP 3362281B2
Authority
JP
Japan
Prior art keywords
color
target value
points
value
color reproduction
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
JP26671892A
Other languages
Japanese (ja)
Other versions
JPH0698161A (en
Inventor
博哲 洪
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.)
Konica Minolta Inc
Original Assignee
Konica Minolta Inc
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 Konica Minolta Inc filed Critical Konica Minolta Inc
Priority to JP26671892A priority Critical patent/JP3362281B2/en
Publication of JPH0698161A publication Critical patent/JPH0698161A/en
Application granted granted Critical
Publication of JP3362281B2 publication Critical patent/JP3362281B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】本発明は、色座標が与えられたと
きの、該色座標に対応した入力色分解画像信号を正確に
推定する場合に適用して好適な色再現装置及び色再現方
法に関し、特にグレー色の階調性の改善を図ったもので
ある。
The present invention relates, when the color coordinates are given, preferred color reproduction device and the color reproduction direction applied to a case of accurately estimating the input color separation image signals corresponding to the color coordinate
Related to the law, and in particular for Improving the gradation of gray color.

【0002】[0002]

【従来の技術】印刷や、感熱転写、インクジェット、電
子写真等によってフルカラープリントする場合、カラー
プリントの出力色としては、一般にY,M,C,Kの4
色が使われることが多い。4色で色を表現するため、こ
のY,M,C,Kの組合せをどのようにするかについて
は、様々な手法がある。例えばこの手法の一つとして、
UCR(Under Color Removal )と呼ばれているものが
ある。
2. Description of the Related Art In the case of full color printing by printing, thermal transfer, ink jet, electrophotography, etc., the output colors of the color print are generally Y, M, C and K.
Color is often used. There are various methods for how to combine the Y, M, C, and K in order to express the colors with four colors. For example, as one of the methods,
There is a so-called UCR (Under Color Removal).

【0003】その具体的な方法として100 %UCR法と
呼ばれているものは、図29に示すように元のY,M,C
(YO , MO , CO )のうちの最小濃度を黒色Kの濃度
に置換する手法であって、同図AのようにYO が最小濃
度であるときには、100 %UCR法で置換すると、同図
Bに示すようになる。置換後の濃度は夫々Yn , Mn,
n , Kn (但しYn =0)となる。
As a concrete method, a method called 100% UCR method is the original Y, M, C as shown in FIG.
A method of replacing the minimum density of (Y O , M O , C O ) with the density of black K, and when Y O is the minimum density as shown in FIG. , As shown in FIG. The concentrations after substitution are Y n , M n , and
C n, the K n (where Y n = 0).

【0004】このような100 %UCR法によると、以下
のような特性を有することが知られている。 ・作業上の信頼性が増す。 ・インキの節約が図れる。 ・多少のインキの不安定性をカバーできる。
It is known that the 100% UCR method has the following characteristics. -Working reliability is increased.・ Ink can be saved.・ Can cover some ink instability.

【0005】・インキ乾燥のエネルギー節約と、乾燥
関する諸問題が減る。これは、Y,M,Cのインキ量の
和に対し、Kのインキ量は大方1/3で対応できるため
である。一方、Y,M,C,Kの色再現特性を求める場
合、Y,M,C,Kの離散的な組み合わせによる少数の
カラーパッチを実際に作成し、このカラーパッチを実際
に測色し、その測色値を補間演算して上述した組み合わ
せ以外のY,M,C,Kの色再現特性、つまり、ある
Y,M,C,Kの組み合わせのときの測色値を推定する
ようにすることが考えられている。
The energy savings of ink drying and the problems associated with drying are reduced. This is because the K ink amount can be roughly ⅓ to the sum of the Y, M, and C ink amounts. On the other hand, when obtaining the color reproduction characteristics of Y, M, C, and K, a small number of color patches by a discrete combination of Y, M, C, and K are actually created, and this color patch is actually measured, The colorimetric values are interpolated to estimate the color reproduction characteristics of Y, M, C, K other than the above-mentioned combinations, that is, the colorimetric values for a certain Y, M, C, K combination are estimated. Is being considered.

【0006】また、これとは逆にある特定の色座標が指
定されたときには、その色座標を示すY,M,C,Kの
組み合わせも、カラーパッチから実際に測色した値を補
間演算して推定することができる。このようにY,M,
C,Kの全ての組み合わせを、Y,M,C,Kの離散的
な組み合わせに基づいて作成されたカラーパッチを用い
て推定すれば、実際の測色値からY,M,C,Kの組み
合わせを推定するものであるから、その推定精度が向上
し、色再現特性がより改善されることになる。
On the contrary, when a specific color coordinate is designated, the combination of Y, M, C, and K indicating the color coordinate is also interpolated by the value actually measured from the color patch. Can be estimated. Thus, Y, M,
If all combinations of C, K are estimated using color patches created based on discrete combinations of Y, M, C, K, the actual colorimetric values of Y, M, C, K Since the combination is estimated, the estimation accuracy is improved and the color reproduction characteristic is further improved.

【0007】なお、Y,M,C,Kの離散的な組み合わ
せによって作成されるカラーパッチの一例を図30に示
す。同図はY,M,C,Kの最大値を夫々255 とし、略
64ステップで5段階に分けて色の濃度を取り出したとき
のカラーパッチである。また、このようにある色座標が
与えられたとき、若しくはその色座標に相当する電気信
号 (例えばR,G,Bの色分解画像信号) が与えられた
ときには、何れの場合でも、これらよりY,M,C,K
の組み合わせを得るための色出力修正手段が必要にな
る。
FIG. 30 shows an example of a color patch created by a discrete combination of Y, M, C and K. In the figure, the maximum values of Y, M, C, and K are set to 255, and
This is a color patch when the density of the color is extracted in 5 steps in 64 steps. In addition, when a certain color coordinate is given, or when an electric signal corresponding to the color coordinate (for example, R, G, B color separation image signals) is given, Y is more than Y in any case. , M, C, K
A color output correction means for obtaining the combination of is required.

【0008】図31はその一例であって、プリンタ12の前
段に色分解画像修正装置の一つであるカラーマスキング
装置10が設けられ、ここで色座標に相当するR,G,B
の色分解画像信号がY,M,C,Kの色信号に変換され
た後、プリンタ12に供給され、その画像がY,M,C,
Kによって記録媒体13上に記録される。ところで、上述
した100 %UCR法では、上述したような利点を有する
ものの、 ・全般に彩度が低下する。
FIG. 31 shows an example thereof, in which a color masking device 10, which is one of color separation image correction devices, is provided in front of the printer 12, where R, G, B corresponding to color coordinates are provided.
After the color separation image signals of Y, M, C, and K are converted into Y, M, C, and K color signals, the images are supplied to the printer 12, and the images thereof are Y, M, C,
It is recorded on the recording medium 13 by K. By the way, the above-mentioned 100% UCR method has the above-mentioned advantages, but the saturation is generally lowered.

【0009】・全体に黒版の割合が増えてくると濃度が
低下する。 などの欠点を有し、色再現特性が劣化してしまう。そこ
で、本願出願人は、100 %UCR法の欠点である色再現
特性を改善すると共に、カラーパッチによる色再現特性
推定方法を適用するに際して、より少ないカラーパッチ
を使用して色再現特性を推定できるようにしたものを先
に提案した(特開平2−136848号公報参照) 。
The density decreases as the proportion of black plates increases. However, the color reproduction characteristic is deteriorated. Therefore, the applicant of the present application can improve the color reproduction characteristic, which is a drawback of the 100% UCR method, and can estimate the color reproduction characteristic using less color patches when applying the color reproduction characteristic estimation method using color patches. The above-mentioned method was previously proposed (see Japanese Patent Laid-Open No. 2-136848).

【0010】このものでは、与えられた色座標を示す
Y,M,C,Kの組み合わせの中、Kの値が最小濃度値
となり、Y,M,C,Kの値が最小濃度値となるY,
M,C,Kの組み合わせを、離散的に与えられた色座標
の各特定点について夫々求めると共に、これら離散的な
組み合わせによるY,M,C,Kのカラーパッチを作成
し、作成されたカラーパッチの測色値から、上記特定点
以外の色座標に対応するY,M,C,Kの濃度値を推定
するようにしたものである。
In this system, among the combinations of Y, M, C, and K indicating the given color coordinates, the value of K is the minimum density value and the value of Y, M, C, and K is the minimum density value. Y,
A combination of M, C, K is obtained for each specific point of discretely given color coordinates, Y, M, C, K color patches are created by these discrete combinations, and the created color is created. The density values of Y, M, C and K corresponding to the color coordinates other than the specific point are estimated from the colorimetric values of the patch.

【0011】任意の色座標が与えられたとき、その色を
表すY,M,C,Kの組み合わせは一般的には無数に存
在する。この無数に存在する組み合わせを唯一のものと
するため、与えられた色座標を示す組み合わせの中、K
の値が最大濃度値をとるという条件が導入される。この
条件は、Y,M,Cの少なくとも一つが0か若しくはK
が最大濃度 (印刷でいうベタ) という条件となる。この
ような「ある色を表すことができる範囲でKの値が最大
濃度値をとる」という条件 (改良型条件) を入れた改良
型によれば、カラーパッチ数を最小限に抑制しつつ色再
現特性を改善できる。
When arbitrary color coordinates are given, there are generally innumerable combinations of Y, M, C and K representing the color. In order to make this innumerable combination unique, K among the combinations showing the given color coordinates
The condition that the value of is the maximum density value is introduced. This condition is that at least one of Y, M and C is 0 or K
Is the maximum density (solid in printing). According to the improved type that includes the condition (improved type condition) that the value of K takes the maximum density value within the range that can express a certain color, the number of color patches is suppressed to the minimum. Reproducibility can be improved.

【0012】[0012]

【発明が解決しようとする課題】しかしながら、かかる
改良型100 %UCR法においては、グレー色の再現の場
合は、前記条件に従うとK一色で再現されることになる
が、その場合、電子写真プリンタやインクジェットプリ
ンタのように階調性が悪いものでは、単一色では階調性
が悪いという欠点が直接的に現れてしまい、良好なグラ
デーションを得ることができないという問題を生じた。
However, in the improved 100% UCR method, in the case of gray color reproduction, K color is reproduced in accordance with the above conditions. In that case, the electrophotographic printer is used. In the case of poor gradation such as an ink jet printer and an ink jet printer, a defect that gradation is poor with a single color directly appears, and a problem that good gradation cannot be obtained occurs.

【0013】本発明は、このような問題点に鑑みなされ
たもので、Y,M,C,Kで作りうる色域をフルに使い
ながら、階調性の悪いプリンタのグラデーションを良好
に出すためのY,M,C,K値の決定方法を提供するこ
とを目的とする。
The present invention has been made in view of the above problems, and in order to produce good gradation of a printer having poor gradation, while making full use of the color gamut that can be created by Y, M, C, and K. It is an object of the present invention to provide a method for determining Y, M, C, K values of.

【0014】[0014]

【課題を解決するための手段】このため本発明は、
(イエロー) 、M (マゼンタ) ,C (シアン) 及びK
(黒) の4色を組み合わせて色再現を行う色再現方法で
あって、与えられた色座標を示すY、M、C、Kの組み
合わせのうち、Y、M、Cの少なくとも一つが最大とな
る場合にKを含めて目標値を求める処理、およびKを0
とした場合にY、M、Cのみで目標値を求める処理を用
いることによって、Kの濃度が最小となるY,M,C,
Kの組み合わせに対応するY、M、C、Kの濃度値を推
定して色再現を行うことを特徴とするものである。
Means for Solving the Problems The present invention for this purpose is, Y
(Yellow), M (Magenta), C (Cyan) and K
A color reproduction method that combines four colors (black) to reproduce colors.
Yes, a set of Y, M, C, K indicating the given color coordinates
Of the combinations, at least one of Y, M, and C is the maximum.
If K is included, the process of obtaining the target value including K, and K is set to 0
In case of, the process of obtaining the target value using only Y, M, and C is used.
, The concentration of K becomes the minimum, Y, M, C,
Estimate the density values of Y, M, C, K corresponding to the combination of K
The feature is that the color reproduction is performed in a fixed manner.

【0015】[0015]

【作用】色座標を示すY,M,C,Kの組み合わせは無
数に存在するが、そのうち、Kの濃度を最小のものを決
める条件は、K=0又はY,M,Cの少なくとも1つが
最大濃度であることと同等である。このように、Kの濃
度を最小とすることで、グレー色をY,M,Cの少なく
とも3色若しくはそれにKを加えた4色を混合して再現
することができるため単一色での階調性の悪さをカバー
できる。
There are innumerable combinations of Y, M, C and K indicating color coordinates, and the condition for determining the minimum concentration of K is that K = 0 or at least one of Y, M and C. It is equivalent to the maximum concentration. As described above, by minimizing the density of K, it is possible to reproduce a gray color by mixing at least three colors of Y, M, and C or four colors obtained by adding K to the gray color. It can cover bad sex.

【0016】また、1色については決められ、残り3色
についての組み合わせを決めればよいため、Kの濃度を
最大にする場合と同様にカラーパッチの数を十分少なく
作成すれば済む。
Since one color can be determined and combinations for the remaining three colors can be determined, the number of color patches can be made sufficiently small as in the case of maximizing the density of K.

【0017】[0017]

【実施例】以下に本発明の実施例を図に基づいて説明す
る。前述したように任意の色座標が与えられたときにそ
の色を表すY,M,C,Kの組み合わせは一般には無数
に存在する。この無数に存在する組み合わせを唯一にす
ると同時に、グレー色を複数色を混合して再現するため
に、与えられた色座標を示す組み合わせのうちKの濃度
を最小とするものを採るという条件を導入する。
Embodiments of the present invention will be described below with reference to the drawings. As described above, there are generally innumerable combinations of Y, M, C, and K that represent a given color when given an arbitrary color coordinate. In order to reproduce this gray color by mixing a plurality of colors while introducing only one of these innumerable combinations, the condition that the density of K is minimized among the combinations showing the given color coordinates is introduced. To do.

【0018】即ち、色座標を示す範囲内でもKは色々な
値を採り得るが、採り得る値のうちで最小値をKの値と
して決めれば、これによってその他のY,M,Cの組み
合わせを決定することができ、かつ、Kを最小とするこ
とで、グレー色をY,M,Cの少なくとも3色若しくは
それにKを加えた4色を混合して再現することができる
ため単一色での階調性の悪さをカバーできるのである。
That is, K can take various values even within the range showing the color coordinates, but if the minimum value of the possible values is determined as the value of K, other combinations of Y, M, and C can be obtained. The gray color can be determined, and the gray color can be reproduced by mixing at least three colors Y, M, and C or four colors obtained by adding K to the gray color. It is possible to cover the poor gradation.

【0019】この条件は、Y,M,Cの何れか一つが最
大濃度 (相対値ではなく絶対量として) であるか、Kが
0である場合、つまり、 ・Yが最大 (255)のとき ・Mが最大 (255)のとき ・Cが最大 (255)のとき ・Kが0のとき ・・・ (1) の少なくとも1つが満たされる条件に他ならない。そし
て、これらによって決定される色立体が色再現範囲とな
る。
This condition is that when any one of Y, M and C is the maximum density (as an absolute amount, not as a relative value) or when K is 0, that is, when Y is maximum (255)・ When M is the maximum (255) ・ When C is the maximum (255) ・ When K is 0 ... It is the condition that at least one of (1) is satisfied. The color solid determined by these becomes the color reproduction range.

【0020】尚、Y,M,Cの何れか一つが最大濃度で
ある条件を有するのは、これらが最大の場合でもKが0
になるとは限らないからである。この方法によってカラ
ーパッチを作成すると図1に示すようになる。同図は、
Y,M,C,Kの離散的なポイント数nを5、最大の量
子化レベルを256 ステップに夫々選択した場合であっ
て、このときには夫々のポイント (0,64, 128,192,
255) の基本色Y,M,C,Kを組み合わせて実際に記
録媒体、例えば印刷紙上にインキにより記録すると図1
のカラーパッチが得られる。
It is to be noted that the condition that any one of Y, M, and C has the maximum density is that K is 0 even when these are maximum.
This is not always the case. When a color patch is created by this method, it becomes as shown in FIG. This figure shows
When the number of discrete points n of Y, M, C, and K is selected to 5, and the maximum quantization level is selected to 256 steps, at this time, the respective points (0, 64, 128, 192,
When the basic colors Y, M, C, K of (255) are actually combined and recorded on the recording medium, for example, printing paper with ink,
To obtain a color patch of.

【0021】この図1のカラーパッチを実際に測色し
て、その測色値 (与えられた色座標)を他の表色系 (例
えばCIELUV表色系) に関する変換式を使用してC
IELUV表色系の値に換算し、これをカラーパッチ毎
にプロットすると、図2のようになる。各カラーパッチ
の測色値が夫々の格子点に対応する。ただし、この図2
は説明の便宜上、彩度と明度との2軸上で表現したもの
で、また、シアンCの値を省略して描いてある。以下に
示す表色系も同様である。
The color patch of FIG. 1 is actually measured and its colorimetric value (given color coordinates) is converted into C using a conversion formula for another color system (for example, CIE LUV color system).
The values are converted into values in the IELUV color system and plotted for each color patch, as shown in FIG. The colorimetric value of each color patch corresponds to each grid point. However, this Figure 2
For convenience of explanation, is expressed on the two axes of saturation and lightness, and the value of cyan C is omitted. The same applies to the color system shown below.

【0022】上述したカラーパッチは図3のようにする
こともできる。この図3によれば、カラーパッチ数が図
1よりさらに減少している。図3のカラーパッチを実際
に測色してCIELUV表色系に写像すると、図4のよ
うになる。黒丸印が実際に測色して得た値である。ここ
で図4の表色系において、Kの値が大きくなるにしたが
って格子間が略直線的に縮小されていくので、Kが大き
いときには白丸格子点は直前の格子点の値に基づいて直
線補間 (例えば内分) して求めても、その誤差は僅少に
なる。
The color patch described above may be arranged as shown in FIG. According to this FIG. 3, the number of color patches is further reduced as compared with FIG. When the color patches of FIG. 3 are actually measured and mapped on the CIE LUV color system, the result is as shown in FIG. The black circles are the values obtained by actual color measurement. Here, in the color system of FIG. 4, the space between the grids is reduced substantially linearly as the value of K increases, so when K is large, the white circle grid points are linearly interpolated based on the values of the immediately preceding grid points. Even if it is calculated by (for example, internal division), the error becomes small.

【0023】白丸格子点の全てを補間すると、図2の表
色系と同一になるがこれはカラーパッチとして図1のカ
ラーパッチを使用したときと同一である。つまり、単純
に (1)式に合致するだけでなく、測色値の性質を利用し
てKの値が増すほどカラーパッチを減らしても補間演算
処理すれば図1と同一になるから、カラーパッチ数を削
減できる。
When all the white circle lattice points are interpolated, they become the same as the color system of FIG. 2, but this is the same as when the color patch of FIG. 1 is used as the color patch. In other words, not only the equation (1) is simply met, but if the color patch is reduced as the value of K is increased by utilizing the property of the colorimetric value, the interpolation calculation processing will be the same as in FIG. The number of patches can be reduced.

【0024】電気的な処理によって格子間隔を上述より
もさらに1/2 にするには、写像された値に基づいて補間
演算すればよい。この場合の補間処理は非線形な補間処
理である。補間処理の一例を図5に示す。図5に示すよ
うに、黒丸●を格子点 (サンプル点) とし、△印と×印
が夫々補間すべき点とすると、×印のように前後に2点
ずつ格子点が存在する場合と、×印のように前後に1点
及び3点ある場合とでは、異なった補間式が使用され
る。
In order to further reduce the lattice spacing by a half by the electrical processing, interpolation calculation may be performed based on the mapped values. The interpolation process in this case is a non-linear interpolation process. An example of the interpolation process is shown in FIG. As shown in FIG. 5, if the black circles ● are grid points (sample points) and the triangles and the crosses are points to be interpolated, respectively, there are two grid points before and after as shown by the crosses. Different interpolation formulas are used in the case where there are one point and three points before and after as in the case of the cross mark.

【0025】補間すべき点の表色系を、CIELUVと
し、各サンプル点の表色系の値をLi’, Ui’, V
i’ (i=1〜4) としたとき、前者の場合は以下のよ
うな補間式によって補間される。 Lm ’=− (1/16) L1’+(9/16)L2’ + (9/16) L3’−(1/16)L4’ Um ’=− (1/16) U1’+(9/16)U2’ + (9/16) U3’−(1/16)U4’ Vm ’=− (1/16) V1’+(9/16)V2’ + (9/16) V3’−(1/16)V4’ 後者の場合には、次の補間式が使用される。
CIE LUV is the color system of the points to be interpolated, and the values of the color system of each sample point are Li ', Ui', V
When i '(i = 1 to 4), the former case is interpolated by the following interpolation formula. L m '=-(1/16) L1' + (9/16) L2 '+ (9/16) L3'- (1/16) L4' U m '=-(1/16) U1' + ( 9/16) U2 '+ (9/16) U3'-(1/16) U4 'V m ' =-(1/16) V1 '+ (9/16) V2' + (9/16) V3 ' -(1/16) V4 'In the latter case, the following interpolation formula is used.

【0026】 Lm ’= (5/16) L1’+(15/16) L2’ − (5/16) L3’− (1/16) L4’ Um ’= (5/16) U1’+(15/16) U2’ − (5/16) U3’− (1/16) U4’ Vm ’= (5/16) V1’+(15/16) V2’ − (5/16) V3’− (1/16) V4’ 補間処理の順序の一例を図6に示す。番号I、II、III
の順序で補間される。
L m '= (5/16) L1' + (15/16) L2'- (5/16) L3'- (1/16) L4 'U m ' = (5/16) U1 '+ (15/16) U2'- (5/16) U3'- (1/16) U4 'V m ' = (5/16) V1 '+ (15/16) V2'- (5/16) V3' -(1/16) V4 'An example of the order of the interpolation processing is shown in FIG. Numbers I, II, III
Are interpolated in the order of.

【0027】このような補間処理によって、実際に測色
した数より多くのカラーパッチ数を使用したのと同様な
表色系の格子点数が得られる。補間処理は、直接近似に
よる内挿処理でもよい。以上のようにすれば、カラーパ
ッチ数を増やさないで、測色値に対応した色座標の値を
得ることができる。その場合の色座標の値は何れも格子
点のデータである。
By such an interpolation process, the number of grid points of the color system similar to the number of color patches actually used can be obtained. The interpolation processing may be interpolation processing by direct approximation. According to the above, the color coordinate values corresponding to the colorimetric values can be obtained without increasing the number of color patches. The values of the color coordinates in that case are data of grid points.

【0028】格子点以外に存在する色座標の値は、以下
のような方法によって推定される。まず、当該推定方法
の原理を説明する。N次元(Nは2以上の整数)の空間
を分割する最も単純な分割空間は、(N+1)個の頂点
を有する分割空間である。例えば、二次元では三角形、
三次元では三角錐である。ここで、図7に示すように空
間1及び空間2を考える。ここでは二次元で考える。こ
の場合、空間1及び空間2は、夫々a〜i及び点a’〜
i’を用いて三角形に分割されており、例えば空間1の
分割空間Δbdeは、空間2の分割空間Δb’d’e’
と対応している。
The values of color coordinates existing at points other than the grid points are estimated by the following method. First, the principle of the estimation method will be described. The simplest division space that divides an N-dimensional (N is an integer of 2 or more) space is a division space having (N + 1) vertices. For example, in two dimensions a triangle,
In three dimensions it is a triangular pyramid. Here, space 1 and space 2 are considered as shown in FIG. Here we consider two dimensions. In this case, the space 1 and the space 2 are a to i and a point a'to respectively.
i ′ is used to divide into triangles. For example, the divided space Δbde of the space 1 is divided into the divided space Δb′d′e ′ of the space 2.
It corresponds to.

【0029】このように対応している分割空間内では線
形に対応していると仮定し、空間1に与えられた点Pは
空間2では点P’に対応させられる。ここで、例えば点
P及び点P’を囲む点を夫々Pi(xi,yi,zi,
・・・)、Pi’(xi’,yi’,zi’,・・・)
(i=1〜N+1)とし、与えられた点Pを(x,
y,z・・・)、求むべき点P’を(x’、y’,
z’,・・・)とすると、次のようにマトリクス形式で
示すことができる。
Assuming that they correspond linearly in the corresponding divided spaces, the point P given to the space 1 is made to correspond to the point P'in the space 2. Here, for example, points surrounding the point P and the point P ′ are respectively Pi (xi, yi, zi,
...), Pi '(xi', yi ', zi', ...)
(I = 1 to N + 1), and the given point P is (x,
y, z ...), and the point P ′ to be obtained is (x ′, y ′,
z ′, ...), the matrix form can be shown as follows.

【0030】[0030]

【数1】 [Equation 1]

【0031】また、求むべき点Pは、この点Pを囲む3
点から重み平均でも求めることができる。図8に示すよ
うに、S(b): S(d): S(e)=S(b’): S
(d’): S(e’)となる。尚、与えられた点Pがど
の分割空間に入るかは、夫々の分割空間の境界線(また
は境界面)のどちら側に入っているかを調べることで特
定できる。
Further, the point P to be obtained is 3 surrounding the point P.
The weighted average can be obtained from the points. As shown in FIG. 8, S (b): S (d): S (e) = S (b '): S
(D '): S (e'). It should be noted that which division space the given point P enters can be specified by checking which side of the boundary line (or boundary surface) of each division space it is.

【0032】本例においては、上述した原理を用いるこ
とにより、ある目標色を示す基本色の組合せを以下のよ
うにして求める。簡単のため、本例においても基本色を
2色(例えば、Y, M)として説明する。前記補間処理
又は内挿処理後のY, M座標系及びCIELUV表色系
は、夫々図9, 図10に示すようになる。同図において白
○の格子点は内挿された点を示している。
In this example, by using the above-mentioned principle, a combination of basic colors showing a certain target color is obtained as follows. For simplification, the basic color will be described as two colors (for example, Y and M) also in this example. The Y, M coordinate system and the CIE LUV color system after the interpolation processing or the interpolation processing are as shown in FIGS. 9 and 10, respectively. In the figure, the grid points of white circles indicate the interpolated points.

【0033】次にY, M座標系の空間及びCIELUV
表色系の空間を、夫々図11及び図12に示すように三角形
に分割する。この結果82 ×2=128 の三角形に分割さ
れる。次に、図13に示すように、CIELUV表色系に
求めたい出力色(目標色)に対応する目標値T’が上述
のように夫々の三角形の境界線のどちら側に入っている
かを調べることで三角形が特定される。
Next, the space of the Y and M coordinate systems and CIE LUV
The color space is divided into triangles as shown in FIGS. 11 and 12, respectively. As a result, it is divided into 8 2 × 2 = 128 triangles. Next, as shown in FIG. 13, it is checked which side of the boundary line of each triangle the target value T'corresponding to the output color (target color) desired in the CIE LUV color system belongs to. This will identify the triangle.

【0034】ここで、目標値T’が、図14に示すように
格子点a’〜c’で形成される三角形は、図13に示すよ
うに格子点a〜cで形成される三角形となる。次に目標
値T’が入っているCIELUV表色系の三角形及び目
標値T’が入っているY, M座標系の三角形の各頂点
(3点ずつ)の座標と、目標値T’を(2) 式に代入し
て、目標値T、したがって求めたい出力色を示す基本色
の組合せを求める。
Here, the triangle whose target value T'is formed by the grid points a'-c 'as shown in FIG. 14 is the triangle formed by the grid points a-c as shown in FIG. . Next , the coordinates of each vertex (three points) of the triangle of the CIE LUV color system containing the target value T'and the triangle of the Y , M coordinate system containing the target value T ', and the target value T' ( 2) Substituting into the equation, the target value T, and thus the basic color combination indicating the desired output color is calculated.

【0035】以下、与えられた目標値T’の夫々で三角
形の特定及び(2) 式の演算を実行して、目標値Tを求め
る。与えられる目標値T’が図15に示すようにCIEL
UV表色系のどの三角形にも入らないとき、つまり色再
現範囲外にあるときには、この目標値T’を色再現範囲
内に移動する必要がある。
Below, the target value T is obtained by specifying the triangle and performing the operation of the equation (2) for each of the given target values T '. The given target value T'is CIEL as shown in FIG.
When it does not fall into any of the triangles of the UV color system, that is, when it is outside the color reproduction range, it is necessary to move this target value T ′ into the color reproduction range.

【0036】この場合には、図16に示すように、目標値
T’を無彩色方向に移動させ、図18に示すように目標値
T’を無彩色方向の直線と色再現範囲の境界との交点の
座標を目標値T’とする。そして、この目標値T’とそ
れが含まれる線分L’を算出すると共に図17に示すよう
にY, M座標系で線分L’に対応する線分Lを算出す
る。これにより、(1) 式をもって目標値Tを算出する。
In this case, the target value T'is moved in the achromatic color direction as shown in FIG. 16, and the target value T'is set as a boundary between the straight line in the achromatic color direction and the color reproduction range as shown in FIG. The coordinate of the intersection point of is the target value T '. Then, the target value T ′ and the line segment L ′ including the target value T ′ are calculated, and the line segment L corresponding to the line segment L ′ in the Y , M coordinate system is calculated as shown in FIG. As a result, the target value T is calculated using the equation (1).

【0037】ところで、上述例は説明の簡単のため、基
本色を2色として説明したが、基本色が3色(Y,M,
C)の場合であっても目標値T(Y, M, C)を同様に
して求めることができる。但し、この場合には空間を三
角錐に分割し、与えられた目標値T’がCIELUV
表色系のどの三角錐に入っているかを調べた後、対応す
るY,M,C座標系の三角錐を決定する。
In the above example, the basic colors are two for the sake of simplicity of explanation, but the basic colors are three (Y, M,
Even in the case of C), the target value T (Y, M, C) can be similarly obtained. However, in this case, the space is divided into triangular pyramids, and the given target value T ′ is CIELUV.
After checking which triangular pyramid of the color system is included, the corresponding triangular pyramid of the Y, M, C coordinate system is determined.

【0038】そして、(2) 式をもって目標値Tを演算す
る。ここで、例えば図19に示すように目標値T (Y,
M,C) が入っているY,M,C座標系の三角錐の4頂
点を、 (Y1, M1, C1) (Y2, M2, C2) (Y3, M3, C3) (Y4, M4, C4) とし、目標値T’ (L, U, V) の入っているCIEL
UV 表色系の三角錐の4頂点を、 (L1, U1, V1) (L2, U2, V2) (L3, U3, V3) (L4, U4, V4) とすると、目標値Tは次式で求めることができる。
Then, the target value T is calculated using the equation (2). Here, for example, as shown in FIG. 19, the target value T (Y,
(Y, M, C1) (Y2, M2, C2) (Y3, M3, C3) (Y4, M4, C4 ) And the target value T '(L, U, V) is entered in the CIEL
If the four vertices of the triangular pyramid of the UV color system are (L1, U1, V1) (L2, U2, V2) (L3, U3, V3) (L4, U4, V4), the target value T is You can ask.

【0039】[0039]

【数2】 [Equation 2]

【0040】以上は、逆行列として4×4のマトリクス
を用いる例を示したが、次のように3×3のマトリクス
を用いても等価であり、繰り返しの演算が少なくなる分
誤差を減少できる。
In the above, an example in which a 4 × 4 matrix is used as an inverse matrix has been shown, but the same is true even if a 3 × 3 matrix is used as follows, and the error can be reduced by the number of repeated operations. .

【0041】[0041]

【数3】 [Equation 3]

【0042】また、目標値T’(L,U,V)が頂点
(L1,U1,V1)、 (L2,U2,V2)、 (L
3,U3,V3)、 (L4,U4,V4) を有した三角
錐の内側にあるか否かの判別は次のようにして行える。
即ち、
Further, the target value T '(L, U, V) is the vertex
(L1, U1, V1), (L2, U2, V2), (L
3, U3, V3), (L4, U4, V4) can be discriminated as to whether or not it is inside the triangular pyramid.
That is,

【0043】[0043]

【数4】 [Equation 4]

【0044】で得られるα,β,γが次式を満たすとき
当該三角錐の内側 (表面上を含む) にあり、満たさない
場合は、外側にある。 α≧0,β≧0,γ≧0,α+β+γ≦1 次に、目標値T’が色再現範囲外にあるときは、前述し
たように目標値T’を色再現範囲内に移動する必要があ
る。三次元の場合について示すと、図20に示すように色
再現範囲外にある点(L5,U5,V5)から無彩色方
向に向かって色再現範囲内の点 (L0,U0,V0)と
を結ぶ直線と、前記三角錐を構成する色再現範囲で最も
外側にある頂点座標を (L1,U1,V1), (L2,
U2,V2), (L3,U3,V3)とする三角形との
交点を(Lp ,Up ,Vp )とすると、該交点は以下の
ようにして求めることができる。まず、
When α, β, and γ obtained in step (4) satisfy the following equation, they are inside (including on the surface) of the triangular pyramid, and when they do not satisfy, they are outside. α ≧ 0, β ≧ 0, γ ≧ 0, α + β + γ ≦ 1 Next, when the target value T ′ is outside the color reproduction range, it is necessary to move the target value T ′ into the color reproduction range as described above. is there. In the three-dimensional case, as shown in FIG. 20, points (L5, U5, V5) outside the color reproduction range and points (L0, U0, V0) within the color reproduction range toward the achromatic color direction are shown. The connecting straight line and the outermost vertex coordinates in the color reproduction range forming the triangular pyramid are (L1, U1, V1), (L2,
Letting (L p , U p , V p ) be the intersection with the triangle defined as (U2, V2), (L3, U3, V3), the intersection can be obtained as follows. First,

【0045】[0045]

【数5】 [Equation 5]

【0046】としてe,f,gを求め、次に、Then, e, f, and g are obtained as

【0047】[0047]

【数6】 [Equation 6]

【0048】としてδを求め、また、Δ is obtained as

【0049】[0049]

【数7】 [Equation 7]

【0050】としてηを求める。 そして、Η is obtained as And

【0051】[0051]

【数8】 [Equation 8]

【0052】尚、この場合も交点(Lp ,Up ,Vp
が頂点座標(L1 ,U1 ,V1 ),(L2 ,U2
2 ),(L3 ,U3 ,V3 )とする三角形の内側にあ
るか否かの判定は以下のように行われる。まず、前記δ
が0でなく、 で得られるα’,β’,γ’の全てが0以上であるとき
は三角形の内側(辺上を含む)にあり、そうでない場合
は外側にある。尚、δ=0の場合は直線と三角形とが平
行で交わらない場合である。
In this case also, the intersection points (L p , U p , V p )
Are vertex coordinates (L 1 , U 1 , V 1 ), (L 2 , U 2 ,
The determination as to whether or not it is inside the triangle defined by V 2 ), (L 3 , U 3 , V 3 ) is performed as follows. First, the δ
Is not 0, When all of α ', β', and γ'obtained in (1) are 0 or more, they are inside the triangle (including on the side), and otherwise they are outside. The case of δ = 0 means that the straight line and the triangle are parallel and do not intersect.

【0053】このようにして表色系と信号値との関係が
判るので、前記したように3変数a,b,cと信号値と
の関係も求められる。尚、空間 (6面体) を三角錐に分
割するには、図21に示すようにしてもよいし、図22に示
すようにしてもよい。これらの図において、空間1及び
空間2の点a〜h及びa’〜h’は互いに対応した点で
ある。
Since the relationship between the color system and the signal value is known in this manner, the relationship between the three variables a, b, c and the signal value can be obtained as described above. The space (hexahedral) may be divided into triangular pyramids as shown in FIG. 21 or as shown in FIG. In these figures, points a to h and a'to h'in the space 1 and the space 2 are points corresponding to each other.

【0054】かかる推定方式を用いて、Kの濃度を最小
とする前記各条件を満たす三次元空間について目標値T
を求めればよい。具体的には、Yが最大 (255)のとき
は、他のM,C,Kの3色について前記推定により目標
値Tを求め、同様にMが最大 (255)のときはY,C,K
について、Cが最大 (255)のときはY,M,Kについ
て、また、Kが0のときはY,M,Cについて目標値T
を求めることができる。
Using this estimation method, the target value T for the three-dimensional space satisfying the above-mentioned conditions that minimizes the density of K is obtained.
You should ask. Specifically, when Y is the maximum (255), the target value T is obtained by the above estimation for the other three colors of M, C, and K. Similarly, when M is the maximum (255), Y, C, and K
, The target value T for Y, M, K when C is maximum (255) and for Y, M, C when K is 0
Can be asked.

【0055】続いて、前記カラーパッチを用いた色再現
特性の推定方法を適用するのに好適なカラーマスキング
装置 (色分解画像修正装置) の一例を説明する。上述の
ようにして算出された目標値 (Y,M,C,Kの組み合
わせ、つまり色修正データ) がメインのルックアップテ
ーブル (MLUT) に予め格納されている。そして、入
力系がカラーCRTの場合には、R, G, Bによって決
まる基本色の色座標系に対応付けられたY,M,C,K
座標系が、与えられた色座標となる。
Next, an example of a color masking device (color separation image correction device) suitable for applying the method for estimating the color reproduction characteristic using the color patch will be described. The target value (combination of Y, M, C, K, that is, color correction data) calculated as described above is stored in advance in the main look-up table (MLUT). When the input system is a color CRT, Y, M, C, K associated with the color coordinate system of the basic color determined by R, G, B
The coordinate system becomes the given color coordinate.

【0056】この色座標は演算によって求められる。
R,G,Bの座標系、例えばB, Gの座標系とY,M,
Kの座標系との関係は図23のようになるから、入力座標
系、例えば点tはt’に、点sはs’に夫々対応付ける
処理が行われる。そのため、これらR,G,B座標系と
Y,M,C,K座標系とを対応させるため、メインルッ
クアップテーブル (MLUT) が用意され、例えば点t
を表す座標系が入力したときには、t’の座標系が参照
されるようになされる。
The color coordinates are calculated.
R, G, B coordinate system, such as B, G coordinate system and Y, M,
Since the relationship between K and the coordinate system is as shown in FIG. 23, the processing of associating the input coordinate system, for example, the point t with t ′ and the point s with s ′ is performed. Therefore, in order to make these R, G, B coordinate systems correspond to the Y, M, C, K coordinate systems, a main lookup table (MLUT) is prepared.
When the coordinate system representing is input, the coordinate system of t'is referred to.

【0057】格子点以外の色修正データは内挿によって
算出される。上述したように、本発明に係る方法におい
ては、入力座標系のデータが予め決まっている。 (a) M,C, K (Yが最大濃度) (b) Y,C,K (Mが最大濃度) (c) Y,M,K (Cが最大濃度) (d) Y,M,C (K=0) ・・・ (3) 説明を簡略化するため、 (d) について説明する。
Color correction data other than grid points are calculated by interpolation. As described above, in the method according to the present invention, the data of the input coordinate system is predetermined. (a) M, C, K (Y is maximum concentration) (b) Y, C, K (M is maximum concentration) (c) Y, M, K (C is maximum concentration) (d) Y, M, C (K = 0) (3) In order to simplify the description, (d) will be described.

【0058】この例では、図24に示すように3つの入力
画像データR,G,Bによって決まる直方体状の空間W
(その対角頂点に内挿点sがある) を含む8つの色修正
データ (Y,M,C ,K (=0)に対応した既知の算出色
修正データP1〜P8) で形成される直方体状の空間領
域Vを定める。空間領域W, VはいずれもP1を基準点
とするものである。そして、各色の、 0, 32, 64, 96, 128, 160, 1
92, 224, 255 の各点における組み合わせの色に対して、上述したよう
な色修正値を持つものとする。このとき、入力画像デー
タR,G,Bが夫々 (100, 130, 150) の値を持っていた場合、以下に示される8点で囲まれる
空間領域の頂点 (格子点) の色修正データを用いて内挿
される。
In this example, as shown in FIG. 24, a rectangular parallelepiped space W determined by three input image data R, G, B.
A rectangular parallelepiped formed from eight color correction data (known calculation color correction data P1 to P8 corresponding to Y, M, C, K (= 0)) including (interpolation point s is located at the diagonal vertex) The spatial area V of the shape is defined. Both the spatial regions W and V use P1 as a reference point. And 0, 32, 64, 96, 128, 160, 1 of each color
It is assumed that the color correction values as described above are applied to the combined colors at the points 92, 224, and 255. At this time, if the input image data R, G, and B each have a value of (100, 130, 150), the color correction data of the vertices (lattice points) of the space area surrounded by 8 points shown below are obtained. Used to interpolate.

【0059】ここに、左辺のPi (i=1〜8 ) は空間
領域Vの各頂点の座標値を示し、右辺はそのときの色修
正データKi, Ci, Mi, Yiを示す。 P1: (0, 96, 128, 128)= (K1, C1, M1, Y
1) P2: (0,128, 128, 128)= (K1, C1, M1, Y
1) P3: (0, 96, 160, 128)= (K1, C1, M1, Y
1) P4: (0,128, 160, 128)= (K1, C1, M1, Y
1) P5: (0, 96, 128, 160)= (K1, C1, M1, Y
1) P6: (0,128, 128, 160)= (K1, C1, M1, Y
1) P7: (0, 96, 160, 160)= (K1, C1, M1, Y
1) P8: (0, 96, 160, 160)= (K1, C1, M1, Y
1) 空間領域Vの各頂点Piに対する重み係数は次のように
して算出される。
Here, Pi (i = 1 to 8) on the left side indicates the coordinate value of each vertex of the spatial region V, and the right side indicates the color correction data Ki, Ci, Mi, Yi at that time. P1: (0, 96, 128, 128) = (K1, C1, M1, Y
1) P2: (0,128, 128, 128) = (K1, C1, M1, Y
1) P3: (0, 96, 160, 128) = (K1, C1, M1, Y)
1) P4: (0,128, 160, 128) = (K1, C1, M1, Y
1) P5: (0, 96, 128, 160) = (K1, C1, M1, Y)
1) P6: (0,128, 128, 160) = (K1, C1, M1, Y
1) P7: (0, 96, 160, 160) = (K1, C1, M1, Y
1) P8: (0, 96, 160, 160) = (K1, C1, M1, Y
1) The weighting coefficient for each vertex Pi of the spatial region V is calculated as follows.

【0060】本例では、求めるべき修正値の点の反対の
頂点と、内挿点sで作られる直方体の空間領域Wの体積
を、求めるべき修正値の心における重み係数Wiとする
ものである。したがって、点P8の重み係数は、P1の
座標とsの座標とを用いて、 (100,130,150)−(96,128,128)=(4,2,22) より、sとP1とで作られる直方体状の空間領域の体積
は、 4×2×22=176 となり、これが点Pの重み係数となる。
In this example, the vertex of the correction value to be obtained and the volume of the spatial area W of the rectangular parallelepiped formed by the interpolation points s are used as the weighting factor Wi at the center of the correction value to be obtained. . Therefore, the weighting factor of the point P8 is a rectangular parallelepiped spatial region formed by s and P1 from (100,130,150)-(96,128,128) = (4,2,22) using the coordinates of P1 and the coordinates of s. Is 4 × 2 × 22 = 176, which is the weighting factor of the point P.

【0061】同様にして、残りの点P1〜P7の重み係
数が算出される。 P1=8400 P2=1200 P3=560 P4=80 P7=1232 P8=178 これら重み係数の和は、立方体状の空間領域の体積Vと
同一となり、この例では、32768 (aとする) となる。
したがって、s点における修正値Ks , Cs ,Ms , Y
s は、 Ks =1/a (P1K1+P2K2+P3K3+P4K4+P5K5+P6K6 +P7K7+P8K8) Cs =1/a (P1C1+P2C2+P3C3+P4C4+P5C5+P6C6 +P7C7+P8C8) Ms =1/a (P1M1+P2M2+P3M3+P4M4+P5M5+P6M6 +P7M7+P8M8) Ms =1/a (P1Y1+P2Y2+P3Y3+P4Y4+P5Y5+P6Y6 +P7Y7+P8Y8) となる。即ち、求めたい点s、それを取り囲む8点の修
正値をKi , Ci , Mi, Yi (これは表色系の内挿値
s ’, Us ’Vs ’に対応したY,M,C,K座標系
の値である) とし、夫々の重み係数をAi とすれば、
Similarly, the weighting factors of the remaining points P1 to P7 are calculated. P1 = 8400 P2 = 1200 P3 = 560 P4 = 80 P7 = 1232 P8 = 178 The sum of these weighting factors is the same as the volume V of the cubic space area, which is 32768 (a) in this example.
Therefore, the correction values K s , C s , M s , Y at the s point
s is a K s = 1 / a (P1K1 + P2K2 + P3K3 + P4K4 + P5K5 + P6K6 + P7K7 + P8K8) C s = 1 / a (P1C1 + P2C2 + P3C3 + P4C4 + P5C5 + P6C6 + P7C7 + P8C8) M s = 1 / a (P1M1 + P2M2 + P3M3 + P4M4 + P5M5 + P6M6 + P7M7 + P8M8) M s = 1 / a (P1Y1 + P2Y2 + P3Y3 + P4Y4 + P5Y5 + P6Y6 + P7Y7 + P8Y8). That is, the correction value of the point s to be obtained and the eight points surrounding it is K i , C i , M i , Y i (this corresponds to the interpolated values L s ', U s ' V s ' of the color system). (Y, M, C, K coordinate system values) and each weighting coefficient is A i ,

【0062】 Ks = (1/ΣAi ) ΣAi is = (1/ΣAi ) ΣAi is = (1/ΣAi ) ΣAi is = (1/ΣAi ) ΣAi i で表すことができる入力座標系R,G,Bのデータが異
なるときには、その入力データに対応して夫々より算出
された重み係数Ai によって内挿される。
K s = (1 / ΣA i ) ΣA i K i C s = (1 / ΣA i ) ΣA i C i M s = (1 / ΣA i ) ΣA i M i Y s = (1 / ΣA i ) When the data of the input coordinate systems R, G, B which can be represented by ΣA i Y i are different, the data are interpolated by the weighting factors A i calculated respectively corresponding to the input data.

【0063】この入力データは既に条件式(3) に合致し
た値が入っているので補間装置では、条件式(3) により
変更することはない。尚、実際にはROMの容量などを
考慮して色修正データの数は、2の巾乗に設定される。
したがって、256 kビットのROMを使用する場合に
は、1色につき32点の色修正データ (3色全体で323
32768 点) を持たせることができる。
Since this input data already contains a value that matches the conditional expression (3), the interpolator does not change it according to the conditional expression (3). The number of color correction data is actually set to the power of 2 in consideration of the capacity of the ROM.
Therefore, when using a 256 k-bit ROM, 32 points of color correction data per color (32 3 =
You can have 32768 points).

【0064】図25は、カラーマスキング装置10の一例で
ある。色修正データ記憶手段20には、各色Y,M,K,
Kに対する色修正データが夫々のMLUT21〜24に格納
されている。25は重み係数記憶手段で、これもLUTと
して構成されている。入力画像データR,G,Bは一旦
アドレス信号形成手段40に供給されて、入力レベルに対
応した上位5ビットのアドレス信号が出力され、これが
色修正データ記憶手段20に供給される。また、これより
出力された下位3ビットの重み係数指定信号が重み係数
記憶手段25に供給される。
FIG. 25 shows an example of the color masking device 10. The color correction data storage means 20 stores each color Y, M, K,
Color correction data for K is stored in each MLUT 21-24. Reference numeral 25 is a weighting coefficient storage means, which is also configured as an LUT. The input image data R, G, B are once supplied to the address signal forming means 40, and the upper 5 bits of the address signal corresponding to the input level are outputted, which are supplied to the color correction data storage means 20. Further, the weighting coefficient designating signal of the lower 3 bits outputted from this is supplied to the weighting coefficient storage means 25.

【0065】入力画像データR,G,Bによって決まる
色座標が図21の格子点から外れているときには、その色
座標を囲むY,M,K,K座標系における4つの格子点
が、後段の色修正データ記憶手段20において指定できる
ように、5ビットのアドレス信号が入力画像データR,
G,Bによって参照されて出力されることになる。アド
レス信号形成手段40も夫々LUT (PLUT) 41〜43で
構成される。LUTとしてはバイポーラROMが好適で
ある。これらPLUT41〜43には、更にコントローラ50
から1ビットの振り分け信号が供給されるが、その詳細
については後述する。
When the color coordinates determined by the input image data R, G, B deviate from the grid points in FIG. 21, four grid points in the Y, M, K, K coordinate system surrounding the color coordinates are located at the latter stage. As can be designated in the color correction data storage means 20, the 5-bit address signal is input image data R,
It is referred to and output by G and B. The address signal forming means 40 also comprises LUTs (PLUTs) 41 to 43, respectively. A bipolar ROM is suitable as the LUT. For these PLUTs 41 to 43, a controller 50
A 1-bit distribution signal is supplied from, but the details will be described later.

【0066】入力画像データの入力レベルに対応したア
ドレス信号によって参照された色修正データ及び重み係
数を示すデータ (以下単に重み係数という) は、計8回
にわたり順次掛算累算手段30側に供給される。掛算累算
手段30は、上述したようにAi i (Bi はY,M,
K,Kの総称)を順次実行すると共に、それらの和を求
めるためのものであって、この例では掛算器31〜34と累
算器35〜38とで構成されている。
The color correction data referred to by the address signal corresponding to the input level of the input image data and the data indicating the weighting coefficient (hereinafter simply referred to as weighting coefficient) are sequentially supplied to the multiplication accumulating means 30 side eight times. It As described above, the multiplication and accumulating means 30 uses A i B i (B i is Y, M,
K) (generic name of K) is sequentially executed and the sum of them is obtained. In this example, the multipliers 31 to 34 and the accumulators 35 to 38 are used.

【0067】従って、各掛算器31〜34は、512 kビット
のROMが使用され、これらには対応する色修正データ
(8ビット) と重み係数Ai とが供給されて、Ai , B
i の乗算処理が実行され、そのうちの上位8ビットの乗
算出力は、後段の累算器 (ALU) 35〜38に供給されて
順次乗算出力が加算処理される。累算器35〜38は色ビッ
トの精度で演算されるが、累算出力 (積和出力) として
はそのうちの上位8ビットが使用される。これによっ
て、累算出力を重み係数Ai で除したと同じ出力が得ら
れることになる。
Therefore, each of the multipliers 31 to 34 uses a ROM of 512 kbit, and the color correction data corresponding to these are used.
(8 bits) and the weighting coefficient A i are supplied, and A i , B
The multiplication process of i is executed, and the multiplication output of the upper 8 bits is supplied to the accumulators (ALU) 35 to 38 in the subsequent stage, and the multiplication outputs are sequentially added. The accumulators 35 to 38 are operated with the accuracy of color bits, and the higher 8 bits are used as the cumulative calculation power (sum of products output). As a result, the same output as that obtained by dividing the cumulative calculation force by the weighting coefficient A i is obtained.

【0068】上位8ビットの累算出力は夫々ラッチ回路
46〜49によってラッチされる。ラッチパルスはコントロ
ーラ50で生成される。各部の構成を更に詳細に説明す
る。色修正データ記憶手段20として使用されるLUT21
〜24は、256 kビット容量のROMを使用した場合入力
画像データの最小レベルから最大レベルまでの間を32点
だけ抽出する。これによって一点につき32点 (したがっ
て、3色では323 =32768 点) の色修正データを格納す
ることができる。
The cumulative calculation powers of the upper 8 bits are respectively latch circuits.
Latched by 46-49. The latch pulse is generated by the controller 50. The configuration of each unit will be described in more detail. LUT 21 used as color correction data storage means 20
In the case of using a ROM having a capacity of 256 kbits, the numbers 24 to 24 extract only 32 points from the minimum level to the maximum level of the input image data. As a result, 32 points of color correction data can be stored for each point (thus, 3 3 = 32 768 points for 3 colors).

【0069】したがって、256 階調の入力レベルである
ときには、32点の配分は、例えば次に示すように、0か
ら順に「8」ずつ区切って、 0,8,16,・・・・240, 248 の合計32個となるように等分に配分し、33点目となる24
9 点以上255 点までは使用しない。若しくは、249 〜25
の点は248 として扱う。
Therefore, when the input level is 256 gradations, the distribution of 32 points is 0, 8, 16, ... 240, It is divided evenly so that a total of 32 pieces of 248 are obtained, and it becomes the 33rd point.
Do not use from 9 points to 255 points. Or 249-25
Points are treated as 248.

【0070】このような各配分点での色修正データが正
確に算出され、算出されたこれら複数の色修正データが
夫々のLUT21〜24に格納されるものである。このよう
に配分点を32点に設定すると、8ビット出力の汎用RO
Mを使用できるから記憶手段にれを安価に構成できるメ
リットがある。重み係数記憶手段4用のLUT25には、
各配分点における重み係数Ai が格納されている。い
ま、上述したように8ビットずつ配分した場合には、8
個の重み係数Ai の総計は、 8×8×8=512 となるが、上述のように出力が8ビットの市販の汎用I
Kを使用しようとするならば、理論値通りの重み係数
(最大512)を持つと素子が増えるため、この例では理論
値を略1/2 に圧縮した近似値が重み係数の実際値として
使用される。
The color correction data at each distribution point is accurately calculated, and the plurality of calculated color correction data are stored in the respective LUTs 21-24. If the distribution points are set to 32 points in this way, 8-bit output general-purpose RO
Since M can be used, there is an advantage that the storage means can be constructed inexpensively. In the LUT 25 for the weighting coefficient storage means 4,
The weighting coefficient A i at each distribution point is stored. Now, when 8 bits are allocated as described above, 8
The total of the individual weighting factors A i is 8 × 8 × 8 = 512, but as described above, the commercially available general-purpose I
If K is to be used, the weighting factor is the theoretical value
If you have (maximum 512), the number of elements will increase, so in this example, the approximate value of the theoretical value compressed to about 1/2 is used as the actual value of the weighting coefficient.

【0071】以下に示す例は、8個の重み係数の和が常
に256 となるように設定し、夫々のうちの最大の重み係
数は、255 とする。こうした場合、例えば図22におい
て、内挿点sがP1と同じ位置にあった場合、P1〜P
8の各重み係数は、 ( ) 内にその理論値で示すよう
に、 となり、重み係数の総和は、256 となる。
In the example shown below, the sum of eight weighting factors is set to always be 256, and the maximum weighting factor of each is 255. In such a case, for example, in FIG. 22, when the interpolation point s is at the same position as P1, P1 to P1
Each weighting factor of 8 is, as shown by its theoretical value in (), Therefore, the sum of the weighting factors is 256.

【0072】また、sがP1とP3との中間で、P1か
ら3 (従って、P3からは5) だけ離れた位置にあった
ときには、P1〜P8の各重み係数は次のようになる。 となり、この場合の重み係数の総和も、256 となるよう
に、各重み係数が適宜選定される。
When s is in the middle of P1 and P3 and is apart from P1 by 3 (thus, 5 from P3), the weighting factors of P1 to P8 are as follows. Therefore, each weighting factor is appropriately selected so that the sum of the weighting factors in this case is also 256.

【0073】同様にして、sがP1〜P4の面から3だ
け離れ、P1, P3 , P5 , P7 の面から1だけ離れ、
そしてP1, P2, P5, P6の面から離れていた場合
には、次のような重み係数P1〜P8となる。 となり、この場合の重み係数の総和も256 となるように
各重み係数が選定される。
Similarly, s is separated from the surfaces of P1 to P4 by 3 and separated from the surfaces of P1, P3, P5 and P7 by 1,
When the surface is away from the planes P1, P2, P5 and P6, the weighting factors P1 to P8 are as follows. Therefore, each weighting factor is selected so that the sum of the weighting factors in this case is also 256.

【0074】上述した1ビットの振り分け信号とは、各
点sを含む前後の色修正データを指定するための制御信
号である。即ち、説明の便宜上、32個の配分点 (格子
点) とそれに対応するアドレス信号との関係を図26に示
すように設定する。今、入力画像データのレベルが100
であったときには、色修正データ記憶手段20からこの入
力レベルを含む前後の色修正データ (96と104)が出力さ
れるようなアドレス信号 (12, 13) を形成する必要があ
る。
The above-mentioned 1-bit distribution signal is a control signal for designating the color correction data before and after including each point s. That is, for convenience of explanation, the relationship between 32 allocation points (lattice points) and the corresponding address signals is set as shown in FIG. Now the input image data level is 100
Then, it is necessary to form the address signal (12, 13) from which the color correction data storage means 20 outputs the color correction data (96 and 104) before and after including this input level.

【0075】そこで、振り向け信号が0のとき、小さい
方の色修正データ (96) が参照されるようなアドレス信
号 (12) が出力され、また振り分け信号が1のとき、大
きい方の色修正データ(104) が参照されるようなアドレ
ス信号 (13) が出力されるようにコントロールされる。
ただし、使用する値の最大値 (この場合は248)のとき
で、振り分け信号が0のときには、それ自身の値の色修
正データを選定し、振り分け信号が1のときには小さい
方の1 (この場合240)を制定する。
Therefore, when the distribution signal is 0, the address signal (12) is output so that the smaller color correction data (96) is referred to. When the distribution signal is 1, the larger color correction data (96) is output. It is controlled so that the address signal (13) for referring to (104) is output.
However, when the distribution signal is 0 at the maximum value to be used (248 in this case), the color correction data of its own value is selected, and when the distribution signal is 1, the smaller one (in this case, 1) 240) is established.

【0076】また、別の補間方法について説明する。い
ま、図27に示すように入力信号値R,G,Bが前記LU
TのR,G,Bの8つの格子点を頂点A〜Hとする立方
体で囲まれているとする。その場合、前記目標値を求め
る場合と同様、該立方体を1点鎖線で示すように6つの
三角錐に分割できる。そして、入力信号値R,G,Bが
どの三角錐に含まれているかが同様にして求めることが
できる。いま、入力信号値R,G,Bの値を(5,1,
2)であるとすると、この点は頂点A,B,C,Gによ
って形成される三角錐に含まれることが判る。この三角
錐に対応するCIELUV表色系の三角錐を図28に示
す。A’〜G’は夫々A〜Gに対応するものとすれば、
補間点P’も三角錐T’内に存在する。
Another interpolation method will be described. Now, as shown in FIG. 27, the input signal values R, G, B are the same as the LU.
It is assumed that the eight grid points of R, G, and B of T are surrounded by a cube having vertices A to H. In that case, as in the case of obtaining the target value, the cube can be divided into six triangular pyramids as shown by the one-dot chain line. Then, which triangular pyramid the input signal values R, G, B are included in can be similarly obtained. Now, the values of the input signal values R, G, B are (5, 1,
2), this point is included in the triangular pyramid formed by the vertices A, B, C, and G. FIG. 28 shows a triangular pyramid of the CIE LUV color system corresponding to this triangular pyramid. If A ′ to G ′ correspond to A to G, respectively,
The interpolation point P ′ also exists within the triangular pyramid T ′.

【0077】この三角錐Tが決定されると、次に同図A
に示すように、補間点Pと頂点A,B,C,Gとが結ば
れて計4個の新たな三角錐が形成され、夫々の体積VBC
GP,VACGP, VABGP, VABCPが求められる。これらの体
積VBCGP, VACGP, VABGP,VABCPと、同図Bの表色系
の頂点A’, B’, C’, G’とから補間点P’が次式
によって算出される。なお、VABGPは三角錐Tの体積で
ある。
When this triangular pyramid T is determined, next, in FIG.
As shown in, the interpolation point P and the vertices A, B, C, and G are connected to form a total of four new triangular pyramids, and each volume VBC
GP, VACGP, VABGP, VABCP are required. From these volumes VBCGP, VACGP, VABGP, VABCP and the vertices A ', B', C ', G'of the color system of FIG. Note that VABGP is the volume of the triangular pyramid T.

【0078】 P’=1/VABGP(VBCGP・A’+VACGP・B’+VABGP・C’+VABCP・G ’) このようにして、入力される任意の信号値R,G,Bに
対してLUTからの検索と、補間演算とにより対応する
CIELUV表色系の値が求められる。尚、本発明にお
いて、与えられた色座標は入力装置に依存する。上述し
たように、入力装置がディスプレーならば、そのR,
G,Bの値から演算によって色座標が求められる。入力
装置が印刷系ならば、そのY,M,C,K4色から演算
で求められる。
P ′ = 1 / VABGP (VBCGP · A ′ + VACGP · B ′ + VABGP · C ′ + VABCP · G ′) In this way, for the arbitrary signal values R, G, B input from the LUT The value of the corresponding CIE LUV color system is obtained by the search and the interpolation calculation. In the present invention, the given color coordinates depend on the input device. As described above, if the input device is a display, its R,
Color coordinates are obtained from the values of G and B by calculation. If the input device is a printing system, it can be calculated from the Y, M, C and K four colors.

【0079】また、入力装置がスキャナ系であるなら
ば、そのR, G, Bの値から色座標が求められることに
なる。
If the input device is a scanner system, color coordinates can be obtained from the R, G, B values.

【0080】[0080]

【発明の効果】以上説明してきたように本発明によれ
ば、与えられた色座標を示すY,M,C,Kの組み合わ
せのうちから、Kの濃度が最小となるものを採る構成と
したため、特定点に対応するカラーパッチの数を小さく
できる一方で、グレー色が多くの色を混合して再現され
るため、階調性の悪い電子写真プリンタやインクジェッ
トプリンタにおける階調性の劣化が目立ちにくくなり、
画質を向上できるものである。
As described above, according to the present invention, the combination of Y, M, C, and K indicating the given color coordinates has the minimum K density. While the number of color patches corresponding to a specific point can be reduced, many gray colors are reproduced by mixing many colors, so deterioration of the gradation is noticeable in electrophotographic printers and inkjet printers with poor gradation. Becomes difficult,
The image quality can be improved.

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

【図1】本発明の実施例に係るカラーパッチの図FIG. 1 is a diagram of a color patch according to an embodiment of the present invention.

【図2】図1で示したカラーパッチ測色値をCIELU
V表色系に写像したときの図
FIG. 2 shows CIELU as the color patch colorimetric values shown in FIG.
Figure when mapped to V color system

【図3】本発明の実施例に係るカラーパッチの別の例を
示す図
FIG. 3 is a diagram showing another example of the color patch according to the embodiment of the invention.

【図4】図3で示したカラーパッチ測色値をCIELU
V表色系に写像したときの図
FIG. 4 shows CIELU as the color patch colorimetric values shown in FIG.
Figure when mapped to V color system

【図5】曲線近似の説明図FIG. 5 is an explanatory diagram of curve approximation.

【図6】図5の曲線近似で得られる目標値拡張の説明図6 is an explanatory diagram of target value expansion obtained by curve approximation in FIG.

【図7】目標値推定の原理を示す図FIG. 7 is a diagram showing the principle of target value estimation.

【図8】重み平均による計算を説明するための図FIG. 8 is a diagram for explaining calculation by weighted average.

【図9】内挿処理をしたY,M座標系の図FIG. 9 is a diagram of a Y and M coordinate system after interpolation processing.

【図10】同じくCIELUV表色系の図FIG. 10 is a diagram of the CIE LUV color system as well.

【図11】YM座標系空間を三角形に分割した図FIG. 11 is a diagram in which the YM coordinate system space is divided into triangles.

【図12】CIELUV座標系空間を三角形に分割した図FIG. 12 is a diagram in which the CIE LUV coordinate system space is divided into triangles.

【図13】目標値T’に対応するY,Mの組み合わせTを
示す図
FIG. 13 is a diagram showing a combination T of Y and M corresponding to a target value T ′.

【図14】目標値T’をCIELUV表色系で示した図FIG. 14 is a diagram showing a target value T ′ in the CIE LUV color system.

【図15】目標値T’が色再現範囲外に存在する場合を示
した図
FIG. 15 is a diagram showing a case where a target value T ′ exists outside the color reproduction range.

【図16】目標値T’を色再現範囲内に移動させる説明図FIG. 16 is an explanatory diagram of moving a target value T ′ into a color reproduction range.

【図17】新目標値T’に対応するY,Mの組み合わせT
の算出
FIG. 17 is a combination T of Y and M corresponding to a new target value T ′.
Calculation of

【図18】新目標値T’のCIELUV表色系の位置を示
す図
FIG. 18 is a diagram showing the position of the CIE LUV color system of the new target value T ′.

【図19】目標値T,T’をY,M,C座標系とCIEL
UV表色系とで示した図
[FIG. 19] Target values T and T ′ are set to Y, M, and C coordinate systems and CIEL
Figure shown with UV color system

【図20】目標値T’がY,M,C座標系の色再現範囲外
に存在する場合を示した図
FIG. 20 is a diagram showing a case where the target value T ′ exists outside the color reproduction range of the Y, M, and C coordinate systems.

【図21】空間を三角形に分割する方法の一例を示す図FIG. 21 is a diagram showing an example of a method of dividing a space into triangles.

【図22】同じく空間を三角形に分割する方法の別の例を
示す図
FIG. 22 is a diagram showing another example of a method of dividing a space into triangles.

【図23】B,G座標系とY,M,K座標系との関係を示
す図
FIG. 23 is a diagram showing a relationship between a B, G coordinate system and a Y, M, K coordinate system.

【図24】内挿処理の説明図FIG. 24 is an explanatory diagram of interpolation processing.

【図25】本発明を適用できるカラーマスキング装置の系
統図
FIG. 25 is a system diagram of a color masking device to which the present invention can be applied.

【図26】同上装置に使用されるアドレスと色修正データ
との関係を示す図
FIG. 26 is a diagram showing a relationship between an address used in the same device and color correction data.

【図27】三次元補間演算を説明するための図FIG. 27 is a diagram for explaining a three-dimensional interpolation calculation.

【図28】同じく補間演算を説明するための図[FIG. 28] Similarly, a diagram for explaining interpolation calculation.

【図29】100 %UCR法の説明図FIG. 29: Illustration of 100% UCR method

【図30】Y,M,C,Kによるカラーパッチの説明図FIG. 30 is an explanatory diagram of color patches using Y, M, C, and K.

【図31】本発明の説明に供するカラー画像形成装置の系
統図
FIG. 31 is a system diagram of a color image forming apparatus used for explaining the present invention.

【符号の説明】[Explanation of symbols]

10 カラーマスキング装置 12 カラープリンタ 13 記録媒体 10 Color masking device 12 color printer 13 recording medium

Claims (2)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】 Y (イエロー) 、M (マゼンタ) ,C
(シアン) 及びK (黒) の4色を組み合わせて色再現を
行う色再現方法であって、 与えられた色座標を示すY、M、C、Kの組み合わせの
うち、Y、M、Cの少なくとも一つが最大となる場合に
Kを含めて目標値を求める処理、およびKを0とした場
合にY、M、Cのみで目標値を求める処理を用いること
によって、Kの濃度が最小となるY,M,C,Kの組み
合わせに対応するY、M、C、Kの濃度値を推定して色
再現を行うことを特徴とする色再現方法。
1. Y (yellow), M (magenta), C
(Cyan) and K (black) are combined for color reproduction
A color reproduction method to be performed, which is a combination of Y, M, C, and K indicating a given color coordinate.
Of these, if at least one of Y, M, and C is maximum
Processing to obtain the target value including K, and when K is 0
In this case, use the process to obtain the target value using only Y, M, and C
The combination of Y, M, C, and K that minimizes the concentration of K
Estimate the density value of Y, M, C, K corresponding to the matching
A color reproduction method characterized by performing reproduction.
【請求項2】 Y (イエロー) 、M (マゼンタ) ,C
(シアン) 及びK (黒) の4色を組み合わせて色再現を
行う色再現方法であって、 与えられた色座標を示すY、M、C、Kの組み合わせの
うち、Y、M、Cの少なくとも一つが最大となる場合に
Kを含めて目標値を求める処理、およびKを0とした場
合にY、M、Cのみで目標値を求める処理を用いること
によって、Kの濃度が最小となるY,M,C,Kの組み
合わせを、離散的に与えられた色座標の各特定点につい
て夫々求めると共に、 これら求めた組み合わせによるY、M、C、Kのカラー
パッチを作成し、 作成されたカラーパッチの測色値から、上記特定点以外
の色座標に対応するY、M、C、Kの濃度値を推定して
色再現を行うことを特徴とする色再現方法。
2. Y (yellow), M (magenta), C
(Cyan) and K (black) are combined for color reproduction
A color reproduction method to be performed, which is a combination of Y, M, C, and K indicating a given color coordinate.
Of these, if at least one of Y, M, and C is maximum
Processing to obtain the target value including K, and when K is 0
In this case, use the process to obtain the target value using only Y, M, and C
The combination of Y, M, C, and K that minimizes the concentration of K
For each specific point of discretely given color coordinates,
And the color of Y, M, C, K according to the combination obtained.
A color characterized by performing a color reproduction by creating a patch, estimating the density values of Y, M, C, and K corresponding to the color coordinates other than the specific point from the colorimetric value of the created color patch. How to reproduce.
JP26671892A 1992-09-09 1992-09-09 Color reproduction device and color reproduction method Expired - Lifetime JP3362281B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP26671892A JP3362281B2 (en) 1992-09-09 1992-09-09 Color reproduction device and color reproduction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP26671892A JP3362281B2 (en) 1992-09-09 1992-09-09 Color reproduction device and color reproduction method

Publications (2)

Publication Number Publication Date
JPH0698161A JPH0698161A (en) 1994-04-08
JP3362281B2 true JP3362281B2 (en) 2003-01-07

Family

ID=17434721

Family Applications (1)

Application Number Title Priority Date Filing Date
JP26671892A Expired - Lifetime JP3362281B2 (en) 1992-09-09 1992-09-09 Color reproduction device and color reproduction method

Country Status (1)

Country Link
JP (1) JP3362281B2 (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3273204B2 (en) * 1993-02-15 2002-04-08 コニカ株式会社 Color separation image processing method and apparatus for four-color printer
JP3912486B2 (en) 2001-10-11 2007-05-09 富士ゼロックス株式会社 Color processing method, storage medium, color processing apparatus, color conversion apparatus, and image forming apparatus
JP3867559B2 (en) 2001-11-22 2007-01-10 ブラザー工業株式会社 Look-up table creation method, look-up table, color conversion method, and color conversion device
US7593133B2 (en) 2004-02-16 2009-09-22 Brother Kogyo Kabushiki Kaisha Inkjet printing device
JP2007329737A (en) * 2006-06-08 2007-12-20 Fujifilm Corp Device and program for creating color conversion definition

Also Published As

Publication number Publication date
JPH0698161A (en) 1994-04-08

Similar Documents

Publication Publication Date Title
US5065234A (en) Method and apparatus for correcting the color of a printed image
US5581376A (en) System for correcting color images using tetrahedral interpolation over a hexagonal lattice
US5121196A (en) Color processing method and apparatus with a color patch
US5162899A (en) Color data correction apparatus ultilizing neural network
US5087126A (en) Method of estimating colors for color image correction
US20030072016A1 (en) Color conversion with hue straightening using multiple look-up tables and interpolation
JPH10136216A (en) Color picture processor
EP1221812B1 (en) Fast interpolation of large color lookup tables
US7304767B2 (en) Color processing method, storage medium, color processing apparatus, color transforming apparatus, and image forming apparatus
JP2952489B2 (en) Image recording device
JP3362281B2 (en) Color reproduction device and color reproduction method
JPH0575848A (en) Interpolation on system and color correcting method
US6031642A (en) Tetrahedral and pruned tetrahedral interpolation
JP2554366B2 (en) Method of estimating color reproduction characteristics using color patches
US5966474A (en) Non-symmetric radial and non-symmetric pruned radial interpolation
JP2575134B2 (en) Color estimation method by color patch
JPH08116460A (en) Color masking parameter determining device
JP2898030B2 (en) Color estimation method
JP3611890B2 (en) Color conversion device
JP2528463B2 (en) Color separation image correction method and apparatus
JPH0946542A (en) Interpolation arithmetic method and data converter
JP2000022975A (en) Device and method for converting data
US6833838B2 (en) Color adjustment method, color adjustment method implementation program and program storage medium
JP4049763B2 (en) Color conversion table, color conversion device, color signal generation device, color conversion table creation method, color conversion method, and color gamut determination method
JPH11220630A (en) Method for generating three-dimensional look up table, image processor executing the same and digital color printer provided with the same

Legal Events

Date Code Title Description
S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20081025

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20081025

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20091025

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101025

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111025

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121025

Year of fee payment: 10

EXPY Cancellation because of completion of term