# Evaluating the Performance of Polynomial Regression Method with Different Parameters during Color Characterization

Authors: Bangyong Sun; Han Liu; Shisheng Zhou; Wenli Li
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

The polynomial regression method is employed to calculate the relationship of device color space and CIE color space for color characterization, and the performance of different expressions with specific parameters is evaluated. Firstly, the polynomial equation for color conversion is established and the computation of polynomial coefficients is analysed. And then different forms of polynomial equations are used to calculate the RGB and CMYK’s CIE color values, while the corresponding color errors are compared. At last, an optimal polynomial expression is obtained by analysing several related parameters during color conversion, including polynomial numbers, the degree of polynomial terms, the selection of CIE visual spaces, and the linearization.


---

## Body

## 1. Introduction

As color electronics often have different imaging characteristics, such as the imaging mechanism, color space, apparatus capability, and material peculiarity [1], the color images always look different in some detail when they are output by different devices. For example, the same image displayed on monitor looks brighter and more colorful than that printed out on paper by printer. Even the same image displayed on two different monitors sometimes produces different visual effects. Thus, for color signal processing system with several color devices, in order to maintain the color consistence of color images, the precision of color signal transmission between different devices must be high enough. Now, for most of the color signal processing systems, the device-connection space is often used, such as the CIE color spaces CIEXYZ and CIELAB [2]. If the CIELAB space, for example, is chosen as the standard connection space, the color transmission process can be divided into two parts, device-to-CIELAB and CIELAB-to-device. Therefore, the color signal processing precision is highly depending on the color conversion algorithms between device colors and CIELAB colors.There are many models which can be used for converting color signals, such as Neugebauer model, neural network, interpolation method, and polynomial regression. The regression model is widely used in color signal processing systems [3], since it can produce the high accuracy by using less sample data and also it can be used in both the device-to-CIELAB and the CIELAB-to-device directions. However, there are still some problems unresolved for this model during processing color signals; for example,(1)
the calculation precision of this model is highly dependent on the number and the degree of polynomial terms [4], so it is important to obtain the optimal polynomial expressions for specific signal processing process;(2)
the selection of device-connection space, such as CIEXYZ and CIELAB spaces, may have some influence on the signal processing precision [5, 6], so it still needs to be analyzed and tested for polynomial regression models;(3)
in some cases, the RGB signals are linearized before processing with CIE colors, but in other cases linearization is not added. Hence, for both the RGB and the CMYK signals, the effect of linearization processing should be analyzed and tested [7, 8], which may reveal whether or not it should be added for specific color devices and polynomials.In this paper, these issues above are analyzed and tested in corresponding experiments. The different polynomials expressions, the different device-connection color spaces, and influence of linearization for signal processing are all tested on RGB and CMYK devices. At last, for the specific RGB and CMYK color signal processing systems, the optimal parameters are obtained with detailed analysis.


## 2. Polynomial Regression Model for Color Signal Processing

Polynomial regression is a form of linear regression in which the relationship between the independent variablex and the dependent variable y is modeled as an nth degree polynomial. Meanwhile, polynomial regression fits a nonlinear relationship between the value of x and the corresponding conditional mean of y values denoted by E(y∣x) and has been used to describe nonlinear phenomena, such as the growth rate of tissues [9], the distribution of carbon isotopes in lake sediments [10], and the progression of disease epidemics [11].Within the polynomial regression model, if the dependent variabley and multiple independent variables x1,x1,…,xM have the linear relationship and there are N groups of sample data:
(1)(xt1,xt2,…,xtM;yt)t=1,2,…,N,
then the relationship between them can be described as [12]
(2)y1=β0+β1x11+β2x12+⋯+βMx1M+ε1,y2=β0+β1x21+β2x22+⋯+βMx2M+ε2,⋮,yN=β0+β1xN1+β2xN2+⋯+βMxNM+εN,
where β0,β1,…βM are the coefficients to be determined and ε1,ε2,…εN are independent random variables.The system of expressions above can be represented using the matrix(3)Y=Xβ+ε,
where
(4)Y=[y1y2⋮yN],X=[1x11x12⋯x1M1x21x22⋯x2M⋯1xN1xN2⋯xNM]β=[β1β2⋮βN],ε=[ε1ε2⋮εN].Ifb0,b1,…,bM are the estimated values by least squares methods for parameter β, then the regression equation is
(5)y^=b0+b1x1+⋯+bMxM.From the principle of least squares [13–15], the coefficients of b0,b1,…,bM should obtain the minimal residuals square sum for all the measured value yt and regression value yt^:
(6)Q=∑i=1N(yt-yt^)2.By using least squares method, the coefficientsb can be resolved as follows:
(7)b=(XTX)-1XTY.In addition, the polynomial regression method can also be used to describe nonlinear problems, in which the dependent variabley is modeled as an nth degree polynomial of independent variables, so this model can be rightly used in color signal processing systems. Taking the RGB monitor as an example, with the CIEXYZ device-connection space, the relationship between RGB and CIEXYZ can be expressed as
(8)X=∑i=0n∑j=0n∑k=0nβXRiGjBk,Y=∑i=0n∑j=0n∑k=0nβYRiGjBk,Z=∑i=0n∑j=0n∑k=0nβZRiGjBk,
where β is the polynomial coefficients, n is the degree of polynomial, and i+j+k≤n, and the expression above can also be represented using matrix
(9)[XYZ]T=B3×L×ρL,
where the B3×L is the coefficients matrix and ρL is the matrix of polynomials, while L represents the number of polynomials.Thus, whenn=1, L=4, the first-degree polynomial matrix ρ4 is shown as follows:
(10)ρ4=[RGB1]T,
when n=2, L=10, the second-degree polynomial matrix ρ10 should be
(11)ρ10=[RBRGGBR2G2B2ρ4]T,
when n=3, L=20, the third-degree polynomial matrix ρ20 is
(12)ρ20=[RGBR3G3B3R2GG2BB2RG2RB2GR2Bρ10]T,
and when n=4, L=35, the fourth-degree polynomial matrix ρ35 is
(13)ρ35=[R4G4B4R2GBG2RBB2RGR3GR3BG3RG3BB3RB3GR2G2R2B2G2B2ρ20]T.
There are also some other forms of polynomials used in color signal processing, and the polynomial matrixes are shown as follows:
(14)ρ3=[RGB]T,ρ6=[RGRBGBρ3]T,ρ8=[RGBρ61]T,ρ14=[R3G3B3RGBρ10]T.When the coefficientsn and L are defined in color signal processing, with the sample color data which consist of device colors and CIEXYZ colors, the coefficient matrix B3×L can be resolved using least square method.


## 3. Study on the Key Parameters during Color Signal Processing

To determine the key parameters within the polynomial regression model, a color signal processing system with several color devices is introduced in the experiment. As the additive primary color isRGB and subtractive primary color is CMYK, an RGB monitor and CMYK printer are chosen as the typical testing color devices. Within the color signal processing system, the device-connection space is CIEXYZ or CIELAB, so the color processing is mainly based on four color spaces.For the purpose of obtaining the relationship between the device colors and the device-connection colors, the training sample data should be gathered in advance. For the IBM monitor, the three primary channels Red, Green, and Blue are all divided evenly into 9 parts, and each value ofR, G, and B colors ranges within [0326496128160192224255]. When all these 93=729 patches are displayed on monitor, the corresponding CIEXYZ and CIELAB colors are measured with Spectrophotometer X-Rite DTP94. These RGB and corresponding CIEXYZ or CIELAB colors form the training sample data. To verify the accuracy of polynomial regression model, the testing sample data should also be collected. Similar to the training sample data, the testing data consists of 73=349 color patches with the single channel ranging within [164880112144176208240].For the CMYK Epson 9880 printer, the single channel is divided into 11 parts with the interval 10, so every color channel ranges within [0102030405060708090100]. Because the subtractive primary colors are Cyan, Magenta, and Yellow (CMY) and the color of Black can be seen as the replacement of a certain amount of CMY, in experiment the device color CMYK is treated as CMY. Thus, when all the 113=1331 CMY color patches are printed out, the corresponding CIEXYZ and CIELAB colors are measured with Spectrophotometer X-Rite 528, and all these CMY and corresponding CIE colors form the training sample data. In addition, the testing data consists of 63=216 patches with the single channel ranging within [52545658595].In general, the regression errors for the training sample reduce data as the number of polynomial terms or the degree of polynomials increases. However, for the color data of entire range, the regression errors will increase when the number of polynomial terms exceeds a certain value. Therefore, it is highly important and necessary to find the optimal polynomial expressions producing least errors, especially determined by the number of polynomial terms or the degree of polynomials. In experiment, to evaluate the different polynomial expressions employed in color signal processing, the error computation is defined by the color difference CIE76 formula [16]:
(15)ΔE=(L2*-L1*)2+(a2*-a1*)2+(b2*-b1*)2,
where (L1*,a1*,b1*) is the regression color and (L2*,a2*,b2*) is the measured color.


### 3.1. Number and the Degree of Polynomial Terms

To find the appropriate polynomial expressions, theRGB monitor is tested and the signal processing errors with different polynomial expressions are compared. Firstly, the training sample data is used to obtain the polynomial coefficients between RGB and CIELAB signals; then for all the RGB colors, the CIELAB values can be simulated by using the obtained coefficients; secondly, in order to analyze the regression precision, the measured CIELAB color values from the training data are used to compute the errors of different polynomial expressions; at last the errors are represented as color differences between measured and simulated CIELAB values as below.From the above result, it can be seen that the regression precision obviously improves as the number of polynomial terms increases, but when the number of terms reaches a certain value the mean error becomes small enough. For example, for the first-degree polynomialρ4, the average error is 10.9106ΔE which exceeds the reproduction error threshold [17, 18], and the maximal error is 37.5866ΔE which is a visually unacceptable error. Additionally its standard deviation and variance are 5.4642ΔE and 29.8576ΔE, respectively, which indicate that the distribution of errors is unsatisfactory.In general, the regression precision can be evaluated mainly from the average error for different polynomials shown in Figure1. The regression errors for polynomials ρ4, ρ6, ρ8, and ρ10 are all exceeding 5ΔE units, while for the other polynomials ρ14, ρ20, and ρ35 their color differences are all acceptable. The figure shows that polynomial ρ35 is the most accurate, but its precision is very close to polynomial ρ20. In addition, too many terms of the polynomial may increase the difficulty of coefficients-solving process [14], so the polynomial ρ20 should be most suitable for RGB color signal processing.Figure 1
Different polynomials’ average error of theRGB training samples.Using the regression coefficients from the training data, the relationship betweenRGB and CIELAB signals can be described. To verify the precision of different polynomials for the whole range of RGB signals, the testing sample data should be used. For the testing data which is outside the range of training data, the RGB’s corresponding simulation CIELAB values can be computed using the relationship obtained above, and the errors can also be calculated by comparing the measured values. For the different polynomials and the testing data, the errors are shown below.From Table2, for the RGB and CIELAB color signals, the signal processing precision of different polynomials for testing data is similar to the training data shown in Table 1, and the acceptable forms of polynomials are ρ14, ρ20, and ρ35, respectively.Table 1
The color differences of RGB training sample data for different polynomials.
Mean ΔE
Min ΔE
Max ΔE
Std. ΔE/var. ΔE
ρ
4
10.9106
1.1670
37.5866
5.4642/29.8576
ρ
6
9.4279
1.1217
35.0208
4.9875/24.8748
ρ
8
8.7390
1.3663
30.4869
4.3973/19.3359
ρ
10
6.5693
0.6068
23.8992
3.1985/10.2302
ρ
14
3.9042
0.1371
20.8066
2.9348/8.6133
ρ
20
2.5073
0.1744
16.2386
1.9063/3.6338
ρ
35
1.8134
0.0788
8.5782
1.2957/1.6787Table 2
The color difference of RGB testing data for different polynomials.
Mean ΔE
Min ΔE
Max ΔE
Std. ΔE/var. ΔE
ρ
4
9.4627
1.2718
29.3782
4.8767/23.7826
ρ
6
8.2973
0.7076
28.1559
4.0859/16.6949
ρ
8
7.7875
1.3483
21.7351
3.5269/12.4392
ρ
10
5.5913
0.2520
17.0864
3.1055/9.6441
ρ
14
3.8443
0.1035
19.2380
3.1141/9.6974
ρ
20
2.6939
0.0751
16.2386
2.3415/5.4824
ρ
35
2.3701
0.2423
13.1631
1.8308/3.3518For the purpose of testing the different polynomials’ performance for CMYK signals, the training sample data of EPSON9880 printer are used to obtain the relationship between CMYK and CIELAB. Because the errors are too large for polynomials with few terms such asρ4 and ρ6, polynomials ρ8, ρ10, ρ14, ρ20, and ρ35 are only tested for CMYK signals. With the obtained polynomial coefficients solved with training data, for the 216 testing patches of testing sample data, the color differences are shown in Table 3.Table 3
The color difference of CMYK testing data for different polynomials.
Mean ΔE
Min ΔE
Max ΔE
Std. ΔE/var. ΔE
ρ
8
2.5541
0.5667
5.7650
1.1044/1.2197
ρ
10
1.8277
0.1946
4.4006
0.7523/0.5659
ρ
14
1.5223
0.4062
3.1525
0.5682/0.3229
ρ
20
1.2485
0.3233
2.6879
0.4690/0.2200
ρ
35
1.1592
0.1275
2.5541
0.4179/0.1746It can be seen that, for the CMYK devices, most of the polynomials perform well with the average error below3ΔE. This is mainly because the color-rendering properties of CMYK printers surpass those of the RGB monitors, such as the regularity of color gamut and color consistence [1, 19]. On the whole the fourth-degree polynomial ρ35 with the largest number of terms has the smallest color difference and its error distribution is also ideal. The performance of the third-degree polynomial ρ20 is close to the ρ35, which indicates that the preferred polynomials for CMYK and CIELAB signal processing should be the ρ35 or ρ20.


### 3.2. Selection of CIE Color Spaces during Color Signal Processing

During the signal processing for different color devices, the use of two CIE color spaces, CIEXYZ and CIELAB spaces, often brings in different color errors. Hence, it is important to find the appropriate CIE color space for the specified color devices and polynomials.In the experiment, the CIEXYZ and CIELAB color spaces are tested, respectively, onRGB monitors and on CMYK printers to compare their color errors. For the IBM monitor, when the CIEXYZ is selected as device-connection space, the errors of different polynomials are shown in Table 4. Corresponding to the color errors listed in Table 2 where the CIELAB is the device-connection space, the mean color differences corresponding to these two CIE spaces are compared in Figure 2.Table 4
The color difference of RGB monitor for different polynomials using CIEXYZ space.
Mean ΔE
Min ΔE
Max ΔE
Std. ΔE/var. ΔE
ρ
8
13.2758
1.3978
93.3629
9.8920/97.8520
ρ
10
4.9488
0.2077
19.3735
3.1563/9.9620
ρ
14
2.6386
0.2299
13.6674
2.0601/4.2441
ρ
20
2.5788
0.1656
17.8186
2.0549/4.2226
ρ
35
2.2450
0.3136
10.8137
1.661/2.7585Figure 2
Comparison of theRGB monitor’s color errors between CIEXYZ and CIELAB spaces.As the figure shows, for theRGB monitors, the polynomials perform somewhat differently with the CIEXYZ and CIELAB color spaces. For the low-degree polynomials, the color error of signal processing with CIEXYZ space is greater than with CIELAB space, and along the increasing of polynomial terms, the influence of device-connection color space becomes very little. So when the polynomials of less than 10 terms are used in RGB signal processing, the CIELAB color space is recommended as the device-connection space, while for the case when the polynomial terms are between 10 and 20, the CIEXYZ space is suggested.To test the influence of CIE color space for CMYK devices, the errors of CIEXYZ and CIELAB color spaces are also compared on Epson9889 printer. The CMYK signal processing errors with CIELAB space have been listed in Table3; for the polynomials of 8 to 35 terms, the color errors with CIEXYZ space are recorded in Table 5.Table 5
The color differences of CMYK printer for different polynomials using CIEXYZ space.
Mean ΔE
Min ΔE
Max ΔE
Std. ΔE/var. ΔE
ρ
8
11.9339
0.9828
37.2081
7.6683/58.8035
ρ
10
4.6879
0.5701
19.8019
3.5922/12.9040
ρ
14
4.3475
0.5852
20.1294
3.5023/12.2659
ρ
20
1.5809
0.3099
4.0116
0.6906/0.4769
ρ
35
1.2455
0.2666
3.0209
0.4968/0.2468It can be seen that, for the second-degree polynomialsρ10 and ρ14, third-degree polynomial ρ20, and fourth-degree polynomial ρ35, all the color errors are acceptable. Taking account of the precision and computing efficiency, the third-degree polynomial ρ20 is the most suitable model for CMYK devices using CIEXYZ color space. In addition, similar to the comparison of RGB signal processing with different device-connection spaces, the CMYK signal processing using CIELAB and CIEXYZ spaces is compared in Figure 3.Figure 3
Comparison of the CMYK printer’s color errors between CIEXYZ and CIELAB spaces.It can be seen that, within the CMYK signal processing based on polynomial regression method, the precision of color conversion using CIELAB color space is higher than that of using CIEXYZ space for a majority of the second, third, and fourth degree of polynomials. Therefore, in CMYK signal processing, the CIELAB color space is preferred as the device-connection space for the polynomial regression model.


### 3.3. Linearization during Color Signal Processing

For theRGB devices, the linearization is often applied to calibrate the color device’s gray balance [20], in which the RGB signals are firstly converted into lightness signals before color conversion process. In the experiment, the linearization is described as follows:
(16)Rl=fR(R),Gl=fG(G),Bl=fB(B),
where Rl, Gl, Bl, respectively, are the lightness signals, respectively, and fR/G/B are the linearization functions which are described as follows:
(17)fC(C)=α0+α1C+α2C2+α3C3,
where C stands for one of the colors R, G, or B and αi(i=0,1,2,3) are the coefficients.Within theRGB color signals processing, the device colors RGB are firstly linearized into the lightness signals (RGB)l, then the relationship between (RGB)l and CIELAB is obtained by polynomial regression model, and at last the color error is calculated by using the measured colors within the testing data. For the IBM monitor in the experiment, the color differences of different linearized polynomials are shown in Table 6.Table 6
The RGB color difference of different polynomials with linearization applied.
Mean ΔE
Min ΔE
Max ΔE
Std. ΔE/var. ΔE
ρ
8
8.0337
0.8814
18.4480
3.0490/9.2965
ρ
10
4.5930
0.4682
17.8630
2.6520/7.0334
ρ
14
3.9435
0.3297
18.0798
2.6400/6.9698
ρ
20
2.7392
0.4485
14.3771
1.8530/3.4335
ρ
35
2.3601
0.2703
9.6480
1.3645/1.8619To test the influence of the linearization on the color signal processing precision, the two groups of color errors in Tables2 and 6 are compared. As shown in Figure 4, the errors are very close for the two processes, so a conclusion is reached that for the RGB devices the linearization has little impact on the signal processing precision especially for the polynomials of degree greater than or equal to three.Figure 4
TheRGB monitor’s color error comparison with linearization.Similarly, to test the linearization for CMYK signal processing, the color errors of testing sample of EPSON printer are recorded in Table7, and the comparison with the signal processing without linearization in Table 3 is described in Figure 5.Table 7
The CMYK color difference of different polynomials with linearization applied.
Mean ΔE
Min ΔE
Max ΔE
Std. ΔE/var. ΔE
ρ
8
4.2971
0.5850
10.9636
2.1104/4.4537
ρ
10
2.1092
0.1526
4.5244
0.8569/0.7343
ρ
14
1.8885
0.6387
4.2186
0.7901/0.6242
ρ
20
1.5002
0.2708
3.3267
0.5553/0.3084
ρ
35
1.1958
0.1047
3.0043
0.5777/0.3338Figure 5
The CMYK printer’s color error comparison with linearization.It can be seen that, for the CMYK printers, the precision of signal processing with linearization is lower than that without linearization for most polynomial models. In some cases the linearization process does not improve the CMYK signal processing precision, so it is not advisable for CMYK devices calibration.


## 3.1. Number and the Degree of Polynomial Terms

To find the appropriate polynomial expressions, theRGB monitor is tested and the signal processing errors with different polynomial expressions are compared. Firstly, the training sample data is used to obtain the polynomial coefficients between RGB and CIELAB signals; then for all the RGB colors, the CIELAB values can be simulated by using the obtained coefficients; secondly, in order to analyze the regression precision, the measured CIELAB color values from the training data are used to compute the errors of different polynomial expressions; at last the errors are represented as color differences between measured and simulated CIELAB values as below.From the above result, it can be seen that the regression precision obviously improves as the number of polynomial terms increases, but when the number of terms reaches a certain value the mean error becomes small enough. For example, for the first-degree polynomialρ4, the average error is 10.9106ΔE which exceeds the reproduction error threshold [17, 18], and the maximal error is 37.5866ΔE which is a visually unacceptable error. Additionally its standard deviation and variance are 5.4642ΔE and 29.8576ΔE, respectively, which indicate that the distribution of errors is unsatisfactory.In general, the regression precision can be evaluated mainly from the average error for different polynomials shown in Figure1. The regression errors for polynomials ρ4, ρ6, ρ8, and ρ10 are all exceeding 5ΔE units, while for the other polynomials ρ14, ρ20, and ρ35 their color differences are all acceptable. The figure shows that polynomial ρ35 is the most accurate, but its precision is very close to polynomial ρ20. In addition, too many terms of the polynomial may increase the difficulty of coefficients-solving process [14], so the polynomial ρ20 should be most suitable for RGB color signal processing.Figure 1
Different polynomials’ average error of theRGB training samples.Using the regression coefficients from the training data, the relationship betweenRGB and CIELAB signals can be described. To verify the precision of different polynomials for the whole range of RGB signals, the testing sample data should be used. For the testing data which is outside the range of training data, the RGB’s corresponding simulation CIELAB values can be computed using the relationship obtained above, and the errors can also be calculated by comparing the measured values. For the different polynomials and the testing data, the errors are shown below.From Table2, for the RGB and CIELAB color signals, the signal processing precision of different polynomials for testing data is similar to the training data shown in Table 1, and the acceptable forms of polynomials are ρ14, ρ20, and ρ35, respectively.Table 1
The color differences of RGB training sample data for different polynomials.
Mean ΔE
Min ΔE
Max ΔE
Std. ΔE/var. ΔE
ρ
4
10.9106
1.1670
37.5866
5.4642/29.8576
ρ
6
9.4279
1.1217
35.0208
4.9875/24.8748
ρ
8
8.7390
1.3663
30.4869
4.3973/19.3359
ρ
10
6.5693
0.6068
23.8992
3.1985/10.2302
ρ
14
3.9042
0.1371
20.8066
2.9348/8.6133
ρ
20
2.5073
0.1744
16.2386
1.9063/3.6338
ρ
35
1.8134
0.0788
8.5782
1.2957/1.6787Table 2
The color difference of RGB testing data for different polynomials.
Mean ΔE
Min ΔE
Max ΔE
Std. ΔE/var. ΔE
ρ
4
9.4627
1.2718
29.3782
4.8767/23.7826
ρ
6
8.2973
0.7076
28.1559
4.0859/16.6949
ρ
8
7.7875
1.3483
21.7351
3.5269/12.4392
ρ
10
5.5913
0.2520
17.0864
3.1055/9.6441
ρ
14
3.8443
0.1035
19.2380
3.1141/9.6974
ρ
20
2.6939
0.0751
16.2386
2.3415/5.4824
ρ
35
2.3701
0.2423
13.1631
1.8308/3.3518For the purpose of testing the different polynomials’ performance for CMYK signals, the training sample data of EPSON9880 printer are used to obtain the relationship between CMYK and CIELAB. Because the errors are too large for polynomials with few terms such asρ4 and ρ6, polynomials ρ8, ρ10, ρ14, ρ20, and ρ35 are only tested for CMYK signals. With the obtained polynomial coefficients solved with training data, for the 216 testing patches of testing sample data, the color differences are shown in Table 3.Table 3
The color difference of CMYK testing data for different polynomials.
Mean ΔE
Min ΔE
Max ΔE
Std. ΔE/var. ΔE
ρ
8
2.5541
0.5667
5.7650
1.1044/1.2197
ρ
10
1.8277
0.1946
4.4006
0.7523/0.5659
ρ
14
1.5223
0.4062
3.1525
0.5682/0.3229
ρ
20
1.2485
0.3233
2.6879
0.4690/0.2200
ρ
35
1.1592
0.1275
2.5541
0.4179/0.1746It can be seen that, for the CMYK devices, most of the polynomials perform well with the average error below3ΔE. This is mainly because the color-rendering properties of CMYK printers surpass those of the RGB monitors, such as the regularity of color gamut and color consistence [1, 19]. On the whole the fourth-degree polynomial ρ35 with the largest number of terms has the smallest color difference and its error distribution is also ideal. The performance of the third-degree polynomial ρ20 is close to the ρ35, which indicates that the preferred polynomials for CMYK and CIELAB signal processing should be the ρ35 or ρ20.


## 3.2. Selection of CIE Color Spaces during Color Signal Processing

During the signal processing for different color devices, the use of two CIE color spaces, CIEXYZ and CIELAB spaces, often brings in different color errors. Hence, it is important to find the appropriate CIE color space for the specified color devices and polynomials.In the experiment, the CIEXYZ and CIELAB color spaces are tested, respectively, onRGB monitors and on CMYK printers to compare their color errors. For the IBM monitor, when the CIEXYZ is selected as device-connection space, the errors of different polynomials are shown in Table 4. Corresponding to the color errors listed in Table 2 where the CIELAB is the device-connection space, the mean color differences corresponding to these two CIE spaces are compared in Figure 2.Table 4
The color difference of RGB monitor for different polynomials using CIEXYZ space.
Mean ΔE
Min ΔE
Max ΔE
Std. ΔE/var. ΔE
ρ
8
13.2758
1.3978
93.3629
9.8920/97.8520
ρ
10
4.9488
0.2077
19.3735
3.1563/9.9620
ρ
14
2.6386
0.2299
13.6674
2.0601/4.2441
ρ
20
2.5788
0.1656
17.8186
2.0549/4.2226
ρ
35
2.2450
0.3136
10.8137
1.661/2.7585Figure 2
Comparison of theRGB monitor’s color errors between CIEXYZ and CIELAB spaces.As the figure shows, for theRGB monitors, the polynomials perform somewhat differently with the CIEXYZ and CIELAB color spaces. For the low-degree polynomials, the color error of signal processing with CIEXYZ space is greater than with CIELAB space, and along the increasing of polynomial terms, the influence of device-connection color space becomes very little. So when the polynomials of less than 10 terms are used in RGB signal processing, the CIELAB color space is recommended as the device-connection space, while for the case when the polynomial terms are between 10 and 20, the CIEXYZ space is suggested.To test the influence of CIE color space for CMYK devices, the errors of CIEXYZ and CIELAB color spaces are also compared on Epson9889 printer. The CMYK signal processing errors with CIELAB space have been listed in Table3; for the polynomials of 8 to 35 terms, the color errors with CIEXYZ space are recorded in Table 5.Table 5
The color differences of CMYK printer for different polynomials using CIEXYZ space.
Mean ΔE
Min ΔE
Max ΔE
Std. ΔE/var. ΔE
ρ
8
11.9339
0.9828
37.2081
7.6683/58.8035
ρ
10
4.6879
0.5701
19.8019
3.5922/12.9040
ρ
14
4.3475
0.5852
20.1294
3.5023/12.2659
ρ
20
1.5809
0.3099
4.0116
0.6906/0.4769
ρ
35
1.2455
0.2666
3.0209
0.4968/0.2468It can be seen that, for the second-degree polynomialsρ10 and ρ14, third-degree polynomial ρ20, and fourth-degree polynomial ρ35, all the color errors are acceptable. Taking account of the precision and computing efficiency, the third-degree polynomial ρ20 is the most suitable model for CMYK devices using CIEXYZ color space. In addition, similar to the comparison of RGB signal processing with different device-connection spaces, the CMYK signal processing using CIELAB and CIEXYZ spaces is compared in Figure 3.Figure 3
Comparison of the CMYK printer’s color errors between CIEXYZ and CIELAB spaces.It can be seen that, within the CMYK signal processing based on polynomial regression method, the precision of color conversion using CIELAB color space is higher than that of using CIEXYZ space for a majority of the second, third, and fourth degree of polynomials. Therefore, in CMYK signal processing, the CIELAB color space is preferred as the device-connection space for the polynomial regression model.


## 3.3. Linearization during Color Signal Processing

For theRGB devices, the linearization is often applied to calibrate the color device’s gray balance [20], in which the RGB signals are firstly converted into lightness signals before color conversion process. In the experiment, the linearization is described as follows:
(16)Rl=fR(R),Gl=fG(G),Bl=fB(B),
where Rl, Gl, Bl, respectively, are the lightness signals, respectively, and fR/G/B are the linearization functions which are described as follows:
(17)fC(C)=α0+α1C+α2C2+α3C3,
where C stands for one of the colors R, G, or B and αi(i=0,1,2,3) are the coefficients.Within theRGB color signals processing, the device colors RGB are firstly linearized into the lightness signals (RGB)l, then the relationship between (RGB)l and CIELAB is obtained by polynomial regression model, and at last the color error is calculated by using the measured colors within the testing data. For the IBM monitor in the experiment, the color differences of different linearized polynomials are shown in Table 6.Table 6
The RGB color difference of different polynomials with linearization applied.
Mean ΔE
Min ΔE
Max ΔE
Std. ΔE/var. ΔE
ρ
8
8.0337
0.8814
18.4480
3.0490/9.2965
ρ
10
4.5930
0.4682
17.8630
2.6520/7.0334
ρ
14
3.9435
0.3297
18.0798
2.6400/6.9698
ρ
20
2.7392
0.4485
14.3771
1.8530/3.4335
ρ
35
2.3601
0.2703
9.6480
1.3645/1.8619To test the influence of the linearization on the color signal processing precision, the two groups of color errors in Tables2 and 6 are compared. As shown in Figure 4, the errors are very close for the two processes, so a conclusion is reached that for the RGB devices the linearization has little impact on the signal processing precision especially for the polynomials of degree greater than or equal to three.Figure 4
TheRGB monitor’s color error comparison with linearization.Similarly, to test the linearization for CMYK signal processing, the color errors of testing sample of EPSON printer are recorded in Table7, and the comparison with the signal processing without linearization in Table 3 is described in Figure 5.Table 7
The CMYK color difference of different polynomials with linearization applied.
Mean ΔE
Min ΔE
Max ΔE
Std. ΔE/var. ΔE
ρ
8
4.2971
0.5850
10.9636
2.1104/4.4537
ρ
10
2.1092
0.1526
4.5244
0.8569/0.7343
ρ
14
1.8885
0.6387
4.2186
0.7901/0.6242
ρ
20
1.5002
0.2708
3.3267
0.5553/0.3084
ρ
35
1.1958
0.1047
3.0043
0.5777/0.3338Figure 5
The CMYK printer’s color error comparison with linearization.It can be seen that, for the CMYK printers, the precision of signal processing with linearization is lower than that without linearization for most polynomial models. In some cases the linearization process does not improve the CMYK signal processing precision, so it is not advisable for CMYK devices calibration.


## 4. Conclusions

In this paper, the polynomial regression model is used for RGB and CMYK color signal processing. For the purpose of improving color signal processing precision, the item number and degree of the polynomials are tested. By comparing the color errors within the color signal processing, the appropriate polynomial expressions forRGB and CMYK color devices are obtained. In addition, the parameters of device-connection color space and linearization are tested. In general, the RGB and CMYK color signal processing by employing the polynomial regression method can be concluded as follows.(1)
During theRGB and CMYK color signal processing with polynomials regression, it is advised to use the third-degree polynomials ρ20 or fourth-degree polynomials ρ35. Taking into account the coefficient solving process and the color errors in experiment, the third-degree ρ20 is the most appropriate model.(2)
When the CIE color space is used as the device-connection space, for theRGB devices, the signal processing precision is higher with CIEXYZ space than with CIELAB space for polynomials including 10 to 20 terms, while in other cases the CIELAB is more precise. For CMYK devices, in most cases the CIELAB color space performs better than the CIEXYZ space.(3)
When the linearization is added in the color signal processing, the precision improves somewhat for parts of polynomials for theRGB devices, while for the CMYK devices, the addition of linearization reduces the signal processing accuracy instead, so the linearization is not recommended to be used within the CMYK signal processing.



---

*Source: 418651-2014-07-03.xml*