Next Article in Journal
Development of a Guidance System for Motor Imagery Enhancement Using the Virtual Hand Illusion
Next Article in Special Issue
Detecting Bacterial Biofilms Using Fluorescence Hyperspectral Imaging and Various Discriminant Analyses
Previous Article in Journal
Towards Precise Positioning and Movement of UAVs for Near-Wall Tasks in GNSS-Denied Environments
Previous Article in Special Issue
The Use of a Micro Near Infrared Portable Instrument to Predict Bioactive Compounds in a Wild Harvested Fruit—Kakadu Plum (Terminalia ferdinandiana)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Minimum Preparation Sampling Strategies for Sugarcane Quality Prediction by vis-NIR Spectroscopy

by
Lucas de Paula Corrêdo
*,
Leonardo Felipe Maldaner
,
Helizani Couto Bazame
and
José Paulo Molin
Precision Agriculture Laboratory, Biosystems Engineering Department, ‘Luiz de Queiroz’ College of Agriculture, University of São Paulo, Av. Pádua Dias 11, 13418-900 Piracicaba, Brazil
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(6), 2195; https://doi.org/10.3390/s21062195
Submission received: 23 January 2021 / Revised: 10 March 2021 / Accepted: 19 March 2021 / Published: 21 March 2021
(This article belongs to the Special Issue Using Vis-NIR Spectroscopy for Predicting Quality Compounds in Foods)

Abstract

:
Proximal sensing for assessing sugarcane quality information during harvest can be affected by various factors, including the type of sample preparation. The objective of this study was to determine the best sugarcane sample type and analyze the spectral response for the prediction of quality parameters of sugarcane from visible and near-infrared (vis-NIR) spectroscopy. The sampling and spectral data acquisition were performed during the analysis of samples by conventional methods in a sugar mill laboratory. Samples of billets were collected and four modes of scanning and sample preparation were evaluated: outer-surface (‘skin’) (SS), cross-sectional scanning (CSS), defibrated cane (DF), and raw juice (RJ) to analyze the parameters soluble solids content (Brix), saccharose (Pol), fibre, pol of cane and total recoverable sugars (TRS). Predictive models based on Partial Least Square Regression (PLSR) were built with the vis-NIR spectral measurements. There was no significant difference (p-value > 0.05) between the accuracy SS and CSS samples compared to DF and RJ samples for all prediction models. However, DF samples presented the best predictive performance values for the main sugarcane quality parameters, and required only minimal sample preparation. The results contribute to advancing the development of on-board quality monitoring in sugarcane, indicating better sampling strategies.

1. Introduction

Near-Infrared (NIR) spectroscopy is a well-established technique to monitor the quality of raw sugarcane received by sugar mills [1], and consequently, for pricing and trading with producers and growers [2]. Crop quality is estimated based on physicochemical parameters related to physiological composition, such as soluble solids content (Brix), water-insoluble solids (Fibre), and the apparent sucrose in the juice (Pol). Furthermore, all other parameters (purity, Pol of cane, reducing sugars, and total recoverable sugars) are calculated based on the former parameters [3], from which total recoverable sugars (TRS) are used for the pricing and trading of the raw material. Sugarcane quality parameters are determined by analytical methods and empirical equations described by the National Council of Sugarcane Producers (CONSECANA), which are based on the International Commission for Uniform Methods of Sugar Analysis (ICUMSA). Using calibration methods, it is now possible to obtain some crop quality parameters by NIR spectroscopy [4].
Different wavelength regions of the electromagnetic spectrum can be used in spectroscopy, such as visible (400 to 750 nm), near-infrared (NIR, 750 to 2500 nm), shortwave near-infrared (SWNIR, 750 to 1100 nm), and visible and near-infrared (vis-NIR, 400 to 2500 nm) [5,6,7,8]. However, there is no consensus on the limits between these regions. Interaction between electromagnetic radiation and matter causes molecular vibrations involving heavy atom (C, N, O, and S) attached to a hydrogen atom [8]. This basic principle has allowed substantial scientific advances to predict organic compounds of agricultural products associated with its quality using vis-NIR spectroscopy as a nondestructive and environmentally friendly analysis technique [5]. Moreover, several studies have shown promising results when using the technique to predict sugar cane quality in the sugar mill [9,10,11,12] and for breeding programs [13,14,15].
Despite the advances in industrial sugarcane quality monitoring, spectral methods are still a distant reality for in-field measurements in line with precision agriculture practices. Current proximal sensing technologies applied at the canopy level allow only monitoring crop yield [16,17]. However, some studies have indicated that vis-NIR could also be a viable technology for acquiring quality data of harvested products in real-time during mechanical harvesting [16]. The monitoring of crop quality parameters across the field is important to adopting precision agriculture (PA) practices, in which quality maps would show the variability of the crop and help guide site-specific management [18]. In this context, mechanical harvest opens a way to obtain a high frequency of sampling and data collection to analyze the quality variability across the field [19]. However, some requirements need to be satisfied to use vis-NIR as proximal sensing technology for this purpose: (i) the location of adaptation in the harvester, (ii) development of a sampler system, and (iii) the type of sampling required for analysis. The first two requirements are fundamentally dependent on the last one.
Nawi et al. (2014) indicated that the ideal place for implementing an on-board sugarcane quality monitor would be in the elevator of the harvester, where the sugarcane material is partially cleaned and processed in the form of billets. In this context, some studies have reported promising results on the prediction of sugarcane Brix from sensor readings made on the outer-surface (‘skin’) [20] or on cross-sections [21] of sugarcane billets. Furthermore, more recent studies have advanced with on-board vis-NIR spectroscopy sensor applications on the elevator of a sugarcane harvester simulator, i.e., analyzing samples at a distance and in motion [2,22]. However, Maraphum et al. [23] and Phuphaphud et al. [24] reported that the waxy material should be removed from the cane surface for maximum accuracy in the spectral data condition, even though this may be impractical for an embedded system. Associated with this fact, Phetpan et al. [2] reported on the need to evaluate the potential of the vis-NIR spectroscopy technique with data sets consisting of a larger number of sugarcane varieties. In addition, despite the advantages of nondestructive measurement, there is a lack of basic studies comparatively evaluating various forms of sampling without and with minimal processing, using extensive numbers of samples obtained over the course of a harvest. Thus, the objective of this study was to compare different sugarcane sample types, including billets, defibrated cane, and raw juice, and to analyze the spectral response of each sampling type for the prediction of quality parameters of sugarcane from vis-NIR spectroscopy.

2. Materials and Methods

2.1. Sampling

Variability of sugarcane quality parameters throughout a harvest occurs due to diverse environmental conditions, mainly temperature and precipitation, during the harvest [25]. Based on this, the data collection procedure occurred on random periods over six months of the 2019 harvest (June to November). We sought to obtain higher variability of the sugarcane quality parameters from this collection procedure throughout the harvest.
The data collection was carried out in the quality laboratory of a sugar mill. Three hundred and two samples were collected, and different levels of processing were applied. Also, all data acquisition was performed in a controlled temperature environment (20 ± 5 °C), minimizing the effects of sugar degradation by microorganisms. In addition, we collected the daily results of conventional analysis performed by sugar mill quality laboratory over the same months in which the samples were collected for spectral analysis.
The sampling procedure for vis-NIR analysis was carried out sequentially and simultaneously to the sample preparation for conventional analysis, as described:
  • An oblique probe collected a sample of sugarcane billets in each truckload to proceed with the conventional analyzes of the sugarcane transported to the sugar mill (Figure 1a). Before the laboratory processes the sample, we took a subsample composed of three random sugarcane billets (Figure 1b);
  • The remaining sample of billets was milled in a mechanical knife crusher and homogenized in a mixer (Figure 1c). Then, samples with 500 g and 50 g of the homogenized defibrated sugarcane were collected (Figure 1d) to proceed with the conventional and spectral analysis, respectively;
  • The 500 g sample was pressed in a hydraulic press under constant pressure, at 24.5 MPa (250 kgf cm−2) for 60 s to obtain the raw juice for conventional analysis (Figure 1e). At the same moment, the third sample composed of extracted raw juice was collected for spectral measurements (Figure 1f).
The sample types composed of three billets, defibrated cane, and raw juice were prepared (Figure 1g,h) and immediately used for spectral measurements at the sugar mill laboratory (Figure 1i). The remains of raw juice and the bagasse without juice (after pressed) were used for conventional analytical analysis (Figure 1j).

2.2. Sugarcane Quality Analysis

All the procedures and equations for sugarcane quality parameters determination followed the standard protocol proposed by the CONSECANA [26]. These protocols agree with the international rules from ICUMSA.
Initially, the soluble solids content (Brix) was determined by pouring raw juice into a refractometer probe (RX-5000α, ATAGO Co Ltd., Tokyo, Japan) with a maximum resolution of 0.1 Brix. Then, 14 g of a mixture composed of a 4:2:1 proportion of Celite (mineral filtering agent), aluminum chloride, and calcium hydroxide, respectively, was added to 200 mL of raw juice homogenized by a magnetic stirrer until the solution was well-mixed. The solution was filtered through filter paper to obtaining clarified juice. A volume of 70 mL was added in a digital polarimeter (Schmidt + Haensch, Polartronic NHZ 8, Berlin, Germany) to the saccharimetric reading. The result was obtained as percentage of apparent sucrose in the juice (Pol).
After the juice extraction, the remaining fibrous cane residue (Figure 1j) was transferred to metal baskets, with holes at the base. The fibrous residue was weighed on a semianalytical balance. Then, samples were maintained in a forced air circulation dryer at a temperature of 105 °C, until constant weight was achieved. The dryer samples were weighed, and the fibre content (insoluble solids) was determined.
The Pol of cane and TRS were calculated for each sample from the previous parameters as described by CONSECANA (2015). The triplicate values obtained from each sample were averaged.
It is important to mention that reducing sugars (fructose and glucose) and purity (apparent sucrose in the soluble solids content) were also determined in the laboratory. However, while these parameters are not the subject of the present study, they were used to calculate TRS by the reference method.

2.3. Acquisition of Spectral Data

The spectral measurements were performed with a Veris vis-NIR spectrometer (Veris Technologies Inc., Salina, KS, USA). This equipment was developed for on-the-go soil measurements mounted on a platform, connected to a three-point hitch, and pulled by a tractor [27]. However, it may also be used in bench mode. The equipment consists of a CCD array spectrometer (USB4000, Ocean optics, Largo, FL, USA), measuring wavelengths between 373 and 1011 nm, and an InGaAs photodiode-array spectrometer (C9914GB, Hamamatsu Photonics, Hamamatsu, Japan), with a spectral range between 1170 and 2222 nm. The system presents a resolution of around 5 nm. Spectral measurement acquisitions were performed via a sapphire window in the lower of the shank using a tungsten halogen lamp as an electromagnetic energy source (Figure 1i). Each spectrum recorded by the equipment software (Veris spectrophotometer software V1.79) corresponded to the average of 20 spectral readings. The inside shutter is operated automatically to obtain dark and reference spectra before each analysis. Four external references with different grey levels were used for spectral calibration of the spectrometer before the analysis. The spectral data were stored as absorbance units.
The three billets of each sugarcane sample were cut transversely at both extremities, and their skin was lightly cleaned with paper to remove residues from harvest (Figure 1b). An ad hoc dark chamber was constructed with PVC pipes and foam to accommodate the billets, and its inside was painted matte black. Also, magnets were placed on the cover of the chamber, next to the hole through which the spectral scans were performed, aiming to fix the chamber on the reading shank of the spectrometer. This device standardizes the distance between the sapphire window of the spectrometer and the sugarcane billets and removes interference from external lighting (Figure 1g).
The spectral scan method on billets was adapted from Nawi [20,21] and Phuphaphud [15]. The spectral scans were performed at three equidistant points (around 120°) on the skin of each billet. Furthermore, the cross-sectional scanning of billets was performed in triplicate in each cross-sectional surface of each billet, only changing the position after each reading. Therefore, each sample type measurement, skin scanning (SS), and cross-sectional scanning (CSS) of billets was represented by an average of nine successive scans.
The spectral measurements of defibrated cane and raw juice were performed in the same manner. A recipient available from the equipment itself with a volume of around 3 mL was used. The recipient was filled with sample (defibrated cane or extracted raw juice, Figure 1h), and the spectral measurements were performed in triplicate. Three replicates were performed for each sample. Thus, the average of nine spectral readings of defibrated cane (DF) and nine spectral readings of raw juice (RJ) were recorded.

2.4. Spectral Preprocessing

Data preprocessing steps were performed to remove or minimize the sources of spectral variabilities, such as noise present in the dataset, which was not related to the analytical signal [8,28].
Firstly, the spectra were preprocessed using standard normal variate (SNV) [29] to eliminate the deviations caused by particle size and scattering, which centers each spectrum on its mean and then scales it by its standard deviation. Also, the second derivative based on the Savitzky-Golay algorithm [30] was applied, with a window size of 11 points and second-order polynomial fitting to minimize hurdles such as baseline shifts drifts and to remove high-frequency noise from a spectrum and improve the signal-to-noise ratio [28]. After the preprocessing of the spectral data, Pareto scaling (PS) was applied to variables, which is the most commonly applied scaling method in infrared data [31]. The method centered all variables at their means, and then divided them by the square root of the standard deviation.

2.5. Multivariate Analysis

Firstly, the spectral data of the four sampling conditions were concatenated. Then, the data was divided into calibration (75%, 227 samples) and external validation (25%, 75 samples) data sets, based on the Kennard-Stone method [19]. This procedure allowed to obtain the same samples for calibration and external validation data set for both sample types evaluated. The spectral measurements were used to build predictive models for sugarcane quality parameters based on Partial least square regression (PLS) [20].
The models were calibrated using the venetian blinds cross-validation method with 10 splits. The optimal PLS models were determined based on the lowest number of latent variables (LV), in which the root mean square error of cross-validation (RMSECV) was not significantly higher than the minimum RMSECV [23]. The root mean square error (RMSE) was calculated as follows:
RMSE =   i = 1 n ( y i y i ^ ) 2 n
where n is the number of samples, yi is the reference measurement of sample i, and y i ^ is the estimated result for sample i.
The outliers were evaluated during the calibration step for the reference lab values and spectral data. The presence of outliers in the spectral data was evaluated by the “influence plots” based on high leverage and unmodeled residuals by Hotelling T2 and Q statistics, respectively [8]. Samples with high values in both cases, at 5% of significance level, were considered outliers and removed from the spectral data set. On the other hand, outliers in reference data were evaluated by the root mean square error in calibration (RMSEC) values. Samples that presented errors in prediction greater than ±3 × RMSEC were considered outliers and removed from the data set [32]. The external validation samples were considered unknown samples. In this way, the outliers were evaluated only for the spectral data set. The process was carried out at most three times in the calibration step, as recommended by ASTM E1655-7 [33].
The model accuracy was evaluated based on the RMSE for calibration, cross-validation, and prediction (RMSEC, RMSECV, RMSEP, respectively). Prediction performance was evaluated based on the determination coefficient (R2) for calibration and prediction (R2c and R2p, respectively), and the ratio of performance to the interquartile range (RPIQ), which is calculated by the ratio between the interquartile difference and the RMSEP. Also, a randomization test [34] with 0.05 significance level of probability was performed. The aim was to compare the accuracy of regression models using different sugarcane sample types in the validation set. The hypothesis evaluated were:
  • Null hypothesis (H0): RMSEPsample type 1 = RMSEPsample type 2 (accuracy is similar);
  • Alternative hypothesis (H1): RMSEPsample type 1 ≠ RMSEPsample type 2 (accuracy is not similar).
An advantage of this test is its simplicity and the fact that assumptions about normality or homoscedasticity of the data are not required (distribution-free) [35]. More details about this test, included an algorithm script, can be found in Olivieri [36].
Moreover, the variable importance for the projection (VIP) was calculated to verify the wavelengths with a more significant impact on the external validate models [27] for each sample type. The VIP was calculated as follows:
VIP j = p k = 1 h [ Z ( w jk | | w k | | ) 2 ] . ( k = 1 h Z ) 1
where VIP is the variable importance for projection (dimensionless), j is a specific wavelength (nm), p is the number of wavelengths (dimensionless), h is the number of latent variables (dimensionless), Z is the fraction of variance in the prediction explained by the latent variable (dimensionless), and w is the loading weight (dimensionless).
All models, routines, and data processing were performed in Matlab R2015a (The MathWorks, Natick, MA, USA) and PLS Toolbox 8.9 (R8.9.1; Eigenvector Research, Wenatchee, WA, USA).

3. Results and Discussion

3.1. Overview of Sugarcane Quality Reference Data and vis-NIR Spectral Measurements of Different Sample Types

From the daily results of analyses performed by conventional methods at the mill, it was possible to characterize the variation of the main parameters determined analytically (Brix, Pol, and Fibre), as well as for TRS, throughout the months in which the experiment was performed (Figure 2).
It is possible to observe an increasing trend in all parameters from June to October. Afterward, there is a tendency to decrease, except for fibre. Weather is highly influential on sucrose storage [25]. In the months corresponding to autumn (June) and winter (June to September), water stress and cooler temperatures contribute to the reduction of vegetative crop growth and favor sucrose storage [37]. With the beginning of spring (September/October) and the beginning of the rainy season, the vegetative growth of the crop is resumed, and the reserves are consumed. The sample acquisition on different periods (vertical bars in Figure 2) throughout the harvest allowed us to obtain data including different stages in this variation. The effect of this variability was reflected in the range of all samples collected during the experiment (Table 1).
On a first view, the Kennard-Stone method provides a representative calibration data set, with external validation data set between its range (Table 1). The sample acquisition method provided a satisfactory variability of data, as expected; TRS varied from 86.94 to 173.80 kg of sugar per Mg of cane.
The distribution of all quality parameter values had wide distribution (Figure 3). Fibre content did not positively or negatively correlate with any other parameters analyzed, with values varying from −0.16 to 0.13 (p < 0.05). On the other hand, the other parameters showed a positive correlation higher than 0.94 (p < 0.05). The highest correlation was observed between Pol of cane and TRS, close to 1.00 (p < 0.05). Higher correlation values are observed between the TRS with parameters analytically determined such as Brix and Pol (0.94 and 0.96, respectively, p < 0.05). The correlation values for these attributes are firstly explained by the composition of the soluble solids content of sugarcane, measured by Brix, in which the largest proportion corresponds to sucrose (about 15–18%), measured by Pol [38,39], and reducing sugars (fructose and glucose) in a smaller proportion (about 0.5%) [40]; note that the determination of reducing sugars was not the objective of the present study.
The vis-NIR raw spectral data obtained for 302 samples of each sugarcane sample type are shown in Figure 4. The spectral data were evaluated to identify possible spectral errors [41]. We observed a noisy aspect in the region corresponding to the visible spectrum (400 to 698 nm), mainly for SS samples. This effect may have been attributed to the influence of skin colors of billet samples, which were obtained for several different sugarcane varieties (Table A1), or soil residues from the harvest present in the RJ samples. Therefore, this spectral region was removed from the data set. Phuphaphud [14] observed the same effect due to the skin color of sugarcane billets. Also, based on the evaluation of the coefficient of variation (CV) obtained for each spectral band, the last spectral bands showed high CV concerning their neighbors and were also removed from the dataset, similar to the procedure performed by Franceschini [27] in a study on the external effects on the spectral reading of vis-NIR of soil samples using the same equipment. Thus, only bands in the spectral range between 699 and 1010 nm and between 1070 and 2153 nm (303 spectral bands) were retained.
A PCA analysis performed an exploratory overview of the data structure. The spectral data were only mean-centered, and the classes were identified by sample type. Two principal components, PC1, and PC2, explained 98.6% of the data variance (Figure 5). The first component explained 96.8% of the data variance. The data structure was different for each sample type, as can be seen from ellipses illustrating the majority of samples (Figure 5). However, a first overview allowed us to verify the greater difficulty in explaining the variance of less processed samples, such as samples obtained by spectral readings in the skin (SS) and cross-sectional (CSS) of the billets, than processed samples, such as raw juice samples (RJ).
Vis-NIR spectroscopy may be used in a number of applications, including the classification of sugarcane varieties, with promising results [42]. The same authors showed that the spectral regions between 650 and 750 nm, corresponding to the visible spectrum, was the most suitable for sugarcane discrimination. The principal component analysis for the four sample types individually showed that the scatter plots were not categorized into groups based on sugarcane varieties (Figure A1). These results indicated that the sample set composed for many different varieties did not affect the spectral characteristics between each variety for both sample types. A similar effect was observed by Phuphaphud et al. (2020) [14] when evaluating the classification of three varieties. Therefore, the present study was conducted for all varieties combined.
The vis-NIR technique principle is based on the detection of compounds and molecules through their molecular vibration states [8]. Different varieties naturally have different concentrations of parameters such as sucrose and fibre according to genetics. Furthermore, for all of them, the plant matrix is essentially composed of water (75–82%), insoluble solids content (Fibre, 10–18%), and soluble solids (Brix, 18–25%), which are composed of nonsugars (1–2%), sucrose (14–24%), and reducing sugars (0–1.5%) [43]. However, the prediction of quality parameters related to chemical compounds of interest should be independent of sugarcane varieties.

3.2. Prediction Performance of Models Based on Different Sugarcane Sample Types

Figure 6 presents scatter plots showing reference versus predicted values of sugarcane quality parameters. There was an underestimation of high values and overestimation of lower values for all attributes and sample types evaluated. However, this effect was more intense for the less prepared sampling condition, such as SS and CSS. Also, overall, the residuals showed no trend (Figure A2).
More LVs were necessary to explain the variance of the data for models constructed from SS samples (between 7 and 10) than those obtained to predict the same parameters from other sample types (Table 2). Also, it could be observed that SS and CSS did not show similar accuracy (p-value < 0.05) for Fibre, Pol of cane, and TRS (Table 3). Moreover, the prediction performance results for these parameters by these sample types were worse than the performance results obtained by DF and RJ samples (Table 3). The RPIQ values for TRS were 40% higher on average than those for SS and CSS samples, for example. Furthermore, the model accuracy observed between sampling methods for all the other conditions was statistically nonsignificant (p-value > 0.05).
The model performance for DF and RJ samples was equivalent for practically all parameters evaluated. There was no significant difference between the model’s accuracy (p-value > 0.05) and very close values of R2p and RPIQ. Moreover, from DF samples, it was possible to obtain a satisfactory performance to predict Fibre content; this was not possible for RJ samples. On the other hand, the models performed for SS samples presented a higher number of LV than for a prepared sample. The model performance for predicting parameters related to sucrose (Brix, Pol, Pol of cane, and TRS) was not satisfactory, with R2p and RPIQ below 0.5 and 2.0, respectively, except for Fibre prediction. The prediction results from SS samples for Fibre were close to those obtained for DF samples, as shown by the values of R2p and RPIQ. However, the results were less promising than those obtained by Phuphaphud et al. (2019) [15], which obtained the following results: maxima of 0.81 for R2p and 0.63 for RMSEP. Although Fibre content is an important attribute for sugarcane quality determination, it is not essential for sucrose estimation. Fibre content has no relation with some important attributes, such as Brix and Pol, and only minimally impacts TRS calculation. The prediction of this parameter is important for producing energy cane and breeding programs, as in work developed by Phuphaphud et al. (2019) [14].
Some models developed for CSS samples were similar to those developed for SS, as for Brix and Pol prediction. However, its predictive performance was lower than those obtained for Fibre, Pol of cane, and TRS predicted by SS samples, with worse results for R2p and RPIQ. In a first investigation, Nawi et al. (2013a) obtained values of 0.87 for R2p. The excellent performance of this index can be explained by the method of data acquisition adopted by the authors, with individualized samples according to the stem portion (lower, middle, and upper portion) and only three varieties of cane. Sucrose accumulation occurs in an ascending manner, with more accumulation in the lower portion and less in the internodes of the upper portion, close to the leaves [38,44,45]. Therefore, samples composed of different sections resulted in more variability in quality parameters. However, if we analyze the characteristics of a sugar cane harvester, after the stems pass through the chopper roll system, the distinction between portions of the sugarcane stem is not viable.
The RPIQ values for the SS method were higher than those obtained by the CSS method. The SS method on billet samples on the harvester conveyor would be the most practical method, due to the better operability of sample acquisition in that portion of the harvester [16]. However, several external factors must be considered to measure quality attributes by the SS method. A critical one is the constitution of the sugarcane skin itself, as various colors depending on the variety, black and white waxy material, and organic compounds may be present [24,46].
There are common waxy materials on the cane surface that affect vis-NIR measurements by the SS method. Maraphum et al. (2018) evaluated the effect of the waxy material on the cane surface to eliminate or avoid getting low accuracy of the models for Pol measurements. They obtained RMSEP values were around 1.20 to 1.50%, i.e., close to those found by the present study. The authors concluded that spectra acquisition by removed-wax samples was convenient for the measurement of Pol. However, other compounds could affect vis-NIR spectroscopy measurements, such as cellulose and lignin [7].
Overall, the predictive performance results of the models (based on the R2p and RPIQ) indicate that DF and RJ samples presented similar performance and provide the best results. However, models built from DF samples require less preparation, i.e., by avoiding juice extraction, making them more attractive for an on-board system. Moreover, CSS samples presented worse performance than all other samples. SS samples presented higher values of R2p and RPIQ than CSS samples for all quality parameters. Furthermore, SS samples showed a nonsignificant difference (p-value > 0.05) of accuracy (RMSEP) with models built from DF samples, but worse results for performance (R2p and RPIQ), except for Fibre. Possibly DF results were satisfactory due to the exposure of the internal constituents, which overlapped concerning the waxy material that becomes visually negligible. On the other hand, the organic compounds in the sugarcane outer-surface may have interfered in the prediction models. Future studies using nonlinear processing methods [47] or advanced filtering methods, with orthogonalization of unwanted signals concerning the compounds of interest [48], may help in increasing the predictive performance of the models, which is more interesting for an on-board system.

3.3. Variable Influence on the Models

VIP scores were used to describe the importance of each wavelength to the prediction of the main sugarcane quality parameters (Figure 7).
As a vibrational spectroscopy technique, the interaction between the vis-NIR electromagnetic radiation and the matter of the sample could be interpreted mainly by overtones and combinations of vibrational modes involving C-H, O-H, and N-H chemical bonds [8,49]. VIP values greater than 1.0 indicate variables with greater influence on the models, and VIP values between 0.8 and 1.0 indicate the moderately influential variables. All variables with VIP smaller than 0.8 are insignificant to the predictions [50]. Some substantial similarities could be observed for different samples and quality parameters on a first overview.
At around 960 nm, there is a small interval with high values of VIP (higher than 2.3), especially related to the prediction of Brix, Pol, Pol of cane, and TRS by CSS and RJ samples, corresponding to the second and third overtone of O-H and C-H stretching, respectively [51]. The region between 980 and 1030 nm can be regarded as an important contributor to quality predictions (VIP higher than 1.0). Between 972 and 1009 nm there is a characteristic signal related to saccharides and the third overtone of O-H [7]. This signal is more expressive for parameters determined for SS samples. This spectral range could be associated with cellulosic fibres, which explained the higher VIP values to SS samples. Similar observations were found by Phuphaphud et al. (2020) to predict commercial cane sugar from growing cane stalks for breeding programs using vis-NIR spectroscopy.
At 1139 nm, there is a small band with VIP values higher than 1.0 for all parameters predicted for four sample types, except for Fibre content. On the interval between 1100 and 1230 nm occurs the second vibrational frequency overtones associated with C-H stretching [51]. Also, at around 1170–1180 nm, there are VIP values higher than 1 for Fibre predicting, mainly by DF and CSS samples. In this region, the third overtone of C-H and unsaturated C=C double bonds are typically associated with fibre, such as lignin [7].
At 1360 nm, there is another expressive region with high VIP values, possibly related to C-H combinations and the O-H first overtone, respectively [49]. Then, at 1600 nm, there are highly similar VIP values possibly associated with to first vibrational frequency overtone of C-H stretching [51]. Another region shows a high contribution, with VIP values higher than 1, for Fibre predicting between 1850 and 1900 nm, mainly for CSS and DF samples. Around 1820 nm occurs the effects of O-H stretching associated with two combinations of C-O stretching commonly associated with Fibre as cellulosic [7]. This interval is lower and with lower VIP values for SS samples, possibly due to the waxy effect on the near-infrared signal [23]. Finally, in the last bands of the spectra, after 2100 nm, the intensity of VIP values is similar for all predicted attributes from any sample type due to O-H bending and C-O stretching combination [7].
The scores of the models for all measured quality parameters are displayed by their first PLS loadings (Figure A3, Figure A4, Figure A5, Figure A6 and Figure A7 see Appendix A), accounting for more than 95% of the data variance. Overall, the most considerable variation occurred in the spectral region between 1300 and 1500 nm and between 1800 and 1950 nm. Other authors have found similar response in these spectral regions for prediction of sugars in other products [52,53]. This effect was similar for all sample types and all parameters evaluated. Therefore, this fact proved the relationships identified by the VIP scores and the key molecular bonds related to the parameters of interest described earlier.
The gap between two spectrometers, starting at 1011 nm until 1070 nm, is not related to significative known vibrational frequency overtones associated with some bands related to sugars or fibres [7]. Therefore, the absence of information in this range would not have significantly affected the development of the models.
Processed samples allowed a more significant interaction of electromagnetic radiation corresponding to vis-NIR bands with matter constituents. This physical effect resulted in more prominent signals from specific vibrational frequency bands related to the chemical constitution of sugarcane quality parameters. Overall, defibrated samples (DF) showed performance prediction results that were close to raw juice samples (RJ). Also, the DF sample allowed us to predict Fibre content as well as other parameters, which is not possible with RJ samples. The prediction of sugarcane quality parameters from less processed samples is a desirable characteristic for mechanization of on-the-go measurements of crops, thus promoting spatial information of crops based on quality. DF samples may partially satisfy this requirement; however, this is a destructive sampling technique.
Improving the predicting performance of sugarcane quality parameters from billets for on-the-go systems may be possible [20]. Some effects, such as waxy and skin organic compounds, need to be considered and minimized [24]. Other data processing techniques such as nonlinear models [47] or advanced filtering methods such as orthogonalization [48], could be investigated to improve the performance aiming to develop reliable models for measuring sugarcane quality using billets of cane.

4. Conclusions

This study demonstrates that vis-NIR spectroscopy could be used as a quick method to assess the abundance of chemical compounds of sugarcane related to its quality. There was no significant difference (p-value > 0.05) in the accuracy (RMSEP) of prediction of whole cane samples when compared to processed samples, such as defibrated cane (DF) and extracted raw juice (RJ), for all evaluated quality parameters. Also, outer-surface measurements of sugarcane billets presented a better accuracy (RMSEP, p-value > 0.05) and performance (R2p and RPIQ) than measurements on the cross-section.
Despite the similar accuracy (p-value > 0.05), DF and RJ sampling presented better performance than outer-surface measurements of sugarcane billets. Moreover, the performance of the models from DF and RJ samples were similar, but DF samples involve less preparation, as they do not require juice extraction of the sample.
The results showed that DF sampling could be used to predict the main sugarcane quality parameters, such as soluble solids content (Brix), saccharose (Pol), Fibre, Pol of cane, and total recoverable sugars (TRS), all of which are used for pricing and trading between mills and sugarcane producers. The DF models presented RMSEP varying between 0.72% and 0.87% for Brix, Pol, Fibre, and Pol of cane, and 6.71 kg Mg−1 for TRS.
The results in this study contribute to advancing the development of on-board quality monitoring in sugarcane. This information shows the spatial variability of crop quality and helps guide site-specific management of sugarcane fields.

Author Contributions

Conceptualization, L.d.P.C., and J.P.M.; methodology, L.d.P.C., and J.P.M.; software, L.d.P.C.; validation, L.d.P.C.; formal analysis, L.d.P.C.; investigation, L.d.P.C., and J.P.M.; resources, L.d.P.C. and J.P.M.; data curation, L.d.P.C.; writing—original draft preparation, L.d.P.C.; writing—review and editing, L.d.P.C., L.F.M., H.C.B., and J.P.M.; visualization, L.d.P.C., L.F.M., H.C.B., and J.P.M.; supervision, J.P.M.; project administration, L.d.P.C. and J.P.M.; funding acquisition, L.d.P.C., and J.P.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by São Paulo Research Foundation (FAPESP), which provided a doctoral scholarship to the first author (grant number 2018/25008-8).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

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

Acknowledgments

The authors would like to thank the São Martinho group for allowed to conduct this study in one of its mills. Also, we would like to thank the researcher André Freitas Colaço, from Commonwealth Scientific and Industrial Research Organization (CSIRO), for his valuable and pertinent contributions on an earlier version of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Number of samples of each Brazilian sugarcane variety used in the study.
Table A1. Number of samples of each Brazilian sugarcane variety used in the study.
Variety# SamplesVariety# SamplesVariety# Samples
CT96-100712CTC 90011RB9669287
CT96-33467CTC 90052RB9752011
CTC 1119CV 66543RB9759522
CTC 1410IACSP95-50001RB9854761
CTC 1512RB8550023SP80-32809
CTC 177RB85515633SP83-284731
CTC 219RB8555364SP83-50732
CTC 2026RB8675159RB9656211
CTC 222RB9356214SP91-10491
CTC 434RB9656211Various26
CTC 71RB96590212
Figure A1. Principal component scores of different sugarcane varieties for spectral data of skin scanning of sugarcane billets (a), cross-sectional scanning of sugarcane billets (b), defibrated cane samples (c), and raw juice samples (d).
Figure A1. Principal component scores of different sugarcane varieties for spectral data of skin scanning of sugarcane billets (a), cross-sectional scanning of sugarcane billets (b), defibrated cane samples (c), and raw juice samples (d).
Sensors 21 02195 g0a1
Figure A2. The plot of standardized residuals versus predicted sugarcane quality values from vis-NIR by skin (SS) and cross-sectional (CSS) scanning of billets, defibrated (DF), and raw juice (RJ) samples. Black dots are calibration samples, and red dots are external validation samples.
Figure A2. The plot of standardized residuals versus predicted sugarcane quality values from vis-NIR by skin (SS) and cross-sectional (CSS) scanning of billets, defibrated (DF), and raw juice (RJ) samples. Black dots are calibration samples, and red dots are external validation samples.
Sensors 21 02195 g0a2
Figure A3. The three first Partial least squares loadings for Brix prediction using near-infrared reflectance spectroscopy from the skin (a) and cross-sectional (b) scanning of billets, defibrated cane (c), and extracted raw juice (d).
Figure A3. The three first Partial least squares loadings for Brix prediction using near-infrared reflectance spectroscopy from the skin (a) and cross-sectional (b) scanning of billets, defibrated cane (c), and extracted raw juice (d).
Sensors 21 02195 g0a3
Figure A4. The three first Partial least squares loadings for Pol prediction using near-infrared reflectance spectroscopy from the skin (a) and cross-sectional (b) scanning of billets, defibrated cane (c), and extracted raw juice (d).
Figure A4. The three first Partial least squares loadings for Pol prediction using near-infrared reflectance spectroscopy from the skin (a) and cross-sectional (b) scanning of billets, defibrated cane (c), and extracted raw juice (d).
Sensors 21 02195 g0a4
Figure A5. The three first Partial least squares loadings for Fibre prediction using near-infrared reflectance spectroscopy from the skin (a) and cross-sectional (b) scanning of billets, and defibrated cane (c).
Figure A5. The three first Partial least squares loadings for Fibre prediction using near-infrared reflectance spectroscopy from the skin (a) and cross-sectional (b) scanning of billets, and defibrated cane (c).
Sensors 21 02195 g0a5
Figure A6. The three first Partial least squares loadings for Pol of cane prediction using near-infrared reflectance spectroscopy from the skin (a) and cross-sectional (b) scanning of billets, defibrated cane (c), and extracted raw juice (d).
Figure A6. The three first Partial least squares loadings for Pol of cane prediction using near-infrared reflectance spectroscopy from the skin (a) and cross-sectional (b) scanning of billets, defibrated cane (c), and extracted raw juice (d).
Sensors 21 02195 g0a6
Figure A7. The three first Partial least squares loadings for total recoverable sugars (TRS) prediction using near-infrared reflectance spectroscopy from the skin (a) and cross-sectional (b) scanning of billets, defibrated cane (c), and extracted raw juice (d).
Figure A7. The three first Partial least squares loadings for total recoverable sugars (TRS) prediction using near-infrared reflectance spectroscopy from the skin (a) and cross-sectional (b) scanning of billets, defibrated cane (c), and extracted raw juice (d).
Sensors 21 02195 g0a7

References

  1. Sexton, J.; Everingham, Y.; Donald, D.; Staunton, S.; White, R. Investigating the identification of atypical sugarcane using NIR analysis of online mill data. Comput. Electron. Agric. 2020, 168, 105111. [Google Scholar] [CrossRef]
  2. Phetpan, K.; Udompetaikul, V.; Sirisomboon, P. An online visible and near-infrared spectroscopic technique for the real-time evaluation of the soluble solids content of sugarcane billets on an elevator conveyor. Comput. Electron. Agric. 2018, 154, 460–466. [Google Scholar] [CrossRef]
  3. Rodrigues, F.A.; Magalhães, P.S.G.; Franco, H.C.J. Soil attributes and leaf nitrogen estimating sugar cane quality parameters: Brix, pol and fibre. Precis. Agric. 2013, 14, 270–289. [Google Scholar] [CrossRef]
  4. Ramírez-Morales, I.; Rivero, D.; Fernández-Blanco, E.; Pazos, A. Optimization of NIR calibration models for multiple processes in the sugar industry. Chemom. Intell. Lab. Syst. 2016, 159, 45–57. [Google Scholar] [CrossRef]
  5. Cortés, V.; Blasco, J.; Aleixos, N.; Cubero, S.; Talens, P. Monitoring strategies for quality control of agricultural products using visible and near-infrared spectroscopy: A review. Trends Food Sci. Technol. 2019, 85, 138–148. [Google Scholar] [CrossRef]
  6. Walsh, K.B.; Blasco, J.; Zude-Sasse, M.; Sun, X. Visible-NIR ‘point’ spectroscopy in postharvest fruit and vegetable assessment: The science behind three decades of commercial use. Postharvest Biol. Technol. 2020, 168, 111246. [Google Scholar] [CrossRef]
  7. Workman, J., Jr.; Weyer, L. Practical Guide and Spectral Atlas for Interpretive Near-Infrared Spectroscopy; CRC Press: Boca Raton, FL, USA, 2012; ISBN 9781439875261. [Google Scholar]
  8. Pasquini, C. Near infrared spectroscopy: A mature analytical technique with new perspectives—A review. Anal. Chim. Acta 2018, 1026, 8–36. [Google Scholar] [CrossRef] [PubMed]
  9. Taira, E.; Ueno, M.; Kawamitsu, Y. Automated quality evaluation system for net and gross sugarcane samples using near infrared spectroscopy. J. Near Infrared Spectrosc. 2010, 18, 209–215. [Google Scholar] [CrossRef]
  10. Sorol, N.; Arancibia, E.; Bortolato, S.A.; Olivieri, A.C. Visible/near infrared-partial least-squares analysis of Brix in sugar cane juice. A test field for variable selection methods. Chemom. Intell. Lab. Syst. 2010, 102, 100–109. [Google Scholar] [CrossRef]
  11. Valderrama, P.; Braga, J.W.B.; Poppi, R.J. Validation of multivariate calibration models in the determination of sugar cane quality parameters by near infrared spectroscopy. J. Braz. Chem. Soc. 2007, 18, 259–266. [Google Scholar] [CrossRef] [Green Version]
  12. Taira, E.; Ueno, M.; Furukawa, N.; Tasaki, A.; Komaki, Y.; Nagai, J.I.; Saengprachatanarug, K. Networking system employing near infrared spectroscopy for sugarcane payment in Japan. J. Near Infrared Spectrosc. 2013, 21, 477–483. [Google Scholar] [CrossRef]
  13. Taira, E.; Ueno, M.; Saengprachatanarug, K.; Kawamitsu, Y. Direct sugar content analysis for whole stalk sugarcane using a portable near infrared instrument. J. Near Infrared Spectrosc. 2013, 21, 281–287. [Google Scholar] [CrossRef]
  14. Phuphaphud, A.; Saengprachatanarug, K.; Posom, J.; Maraphum, K.; Taira, E. Non-destructive and rapid measurement of sugar content in growing cane stalks for breeding programmes using visible-near infrared spectroscopy. Biosyst. Eng. 2020, 197, 76–90. [Google Scholar] [CrossRef]
  15. Phuphaphud, A.; Saengprachatanarug, K.; Posom, J.; Maraphum, K.; Taira, E. Prediction of the fibre content of sugarcane stalk by direct scanning using visible-shortwave near infrared spectroscopy. Vib. Spectrosc. 2019, 101, 71–80. [Google Scholar] [CrossRef]
  16. Nawi, N.M.; Chen, G.; Jensen, T. In-field measurement and sampling technologies for monitoring quality in the sugarcane industry: A review. Precis. Agric. 2014, 15, 684–703. [Google Scholar] [CrossRef] [Green Version]
  17. Bramley, R.G.V. Lessons from nearly 20 years of Precision Agriculture research, development, and adoption as a guide to its appropriate application. Crop Pasture Sci. 2009, 60, 197–217. [Google Scholar] [CrossRef]
  18. Bramley, R.G.V.; Panitz, J.H.; Jensen, T.A.; Baillie, C.P. Within block spatial variation in CCS—Another potentially important consideration in the application of precision agriculture to sugarcane production. Int. Sugar J. 2013, 115, 1–8. [Google Scholar]
  19. Corrêdo, L.d.P.; Canata, T.F.; Maldaner, L.F.; de Lima, J.J.A.; Molin, J.P. Sugarcane Harvester for In-field Data Collection: State of the Art, Its Applicability and Future Perspectives. Sugar Tech 2020, 22, 1–14. [Google Scholar] [CrossRef]
  20. Nawi, N.M.; Chen, G.; Jensen, T.; Mehdizadeh, S.A. Prediction and classification of sugar content of sugarcane based on skin scanning using visible and shortwave near infrared. Biosyst. Eng. 2013, 115, 154–161. [Google Scholar] [CrossRef]
  21. Nawi, N.M.; Chen, G.; Jensen, T. Visible and shortwave near infrared spectroscopy for predicting sugar content of sugarcane based on a cross-sectional scanning method. J. Near Infrared Spectrosc. 2013, 21, 289–297. [Google Scholar] [CrossRef] [Green Version]
  22. Udompetaikul, V.; Phetpan, K.; Sirisomboon, P. Development of the partial least-squares model to determine the soluble solids content of sugarcane billets on an elevator conveyor. Meas. J. Int. Meas. Confed. 2021, 167, 107898. [Google Scholar] [CrossRef]
  23. Maraphum, K.; Chuan-Udom, S.; Saengprachatanarug, K.; Wongpichet, S.; Posom, J.; Phuphaphud, A.; Taira, E. Effect of waxy material and measurement position of a sugarcane stalk on the rapid determination of Pol value using a portable near infrared instrument. J. Near Infrared Spectrosc. 2018, 26, 287–296. [Google Scholar] [CrossRef]
  24. Phuphaphud, A.; Saengprachatanarug, K.; Posom, J.; Wongpichet, S.; Maraphum, K.; Taira, E. Effects of waxy types of a sugarcane stalk surface on the spectral characteristics of visible-shortwave near infrared measurement. Eng. J. 2019, 23, 13–24. [Google Scholar] [CrossRef]
  25. Cardozo, N.P.; Sentelhas, P.C. Climatic effects on sugarcane ripening under the influence of cultivars and crop age. Sci. Agric. 2013, 70, 449–456. [Google Scholar] [CrossRef] [Green Version]
  26. CONSECANA National Concill of Sugarcane Producers of São Paulo State. Instruction Manual, 6th ed.; CONSECANA-SP: Piracicaba, SP, Brazil, 2015; 81p. (In Portuguese) [Google Scholar]
  27. Franceschini, M.H.D.; Demattê, J.A.M.; Kooistra, L.; Bartholomeus, H.; Rizzo, R.; Fongaro, C.T.; Molin, J.P. Effects of external factors on soil reflectance measured on-the-go and assessment of potential spectral correction through orthogonalisation and standardisation procedures. Soil Tillage Res. 2018, 177, 19–36. [Google Scholar] [CrossRef]
  28. Oliveri, P.; Malegori, C.; Simonetti, R.; Casale, M. The impact of signal pre-processing on the final interpretation of analytical outcomes—A tutorial. Anal. Chim. Acta 2019, 1058, 9–17. [Google Scholar] [CrossRef] [PubMed]
  29. Barnes, R.J.; Dhanoa, M.S.; Lister, S.J. Standard normal variate transformation and de-trending of near-infrared diffuse reflectance spectra. Appl. Spectrosc. 1989, 43, 772–777. [Google Scholar] [CrossRef]
  30. Savitzky, A.; Golay, M.J.E. Smoothing and Differentiation of Data by Simplified Least Squares Procedures. Anal. Chem. 1964, 36, 1627–1639. [Google Scholar] [CrossRef]
  31. Gerretzen, J.; Szymańska, E.; Jansen, J.J.; Bart, J.; Van Manen, H.J.; Van Den Heuvel, E.R.; Buydens, L.M.C. Simple and Effective Way for Data Preprocessing Selection Based on Design of Experiments. Anal. Chem. 2015, 87, 12096–12103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Valderrama, P.; Braga, J.W.B.; Poppi, R.J. Variable selection, outlier detection, and figures of merit estimation in a partial least-squares regression multivariate calibration model. A case study for the determination of quality parameters in the alcohol industry by near-infrared spectroscopy. J. Agric. Food Chem. 2007, 55, 8331–8338. [Google Scholar] [CrossRef]
  33. The American Society for Testing and Materials (ASTM). ASTM E1655-17 Standard Practices for Infrared Multivariate Quantitative Analysis; ASTM Annual Book of Standards: West Conshohocken, PA, USA, 2017; Volume 5, 30p. [Google Scholar]
  34. Van der Voet, H. Comparing the predictive accuracy of models using a simple randomization test. Chemom. Intell. Lab. Syst. 1994, 25, 313–323. [Google Scholar] [CrossRef]
  35. Santana, F.B.; Giuseppe, L.O.; Souza, A.M.; Poppi, R.J. Removing the moisture effect in soil organic matter determination using NIR spectroscopy and PLSR with external parameter orthogonalization. Microchem. J. 2019, 145, 1094–1101. [Google Scholar] [CrossRef]
  36. Olivieri, A.C. Practical guidelines for reporting results in single- and multi-component analytical calibration: A tutorial. Anal. Chim. Acta 2015, 868, 10–22. [Google Scholar] [CrossRef] [PubMed]
  37. Pagani, V.; Stella, T.; Guarneri, T.; Finotto, G.; van den Berg, M.; Marin, F.R.; Acutis, M.; Confalonieri, R. Forecasting sugarcane yields using agro-climatic indicators and Canegro model: A case study in the main production region in Brazil. Agric. Syst. 2017, 154, 45–52. [Google Scholar] [CrossRef]
  38. Wang, J.; Zhao, T.; Yang, B.; Zhang, S. Sucrose Metabolism and Regulation in Sugarcane. J. Plant Physiol. Pathol. 2017, 5, 2. [Google Scholar] [CrossRef] [Green Version]
  39. Mancini, M.C.; Leite, D.C.; Perecin, D.; Bidóia, M.A.P.; Xavier, M.A.; Landell, M.G.A.; Pinto, L.R. Characterization of the Genetic Variability of a Sugarcane Commercial Cross Through Yield Components and Quality Parameters. Sugar Tech 2012, 14, 119–125. [Google Scholar] [CrossRef]
  40. Solomon, S. Post-harvest deterioration of sugarcane. Sugar Tech 2009, 11, 109–123. [Google Scholar] [CrossRef]
  41. Bazar, G.; Kovacs, Z.; Tsenkova, R. Evaluating spectral signals to identify spectral error. PLoS ONE 2016, 11, e0146249. [Google Scholar] [CrossRef] [PubMed]
  42. Steidle Neto, A.J.; Lopes, D.C.; Toledo, J.V.; Zolnier, S.; Silva, T.G.F. Classification of sugarcane varieties using visible/near infrared spectral reflectance of stalks and multivariate methods. J. Agric. Sci. 2018, 156, 537–546. [Google Scholar] [CrossRef]
  43. da Costa, M.V.A.; Fontes, C.H.; Carvalho, G.; de Moraes Júnior, E.C. UltraBrix: A Device for Measuring the Soluble Solids Content in Sugarcane. Sustainability 2021, 13, 1227. [Google Scholar] [CrossRef]
  44. Uys, L.; Botha, F.C.; Hofmeyr, J.H.S.; Rohwer, J.M. Kinetic model of sucrose accumulation in maturing sugarcane culm tissue. Phytochemistry 2007, 68, 2375–2392. [Google Scholar] [CrossRef]
  45. Wang, J.; Nayak, S.; Koch, K.; Ming, R. Carbon partitioning in sugarcane (Saccharum species). Front. Plant Sci. 2013, 4, 201. [Google Scholar] [CrossRef] [Green Version]
  46. Inarkar, M.B.; Lele, S.S. Extraction and Characterization of Sugarcane Peel Wax. ISRN Agron. 2012, 2012, 1–6. [Google Scholar] [CrossRef] [Green Version]
  47. Sexton, J.; Everingham, Y.; Donald, D.; Staunton, S.; White, R. A comparison of non-linear regression methods for improved on-line near infrared spectroscopic analysis of a sugarcane quality measure. J. Near Infrared Spectrosc. 2018, 26, 297–310. [Google Scholar] [CrossRef]
  48. Roger, J.M.; Boulet, J.C. A review of orthogonal projections for calibration. J. Chemom. 2018, 32, 1–16. [Google Scholar]
  49. Osborne, B.G. Near-Infrared Spectroscopy in Food Analysis. In Encyclopedia of Analytical Chemistry; John Wiley & Sons, Ltd.: Chichester, UK, 2006. [Google Scholar]
  50. Steidle Neto, A.J.; Toledo, J.V.; Zolnier, S.; Lopes, D.C.; Pires, C.V.; Silva, T.G.F. Prediction of mineral contents in sugarcane cultivated under saline conditions based on stalk scanning by Vis/NIR spectral reflectance. Biosyst. Eng. 2017, 156, 17–26. [Google Scholar] [CrossRef]
  51. Golic, M.; Walsh, K.; Lawson, P. Short-wavelength near-infrared spectra of sucrose, glucose, and fructose with respect to sugar concentration and temperature. Appl. Spectrosc. 2003, 57, 139–145. [Google Scholar] [CrossRef]
  52. Pomares-Viciana, T.; Martínez-Valdivieso, D.; Font, R.; Gómez, P.; del Río-Celestino, M. Characterisation and prediction of carbohydrate content in zucchini fruit using near infrared spectroscopy. J. Sci. Food Agric. 2018, 98, 1703–1711. [Google Scholar] [CrossRef]
  53. Bázár, G.; Romvári, R.; Szabó, A.; Somogyi, T.; Éles, V.; Tsenkova, R. NIR detection of honey adulteration reveals differences in water spectral pattern. Food Chem. 2016, 194, 873–880. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The sequence samples collection and spectral measurements during the preparation steps of samples for conventional analysis. (a) Sampling of sugarcane billets by an oblique probe in the cargo truck; (b) sugarcane billets for skin and cross-sectional scanning measurements; (c) milling and homogenization of sugarcane to defibrated sample; (d) defibration sample; (e) pressing of defibrated sample to extracting of juice; (f) extracted raw juice; prepared samples for vis-NIR spectral measurements: (g) cross-sectional and skin of billets inside pipeline chambers, (h) defibrated cane and raw juice; (i) Veris vis-NIR spectrometer and internal configuration scheme of the measurement shank; (j) fibrous cane residue and extracted raw juice for conventional analysis.
Figure 1. The sequence samples collection and spectral measurements during the preparation steps of samples for conventional analysis. (a) Sampling of sugarcane billets by an oblique probe in the cargo truck; (b) sugarcane billets for skin and cross-sectional scanning measurements; (c) milling and homogenization of sugarcane to defibrated sample; (d) defibration sample; (e) pressing of defibrated sample to extracting of juice; (f) extracted raw juice; prepared samples for vis-NIR spectral measurements: (g) cross-sectional and skin of billets inside pipeline chambers, (h) defibrated cane and raw juice; (i) Veris vis-NIR spectrometer and internal configuration scheme of the measurement shank; (j) fibrous cane residue and extracted raw juice for conventional analysis.
Sensors 21 02195 g001
Figure 2. Mean and standard deviation (SD) of annual variation of sugarcane quality parameters, Brix, Pol, Fibre, and total recoverable sugars (TRS); and spectral data collection periods (vertical bars).
Figure 2. Mean and standard deviation (SD) of annual variation of sugarcane quality parameters, Brix, Pol, Fibre, and total recoverable sugars (TRS); and spectral data collection periods (vertical bars).
Sensors 21 02195 g002
Figure 3. Correlogram of sugarcane quality parameters with frequency distributions on the principal diagonal, Pearson’s correlation coefficient and respective p-values, and the correlation trend line.
Figure 3. Correlogram of sugarcane quality parameters with frequency distributions on the principal diagonal, Pearson’s correlation coefficient and respective p-values, and the correlation trend line.
Sensors 21 02195 g003
Figure 4. vis-NIR mean spectra and standard deviation (SD) of all 302 sugarcane samples for (a) skin and (b) cross-sectional scanning of billets, (c) defibrated, and (d) raw juice samples.
Figure 4. vis-NIR mean spectra and standard deviation (SD) of all 302 sugarcane samples for (a) skin and (b) cross-sectional scanning of billets, (c) defibrated, and (d) raw juice samples.
Sensors 21 02195 g004
Figure 5. PCA score plot for the sugarcane sample types analyzed. SS-skin scanning of billets; CSS-cross-sectional scanning of billets; DF-defibrated samples; RJ-raw juice samples.
Figure 5. PCA score plot for the sugarcane sample types analyzed. SS-skin scanning of billets; CSS-cross-sectional scanning of billets; DF-defibrated samples; RJ-raw juice samples.
Sensors 21 02195 g005
Figure 6. Plots of observed values versus predicted sugarcane quality values from vis-NIR by skin (SS) and cross-sectional (CSS) scanning of billets, defibrated (DF), and raw juice (RJ) samples. Brix, Pol, Fibre, and Pol of cane are in percentage, and TRS values are in kg Mg−1.
Figure 6. Plots of observed values versus predicted sugarcane quality values from vis-NIR by skin (SS) and cross-sectional (CSS) scanning of billets, defibrated (DF), and raw juice (RJ) samples. Brix, Pol, Fibre, and Pol of cane are in percentage, and TRS values are in kg Mg−1.
Sensors 21 02195 g006
Figure 7. Heatmap of Variable Importance in Projection (VIP) for models used to predict Brix, Pol, Fibre, Pol of cane, and TRS based on different spectral sample types datasets. SS: skin scanning of billets; CSS: cross-sectional scanning of billets; DF: defibrated samples; RJ: raw juice samples; TRS–total recoverable sugars.
Figure 7. Heatmap of Variable Importance in Projection (VIP) for models used to predict Brix, Pol, Fibre, Pol of cane, and TRS based on different spectral sample types datasets. SS: skin scanning of billets; CSS: cross-sectional scanning of billets; DF: defibrated samples; RJ: raw juice samples; TRS–total recoverable sugars.
Sensors 21 02195 g007
Table 1. Descriptive statistics of the reference results for the sugarcane quality attributes of all samples, calibration, and external validation data sets.
Table 1. Descriptive statistics of the reference results for the sugarcane quality attributes of all samples, calibration, and external validation data sets.
ParameterunitMean ± SDMedianRangep25p75SEL
All samples (n = 302)
Brix%18.95 ± 1.7118.9913.08–23.4217.8020.010.03
Pol%16.67 ± 1.9016.6610.78–21.2015.4117.950.01
Fibre%13.29 ± 1.7912.907.22–20.0812.0714.330.07
Pol of cane%13.80 ± 1.5613.918.40–17.5612.9214.780.01
TRSkg Mg−1137.66 ± 14.48138.6686.94–173.80129.75146.841.12
Calibration set (n = 227)
Brix%18.86 ± 1.6618.8013.08–23.4217.7919.98-
Pol%16.54 ± 1.8616.5510.78–21.2015.3817.78-
Fibre%13.31 ± 1.8912.837.22–20.0812.0514.41-
Pol of cane%13.69 ± 1.5213.798.40–17.5612.8314.60-
TRSkg Mg−1136.65 ± 14.07137.0186.94–173.80128.99145.19-
Validation set (n = 75)
Brix%19.24 ± 1.8519.5913.55–23.0518.0620.61-
Pol%17.06 ± 1.9817.3311.24–20.9015.7318.54-
Fibre%13.23 ± 1.4413.0210.49–17.1512.2414.16-
Pol of cane%14.14 ± 1.6414.348.96–17.1413.2715.47-
TRSkg Mg−1140.76 ± 15.36142.5992.16–169.02132.35152.44-
SD: standard deviation; p25: lower quartile; p75: upper quartile; SEL: standard error of laboratory; TRS: total recoverable sugar.
Table 2. Figures of merit for the PLSR models for all studied sugarcane quality attributes and sample types.
Table 2. Figures of merit for the PLSR models for all studied sugarcane quality attributes and sample types.
AttributeSample TypeLVRMSEC aRMSECV aRMSEP aR2cR2pRPIQ
BrixSS90.921.101.290.640.481.98
CSS60.951.041.380.620.411.85
DF70.670.750.840.810.803.05
RJ80.640.830.750.850.853.39
PolSS81.091.261.420.600.481.98
CSS61.091.191.440.610.441.95
DF70.820.930.870.790.833.24
RJ70.800.970.900.820.813.12
FibreSS101.021.290.870.590.652.22
CSS41.451.501.270.240.231.51
DF50.931.040.820.690.692.36
RJ b-------
Pol of caneSS70.951.071.130.520.461.94
CSS51.011.091.270.520.311.73
DF70.730.840.720.760.813.04
RJ70.710.850.720.780.813.07
TRSSS98.5710.2710.860.600.501.85
CSS59.4910.1711.860.500.321.69
DF76.507.516.710.760.822.99
RJ76.387.956.790.780.812.96
SS: skin scanning of billets samples; CSS: cross-sectional scanning of billets samples; DF: defibrated samples; RJ: raw juice samples. LV: latent variable. RMSEC: root mean square error of calibration. RMSECV: Root Mean Square Error of Cross-Validation. RMSEP: Root Mean Square Error of Prediction. R2c: calibration coefficient of determination. R2p: prediction coefficient of determination. RPIQ: Ratio of performance to interquartile distance. a values for Brix, Pol, Fibre, and Pol of cane are in percentage and TRS in kg Mg−1. b The fibre content was not determined from raw juice samples.
Table 3. p-Values of randomization test of external validation set for all compared sugarcane sample types.
Table 3. p-Values of randomization test of external validation set for all compared sugarcane sample types.
Binary Combination
(Sample Types)
Sugarcane Quality Parameters
BrixPolFibrePol of caneTRS
SS vs. CSS0.1040.116<0.0010.0360.008
SS vs. DF1.001.000.6671.001.00
SS vs. RJ1.001.00-1.001.00
CSS vs. DF1.001.001.001.001.00
CSS vs. RJ1.001.00-1.001.00
DF vs. RJ0.8790.344-0.6060.502
SS: skin scanning of billets; CSS: cross-sectional scanning of billets; DF: defibrated samples; RJ: raw juice samples; TRS: total recoverable sugar.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Corrêdo, L.d.P.; Maldaner, L.F.; Bazame, H.C.; Molin, J.P. Evaluation of Minimum Preparation Sampling Strategies for Sugarcane Quality Prediction by vis-NIR Spectroscopy. Sensors 2021, 21, 2195. https://doi.org/10.3390/s21062195

AMA Style

Corrêdo LdP, Maldaner LF, Bazame HC, Molin JP. Evaluation of Minimum Preparation Sampling Strategies for Sugarcane Quality Prediction by vis-NIR Spectroscopy. Sensors. 2021; 21(6):2195. https://doi.org/10.3390/s21062195

Chicago/Turabian Style

Corrêdo, Lucas de Paula, Leonardo Felipe Maldaner, Helizani Couto Bazame, and José Paulo Molin. 2021. "Evaluation of Minimum Preparation Sampling Strategies for Sugarcane Quality Prediction by vis-NIR Spectroscopy" Sensors 21, no. 6: 2195. https://doi.org/10.3390/s21062195

APA Style

Corrêdo, L. d. P., Maldaner, L. F., Bazame, H. C., & Molin, J. P. (2021). Evaluation of Minimum Preparation Sampling Strategies for Sugarcane Quality Prediction by vis-NIR Spectroscopy. Sensors, 21(6), 2195. https://doi.org/10.3390/s21062195

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