Spectral Reflectance Recovery from Tristimulus Values under Multi-Illuminants.
The object's spectral reflectance is almost the definition of "fingerprinting" that predicts accurately the object appearance under arbitrary illuminants and observers. The most accurate and effective information to represent an object color by its spectral reflectance is highly desirable, since it is essential in common application scenarios, for example, printing inspection, disease diagnosis, textile color matching, and computer vision [1-4]. Normally, the spectral reflectance could be directed to acquire from the spectrophotometers and spectral cameras. Unfortunately, portability, complexity, and expensiveness of these devices limit their own availability. The most used colorimetric values, instead, can be easily procurable by digital still cameras, smartphone, and colorimeters.
The colorimetric values such as RGB or CIE XYZ tristimulus values, obtaining from three channels, only record color information under fixed viewing conditions. There are some cases in which this representation is not enough, because it is heavily dependent on the illuminant. So, the colorimetric values still vary significantly with the illuminants. Although the computation of colorimetric values is uniquely made from the object's spectral reflectance, the problem in the spectral reflectance recovery calculating from the colorimetric values is usually ill-posed. Many different mathematical methods are still widely studied to recover the spectrum, for example, the pseudoinverse method (PI) , principal component analysis (PCA) , matrix R method , non-negative matrix transformation (NNMF) , simulated annealing , compressive sensing , simplex method , and so on. Of all these methods, the PI method is an uncomplicated and straightforward solution for the ill-posed inverse problem, which shows more clearly the connection between colorimetric values and the corresponding spectrum.
In fact, spectral recovery from a set of the corresponding colorimetric values under a specified illuminant-observer condition limits only three available dimensions, which leads to a noticeable insufficiency for spectral reflectance recovery. The colorimetric values under multi-illuminants have inspired researchers to start using several different sets of the corresponding values for the spectral recovery process. Schettini and Zuffi  exploited genetic algorithms to spectrum recovery from the CIE XYZ tristimulus values under one or multi-illuminants on the premise that basis functions and cardinalities were investigated. Abed et al.  proposed the different color spaces of lookup tables (LUTs) using a novel interpolation strategy for spectrum recovery under illuminate D65 and illuminate A; they stated that the more color coordinates reconstructed LUTs, the better the spectrum recovery results could be obtained. Harifi et al.  initially applied the nonlinear regression method to virtually increase the number of tristmulus values and then used six eigenvectors in the spectral recovery process. Amiri and Amirshahi  increased virtually the CIE XYZ tristimulus values under another illuminant and then recovered the spectral reflectance using of six eigenvectors by adopting the PCA or NNMF method, respectively. Zhang et al.  initially predicted CIE XYZ tristimulus values from camera responses values according to different illuminants and then recovered the spectral reflectance through the PI method. However, the aforementioned methods are available for spectral reflectance recovery from tristimulus values under predefined reference illuminants (such as illuminant A and illuminant D65) and treat equally each of the training samples to effect on the spectral recovery process under multireference illuminants.
This study proposes a more accurate recovering spectral reflectance method from tristimulus values under multi-illuminants that solves the problems mentioned above.
The novelty of our proposed method is concerned with the reference illuminant selection strategy and local sample weighted strategy under the multi-illuminants condition. The experimental results of the proposed and traditional methods are compared to evaluate the spectral and colorimetric accuracy between the actual and estimated spectra.
2. Mathematic Background and Method
The CIE XYZ tristimulus values under multi-illuminants are simply computed by the following equations:
[mathematical expression not reproducible] (1)
K = 100 / [integral] [I.sub.i] ([lambda])[bar.y]d[lambda] (2)
where r([lambda]) denotes the object's reflectance spectrum; K represents the factor of standardization; [bar.x]([lambda]), [bar.y]([lambda]), and [bar.z]([lambda]) denote the color matching functions of the CIE standard colorimetric observer; [I.sub.i] ([lambda]) denotes the reference illuminants of index i(i = 1, ..., m); [t.sub.i] = [[[X.sub.i], [Y.sub.i], [Z.sub.i]].sup.T] denotes the object color tristimulus values under a given reference illuminant [I.sub.i] ([lambda]); and the superscript T denotes the matrix transpose. Equation (1) can then be represented in matrix notation as follows:
t = [A.sup.T]r, (3)
where t = [[[t.sub.1], ..., [t.sub.n]].sup.T] is the combining matrix of tristimulus values under multi-illuminants and [A.sup.T] is the coefficient matrix involving the reference illuminants and the CIE color matching functions. For arbitrary color of the combining matrix [??] from tristimulus values under multi-illuminants, the recovered spectra [??] can be calculated to obtain from the inverse matrix [A.sup.T] directly. So, this spectral recovery process immediately implements transformation from tristimulus value space to spectral reflectance space called the direct spectral recovery method . Since [A.sup.T] is an underdetermined matrix, the recovered spectra are obtained from the pseudoinverse matrix [([A.sup.T]).sup.+], which inevitably causes large calculation error on ill-posed inverse problems. The recovery of spectral reflectance [??] using the PI method is simply calculated by the following equation:
[mathematical expression not reproducible] (4)
where the superscript "+" is the matrix pseudoinverse.
Equation (4) creates a linear conversion between tristimulus values under multi-illuminants and the spectral reflectance, so it is assumed to be a straightforward and resultful solution. When the transformation matrix [([A.sup.T]).sup.+] is calculated, the recovery process directly yields a single linear relation. It must be accepted that for the spectral recovery process by the implementation of the PI method, each of the training samples has made an equal difference in the formation of [([A.sup.T]).sup.+], so it leads to optimization for the samples fully, but not for each individual sample. Hence, the standard pseudoinverse method leads to large calculation error between the actual and estimated spectra and generates an imprecise result. Obviously, the more similar the spectral reflectance between the testing samples and the training samples is, the more accurate the results can be generated due to the more linearly relation [16, 17]. Theoretically speaking, if the matrix [([A.sup.T]).sup.+] can be formed by the testing sample's own characteristics, the performance of the recovered spectra would be calculated optimally. In this work, to obtain the matrix [([A.sup.T]).sup.+] adapting with the characteristics of the testing sample, the local sample weighted strategy is proposed for calculating the spectral reflectance recovery. The local sample weighted strategy involves selection of the optimal training samples and weighting of the local selected samples. The local optimal training samples based on the testing sample's own characteristics are selected, which use the multicolor space Euclidean distance between the testing sample and the training samples. Specifically, the multicolor space Euclidean distance between each testing sample and the training samples can be calculated as follows:
[mathematical expression not reproducible] (5)
where [X.sub.i,test], [Y.sub.i,test], and [Z.sub.i,test] denote the tristimulus values of the testing sample under the ith reference illuminant; n demonstrates the number of the training samples; [X.sub.i,j], [Y.sub.i,j], and [Z.sub.i,j] denote the jth training sample's tristimulus values under the ith reference illuminant; and [d.sub.j] refers to the multicolor space Euclidean distance between the testing sample and the jth training samples under the multi-illuminants condition. After that the training samples arrange in the increasing order based on the [d.sub.j] values. The local optimal training samples have been obtain to select the p (1 [less than or equal to] p [less than or equal to] n) neighboring samples from the whole training samples. Meanwhile, different local training samples should choose different weighting efficients since the larger the weighting efficient a certain training sample chooses, the stronger the influence the matrix [([A.sup.T]).sup.+] determinates. So, the weight coefficient, depending on the inverse multicolor space Euclidean distance, can be calculated as
[[omega].sub.k] = 1 / [d.sub.k] + [epsilon], k = 1, 2, ..., p, (6)
where the subscript k denotes the kth local optimal training samples; [d.sub.k] denotes the multicolor space Euclidean distance between the testing sample and the kth sample of the local optimal training samples; and [epsilon] = 0.001 is used in this study. Clearly, more neighboring training samples to the testing sample would generate the larger [w.sub.k] value. This mathematical symbol wk is shown in the galley proof. The weighting matrix W is a diagonal matrix defined as
[mathematical expression not reproducible] (7)
[([A.sup.T]).sup.+] = rW[(tW).sup.+]. (8)
Indeed, the designing method has selectively controlled the influence on the formation of the matrix [([A.sup.T]).sup.+] based on the characteristics of the testing sample, which varies from one testing sample to the other one. Each of the testing samples should have its own the transformation matrix rather than a unique transformation matrix for the whole testing samples, which generates a more precise estimation result.
3. Experiment and Procedure
In this study, three different datasets were selected that include the Munsell Chips , ColorChecker SG, and Vrhel dataset . The Munsell Chips consist 1269 spectral color chips in Munsell Matte Color Book. The 140 chips of the ColorChecher SG use a X-rite i1 pro spectrophotometer to measure personally. The Vrhel dataset contains 354 samples of reflectance spectra. The spectral reflectance function is proved to sample at 10 nm intervals without impacting greatly the mathematical precision . So, the spectral reflectance of three different dataset samples range from 400 nm to 700 nm at 10 nm intervals. The Munsell Chips were selected as the training samples in this study, and their transformation matrix was used to recover three different testing samples.
To illuminate the advantages of the proposed method, the standard pseudoinverse method (method PI) and principal component analysis embedded weight regression technique suggested by Amiri and Amirshahi  were compared with the results of the spectral recovery from tristimulus values under multireference illuminants. The root mean square error (RMSE), goodness of fit coefficient (GFC), and CIELAB color difference ([DELTA][E.sub.ab]) were selected as the evaluation of the spectral recovery accuracy , while all the methods of the spectral reflectance recovery perform experimental simulation by the Matlab software.
The dominant spectral power distributions (SPDs) of the selected reference illuminants and light sources were of various distributions and were relatively smooth, which did not have considerable spiky radiance. So, the CIE illuminants (A, B, C, D50, D55, and D65) and the two actual light-emitting diode (LED) light sources (LED1: Sylvania Concord 2048794 and LED2: Photon Star CS5)  were selected as the reference illuminants so as to evaluate spectral recovery accuracy in this study. All the CIE illuminants and LEDs were sampled at a range from 400 nm to 700 nm at 10 nm intervals. To evaluate the performance of the proposed method, it is critical to determine the type of the reference illuminants, that is, what is the optimal type of reference illuminants to improve the spectral reflectance recovery accuracy. As Zhang et al.  discussed previously, a trend can be observed: the more similar the spectral power distributions (SPDs) of two selected reference illuminants are, the lower the spectral recovery accuracy is. However, this point lacks of a clear numerical explanation to verify. This study adopted the spectral angle mapper (SAM) algorithm [23, 24] to calculate the similarity of the reference illuminants as the reference illuminants selection strategy, shown as
[S.sub.[alpha],[beta]] = accos [I.sub.[alpha]] * [I.sub.[epsilon]] / [absolute value of [I.sub.[alpha]]] * [absolute value of [I.sub.[epsilon]]], (9)
where [I.sub.[alpha]] and [I.sub.[epsilon]] denote the SPDs value vectors of the reference illuminant of index [alpha] ([alpha] = 1, ..., m) and the reference illuminant of index [epsilon]([epsilon] = 1, ..., m), respectively. This method determines the SPDs similarity to treat the two reference illuminants as vectors with treating equally to each dimension and calculate the angle between the two reference illuminants, which do not adopt the vector length but rather the vector direction .
4. Results and Discussion
The proposed spectral reflectance recovery method mainly involves two optimal parameters, namely, optimal parameter of the reference illuminants and optimal parameter of the local sample weighted method. This study first investigated two optimal parameters of the spectral recovery process and finally evaluated the spectral reflectance recovery accuracy of the proposed compared with the traditional methods.
4.1. Optimal Parameter of the Reference Illuminants. The spectral reflectance recovery accuracy can be affected by the number and type of selected reference illuminants. Theoretically, increasing the number of reference illuminants can increase the spectral recovery accuracy, but obtaining extra tristimulus values by the colorimetric device would be not easily procurable. Meanwhile, according to the conclusion made by Schettini and Zuffi , the result of the spectral recovery under the three reference illuminants is less to improve the computing accuracy than under the two illuminants. In this study, the number of reference illuminants determined to two illuminants. The type of reference illuminants selected optimally are the relatively smooth illuminants rather than the considerable spiky illuminants; this is because spectral recovery accuracy was low when one or two spiky illuminants were adopted as reference illuminants . First, all the samples of the Munsell Chips, ColorChecker SG, and Vrhel dataset were calculated numerically to obtain the CIE 1964 XYZ tristimulus values under the CIE illuminants (A, B, C, D50, D55, and D65) and the two actual light-emitting diode (LED) light sources ([LED.sub.1]: Sylvania Concord 2048794, [LED.sub.2]: Photon Star CS5), which illustrated to select optimally the type of reference illuminants. In this study, we considered the Munsell Chips as the training samples to estimate the spectral reflectance of the Munsell Chips, ColorChecker SG, and Vrhel dataset under two different illuminants with the corresponding CIE 1964 XYZ tristimulus values. Table 1 shows the mean root-mean-square error (RMSE) between the actual and reconstructed spectra for the Munsell Chips, ColorChecker SG, and Vrhel dataset under two illuminant combinations. First, it is noted that the CIE 1964 XYZ tristimulus values under two different illuminants have a very significant influence on the spectrally recovered result. Second, Table 1 illustrates the spectral reflectance accuracy for the Munsell Chips presented better than for the ColorChecker SG and Vrhel dataset. As Babaei et al.  proposed in a previous study, the best optimum condition would be obtained when the testing sample is a number of the training samples. Third, the similarity of the two reference illuminants also affects spectral estimation accuracy. Table 2 shows the spectral angle mapper (SAM) statistics for two reference illuminants. As the spectral curve shape of the two reference illuminants is more dissimilar between each other, the spectral recovery accuracy will be better, as shown in Tables 1 and 2. Combining all the experimental results of Tables 1 and 2, this study selects finally the two illuminants D65 and A as the reference illuminants.
4.2. Optimal Parameter of the Local Sample Weighted Method. The number of the selected local training samples and the weighted local samples can also influence the accuracy of the spectral reflectance recovery. To analyse these issues, the spectral reflectance of the Munsell Chips, ColorChecker SG, and Vrhel dataset was recovered by using the multicolor space Euclidean distance [d.sub.j] and the weighting matrix W under the reference illuminants A and D65 with different numbers of the local Munsell Chips as the training samples. The mean RMSE between the actual and recovered spectra under different numbers of the local Munsell Chips was computed. It is noted that in selecting, dynamically, the suitable local samples, we always adopt the strategy that if the testing sample itself is included in the Munsell Chips, then this reflectance is removed from the training samples. The relationship between the mean RMSE of the Munsell Chips, ColorChecker SG, and Vrhel dataset and the number of local training samples is shown in Figure 1. It is found that the mean RMSE initially decreases with the increase in the number of the local training samples and trends basically stable in the end. To achieve the better spectral reflectance accuracy, sufficient training samples can be selected. So, we selected adaptively 100 local training samples from the Munsell Chips for the spectral recovery process.
4.3. Accuracy Evaluation of the Proposed Method. To evaluate the colorimetric and spectral performance of the proposed method, the proposed method was implemented to compare with the pseudoinverse method (method PI) and the Morteza Maali Amiri's method (method WRPCA) . The statistic comparison results among PI, WRPCA, and the proposed method are summarized in Table 3. First, it is easy find that the proposed method obviously outperforms the PI and WRPCA for the three testing datasets. Second, results in Table 3 suggest that all the methods for Munsell Chips implemented better than for the ColorChecker SG and Vrhel dataset. This finding shows that the training sample has the better performance to recover itself than other testing samples. Second, Table 3 indicates that the proposed method presents the spectral reflectance with the highest accuracy compared with the other methods.
To further accurately evaluate the performance of the proposed method, the average spectral residuals between recovered and measured spectra for the three datasets are shown in Figures 2-4. In each figure, PI, WRPCA, and the proposed methods are contrasted for the three different datasets: the Munsell Chips (Figure 2), the ColorChecker SG (Figure 3), and the Vrhel dataset (Figure 4). By the analysis of Figures 2-4, it is easy to note that the average spectral residuals are more accurate to the middle of the wavelengths than both ends. The phenomenon is consistent with the conclusion made by Wu et al. and is attributed mainly to the human visual system that obtains the XYZ tristimulus values for three methods being calculated by the CIE 1964 standard observer .
In order to illustrate the performance of the proposed method, the recovered results of the spectral reflectance for two randomly selected samples from the Munsell Chips are shown in Figures 5(a) and 5(b), ColorChecker SG are shown in Figures 5(c) and 5(d), and Vrhel dataset are shown in Figures 5(e) and 5(f). Compared with the PI and WRPCA methods, the performance of the proposed method is more accurately illustrated in Figure 5.
A method for the spectral reflectance recovery from tristimulus values under multi-illuminants was presented by adapting with the characteristics of the testing sample to obtain the transformation matrix of pseudoinverse. To obtain the optimal transformation matrix, the proposed method for spectral reflectance recovery mainly involves two optimal parameters: optimal parameter of the reference illuminants and optimal parameter of the local sample weighted method. Selecting the two illuminants A and D65 as the reference illuminants are based on the result of the spectral angle mapper (SAM) statistics for two reference illuminants. The number of the selected local training samples and the weighted local samples can be also determined for the spectral reflectance recovery.
The reflectance of the testing samples from the Munsell Chips, ColorChecker SG, and Vrhel dataset were used to evaluate the spectral reflectance accuracy for different methods in this study. The Munsell Chips were selected as the training samples. Meanwhile, the performance of three different methods, namely, PI, WRPCA, and the proposed method was assessed by the root-mean-square (RMSE), the goodness of fit coefficient (GFC), and CIE LAB color differences under illuminants F2 and F11. The spectral and colorimetrical recovery accuracy of the proposed method were compared with those of PI and WRPCA, and the results showed the proposed method is an effective spectral reflectance recovery method.
No data were used to support this study.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
This study was supported by the Shandong Provincial Natural Science Foundation, China (ZR2017LF027, ZR2017PC020); Key Research and Development Program of Shandong Province (2018GGX106009, 2019GGX105016); Project of Shandong Province Higher Educational Science and Technology Program (J17KA178, J18KA332); Open Fund of National-Local Joint Engineering Laboratory for Digitalize Electrical Design Technology (NELDED2017K003); National Natural Science Foundation of China (NSFC) (11704162); Foundation of Xuzhou City (KC18002); and the Key Laboratory of Auxiliary Chemistry and Technology for Chemical Industry, Ministry of Education (KFKT2019-03).
 E. M. Valero, Y. Hu, J. Hernandez-Andres et al., "Comparative performance analysis of spectral estimation algorithms and computational optimization of a multispectral imaging system for print inspection," Color Research & Application, vol. 39, no. 1, pp. 16-27, 2014.
 G. Wu, "Statistical characterization of skin color spectrum and its application to dermatologic diagnosis," Basic & Clinical Pharmacology & Toxicology, vol. 124, pp. 69-70, 2019.
 J. Mohtasham, A. S. Nateri, and H. Khalili, "Textile colour matching using linear and exponential weighted principal component analysis," Coloration Technology, vol. 128, no. 3, pp. 199-203, 2012.
 G. Wu, "Charamer mismatch-based spectral gamut mapping," Laser Physics Letters, vol. 16, no. 9, Article ID 095206, 2019.
 V. Babaei, S. H. Amirshahi, and F. Agahian, "Using weighted pseudo-inverse method for reconstruction of reflectance spectra and analyzing the dataset in terms of normality," Color Research & Application, vol. 36, no. 4, pp. 295-305, 2011.
 G. Wu, X. Shen, and Z. Liu, "Wavelength-sensitive-function-based spectral reconstruction using segmented principal component analysis," Optica Applicata, vol. 46, no. 3, pp. 365-374, 2016.
 Y. Zhao and R. S. Berns, "Image-based spectral reflectance reconstruction using the matrix R method," Color Research & Application, vol. 32, no. 5, pp. 343-351, 2007.
 M. M. Amiri and S. H. Amirshahi, "A hybrid of weighted regression and linear models for extraction of reflectance spectra from CIEXYZ tristimulus values," Optical Review, vol. 21, no. 6, pp. 816-825, 2014.
 D. Dupont, "Study of the reconstruction of reflectance curves based on tristimulus values: comparison of methods of optimization," Color Research & Application, vol. 27, no. 2, pp. 88-99, 2002.
 G. Wu, "Reflectance spectra recovery from a single RGB image by adaptive compressive sensing," Laser Physics Letters, vol. 16, no. 8, Article ID 085208, 2019.
 R. Schettini and S. Zuffi, "A computational strategy exploiting genetic algorithms to recover color surface reflectance functions," Neural Computing and Applications, vol. 16, no. 1, pp. 69-79, 2007.
 F. M. Abed, S. H. Amirshahi, and M. R. M. Abed, "Reconstruction of reflectance data using an interpolation technique," Journal of the Optical Society of America A, vol. 26, no. 3, pp. 613-624, 2009.
 T. Harifi, S. H. Amirshahi, and F. Agahian, "Recovery of reflectance spectra from colorimetric data using principal component analysis embedded regression technique," Optical Review, vol. 15, no. 6, pp. 302-308, 2008.
 X. Zhang, Q. Wang, J. Li, X. Zhou, Y. Yang, and H. Xu, "Estimating spectral reflectance from camera responses based on CIEXYZtristimulus values under multi-illuminants," Color Research & Application, vol. 42, no. 1, pp. 68-77, 2016.
 G. Wu, X. Shen, Z. Liu, S. Yang, and M. Zhu, "Reflectance spectra recovery from tristimulus values by extraction of color feature match," Optical and Quantum Electronics, vol. 48, no. 1, 64 pages, 2016.
 J. Liang and X. Wan, "Optimized method for spectral reflectance reconstruction from camera responses," Optics Express, vol. 25, no. 23, pp. 28273-28287, 2017.
 J. Liang, K. Xiao, M. R. Pointer, X. Wan, and C. Li, "Spectra estimation from raw camera responses based on adaptive local-weighted linear regression," Optics Express, vol. 27, no. 4, pp. 5165-5180, 2019.
 Spectral Database, University of Eastern Finland, Kuopio, Finland, http://www2.uef.fi/fi/spectral/8.
 M. J. Vrhel, R. Gershon, and L. S. Iwan, "Measurement and analysis of object reflectance spectra," Color Research & Application, vol. 19, no. 1, pp. 4-9, 1994.
 G. Wu, Z. Liu, E. Fang, and H. Yu, "Reconstruction of spectral color information using weighted principal component analysis," Optik-International Journal for Light and Electron Optics, vol. 126, no. 11-12, pp. 1249-1253, 2015.
 G. Wu, Z. Liu, and J. Zhang, "Spectral color reproduction from CIE tristimulus values using a node address array selection technique," International Journal of Signal Processing, Image Processing and Pattern Recognition, vol. 8, no. 9, pp. 141-150, 2015.
 Spectral Power Distribution (SPD) Curves, National Gallery, London, UK, http://research.ng-london.org.uk/scientific/spd/.
 F. A. Kruse, A. B. Lefkoff, J. W. Boardman et al., "The spectral image processing system (SIPS)-interactive visualization and analysis of imaging spectrometer data," Remote Sensing of Environment, vol. 44, no. 2-3, pp. 145-163, 1993.
 A. Pelagotti, A. Mastio, A. Rosa, and A. Piva, "Multispectral imaging of paintings," IEEE Signal Processing Magazine, vol. 25, no. 4, pp. 27-36, 2008.
Guangyuan Wu [ID], (1,2) Linyong Qian, (3) Guichun Hu, (1,4) and Xiaozhou Li (1)
(1) School of Light Industry and Engineering, Qilu University of Technology (Shandong Academy of Sciences), Jinan 250353, China
(2) National-Local Joint Engineering Laboratory for Digitalized Electrical Design Technology, Wenzhou 325035, China
(3) School of Physics and Electronic Engineering, Jiangsu Normal University, Xuzhou 221116, China
(4) Key Laboratory of Auxiliary Chemistry and Technology for Chemical Industry, Ministry of Education,
Shaanxi University of Science and Technology, Xi'an 710021, China
Correspondence should be addressed to Guangyuan Wu; firstname.lastname@example.org
Received 23 February 2019; Accepted 9 October 2019; Published 3 November 2019
Academic Editor: Khalique Ahmed
Caption: Figure 1: Relationship between the number of local training samples and mean RMSE of the Munsell chips, ColorChecker SG, and Vrhel dataset.
Caption: Figure 2: Average spectral residuals between recovered and measured spectra on the Munsell chips.
Caption: Figure 3: Average spectral residuals between recovered and measured spectra on the ColorChecker SG.
Caption: Figure 4: Average spectral residuals between recovered and measured spectra on the Vrhel dataset.
Caption: Figure 5: The measured spectral reflectance curve of two randomly selected samples from (a) and (b) the Munsell Chips, (c) and (d) the ColorChecker SG, and (e) and (f) the Vrhel dataset by using PI, WRPCA, and the proposed method.
Table 1: Mean RMSE of the Munsell chips, ColorChecker SG, and Vrhel dataset under two illuminant combinations. A B C Munsell chips A 0.0218 0.0112 0.0113 B 0.0112 0.0227 0.0115 C 0.0113 0.0115 0.0232 D50 0.0097 0.0126 0.0128 D55 0.0096 0.0134 0.0136 D65 0.0095 0.0113 0.0149 [LED.sub.1] 0.0140 0.0123 0.0121 [LED.sub.2] 0.0123 0.0129 0.0125 ColorChecker SG A 0.0348 0.0188 0.0191 B 0.0188 0.0369 0.0194 C 0.0191 0.0194 0.0380 D50 0.0175 0.0213 0.0228 D55 0.0171 0.0241 0.0250 D65 0.0165 0.0228 0.0267 [LED.sub.1] 0.0250 0.0197 0.0195 [LED.sub.2] 0.0178 0.0211 0.0206 Vrhel dataset A 0.0398 0.0201 0.0207 B 0.0201 0.0425 0.0215 C 0.0207 0.0215 0.0439 D50 0.0186 0.0251 0.0266 D55 0.0179 0.0264 0.0287 D65 0.0170 0.0237 0.0273 [LED.sub.1] 0.0293 0.0251 0.0244 [LED.sub.2] 0.0195 0.0277 0.0265 D50 D55 D65 Munsell chips A 0.0097 0.0096 0.0095 B 0.0126 0.0134 0.0113 C 0.0128 0.0136 0.0149 D50 0.0228 0.0102 0.0102 D55 0.0102 0.0230 0.0101 D65 0.0102 0.0101 0.0233 [LED.sub.1] 0.0119 0.0116 0.0113 [LED.sub.2] 0.0128 0.0112 0.0117 ColorChecker SG A 0.0175 0.0171 0.0165 B 0.0213 0.0241 0.0228 C 0.0228 0.0250 0.0267 D50 0.0370 0.0160 0.0159 D55 0.0160 0.0374 0.0158 D65 0.0159 0.0158 0.0380 [LED.sub.1] 0.0181 0.0177 0.0172 [LED.sub.2] 0.0208 0.0195 0.0182 Vrhel dataset A 0.0186 0.0179 0.0170 B 0.0251 0.0264 0.0237 C 0.0266 0.0287 0.0273 D50 0.0427 0.0175 0.0173 D55 0.0175 0.0432 0.0173 D65 0.0173 0.0173 0.0440 [LED.sub.1] 0.0235 0.0226 0.0217 [LED.sub.2] 0.0277 0.0256 0.0234 [LED.sub.1] [LED.sub.2] Munsell chips A 0.0140 0.0123 B 0.0123 0.0129 C 0.0121 0.0125 D50 0.0119 0.0128 D55 0.0116 0.0112 D65 0.0113 0.0117 [LED.sub.1] 0.0221 0.0107 [LED.sub.2] 0.0107 0.0221 ColorChecker SG A 0.0250 0.0178 B 0.0197 0.0211 C 0.0195 0.0206 D50 0.0181 0.0208 D55 0.0177 0.0195 D65 0.0172 0.0182 [LED.sub.1] 0.0355 0.0176 [LED.sub.2] 0.0176 0.0356 Vrhel dataset A 0.0293 0.0195 B 0.0251 0.0277 C 0.0244 0.0265 D50 0.0235 0.0277 D55 0.0226 0.0256 D65 0.0217 0.0234 [LED.sub.1] 0.0404 0.0196 [LED.sub.2] 0.0196 0.0405 Table 2: Spectral angle mapper (SAM) statistics for two reference illuminants. A B C D50 SAM A 0.0000 0.4276 0.6339 0.4552 B 0.4276 0.0000 0.2141 0.0620 C 0.6339 0.2141 0.0000 0.2087 D50 0.4552 0.0620 0.2087 0.0000 D55 0.5280 0.1137 0.1417 0.0774 D65 0.6406 0.2249 0.0831 0.1986 [LED.sub.1] 0.3053 0.3891 0.5748 0.3990 [LED.sub.2] 0.3185 0.3558 0.5392 0.3659 D55 D65 [LED.sub.1] [LED.sub.2] SAM A 0.5280 0.6406 0.3053 0.3185 B 0.1137 0.2249 0.3891 0.3558 C 0.1417 0.0831 0.5748 0.5392 D50 0.0774 0.1986 0.3990 0.3659 D55 0.0000 0.1215 0.4653 0.4313 D65 0.1215 0.0000 0.5751 0.5410 [LED.sub.1] 0.4653 0.5751 0.0000 0.0619 [LED.sub.2] 0.4313 0.5410 0.0619 0.0000 Table 3: Statistics of the comparison results among PI, WRPCA, and the proposed method. Method Munsell chips PI WRPCA Proposed Mean RMSE 0.0095 0.0112 0.0034 Max RMSE 0.0714 0.0637 0.0370 Min GFC 0.9239 0.9088 0.9850 Mean GFC 0.9989 0.9983 0.9998 Max GFC 1.0000 1.0000 1.0000 CIELAB color difference under testing illuminant F2 Mean [DELTA][E.sub.ab] 0.2520 0.3719 0.0892 Var. (a) 0.0723 0.2046 0.0054 Max [DELTA][E.sub.ab] 2.8043 4.9694 0.6701 %[DELTA][E.sub.ab] > 3 (b) 0.0000 0.6304 0.0000 CIELAB color difference under testing illuminant F11 Mean [DELTA][E.sub.ab] 0.8966 1.2880 0.3095 Var. 0.9990 3.2857 0.0836 Max [DELTA][E.sub.ab] 11.5112 28.5594 2.4955 %[DELTA][E.sub.ab] > 3 3.9401 8.7470 0.0000 Testing samples Method ColorChecker SG PI WRPCA Proposed Mean RMSE 0.0165 0.0296 0.0120 Max RMSE 0.0536 0.1361 0.0373 Min GFC 0.9708 0.7997 0.9803 Mean GFC 0.9980 0.9731 0.9988 Max GFC 0.9999 0.9998 1.0000 CIELAB color difference under testing illuminant F2 Mean [DELTA][E.sub.ab] 0.4135 1.2824 0.2648 Var. (a) 0.1962 1.2150 0.0403 Max [DELTA][E.sub.ab] 2.7253 5.8017 1.2775 %[DELTA][E.sub.ab] > 3 (b) 0.0000 9.2857 0.0000 CIELAB color difference under testing illuminant F11 Mean [DELTA][E.sub.ab] 1.3238 2.7523 0.7869 Var. 3.1423 4.9938 0.5667 Max [DELTA][E.sub.ab] 12.7500 9.7594 3.7572 %[DELTA][E.sub.ab] > 3 14.2857 42.1429 1.4286 Method Vrhel dataset PI WRPCA Proposed Mean RMSE 0.0170 0.0365 0.0152 Max RMSE 0.1278 0.1690 0.1454 Min GFC 0.9246 0.6299 0.9061 Mean GFC 0.9955 0.9831 0.9959 Max GFC 0.9999 0.9999 1.0000 CIELAB color difference under testing illuminant F2 Mean [DELTA][E.sub.ab] 0.5012 1.4597 0.3547 Var. (a) 0.2096 1.4668 0.0817 Max [DELTA][E.sub.ab] 3.0989 5.3056 1.8061 %[DELTA][E.sub.ab] > 3 (b) 0.5650 13.2768 0.0000 CIELAB color difference under testing illuminant F11 Mean [DELTA][E.sub.ab] 1.4677 2.7012 1.2836 Var. 1.9217 5.3176 1.4947 Max [DELTA][E.sub.ab] 15.8069 13.1744 6.7993 %[DELTA][E.sub.ab] > 3 10.4520 33.6158 9.0395 (a) "Var." means the variance of the the CIELAB color difference. (b) "%[DELTA][E.sub.ab] > 3" means the percentage of testing samples with a color difference greater than 3 CIELAB units.
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Wu, Guangyuan; Qian, Linyong; Hu, Guichun; Li, Xiaozhou|
|Publication:||Journal of Spectroscopy|
|Date:||Nov 1, 2019|
|Previous Article:||Improving Neutron-Gamma Discrimination with Stilbene Organic Scintillation Detector Using Blind Nonnegative Matrix and Tensor Factorization Methods.|
|Next Article:||Superparamagnetic Clustering of Diabetes Patients Raman Spectra.|