Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Fourier based partial least squares algorithm: new insight into influence of spectral shift in “frequency domain”

Open Access Open Access

Abstract

Developments in analytical chemistry technology, especially the combination between the partial least squares and spectroscopy, have contributed significantly to predicting the chemical concentrations and discriminating similar chemical analytes. However, spectral shift is an unwanted but inevitable factor for the spectroscopic analyzer, especially in practical application, which decreases the method’s accuracy and stability. To remove the term of spectral shift completely and increase the robustness of spectroscopic analysis method, Fourier transform based partial least squares method was proposed. The approach used Fourier transform first to transform the spectral shift in the “time domain” to the phase term in the “frequency domain.” The module of the Fourier transformed spectra was then calculated. As a result, the phase term was removed (the module of the phase term is 1), which means the spectral shift term was removed completely. Finally, the spectra modules were used to build the model and validate. The approach’s advantages are: (i) that the approach provides a new insight to treat the spectral shift in spectroscopic analyzer; (ii) that the model is insensitive to spectral shift; (iii) that the approach makes partial least squares combined with spectroscopy more suitable for practical application, rather than lab experiment, because spectral shift is permitted, which means the decreased requirements of measure environment. As an example, blood species discrimination, using Raman spectroscopy, was used in order to demonstrate this approach’s effectiveness.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Partial least squares discriminant analysis is a traditional chemometrics analysis method [1–4]. Recent years, it is combined with spectroscopic technology and stands out in the analysis of chemical components because of its advantages of non-invasive, non-contact and high sensitivity [5–8]. In analytical chemistry, PLS models are nowadays widely used to predict analyte concentrations and discriminate similar analytes from spectral measurements [9,10]. For example, C. J. Meunier used PLS model to discriminate Hydrogen peroxide (H2O2) and ΔpHsignals [11]. H. Khajehsharifi used PLS to simultaneously determine the concentration of ascorbic acid, dopamine and uric acid and satisfactory results were achieved [12]. In addition, it has also been widely used by chemometricians in many other applications such as forensic science [13,14], food safety [15,16], clinical diagnosing [17,18] and public safety [19,20]. The most common approach of the detection is: firstly, a series of spectra is first measured of the sample; secondly, signal pre-processing algorithm is used to select the proper training set; finally, PLS model is built with the training data and is validated by cross-validation.

Even though such great advances have been achieved in lab experiments by many research groups, the practical application of this approach is not satisfied because of spectral shift induced by the change of the environment, the shift of the laser wavelength and the performance of the spectroscopy, which has been recognized by many researchers [21-24]. To decrease the influence of spectral shift, the simplest method is to strictly control the environment such as the temperature and humidity and using high quality laser. However, these requirements mean the increasing cost and inconvenience for the end user.

To utilize the full power of spectroscopic analyzer, many efforts have been made to avoid the influence of the spectral shift. Among these reported methods, the methods can be divided into two kinds: (i) to calibrate the original spectra; Brown and Stoyanova proposed using principal-component analysis (PCA) as an automatic method of spectral quantification to determine the frequency misalignment and phase misadjustment in a single resonant peak across a series of spectra [25]. Their results demonstrated that it has been possible to use PCA methodology to determine and correct the spectra. Westad et al. proposed using Horn and Schunck algorithm to give estimates that can move the spectra with respect to some selected reference spectrum back to reference position [26]. (ii) to derived the relationship between the spectral shift and prediction value and then compensate the prediction value; Bian et al. analyzed the relationship between the spectral shift and the prediction errors and the formula was derived [27]. When the spectral shift is known, the prediction error can be calculated with the formula. The disadvantages of these methods include: (i) all these methods can only decrease the influence of spectral shift and most of them will no longer valid when the spectral shift is large enough; (ii) the parameters to decrease the influence of spectral shift need to be modified according to the value of spectral shift which is difficult to end user; (iii) it is trouble and time-consuming to select the spectra without spectral shift as the training data. All these disadvantages are fatal and unwanted to end users. Actually, a perfect instrument for end user should be easily operated for everyone who has no background about the principle of the instrument.

Here, we investigated the possibility to remove the term of spectral shift in the spectra completely using Fourier partial least squares method. The spectra were first done Fourier transform and as a result, the spectra shift in “time domain” is transferred to the change of the phase in “frequency domain” which means that the modules of the Fourier transformed spectra in “frequency domain” have no relationship to the spectral shift. When the modules are used to build the PLS model, the model will be insensitive to the spectral shift. Different from all the methods until now, the method proposed here provides a new insight to the influence of spectral shift instead of compensation of the spectra or the prediction error. The discrimination of blood species based on Raman spectroscopy was used to demonstrate the effectiveness of FPLS. By validating 33 spectra with different spectral shift, the results indicated that FPLS compared with PLS actually removed the term of spectral shift completely and the algorithm is really simple.

2. Principle

The experimental setup and the preparation of the blood samples are similar to that in our previous work [27,28]. 31 Raman spectra of human blood and 41 Raman spectra of nonhuman blood originating from chicken, duck, dove, rabbit, rat, monkey, cat, dog, sheep and pig were used to build the PLS model. 8 Raman spectra of human blood and 25 Raman spectra of nonhuman blood originating from the same species were used to validate the PLS model. The Raman spectra of the same sample with a spectral shift of 2.82, 6.36, 9.18 and 19.77 cm−1 were obtained and were used to demonstrate that FPLS is insensitive to spectral shift.

The Raman shift k is defined as

k=1λexcitation1λ
where λexcitationis the wavelength of the laser, λis the wavelength of the scattered light.

From the above equation, we can see if the laser wavelength is changed, the wavelength of scattered light will change accordingly since Raman shift k is a constant value for a given functional group. For a given spectroscopy, each pixel is corresponding to a particular wavelength, as a result, the Raman peaks will have a shift in pixel number domain as shown in Fig. 1. When the spectral shift is smaller, the error of the prediction value can be calculated with the following equation [27]

 figure: Fig. 1

Fig. 1 The Raman spectra of blood samples with and without spectral shift.

Download Full Size | PDF

Δy=(a(k1)I'(k1)+a(k2)I'(k2)+...+a(kn)I'(kn))δk+(a(k1)I'(k1)2!+a(k2)I'(k2)2!+...+a(kn)I'(kn)2!)δk2+...+(a(k1)I'(k1)n!+a(k2)I'(k2)n!+...+a(kn)I'(kn)n!)δkn

However, when the spectral shift is large, the equation is not suitable to calculate the error. As a result, these Raman spectra with large spectral shift should be treated as outlier data for the PLS model. This means the repeat measurements should be performed after the maintenance of Raman spectroscopy or the environment. To make the measurements easier and time-saving, here, Fourier based PLS method is proposed.

Assume that the Raman spectrum for blood is I(k), the actual Raman spectrum sampled by CCD is

Isampled(k)=I(k)Πδξ/2(ξ)
where Π is the rectangular function (Rect function) which is defined as Πδξ/2(ξ)={1,-δξ/2<ξ<δξ/20,others,δξ is the size of the pixel of CCD.ξ is the coordinate in the pixel direction.

To remove the term of the spectral shift, the Fourier transform of the Raman spectrum is first performed. Assume the Fourier transform of the theoretical Raman spectrum I(k) is G(ω) and considering that the Fourier transform of Rect function is Sinc function, the Fourier transform of the sampled Raman spectrum can be calculated as

FFT(Isampled(k))=FFT(I(k)Πδξ/2(ξ))=FFT(I(k))×FFT(Πδξ/2(ξ))=G(ω)×sinc(ωδξ/2)

The module of the Fourier transformed Raman spectrum is

|FFT(Isampled(k))|=|G(ω)×sinc(ωδξ/2)|=A(ω)×|sinc(ωδξ/2)|
where A(ω) is the amplitude-frequency characteristic

For the Raman spectrum with a spectral shift δk, the module of the Fourier transform is

|FFT(Isampled(k+δk))|=|FFT(I(k+δk)Πδξ/2(ξ))|=|FFT(I(k+δk))×FFT(Πδξ/2(ξ))|=|G(ω)exp(jωδk)×sinc(ωδξ/2)|=A(ω)×|exp(jωδk)|×|sinc(ωδξ/2)|=A(ω)×|sinc(ωδξ/2)|

In Eq. (6), we find that the term of spectral shift is removed by calculating the module. Compared Eq. (5) and (6), we can find that

|FFT(I(k))|=|FFT(I(k+δk))|=A(ω)×|sinc(ωδξ/2)|

It can be seen that the modules of the Fourier transform of Raman spectra has no relationship to the spectral shift. To demonstrate the result, we do Fourier transform of the Raman spectra in Fig. 1 and the results were shown in Fig. 2. The three lines in Fig. 2 were completely overlapped which demonstrated the accuracy of Eq. (7). However, there is small difference of the amplitude of the Fourier transformed Raman spectra. The main reason is the non-uniform sampling in the wave-number domain which means the spectral shift δk is not a constant value for each Raman peaks. Assume that the Raman spectrum has n Raman peaks, then

 figure: Fig. 2

Fig. 2 The modules of the Fourier transform of the Raman spectra in Fig. 1.

Download Full Size | PDF

Isampled(k)=i=1nIsampled(ki)

Thus, Eq. (5) should be revised as

|FFT(Isampled(k))|=|i=1nGi(ω)×sinc(ωδξ/2)|=|i=1nGi(ω)|×|sinc(ωδξ/2)|

For the Raman spectrum with spectral shift, Eq. (6) should be revised as

|FFT(Isampled'(k))|=|FFT(i=1nIsampled(ki+δki)Πδξ/2(ξ))|=|i=1n(Gi(ω)exp(jωδki))×sinc(ωδξ/2)|=|i=1n(Gi(ω)exp(jωδki))|×|sinc(ωδξ/2)|
When δki=constant or δki=0, i=1nGi(ω)=i=1n(Gi(ω)exp(jωδki))which means that the value calculated by Eq. (9) and (10) is equal. Even though the amplitude of the Fourier transformed Raman spectrum is not equal for the Raman spectral with and without spectral shift, the influence of the non-uniform sampling to PLS model is small because the amplitude will be normalized first and the difference of spectral shift induced by non-uniform sampling is small.

3. Results

To avoid the influence of different measurement which may be induced by the temperature, the mixture of the blood or the focus position, we use the same Raman spectra to validate the performance of FPLS model. The spectral shift was introduced by shifting the pixel mathematically. Figure 3 presented the results obtained by PLS and FPLS for the same Raman spectra. The prediction value obtained by PLS distributed discretely and when the shift of the pixel is larger, the prediction value which is much larger than 1.5 is meaningless. However, the prediction value obtained by FPLS was distributed around 0.9 which means FPLS is insensitive to the shift of the pixel as analyzed in the theory section. The prediction value was not a constant value as analyzed because of the influence of noise. Nevertheless, the influence of the noise can be neglected even when the shift of the pixel is 100. In practical application, the Raman spectra of different measurements are more complex.

 figure: Fig. 3

Fig. 3 The prediction value. (a) is obtained by PLS when the Raman spectrum was shifted the pixel mathematically; (b) is obtained by FPLS when the Raman spectrum was shifted the pixel mathematically.

Download Full Size | PDF

To demonstrate the effectiveness of FPLS in the practical application, we measured the Raman spectra of the same blood samples with different spectral shift. The spectral shift was realized by displacing the CCD in the wavelength direction which can be easily controlled using a motor. The Raman spectra with a spectral shift of 2.82, 6.36, 9.18 and 19.77 cm−1 were obtained. The training data set included 72 Raman spectra without spectral shift. The prediction value for the Raman spectral with different spectral shift was shown in Fig. 4. The circular sign in Fig. 4 represented human blood, while the asterisk sign represented nonhuman blood. Since in the model 1 represents human blood and 2 represents nonhuman blood, the prediction value of human and nonhuman blood should be distributed around 1 and 2. Usually, to achieve the higher discrimination accuracy, algorithm such as three sigma [29] and posterioti probability [30–32] can be used to establish the threshold value to discriminate blood species. In this paper, to simply explain the problem, the threshold value was set as 1.5 [33,34]. The discrimination accuracy of the PLS model is 100% for the Raman spectra without spectral shift. When small spectral shift occurred (Fig. 4(b)), the prediction value was still distributed around 1 and 2. But the discrimination accuracy decreased and this small spectral shift was usually inconspicuous for the end user. This is a disaster for them because they will obtain the wrong results while they are not aware of it. Compared with small spectral shift, the large spectral shift can be easily observed. Since the prediction value of the Raman spectra with large spectral shift was obviously larger than 1 and 2 for human and nonhuman blood. The remedial method such as compensation or repeat measurement can be done to achieve the right results. However, no matter the compensation or the repeat measurement is undesired for the end user. Figure 5 is the average prediction value for the Raman spectra with different spectral shift. In Fig. 3, the prediction value became smaller when the shift of the pixel (in the range between 1 and 10) increased (each pixel represents 2.6 cm−1 if the nonlinearity between the Raman shift and pixel number is not considered). While in Fig. 5, the prediction value became larger when the spectral shift increased. The difference was induced by the spectral shift direction (for Fig. 3, the wavelength was shifted in the direction of short wavelength, while that for Fig. 5 was shifted in the direction of long wavelength).

 figure: Fig. 4

Fig. 4 The prediction value obtained by PLS when the Raman spectra have a spectral shift. (a) shown prediction value of the Raman spectra without spectral shift; (b)-(e) were corresponding to prediction value of the Raman spectra with a shift of 2.82 cm−1, 6.36 cm−1, 9.18 cm−1 and 19.77 cm−1.

Download Full Size | PDF

 figure: Fig. 5

Fig. 5 The average prediction value with different spectral shift when PLS model was used.

Download Full Size | PDF

The same training data set and validation data set were used to build FPLS and validate the FPLS model. The results were shown in Fig. 6. The prediction value for the Raman spectra with and without spectral shift were all distributed around 1 and 2. The discrimination accuracy for the Raman spectra without spectral shift was 100% which was the same as that of the PLS model. The discrimination accuracy is 94% for the Raman spectra with spectral shifts of both 2.82 and 6.36 cm−1. For the Raman spectra with a spectral shift of 2.82 cm−1, one human blood sample and one nonhuman blood sample were mis-discriminated while for the Raman spectra with a spectral shift of 6.36 cm−1, two nonhuman blood samples were mis-discriminated. For the Raman spectra with a spectral shift of 9.18 cm−1, only one nonhuman blood sample was mis-discriminated and the discrimination accuracy was 97%. For the Raman spectra with a spectral shift of 19.77 cm−1, the discrimination accuracy was 100%. The difference for the results of different spectral shifts was mainly induced by different focus position and the homogeneity of the blood (the blood would gradually subside with the time) in the repeat measurements. Comparing the results in Fig. 4 with that in Fig. 6, it seems that spectral shift have greater influence on the prediction value than other factors such as noise, intensity or the homogeneity of the blood. Even though the discrimination accuracy decreases, it is still larger than 94% which may be acceptable in practical application. The discrimination accuracy can be improved by establishing the threshold value with three sigma method or using dual-model method to avoid the influence of the human-like nonhuman blood samples [28]. Figure 7 was the averaged prediction value, compared with that in Fig. 5, the averaged prediction value of FPLS was found more stable and not change with the spectral shift which means that FPLS method is insensitive to the spectral shift.

 figure: Fig. 6

Fig. 6 The prediction value obtained by FPLS when the Raman spectra have a spectral shift. (a) shown prediction value of the Raman spectra without spectral shift; (b)-(e) were corresponding to prediction value of the Raman spectra with a shift of 2.82 cm−1, 6.36 cm−1, 9.18 cm−1 and 19.77 cm−1

Download Full Size | PDF

 figure: Fig. 7

Fig. 7 The average prediction value with different spectral shift when FPLS model was used

Download Full Size | PDF

Figure 8 compared the results obtained by PLS and FPLS model when the Raman spectra in the training data set have spectral shift. PLS model cannot discriminate blood species any more, while FPLS model’s discrimination accuracy is 100%. The spectral shift in the training data set made the PLS model more complex. As the spectral shift is one of the loading factors, the Raman spectra would be predicted by the combination of the Raman peaks and spectral shift. In Fig. 8 the prediction value was distributed around 2, because the spectral shift in the training data set is the same as that in the validation set. However, in the practical application, the training data set is impossible to include all the spectral shift. The result may be worse than that obtained in Fig. 8.

 figure: Fig. 8

Fig. 8 The prediction value. (a) is obtained by PLS when the Raman spectra in the training data set have spectral shift; (b) is obtained by FPLS model when the Raman spectra in the training data set have spectral shift

Download Full Size | PDF

4. Conclusion

A new approach aims to remove the term of spectral shift, instead of compensating the spectral shift, based on the combination of Fourier transform and PLS algorithm has been illustrated in this study. This method treated the influence of the spectral shift in a new direction-“frequency domain” which made the problem easier. To demonstrate the effectiveness of FPLS method, a total of 237 Raman spectra of blood originating from 105 human or nonhuman samples were collected. Among 237 Raman spectra, 132 Raman spectra have spectral shift including 33 Raman spectra with a spectral shift of 2.82 cm−1, 6.36 cm−1, 9.18 cm−1 and 19.77 cm−1. PLS and FPLS methods were compared in three aspects. First, the mathematically shifted Raman spectra which aim to decrease the influence of the noise and repeat measurements were predicted by the two models while the Raman spectra in the training data set had no spectral shift. The results indicated that FPLS method was insensitive to spectral shift; Second, the Raman spectra with a spectral shift of 2.82 cm−1, 6.36 cm−1, 9.18 cm−1 and 19.77 cm−1 induced by displacing CCD were predicted by the two models while the Raman spectra in the training data set had no spectral shift. The results indicated two facts: (1) the influence of the spectral shift for PLS models is more obvious compared with the factors such as noise and intensity which means that it is necessary to avoid the influence of the spectral shift; (2) the Fourier transform based on partial least squares is indeed insensitive to spectral shift in the validating data set; Third, the Raman spectra in training data set and validating data set both had a random spectral shift, and the results also indicated FPLS method was insensitive to spectral shift.

The three experiments demonstrated the effectiveness of FPLS method. Compared with other methods reported until now, FPLS method is simpler and can solve the influence of spectral shift once and for all by removing the term of spectral shift in “frequency domain”. This method is important for Raman or infrared spectroscopy in practical application and readers in many applications such as forensic science, food safety, drug testing, chemistry analysis or clinical diagnosis would be interested in this method.

Funding

National Key R&D Program of China (2018YFF01011104, 2016YFC1000701); National High Technology Research and Development Program of China (863 Program 2015AA021105); the Natural Science Foundation of Jiangsu Province (BK20180220); Key Project of Jiangsu Province (Grant No. BE2016090 and BE2016005-2); Jiangsu Postdoctoral Research Foundation (Grant No. 1701045B, 1701046B).

Acknowledgments

The authors gratefully acknowledge Y. B. Tian, J. Wang and N. Wang for the help of contacting the supplier of the blood samples and Wanrong Gao for improving the English of the manuscript.

References

1. J. Gottfries, K. Blennow, A. Wallin, and C. G. Gottfries, “Diagnosis of dementias using partial least squares discriminant analysis,” Dementia 6(2), 83–88 (1995). [PubMed]  

2. S. Kasemsumran, N. Suttiwijitpukdee, and V. Keeratinijakal, “Rapid classification of turmeric based on DNA fingerprint by near-infrared spectroscopy combined with moving window partial least squares-discrimination analysis,” Anal. Sci. 33(1), 111–115 (2017). [CrossRef]   [PubMed]  

3. I. Marquetti, J. V. Link, A. L. G. Lemes, M. B. S. Scholz, P. Valderrama, and E. Bona, “Partial least square with discriminant analysis and near infrared spectroscopy for evaluation of geographic and genotypic origin of arabica coffee,” Comput. Electron. Agric. 121, 313–319 (2016). [CrossRef]  

4. F. Allegrini and A. C. Olivieri, “IUPAC-consistent approach to the limit of detection in partial least-squares calibration,” Anal. Chem. 86(15), 7858–7866 (2014). [CrossRef]   [PubMed]  

5. X. X. Zhang, J. H. Yin, Z. H. Mao, and Y. Xia, “Discrimination of healthy and osteoarthritic articular cartilages by Fourier transform infrared imaging and partial least squares-discriminant analysis,” J. Biomed. Opt. 20(6), 060501 (2015). [CrossRef]   [PubMed]  

6. K. C. Doty and I. K. Lednev, “Differentiation of human blood from animal blood using Raman spectroscopy: a survey of forensically relevant species,” Forensic Sci. Int. 282, 204–210 (2018). [CrossRef]   [PubMed]  

7. K. H. Wong, V. Razmovski‐Naumovski, K. M. Li, G. Q. Li, and K. Chan, “The quality control of two Pueraria species using Raman spectroscopy coupled with partial least squares analysis,” J. Raman Spectrosc. 46(4), 361–368 (2015). [CrossRef]  

8. L. Zhang, Q. Li, W. Tao, B. Yu, and Y. Du, “Quantitative analysis of thymine with surface-enhanced Raman spectroscopy and partial least squares (PLS) regression,” Anal. Bioanal. Chem. 398(4), 1827–1832 (2010). [CrossRef]   [PubMed]  

9. D. M. Haaland and E. V. Thomas, “Partial least-squares methods for spectral analyses. 1. Relation to other quantitative calibration methods and the extraction of qualitative information,” Anal. Chem. 60(11), 1193–1202 (1988). [CrossRef]  

10. M. Singh, V. Karki, R. K. Mishra, A. Kumar, C. P. Kaushik, X. Mao, R. E. Russo, and A. Sarkar, “Analytical spectral dependent partial least squares regression: a study of nuclear waste glass from thorium based fuel using LIBS,” J. Anal. At. Spectrom. 30(12), 2507–2515 (2015). [CrossRef]  

11. C. J. Meunier, E. C. Mitchell, J. G. Roberts, J. V. Toups, G. S. McCarty, and L. A. Sombers, “Electrochemical selectivity achieved using a double voltammetric waveform and partial least squares regression: differentiating endogenous hydrogen peroxide fluctuations from shifts in pH,” Anal. Chem. 90(3), 1767–1776 (2018). [CrossRef]   [PubMed]  

12. H. Khajehsharifi, E. Pourbasheer, H. Tavallali, S. Sarvi, and M. Sadeghi, “The comparison of partial least squares and principal component regression in simultaneous spectrophotometric determination of ascorbic acid, dopamine and uric acid in real samples,” Arab. J. Chem. 10, S3451–S3458 (2017). [CrossRef]  

13. V. Sikirzhytski, K. Virkler, and I. K. Lednev, “Discriminant analysis of Raman spectra for body fluid identification for forensic purposes,” Sensors (Basel) 10(4), 2869–2884 (2010). [CrossRef]   [PubMed]  

14. A. G. Ryder, G. M. O’Connor, and T. J. Glynn, “Identifications and quantitative measurements of narcotics in solid mixtures using near-IR Raman spectroscopy and multivariate analysis,” J. Forensic Sci. 44(5), 12031J (1999). [CrossRef]  

15. J. Lim, G. Kim, C. Mo, M. S. Kim, K. Chao, J. Qin, X. Fu, I. Baek, and B. K. Cho, “Detection of melamine in milk powders using near-infrared hyperspectral imaging combined with regression coefficient of partial least square regression model,” Talanta 151, 183–191 (2016). [CrossRef]   [PubMed]  

16. S. He, W. Xie, W. Zhang, L. Zhang, Y. Wang, X. Liu, Y. Liu, and C. Du, “Multivariate qualitative analysis of banned additives in food safety using surface enhanced Raman scattering spectroscopy,” Spectrochim. Acta A Mol. Biomol. Spectrosc. 137, 1092–1099 (2015). [CrossRef]   [PubMed]  

17. H. Nawaz, N. Rashid, M. Saleem, M. Asif Hanif, M. Irfan Majeed, I. Amin, M. Iqbal, M. Rahman, O. Ibrahim, S. M. Baig, M. Ahmed, F. Bonnier, and H. J. Byrne, “Prediction of viral loads for diagnosis of Hepatitis C infection in human plasma samples using Raman spectroscopy coupled with partial least squares regression analysis,” J. Raman Spectrosc. 48(5), 697–704 (2017). [CrossRef]  

18. Q. Yang, S. S. Lin, J. T. Yang, L. J. Tang, and R. Q. Yu, “Detection of inborn errors of metabolism utilizing GC-MS urinary metabolomics coupled with a modified orthogonal partial least squares discriminant analysis,” Talanta 165, 545–552 (2017). [CrossRef]   [PubMed]  

19. T. Lu, Y. Yuan, X. He, M. Li, X. Pu, T. Xu, and Z. Wen, “Simultaneous determination of multiple components in explosives using ultraviolet spectrophotometry and a partial least squares method,” RSC Advances 5(17), 13021–13027 (2015). [CrossRef]  

20. J. B. Sleiman, B. Bousquet, N. Palka, and P. Mounaix, “Quantitative Analysis of Hexahydro-1,3,5-trinitro-1,3,5, Triazine/Pentaerythritol Tetranitrate (RDX-PETN) Mixtures by Terahertz Time Domain Spectroscopy,” Appl. Spectrosc. 69(12), 1464–1471 (2015). [CrossRef]   [PubMed]  

21. T. Brekke, O. M. Kvalheim, and E. Sletten, “Prediction of physical properties of hydrocarbon mixtures by partial-least-squares calibration of carbon-13 nuclear magnetic resonance data,” Anal. Chim. Acta 223, 123–134 (1989). [CrossRef]  

22. J. Vogels, A. C. Tas, J. Venekamp, and J. Van Der Greef, “Partial linear fit: a new NMR spectroscopy preprocessing tool for pattern recognition applications,” J. Chemometr. 10(5–6), 425–438 (1996). [CrossRef]  

23. M. Spraul, P. Neidig, U. Klauck, P. Kessler, E. Holmes, J. K. Nicholson, B. C. Sweatman, S. R. Salman, R. D. Farrant, E. Rahr, C. R. Beddell, and J. C. Lindon, “Automatic reduction of NMR spectroscopic data for statistical and pattern recognition classification of samples,” J. Pharm. Biomed. Anal. 12(10), 1215–1225 (1994). [CrossRef]   [PubMed]  

24. F. Vogt and K. Booksh, “Influence of wavelength-shifted calibration spectra on multivariate calibration models,” Appl. Spectrosc. 58(5), 624–635 (2004). [CrossRef]   [PubMed]  

25. T. R. Brown and R. Stoyanova, “NMR spectral quantitation by principal-component analysis. II. Determination of frequency and phase shifts,” J. Magn. Reson. B. 112(1), 32–43 (1996). [CrossRef]   [PubMed]  

26. F. Westad and H. Martens, “Shift and intensity modeling in spectroscopy—general concept and applications,” Chemometr. Intel. Lab. 45(1–2), 361–370 (1999). [CrossRef]  

27. H. Bian and J. Gao, “Error analysis of the spectral shift for partial least squares models in Raman spectroscopy,” Opt. Express 26(7), 8016–8027 (2018). [CrossRef]   [PubMed]  

28. H. Bian, P. Wang, N. Wang, Y. Tian, P. Bai, H. Jiang, and J. Gao, “Dual-model analysis for improving the discrimination performance of human and nonhuman blood based on Raman spectroscopy,” Biomed. Opt. Express 9(8), 3512–3522 (2018). [CrossRef]   [PubMed]  

29. P. Bai, J. Wang, H. Yin, Y. Tian, W. Yao, and J. Gao, “Discrimination of human and nonhuman blood by Raman spectroscopy and partial least squares discriminant analysis,” Anal. Lett. 50(2), 379–388 (2017). [CrossRef]  

30. C. Y. Wang, C. T. Chen, C. P. Chiang, S. T. Young, S. N. Chow, and H. K. Chiang, “A probability-based multivariate statistical algorithm for autofluorescence spectroscopic identification of oral carcinogenesis,” Photochem. Photobiol. 69(4), 471–477 (1999). [CrossRef]   [PubMed]  

31. C. Botella, J. Ferré, and R. Boqué, “Classification from microarray data using probabilistic discriminant partial least squares with reject option,” Talanta 80(1), 321–328 (2009). [CrossRef]   [PubMed]  

32. G. Aliakbarzadeh, H. Parastar, and H. Sereshti, “Classification of gas chromatographic fingerprints of saffron using partial least squares discriminant analysis together with different variable selection methods,” Chemometr. Intell. Lab. 158, 165–173 (2016). [CrossRef]  

33. Å. Eriksson, K. Persson Waller, K. Svennersten-Sjaunja, J.-E. Haugen, F. Lundby, and O. Lind, “Detection of mastitic milk using a gas-sensor array system (electronic nose),” Int. Dairy J. 15(12), 1193–1201 (2005). [CrossRef]  

34. D. Cozzolino, A. Chree, J. R. Scaife, and I. Murray, “Usefulness of near-infrared reflectance (NIR) spectroscopy and chemometrics to discriminate fishmeal batches made with different fish species,” J. Agric. Food Chem. 53(11), 4459–4463 (2005). [CrossRef]   [PubMed]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (8)

Fig. 1
Fig. 1 The Raman spectra of blood samples with and without spectral shift.
Fig. 2
Fig. 2 The modules of the Fourier transform of the Raman spectra in Fig. 1.
Fig. 3
Fig. 3 The prediction value. (a) is obtained by PLS when the Raman spectrum was shifted the pixel mathematically; (b) is obtained by FPLS when the Raman spectrum was shifted the pixel mathematically.
Fig. 4
Fig. 4 The prediction value obtained by PLS when the Raman spectra have a spectral shift. (a) shown prediction value of the Raman spectra without spectral shift; (b)-(e) were corresponding to prediction value of the Raman spectra with a shift of 2.82 cm−1, 6.36 cm−1, 9.18 cm−1 and 19.77 cm−1.
Fig. 5
Fig. 5 The average prediction value with different spectral shift when PLS model was used.
Fig. 6
Fig. 6 The prediction value obtained by FPLS when the Raman spectra have a spectral shift. (a) shown prediction value of the Raman spectra without spectral shift; (b)-(e) were corresponding to prediction value of the Raman spectra with a shift of 2.82 cm−1, 6.36 cm−1, 9.18 cm−1 and 19.77 cm−1
Fig. 7
Fig. 7 The average prediction value with different spectral shift when FPLS model was used
Fig. 8
Fig. 8 The prediction value. (a) is obtained by PLS when the Raman spectra in the training data set have spectral shift; (b) is obtained by FPLS model when the Raman spectra in the training data set have spectral shift

Equations (10)

Equations on this page are rendered with MathJax. Learn more.

k = 1 λ e x c i t a t i o n 1 λ
Δ y = ( a ( k 1 ) I ' ( k 1 ) + a ( k 2 ) I ' ( k 2 ) + ... + a ( k n ) I ' ( k n ) ) δ k + ( a ( k 1 ) I ' ( k 1 ) 2 ! + a ( k 2 ) I ' ( k 2 ) 2 ! + ... + a ( k n ) I ' ( k n ) 2 ! ) δ k 2 + ... + ( a ( k 1 ) I ' ( k 1 ) n ! + a ( k 2 ) I ' ( k 2 ) n ! + ... + a ( k n ) I ' ( k n ) n ! ) δ k n
I s a m p l e d ( k ) = I ( k ) Π δ ξ / 2 ( ξ )
F F T ( I s a m p l e d ( k ) ) = F F T ( I ( k ) Π δ ξ / 2 ( ξ ) ) = F F T ( I ( k ) ) × F F T ( Π δ ξ / 2 ( ξ ) ) = G ( ω ) × sin c ( ω δ ξ / 2 )
| F F T ( I s a m p l e d ( k ) ) | = | G ( ω ) × sin c ( ω δ ξ / 2 ) | = A ( ω ) × | sin c ( ω δ ξ / 2 ) |
| F F T ( I s a m p l e d ( k + δ k ) ) | = | F F T ( I ( k + δ k ) Π δ ξ / 2 ( ξ ) ) | = | F F T ( I ( k + δ k ) ) × F F T ( Π δ ξ / 2 ( ξ ) ) | = | G ( ω ) exp ( j ω δ k ) × sin c ( ω δ ξ / 2 ) | = A ( ω ) × | exp ( j ω δ k ) | × | sin c ( ω δ ξ / 2 ) | = A ( ω ) × | sin c ( ω δ ξ / 2 ) |
| F F T ( I ( k ) ) | = | F F T ( I ( k + δ k ) ) | = A ( ω ) × | sin c ( ω δ ξ / 2 ) |
I s a m p l e d ( k ) = i = 1 n I s a m p l e d ( k i )
| F F T ( I s a m p l e d ( k ) ) | = | i = 1 n G i ( ω ) × sin c ( ω δ ξ / 2 ) | = | i = 1 n G i ( ω ) | × | sin c ( ω δ ξ / 2 ) |
| F F T ( I s a m p l e d ' ( k ) ) | = | F F T ( i = 1 n I s a m p l e d ( k i + δ k i ) Π δ ξ / 2 ( ξ ) ) | = | i = 1 n ( G i ( ω ) exp ( j ω δ k i ) ) × sin c ( ω δ ξ / 2 ) | = | i = 1 n ( G i ( ω ) exp ( j ω δ k i ) ) | × | sin c ( ω δ ξ / 2 ) |
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.