A Modeling Approach on the Correction Model of the Chromatic Aberration of Scanned Wood Grain Images

Abstract

There always exists subjective and objective color differences between digital wood grain and real wood grain, making it difficult to replicate the color of natural timber. Therefore, we described a novel method of correcting the chromatic aberration of scanned wood grain to maximally restore the objective color information of the real wood grain. A point-to-point correction model of chromatic aberration between the scanned wood grain and the measured wood grain was established based on Circle 1 by adjusting the three channels ($sR$, $sG$, and $sB$) of the scanned images. A conversion of the color space was conducted using the mutual conversion formulas. The color change of the scanned images before and after the correction was evaluated through the $L^{\ast }{a}^{\ast }{b}^{\ast }$ color-mode-based $\Delta E^{\ast }$. and the lαβ color-model-based CIQI (Color Image Quality Index) and CQE (Color Quality Enhancement). The experimental results showed that the chromatic aberration $\Delta E^{\ast }$ between the scanned wood grain and the measured wood grain decreased and the colorfulness index CIQI of the scanned wood grain increased for most wood specimens after the correction. The values of $\Delta E^{\ast }$ of the twenty kinds of wood specimens decreased by an average of 3.1 in Circle 1 and 2.3 in Circle 2, thus the correction model established based on Circle 1 was effective. The color of the scanned wood grain was more consistent with that of the originals after the correction, which would provide a more accurate color information for the reproductions of wood grain and had an important practical significance.

1. Introduction

The color of wood or wood grain depends on the absorption and reflection properties of wood components to different wavelengths in the visible light. It is known that color is one of the most important parameters of the aesthetic performance of wood surface, and yet wood components are vulnerable to weathering factors, such as UV radiation and humidity, leading to the discoloration [1,2,3]. In order to imitate the color of natural timber and reduce the influence of environmental factors on the color of wood surface, technologies such as wood staining [4], decorative paper, and 3D-printed wood grain [5] have been applied. Most of these surface decoration methods require the image acquisition of the natural wood grain to obtain the digital grain, and they are then applied to mass reproductions [6]. However, there always exists subjective and objective color differences between the digital wood grain and the real wood grain, making it difficult to replicate the color of natural timber.

Sang et al. [5] printed the scanned images of a real wood surface by UV-curable inkjet technology, and the replicated wood grain on wood-based panels presented a very vivid and realistic effect in color, texture, and touch under natural light. However, it was not mentioned whether there was any chromatic aberration between the real wood surface and the replicated surface. Nordvik et al. [7] studied the experienced difference between wood seen physically and the picture of it on a computer screen. Forsthuber et al. [8] proved that there was a strong correlation between the values of the relative color change obtained by a spectrophotometer and by a calibrated photo-scanner, while only a moderate correlation was obtained comparing the absolute colorimetric values. In fact, it is because the median filter denoising conducted by the authors during the image analysis and the denoising mainly works on the blue channel. As the blue filter is the least efficient in terms of allowing light to pass through, the scanner noise is typically strongest in the blue channel [9]. To obtain the scanned wood grain whose color is maximally consistent with that of the real grain, image preprocessing such as denoising should be avoided during the image acquisition. In addition, the color-management software compares the known colors of the target (CIE-LAB values in the reference file) with the corresponding RGB values of the scanned target file. From this comparison, the software builds a profile to help applications such as Photoshop or Lightroom understand the input device’s RGB values; that is, calibrating and profiling input devices helps us to better understand the colors of our originals [10], and it apparently has no effects on the scanned colors. Feng et al. [11] adjusted the $L^{\ast }{a}^{\ast }{b}^{\ast }$ values in Photoshop CC 2018 (Adobe Inc., San Jose, CA, USA) to decrease the chromatic aberration of the scanned images. The $L^{\ast }{a}^{\ast }{b}^{\ast }$ differences between the real wood grain and the digital grain reduced, and yet it was a subjective operation process requiring constant adjustments to achieve the ideal result.

From the above, comparing the chromatic aberration between the real wood grain and the digital grain has been proven to be an interesting exploration. However, color correction is mostly used in the distorted images caused by shooting environment [12,13] or in the chromatic aberration of optics themselves. Han et al. [14] proposed an RGBN (Red-green-blue and Near-infrared) camera color correction pipeline to solve color bias and the multicollinearity of the imaging devices. Kozubek et al. [15] invented a novel method of chromatic aberration measurement and software correction using sub-resolution bead imaging and computer image analysis. Additionally, in terms of studies on real wood color, Yao et al. [16] investigated the effect of heat on color and surface chemistry of wood flour (WF). Forsthuber et al. [17] developed a novel color measurement procedure to estimate the color of early- and latewood in the annual rings of wood. Nasir et al. [18] used the $L^{\ast }{a}^{\ast }{b}^{\ast }$ color measuring system to develop a machine learning model for predicting the mechanical properties of artificially weathered fir, alder, oak, and poplar wood. Studies on the digital wood color such as the color analysis of digital Birch wood after the heat treatment [19] and those about the color image itself, such as the perceptually optimized enhancement of color in images using JND (Just Noticeable Difference) transform and color constancy [20] and the novel method on color image enhancement [21] have also been widely reported.

Overall, only a few studies focused on the correction of chromatic aberration between real wood grain and digital grain. In order to maximally restore the objective color information of the real wood grain by digitalizing, we reduced the chromatic aberration by establishing a point-to-point correction model of $sR$, $sG$, and $sB$ between the measured wood grain and the scanned wood grain. The color information of digital images and that of real objects was connected through this method. The $L^{\ast }{a}^{\ast }{b}^{\ast }$ values of the scanned wood grain after the correction were more consistent with those of the real grain and they could be saved as wood material stock images. Most importantly, a more accurate color information of the natural timber would be provided for the reproductions of wood grain by color correcting.

2. Materials and Methods

2.1. Pretreatment and Scanning of Wood Grain

Wood specimens from twenty kinds of tree species were selected as the research objects. The specimens free from knots were cut to the size of 65 mm × 40 mm × 10 mm (length × width × height) and polished to prepare a clear and flat surface. All these pieces were conditioned for 3 months in a climate room at 20 °C and 65% RH to constant weight.

Flatbed scanners are commonly available at low cost [8]. Prior to the scanning, two identical circles with a diameter of about 32 mm were drawn on the wood surface according to the area covered by the colorimeter when measuring colors. The twenty pieces of wood were then scanned centrally at one time by HP Smart Tank (Hewlett-Packard Co., Beijing, China) scanner with a charge coupled device CIS at 8 bits per channel and 300 ppi (pixels per inch) resolution. As mentioned above, color management had no effects on the scanned colors, and calibrating and profiling the scanner was canceled. Automatic color adjustments in the scanner settings were turned off and no denoising or sharpening options were enabled in order to maximally reproduce the colors of originals.

The scanned wood grain was stored in a lossless TIFF format and cropped to the size of 60 mm × 35 mm (length × width). Hence, the size of each wood image is 709 × 413 pixels with 256 gray levels. The twenty kinds of wood specimens are shown in Figure 1 and

Table 1. Basic information of the twenty tree species.

Serial Numbers	Wood Name	Tree Species	Production Area
1	Buxus	B. sinica	China
2	Guajacum	G. officinale	China
3	Juglans spp.	J. nigra	China
4	Guibourtia spp.	G. echie	Africa
5	Diospyros spp.	D. melanoxylon	Africa
6	Millettia spp.	M. laurentii	Africa
7	Dalbergia	D. melanoxylon	Myanmar
8	Acer spp.	A. nigrum	North America
9	Quercus spp.	Q. mongolica	Canada
10	Zelkova spp.	Z. schneideriana	China
11	Guibourtia spp.	G. tessmannii	Africa
12	Parinari	P. campestris	Brazil
13	Colophospermum	C. mopane	Africa
14	Pterocarpus spp.	P. erinaceus	Vietnam
15	Phoebe	P. zhennan	China
16	Myroxylon spp.	M. balsamum	Brazil
17	Pterocarpus spp.	P. macrocarpus	Vietnam
18	Pterocarpus spp.	P. tinctorius	Africa
19	Swartzia spp.	S. madagascariensis	Africa
20	Pterocarpus spp.	P. soyauxii	Africa

.

2.2. Color Measurement and Calculation of Wood Grain

2.2.1. Conversion of Color Space

The images we acquire are usually the sRGB format of nonlinear color space [22], and it is necessary to convert the sRGB to the RGB of linear color space, i.e., gamma correction [23], then to XYZ and $L^{\ast }{a}^{\ast }{b}^{\ast }$. However, roundtripping from RGB to $L^{\ast }{a}^{\ast }{b}^{\ast }$ to RGB is a destructive process on 8-bit/channel images [9]. In consequence, the original sRGB format of the scanned images was kept, and the colorimetric values of sRGB and $L^{\ast }{a}^{\ast }{b}^{\ast }$ were calculated using the mutual conversion formulas.

(1)

sRGB to $L^{\ast }{a}^{\ast }{b}^{\ast }$

The conversion of gamma correction is as Equation (1). The conversion of RGB to XYZ is as Equation (2). The conversion of XYZ to $L^{\ast }{a}^{\ast }{b}^{\ast }$ is as Equations (3) and (4). As the conversion matrix under D65 lighting condition and 10° field angle was used, the corresponding stimulus values of white point $X_0$, $Y_0$, $Z_0$ are 94.811, 100, and 107.304, respectively.

$$f\left(k\right)=\{\begin{array}{r}{(\frac{k+0.055}{1.055})}^{2.4},k>0.04045\\ \frac{k}{12.92},k\le 0.04045\end{array}$$

(1)

where the value of $k$ is $\frac{sR}{255}$, $\frac{sR}{255}$, and $\frac{sR}{255}$. The values of $R$, $G$, and $B$ are $f\left(\frac{sR}{255}\right)$, $f\left(\frac{sR}{255}\right)$, and $f\left(\frac{sR}{255}\right)$, respectively.

$$\left[\begin{array}{c}X\\ Y\\ Z\end{array}\right]=\left[\begin{array}{ccc}0.4124564& 0.3575761& 0.1804375\\ 0.2126729& 0.7151522& 0.0721750\\ 0.0193339& 0.1191920& 0.9503041\end{array}\right]\left[\begin{array}{c}R\\ G\\ B\end{array}\right]$$

(2)

$$f\left(t\right)=\{\begin{array}{r}{(t)}^{1/3},t>{\left(\frac{6}{29}\right)}^3\\ \frac{1}{3}{\left(\frac{29}{6}\right)}^{2}t+\frac{16}{116},t\le {\left(\frac{6}{29}\right)}^3\end{array}$$

(3)

where the value of $t$ is $\frac{X}{X_0}$, $\frac{Y}{Y_0}$, and $\frac{Z}{Z_0}$. The values of $L^{\ast }$, $a^{\ast }$, and $b^{\ast }$ are obtained by substituting the values of $f\left(t\right)$ into Equation (4).

$$\{\begin{array}{c}{L}^{\ast }=116f\left(\frac{Y}{Y_0}\right)-16\\ a^{\ast }=500[f\left(\frac{X}{X_0}\right)-f\left(\frac{Y}{Y_0}\right)]\\ b^{\ast }=200[f\left(\frac{Y}{Y_0}\right)-f\left(\frac{Z}{Z_0}\right)]\end{array}$$

(4)

(2)

$L^{\ast }{a}^{\ast }{b}^{\ast }$ to sRGB

Similarly, the conversion of $L^{\ast }{a}^{\ast }{b}^{\ast }$ to sRGB is as Equations (5)–(8).

$$\{\begin{array}{c}\frac{Y}{Y_0}=\frac{L^{\ast }+16}{116}\\ \frac{X}{X_0}=\frac{a^{\ast }}{500}+\frac{Y}{Y_0}\\ \frac{Z}{Z_0}=\frac{Y}{Y_0}-\frac{b^{\ast }}{200}\end{array}$$

(5)

$$f\left(t\right)=\{\begin{array}{c}{t}^3,t>{\left(\frac{6}{29}\right)}^3\\ \frac{t-\frac{16}{116}}{\frac{1}{3}{\left(\frac{29}{6}\right)}^2},t\le {\left(\frac{6}{29}\right)}^3\end{array}$$

(6)

where the value of $t$ is $\frac{X}{X_0}$, $\frac{Y}{Y_0}$, and $\frac{Z}{Z_0}$. The values of $X$, $Y$, and $Z$ are $X_{0}f\left(\frac{X}{X_0}\right)$, $Y_{0}f\left(\frac{Y}{Y_0}\right)$, and $Z_{0}f\left(\frac{Z}{Z_0}\right)$, respectively.

$$\left[\begin{array}{c}R\\ G\\ B\end{array}\right]=\left[\begin{array}{ccc}3.2404542& -1.5371385& -0.4985314\\ -0.9692660& 1.8760108& 0.0415560\\ 0.0556434& -0.2040259& 1.0572252\end{array}\right]\left[\begin{array}{c}X\\ Y\\ Z\end{array}\right]$$

(7)

$$f\left(k\right)=\{\begin{array}{c}1.055k^{\frac{1}{2.4}}-0.055,k>0.0031308\\ 12.92k,k\le 0.0031308\end{array}$$

(8)

where the value of $k$ is $R$, $G$, and $B$. The values of $sR$, $sG$, and $sB$ are $255f\left(R\right)$, $255f\left(G\right)$, and $255f\left(B\right)$, respectively.

2.2.2. Color Measurement of Original Wood Grain

The color measurement of the wood specimens was conducted using a Color Meter Pro (CHNSpec Technology Co., Hangzhou, China) with D/8-SCI at D65 lighting condition and 10° field angle. As the $L^{\ast }{a}^{\ast }{b}^{\ast }$ color space has been widely used to represent surface colors [24,25,26], values of lightness $L^{\ast }$ and chromaticity $a^{\ast }$ and $b^{\ast }$ were measured.

The colorimeter was then placed on the two circles of the wood specimens. Each circle was measured three times and the average value was taken. This value represented the mean colorimetric value within the circle with a diameter of 8 mm on the wood surface due to the 8 mm measuring aperture. The 8 mm circle was also concentric with the 32 mm drawn one.

2.2.3. Color Calculation of Digital Wood Grain

As shown in Figure 2, a concentric circle with a diameter of 8 mm was drawn in the 32 mm circle, and the 8 mm circle was cut out and saved in a PNG format at 8 bits per channel and 300 ppi resolution.

The mean values of sRGB of the 8 mm circles were obtained in Photoshop CC2018, and color conversion was conducted using Equations (1)–(4) to calculate the mean values of $L^{\ast }{a}^{\ast }{b}^{\ast }$.

2.3. Correction Model of Chromatic Aberration

Since there might be a strong correlation between the values of the relative color change obtained by a spectrophotometer and by an uncalibrated scanner, it was assumed that the colorimetric values between the measured point on the original wood grain and the corresponding point on the digital grain had a certain proportion. Based on this assumption, the colorimetric values of Circle 1 were used to establish the correction model, and those of Circle 2 were used to further confirm whether the assumption was tenable.

Initially, we converted the measured $L^{\ast }{a}^{\ast }{b}^{\ast }$ values of the 8 mm circles in Circle 1 to the sRGB values. Secondly, we calculated the chromatic aberration between this converted sRGB values and the sRGB values of the corresponding point on the scanned wood grain. The correction model of the chromatic aberration of $sR$, $sG$, and $sB$ are as Equations (9)–(11), respectively.

$$\{\begin{array}{c}{I^{\prime }}_r=I_r+\frac{{\mu }_r}{s{R}_d}\times \left(s{R}_o-s{R}_d\right),s{R}_o\ge s{R}_d\\ {I^{\prime }}_r=I_r-\frac{{\mu }_r}{s{R}_d}\times \left(s{R}_d-s{R}_o\right),s{R}_o<s{R}_d\end{array},$$

(9)

$$\{\begin{array}{c}{I^{\prime }}_g=I_g+\frac{{\mu }_g}{s{G}_d}\times \left(s{G}_o-s{G}_d\right),s{G}_o\ge s{G}_d\\ {I^{\prime }}_g=I_g-\frac{{\mu }_g}{s{G}_d}\times \left(s{G}_d-s{G}_o\right),s{G}_o<s{G}_d\end{array},$$

(10)

$$\{\begin{array}{c}{I^{\prime }}_b=I_b+\frac{{\mu }_b}{s{B}_d}\times \left(s{B}_o-s{B}_d\right),s{B}_o\ge s{B}_d\\ {I^{\prime }}_b=I_b-\frac{{\mu }_b}{s{B}_d}\times \left(s{B}_d-s{B}_o\right),s{B}_o<s{B}_d\end{array},$$

(11)

where $s{R}_{o}s{G}_{o}s{B}_o$ are the converted values of the 8 mm circles on the original wood grain using Equations (5)–(8), and $s{R}_{d}s{G}_{d}s{B}_d$ are the values of the corresponding point on the digital wood grain. $\mu$ is the average values of all pixels in the scanned image [27]. $I\left(i,j\right)$ is the pixel value of one point at the image coordinate, and $I^{\prime }\left(i,j\right)$ is that after the chromatic aberration correction.

2.4. Calculation of Chromatic Aberration

The total color difference of the 8 mm circles between the digital wood grain and the original wood grain was calculated using the Equation (12) [8,26]:

$$\Delta E^{\ast }=\sqrt{\Delta L^{\ast }+\Delta a^{\ast }+\Delta b^{\ast }}$$

(12)

where $\Delta E^{\ast }$ is the comprehensive evaluation index of color difference, $\Delta L^{\ast }$ is the lightness difference, $\Delta a^{\ast }$ and $\Delta b^{\ast }$ are the chromaticity differences.

2.5. Evaluation of Colorfulness Index

Colorfulness is the attribute of chrominance information humans perceive. Panetta et al. [28] demonstrate that colorfulness can be represented effectively with combinations of image statistics. The lαβ color model proposed by Ruderman et al. [29] is based on human vision system, and it will not be affected by environment or device [30]. Therefore, to further investigate the color change of scanned wood grain before and after the correction, two metrics using the lαβ color model were introduced to evaluate the colorfulness.

The CIQI (Color Image Quality Index) colorfulness metric [31] is shown in Equation (13), where $\alpha =R-G$, $\beta =0.5\times \left(R+G\right)-B$, ${\sigma }^2$ and $\mu$ represent the variance and mean values of the two-color axes.

$$CIQ{I}_{\_Colorfulness}=\frac{\sqrt{{\sigma }_{\alpha }{}^2+{\sigma }_{\beta }{}^2}+0.3\sqrt{{\mu }_{\alpha }{}^2+{\mu }_{\beta }{}^2}}{85.59}$$

(13)

Zhang et al. [32] introduced the logarithmic operation into the colorfulness measure because the subjective sensation was proportional to the logarithm of the stimulus. As shown in Equation (14), the CQE (Color Quality Enhancement) colorfulness metric was proposed by combinations of a subset of color image quantities in the logarithmic domain [28].

$$CQ{E}_{\_Colorfulness}=0.02\times \frac{log{\sigma }_{\alpha }{}^2\times log{\sigma }_{\beta }{}^2}{\log \left({\sigma }_{\alpha }{}^2+{\sigma }_{\beta }{}^2\right)}\times \frac{log{\mu }_{\alpha }{}^2\times log{\mu }_{\beta }{}^2}{\log {\left({\mu }_{\alpha }+{\mu }_{\beta }\right)}^2}$$

(14)

3. Results and Discussion

3.1. Statistical Analyses before the Chromatic Aberration Correction

The colorimetric values of the 8 mm circles obtained by measuring and scanning were compared. The correlation diagrams between them are shown in Figure 3. Clearly there was a highly linear positive correlation of the $sR$, $sG$, and $sB$ values obtained through the two methods, and the difference between Circle 1 and Circle 2 was small. As a result, the assumption mentioned in 2.3 was valid, and the colorimetric values of Circle 1 were used to establish the correction model.

The values of $\Delta E^{\ast }$ of the 8 mm circles between the measured wood grain and the scanned wood grain before the chromatic aberration correction were calculated using Equation (12). The results are shown in

Table 2. Chromatic aberration $\Delta E^{\ast }$ of the 8 mm circles in Circle 1 and Circle 2 before the correction.

No.		1	2	3	4	5	6	7	8	9	10
$\Delta E^{\ast }$	Circle 1	6.34	4.09	1.48	2.49	2.85	1.50	2.59	4.15	3.14	3.79
	Circle 2	6.00	3.10	1.14	2.35	2.74	1.35	1.20	4.11	3.38	4.03
No.		11	12	13	14	15	16	17	18	19	20
$\Delta E^{\ast }$	Circle 1	3.00	3.05	1.16	2.28	1.89	3.89	5.18	5.23	4.46	6.09
	Circle 2	3.89	2.34	1.89	2.44	2.47	5.96	4.39	4.03	6.10	5.94

, and the corresponding cluster analysis diagram is shown in Figure 4.

Figure 4 shows that the values of $\Delta E^{\ast }$ are closer for most wood specimens sharing a similar color, but the color difference between the No.1 and No.20 specimens is conspicuous (see Figure 1). This demonstrated that a closer value of $\Delta E^{\ast }$ could not represent a similar color of the originals, and the chromatic aberration $\Delta E^{\ast }$ between the original wood grain and the digital wood grain depended on the scanner’s ability of reproducing different values of $sR$, $sG$, and $sB$. Consequently, it was not recommended to correct the chromatic aberration of different kinds of wood specimens by increasing or decreasing the same sRGB values even if their values of $\Delta E^{\ast }$ were close.

3.2. Color Comparison after the Chromatic Aberration Correction

The measured $L^{\ast }{a}^{\ast }{b}^{\ast }$ values of the 8 mm circles in Circle 1 were converted to sRGB values using Equations (5)–(8), and the sRGB values of the corresponding point on the scanned wood grain were obtained in Photoshop. The correction model was then established on the basis of the difference of sRGB values using Equations (9)–(11).

The values of $\Delta E^{\ast }$ of the 8 mm circles between the measured wood grain and the scanned wood grain after the chromatic aberration correction were calculated using Equation (12). The results are shown in

Table 3. Chromatic aberration $\Delta E^{\ast }$ of the 8 mm circles in Circle 1 and Circle 2 after the correction.

No.		1	2	3	4	5	6	7	8	9	10
$\Delta E^{\ast }$	Circle 1	0.17	0.29	0.42	0.54	0.76	0.18	0.40	0.07	0.36	0.17
	Circle 2	0.90	0.96	0.29	1.34	1.84	0.59	1.72	0.75	0.92	0.41
No.		11	12	13	14	15	16	17	18	19	20
$\Delta E^{\ast }$	Circle 1	0.31	0.25	0.32	0.24	0.17	0.20	0.33	0.31	0.23	0.54
	Circle 2	0.80	1.67	1.03	0.69	0.75	2.41	1.18	1.43	2.40	1.54

, and the changing trends of $\Delta E^{\ast }$ before and after the correction are shown in Figure 5.

It can be seen from Figure 5 that the values of $\Delta E^{\ast }$ went down sharply for most wood specimens after the chromatic aberration correction. However, the No.7 value of $\Delta E^{\ast }$ of the 8 mm circle in Circle 2 went up, and the whole values of $\Delta E^{\ast }$ of the 8 mm circles in Circle 1 showed a greater drop compared with those in Circle 2. The values of $\Delta E^{\ast }$ of the twenty kinds of wood specimens decreased by an average of 3.1 in Circle 1 and 2.3 in Circle 2.

This phenomenon might be explained by the fact that the correction model of the chromatic aberration was established based on Circle 1. It could be further summarized that if one wood surface had distinct colors, multi-point measurements were required to establish different correction models by dividing color regions, which would inevitably cause uncoordinated color block areas in the scanned images after reorganization. Also, it might be very difficult to determine suitably chromatic aberration regions in the images [33]. If the falling range of $\Delta E^{\ast }$ was required to be consistent on any point of the scanned images, the better solution was to maximally find the corresponding points between the digital wood grain and the original wood grain, and then to substitute the mean values of the scanned points and the measured points, respectively, into the correction model.

In this paper, the surface color of each wood specimen was relatively uniform. Additionally, as the values of $\Delta E^{\ast }$ of the 8 mm circles in both Circle 1 and Circle 2 decreased, the correction model established based on Circle 1 was effective. The comparison of the scanned wood grain images before and after the correction is shown in Figure 6.

According to the assessment guidelines on the classification of the overall color change used by Hrčková et al. [19], the color change is visible by medium quality filter when the value of $\Delta E^{\ast }$ is between 3 and 6. Since less than half of the wood specimens showed a change over 3 of $\Delta E^{\ast }$ after the correction, the color difference in Figure 6 was not obvious for most of them.

3.3. Colorfulness Evaluation of the Scanned Wood Grain Images

The colorfulness of the scanned wood grain images before and after the chromatic aberration correction was calculated using Equations (13) and (14), respectively. The changes of the evaluation values of CIQI and CQE after the correction are shown in

Table 4. Changes of the evaluation values of CIQI and CQE after the chromatic aberration correction.

No.	1	2	3	4	5	6	7	8	9	10
CIQI	0.050	0.021	0.006	−0.005	0.001	0.005	0.012	0.038	0.029	0.032
CQE	0.028	0.013	0.004	−0.012	0.002	0.013	0.042	0.030	0.024	0.024
No.	11	12	13	14	15	16	17	18	19	20
CIQI	0.013	−0.006	−0.007	0.005	0.013	−0.002	0.018	−0.019	−0.023	−0.007
CQE	0.012	−0.007	−0.009	0.005	0.014	0.004	0.026	−0.016	−0.018	0.008

, and the corresponding scatter diagram before and after the correction is shown in Figure 7.

In general, because the color change of scanned wood grain before and after the correction was not obvious, both the colorfulness metrics showed a low sensitivity. It could be observed from Figure 6 that the saturation of No.18 and No.19 wood grain reduced after the correction, and correspondingly, the values of CIQI and CQE of No.18 and No.19 decreased.

Notably, though the color change of most wood specimens was not obvious, it was found from Figure 6 that the saturation of the No.20 wood specimen reduced while its CQE value rose fractionally (see Table 4 and Figure 7), which indicated that the CQE method proposed by Panetta et al. [28] might not be suitable to evaluate the colorfulness of the scanned wood grain images. Accordingly, only the colorfulness index CIQI could accurately evaluate the color change of the scanned wood grain after the chromatic aberration correction of sRGB.

4. Conclusions

Even the best flatbed scanners available on the market exhibit chromatic aberration to some extent. In this paper, a point-to-point correction model of the chromatic aberration between the scanned wood grain and the measured wood grain was established by adjusting the three channels ($sR$, $sG$, and $sB$) of the scanned images to maximally restore the objective color information of the originals. The conversion of the color space was conducted using the mutual conversion formulas. The color change of scanned wood grain before and after the correction was evaluated through the $L^{\ast }{a}^{\ast }{b}^{\ast }$ color-mode-based $\Delta E^{\ast }$ and the lαβ color-model-based CIQI and CQE.

The experimental results showed that the chromatic aberration $\Delta E^{\ast }$ decreased and the colorfulness index CIQI increased for most wood specimens after the correction, while CQE had a few errors in the evaluation. The values of $\Delta E^{\ast }$ of the twenty kinds of wood specimens decreased by an average of 3.1 in Circle 1 and 2.3 in Circle 2, thus the correction model established based on Circle 1 was effective. Moreover, both the $\Delta E^{\ast }$ and CIQI showed a slight change, indicating that this scanner had already retained the colors of the originals to a large extent.

Multi-point measurements between the digital wood grain and the original wood grain were suggested on one wood surface to ensure that the overall colors of the scanned surface share a consistent change after the correction. The method was implemented by substituting the mean values of the scanned points and the measured points, respectively, into the correction model whether the surface had distinct colors or uniform colors.

Due to inevitable human errors, the measured point on the original wood grain and the corresponding point on the digital grain could not be completely coincident. Furthermore, although the influence of environmental factors on the color of wood surface could be avoided by digitalizing, the reading accuracy of the color information relied on the performance and stability of the measuring and scanning equipment. The influence of these factors on the correction model should be further studied.