Next Article in Journal
Data-Driven Artificial Intelligence Model of Meteorological Elements Influence on Vegetation Coverage in North China
Previous Article in Journal
Susceptibility Prediction of Post-Fire Debris Flows in Xichang, China, Using a Logistic Regression Model from a Spatiotemporal Perspective
Previous Article in Special Issue
A NIRS-Aided Methodology to Elucidate the Nutrition of the Endangered Mountain Gazelle (Gazella gazella) Using Samples of Rumen Contents from Roadkills
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Near-Infrared Spectroscopy and Mode Cloning (NIR-MC) for In-Situ Analysis of Crude Protein in Bamboo

1
Department of Biochemistry, Molecular Biology, Entomology, and Plant Pathology, Mississippi State University, Starkville, MS 39762, USA
2
Department of Wildlife, Fisheries & Aquaculture, Mississippi State University, Starkville, MS 39762, USA
*
Author to whom correspondence should be addressed.
Submission received: 1 January 2022 / Revised: 22 February 2022 / Accepted: 4 March 2022 / Published: 8 March 2022
(This article belongs to the Special Issue Near Infrared Spectroscopy in Animal Ecophysiology)

Abstract

:
This study develops Near-Infrared Spectroscopy (NIRS) and Mode-Cloning (MC) for the rapid assessment of the nutritional quality of bamboo leaves, the primary diet of giant pandas (Ailuropoda melanoleuca) and red pandas (Ailurus fulgens). To test the NIR-MC approach, we evaluated three species of bamboo (Phyllostachys bissetii, Phyllostachys rubromarginata, Phyllostachys aureosulcata). Mode-Cloning incorporated a Slope and Bias Correction (SBC) transform to crude protein prediction models built with NIR spectra taken from Fine–Ground leaves (master mode). The modified models were then applied to spectra from leaves in the satellite minimal processing modes (Course–Ground, Dry–Whole, and Fresh–Whole). The NIR-MC using the SBC yielded a residual prediction deviation (RPD) = 2.73 and 1.84 for Course–Ground and Dry–Whole sample modes, respectively, indicating a good quantitative prediction of crude protein for minimally processed samples that could be easily acquired under field conditions using a portable drier and grinder. The NIR-MC approach also improved the model of crude protein for spectra collected from Fresh–Whole bamboo leaves in the field. Thus, NIR-MC has the potential to provide a real-time prediction of the macronutrient distribution in bamboo in situ, which affects the foraging behavior and dispersion of giant and red pandas in their natural habitats.

1. Introduction

Forage grasses and legumes supply energy and nutrients to various ruminants and non-ruminants by providing the protein, vitamins, and minerals necessary for growth and metabolism [1]. Near-infrared spectroscopy (NIRS) has been used widely in assessing the forage quality of livestock, and in recent years, an increase in studies of forage nutrition for wild species [2,3,4]. Importantly, the determination of the nutritional value of forage grasses can be used to understand the nutritional drivers of wild-animal movement in situ [5]. For such studies, miniaturized NIR spectrometers designed for in-field spectral collection are more practical for real-time assessment than bench-top systems [6,7,8].
NIRS is a cost-effective and rapid technique for measuring the nutritional quality of many forage grasses, including bamboo [9,10]. The quantitative nutrient composition of forages can be obtained using NIRS combined with mathematical transformations to generate calibration models for nutrient prediction [11]. However, sample processing emphasizes different contributions to the spectral profile of the sample, which can then affect the model performance and prediction results. In standard practice, NIR spectra are collected from forage samples that are entirely dried at temperatures from 55 to 70 °C because samples with high moisture generate a strong water absorption signal that overlaps and obscures other spectral features, which causes peak shifts and results in non-linear responses [1,12,13]. Several studies in soil and forage analysis have found that NIR results from wet samples have higher standard errors of calibration compared with dried samples and hence have lower prediction accuracy [13,14]. Moreover, the size, shape, and uniformity of particles in a sample can change the surface light scattering and affect the slope of the NIR spectra [15], as well as cause absorption flattening [16,17]. Under laboratory conditions, isotropic and uniform sized particles can be obtained by grinding in a cyclone mill to minimize scattering noise in the spectrum. The grinding process also increases homogenization and provides the highest accuracy in prediction equations from spectral data [18].
Although NIR spectra and model prediction of nutrient quantities from dried, highly ground samples can be more accurate, the sample processing steps require time, labor, and unwieldy equipment that is unsuitable for in situ applications. Field studies require a simplified procedure for minimal or unprocessed fresh samples. The calibration transfer of a master prediction model generated for fine-ground samples may be applicable to non-ground or even wet samples [19]. However, in some cases, the model cannot be applied directly to new datasets because the gaps between the existing and new types of samples cause variation between the spectra, leading to biased predictions. For example, different processing methods of the same sample can cause a change in the physical properties of the sample, such as particle size or surface texture [20]. Here, we propose an NIR-Mode Cloning (NIR-MC) approach that consists of calibration transfer based on mathematical and statistical techniques to correct for differences caused by various processing modes [21]. One method commonly used for correcting predictions is the simple univariate Slope and Bias Correction (SBC), in which the prediction model is developed based on the Master mode and can then be modified to correct predicted values in a satellite mode [22,23]. The SBC model transfer was applied to the leaves of bamboo, which is a wild forage grass for numerous species. The predictive performance of the SBC model transfer has been shown to be similar to that of other spectral correction methods such as direct standardization (DS), piecewise direct standardization (PDS), and external parameter orthogonalization (EPO) corrections [24,25,26].
Bamboo is the primary food source for the giant panda (Ailuropoda melanoleuca) and red panda (Ailurus fulgens), representing 99% and 91.4–99.1% of their diet, respectively [27,28,29]. Bamboos are distributed throughout the panda habitat and are cold-tolerant and evergreen [30,31], serving as the primary food source throughout the year, especially in the winter. The amount of protein and fat that either panda species can acquire from a diet of bamboo is limited because bamboo is a low-nutrition/energy food source with crude protein composing only 7.0–21.6% of total dry matter and crude fat only 1.6–4.2%, with the remaining content composed of complex carbohydrates and digestible sugars [9,32]. Giant pandas and red pandas have a relatively short gastrointestinal tract (GIT), typical of carnivores, and the low abundance of cellulases and endo-hemicellulases in the carnivore digestive system make it difficult to digest cellulose into usable nutrients or to extract the protein and fat constituents [31,32,33]. The digestive efficiency of bamboo leaves is less than 30% in giant pandas and red pandas [34,35], resulting in the need to consume large quantities of bamboo each day.
The variation in nutrient content across different bamboo species and habitats can affect giant and red panda diet quality and dietary selection. For example, studies have shown that captive giant pandas prefer culm over leaves in spring due to the seasonal changes in carbohydrate distribution within different parts of the bamboo plant [36]. Moreover, giant pandas shift between consuming bamboo leaves and shoots of varying bamboo species throughout the year, as the concentrations of primary energy constituents and critical elements such as calcium, phosphorus, and nitrogen shift within the plant structure, potentially influencing annual patterns of migration [37]. Red pandas also prefer bamboo leaves in the summer and fall, as leaves contain the highest levels of crude protein and fat during these seasons [32]. A better understanding of the bamboo foodscape can assist wild-population management through habitat conservation and the establishment of protected areas to account for the seasonal foraging behavior, migration, and population distribution changes. The direct assessment of forage quality in situ for giant and red panda population management is possible by using portable NIRS and sampling processing devices. To test the NIR-MC approach, we evaluated three species of bamboo from the genus Phyllostachys that are endemic to China and that serve as the primary forage of both panda species. Once the calibration models for predicting nutritional parameters were built using spectra from highly processed samples, the NIR-MC method was applied to fit spectra from minimally processed bamboo leaves with calibration transfer and mathematical correction. Building predictive models using SBC transforms that adjust for sample processing modes and are transferable to spectra collected from bamboo leaves in situ willenable real-time assessment of the bamboo nutritional landscape of panda habitat.

2. Materials and Methods

2.1. Bamboo Sampling and NIR Spectra Collection

The leaves from three bamboo species, Phyllostachys bissetii (Bissetii), Phyllostachys rubromarginata (Rubro), and Phyllostachys aureosulcata (Sulcata), were selected and collected bi-weekly over 20 months from established cultivars at the Memphis Zoo Bamboo Farm (Agricenter, Memphis, TN, USA). Fresh live leaves (n = 10) were randomly selected from different rows throughout the plot for each species and NIR spectra (n = 10/ scans per leaf) were taken in the field from both sides of each leaf to create the Fresh–Whole mode dataset (Figure 1A) using an ASD FieldSpec®3 Vis-NIRS portable spectrometer (350–2500 nm, 50 scans with a 34 ms integration time, approximately 1 nm resolution) and a 20 mm in diameter ASD low-intensity ‘plant probe’ specially designed for heat-sensitive targets such as plant tissues. The selected leaves were collected and dried in an air-forced oven at 60 °C to constant weight for 48 h and pressed for NIR spectra (n = 10 scans/per leaf) to generate the Dry–Whole leaf mode dataset (Figure 1B). Next, dried leaves were coarsely ground with a Wiley mill (Figure 1C) and then finely ground into 0.1 mm isotropic particles with a cyclone mill (Figure 1D, Udy Corporation, Fort Collins, CO, USA). The NIR spectra of dried and ground samples were then collected (n = 3 replicates/ sample). In the standard procedures of this study, each ground sample mode was scanned three times as replicates, as those samples were homogeneous and in small volumes. However, extra scans were collected from different parts of the whole leaves of Dry–Whole and Fresh–Whole sample modes (n = 10 replicates/sample) to cover the structural variation of the tissue in the whole bamboo leaf. All the samples and spectra were classified into four processing modes: the Fine–Ground samples are considered the Master mode, while the minimally processed Coarse–Ground, Dry–Whole, and Fresh–Whole are considered satellite modes. The Fresh–Whole samples are unprocessed, while Dry–Whole and Coarse–Ground samples can be simply processed by a portable dryer and coffee grinder, which are practical in situ (Figure 2). Cyclone milling to obtain the Fine–Ground samples is not feasible for in-field applications as the equipment is not portable, and the samples are easily contaminated or lost during processing.

2.2. Chemical Analysis

The ground samples were sent for proximate analysis to generate the chemical reference data needed to calibrate the NIR Master mode (Forage Labs, Cumberland, MD, USA). The dry matter (DM) content was determined after the ground subsamples were weighed and dried at 100 °C for 24 h. Crude protein (CP) is the total protein and was determined using the micro-Kjeldahl method [38]. The percentage of nitrogen was calculated using a Leco Nitrogen Analyzer (model FP-2000, Leco Corporation, St. Joseph, MI, USA) (AOAC, 1995), and the percentage of CP was obtained by multiplying the percentage of nitrogen by a standardized factor of 6.25 as protein averages approximately 16 % nitrogen. In addition, a randomized set of samples (n = 10) were split and measured separately to calculate the standard error of the laboratory values (SEL) [39].

2.3. Chemometric Analysis

2.3.1. Calibration Models

Before the application of chemometrics, an outlier detection analysis was performed, in which outlier samples were removed from the database. The spectra collected from each leaf under the conditions of each sample processing step were averaged and arranged into four balanced databases (Table 1). The databases were equally weighted by random sampling within each species for calibration and internal validation. Spectra in Database 1 contained equal numbers of samples from the three bamboo species and were randomly divided into calibration (23 samples for each bamboo species) and internal validation sets (6 samples for each species) comprising 80% and 20% of the original samples, respectively. In databases 2–4, spectra from two bamboo species were combined such that 80% of samples were randomly selected to establish calibration models (n = 23 for Rubro and Sulcata, n = 31 for Rubro and Bissetii, n = 23 for Bissetii and Sulcata), and the remaining 20% of samples were used as internal validation sets (n = 47 for Bissetii, n = 29 for Sulcata, n = 39 for Rubro). Spectra collected from the third species were used as the external validation set, with 47 samples from Bissetii, 29 samples from Sulcata, and 39 samples from Rubro in Databases 2–4, respectively. The mathematical pre-treatments of Standard Normal Variate (SNV) with de-trending and a second derivative (Savitzky–Golay smoothing, symmetric points = 16, smoothing points = 33) were applied on the full data spectrum (350–2500 nm). Principal Component Analysis (PCA) and wavelength-range evaluation were applied to the databases to determine dominant peaks in the loadings using the Unscrambler® X v.11, Aspen Technology Inc, MA, USA. The Partial Least Squares Regression (PLSR) algorithm for each processing mode was applied on the transformed spectra associated with the crude protein information (1000–2500 nm) from the balanced databases to predict the crude protein values. The models were calibrated and validated using the calibration set and internal validation set, respectively, while the prediction results were obtained from the external validation set.

2.3.2. NIR Mode Cloning (NIR-MC) Using Slope and Bias Correction (SBC)

The NIR-MC was carried out in Unscrambler® X v.11 (Aspen Technology Inc, MA, USA). The best calibration model for the master (Fine–Ground) mode to predict the crude protein values was applied to the three satellite sampling modes. Additionally, models built with samples from the Coarse–Ground mode were applied to the satellite modes Dry–Whole and Fresh–Whole, and the best-fit calibration model for the Dry–Whole mode was applied to the Fresh–Whole samples (Figure 2). The SBC was performed on predicted Y-values only, without correcting the spectral matrix. The correction of the predicted values and model evaluation was performed using R studio (Version 1.2.5019, RStudio, Inc., Boston, MA, USA). With the SBC, it is assumed that a linear relationship exists between the predictions of spectra measured from the master (Fine–Ground) mode samples and the predictions obtained from the satellite (Coarse–Ground) samples using the same PLSR model built on the Fine–Ground mode. The SBC was conducted using the external validation set in Database 2 containing the predicted CP values of Fine–Ground and Coarse–Ground samples using the same Fine–Ground PLSR model. Database 2 was divided into two subsets, a training set and a test set, comprised of 70% and 30% of the dataset, respectively. With ordinary least squares, the training set was used to build a linear regression between the predictions of the Fine–Ground samples and the predictions of the Coarse–Ground samples using the same Fine–Ground model:
yfg = bias + slope ∗ ycg
where yfg is the predicted CP values of Fine–Ground samples using the Fine–Ground calibration model, and ycg is the predicted CP values of Coarse–Ground samples using the Fine–Ground calibration model. Slope and bias coefficients were computed using a linear regression between those two sets of predictions.
Then, in the test set, the predicted values of Coarse–Ground samples were corrected by the slope and bias (intercept) based on the linear equation:
ycg. corr = bias + slope ∗ ycg
where ycg. corr is the corrected predictions of Coarse–Ground samples.
The same procedures were applied to the mode cloning from the Fine-Ground model to Dry–Whole and Fresh–Whole samples, from the Coarse–Ground model to Dry–Whole and Fresh–Whole samples, and from the Dry–Whole model to Fresh–Whole samples (Figure 2). The root mean square error of prediction (RMSEP), standard error of prediction (SEP), bias, the coefficient of determination (R2), and the residual prediction deviation (RPD) were used to summarize the accuracy of the prediction. The RPD is the standard deviation (SD) of the sample’s reference values divided by the root mean square error of prediction [40]. The RPD value is used as the criteria for evaluating model performance, a model with an RPD value higher than 1.8 is considered a good model for quantitative prediction of a constituent such as crude protein, and RPD > 2.5 indicates an excellent model [41].

3. Results

3.1. Sample Statistics

The CP content of all the bamboo samples varied from 11.6% to 22.8%, with an average of 17.1% and an SD of 2.1%, determined by the chemical analyses. The frequency distribution plots of CP values are shown in Figure 3. The histogram shows a normal distribution of CP values overlapping across the three bamboo species. The histogram also shows a broad range of CP values was available to build reliable calibration models for predicting CP values in bamboo samples. The Tukey’s test (p > 0.05) showed no significant difference between species in the means. The standard error of laboratory (SEL) was 0.49, calculated using the technical duplicates of the same samples.

3.2. Spectral Characteristics

Figure 4 shows the averaged raw and transformed spectra collected from the four evaluated processing modes, grouped by species. Several major bands can be observed in the raw spectra. For Fresh–Whole samples, the two prominent peaks appeared near the wavelength ranges of 1460–1500 nm and 1850–1950 nm are the dominant and broad near-IR absorption bands of water [42,43]. The major absorption bands of protein can be observed near the wavelength ranges of 1640–1680 nm, 2050–2100 nm, and 2290–2390 nm, which are assigned to the amide A-amide II combination and high-frequency aliphatic CH stretching bands [44,45].
The transformation helped accentuate the spectral response to chemical changes in the samples. Transformed spectra show that the absorbance at around 1425 (first overtone of O-H and N-H in amino and amide groups), 1650–1680 (CH), 1900 (NH), 2050–2060 (ROH/NH), and 2060 (carbonyl stretch of the primary amide) nm have more substantial peaks as the protein content increased [45,46]. These protein signals would be an important factor in building the calibration models for determining CP content in bamboo samples but are masked by the water signal in Fresh-Whole samples. Thus, the overlapping water signals would have to be factored out in transferring calibration models from dry samples to fresh samples.
The raw and transformed spectra of the three dried modes (Fine–Ground, Coarse–Ground, and Dry–Whole) shared similar major absorption bands. Considering that there is no significant difference between bamboo species in their mean reference CP values suggests that it is possible that SBC without spectra correction between processing modes can provide acceptable results in mode cloning.

3.3. NIR Modeling across Species

The PCA scores plots and loadings for all the processing modes are shown in Figure 5 and Figure 6. In the Fine–Ground (Figure 5A), the first two PCs account for 73% of the variance when the three species were not well distinguished, indicating different bamboo species exhibit overlapping confidence intervals. In this case, the PCA loadings (Figure 5B) showing the prominent peaks that explained the trends in the scores plot accounted for PC-1 = 59%, PC-2 = 14% of the variance Similarly, in the Coarse–Ground plot (Figure 5C), the three species of bamboo overlap with PCA loadings explaining 58% and 14% of the variation for PC-1 and PC-2, respectively (Figure 5D). In these two ground modes, the PCA loadings plots (Figure 5B–D) show that wavelengths from the regions 1400–1425, 1900–1950, and 2310–2340 nm are essential for the first and second PC loading values, which are the dominant peaks influencing the trends in the prediction of crude protein. Likewise, in the two whole modes, the scores from the three bamboo species overlap (Figure 6A–C), and the regions 1900–1950 and 2310–2340 nm in the loadings are shown as crucial for PC-1 and PC-2 values (Figure 6B–D).
The statistical summary of the calibration models for predicting CP across species is shown in Table 2, and contains RMSEC, RMSECV, RMSEP, SEC, SECV, SEP, R2, and RPD values. The main goal was to test if the calibration models based on samples from two species could be directly transferred to the third species of the same genus. For the Fine-Ground, all three calibration models (Databases 2–4) transferred to the third species well, with R2 ranging from 0.82 to 0.99, and RPD ranging from 2.14 to 12.8 in the calibration, cross-validation, validation, and external validation sets, indicating an excellent quantitative predictive performance. The Dry–Whole and Fresh–Whole in Database 4 provided a very poor prediction of CP values, with RPD values lower than 1.0. When Database 4 was excluded, it showed that Coarse–Ground models had good performance in quantitative prediction of CP contents (0.76 < R2 < 0.98, 1.86 < RPD < 9.10) in the calibration, cross-validation, validation, and external validation sets. Though the Dry–Whole and Fresh–Whole models still performed well in the calibration and cross-validation set (0.71 < R2 < 0.92, 1.78 < RPD < 3.53; 0.83 < R2 < 0.97, 2.50 < RPD < 6.19, respectively), they failed to provide good performance in the validation and external validation set. The Dry–Whole models had fair performance in predicting CP values (0.25 < R2 < 0.75, 1.04 < RPD < 2.00) in the validation and external validation sets, whereas the Fresh–Whole models had less accuracy (0.21 < R2 < 0.39, 1.09 < RPD < 1.13). The ratio of RMSEP and RMSECV (RMSEP/RMSECV) in databases 1–3 ranged from 0.46 to 1.60, 0.78 to 3.13, 0.94 to 2.13, and 1.38 to 3.06 in Fine–Ground, Coarse–Ground, Dry–Whole, and Fresh–Whole mode, respectively. An RMSEP/RMSECV ratio larger than 1.2 might suggest that the acquired calibration set is not large enough to properly capture the variance of the population [47].
Thus, it can be concluded that, for the three bamboo species used in this study, the cross-species model transfer can be applied directly without performing any standardization procedure in Databases 2 and 3. However, when transferring the calibration model built for Bissetii and Sulcata samples to Rubro samples in Database 4, only the Fine–Ground and Coarse–Ground modes performed well enough for quantitative prediction of CP, while the two whole-leaf modes all had poor prediction results in the validation set. For further testing of NIR-MC across processing modes, Database 4 was excluded and since Database 2 had higher RPD values in calibration, cross-validation, and validation sets in general, compared with Database 3, it was selected to apply NIR-MC in the following steps.

3.4. NIR-MC and SBC across Processing Modes

The prediction statistics of CP values after applying NIR-MC and SBC on Database 2 are shown in Table 3 and Figure 7 and Figure 8. Figure 7 shows that all the transferred models after SBC had lower RMSEP values and increased RPD values than models without SBC. The improved model performances after SBC indicated the necessity of standardizing predicted values in mode cloning. Table 3 shows that, when applying the Fine–Ground Master model to Coarse–Ground bamboo samples, the prediction results were excellent before (RMSEP = 0.83, RPD = 2.62) and after SBC (RMSEP = 0.79, RPD = 2.73). The RPD values increased by 0.74 and 1.07 after SBC when applying the Fine–Ground model to Dry–Whole and Fresh–Whole samples. When transferring the Coarse–Ground model to Dry–Whole and Fresh–Whole samples, RPD values increased after SBC. When applying the Dry–Whole model to Fresh-Whole samples, the prediction results also improved, with the RPD value increasing after SBC.
By comparing the performances of the original model generated for each processing mode and the models using the NIR-MC approach (Figure 8 and Table 3), we found the best-fit model for each sampling process. For Coarse–Ground samples, the transferred Fine–Ground model after SBC (RPD = 2.73) had a better performance than the original Coarse–Ground model (RPD = 2.65). For Dry–Whole samples, the transferred Fine–Ground model after SBC (RPD = 1.84) provided better prediction than either the original Dry–Whole model or the transferred Coarse–Ground model after SBC (RPD = 1.73, 1.32). For Fresh–Whole samples, the transferred Dry–Whole model after SBC (RPD = 1.49) was better than the original Fresh–Whole model, the transferred Coarse–Ground model after SBC, and the transferred Fine–Ground model after SBC (RPD = 1.13, 1.16, 1.22). Thus, the transferred Fine–Ground models after SBC were the best-fit models to predict CP values of samples of Coarse–Ground and Dry–Whole samples, while the transferred Dry–Whole models after SBC had the best performance in Fresh–Whole samples.

4. Discussion

The NIR-MC approach was used across bamboo species and sample processing modes to determine the CP content in bamboo leaves. The model transfer, and Slope and Bias Correction (SBC), were applicable to a model generated from Fine–Ground samples to predict CP values of less-processed bamboo samples across species. Constructing a master calibration model based on an extensive database containing different bamboo species and Fine–Ground samples under laboratory conditions that are transferable through NIR-MC allows a fast and accurate CP prediction of new bamboo samples under various circumstances and eliminates the need for creating additional models specific to sample processing states. Currently, NIR mode cloning or calibration transfer techniques are primarily applied to address variation in sample physical or chemical composition due to batch differences, instrumental changes, or environment change over time [22]. Only a few investigators have performed mode cloning from ground samples to non-milled samples using spectral correction methods [19]. This study demonstrates that the simple SBC, without complex spectra correction, was capable of achieving good prediction accuracy when transferring the calibration model from Fine–Ground to Coarse–Ground and even Dry–Whole forage samples. The target subject matter of bamboo leaves makes our study unique and applicable to ecophysiology as this is the first time the NIR-Mode cloning method has been used as a rapid tool in assessing the CP contents of the primary food sources of the giant panda and red panda.
We showed that NIR could be applied directly across bamboo species from the same genus. Phyllostachys is one of the easiest bamboo genera to identify and can be defined by its taxonomic and morphological characteristics. For example, compared with the tropical Bambusa group, Phyllostachys species are less fibrous [48]. In leaf anatomy, they lack fusoid cells, adaxial ribs, expanded leaf margins, and well-developed abaxial papillae and prickles [49]. Due to the similarities in leaf morphology and chemical composition, the spectra patterns between species are similar, showing no significant separation in the PCA. Furthermore, the distribution and means of CP values do not show significant variance. The similarities in spectra and reference data enable us to use calibration models across species, even in less-processed samples. However, when applying the same models to Pseudosasa japonica, which is from a different genus, all the transferred models showed poor performance in predicting CP values (RPD ≤ 1.54, R2 ≤ 0.80). We found that the mean amount of CP in Pseudosasa japonica (15.32% DM) was significantly different (p < 0.05) than in the Phyllostachys species (16.35, 17.19, 17.31% DM), which would cause an error in the calibration if this species were added into the spectral databases. A micromorphological study has revealed distinct characteristic variations in the leaf surface between the two genera as well. For example, there are no papillae around the stomata in Phyllostachys species, while in Pseudosasa japonica, the stomata are surrounded by 4–10 elongated papillae [50]. The differences in micromorphological features might lead to spectra variation between the two genera, affecting the model performance for whole leaves. Therefore, we need a different calibration model for Pseudosasa or other genera in future studies.
When considering the RPD value as the criterium for evaluating model performance, a model with an RPD value higher than 1.8 indicates a good model for quantitative prediction of a constituent such as crude protein, and RPD > 2.5 indicates an excellent model [41]. We found that, with the NIR-MC method, the master Fine–Ground model is the best-fit model for all three satellite modes. It can also provide more-accurate predictions than using the original satellite models directly. The RPD values were higher than 2.5 when applying the Fine–Ground model to both Fine–Ground and Coarse-Ground samples with NIR-MC after SBC. This demonstrated that the cloned model has a similar prediction ability in predicting CP contents using Coarse–Ground samples as the laboratory-standard Fine–Ground samples. The Fine–Ground model transferred and standardized to Dry–Whole samples had an RPD value larger than 1.8, indicating that this model met the standard for quantitative prediction. None of the models applied to Fresh–Whole samples were good enough for quantitative predictions, with RPD < 1.8 (Figure 7).
In this study, we have established a database containing three bamboo species and four levels of processing. Although the Fine–Ground model applied to Fine–Ground samples yields the best accuracy in predicting protein content, it is time-consuming and requires a laboratory cyclone mill grinder, which is not amenable to field studies. We also show that the Fine–Ground Master model has good performance in predicting protein content of the satellite modes Coarse–Ground and Dry–Whole. Therefore, the Fine–-Ground model obtained from the existing databases built in the laboratory can be applied to new field samples of satellite processing modes by NIR-MC using SBC. According to various environmental and practical situations, we can choose between processing modes to obtain predictive results of CP content. The drying process is required to obtain good predictions of CP values as the water signal overlaps that of protein. However, Dry–Whole bamboo leaves can be acquired by drying with a portable drier or natural sunlight. To obtain higher-quality predictions of CP values with similar prediction ability as the laboratory standard, grinding the dried samples with a portable coffee grinder is necessary. This Coarse–Ground processing can be performed in situ using dried samples.
The differences between the evaluated processing modes affected the performance of the NIR-MC. The calibration model transferred from Fine–Ground to Coarse–Ground mode had the best performance of predicting CP values because Coarse–Ground samples were similar to Fine–Ground samples in particle size and distribution. The differences between Dry–Whole and Fine–Ground modes were more significant as the Dry–Whole samples were not homogenous and still maintained the original physical structure of the bamboo leaves, which can cause more significant variation in spectra dependent on where and how the spectra are collected from the leaves. However, all the dry modes (Fine–Ground, Coarse–Ground, and Dry–Whole) shared a similar spectral pattern and similar major absorption bands. The similarity of spectra patterns presented in the three dried modes indicated that a simple univariate correction of slope and bias is enough without the spectral correction. The moisture in Fresh–Whole samples caused the variation of spectra, particularly in the shifts of some absorption bands related to water, thus limiting the application of NIR-MC by having strongly biased predictions. To remove the effects of moisture from NIR spectra, a mathematical algorithm is needed to project all the spectra of wet samples orthogonally to the space of unwanted variation [51]. We propose that external parameter orthogonalization (EPO) or a model transfer tool of combining direct orthogonal signal correction with slope and bias correction (DOSC–SBC) [52] can be applied on Fresh–Whole samples to remove the spectral variation caused by water before using the NIR-MC.
We chose to analyze the CP of bamboo leaves as a nutritional parameter of the diets of giant and red pandas because it affects their forage behavior and seasonal movement. Both species maintain a carnivorous digestive tract but have evolutionarily adapted to a diet of bamboo; as such, these two foragers depend on the changes in the nutrient content distribution in the bamboo plant based on species and season [32,36]. Dietary proteins that giant and red pandas can obtain from bamboo leaves are limited because the bamboo leaf is a high-fiber and low-nutrition food. To meet this nutritional requirement, during the adaptation to becoming bamboo specialists, specific genes helped by increasing the efficiency of releasing lysine and arginine from dietary proteins and amino acid recycling, thus offsetting the limited protein content in the bamboo diet [53]. Despite the low protein content (7.0% to 17.9%) in bamboo leaves, energy utilization of giant pandas is sourced primarily from protein (61%), with the remainder split between carbohydrates (23%) and fats (16%) [54]. Red pandas also have a foraging behavior to increase protein intake by selecting certain bamboo species with higher protein content [55]. The macronutrient composition of giant and red pandas’ diets with high dietary protein-energy is similar to other carnivore diets, showing the importance of assessing the CP values of bamboo as an indicator of their diet quality. The NIR-MC technique provides real-time assessment of protein values without the time-consuming procedures of chemical analysis or sample processing, potentially allowing rapid in-field decisions for the management of panda habitat based on forage structure and nutrient patterns.
Crude protein is not the only necessary nutrient for giant and red pandas. Previous studies found that the changes in carbohydrates (culm starch and bound glucose) of bamboo caused the changes in giant pandas’ foraging behavior during spring [36]. The changes in calcium, phosphorus, and nitrogen concentrations also affected giant panda annual migration patterns as they switch their diet between the leaves and shoots of bamboo species. The NIR-MC is not limited to crude protein analysis, nor is it limited to analysis of bamboo, but has potential for application of other nutrient parameters and other forage systems, such as that of the eucalyptus and koala [56].

5. Conclusions

In this study, the viability of the NIR-MC approach was examined in predicting the crude protein values of bamboo leaves. With databases containing different bamboo species and sample processing modes, we demonstrated the possibility of using NIR across different bamboo species sharing similar leaf morphology. We also demonstrated the potential of NIR-MC in transferring the Master mode to less-processed satellite modes. For future studies in the evaluation of bamboo nutrition, with the help of NIR-MC and SBC, new spectra of field samples processed with a portable dryer and grinder could be fitted into the existing database of the master model, allowing a real-time prediction of diet quality available to the giant panda and red panda.

Author Contributions

Conceptualization, C.K.V.; methodology, Q.S., X.O. and C.K.V.; formal analysis, Q.S., C.K.V. and M.S.-R.; validation, Q.S.; investigation, Q.S., X.O. and C.K.V.; resources, A.J.K. and C.K.V.; data curation, Q.S., M.S.-R., X.O. and C.K.V.; writing—original draft preparation, Q.S.; writing—review and editing, Q.S., M.S.-R., A.J.K. and C.K.V.; visualization, Q.S. and C.K.V.; supervision, A.J.K. and C.K.V.; project administration, C.K.V.; funding acquisition, A.J.K. and C.K.V. All authors have read and agreed to the published version of the manuscript.

Funding

This project was supported by a USDA U.S. Forest Service International Programs cooperative agreement, grant # 18-DG-11132762-248, on ‘Coupled Forest-Human Ecosystems: Connecting Forest Management with Conservation Initiatives’, the U.S. Department of Agriculture, Agricultural Research Service, Biophotonics project # 6066-31000-015-00D, and the Mississippi Agricultural and Forestry Experiment Station, the National Institute of Food and Agriculture, U.S. Department of Agriculture, Hatch project under accession number W3173.

Data Availability Statement

Data are available upon request from Carrie K. Vance.

Acknowledgments

We thank the giant panda staff of the Memphis Zoo for providing information about the bamboo species in the giant pandas’ diet.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Van Soest, P.J. Nutritional Ecology of the Ruminant; Cornell University Press: Ithaca, NY, USA, 1994; Volume 32, p. 7. [Google Scholar]
  2. Showers, A.; Scott, E.; Douglas, R.; Jerry, W.; James, C.; Ben, H.; Showers, S.E.; Tolleson, D.R.; Stuth, J.W.; Kroll, J.C.; et al. Predicting Diet Quality of White-Tailed Deer via NIRS Fecal Profiling Predicting Diet Quality of White-Tailed Deer via NIRS Fecal Profiling. Rangel. Ecol. Manag. 2021, 59, 300–307. [Google Scholar] [CrossRef]
  3. Dixon, R.; Coates, D. Near infrared spectroscopy of faeces to evaluate the nutrition and physiology of herbivores Review. J. Near Infrared Spectrosc. 2009, 17, 1–31. [Google Scholar] [CrossRef]
  4. Rothman, J.M.; Chapman, C.A.; Hansen, J.L.; Cherney, D.J.R.; Pell, A.N. Rapid assessment of the nutritional value of foods eaten by mountain gorillas: Applying near-infrared reflectance spectroscopy to primatology. Int. J. Primatol. 2009, 30, 729–742. [Google Scholar] [CrossRef]
  5. Givens, D.I.; De Boever, J.L.; Deaville, E.R. The principles, practices and some future applications of near infrared spectroscopy for predicting the nutritive value of foods for animals and humans. Nutr. Res. Rev. 1997, 10, 83–114. [Google Scholar] [CrossRef]
  6. Beć, K.B.; Grabska, J.; Huck, C.W. Principles and Applications of Miniaturized Near-Infrared (NIR) Spectrometers. Chem.—A Eur. J. 2021, 27, 1514–1532. [Google Scholar] [CrossRef]
  7. Goi, A.; Simoni, M.; Righi, F.; Visentin, G.; De Marchi, M. Application of a handheld near-infrared spectrometer to predict gelatinized starch, fiber fractions, and mineral content of ground and intact extruded dry dog food. Animals 2020, 10, 1660. [Google Scholar] [CrossRef]
  8. Beć, K.B.; Grabska, J.; Plewka, N.; Huck, C.W. Insect protein content analysis in handcrafted fitness bars by NIR spectroscopy. Gaussian process regression and data fusion for performance enhancement of miniaturized cost-effective consumer-grade sensors. Molecules 2021, 26, 6390. [Google Scholar] [CrossRef]
  9. Wiedower, E.; Hansen, R.; Bissell, H.; Ouellette, R.; Kouba, A.; Stuth, J.; Rude, B.; Tolleson, D. Use of near infrared spectroscopy to discriminate between and predict the nutrient composition of different species and parts of bamboo: Application for studying giant panda foraging ecology. J. Near Infrared Spectrosc. 2009, 17, 265–273. [Google Scholar] [CrossRef]
  10. García, J.; Cozzolino, D. Use of near infrared reflectance (NIR) spectroscopy to predict chemical composition of forages in broad-based calibration models Uso de la espectroscopía de reflectancia en el infrarrojo cercano (NIR) para predecir la composición química de forrajes e. Agric. Tec. 2006, 66, 41–47. [Google Scholar]
  11. Deaville, E.R.; Flinn, P.C. Near infrared (NIR) spectroscopy: An alternative approach for the estimation of forage quality and voluntary intake. In Forage Evaluation in Ruminant Nutrition; CABI Publishing: Wallingford, UK, 2000; pp. 301–320. ISBN 0851993443. [Google Scholar]
  12. Alomar, D.; Fuchslocher, R.; De Pablo, M. Effect of preparation method on composition and NIR spectra of forage samples. Anim. Feed Sci. Technol. 2003, 107, 191–200. [Google Scholar] [CrossRef]
  13. Reeves, J.B., III. Use of Near Infrared Reflectance Spectroscopy. In Farm Animal Metabolism and Nutrition; CABI Publishing: Wallingford, UK, 2009; pp. 184–209. ISBN 0851993788. [Google Scholar]
  14. Bricklemyer, R.S.; Brown, D.J. On-the-go VisNIR: Potential and limitations for mapping soil clay and organic carbon. Comput. Electron. Agric. 2010, 70, 209–216. [Google Scholar] [CrossRef]
  15. Ely, D.R.; Thommes, M.; Carvajal, M.T. Analysis of the effects of particle size and densification on NIR spectra. Colloids Surfaces A Physicochem. Eng. Asp. 2008, 331, 63–67. [Google Scholar] [CrossRef]
  16. Duyens, L.N.M. The flattering of the absorption spectrum of suspensions, as compared to that of solutions. Biochim. Biophys. Acta 1956, 19, 1–12. [Google Scholar] [CrossRef]
  17. Merzlyak, M.N.; Chivkunova, O.B.; Zhigalova, T.V.; Naqvi, K.R. Light absorption by isolated chloroplasts and leaves: Effects of scattering and “packing”. Photosynth. Res. 2009, 102, 31–41. [Google Scholar] [CrossRef]
  18. Barthès, B.G.; Brunet, D.; Ferrer, H.; Chotte, J.L.; Feller, C. Determination of total carbon and nitrogen content in a range of tropical soils using near infrared spectroscopy: Influence of replication and sample grinding and drying. J. Near Infrared Spectrosc. 2006, 14, 341–348. [Google Scholar] [CrossRef]
  19. De La Haba, M.J.; Garrido-Varo, A.; Núñez, N.; Pérez-Marín, D.C.; Guerrero, J.E. Using calibrations developed for fine milled meat and bone meal on spectra measured on non-milled samples. J. Near Infrared Spectrosc. 2008, 16, 275–279. [Google Scholar] [CrossRef]
  20. Feudale, R.N.; Woody, N.A.; Tan, H.; Myles, A.J.; Brown, S.D.; Ferré, J. Transfer of multivariate calibration models: A review. Chemom. Intell. Lab. Syst. 2002, 64, 181–192. [Google Scholar] [CrossRef]
  21. Shenk, J.S.; Westerhaus, M.O. Optical Instrument Calibration System. U.S. Patent No. 4,866,644, 12 September 1989. [Google Scholar]
  22. Bouveresse, E.; Massart, D.L. Standardisation of near-infrared spectrometric instruments: A review. Vib. Spectrosc. 1996, 11, 3–15. [Google Scholar] [CrossRef]
  23. Salguero-Chaparro, L.; Palagos, B.; Peña-Rodríguez, F.; Roger, J.M. Calibration transfer of intact olive NIR spectra between a pre-dispersive instrument and a portable spectrometer. Comput. Electron. Agric. 2013, 96, 202–208. [Google Scholar] [CrossRef] [Green Version]
  24. Dong, X.; Dong, J.; Li, Y.; Xu, H.; Tang, X. Maintaining the predictive abilities of egg freshness models on new variety based on VIS-NIR spectroscopy technique. Comput. Electron. Agric. 2019, 156, 669–676. [Google Scholar] [CrossRef]
  25. Abdelkader, M.F.; Cooper, J.B.; Larkin, C.M. Calibration transfer of partial least squares jet fuel property models using a segmented virtual standards slope-bias correction method. Chemom. Intell. Lab. Syst. 2012, 110, 64–73. [Google Scholar] [CrossRef]
  26. Shan, P.; Zhao, Y.; Wang, Q.; Wang, S.; Ying, Y.; Peng, S. A nonlinear calibration transfer method based on joint kernel subspace. Chemom. Intell. Lab. Syst. 2021, 210, 104247. [Google Scholar] [CrossRef]
  27. Schaller, G.B.; Hu, J.; Pan, W.; Zhu, J. The Giant Pandas of Wolong; The University of Chicago Press: Chicago, IL, USA, 1985. [Google Scholar]
  28. Johnson, K.G.; Schaller, G.B.; Jinchu, H. Comparative Behavior of Red and Giant Pandas in the Wolong Reserve, China. J. Mammal. 1988, 69, 552–564. [Google Scholar] [CrossRef]
  29. Wei, F.; Wang, W.; Zhou, A.; Hu, J.; Wei, Y. Preliminary study on food selection and feeding strategy of red pandas. Acta Theriol. Sin. 1995, 15, 259–266. [Google Scholar]
  30. Zhang, J.; Hull, V.; Huang, J.; Yang, W.; Zhou, S.; Xu, W.; Huang, Y.; Ouyang, Z.; Zhang, H.; Liu, J. Natural recovery and restoration in giant panda habitat after the Wenchuan earthquake. For. Ecol. Manag. 2014, 319, 1–9. [Google Scholar] [CrossRef]
  31. Roberts, B.M.S.; Gittleman, J.L. Ailurus fulgens. Mamm. Species 1984, Issue 222. 1–8. [Google Scholar] [CrossRef]
  32. Wei, F.; Feng, Z.; Wang, Z.; Zhou, A.; Hu, J. Use of the nutrients in bamboo by the red panda (Ailurus fulgens). J. Zool. 1999, 248, 535–541. [Google Scholar] [CrossRef]
  33. Dierenfeld, E.S.; Hintz, H.F.; Robertson, J.B.; Van Soest, P.J.; Oftedal, O.T. Utilization of bamboo by the giant panda. J. Nutr. 1982, 112, 636–641. [Google Scholar] [CrossRef]
  34. Finley, T.G.; Sikes, R.S.; Parsons, J.L.; Rude, B.J.; Bissell, H.A.; Ouellette, J.R. Energy digestibility of giant pandas on bamboo-only and on supplemented diets. Zoo Biol. 2011, 30, 121–133. [Google Scholar] [CrossRef]
  35. Wei, F.; Wang, Z.; Feng, Z.; Li, M.; Zhou, A. Seasonal energy utilization in bamboo by the red panda (Ailurus fulgens). Zoo Biol. Publ. Affil. Am. Zoo Aquar. Assoc. 2000, 19, 27–33. [Google Scholar] [CrossRef]
  36. Knott, K.K.; Christian, A.L.; Falcone, J.F.; Vance, C.K.; Bauer, L.L.; Fahey, G.C.; Kouba, A.J. Phenological changes in bamboo carbohydrates explain the preference for culm over leaves by giant pandas (Ailuropoda melanoleuca) during spring. PLoS ONE 2017, 12, e0177582. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Nie, Y.; Zhang, Z.; Raubenheimer, D.; Elser, J.J.; Wei, W.; Wei, F. Obligate herbivory in an ancestrally carnivorous lineage: The giant panda and bamboo from the perspective of nutritional geometry. Funct. Ecol. 2015, 29, 26–34. [Google Scholar] [CrossRef]
  38. Dairy, N. Kjeldahl Nitrogen Analysis as a Reference Method for Protein Determination in Dairy Products. J. AOAC Int. 1999, 82, 1389–1398. [Google Scholar]
  39. Bernhard, T.; Truberg, B.; Friedt, W.; Snowdon, R.; Wittkop, B. Development of Near-Infrared Reflection Spectroscopy Calibrations for Crude Protein and Dry Matter Content in Fresh and Dried Potato Tuber Samples. Potato Res. 2016, 59, 149–165. [Google Scholar] [CrossRef]
  40. Williams, P.C. Implementation of near-infrared technology. In Near-Infrared Technology in the Agricultural and Food Industries; Williams, P., Norris, K., Eds.; American Association of Cereal Chemists Inc.: Saint Paul, MN, USA, 1987; pp. 145–169. [Google Scholar]
  41. Kuang, B.; Mouazen, A.M. Non-biased prediction of soil organic carbon and total nitrogen with vis-NIR spectroscopy, as affected by soil moisture content and texture. Biosyst. Eng. 2013, 114, 249–258. [Google Scholar] [CrossRef] [Green Version]
  42. Büning-Pfaue, H. Analysis of water in food by near infrared spectroscopy. Food Chem. 2003, 82, 107–115. [Google Scholar] [CrossRef]
  43. Manley, M. Near-infrared spectroscopy and hyperspectral imaging: Non-destructive analysis of biological materials. Chem. Soc. Rev. 2014, 43, 8200–8214. [Google Scholar] [CrossRef] [Green Version]
  44. Bai, S.; Nayar, R.; Carpenter, J.F.; Manning, M.C. Noninvasive determination of protein conformation in the solid state using near infrared (NIR) spectroscopy. J. Pharm. Sci. 2005, 94, 2030–2038. [Google Scholar] [CrossRef]
  45. Shenk, J.S.; Workman, J.J., Jr.; Westerhaus, M.O. Application of NIR Spectroscopy to Agricultural Products. In Handbook of Near-Infrared Analysis; Burns, D.A., Ed.; CRC Press: Boca Raton, FL, USA, 2008; pp. 347–386. [Google Scholar]
  46. Chen, J.; Ren, X.; Zhang, Q.; Diao, X.; Shen, Q. Determination of protein, total carbohydrates and crude fat contents of foxtail millet using effective wavelengths in NIR spectroscopy. J. Cereal Sci. 2013, 58, 241–247. [Google Scholar] [CrossRef]
  47. Lu, Y.; Du, C.; Yu, C.; Zhou, J. Fast and nondestructive determination of protein content in rapeseeds (Brassica napus L.) using Fourier transform infrared photoacoustic spectroscopy (FTIR-PAS). J. Sci. Food Agric. 2014, 94, 2239–2245. [Google Scholar] [CrossRef]
  48. Halvorson, J.J.; Cassida, K.A.; Turner, K.E.; Belesky, D.P. Nutritive value of bamboo as browse for livestock. Renew. Agric. Food Syst. 2010, 26, 161–170. [Google Scholar] [CrossRef]
  49. Hodkinson, T.R.; Renvoize, S.A.; Chonghaile, G.N.; Stapleton, C.M.A.; Chase, M.W. A comparison of ITS nuclear rDNA sequence data and AFLP markers for phylogenetic studies in Phyllostachys (Bambusoideae, Poaceae). J. Plant Res. 2000, 113, 259–269. [Google Scholar] [CrossRef]
  50. Zhang, Y.; Zeng, C.; Li, D. Scanning electron microscopy of the leaf epidermis in Arundinarieae (Poaceae: Bambusoideae): Evolutionary implications of selected micromorphological features. Bot. J. Linn. 2014, 176, 46–65. [Google Scholar] [CrossRef] [Green Version]
  51. Minasny, B.; McBratney, A.B.; Bellon-Maurel, V.; Roger, J.M.; Gobrecht, A.; Ferrand, L.; Joalland, S. Removing the effect of soil moisture from NIR diffuse reflectance spectra for the prediction of soil organic carbon. Geoderma 2011, 167–168, 118–124. [Google Scholar] [CrossRef] [Green Version]
  52. Wang, A.; Yang, P.; Chen, J.; Wu, Z.; Jia, Y.; Ma, C.; Zhan, X. A new calibration model transferring strategy maintaining the predictive abilities of NIR multivariate calibration model applied in different batches process of extraction. Infrared Phys. Technol. 2019, 103, 103046. [Google Scholar] [CrossRef]
  53. Hu, Y.; Wu, Q.; Ma, S.; Ma, T.; Shan, L.; Wang, X.; Nie, Y.; Ning, Z.; Yan, L.; Xiu, Y.; et al. Comparative genomics reveals convergent evolution between the bamboo-eating giant and red pandas. Proc. Natl. Acad. Sci. USA 2017, 114, 1081–1086. [Google Scholar] [CrossRef] [Green Version]
  54. Nie, Y.; Wei, F.; Zhou, W.; Hu, Y.; Senior, A.M.; Wu, Q.; Yan, L.; Raubenheimer, D. Giant Pandas are Macronutritional Carnivores. Curr. Biol. 2019, 29, 1677–1682.e2. [Google Scholar] [CrossRef]
  55. Zhang, Z.; Hu, J.; Yang, J.; Li, M.; Wei, F. Food habits and space-use of red pandas Ailurus fulgens in the Fengtongzhai Nature Reserve, China: Food effects and behavioural responses. Acta Theriol. 2009, 54, 225–234. [Google Scholar] [CrossRef]
  56. Morgan, L.R.; Marsh, K.J.; Tolleson, D.R.; Youngentob, K.N. The application of nirs to determine animal physiological traits for wildlife management and conservation. Remote Sens. 2021, 13, 3699. [Google Scholar] [CrossRef]
Figure 1. Four processing modes of forage samples: (A) Fresh–Whole, (B) Dry–Whole, (C) Coarse–Ground, and (D) Fine–Ground.
Figure 1. Four processing modes of forage samples: (A) Fresh–Whole, (B) Dry–Whole, (C) Coarse–Ground, and (D) Fine–Ground.
Remotesensing 14 01302 g001
Figure 2. Flow chart of the sample processing for NIR data collection from field samples to the most processed for obtaining the best NIR spectra. The sample processing steps eliminate complicating factors in the NIR spectrum, such as overlapping water signals, mixture inhomogeneity, and scattering from variable particle size and shape.
Figure 2. Flow chart of the sample processing for NIR data collection from field samples to the most processed for obtaining the best NIR spectra. The sample processing steps eliminate complicating factors in the NIR spectrum, such as overlapping water signals, mixture inhomogeneity, and scattering from variable particle size and shape.
Remotesensing 14 01302 g002
Figure 3. Distribution of crude protein of bamboo samples grouped by species presented by the histogram (A) and boxplot (B).
Figure 3. Distribution of crude protein of bamboo samples grouped by species presented by the histogram (A) and boxplot (B).
Remotesensing 14 01302 g003
Figure 4. Raw (up) and transformed (down) spectra from four processing modes, grouped by species. Spectra collected by each processing mode from each species were averaged. Fresh–Whole (FW), Dry–Whole (DW), Coarse–Ground (CG), Fine–Ground (FG).
Figure 4. Raw (up) and transformed (down) spectra from four processing modes, grouped by species. Spectra collected by each processing mode from each species were averaged. Fresh–Whole (FW), Dry–Whole (DW), Coarse–Ground (CG), Fine–Ground (FG).
Remotesensing 14 01302 g004
Figure 5. PCA scores plot (A) and loadings plot (B) for spectra collected from Fine–Ground bamboo samples. PCA scores plot (C) and loadings plot (D) for spectra collected from Coarse–Ground bamboo samples.
Figure 5. PCA scores plot (A) and loadings plot (B) for spectra collected from Fine–Ground bamboo samples. PCA scores plot (C) and loadings plot (D) for spectra collected from Coarse–Ground bamboo samples.
Remotesensing 14 01302 g005
Figure 6. PCA scores plot (A) and loadings plot (B) for spectra collected from Dry–Whole bamboo samples. PCA scores plot (C) and loadings plot (D) for spectra collected from Fresh–Whole bamboo samples.
Figure 6. PCA scores plot (A) and loadings plot (B) for spectra collected from Dry–Whole bamboo samples. PCA scores plot (C) and loadings plot (D) for spectra collected from Fresh–Whole bamboo samples.
Remotesensing 14 01302 g006
Figure 7. The root mean square error of prediction (RMSEP), standard error of prediction (SEP), and the residual prediction deviation (RPD) values before and after slope and bias correction (SBC), predicting samples in satellite modes using Fine–Ground (FG), Coarse–Ground (CG), and Dry–Whole (DW) models. RMSEP and SEP values decrease after SBC, and RPD values increase after SBC.
Figure 7. The root mean square error of prediction (RMSEP), standard error of prediction (SEP), and the residual prediction deviation (RPD) values before and after slope and bias correction (SBC), predicting samples in satellite modes using Fine–Ground (FG), Coarse–Ground (CG), and Dry–Whole (DW) models. RMSEP and SEP values decrease after SBC, and RPD values increase after SBC.
Remotesensing 14 01302 g007
Figure 8. Predicted vs. Reference values of CP. Fine–Ground samples using Fine–Ground model (A), Coarse–Ground samples using Fine–Ground and Coarse–Ground models (B), Dry–Whole samples using Fine–Ground and Coarse–Ground models (C), and Fresh–Whole samples using Fine–Ground, Coarse–Ground, and Dry–Whole models (D). SBC was used in each mode cloning.
Figure 8. Predicted vs. Reference values of CP. Fine–Ground samples using Fine–Ground model (A), Coarse–Ground samples using Fine–Ground and Coarse–Ground models (B), Dry–Whole samples using Fine–Ground and Coarse–Ground models (C), and Fresh–Whole samples using Fine–Ground, Coarse–Ground, and Dry–Whole models (D). SBC was used in each mode cloning.
Remotesensing 14 01302 g008
Table 1. Balanced databases containing NIR spectra collected from each bamboo species under each processing mode. Databases 2–4 are generated using spectra from two bamboo species and the calibration models were tested on spectra collected from the third bamboo species.
Table 1. Balanced databases containing NIR spectra collected from each bamboo species under each processing mode. Databases 2–4 are generated using spectra from two bamboo species and the calibration models were tested on spectra collected from the third bamboo species.
DatabaseSpeciesCalibration for
Each Processing Mode
Validation for Each
Processing Mode
External Validation for
Each Processing Mode
1Rubro2360
Bissetii2360
Sulcata2360
2Rubro2360
Sulcata2360
Bissetii0047
3Rubro3180
Bissetii3180
Sulcata0029
4Bissetii2360
Sulcata2360
Rubro0039
Table 2. Prediction results of the PLSR models for predicting the quantity of crude protein across bamboo species. The calibration model of each processing mode was applied to its corresponding sample type in all databases. In Database 1, spectra from all bamboo species were combined to generate calibration models. Calibration models for crude protein were built from Sulcata and Rubro and used to predict Bissetii in Database 2; models were built from Bissetii and Rubro and used to predict crude protein in Sulcata in Database 3, and models built from Bissetii and Sulcata were used to predict crude protein in Rubro in Database 4.
Table 2. Prediction results of the PLSR models for predicting the quantity of crude protein across bamboo species. The calibration model of each processing mode was applied to its corresponding sample type in all databases. In Database 1, spectra from all bamboo species were combined to generate calibration models. Calibration models for crude protein were built from Sulcata and Rubro and used to predict Bissetii in Database 2; models were built from Bissetii and Rubro and used to predict crude protein in Sulcata in Database 3, and models built from Bissetii and Sulcata were used to predict crude protein in Rubro in Database 4.
Predicted SpeciesCalibration ModelPCsCalibrationCross-ValidationValidationExternal-Validation
R2 RMSECSECRPDR2 RMSECVSECVRPDR2 RMSEPSEPRPDR2 RMSEPSEPRPD
Database 1
—3 species combined
FG model60.970.340.346.290.960.430.435.000.880.690.702.89
CG model90.970.380.385.700.940.570.573.790.820.880.892.26
DW model90.900.630.643.250.800.940.942.200.421.481.521.34
FW model90.960.360.365.280.830.750.762.500.221.781.841.11
Database 2
—Bissetti
FG model70.980.250.258.150.960.400.415.090.920.590.613.600.910.650.653.37
CG model90.980.210.219.100.960.470.473.980.810.900.932.360.910.820.702.65
DW model60.850.770.782.580.711.121.131.780.751.061.112.000.681.261.161.73
FW model90.960.340.345.030.890.650.662.610.211.962.051.090.391.931.821.13
Database 3
—Sulcata
FG model70.960.290.294.920.930.390.393.640.930.470.463.580.960.700.703.59
CG mode90.970.270.275.860.920.450.463.500.950.360.384.590.761.341.371.86
DW model70.920.420.423.530.720.790.791.870.251.611.651.040.611.531.561.63
FW model90.970.210.216.190.880.470.482.780.371.471.401.130.292.302.331.09
Database 4
—Rubro
FG model60.970.400.415.750.930.620.633.700.990.290.2612.800.820.650.652.14
CG model60.970.410.415.860.940.590.604.060.761.881.911.950.841.200.561.16
DW model60.910.570.583.330.611.261.291.510.442.752.831.330.301.431.450.97
FW model110.960.460.465.050.581.441.461.610.802.012.101.820.341.441.380.97
The root mean square error of calibration (RMSEC), standard error of calibration (SEC), the coefficient of determination (R2), and the residual prediction deviation (RPD). The root mean square error of cross-validation (RMSECV), standard error of cross-validation (SECV), the coefficient of determination (R2), and the residual prediction deviation (RPD). The root mean square error of prediction (RMSEP), standard error of prediction (SEP), the coefficient of determination (R2), and the residual prediction deviation (RPD). Fresh–Whole (FW), Dry–Whole (DW), Coarse–Ground (CG), Fine–Ground (FG).
Table 3. Prediction results of the PLSR models for crude protein transferred across processing modes before and after the slope and bias correction (SBC) were applied to Database 2. The calibration model from the Master mode (Fine–Ground) was applied to each of the three satellite modes. The calibration model from the Coarse–Ground mode was applied to the Dry–Whole and Fresh–Whole modes, and the calibration from the Dry–Whole mode was applied to the Fresh–Whole mode.
Table 3. Prediction results of the PLSR models for crude protein transferred across processing modes before and after the slope and bias correction (SBC) were applied to Database 2. The calibration model from the Master mode (Fine–Ground) was applied to each of the three satellite modes. The calibration model from the Coarse–Ground mode was applied to the Dry–Whole and Fresh–Whole modes, and the calibration from the Dry–Whole mode was applied to the Fresh–Whole mode.
Calibration ModelSlopeOffsetCorrelationR2 RMSEPSEPBiasRPD
FG sample–FG model0.852.590.960.910.650.650.023.37
CG sample−FG model—before SBC0.960.510.930.870.830.82−0.232.62
CG sample–FG model—after SBC0.823.190.930.870.790.800.162.73
CG sample–CG model0.803.930.950.910.820.700.442.65
DW sample–FG model—before SBC0.608.430.850.731.991.191.621.10
DW sample–FG model—after SBC0.587.260.850.731.191.200.201.84
DW sample–CG model—before SBC0.4112.510.720.512.671.532.210.80
DW sample–CG model—after SBC0.3112.400.720.511.621.620.381.32
DW sample–DW model0.577.270.830.681.261.26−0.161.73
FW sample–FG model—before SBC0.70−9.300.660.4414.421.77−14.320.15
FW sample–FG model—after SBC0.2612.960.660.441.721.670.561.22
FW sample–CG model—before SBC0.68−13.330.530.2819.052.45−18.910.11
FW sample–CG model—after SBC0.1914.410.530.281.851.860.361.16
FW sample–DW model—before SBC0.50−6.150.770.5914.591.42−14.520.15
FW sample–DW model—after SBC0.4010.220.770.591.461.500.091.49
FW sample–FW model0.576.680.620.391.931.82−0.681.13
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sheng, Q.; Santos-Rivera, M.; Ouyang, X.; Kouba, A.J.; Vance, C.K. Near-Infrared Spectroscopy and Mode Cloning (NIR-MC) for In-Situ Analysis of Crude Protein in Bamboo. Remote Sens. 2022, 14, 1302. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14061302

AMA Style

Sheng Q, Santos-Rivera M, Ouyang X, Kouba AJ, Vance CK. Near-Infrared Spectroscopy and Mode Cloning (NIR-MC) for In-Situ Analysis of Crude Protein in Bamboo. Remote Sensing. 2022; 14(6):1302. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14061302

Chicago/Turabian Style

Sheng, Qingyu, Mariana Santos-Rivera, Xiaoguang Ouyang, Andrew J. Kouba, and Carrie K. Vance. 2022. "Near-Infrared Spectroscopy and Mode Cloning (NIR-MC) for In-Situ Analysis of Crude Protein in Bamboo" Remote Sensing 14, no. 6: 1302. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14061302

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