Next Article in Journal
Tartary Buckwheat Starch Modified with Octenyl Succinic Anhydride for Stabilization of Pickering Nanoemulsions
Next Article in Special Issue
Current Application of Advancing Spectroscopy Techniques in Food Analysis: Data Handling with Chemometric Approaches
Previous Article in Journal
Uncovering the Biotechnological Importance of Geotrichum candidum
Previous Article in Special Issue
Phenotypic Analysis of Fourier-Transform Infrared Milk Spectra in Dairy Goats
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Heavy Metal Detection in Fritillaria thunbergii Using Laser-Induced Breakdown Spectroscopy Coupled with Variable Selection Algorithm and Chemometrics

1
College of Biosystems Engineering and Food Science, Zhejiang University, 866 Yuhangtang Road, Hangzhou 310058, China
2
Department of Agricultural and Bio-Resource Engineering, Abubakar Tafawa Balewa University, Bauchi PMB 0248, Nigeria
3
College of Mathematics and Computer Science, Zhejiang A&F University, Hangzhou 311300, China
4
Key Laboratory of Spectroscopy Sensing, Ministry of Agriculture and Rural Affairs, Hangzhou 310058, China
*
Author to whom correspondence should be addressed.
Foods 2023, 12(6), 1125; https://doi.org/10.3390/foods12061125
Submission received: 29 January 2023 / Revised: 3 March 2023 / Accepted: 5 March 2023 / Published: 7 March 2023

Abstract

:
Environmental and health risks associated with heavy metal pollution are serious. Human health can be adversely affected by the smallest amount of heavy metals. Modeling spectrum requires the careful selection of variables. Hence, simple variables that have a low level of interference and a high degree of precision are required for fast analysis and online detection. This study used laser-induced breakdown spectroscopy coupled with variable selection and chemometrics to simultaneously analyze heavy metals (Cd, Cu and Pb) in Fritillaria thunbergii. A total of three machine learning algorithms were utilized, including a gradient boosting machine (GBM), partial least squares regression (PLSR) and support vector regression (SVR). Three promising wavelength selection methods were evaluated for comparison, namely, a competitive adaptive reweighted sampling method (CARS), a random frog method (RF), and an uninformative variable elimination method (UVE). Compared to full wavelengths, the selected wavelengths produced excellent results. Overall, RC2, RV2, RP2, RSMEC, RSMEV and RSMEP for the selected variables are as follows: 0.9967, 0.8899, 0.9403, 1.9853 mg kg−1, 11.3934 mg kg−1, 8.5354 mg kg−1; 0.9933, 0.9316, 0.9665, 5.9332 mg kg−1, 18.3779 mg kg−1, 11.9356 mg kg−1; 0.9992, 0.9736, 0.9686, 1.6707 mg kg−1, 10.2323 mg kg−1, 10.1224 mg kg−1 were obtained for Cd Cu and Pb, respectively. Experimental results showed that all three methods could perform variable selection effectively, with GBM-UVE for Cd, SVR-RF for Pb, and GBM-CARS for Cu providing the best results. The results of the study suggest that LIBS coupled with wavelength selection can be used to detect heavy metals rapidly and accurately in Fritillaria by extracting only a few variables that contain useful information and eliminating non-informative variables.

Graphical Abstract

1. Introduction

The World Health Organization (WHO) reports that herbal medicines remain the primary treatment for a number of diseases in developing countries [1]. The use of nutraceuticals and medicinal products derived from medicinal herbs is increasing even in developed countries [1]. The consumption of healthy herbs is currently receiving considerable attention, and there is a focus on consuming products that are as natural as possible and have as little contamination as possible. The demand for nutraceuticals and herbal dietary supplements has increased significantly in recent years. Plants and other natural materials are used to make these products. There is, therefore, a high probability of heavy metal contamination [2]. The food chain may introduce heavy metals to humans through the accumulation of heavy metals in the environment [3,4,5,6]. There are numerous anthropogenic sources and activities that lead to heavy metal contamination of the environment, such as mining, traffic, agriculture, and industrial processes. However, heavy metals are naturally present in the Earth’s crust [6,7]. Risks associated with heavy metals can be represented in a variety of ways [8]. They can be found in the air we breathe, the food we eat, and the water we drink. They can also be found in soil and dust, which can be inhaled or ingested. As a result, they can enter the body and cause health problems.
In traditional Chinese medicine (TCM), heavy metals such as arsenic (As), mercury (Hg), lead (Pb), copper (Cu) and cadmium (Cd) are of particular concern [2]. Human health is at risk from these pollutants. The intake of excessive amounts of heavy metals is detrimental to the human body because it results in neurotoxicity, organic damage, and diseases of the skin and blood [9,10].
Fritillaria thunbergii Miq is a perennial herbaceous plant that is native to the provinces of Zhejiang, Jiangsu, and Anhui in China [11]. A dry bulb of Thunbergii fritillaria (Zhebeimu), a plant from the family of Thunbergii, is frequently used in Chinese medical clinical practice for treating coughs caused by wind and phlegm heats, as well as bronchitis, inflammation, hypertension, gastric ulcer, diarrhea, and bacterial infections [12]. Furthermore, Zhebeimu is extensively used to treat leukemia that is resistant to drugs [13]. Heavy metals have become one of the most serious safety concerns due to increasingly stringent TCM regulations [14,15,16]. Thus, accurate detection of heavy metal concentrations in TCM is crucial. Among the methods commonly used for heavy metal detection are atomic absorption spectroscopy (AAS), atomic fluorescence spectroscopy (AFS), X-ray fluorescence spectroscopy (XRFS), inductively coupled plasma optical emission spectroscopy (ICP-OES) [17,18], and electrochemical methods, especially stripping and cyclic voltammetry, which are commonly used methods for detection. Voltammetry is a sensitive electrochemical method that is widely used for heavy metal detection [19,20,21]. A traditional heavy metal detection method involves sampling, pretreatment, and laboratory chemical analysis; all of which are time-consuming, costly, and require extensive preparation [22,23].
Multi-elemental detection can be achieved through laser-induced breakdown spectroscopy (LIBS) [24]. LIBS uses pulsed laser ablation to create plasma on a sample and then detects and analyzes the emission light emanating from the plasma. With LIBS, a sample does not need to be prepared prior to analysis, thus allowing for rapid results [25], minimum requirements for small samples [26], and cost-effective instrumentation [27]. LIBS is widely used in a wide variety of industries [28], including mining [29], plastics [30], biomedicine [31,32], food [33], and the environment [34]. LIBS can provide information on the composition of samples within a short period of time, as well as the element content of samples. In comparison with other detection technologies, LIBS has many advantages, such as the requirement for fewer samples, the lack of complex pretreatment, the ability to measure multiple elements simultaneously, and the possibility of rapid implementation [35].
A number of studies have been conducted in recent years that focus on using LIBS to detect heavy metal pollution. Wang et al. [36] applied LIBS to detect cadmium content in rice. The LIBS analysis of rice stems demonstrated that it is an effective method for detecting cadmium. Su et al. [37] simultaneously and quantitatively analyzed the heavy metals in Sargassum fusiforme by using laser-induced breakdown spectroscopy. Liu et al. (2020) [38] analyzed cosmetics for trace lead and cadmium through laser-induced breakdown spectroscopy and ultrasound-assisted extraction. Rehan et al. [39] analyzed henna paste, and fresh leaves and soils were tested with LIBS to detect lead and nutrients. Wang et al. [40] conducted an analysis of lead and copper in Ligusticum wallichii using LIBS. Lead (Pb) levels in soil were quantitatively analyzed by Zhao et al. [41]. It was demonstrated that the dual-pulse laser-induced breakdown spectroscopy (DP-LIBS) was an efficient spectroscopic tool for improving the quantitative analysis of Pb heavy metal in soil. Zhu et al. [42] performed an analysis of the content of arsenic in traditional Chinese medicine using laser-induced breakdown spectroscopy (LIBS). Rehan et al. [43] used LIBS to assess the amount of toxic heavy metals (Pb, Cr, Ni) present in different brands of face foundation powder. Zhu et al. [44] detected lead in rhododendron leaves using laser-induced breakdown spectroscopy assisted by laser-induced fluorescence. Peng et al. [45] analyzed rice leaves using collinear DP-LIBS to determine their chromium content. Yang et al. [46] determined the content of lead and cadmium in rice using LIBS. These studies are indicative of the increasing consolidation of LIBS associated with chemometrics methods for the analysis of heavy metals.
However, despite the increase of fast and clean methods for TCM analysis [47], the application of gradient-boosting machine learning algorithms has rarely been used for heavy metal prediction in TCM. Thus, detecting heavy metals in different varieties of Fritillaria based on LIBS technology combined with a gradient boosting machine is unique and important for monitoring human exposure and establishing effective environmental control strategies.
Considering the potential risks associated with heavy metals in traditional Chinese medicine (TCM), the present work investigated the feasibility of LIBS combined with chemometrics in measuring cadmium (Cd), copper (Cu) and lead (Pb) simultaneously in twelve (12) different varieties of Fritillaria thunbergii. However, the specific objectives were as follows: (a) to test the feasibility of the gradient boosting machine (GBM) as a method of measuring heavy metals in different varieties of Fritillaria using both full and extracted variables, (b) to verify the effectiveness of using three variable selection methods, namely, competitive adaptive reweighed sampling (CARS), random frog (RF) and uninformative variable elimination (UVE) by comparing model (GBM, SVR and PLSR) performances, and (c) to establish a quantitative analysis model for heavy metals based on full and extracted variables.

2. Materials and Methods

2.1. Sample Collection and Preparation

A total of twelve (12) different varieties of Fritillaria thunbergia were used in the experiment. The varieties were provided by the Faculty of Biosystems Engineering and Food Science (Zhejiang University, Hangzhou, China). Different copper (Cu), cadmium (Cd) and lead (Pb) samples were prepared in the laboratory using Cu(NO3)2, Cd(NO3)2. 4H2O and Pb (NO3)2 [40,44,46,48], respectively. Fritillaria varieties were randomly divided into eight groups to obtain samples with varying levels of Cu, Cd, and Pb. The first group was designated as a control group, whereas the remaining seven groups were designated as treatment groups. In order to accurately quantify Cu, Cd, and Pb concentrations in Fritillaria samples, the solution was artificially contaminated for 48 h at 4 °C and rinsed three times with super-pure water. This was done to simulate the effects of heavy metal pollution in the environment and to determine the best course of action for mitigating the environmental damage caused by these pollutants. Following drying at room temperature, all samples were milled at high speeds using a high-speed pulverizer (FW100, TAISITE, Tianjin, China). In order to produce pellets from these ground samples, they were compressed using a tablet compressor (FY-24, SCJS, Tianjin, China), of 1.5 cm in diameter, at a pressure of 25 kN for the duration of one minute. In total, 288 pellets were prepared.

2.2. Experimental Setup

A LIBS self-assembled schematic diagram is available in [49]; Figure 1 illustrates the methods used in this experiment. Laser pulses were generated at 532 nm with a maximum energy of 200 mJ and a pulse width of 8 ns using a Q-switched Nd: YAG pulse laser (Vlite 200, Beamtech, Beijing, China). A planoconvex lens (f = 100 mm) finally focused the laser onto the sample surface after passing through the optical system. The laser ablation generated plasma which emitted electromagnetic waves that diffused outward. In order to measure the waves, a light collector was used, and the waves were captured by a spectrometer (SR-500i-A-R, Andor Technology, Belfast, UK) along with an intensified charge-coupled device (ICCD) camera (DH334T-18F-03, Andor Technology, Belfast, UK). The laser Q-switch and ICCD camera were delayed using a delay generator (DG645, Stanford Research Systems, Sunnyvale, CA, USA). Several parameters were optimized before the experiment, including the laser energy of 60 mJ, the delay time of 1.5 µs, and the gate width of 10µs. The Fritillaria pellets were automatically placed, and the laser ablation path was controlled using a 4 × 4 array of craters designed using automatic x-y-z translation. Laser pulses accumulated five times faster in each crater. For each sample, an average of 80 spectra (4 × 4 × 5) were taken in order to reduce fluctuation between the laser points. Approximately one minute was required to collect LIBS information for one sample.

2.3. Determination of Heavy-Metals Reference Value

Inductively coupled plasma mass spectrometry (ICP-MS) was used to determine the contents of three (3) heavy metals (Cu, Cd, and Pb) in Fritillaria thunbergii. Table 1 lists the statistics of heavy metal contents in different samples. A number of procedures are involved in the ICP-MS analysis of the samples, including digestion of the sample, filtration and purification of the digestion solution, and detection of the digestion solution using ICP-MS (ELAN DRC-e, Perkin Elmer, USA). The pH value as well as all steps in the ICP-MS analysis were carried out by experimental technicians at the Zhejiang College of Biosystems Engineering and Food Science, Zhejiang University. This was a preliminary attempt to detect multiple heavy metals (Cd, Cu and Pb) in Fritillaria thubergii using the LIBS technique. The pH value was determined by measuring the electrical conductivity of the sample solution with a pH meter. The ICP-MS analysis was performed using a quadrupole-based ICP-MS system which allowed the technicians to accurately measure the concentrations of elements in the sample solution. Similar steps are described by Su et al. [37]. Unscrambler X, version 10.1 (CAMO Software AS, Oslo, Norway, 2011) was used for the descriptive statistics (file imported in MATLAB format).

2.4. Data Analysis

2.4.1. Spectral Modeling

The gradient boosting technique (GB) is a machine-learning method used in regression and classification problems. As a result, a prediction model is produced as an ensemble of weak prediction models. Every step evaluates the model values at each training data point, using the residuals of previous steps to minimize the loss function [50]. A GBM utilizes the best practices to avoid overfitting the classification machine. A subsample of the training data is randomly selected (without replacement) from the full dataset for each iteration in order to fit the base learner for that iteration. Figure 2 illustrates the main processes involved in gradient boosting through a simplified flow chart.
In comparison to other machine learning methods, gradient boosting machines (GBMs) have several advantages. GBM, in addition to its complex classification capability, allows soft classification, which entails calculating the probability of each sample being a member of each class, rather than labeling every sample as part of a single category (hard classification). It facilitates the assessment of the reliability of the statistical model (the potential for overfitting) and the study of the chemical-physical properties of the model, thus fostering the development of further qualitative and quantitative research studies.
In addition, it provides a natural measure of how significant each spectral feature is for classifying data, something that is usually lost in the black-box nature of many machine learning algorithms, which include artificial neural networks. The relative importance of variables is computed in GBM. A refinement of the split criterion is computed at each split in each tree (MSE for regression). Averaging the improvements made by each variable over every tree that uses that variable follows. In the split criterion, the variables with the greatest average decrease are listed as the most significant [50].When it comes to GBM modeling, there are a variety of tuning parameters available. The following variables were used in this study: (boosting_type = ‘gbdt’, num_leaves = 31, max_depth = −1, learning_rate = 0.1, n_estimators = 100). In order to forecast discrete values, SVR is an approach that uses supervised learning. By comparison, it aims to determine the hyperplane with the most points, or the line of best fit [51,52]. PLS uses projections to build linear regression models using variables and observables [53,54,55]. A significant amount of collinearity can be analyzed by this algorithm, unlike previous algorithms.

2.4.2. Variable Selection Methods

Generally, LIBS data display high covariance due to the capability of the technique to measure multiple emission lines associated with the same element or species. Moreover, each peak is a result of a combination of many factors. It is therefore possible to reduce the covariance and complexity of the model by selecting variables. Selection of more explanatory variables improves the understanding of the multivariate system [56]. The removal of noisy areas, such as the extremes of each spectrometer and variables without analytical information, is also likely to increase the explained variance and enhance the accuracy of the model [56]. To select the optimal variable, variables should be selected and eliminated. In this study, three methods of variable selection were used to simplify the calculation process and improve model performance (CARS, RF, and UVE).
In CARS, wavelengths are selected using a survival-of-the-fittest principle [57]. Firstly, wavelengths with small regression coefficients are removed by using an exponentially decreasing function (EDF). An EDF equation is then used to calculate the wavelength ratio. The following steps are involved in each sampling run: (a) the Monte Carlo (MC) principle is used to sample models; (b) EDF-based wavelength selection is performed; (c) adaptive reweighted sampling is used for competitive wavelength selection and (d) cross-validation evaluation of the subset is performed [57]. A subset of wavelengths with the lowest root mean squared error of cross-validation (RMSECV) is retained as the effective wavelength, and wavelengths with little or no effective information are eliminated [58].
The RF method is iterative in nature. There are primarily three steps in the random frog algorithm: (1) A random selection of features is used to create an initialized feature subset. (2) Iteration is performed until a candidate feature subset is selected. This is accepted with a certain probability, then replaced, and this step is looped until the desired number of iterations is achieved. (3) As a measure of feature importance, the selection probability of each feature is calculated [59]. The total number of feature subsets can be determined after multiple iterations. Its selection probability can be calculated as
P j = N j N ,   j = 1 ,   2 ,   n
Nj denotes the frequency of the jth feature, j = 1, 2, …, n, selected from the features. For each feature, Nj is the feature subsets after iteration, and N is the feature frequency.
Accordingly, the more optimal a feature, the more likely it is to be selected for inclusion in these subsets of features. As a result, features can be selected according to feature importance. The UVE method uses regression coefficients from a PLS model to select variables. This method is useful in eliminating non-informative variables, and the remaining variables can be used to analyze and classify chemicals [60,61].

2.4.3. Model Evaluation and Calculation

In order to assess the performance of the calibration model, several evaluation indices were calculated. The point-to-point fluctuations in the spectra were reduced by area normalization [45]. The data were all normalized prior to modeling. The accuracy of the model was assessed by determining the root mean square error of calibration (RMSEC), the root mean square error of validation (RMSEV), the root mean square error of prediction (RMSEP) and the coefficient of calibration (RC2), coefficient of validation (RV2) and coefficient of prediction (RP2) based on the predicted results. In summary, a good calibration model should have a small RMSEC, and RMSEP, as well as large values for (RC2) and (RP2) [37,48,62,63]. The calculations were performed using Python with Scikit-Learn and the figures were generated using Origin 2022.

3. Results and Discussion

3.1. Spectra Analysis

In LIBS, atoms and ions are expelled from a generated plasma as a result of their excited state [64]. Figure 3 illustrates the normalized spectra of 12 Fritillaria varieties. In accordance with the National Institute of Standards and Technology, USA, Electronic Database, characteristic lines for Cd, Cu, and Pb were identified. In this study, however, the purpose is to quantitatively analyze the content of three heavy metals (Cd, Cu and Pb), and because the LIBS spectra of the different varieties have similar curves, it is difficult to observe the LIBS spectra simultaneously and quantitatively to analyze the heavy metal content of Fritillaria simultaneously and quantitatively. Hence, the LIBS data must be further analyzed using chemometric methods.

3.2. Heavy Metals Prediction Using Full and Selected Variables

We performed LIBS analyses in accordance with best practices [65]. The full spectrum was divided into three parts: calibration (60%), validation (28%), and prediction (12%). A summary of the GBM results for the full wavelength and the selected wavelength is presented in Table 1 and Table 2. The peaks in full LIBS data usually exhibit high covariance due to the capability of measuring multiple emission lines of the same element or species and the fact that a number of x variables are required to generate a single peak. A useful range of wavelengths should be selected from the entire range of wavelengths in order to simplify the LIBS calibration models and improve prediction accuracy. As a result, three methods of variable selection (CARS, RF, and UVE) were used to select the informative variables from the whole wavelength range to simplify the calculation process and to improve the performance of the model.

3.2.1. Cd Content Prediction Using Full and Selected Variables

All three variable selection methods (CARS, RF and UVE) were computed with SVR, PLSR and GBM to predict Cd (Table 2), and the best result was achieved with UVE-GBM. The calibration RC2 achieved by the model was 0.9967, the RMSEC was 1.9853 mg kg−1, the validation set RV2 achieved 0.8899, the RMSEV was 11.3934 mg kg−1, the prediction RP2 achieved 0.9403 and the RMSEP was 8.5344 mg kg−1. Figure 4b displays the scatter plots of the reference value and the prediction value for the Cd using the UVE-GBM.
RF-SVR and RF-PLSR also produced comparatively good results. For the RF-SVR calibration set, RC2 = 0.9999, RMSEC = 0.1000 mg kg−1, the validation set RV2 achieved 0.9287, the RMSEV was 9.1671 mg kg−1 and the prediction set was RP2 = 0.9322, RMSEP = 9.0933 mg kg−1. For the RF-PLSR calibration set RC2 = 0.9825, RMSEC = 4.6313 mg kg−1, the validation set RV2 achieved 0.9313, the RMSEV was 9.0014 mg kg−1 and the prediction set was RP2 = 0.9182, RMSEP = 9.9903 mg kg−1. On the other hand, CARS-SVR, UVE-SVR and CARS-GBM, RF-GBM produced lower prediction outcomes (Table 3). Meanwhile, for the full spectra, GBM was the best (Table 2); for the calibration set, RC2 was 0.9984, RMSEC = 1.3693 mg kg−1, the validation set RV2 achieved 0.8924, the RMSEV was 11.2625 mgkg−1, and the prediction set was RP2 = 0.9139, RMSEP = 10.2494 mg kg−1, as shown in Figure 4a,b below. GBM combined with UVE offers superior predictions than full wavelength GBM alone. GBM-UVE performs better with Cd because it is able to handle the large number of variables and select the most important ones.

3.2.2. Cu Content Prediction Using Full and Selected Variables

For the prediction of Cu content, GBM combined with CARS achieved the best outcome with the following calibration set: RC2 = 0.9933, RMSEC = 5.9332 mg kg−1, the validation set RV2 achieved 0.9316, the RMSEV was 18.3779 mg kg−1, and the prediction set was RP2 = 0.9665, RMSEP = 11.9356 mg kg−1. This was followed by RF-SVR: for the calibration set RC2 = 0.9860, RMSEC = 8.5699 mg kg−1, the validation set RV2 achieved 0.9640, the RMSEV was 13.3319 mg kg−1, and the prediction set was RP2 = 0.9658, RMSEP = 12.0658 mg kg−1; UVE-SVR: for the calibration set RC2 = 0.9715, RMSEC = 12.2395 mg kg−1, the validation set RV2 achieved 0.9412, the RMSEV was 17.0324 mg kg−1, and the prediction set was RP2 = 0.9648, RMSEP = 12.2429 mg kg−1 and UVE-GBM: for the calibration set (RC2 = 0.9964, RMSEC = 4.3168 mg kg−1), the validation set RV2 achieved 0.9304, the RMSEV was 18.5337 mg kg−1, and the prediction set was RP2 = 0.9648, RMSEP = 12.2323 mg kg−1. Figure 5 shows a scatter plot of the reference value and prediction value for Cu content using the CARS-GBM model. Although for the full spectra, GBM was the best (Table 2), for the calibration set RC2 = 0.9979, RMSEC = 3.2885 mg kg−1, the validation set RV2 achieved 0.9308, the RMSEV was 18.4815 mg kg−1, and the prediction set was RP2 = 0.9596, RMSEP = 13.1156 mg kg−1. It can be seen in Figure 5a,b that the prediction result from GBM combined with CARS is better than the result from the full wavelength. GBM-CARS is better for Cu because it is able to capture interactions between variables.

3.2.3. Pb Content Prediction Using Full and Selected Variables

As shown in Table 3, RF combined with SVR produced the best result for the variable selection methods for Pb content prediction with the following calibration set: RC2 = 0.9992, RMSEC = 1.6707 mg kg−1; the validation set RV2 achieved 0.9736, the RMSEV was 10.2323 mg kg−1, and the prediction set was RP2 = 0.9686, RMSEP = 10.1224 mg kg−1. This was followed by UVE-GBM: calibration set RC2 = 0.9992, RMSEC = 1.7248 mg kg−1, the validation set RV2 achieved 0.9562, the RMSEV was 13.1967 mg kg−1, and the prediction set was RP2 = 0.9609, RMSEP = 11.3113 mg kg−1 and CARS-PLSR: calibration set RC2 = 0.9683, RMSEC = 11.2039 mg kg−1, the validation set RV2 achieved 0.9718, the RMSEV was 10.5798 mgkg−1, and the prediction set was RP2 = 0.9599, RMSEP = 11.4563 mg kg−1. Although for the full spectra, GBM was also the best (Table 2), for the calibration set RC2 = 0.9998, RMSEC = 0.8488 mg kg−1, the validation set RV2 achieved 0.9673, the RMSEV was 11.3933 mg kg−1, and the prediction set was RP2 = 0.9635, RMSEP = 10.9220 mg kg−1. It can be noticed that the prediction result for SVR combined with RF, not GBM, as it is with Cd and Cu, as shown in the Figure 6a,b, it is slightly better than that for the full wavelength. SVR-RF is better for Pb because it is able to capture non-linear relationships between variables.

4. Discussion

Among the perennial herbaceous plants found mainly in the Zhejiang, Jiangsu, and Anhui provinces of China, Fritillaria thunbergii Miq. [11] was selected as the research object for collecting LIBS data and assessing its heavy metal contents (Cd, Cu, and Pb). The presence of heavy metals in excess of the standard poses a significant health risk. Due to increasing regulation of (TCM), heavy metals have become a priority pollutant in TCM and a serious safety concern [14,15,16]. Thunbergii fritillaria bulbs (Zhebeimu) are most commonly used in Chinese medical clinical practice to treat coughs caused by wind-heat and phlegm-heat, bronchitis, inflammation, hypertension, gastric ulcers, diarrhea, and bacterial infections [12].
In addition, Zhebeimu is widely used in the treatment of leukemia that is resistant to drugs [13]. Despite the development of fast and clean methods for TCM analysis [47]. Using LIBS to detect heavy metals in Fritillaria Thunbergii has still not been reported. Therefore, the ability to detect heavy metals in different varieties of Fritillaria using LIBS technology is a unique and vital process for establishing effective environmental control strategies and monitoring human exposure to heavy metals. In light of the above considerations and the advantages of the LIBS technique, in this study, a quantitative and simultaneous analysis of the contents of three heavy metals in Fritillaria thunbergii was performed. LIBS data, however, typically show a high degree of covariance due to the ability of LIBS to measure multiple emission lines from the same species or element; besides, several variables are responsible for creating each peak. Therefore, selecting variables may result in a reduction in the covariance and complexity of the model. It is also helpful to select more explanatory variables because this allows a better understanding of the multivariate system in terms of its chemical characteristics [56]. As well as improving the model fit and increasing the explained variance, removing noisy areas, such as extreme regions of the spectrometer and x variables with no analytical information, can also help to remove noisy regions [56].
In multivariate analysis, the matrix effect and fluctuations in the LIBS spectrum can be taken into account, in addition to the fluctuation in the LIBS spectrum from shot to shot. Recent years have seen an extensive use of chemometric methods, such as partial least squares (PLSR) and support vector machines (SVM) in the analysis of LIBS spectra for multivariable analysis [66,67]. An analysis of three variable selection methods (CARS, RF, and UVE) was conducted in order to simplify the calculation process and improve model performance. PLSR, SVM, and GBM models were run on the full spectra and spectral variables from CARS, RF, and UVE, respectively. A comparison of the performance of PLSR, SVM and GBM models was made by determining the root mean square error of calibration and prediction (RMSEC and RMSEP), as well as the greatest correlation coefficient square (R2) between calibration and prediction sets (RC2 and RP2). It has been found that the combination of SVR and CARS yielded the smallest number of variables selected.
The combination of GBM with UVE yielded the greatest number of variables for Cd and Cu, and for Pb, GBM with RF yielded a limited number of variables (Table 3). The results of the different feature selection methods differ when compared when calibration, validation, and prediction sets are analyzed. Different feature selection methods were best for the different heavy metals. For Cd: GBM combined with UVE obtained the best performance; for Cu: GBM combined with CARS obtained the best performance, whereas for Pb: SVR combined with RF obtained the best performance. Several studies have demonstrated that there is no single technique for selecting features that is universally optimal [68], and multiple subsets of features are usually equally effective in predicting the data [69,70,71]. There is no doubt that full spectra contain essential information for elemental analysis but inevitably contain irrelevant information and noise, which weaken the model’s capability [63,72,73]. The data analysis step in LIBS, as in other fields of spectroscopic analysis, is heavily constrained both by the high-dimensional input spaces and their inherent sparsity. In addition to reducing the measurement and storage requirements for LIBS data, properly selecting spectral features can facilitate the visualization and understanding of the data and enable more timely and cost-effective classification methods to be developed. The selection of feature variables needs to be applied to reduce computational complexity. As described above, it is worth noting that the number of variables selected by different feature selection methods varies widely (Table 3). For Cd, CARS showed the least variables (37), followed by RF (151) and UVE (192) and for Cu, CARS also showed the least variables (66), followed by RF (120), UVE (231) and finally for Pb, CARS showed the least variables (17), followed by UVE (93), RF (124), as can be seen from Figure 7. All methods select the informative region around 200–1000 nm, which is consistent with the Cd, Cu and Pb results in previous literature [1,37,39,40,48,63,74,75], indicating that these variable intervals are the informative variables.
The different feature selection algorithms indicate some essential variables, so the variables in those intervals must be supported to improve the model’s prediction ability, indicating those regions are informative variable intervals. Table 3 summarizes the results of CARS, RF, and UVE variable selection methods. Based on Figure 5, all selection methods perform significantly better in the test set than when compared with the full spectrum. The spectral matrix was then transformed by reducing the spectral information of the LIBS measurements to the most relevant variables that contained the most relevant spectral information of the respective heavy metals in Fritillaria. Multivariate analysis could address laser-to-sample interaction, experimental parameter variance, and matrices, among other factors [76].

5. Conclusions

The study indicated that when using the wavelength selection method, only a limited number of useful variables were extracted, and non-informative variables were eliminated. Therefore, the study explored effective variables by variable selection methods. The experimental results showed that all three methods applied could accomplish variable selection effectively, among which GBM-UVE for Cd, SVR-RF for Pb, and GBM-CARS for Cu produced the best results. Results of some recent heavy metal detection using LIBS are presented in Table 4 for comparison with the current work. Table 4 provides a comprehensive comparison of the results from the current work with those obtained from recent LIBS studies. It allows for a direct comparison of the accuracy of the detection methods. This study demonstrated the potential of LIBS coupled with variable selection and chemometrics as a tool for the rapid detection of heavy metals in varieties of Fritillaria thunbergii. It is essential to select variables that are associated with spectral information, so that subsequent modeling analysis can be based on more concise and effective spectral data. Additionally, the wavelengths selected could provide a theoretical foundation for the development of new instruments.

Author Contributions

Conceptualization, F.L. and M.H.K.; data curation, M.H.K., X.L. and R.C.; formal analysis, M.H.K. and R.C.; funding acquisition, F.L.; investigation, M.H.K.; methodology, M.H.K.; project administration, F.L.; resources, F.L. and W.K.; software, M.H.K. and M.L.G.; validation, M.H.K. and F.L.; visualization, M.H.K. and F.L.; writing—original draft, M.H.K.; writing—review and editing, F.L., M.H.K., M.L.G., R.C., X.L. and W.K. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by National Natural Science Foundation of China (61975174) and Science and Technology Department of Zhejiang Province (2021C02023).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors wish to acknowledge the funding agency stated above.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jiang, Y.; Kang, J.; Wang, Y.; Chen, Y.; Li, R. Rapid and sensitive analysis of trace leads in medicinal herbs using laser-lnduced breakdown spectroscopy-laser-induced fluorescence (libs-lif). Appl. Spectrosc. 2019, 73, 1284–1291. [Google Scholar] [PubMed]
  2. Gyamfi, E.T. Metals and metalloids in traditional medicines (ayurvedic medicines, nutraceuticals and traditional chinese medicines). Environ. Sci. Pollut. Res. Int. 2019, 26, 15767–15778. [Google Scholar] [CrossRef] [PubMed]
  3. Dudka, S.; Miller, W.P. Accumulation of potentially toxic elements in plants and their transfer to human food chain. J. Environ. Sci. Health Part B 1999, 34, 681–708. [Google Scholar] [CrossRef] [PubMed]
  4. Dube, A.; Zbytniewski, R.; Kowalkowski, T.; Cukrowska, E.; Buszewski, B. Adsorption and migration of heavy metals in soil. Pol. J. Environ. Stud. 2001, 10, 1–10. [Google Scholar]
  5. Ercal, N.; Gurer-Orhan, H.; Aykin-Burns, N. Toxic metals and oxidative stress part i: Mechanisms involved in metal-induced oxidative damage. Curr. Top Med. Chem. 2001, 1, 529–539. [Google Scholar] [CrossRef]
  6. Duruibe, J.; Ogwuegbu, M.O.C.; Egwurugwu, J. Heavy metal pollution and human biotoxic effects. Int. J. Phys. Sci. 2007, 2, 112–118. [Google Scholar]
  7. Wei, B.; Yang, L. A review of heavy metal contaminations in urban soils, urban road dusts and agricultural soils from china. Microchem. J. 2010, 94, 99–107. [Google Scholar] [CrossRef]
  8. Orecchio, S.; Amorello, D.; Barreca, S.; Valenti, A. Wood pellets for home heating can be considered environmentally friendly fuels? polycyclic aromatic hydrocarbons (pahs) in their ashes. Microchem. J. 2016, 124, 267–271. [Google Scholar] [CrossRef]
  9. Tang, Y.; Wang, X.; Jia, J. Mercury poisoning presenting as sporadic creutzfeldt-jakob disease: A case report. Ann. Intern. Med. 2015, 162, 462–463. [Google Scholar] [CrossRef]
  10. Chan, N.C.N.; Chan, K.P. Coarse basophilic stippling in lead poisoning. Blood 2017, 129, 3270. [Google Scholar] [CrossRef] [Green Version]
  11. Hao, D.-C.; Gu, X.-J.; Xiao, P.-G.; Peng, Y. Phytochemical and biological research of fritillaria medicine resources. Chin. J. Nat. Med. 2013, 11, 330–344. [Google Scholar] [CrossRef]
  12. Hill, C.B.; Roessner, U. Metabolic Profiling of Plants by GC–MS. In The Handbook of Plant Metabolomics; John Wiley & Sons: Hoboken, NJ, USA, 2013. [Google Scholar]
  13. Ju, W. Clinical observation on effects of adjuvant chemotherapy with fritillary thunbergii bulb complex granules on symptoms of refractory acute leukemia patients. Beijing J. Tradit. Chin. Med. 2013, 32, 833–835. [Google Scholar]
  14. Tian, M.; Fang, L.; Yan, X.; Xiao, W.; Row, K.H. Determination of heavy metal ions and organic pollutants in water samples using ionic liquids and ionic liquid-modified sorbents. J. Anal. Methods Chem. 2019, 2019, 1948965. [Google Scholar] [CrossRef]
  15. Yin, H.; Truskewycz, A.; Cole, I. Quantum dot (qd)-based probes for multiplexed determination of heavy metal ions. Microchim. Acta 2020, 187, 336. [Google Scholar] [CrossRef]
  16. Fu, F.; Wang, Q. Removal of heavy metal ions from wastewaters: A review. J. Environ. Manag. 2011, 92, 407–418. [Google Scholar] [CrossRef]
  17. Fortes, F.J.; Moros, J.; Lucena, P.; Cabalín, L.M.; Laserna, J.J. Laser-induced breakdown spectroscopy. Anal. Chem. 2013, 85, 640–669. [Google Scholar] [CrossRef]
  18. Kim, G.; Kwak, J.; Kim, K.-R.; Lee, H.; Kim, K.-W.; Yang, H.; Park, K. Rapid detection of soils contaminated with heavy metals and oils by laser induced breakdown spectroscopy (libs). J. Hazard. Mater. 2013, 263, 754–760. [Google Scholar] [CrossRef]
  19. Sawan, S.; Maalouf, R.; Errachid, A.; Jaffrezic-Renault, N. Metal and metal oxide nanoparticles in the voltammetric detection of heavy metals: A review. TrAC Trends Anal. Chem. 2020, 131, 116014. [Google Scholar] [CrossRef]
  20. Malik, L.A.; Bashir, A.; Qureashi, A.; Pandith, A.H. Detection and removal of heavy metal ions: A review. Environ. Chem. Lett. 2019, 17, 1495–1521. [Google Scholar] [CrossRef]
  21. Amorello, D.; Barreca, S.; Gulli, E.; Orecchio, S. Platinum and rhodium in wine samples by using voltammetric techniques. Microchem. J. 2017, 130, 229–235. [Google Scholar] [CrossRef]
  22. Rai, N.K.; Rai, A.K. Libs—An efficient approach for the determination of cr in industrial wastewater. J. Hazard. Mater. 2008, 150, 835–838. [Google Scholar] [CrossRef] [PubMed]
  23. Gondal, M.; Seddigi, Z.; Nasr, M.; Gondal, B. Spectroscopic detection of health hazardous contaminants in lipstick using laser induced breakdown spectroscopy. J. Hazard. Mater. 2010, 175, 726–732. [Google Scholar] [CrossRef]
  24. Hahn, D.W.; Omenetto, N. Laser-induced breakdown spectroscopy (libs), part i: Review of basic diagnostics and plasma–particle interactions: Still-challenging issues within the analytical plasma community. Appl. Spectrosc. 2010, 64, 335A–366A. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Ma, S.; Tang, Y.; Ma, Y.; Chu, Y.; Chen, F.; Hu, Z.; Zhu, Z.; Guo, L.; Zeng, X.; Lu, Y. Determination of trace heavy metal elements in aqueous solution using surface-enhanced laser-induced breakdown spectroscopy. Opt. Express 2019, 27, 15091–15099. [Google Scholar] [CrossRef] [PubMed]
  26. Ma, S.; Tang, Y.; Ma, Y.; Dong, D.; Guo, L.; Zhu, H.; Liu, J.; Lu, Y. The pH effect on the detection of heavy metals in wastewater by laser-induced breakdown spectroscopy coupled with a phase transformation method. J. Anal. At. Spectrom. 2020, 35, 198–203. [Google Scholar] [CrossRef]
  27. Hahn, D.W.; Omenetto, N. Laser-induced breakdown spectroscopy (libs), part ii: Review of instrumental and methodological approaches to material analysis and applications to different fields. Appl. Spectrosc. 2012, 66, 347–419. [Google Scholar] [CrossRef]
  28. Noll, R.; Fricke-Begemann, C.; Connemann, S.; Meinhardt, C.; Sturm, V. Libs analyses for industrial applications—An overview of developments from 2014 to 2018. J. Anal. At. Spectrom. 2018, 33, 945–956. [Google Scholar] [CrossRef] [Green Version]
  29. Zou, X.H.; Guo, L.B.; Shen, M.; Li, X.Y.; Hao, Z.Q.; Zeng, Q.D.; Lu, Y.F.; Wang, Z.M.; Zeng, X.Y. Accuracy improvement of quantitative analysis in laser-induced breakdown spectroscopy using modified wavelet transform. Opt. Express 2014, 22, 10233–10238. [Google Scholar] [CrossRef]
  30. Gondal, M.A.; Siddiqui, M.N. Identification of different kinds of plastics using laser-induced breakdown spectroscopy for waste management. J. Environ. Sci. Health Part A 2007, 42, 1989–1997. [Google Scholar] [CrossRef]
  31. Unnikrishnan, V.K.; Choudhari, K.S.; Kulkarni, S.D.; Nayak, R.; Kartha, V.B.; Santhosh, C.; Suri, B.M. Biomedical and environmental applicationsof laser-induced breakdown spectroscopy. Pramana 2014, 82, 397–401. [Google Scholar] [CrossRef]
  32. Ahmed, I.; Yang, J.; Law, A.W.L.; Manno, F.A.M.; Ahmed, R.; Zhang, Y.; Lau, C. Rapid and in situ optical detection of trace lithium in tissues. Biomed. Opt. Express 2018, 9, 4459–4471. [Google Scholar] [CrossRef]
  33. Tang, S.S.; Hao, Z.Q.; Ma, Y.Y.; Ma, S.X.; Guo, L.B.; Guo, Y.M.; Zeng, X.Y.; Chu, Y.W.; Lu, Y.F. Accuracy and stability improvement for meat species identification using multiplicative scatter correction and laser-induced breakdown spectroscopy. Opt. Express 2018, 26, 10119. [Google Scholar]
  34. Pandhija, S.; Rai, N.K.; Rai, A.K.; Thakur, S.N. Contaminant concentration in environmental samples using libs and cf-libs. Appl. Phys. B 2010, 98, 231–241. [Google Scholar] [CrossRef]
  35. Gondal, M.A.; Hussain, T.; Yamani, Z.H.; Baig, M.A. On-line monitoring of remediation process of chromium polluted soil using libs. J. Hazard. Mater. 2009, 163, 1265–1271. [Google Scholar] [CrossRef]
  36. Wang, W.; Kong, W.; Shen, T.; Man, Z.; Zhu, W.; He, Y.; Liu, F.; Liu, Y. Application of laser-induced breakdown spectroscopy in detection of cadmium content in rice stems. Front. Plant Sci. 2020, 11, 599616. [Google Scholar] [CrossRef]
  37. Su, L.; Shi, W.; Chen, X.; Meng, L.; Yuan, L.; Chen, X.; Huang, G. Simultaneously and quantitatively analyze the heavy metals in Sargassum fusiforme by laser-induced breakdown spectroscopy. Food Chem. 2021, 338, 127797. [Google Scholar] [CrossRef]
  38. Liu, Y.; Chu, Y.; Hu, Z.; Zhang, S.; Ma, S.; Khan, M.S.; Chen, F.; Zhang, D.; Guo, L.; Lau, C. High-sensitivity determination of trace lead and cadmium in cosmetics using laser-induced breakdown spectroscopy with ultrasound-assisted extraction. Microchem. J. 2020, 158, 105322. [Google Scholar] [CrossRef]
  39. Rehan, I.; Gondal, M.A.; Rehan, K.; Sultana, S.; Dastageer, M.A.; Al-Adel, F.F. Libs for the detection of lead in ready to use henna paste and nutrients in fresh henna leaves and cultivated soils. Talanta 2019, 199, 203–211. [Google Scholar] [CrossRef]
  40. Wang, J.; Xue, S.; Zheng, P.; Chen, Y.; Peng, R. Determination of lead and copper in ligusticum wallichii by laser-induced breakdown spectroscopy. Anal. Lett. 2017, 50, 2000–2011. [Google Scholar] [CrossRef]
  41. Zhao, S.; Song, C.; Gao, X.; Lin, J. Quantitative analysis of pb in soil by femtosecond-nanosecond double-pulse laser-induced breakdown spectroscopy. Results Phys. 2019, 15, 102736. [Google Scholar] [CrossRef]
  42. Zhu, C.; Lv, J.; Liu, K.; Chen, J.; Liu, K.; Li, G.; Lu, B.; Li, X. Rapid determination of arsenic in traditional chinese medicine by laser-induced breakdown spectroscopy (libs). Anal. Lett. 2022, 55, 2531–2541. [Google Scholar] [CrossRef]
  43. Rehan, I.; Gondal, M.A.; Rehan, K.; Sultana, S. Spectral diagnosis of health hazardous toxins in face foundation powders using laser induced breakdown spectroscopy and inductively coupled plasma-optical emission spectroscopy (icp-oes). Talanta 2020, 217, 121007. [Google Scholar] [CrossRef] [PubMed]
  44. Zhu, C.; Tang, Z.; Li, Q.; Zhou, R.; Lv, J.; Zhang, W.; Zhan, K.; Li, X.; Zeng, X. Lead of detection in rhododendron leaves using laser-induced breakdown spectroscopy assisted by laser-induced fluorescence. Sci. Total Environ. 2020, 738, 139402. [Google Scholar] [CrossRef] [PubMed]
  45. Peng, J.; He, Y.; Jiang, J.; Zhao, Z.; Zhou, F.; Liu, F. High-accuracy and fast determination of chromium content in rice leaves based on collinear dual-pulse laser-induced breakdown spectroscopy and chemometric methods. Food Chem. 2019, 295, 327–333. [Google Scholar] [CrossRef] [PubMed]
  46. Yang, P.; Zhou, R.; Zhang, W.; Yi, R.; Tang, S.; Guo, L.; Hao, Z.; Li, X.; Lu, Y.; Zeng, X. High-sensitivity determination of cadmium and lead in rice using laser-induced breakdown spectroscopy. Food Chem. 2019, 272, 323–328. [Google Scholar] [CrossRef] [PubMed]
  47. Li, K.; Yang, H.; Yuan, X.; Zhang, M. Recent developments of heavy metals detection in traditional chinese medicine by atomic spectrometry. Microchem. J. 2021, 160, 105726. [Google Scholar] [CrossRef]
  48. Liu, F.; Ye, L.; Peng, J.; Song, K.; Shen, T.; Zhang, C.; He, Y. Fast detection of copper content in rice by laser-induced breakdown spectroscopy with uni-and multivariate analysis. Sensors 2018, 18, 705. [Google Scholar] [CrossRef] [Green Version]
  49. Kabir, M.H.; Guindo, M.L.; Chen, R.; Sanaeifar, A.; Liu, F. Application of laser-induced breakdown spectroscopy and chemometrics for the quality evaluation of foods with medicinal properties: A review. Foods 2022, 11, 2051. [Google Scholar] [CrossRef]
  50. Huffman, C.; Sobral, H.; Terán-Hinojosa, E. Laser-induced breakdown spectroscopy spectral feature selection to enhance classification capabilities: A t-test filter approach. Spectrochim. Acta Part B At. Spectrosc. 2019, 162, 105721. [Google Scholar] [CrossRef]
  51. Gholizadeh, A.; Žižala, D.; Saberioon, M.; Borůvka, L. Soil organic carbon and texture retrieving and mapping using proximal, airborne and Sentinel-2 spectral imaging. Remote Sens. Environ. 2018, 218, 89–103. [Google Scholar] [CrossRef]
  52. Guindo, M.L.; Kabir, M.H.; Chen, R.; Liu, F. Particle swarm optimization and multiple stacked generalizations to detect nitrogen and organic-matter in organic-fertilizer using vis-nir. Sensors 2021, 21, 4882. [Google Scholar] [CrossRef]
  53. Douglas, R.K.; Nawar, S.; Cipullo, S.; Alamar, M.C.; Coulon, F.; Mouazen, A.M. Evaluation of vis-nir reflectance spectroscopy sensitivity to weathering for enhanced assessment of oil contaminated soils. Sci. Total Environ. 2018, 626, 1108–1120. [Google Scholar] [CrossRef] [Green Version]
  54. Shi, T.; Chen, Y.; Liu, H.; Wang, J.; Wu, G. Soil organic carbon content estimation with laboratory-based visible-near-infrared reflectance spectroscopy: Feature selection. Appl. Spectrosc. 2014, 68, 831–837. [Google Scholar] [CrossRef]
  55. Rossel, R.A.V.; Behrens, T. Using data mining to model and interpret soil diffuse reflectance spectra. Geoderma 2010, 158, 46–54. [Google Scholar] [CrossRef]
  56. Xiaobo, Z.; Jiewen, Z.; Povey, M.J.W.; Holmes, M.; Hanpin, M. Variables selection methods in near-infrared spectroscopy. Anal. Chim. Acta 2010, 667, 14–32. [Google Scholar] [CrossRef]
  57. Li, H.; Liang, Y.; Xu, Q.; Cao, D. Key wavelengths screening using competitive adaptive reweighted sampling method for multivariate calibration. Anal. Chim. Acta 2009, 648, 77–84. [Google Scholar] [CrossRef]
  58. Vohland, M.; Ludwig, M.; Thiele-Bruhn, S.; Ludwig, B. Determination of soil properties with visible to near- and mid-infrared spectroscopy: Effects of spectral variable selection. Geoderma 2014, 223–225, 88–96. [Google Scholar] [CrossRef]
  59. Wang, S.-L.; Li, J.; Fang, J. Predicting progression of als disease with random frog and support vector regression method. In Intelligent Computing Methodologies; Springer International Publishing: Cham, Switzerland, 2016; pp. 160–170. [Google Scholar]
  60. Pan, L.; Lu, R.; Zhu, Q.; Tu, K.; Cen, H. Predict compositions and mechanical properties of sugar beet using hyperspectral scattering. Food Bioprocess Technol. 2016, 9, 1177–1186. [Google Scholar] [CrossRef]
  61. Dong, J.; Guo, W.; Wang, Z.; Liu, D.; Zhao, F. Nondestructive determination of soluble solids content of ‘fuji’ apples produced in different areas and bagged with different materials during ripening. Food Anal. Methods 2016, 9, 1087–1095. [Google Scholar] [CrossRef]
  62. Guindo, M.L.; Kabir, M.H.; Chen, R.; Liu, F. Potential of vis-nir to measure heavy metals in different varieties of organic-fertilizers using boruta and deep belief network. Ecotoxicol. Env. Saf. 2021, 228, 112996. [Google Scholar] [CrossRef]
  63. Shen, T.; Kong, W.; Liu, F.; Chen, Z.; Yao, J.; Wang, W.; Peng, J.; Chen, H.; He, Y. Rapid determination of cadmium contamination in lettuce using laser-induced breakdown spectroscopy. Molecules 2018, 23, 2930. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Tawfik, W.; Askar, A. Study of the matrix effect on the plasma characterization of heavy elements in soil sediments. Prog. Phys. 2007, 3, 46–52. [Google Scholar]
  65. El Haddad, J.; Canioni, L.; Bousquet, B. Good practices in libs analysis: Review and advices. Spectrochim. Acta Part B At. Spectrosc. 2014, 101, 171–182. [Google Scholar] [CrossRef] [Green Version]
  66. Zhang, T.; Wu, S.; Dong, J.; Wei, J.; Wang, K.; Tang, H.; Yang, X.; Li, H. Quantitative and classification analysis of slag samples by laser induced breakdown spectroscopy (libs) coupled with support vector machine (svm) and partial least square (pls) methods. J. Anal. At. Spectrom. 2015, 30, 368–374. [Google Scholar] [CrossRef]
  67. Yan, C.; Qi, J.; Ma, J.; Tang, H.; Zhang, T.; Li, H. Determination of carbon and sulfur content in coal by laser induced breakdown spectroscopy combined with kernel-based extreme learning machine. Chemom. Intell. Lab. Syst. 2017, 167, 226–231. [Google Scholar] [CrossRef]
  68. Pořízka, P.; Klus, J.; Hrdlička, A.; Vrábel, J.; Škarková, P.; Prochazka, D.; Novotný, J.; Novotný, K.; Kaiser, J. Impact of laser-induced breakdown spectroscopy data normalization on multivariate classification accuracy. J. Anal. At. Spectrom. 2017, 32, 277–288. [Google Scholar] [CrossRef]
  69. Myakalwar, A.K.; Spegazzini, N.; Zhang, C.; Anubham, S.K.; Dasari, R.R.; Barman, I.; Gundawar, M.K. Less is more: Avoiding the libs dimensionality curse through judicious feature selection for explosive detection. Sci. Rep. 2015, 5, 13169. [Google Scholar] [CrossRef] [Green Version]
  70. Wang, Q.; Teng, G.; Qiao, X.; Zhao, Y.; Kong, J.; Dong, L.; Cui, X. Importance evaluation of spectral lines in laser-induced breakdown spectroscopy for classification of pathogenic bacteria. Biomed. Opt. Express 2018, 9, 5837–5850. [Google Scholar] [CrossRef]
  71. Menze, B.H.; Kelm, B.M.; Masuch, R.; Himmelreich, U.; Bachert, P.; Petrich, W.; Hamprecht, F.A. A comparison of random forest and its gini importance with standard chemometric methods for the feature selection and classification of spectral data. BMC Bioinform. 2009, 10, 213. [Google Scholar] [CrossRef] [Green Version]
  72. Duan, F.; Fu, X.; Jiang, J.; Huang, T.; Ma, L.; Zhang, C. Automatic variable selection method and a comparison for quantitative analysis in laser-induced breakdown spectroscopy. Spectrochim. Acta Part B At. Spectrosc. 2018, 143, 12–17. [Google Scholar] [CrossRef]
  73. Zhang, C.; Wang, Q.; Liu, F.; He, Y.; Xiao, Y. Rapid and non-destructive measurement of spinach pigments content during storage using hyperspectral imaging with chemometrics. Measurement 2017, 97, 149–155. [Google Scholar] [CrossRef]
  74. Sun, D.; Yang, F.; Su, M.; Han, W.; Dong, C. Quantitative analysis of cu in traditional chinese medicinal materials using laser-induced breakdown spectroscopy. Microw. Opt. Technol. Lett. 2022. [Google Scholar] [CrossRef]
  75. Zhu, C.; Lv, J.; Liu, K.; Li, Q.; Tang, Z.; Zhou, R.; Zhang, W.; Chen, J.; Liu, K.; Li, X. Fast detection of harmful trace elements in glycyrrhiza using standard addition and internal standard method–Laser-induced breakdown spectroscopy (sais-libs). Microchem. J. 2021, 168, 106408. [Google Scholar] [CrossRef]
  76. Liu, F.; Shen, T.; Kong, W.; Peng, J.; Zhang, C.; Song, K.; Wang, W.; Zhang, C.; He, Y. Quantitative analysis of cadmium in tobacco roots using laser-induced breakdown spectroscopy with variable index and chemometrics. Front. Plant Sci. 2018, 9, 1316. [Google Scholar] [CrossRef] [Green Version]
  77. Gondal, M.A.; Nasr, M.M.; Ahmed, M.M.; Yamani, Z.H.; Alsalhi, M. Detection of lead in paint samples synthesized locally using-laser-induced breakdown spectroscopy. J. Environ. Sci. Health Part A 2011, 46, 42–49. [Google Scholar] [CrossRef]
Figure 1. An illustration of the LIBS experimental setup.
Figure 1. An illustration of the LIBS experimental setup.
Foods 12 01125 g001
Figure 2. Gradient Boosting Machine Flow Chart.
Figure 2. Gradient Boosting Machine Flow Chart.
Foods 12 01125 g002
Figure 3. Normalized spectra of 12 Fritillaria varieties.
Figure 3. Normalized spectra of 12 Fritillaria varieties.
Foods 12 01125 g003
Figure 4. (a) Scatter plots of the full variables and reference value (b) Scatter plots of GBM combined with UVE feature selection method and reference value for the content of Cd.
Figure 4. (a) Scatter plots of the full variables and reference value (b) Scatter plots of GBM combined with UVE feature selection method and reference value for the content of Cd.
Foods 12 01125 g004
Figure 5. (a) Scatter plots of the full variables and reference value (b) Scatter plots of GBM combined with CARS feature selection method and reference value for the content of Cu.
Figure 5. (a) Scatter plots of the full variables and reference value (b) Scatter plots of GBM combined with CARS feature selection method and reference value for the content of Cu.
Foods 12 01125 g005
Figure 6. (a) Scatter plots of the full variables and reference value (b) Scatter plots of SVR combined with RF feature selection method and reference value for the content of Pb.
Figure 6. (a) Scatter plots of the full variables and reference value (b) Scatter plots of SVR combined with RF feature selection method and reference value for the content of Pb.
Foods 12 01125 g006
Figure 7. (a) Distributions of selected wave number by CARS (b) Distributions of selected wave number by RF (c) Distributions of selected wave number by UVE.
Figure 7. (a) Distributions of selected wave number by CARS (b) Distributions of selected wave number by RF (c) Distributions of selected wave number by UVE.
Foods 12 01125 g007
Table 1. Heavy metals (Cu, Cd and Pb) contents of Fritillaria obtained by ICP-MS.
Table 1. Heavy metals (Cu, Cd and Pb) contents of Fritillaria obtained by ICP-MS.
Heavy MetalGroupsCK1234567
Number3636363636363636
CuMin.1.5013.8033.0142.5153.7666.49155.41169.66
Max.4.8827.9955.9761.2869.3382.97242.04245.36
Mean2.5020.0637.3152.4865.9579.66172.57215.06
Range3.3714.1822.9518.7615.5716.4786.6375.69
Var.1.2013.1436.7137.6716.8720.33517.21773.97
Std1.093.626.056.134.104.5022.7427.82
CdMin.0.175.219.6420.8124.7344.8563.0683.66
Max.1.197.0811.7025.3429.3897.26117.89100.65
Mean0.435.8610.6422.8726.3960.3282.4893.91
Range1.021.872.054.524.6552.4154.8216.98
Var.0.070.210.312.341.72394.02175.3735.38
Std.0.280.460.561.531.3119.8513.245.94
PbMin.0.134.1221.2843.1358.6163.14102.16143.12
Max.0.767.3323.90102.0371.1590.51112.26219.07
Mean0.296.0622.7551.5865.6285.00112.26199.24
Range0.633.212.6158.9012.5327.3630.3375.95
Var.0.030.560.64265.5114.5559.5463.84604.43
Std.0.190.740.8016.293.817.717.9924.58
CK: Control group.
Table 2. Prediction results of the different models using full variables.
Table 2. Prediction results of the different models using full variables.
ModelsVariablesRC2RMSECRV2RMSEVRP2RMSEP
CdSVRFull variables0.99990.09950.895711.08740.827614.5054
PLSRFull variables0.99233.06370.679119.45390.549923.4405
GBMFull variables0.99841.36930.892411.26250.913910.2494
CuSVRFull variables0.99990.09950.960613.93590.928017.5096
PLSRFull variables0.99365.76910.690239.11380.562243.1827
GBMFull variables0.99793.28850.930818.48150.959613.1156
PbSVRFull variables0.99990.09960.930416.63080.887619.1794
PLSRFull variables0.99384.94430.692534.97910.627234.9311
GBMFull variables0.99980.84880.967311.39330.963510.9220
RC2: Coefficient of determination for calibration; RV2: Coefficient of determination for validation; RP2: Coefficient of determination for prediction; RMSEC: Root mean square error for calibration; RMSEV: Root mean square error for validation; RMSEP: Root mean square error for prediction.
Table 3. Prediction results of the different models using three variable selection (CARS, RF and UVE) methods.
Table 3. Prediction results of the different models using three variable selection (CARS, RF and UVE) methods.
ElementsMethodsNRC2RMSECRV2RMSEVRP2RMSEP
SVRCdCARS370.97075.99350.92359.49840.813915.0703
RF1510.99990.10000.92879.16710.93229.0933
UVE1920.95257.62770.906610.49290.911610.3842
CuCARS660.980310.18230.962213.64760.943015.5721
RF1200.98608.56990.964013.33190.965812.0658
UVE2310.971512.23950.941217.03240.964812.2429
PbCARS170.98697.19470.971010.72510.934114.6841
RF1240.99921.67070.973610.23230.968610.1224
UVE930.97998.91120.958512.83550.939514.0633
PLSRCdCARS370.93159.16590.909810.31320.906010.7075
RF1510.98254.63130.93139.00140.91829.9903
UVE1920.96596.46580.891011.33860.907210.6405
CuCARS660.975811.27460.952015.38410.956813.5625
RF1200.976911.02910.938517.42000.945715.1970
UVE2310.98568.70400.941117.04620.957513.4396
PbCARS170.968311.20390.971810.57980.959911.4563
RF1240.98447.84360.930316.65190.912916.8794
UVE930.972010.51750.955613.28320.938114.2331
GBMCdCARS370.99073.36490.900910.81130.914610.2061
RF1510.99821.47530.847413.41340.890911.5370
UVE1920.99671.98530.889911.39340.94038.5344
CuCARS660.99335.93320.931618.37790.966511.9356
RF1200.99296.09520.937117.62350.954513.9099
UVE2310.99644.31680.930418.53370.964812.2323
PbCARS170.99703.44340.946914.53350.942913.6631
RF1240.99822.61030.932916.33400.913616.8175
UVE930.99921.72480.956213.19670.960911.3113
N: Number of features selected; RC2: Coefficient of determination for calibration; RV2: Coefficient of determination for validation; RP2: Coefficient of determination for prediction; RMSEC: Root mean square error for calibration; RMSEV: Root mean square error for validation; RMSEP: Root mean square error for prediction.
Table 4. Results of some recent heavy metals (Pb, Cd and Cu) detection using LIBS of various samples and current work.
Table 4. Results of some recent heavy metals (Pb, Cd and Cu) detection using LIBS of various samples and current work.
Heavy-MetalSampleSpectral Line (nm)Reference
CdLettuce214.44, 226.50, 228.80[63]
Sargassum fusiforme441.56, 643.85[37]
Lipstick467.9, 573.80[23]
This work214.44, 226.50, 441.56, 467.90, 573.80, 643.85
CuTraditional Chinese medicinal materials324.79, 327.35[74]
Glycyrrhiza324.70[75]
Ligusticum wallichii324.46, 327.09[40]
Sargassum fusiforme324.75, 327.40[37]
Rice324.754, 327.396[48]
This work324.09, 324.79
PbPaint samples405.70[77]
Henna paste405.78[39]
Ligusticum wallichii405.80[40]
Medicinal herbs405.78, 404.00[1]
This work280.00, 404.00, 405.70
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kabir, M.H.; Guindo, M.L.; Chen, R.; Luo, X.; Kong, W.; Liu, F. Heavy Metal Detection in Fritillaria thunbergii Using Laser-Induced Breakdown Spectroscopy Coupled with Variable Selection Algorithm and Chemometrics. Foods 2023, 12, 1125. https://doi.org/10.3390/foods12061125

AMA Style

Kabir MH, Guindo ML, Chen R, Luo X, Kong W, Liu F. Heavy Metal Detection in Fritillaria thunbergii Using Laser-Induced Breakdown Spectroscopy Coupled with Variable Selection Algorithm and Chemometrics. Foods. 2023; 12(6):1125. https://doi.org/10.3390/foods12061125

Chicago/Turabian Style

Kabir, Muhammad Hilal, Mahamed Lamine Guindo, Rongqin Chen, Xinmeng Luo, Wenwen Kong, and Fei Liu. 2023. "Heavy Metal Detection in Fritillaria thunbergii Using Laser-Induced Breakdown Spectroscopy Coupled with Variable Selection Algorithm and Chemometrics" Foods 12, no. 6: 1125. https://doi.org/10.3390/foods12061125

APA Style

Kabir, M. H., Guindo, M. L., Chen, R., Luo, X., Kong, W., & Liu, F. (2023). Heavy Metal Detection in Fritillaria thunbergii Using Laser-Induced Breakdown Spectroscopy Coupled with Variable Selection Algorithm and Chemometrics. Foods, 12(6), 1125. https://doi.org/10.3390/foods12061125

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