Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
An Edge Computing Based Smart Healthcare Framework for Resource Management
Previous Article in Journal
Implementation and Analysis of Tightly Coupled Global Navigation Satellite System Precise Point Positioning/Inertial Navigation System (GNSS PPP/INS) with Insufficient Satellites for Land Vehicle Navigation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lifting Wavelet Transform De-noising for Model Optimization of Vis-NIR Spectroscopy to Predict Wood Tracheid Length in Trees

1
College of Engineering and Technology, Northeast Forestry University, Harbin 150040, China
2
Forest Products Development Center, SFWS, Auburn University, Auburn, AL 36849, USA
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(12), 4306; https://doi.org/10.3390/s18124306
Submission received: 25 October 2018 / Revised: 26 November 2018 / Accepted: 4 December 2018 / Published: 6 December 2018
(This article belongs to the Section Remote Sensors)

Abstract

:
The data analysis of visible-near infrared (Vis-NIR) spectroscopy is critical for precise information extraction and prediction of fiber morphology. The objectives of this study were to discuss the de-noising of Vis-NIR spectra, taken from wood, to improve the prediction accuracy of tracheid length in Dahurian larch wood. Methods based on lifting wavelet transform (LWT) and local correlation maximization (LCM) algorithms were developed for optimal de-noising parameters and partial least squares (PLS) was employed as the prediction method. The results showed that: (1) The values of tracheid length in the study were generally high and had a great positive linear correlation with annual rings (R = 0.881), (2) the optimal de-noising parameters for larch wood based Vis-NIR spectra were Daubechies-2 (db2) mother wavelet with 4 decomposition levels while using a global fixed hard threshold based on LWT, and (3) the Vis-NIR model based on the optimal LWT de-noising parameters ( R c 2 = 0.834, RMSEC = 0.262, RPD c = 2.454) outperformed those based on the LWT coupled with LCM algorithm (LWT-LCM) ( R c 2 = 0.816, RMSEC = 0.276, RPD c = 2.331) and raw spectra ( R c 2 = 0.822, RMSEC = 0.271, RPD c = 2.370). Thus, the selection of appropriate LWT de-noising parameters could aid in extracting a useful signal for better prediction accuracy of tracheid length.

1. Introduction

Visible and near-infrared (Vis-NIR) spectroscopy has been widely applied to agriculture, life science, polymer materials, and other fields [1,2,3] because of its many advantages, including simplicity of operation, no need for sample pretreatment, real-time online monitoring, and great adaptability of the model. It has been shown that Vis-NIR spectroscopy can be used in wood science for the prediction of wood’s physical, chemical, anatomical, and mechanical properties such as moisture content, microfibril angle, bending strength, etc. [4,5,6,7].
As an indirect and nondestructive method, Vis-NIR spectroscopy has the main bands of spectra that correspond to the stretching vibrations of C-H, N-H, O-H, and S-H in almost all biological samples [8]. Therefore, Vis-NIR could well characterize the major properties or compositions of the corresponding samples. Prediction of samples’ properties using Vis-NIR spectroscopy depends on chemometric analysis that explains the relationship between spectra and the corresponding properties. Both linear and nonlinear methods are used successfully for modeling algorithm, such as partial least squares (PLS) regression, support vector machine (SVM), least squares-support vector machine (LS-SVM), and artificial neural network (ANN) [9,10].
Wood, a complex natural organic compound, is mainly composed of cellulose, hemicellulose, and lignin. Wood microstructure parameters mainly include the microfibril angle, fiber length, tracheid length, crystallinity, or any fiber morphology that influence the wood physical, chemical, or mechanical properties [11]. These parameters are of importance for the forest genetics, trees-oriented breeding or the processing and use of wood. However, compared to the traditional methods of physical and chemical analysis, the determination of anatomical features is more complex, time-consuming, and places higher requirements upon forestry researchers. These determination methods cannot be used to realize the assessment of wood traits on a large-scale, multispecies and multisite.
Vis-NIR spectroscopy presents wood composition information through the interaction between Vis-NIR light and internal molecules. The goal for the prediction of wood properties is to develop a calibration model that relates Vis-NIR spectra to cell wall composition and morphology with multivariate statistics. Therefore, spectral quality is an important factor affecting the prediction precision for Vis-NIR models. However, environmental factors, spectral devices, surface roughness, and other various factors induce noise in spectra during acquisition [12,13], such as optical noise, electrical noise, and their interactions, which reduce prediction accuracy and precision. Determining how to remove noise while retaining as much useful information as possible is a challenge for forest researchers.
In recent years, spectral pretreatments including derivatives, multiplication scatter correction (MSC), standard normalized variate (SNV), Savitzky-Golay (SG) smoothing and wavelet transform (WT) [14,15,16] have been used widely for removing noise or other irrelevant information from the signal to improve model performance. While WT is superior in de-noising and compressing Vis-NIR data, removing the noise without losing useful information is still a challenge because the noise in the Vis-NIR spectra is mostly a random signal. Additionally, the traditional WT is time-consuming, and sometimes a large amount of calculations with advanced software is needed. In light of the above deficiencies, as the second-generation wavelet transform, lifting wavelet transform (LWT), with its advantages of rapid computing speed, simple algorithms, and lack of dependence on Fourier transform can be directly refactored into traditional wavelets in time domain and overcome the weakness of the boundary problem [17].
Dahurian larch is one of the major commercial tree species in northeastern China. Tracheid length could be a good indication for pulping properties and tensile strength of wood. Hence, the rapid determination of tracheid length is the key for timber quality classification and evaluation for fiber. In this study, Dahurian larch wood samples were used for Vis-NIR spectra collection and tracheid length estimation. Spectral data were processed by lifting wavelet transform (LWT) and the combination of local correlation maximization (LCM) algorithm. Wavelet transform (WT), moving average, loess, Savitzky-Golay and lowess were compared in this study, as well. The preprocessed Vis-NIR data were used for model development with the PLS method, which provides a theoretical basis and technical support for the rapid and accurate de-noising of wood Vis-NIR spectra.

2. Materials and Methods

2.1. Sample Preparation

Seven Dahurian larch (Larix gmelinii) trees ranging between 39 and 43 years in age were harvested from a larch plantation located at 45°44′00″–45°53′02″ N and 131°8′40″–131°21′23″ E, in Heilongjiang Province, China. Tree heights ranged from 18.8 to 22.5 m, with an average of 21.0 m and the crown width ranged from 2.8 to 4.6 m. The trees were removed from different parts to better represent the range of natural variability in the region. Five-centimeter-thick (5 cm) disks were removed from each tree near breast height (1.3 m) with a total of seven discs (Figure 1) used for model calibration. The disks were left to air dry in an environmental-controlled laboratory for three months. After removing bark from the edge of each disk, slab samples from pith to bark 2 × 2 × 5 cm3 in tangential, radial, and longitudinal dimensions were prepared for Vis-NIR spectra collection and tracheid length measurement with a total of 83 slab samples.

2.2. Vis-NIR Spectroscopy

The Vis-NIR spectra were collected from the rings of each sample using a LabSpec Pro FR/A114260 (Analytical Spectral Devices, Inc., Boulder, CO, USA). Thirty scans were obtained at 1 nm intervals, averaged into one spectrum, and saved as log (1/R). Because the ring width in the mature wood zone was smaller than the spot size of the fiber-optic probe, which was approximately 5 mm, the spectra from latewood portion were obtained [18]. Additionally, to reduce the noise from the spectrometer system, spectra data were collected for each ring from the 1st to 20th annual ring and in two ring intervals after the 21st annual ring. Wavelength range was limited to 350–2397 nm (wavelengths = 2048).

2.3. Tracheid Length Measurement and Calibrations

After Vis-NIR spectra were collected, the tracheid lengths of the latewood were detected by the isolation method [19], macerated with Jeffrey’s solution, and measured using a universal projector. Each tracheid length was recorded by averaging 25 measurements per sample.
The partial least squares (PLS) regression was conducted in the multivariate statistical analysis software called Unscrambler V10.4 (CAMO Software AS, Oslo, Norway). The models’ performance was assessed with the following prediction diagnostics [20], including the coefficient of determination (R2), root mean square error (RMSE), mean absolute percentage error (MAPE), standard error of calibration in the laboratory (SEC), and ratio of performance to standard deviation (RPD). The computation equations of these criteria are shown as follows:
R 2 = 1 i = 1 n ( y i y i ^ ) 2 / i = 1 n ( y i y ¯ ) 2
RMSE = i = 1 n ( y i y i ^ ) 2 / n
MAPE = 1 n i = 1 n | ( y i y i ^ ) / y i | × 100 %
SEC = i = 1 n ( y i y i ^ 1 n i = 1 n ( y i y i ^ ) ) 2 / n 1
RPD = SD / i = 1 n ( y i y i ^ 1 n i = 1 n ( y i y i ^ ) ) 2 / n 1
where y i represents the tracheid length value, y i ^ and y ¯ are the predicted value and the mean of y i , respectively. n is the number of samples. When n is the number of samples of the calibration set, coefficient of determination, ratio of performance to standard deviation and error are named R c 2 , RPD c , RMSEC and MAPE c , respectively. When n is the number of samples of the prediction set, coefficient of determination, ratio of performance to standard deviation and error are named R p 2 , RPD p and RMSEP , respectively.

2.4. Vis-NIR Data Processing

2.4.1. Lifting Wavelet Transform Analysis

The lifting wavelet transform (LWT), a second-generation wavelet, inherits the advantages of the multi-resolution of the traditional wavelet transform and has been recognized as a powerful signal processing tool with fast speed and small memory requirements. The basic steps of LWT are split, predict, and update [21,22]. The aim of LWT is to process lifting transform coefficients obtained by the decomposition of spectra based on the lifting scheme using the fixed threshold and then the reconstruction of the Vis-NIR spectra.
In this study, in order to ensure optimal de-noising parameters, four different mother wavelets (symN, biorN, rbioN, and dbN) were first used with the hypotheses of 5 decomposition layer (k = 5) and order (N) of 5 (N = 5.5 for bior and rbio wavelet) for the optimal wavelet, and then the spectral data were pretreated with a different order. After the optimal order was determined, the processed spectra were analyzed by different decomposition levels.
The LWT was implemented using MatlabR2014b (MathWorks, Natick, MA, USA), and the R c 2 , RMSEC , and MAPE c were used for measuring the optimal de-noising parameters including mother wavelet, decomposition layer (k), and order (N). Generally, higher R2 and lower RMSE and MAPE indicate a better prediction result [23].

2.4.2. Local Correlation Maximization Algorithm

In view of the phenomenon that the noise removal of some wavelengths is insufficient or eliminates useful signals, the local correlation maximization (LCM) algorithm was operated under the optimal LWT de-noising parameters. It was used for construction of the spectra through selecting the absorbance with the highest correlation between decomposition levels (one to the optimal level) and tracheid length in the wavelength range of 350–2397 nm. The specific steps are as follows [24], and it was implemented by using Matlab R2014b.
  • Correlation analysis: the correlation coefficient (r) between Vis-NIR spectra under different decomposition levels and tracheid length were obtained.
  • Judgment analysis: for the wavelength range of 350–2397 nm, each wavelength corresponds to multiple correlation coefficients, and the decomposition layer with the largest r in all decomposition layers was selected as the decomposition layer for this wavelength.
  • Construction of spectra: the spectra were constructed by the absorbance of the decomposition level with the largest r.

2.4.3. Comparison with Basic De-Noising Methods

Environmental factors, spectrometer, human operation, and other various factors induce noise in Vis-NIR spectra, such as baseline drift and light scattering. To eliminate the influence of this noise and analyze the effect of LWT and LCM, wavelet transform (WT), moving average, loess, Savitzky-Golay and lowess were compared. The parameters of wavelet de-noising were selected according to the optimal parameters of LWT. Additionally, the segment size for moving average, loess, Savitzky-Golay, and lowess were discussed.

2.5. Overview of Vis-NIR Spectra De-Noising Processing

A total of 164 samples were analyzed, which included a calibration set with 117 samples and a validation set with 47 samples. Firstly, Vis-NIR spectra of the calibration set were pretreated by four mother wavelets (sym5, bior5.5, rbio5.5 and db5) with five decomposition levels to ensure the optimal wavelet. The optimal order (N) and decomposition level (k) were then determined under the optimal wavelet. After the optimal LWT de-noising parameters were obtained, LCM algorithm, wavelet transform (WT), moving average, loess, Savitzky-Golay, and lowess were used for further de-noising. Finally, the processed spectral data were modeled by PLS. The total treatment is depicted in Figure 2.

3. Results and Analysis

3.1. Statistical Characteristics of Wood Tracheid Length

The total samples were randomly divided into a calibration set (117 samples) and a prediction set (47 samples). As illustrated in Table 1, the tracheid length ranged from 1.626 mm to 4.618 mm, with an average value of 3.320 mm, indicating that the larch tracheid length in this area was generally large relative to the literature. The distribution exhibited low negative skewness and kurtosis, indicating low scatter distribution.

3.2. Radial Development of Wood Tracheid Length in Annual Rings

A strong positive linear correlation between tracheid length and annual ring was found as shown in Figure 3, and the correlation coefficient was larger than 0.80. The tracheid length near the pith was relatively small and increased rapidly before 18 years and levelled off between the 23rd and 35th annual rings. The tracheid length reached the maximum value at the 40 annual rings.

3.3. Selection of Optimal LWT De-Noising Parameters

As the first step of LWT de-noising, the choice of mother wavelets greatly influences the effects of de-noising, but there is no definitive guideline on how to select appropriate wavelets, so flexibility in selecting parameters is needed. In this study, symN, biorN, rbioN, and dbN wavelets were analyzed, and the assumed decomposition levels and order are all 5 (k = N = 5) under the global fixed threshold.
As shown in Table 2, the modeling results from bior5.5 and rbio5.5 were inferior to those from db5 and sym5, while the spectral data processed by db5 with the fifth decomposition level under the global fixed hard threshold had a higher accuracy, which may be due to the property of the biorthogonal compactly supported for the db wavelet [25]. Thus, the Vis-NIR spectra pretreated by the dbN wavelet family with fifth decomposition level under the global fixed hard threshold were used for selecting optimal order (N).
It is difficult to apply the dbN wavelet family without selecting an appropriate choice of order. To effectively remove noise, the effect of de-noising of different orders, (N = 1–8) for the db wavelet, are demonstrated in Table 3, and it can be seen that the results were influenced by wavelet orders. Among these models with dbN preprocessing, the performance of the db8 wavelet was the worst for the data set in which the R c 2 was less than 0.50. The model with the db2 wavelet obtained the best results with a higher R c 2 and lower RMSEC and MAPE c ( R c 2 = 0.818, RMSEC = 0.274, MAPE c = 7.443). Therefore, the optimal order for Vis-NIR spectra of larch wood is 2 (N = 2), and the spectra pretreated with db2 was used for further analysis.
It is a challenge to select appropriate decomposition levels for wavelet de-noising, but the levels between second to fifth were randomly selected in most studies [26]. After the optimal wavelet and order were selected, decomposition levels from first to fifth were analyzed.
As shown in Figure 4, the R c 2 increased gradually first and then declined with increasing decomposition levels, and similar conclusions were also deduced by Cai [27] in the study on soil moisture content (SMC). Among models with first to fifth levels of the db2 wavelet, the optimal decomposition level is fourth because the R c 2 was relatively higher than that of other models indicating that the noise or irrelative signal was removed. Overall, the optimal de-noising parameters of LWT for larch Vis-NIR spectra is db2 with fourth level (k = 4, N = 2).

3.4. Comparison with Basic De-Noising Methods

After the optimal parameters were obtained, moving average, loess, Savitzky-Golay, and lowess were used for further de-noising analysis. The selection of an appropriate segment size parameter is a key point. Figure 5 shows the de-noising effects of different segment size for moving average, loess, Savitzky-Golay, and lowess.
As seen in Figure 5, R c 2 and RPD c decreased gradually when segment size increased from 3 to 25 for the four pretreatment methods, while RMSEC and SEC increased with increase of segment size. For pretreatments, the best performance was obtained with three segment size. Except for moving average ( R c 2 = 0.809 ), the R c 2 of the best model was the same ( R c 2 = 0.822 ). The performance of the model with 25 segment size was the worst, while loess obtained the worst performance with the segment size of 25. However, the R c 2 of the worst models were all larger than 0.78.

3.5. Establishment of Vis-NIR Models for Wood Tracheid Length

The spectral data pretreated by LWT, LCM algorithm, and wavelet transform (db2 with fourth level) were then modeled by PLS. In terms of the results (Table 4), compared to the raw model, the modeling results pretreated by LWT-LCM and WT were inferior to. However, the performances of calibration and validation models were improved after the LWT pretreatment. As for the calibration model, the R c 2 and RPD c were improved by 2% and 4%, respectively, while RMSEC and SEC were decreased by 3% and 3%, respectively. Additionally, the performance of LWT was better than loess, Savitzky-Golay, and lowess.
Even though WT had the same R c 2 as LWT-LCM, the spectra of LWT-LCM model were constructed by the absorbance of the decomposition level with the largest r. Therefore, the prediction accuracies of the other three models for prediction sets with various pretreatments were analyzed. As shown in Figure 6, regardless of various preprocessing methods for the prediction sets, the model whose spectral data were processed by LWT exhibited a good fit for three prediction sets with an R p 2 larger than 0.700. However, the performance of LWT-LCM coupled model was lower, which may be attributed to more error in the calibration model. When different prediction sets were used, the following three models obtained better results in untreated prediction set than other prediction sets. For the LWT-LCM coupled model, the R p 2 was higher for prediction set with LWT preprocessing than that for LWT-LCM preprocessing. However, the situation was reversed for the raw model and LWT model.

4. Discussion

This study demonstrated that the tracheid length had a positive linear correlation with annual ring, which is consistent with the results of Boruszewski et al. [28]. For different pretreatment methods, the spectral data pretreated with the optimal LWT de-noising parameters (db2 wavelet with the fourth decomposition level under the global fixed hard threshold) can be used to better remove noise and extract characteristic signal, and thus obtain better prediction results for tracheid length ( R c 2 = 0.834) than raw model ( R c 2 = 0.822). With comparisons of previous studies on the prediction of tracheid length based on NIR technology, the accuracy of model in this study is inferior to the results in Schimleck et al. [29] ( R c 2 = 0.88), even though the samples are the same (Pinus), different geography and environment could induce different results. Additionally, different wavelengths for calibration were used. For the prediction of tracheid length in other species, although the accuracy of raw model in this study is a little lower than the results of Norway spruce [30] (0.828), the results of LWT pretreatment is better than spruce, which demonstrated that the feasibility of de-noising based on LWT.
The ability to effectively predict tracheid length was due to the change of Vis-NIR spectra with annual ring and there was a significant correlation between ring, chemical properties (lignin and cellulose as associated with 1415, 1965, and 2315 nm), and tracheid length [31].
However, after LWT coupled with the LCM algorithm pretreatment, model accuracy was lower than the raw model, and Figure 6 shows that the R p 2 of the LWT-LCM model was lower than others (excluding the prediction set with LWT). This could be caused by the unnecessary wavelengths embedded in the reorganized data matrix which decrease the correlation in the characteristic band and results in more prediction errors. To further verify this speculation, correlation analysis between tracheid length and corresponding pretreatment spectra were studied.
As shown in Figure 7a, there were four main peaks (maximum value of correlation coefficient) at around 438 nm, 1114 nm, 1306 nm, and 1670 nm, respectively, and four main valleys (minimum value of correlation coefficient) were present near 734 nm, 1212 nm, 1470 nm, and 1944 nm, respectively. The correlation coefficient curves for raw spectra and spectra pretreated by LWT exhibited a uniform pattern, and there was little change compared to that of the raw spectra. However, in terms of the results for spectra pretreated by the LWT-LCM algorithm (Figure 7b), in the range of 1409–2397 nm, the unstable correlation coefficient results in a large fluctuation. Additionally, in the wavelengths of peaks and valleys, the correlation of spectra pretreated by the LWT-LCM algorithm was slightly lower than that of the raw spectra. This indicates that the speculation is reasonable.

5. Conclusions

The tracheid length values were generally high and exhibited a strong positive linear correlation with annual rings (ring from pith). The modeling results with different de-noising methods were different. The optimal parameters for larch wood spectra were db2 wavelet with decomposition layer four based on the lifting scheme. The LWT model with optimal de-noising parameters had relatively better accuracy and precision. However, the LWT-LCM model was not a good fit, which may be due to the unnecessary wavelengths embedded in the reorganized data matrix that decrease the correlation in characteristic bands and result in more prediction errors. Regarding models with different pretreatments, better prediction results were obtained for prediction sets with raw spectra than those of other sets. This study mainly discussed the de-noising of wood Vis-NIR spectra for improving the prediction accuracy of wood tracheid length based on different de-noising methods. In a follow-up study, more applications for prediction of wood’s properties based on smoothed data will be explored to provide a theoretical basis and technical support for the optimization of wood Vis-NIR spectra.

Author Contributions

Conceptualization, Y.L. and X.L.; methodology, B.K.V., Z.C., X.L.; software, B.K.V., Z.C. and X.L.; validation, Z.C., X.L.; formal analysis, B.K.V.; investigation, Y.L., X.L.; resources, Y.L., X.L.; data curation, Y.L.; writing—original draft preparation, Y.L.; writing—review and editing, B.K.V., Z.C. and X.L.; visualization, Y.L., B.K.V., Z.C. and X.L.; supervision, Y.L., B.K.V., Z.C. and X.L.; project administration, Y.L., B.K.V., Z.C. and X.L.; funding acquisition, X.L.

Funding

This research was funded by the China National Key Research and Development Program, grant number 2017YFC0504103 and the Fundamental Research Funds for Central Universities, grant number 2572017AB04. And the APC was funded by Yaoxiang Li.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Miloš, B.; Bensa, A. Estimation of calcium carbonate in anthropogenic soils on flysch deposits from Dalmatia (Croatia) using Vis-NIR spectroscopy. Agric. For. 2017, 63, 121–130. [Google Scholar] [CrossRef]
  2. Verbeek, F.P.R.; Troyan, S.L.; Mieog, J.S.; Liefers, G.J.; Moffitt, L.A.; Rosenberg, M.; Hirshfield-Bartek, J.; Gioux, S.; Velde, C.J.H.; Vahrmeijer, A.L.; et al. Near-infrared fluorescence sentinel lymph node mapping in breast cancer: A multicenter experience. Breast Cancer Res. Treat. 2014, 143, 333–342. [Google Scholar] [CrossRef] [PubMed]
  3. Duconseille, A.; Andueza, D.; Picard, F.; Santé-Lhoutellier, V.; Astruc, T. Molecular changes in gelatin aging observed by NIR and fluorescence spectroscopy. Food Hydrocolloids 2016, 61, 496–503. [Google Scholar] [CrossRef]
  4. Zanuncio, A.J.V.; Hein, P.R.G.; Carvalho, A.G.; Rocha, M.F.V.; Carneiro, A.C.O. Determination of heat-treated eucalyptus and pinus wood properties using NIR spectroscopy. J. Trop. For. Sci. 2018, 30, 117–125. [Google Scholar] [CrossRef]
  5. Via, B.K.; Zhou, C.F.; Acquah, G.; Jiang, W.; Eckhardt, L. Near infrared spectroscopy calibration for wood chemistry: Which chemometric technique is best for prediction and interpretation. Sensors Basel 2014, 14, 13532–13547. [Google Scholar] [CrossRef] [PubMed]
  6. Schimleck, L.R.; Matos, J.L.M.; Trianoski, R.; Prata, J.G. Comparison of methods for estimating mechanical properties of wood by NIR spectroscopy. J. Spectrosc. 2018, 2018, 1–10. [Google Scholar] [CrossRef]
  7. Hein, P.R.G.; Pakkanen, H.K.; Dos Santos, A.A. Challenges in the use of near infrared spectroscopy for improving wood quality: A review. For. Syst. 2017, 26, eR03. [Google Scholar] [CrossRef]
  8. Wu, Z.S.; Ouyang, G.Q.; Shi, X.Y.; Ma, Q.; Wan, G.; Qiao, Y.J. Absorption and quantitative characteristics of C-H bond and O-H bond of NIR. Opt. Spectrosc. 2014, 117, 703–709. [Google Scholar] [CrossRef]
  9. Tenenhaus, M.; Esposito Vinzi, V.; Chatelinc, Y.M.; Lauro, C. PLS path modeling. Comput. Stat. Data Anal. 2005, 48, 159–205. [Google Scholar] [CrossRef]
  10. Balabin, R.M.; Lomakina, E.I.; Safieva, R.Z. Neural network (ANN) approach to biodiesel analysis: Analysis of biodiesel density, kinematic viscosity, methanol and water contents using near infrared (NIR) spectroscopy. Fuel 2011, 90, 2007–2015. [Google Scholar] [CrossRef]
  11. Cheng, Q.; Zhou, C.F.; Jiang, W.; Zhao, X.P.; Via, B.; Wan, H. Mechanical and physical properties of OSB exposed to high temperature and relative humidity and coupled with NIR modeling. For. Prod. J. 2018, 68, 78–85. [Google Scholar] [CrossRef]
  12. Kurata, Y. Nondestructive classification analysis of wood soaked in seawater by using near-infrared spectroscopy. For. Prod. J. 2017, 67, 63–68. [Google Scholar] [CrossRef]
  13. Rossel, R.A.; Webster, R. Predicting soil properties from the Australian soil visible-near infrared spectroscopic database. Eur. J. Soil Sci. 2012, 63, 848–860. [Google Scholar] [CrossRef]
  14. Li, Y.; Li, Y.X.; Li, W.B.; Jiang, L.C. Model optimization of wood property and quality tracing based on wavelet transform and NIR spectroscopy. Spectrosc. Spect. Anal. 2018, 38, 1384–1392. [Google Scholar] [CrossRef]
  15. Liu, J.L.; Sun, B.L.; Yang, Z. Estimation of the physical and mechanical properties of Neosinocalamus Affinins using near infrared spectroscopy. Spectrosc. Spect. Anal. 2011, 31, 647–651. [Google Scholar] [CrossRef]
  16. Li, Y.; Li, Y.X.; Xu, H.K.; Jiang, L.C. Model calibrating for NIRS-based oak wood air-dry density prediction with de-noising pretreatment. J. Nanjing For. Univ. 2016, 40, 148–156. [Google Scholar] [CrossRef]
  17. Yang, X.F.; Zhang, W.; Yang, Y.X. Denoising technology of radar life signal based on lifting wavelet transform. Acta Opt. Sin. 2014, 34, 300–305. [Google Scholar] [CrossRef]
  18. Via, B.K.; So, C.L.; Groom, L.H.; Shupe, T.F.; Stine, M.; Wikaira, J. Within tree variation of lignin, extractives, and microfibril angle coupled with the theoretical and near infrared modeling of microfibril angle. IAWA J. 2007, 28, 189–209. [Google Scholar] [CrossRef]
  19. Sadegh, A.N.; Kiaei, M. Formation of juvenile/mature wood in Pinus eldarica medw and related wood properties. World Appl. Sci. J. 2011, 12, 460–464. [Google Scholar]
  20. Zhang, X.C.; Wu, J.Z.; Xu, Y. Near Infrared Spectroscopy Technique and Applications in Modern Agriculture; Publishing House of Electronics Industry: Beijing, China, 2012; pp. 48–49. [Google Scholar]
  21. Wang, Y.T.; Yang, Z.; Hou, P.G.; Cheng, P.F.; Cao, L.F. Research on spectroscopy spectrum de-noising of mineral oil based on lifting scheme wavelet transform. Spectrosc. Spect. Anal. 2016, 36, 2144–2147. [Google Scholar] [CrossRef]
  22. Zhou, F.B.; Li, C.G.; Zhu, H.Q. Research on threshold improved denoising algorithm based on lifting wavelet transform in UV-Vis spectrum. Spectrosc. Spect. Anal. 2018, 38, 506–510. [Google Scholar] [CrossRef]
  23. Yan, Y.L.; Chen, B.; Zhu, D.Z. Near Infrared Spectroscopy-Principle/Technology and Application; China Light Industry Press: Beijing, China, 2013; pp. 172–174. [Google Scholar]
  24. Zhang, R.; Li, Z.F.; Pan, J.J. Coupling discrete wavelet packet transformation and local correlation maximization improving prediction accuracy of soil organic carbon based on hyperspectral reflectance. Trans. Chin. Soc. Agric. Eng. 2017, 33, 175–181. [Google Scholar] [CrossRef]
  25. Boggess, A.; Narcowich, F.J. A First Course in Wavelets with Fourier Analysis, 2nd ed.; Publishing House of Electronics Industry: Beijing, China, 2010; pp. 211–217. [Google Scholar]
  26. Srivastava, M.; Lindsay Anderson, C.; Freed, J.H. A new wavelet denoising method for selecting decomposition levels and noise thresholds. IEEE Access 2016, 4, 3862–3877. [Google Scholar] [CrossRef] [PubMed]
  27. Cai, L.H.; Ding, J.L. Inversion of soil moisture content based on hyperspectral multi-scale decomposition. Laser Optoelectron. Progress 2018, 55, 013001. [Google Scholar] [CrossRef]
  28. Boruszewski, P.; Jankowska, A.; Kurowska, A. Comparison of the structure of juvenile and mature wood of Larix decidua mill from fast-growing plantations in Poland. Bioresources 2017, 12, 1813–1825. [Google Scholar] [CrossRef]
  29. Schimleck, L.R.; Jones, P.D.; Peter, G.F.; Daniels, R.F.; ClarkllI, A. Nondestructive estimation of tracheid length from sections of radial wood strips by near infrared spectroscopy. Holzforschung 2004, 58, 375–381. [Google Scholar] [CrossRef]
  30. Mäkinen, H.; Jaakkola, T.; Piispanen, R.; Saranpää, P. Predicting wood and tracheid properties of Norway spruce. For. Ecol. Manag. 2007, 24, 1175–1188. [Google Scholar] [CrossRef]
  31. Via, B.K.; Shupe, T.F.; Stine, M.; So, C.L.; Groom, L.H. Tracheid length prediction in Pinus palustris by means of near infrared spectroscopy: The influence of age. Eur. J. Wood Prod. 2005, 63, 231–236. [Google Scholar] [CrossRef]
Figure 1. Dahurian larch disc sample.
Figure 1. Dahurian larch disc sample.
Sensors 18 04306 g001
Figure 2. An overview of steps in Vis-NIR spectra processing.
Figure 2. An overview of steps in Vis-NIR spectra processing.
Sensors 18 04306 g002
Figure 3. Radial variation in larch tracheid length. Data are averages (±SD) for the ring from 1 to 40.
Figure 3. Radial variation in larch tracheid length. Data are averages (±SD) for the ring from 1 to 40.
Sensors 18 04306 g003
Figure 4. Results of PLS model for various decomposition levels of db2 wavelet.
Figure 4. Results of PLS model for various decomposition levels of db2 wavelet.
Sensors 18 04306 g004
Figure 5. Results of PLS model for various segment sizes. MA: moving average; SG: Savitzky-Golay.
Figure 5. Results of PLS model for various segment sizes. MA: moving average; SG: Savitzky-Golay.
Sensors 18 04306 g005
Figure 6. Predicted versus measured wood tracheid length for various prediction sets.
Figure 6. Predicted versus measured wood tracheid length for various prediction sets.
Sensors 18 04306 g006
Figure 7. Correlation between corresponding spectra and tracheid length values. (a) The comparison of correlation coefficient between raw spectra and spectra pretreated by LWT; (b) the comparison of correlation coefficient between raw spectra and spectra pretreated by LWT-LCM.
Figure 7. Correlation between corresponding spectra and tracheid length values. (a) The comparison of correlation coefficient between raw spectra and spectra pretreated by LWT; (b) the comparison of correlation coefficient between raw spectra and spectra pretreated by LWT-LCM.
Sensors 18 04306 g007
Table 1. Summary statistics on tracheid length of larch in each data set.
Table 1. Summary statistics on tracheid length of larch in each data set.
Sample SetNo. of SamplesMax (mm)Min (mm)Avg. (mm)SD (mm)SkewnessKurtosis
Calibration Set1174.1691.6263.2340.645−0.715−0.539
Prediction Set474.6181.9513.5340.678−0.614−0.285
Total1644.6181.6263.3200.667−0.603−0.421
Table 2. Results of PLS model of various wavelet with fifth decomposition level.
Table 2. Results of PLS model of various wavelet with fifth decomposition level.
WaveletPCs R c 2 RMSECMAPEc (%)
sym570.7830.3008.185
bior5.550.3950.50013.460
rbio5.570.5030.45312.407
db570.8110.2797.639
Table 3. Results of PLS model of dbN wavelet family with fifth decomposition level.
Table 3. Results of PLS model of dbN wavelet family with fifth decomposition level.
dbN R c 2 RMSECMAPEc (%)
db10.8070.2827.782
db20.8180.2747.443
db30.7630.3138.738
db40.8090.2817.839
db50.8110.2797.639
db60.6000.40610.971
db70.7890.2958.304
db80.4440.47913.384
Table 4. Model statistics of wood tracheid length. Raw: raw model; LWT: model with LWT pretreatment; LWT-LCM: model with LWT coupled with LCM pretreatment; WT: model with WT pretreatment.
Table 4. Model statistics of wood tracheid length. Raw: raw model; LWT: model with LWT pretreatment; LWT-LCM: model with LWT coupled with LCM pretreatment; WT: model with WT pretreatment.
ModelPCsCalibration SetValidation Set
R c 2 RMSECSEC RPD c R p 2 RMSEPSEP RPD p
Raw70.8220.2710.2722.3700.7140.3470.3491.870
LWT70.8340.2620.2632.4540.7220.3440.3451.897
LWT-LCM70.8160.2760.2772.3310.6830.3650.3671.776
WT70.8160.2750.2772.3310.7170.3470.3461.880

Share and Cite

MDPI and ACS Style

Li, Y.; Via, B.K.; Cheng, Q.; Li, Y. Lifting Wavelet Transform De-noising for Model Optimization of Vis-NIR Spectroscopy to Predict Wood Tracheid Length in Trees. Sensors 2018, 18, 4306. https://doi.org/10.3390/s18124306

AMA Style

Li Y, Via BK, Cheng Q, Li Y. Lifting Wavelet Transform De-noising for Model Optimization of Vis-NIR Spectroscopy to Predict Wood Tracheid Length in Trees. Sensors. 2018; 18(12):4306. https://doi.org/10.3390/s18124306

Chicago/Turabian Style

Li, Ying, Brian K. Via, Qingzheng Cheng, and Yaoxiang Li. 2018. "Lifting Wavelet Transform De-noising for Model Optimization of Vis-NIR Spectroscopy to Predict Wood Tracheid Length in Trees" Sensors 18, no. 12: 4306. https://doi.org/10.3390/s18124306

APA Style

Li, Y., Via, B. K., Cheng, Q., & Li, Y. (2018). Lifting Wavelet Transform De-noising for Model Optimization of Vis-NIR Spectroscopy to Predict Wood Tracheid Length in Trees. Sensors, 18(12), 4306. https://doi.org/10.3390/s18124306

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop