Next Article in Journal
Comparative Analysis of Chinese HJ-1 CCD, GF-1 WFV and ZY-3 MUX Sensor Data for Leaf Area Index Estimations for Maize
Next Article in Special Issue
Recognition of Wheat Spike from Field Based Phenotype Platform Using Multi-Sensor Fusion and Improved Maximum Entropy Segmentation Algorithms
Previous Article in Journal
Preliminary Study of Soil Available Nutrient Simulation Using a Modified WOFOST Model and Time-Series Remote Sensing Observations
Previous Article in Special Issue
Estimation of Wheat LAI at Middle to High Levels Using Unmanned Aerial Vehicle Narrowband Multispectral Imagery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Comparison of Regression Techniques for Estimation of Above-Ground Winter Wheat Biomass Using Near-Surface Spectroscopy

1
Key Laboratory of Quantitative Remote Sensing in Agriculture of Ministry of Agriculture China, Beijing Research Center for Information Technology in Agriculture, Beijing 100097, China
2
International Institute for Earth System Science, Nanjing University, Nanjing 210023, China
3
National Engineering Research Center for Information Technology in Agriculture, Beijing 100097, China
4
Beijing Engineering Research Center for Agriculture Internet of Things, Beijing 100097, China
*
Authors to whom correspondence should be addressed.
Both authors contributed equally to this work.
Remote Sens. 2018, 10(1), 66; https://doi.org/10.3390/rs10010066
Submission received: 16 November 2017 / Revised: 25 December 2017 / Accepted: 2 January 2018 / Published: 5 January 2018

Abstract

:
Above-ground biomass (AGB) provides a vital link between solar energy consumption and yield, so its correct estimation is crucial to accurately monitor crop growth and predict yield. In this work, we estimate AGB by using 54 vegetation indexes (e.g., Normalized Difference Vegetation Index, Soil-Adjusted Vegetation Index) and eight statistical regression techniques: artificial neural network (ANN), multivariable linear regression (MLR), decision-tree regression (DT), boosted binary regression tree (BBRT), partial least squares regression (PLSR), random forest regression (RF), support vector machine regression (SVM), and principal component regression (PCR), which are used to analyze hyperspectral data acquired by using a field spectrophotometer. The vegetation indexes (VIs) determined from the spectra were first used to train regression techniques for modeling and validation to select the best VI input, and then summed with white Gaussian noise to study how remote sensing errors affect the regression techniques. Next, the VIs were divided into groups of different sizes by using various sampling methods for modeling and validation to test the stability of the techniques. Finally, the AGB was estimated by using a leave-one-out cross validation with these powerful techniques. The results of the study demonstrate that, of the eight techniques investigated, PLSR and MLR perform best in terms of stability and are most suitable when high-accuracy and stable estimates are required from relatively few samples. In addition, RF is extremely robust against noise and is best suited to deal with repeated observations involving remote-sensing data (i.e., data affected by atmosphere, clouds, observation times, and/or sensor noise). Finally, the leave-one-out cross-validation method indicates that PLSR provides the highest accuracy (R2 = 0.89, RMSE = 1.20 t/ha, MAE = 0.90 t/ha, NRMSE = 0.07, CV (RMSE) = 0.18); thus, PLSR is best suited for works requiring high-accuracy estimation models. The results indicate that all these techniques provide impressive accuracy. The comparison and analysis provided herein thus reveals the advantages and disadvantages of the ANN, MLR, DT, BBRT, PLSR, RF, SVM, and PCR techniques and can help researchers to build efficient AGB-estimation models.

Graphical Abstract

1. Introduction

Accurate estimates of crop biophysical variables are crucial for monitoring vegetation growth and for analyzing important physiological parameters during the crop growth cycle [1,2]. One such variable, above-ground biomass (AGB), plays an important role in plant functioning because it reflects the status of crop growth and is related to solar-energy consumption, yield, and grain quality [3,4]. Therefore, AGB is considered as one of the most important crop biophysical parameters, and its accurate estimation can help improve crop monitoring and yield prediction [5]. Traditional AGB estimates are based on destructive measurements, which are not only time and labor consuming, but more importantly, are difficult to apply over large areas [6]. In recent years, Hyperspectral remote-sensing data acquired from the ground [7,8], unmanned aerial vehicles [9], airborne platforms [10,11,12], and satellite platforms [13] have been able to capture crop canopy spectra in narrow bands and thereby provide information on the biochemical composition of the canopy. Crop physiology research shows that spectral absorption by plant leaves is mainly due to the leaf pigments, especially chlorophyll content (Chl) [14]. The reflectance is low in both the blue and red regions of the spectrum, due to absorption by chlorophyll for photosynthesis; it has a peak at the green region which gives rise to the green color of vegetation [15]. In the near-infrared region, the reflectance is much higher than that in the visible band due to the cellular structure in the leaves [16]. Previous studies have shown that near-infrared- and red-band vegetation indexes (VIs) are effective for estimating AGB [8,9,11]. However, during the reproductive growth of crops, with the senescence of leaves, the effectiveness of photosynthesis is reduced [14,17]. With clear decreases in both photosynthesis and the near-infrared reflectance, the correlation between AGB and the red- and near-infrared-based VIs reduced. Therefore, hyperspectral remote sensing of AGB has received increasing attention as an efficient and precise method for nondestructive monitoring in agricultural research [18].
Physically based models and empirical regression techniques are two essential approaches for estimating vegetation characteristics from hyperspectral measurements [19]. Physically based models were founded on physical principles. The two main examples of this approach are radiative transfer (RT) models and geometric optical models [19]. Because vegetation canopy reflectance depends on a number of factors [20] (e.g., leaf-area index, Chl, water content, matter content, soil reflectance, and bidirectional reflectance distribution function), physically based models require canopy biophysical parameters, soil parameters, and some external parameters to simulate canopy reflectance, and these are often not readily available. In contrast, empirical regression techniques require a large number of ground measurements, and offer a direct relationship between spectral features and vegetation parameters. Previous research has used many powerful empirical regression techniques that make full use of the narrow hyperspectral bands, VIs, and even different types of sensor data [21]. These techniques essentially fall into two categories: (i) machine-learning techniques such as artificial neural network (ANN) [22], decision tree regression (DT) [23], boosted binary regression tree (BBRT) [24], random forest regression (RF) [25], support vector machine regression (SVM) [26], and (ii) conventional regression techniques such as multivariable linear regression (MLR) [26,27], partial least squares regression (PLSR) [7,8,22], and principal component regression (PCR) [7]. Many studies have obtained promising results by using these techniques [8,9,10,11,26]. However, hyperspectral data redundancy is a big problem because of the high spectral dimensions and large number of bands [28]. In addition, the correlation between the spectral and AGB vary with the crop growth period, which is related to the physiological state of the crop [17]. To address this problem, many researchers have tried to extract features from narrow hyperspectral bands first, and many methods to do this have been proposed; for example, correlation analysis, continuum removal [29], red-edge position [30], gray relational analysis [31], and out-of-bag analysis [21]. Spectral vegetation indexes (VIs) have been widely used for decades, and more than 60 VIs [32] have been proposed for estimating biophysical variables [33].
Conventional regression techniques are more suitable for data that have a clear linear or exponential relationship with a distinct estimation equation, whereas machine-learning techniques are typically better able to cope with the strong nonlinearity between the biophysical and biochemical parameters and the reflection spectra [34]. However, many studies indicate that empirical regression techniques are rarely transferable to other sites with different vegetation, or to data acquired from other types of sensors or under different acquisition conditions. Despite this, empirical regression techniques still have some advantages, such as fewer input variables, less computation, and ease of application, which have resulted in their widespread use under many conditions.
Numerous studies have used hyperspectral remote-sensing data and empirical regression techniques to estimate AGB [26], and some analyses of the performance of these techniques have also been carried out, although they focus mostly on comparing the estimation accuracy. No comprehensive study is available as yet that evaluates these regression techniques for estimating AGB, and no studies have evaluated the different statistical techniques to better understand their respective advantages and disadvantages.
The main objective of the present study is to evaluate the performance (in particular, data selection, sampling methods, noise immunity) of eight regression techniques for estimating AGB. The following four tests were applied:
(1)
VIs were used to train regression models, which were validated to select best VI input (Section 4.1).
(2)
The noise immunities of these techniques were compared by simulating remote-sensing errors by adding white Gaussian noise (Section 4.2).
(3)
The stability of these techniques was examined by using samples of varying sizes and different sampling methods for modeling and validation (Section 4.3).
(4)
Leave-one-out cross validation was used to evaluate the accuracy of the AGB estimation of these techniques (Section 4.4).
We discuss the performance of eight AGB estimation techniques and the advantages and disadvantages of each technique (Section 5), then summarize the optimal conditions for using these techniques.

2. Materials

2.1. Study Area

The study area was situated in Changping District, which is located in the northwest part of Beijing City, China (see Figure 1). Experiments were conducted at the National Precision Agriculture Research Center of China (116°26′36″E, 40°10′44″N). Changping District has an average altitude of 36 m, its total area is about 1352 km2, and it has a warm temperate semi-humid continental monsoon climate, with an average rainfall of 450 mm, an average low temperature of −10 to 7.5 °C, and an average high temperature of 35 to 40 °C.
The aim of the agronomy experiment was to increase the difference in AGB by using two crop varieties, three water treatments, and four nitrogen treatments. The AGB was measured by using ground-based techniques. The experiments involved two winter wheat cultivars, J9843 and ZM175, which are the main winter wheat varieties grown in northern China. The irrigation treatment included rainfall only (W0, see Figure 1), rainfall plus normal irrigation (W1, 100 mm), and rainfall plus double the normal irrigation (W2, 200 mm). The nitrogen fertilizer treatment included no fertilizer (N0), one-half the normal fertilization (N1, 195 kg/ha), normal fertilization (N2, 390 kg/ha), and twice the normal fertilization (N3, 780 kg/ha).

2.2. Measurement of Data

A 5000 m2 square area was selected as the experimental field (see Figure 1) and divided into 48 plots each of size 6 m × 8 m. In each plot, winter wheat near the center of the given plot was selected for spectral, physiological, and biochemical measurements and analyses. AGB, Chl and canopy spectral measurements were made at four growth stages: the winter wheat jointing stage (13 and 14 April 2015), the flag leaf stage (26 and 27 April 2015), the flowering period (12 to 14 May 2015), and the filling period (25 to 27 May 2015). Four ground-based measurements allowed 192 sets of Chl, winter wheat biomass, and canopy hyperspectral data to be collected.

2.2.1. Measurements of Winter Wheat Canopy Reflectance

Canopy hyperspectral reflectance was acquired by using an ASD FieldSpec 3 spectrometer (FieldSpec 3 spectrometer, Analytical Spectral Devices, Boulder, Colorado, CO, USA) from 10:00~14:00 (Beijing time, UTC/GMT+08:00) in windless and cloudless conditions. We calibrated the field spectrometer based on the reflectance from a 40 cm × 40 cm BaSO4 white board, and the vertical height from the canopy is 1.3 m. The winter wheat canopy reflectance was measured 10 times (the scanning time was 0.2 s) at the center of each plot, and the average reflectance was recorded. To reduce the influence of sky and field conditions on the spectral measurements, each plot was measured three times, and the mean value was used as the canopy reflectance for the given experimental plot. Figure 2a shows the average hyperspectral reflectance spectrum for the four growing stages.

2.2.2. Measurements of Winter Wheat Chlorophyll and Above-Ground Biomass

During measurements, the planting density of winter wheat (row spacing 15 cm) was investigated, and 20 stems were collected near the center of each plot. Chl was measured from the first and second uppermost leaves by using a Dualex 4 (Dualex Scientific Portable Sensor for Leaf Measurements, Force-a, Université Paris Sud, Orsay, France) and the average values were processed (see Figure 2b).
After ground measurements, the winter wheat organs were processed in the laboratory. They were first put into paper bags and dried at 80 °C to remove moisture, then, once the sample weight became constant (about 24 h), they were weighed by using a balance with an accuracy of 0.001 g. Finally, the biomass per unit area was calculated based on the measured planting density and sample dry weight. The winter wheat AGB was calculated by using
AGB = m × n 20 × l C
where m is the dry weight of the sample, n is the number of winter wheat ears per unit area, and l is the row spacing. The statistics of the AGB measurement for different growing periods is shown in Table 1.

3. Methods

Data selection, sampling methods, noise immunity, and prediction performance were analyzed by using a series of VIs and ground-based measurements of the AGB. The flowchart in Figure 3 illustrates the process.

3.1. Regression Techniques

3.1.1. Machine Learning Techniques

(1) Artificial neural networks have represented a hot research topic in artificial intelligence since the 1980s; this method is very powerful in dealing with nonlinear relationships [35]. An ANN is based on a collection of connected units called artificial neurons, and each neuron can transmit a signal to other neurons. ANN is composed of a large number of neurons, with each neuron representing a particular output function. The connection between two neurons represents the weighted value of the signal through the connection. The network outputs different weighted values and the incentive function, according to the network connection mode. (2) Support vector machines were proposed by Cortes and Vapnik [36] in 1995 and offer many unique advantages for dealing with complex multidimensional data. A SVM constitutes a supervised learning model with associated learning algorithms that analyze data used for classification and regression analysis. A SVM can be used as a regression method, maintaining all the main features that characterize the algorithm (i.e., maximal margin). Support vector regression (SVR) uses the same principles as a SVM for classification, with only a few minor differences. Herein, we use LIBSVM (LIBSVM 3.1.2—A Library for Support Vector Machines, Version 3.12 [37]) for the tests. (3) A decision tree is a tree structure in which each internal node represents a test of an attribute, each branch represents a test output, and each leaf node represents a category [38]. A decision node has two or more branches, each representing values for the attribute tested. Breaking down a dataset into smaller and smaller subsets incrementally develops an associated decision tree. The final result is a tree with decision nodes and leaf nodes. A leaf node represents a decision for the numerical target. (4) A boosted binary regression tree is a powerful regression method proposed by Friedman [24] in 2001. Boosted binary regression trees combine binary regression trees by using a gradient-boosting technique [39]. (5) Random forest regression is a data analysis and statistical method that is widely used in machine-learning research. It was proposed by Breiman and Cutler [25] in 2001 and is an ensemble learning method for classification and regression that operates by constructing a multitude of decision trees at training time and outputting the class that is the mode of the classes (classification) or the mean prediction (regression) of the individual trees. RF has a higher accuracy, better tolerance to outliers and noise, and makes excellent use of the full spectral information.

3.1.2. Conventional Regression Techniques

(1) Multiple linear regression is a regression method in which two or more independent variables are used to analyze a dependent variable. The regression equation is used to calculate the parameters by using the least squares method in which the sum of the errors squared is minimized. (2) Partial least squares regression is a data analysis method proposed by Wold [40] in 1966. PLSR has also been widely used in studies of vegetation because it provides an efficient way to make full use of hyperspectral information. Previous studies [8,9,10] indicate that PLSR makes excellent use of the full spectral information and is a flexible method for monitoring agricultural crop parameters. (3) Principal component analysis (PCA) is a technique to simplify data sets based on a linear transformation of data into a new coordinate system. After that transformation, the largest variance in the data projection appears in the first coordinate (called the first principal component), the second largest variance appears in second coordinate (second principal component), and so on. PCA often reduces the dimensionality of data sets. This method can reduce the dimensionality of hyperspectral data, thus avoiding the problem of collinear variables that can occur in PLSR and MLR regression [7].
In the present work, we analyze the ANN, SVM, RF, BBRT, DT, MLR, PLSR, and PCR regression by using Matlab2014a (Matrix Laboratory 2014a, MathWorks, Inc., Natick, MA, USA) on a Microsoft Windows platform.

3.2. Selection of Vegetation Indexes

A VI is a combination of two or more characteristic spectra acquired by multispectral or hyperspectral remote sensing. It is a simple, effective, and empirical measure of the surface vegetation status. VIs are widely used to classify vegetation and environmental changes, determine crop and forage yield, monitor droughts, etc. After many years of research on narrow-band hyperspectral spectra, incomplete statistics show that dozens of VIs exist that can be used to estimate biophysical parameters [13,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74].
Data redundancy and multi-collinearity can seriously affect regression performance. By using selected 54 VIs (Table 2) [32], the abilities of eight techniques to solve the multi-collinearity problem can be analyzed. Section 4.1 gives the best input VIs (we set seven levels: 5, 10, 15, 20, 30, 40 and 54) based on eight techniques by analyzing the modeling and validation as a function of the VIs used as the input.

3.3. Noise Simulation

Many error sources exist in remote-sensing imaging and sensor systems (see Figure 4), including radiation errors caused by the atmosphere, topography, or other geometric errors and systematic errors related to the charge-coupled device (CCD) sensor [78,79]. Although radiometric calibration and radiometric correction are applied to correct for sensor degradation and atmospheric effects, the noise cannot be completely removed. Noise such as shot noise due to the quantum properties of light and readout noise generated by the output amplifier remain, and both follow a Poisson distribution [78]. In addition, the dark-current noise and thermal noise are present and are proportional to the CCD temperature; these follow a Gaussian distribution [78,79].
To evaluate how sensor noise and other uncertainty sources affect these data-analysis techniques, we simulate the internal noise (dark current, random noise) of a CCD used for remote sensing. We add random noise to the validation VIs to analyze how remote-sensing noise affects the stability of these techniques with the help of white Gaussian noise [7]. These models to estimate the AGB use the original VIs and are validated by using VIs with noise added. In these tests, the signal-to-noise ratio (SNR) is set to 5, 10, 15, 20, 30, and 50 (the noise increases with decreasing SNR). We also compare the results to those obtained with noiseless VIs.

3.4. Modeling Parameters and Sampling Methods

Before the test, the ability to deal with multiple VIs (best input VIs) and to obtain the modeling parameters must be analyzed. In this test, the optimal modeling parameters are determined according to the accuracy of validation by using an exhaustive method within limits (Table 3). In the next modeling test, these optimal parameters were selected for testing. Through thousands of rough modeling and verification processes (Note in Table 3), we obtained the best optimal model parameters under the optimal VIs input.
Leave one sampling (LOS) was used to evaluate the performance of each technique, and global random sampling (GRS) and growth-period sampling (GPS) were used to evaluate the performance and stability of each technique with different sampling methods [81,82]. Global random sampling represents random samples from all samples, with a total of three samplings taken and denoted GRS1/3 (64 samples for modeling, the remaining 128 samples for validation), GRS1/2 (96 samples for modeling, the remaining 96 samples for validation), and GRS2/3 (128 samples for modeling, the remaining 64 samples for validation). Growth-period sampling uses random samples from each period and ensures an equal number of samples per birth period. Therefore, all samples from the whole growth period were divided into four layers, with each layer including one growth period. Again a total of three samplings were taken and are likewise denoted GPS1/3, GPS1/2, and GPS2/3. Leave one sampling represents leave-one-out cross-validation [83] in which only one sample is selected for verification, with all other samples taken as training samples (191 samples for modeling, the remaining one sample used for validation).

3.5. Precision Evaluation

We use the coefficient of determination R2, the root mean square error (RMSE), the mean absolute error (MAE), and normalized root mean square error (NRMSE) and coefficient of variance of the root mean square error [CV(RMSE)] to evaluate the accuracy of each technique. A larger R2 value corresponds to a smaller RMSE, MAE, NRMSE, and CV(RMSE) and greater model accuracy. R2, RMSE, MAE, NRMSE, and CV(RMSE) are calculated as follows:
R 2 = 1 i = 1 n ( y i x i ) 2 i = 1 n ( y i y ¯ ) 2
RMSE = i = 1 n ( x i y i ) 2 n
MAE = i = 1 n | x i y i | n
NRMSE = i = 1 n ( x i y i ) 2 n y max y min
CV ( RMSE ) = i = 1 n ( x i y i ) 2 n y ¯
where xi and yi are the estimated and measured AGB values, respectively, ymax and ymin are the maximum and minimum measured values, respectively, x ¯ and y ¯ are the average estimated and measured values, respectively, and n is the sample number.

4. Results

The correlation coefficients r between the AGB and the VIs are shown in Table 4. The results show that all the measured VIs are correlated with biomass to varying degrees. Of the 54 VIs investigated, NPQI performs the best (r = 0.757). The correlation with AGB of the red-, green-, and blue-band-based spectral VIs [e.g., NPQI (0.757), BGI (0.555), TCARI (0.188), NPCI (0.474), BRI (0.519), and MCARI (0.237)] is greater than that of red- and near-infrared-band-based spectral VIs [e.g., MND680(0.193), MND705(0.086), NDVI (0.039), SR (0.029), EVI2 (0.182), OSAVI (0.088), and EVI (0.149)].
The correlation coefficients r among the 54 VIs are shown in Figure 5. For each VI, there are 53 colors that represent different correlation coefficient values. The results (Figure 5) show that complex correlations exist among these 54 VIs. NPQI (first), BGI (second), ARI (16th) and TCARI (22nd) are low correlated with other VIs. Among the top 22 VIs, a low correlation (zone a) was observed, but the remaining 32 VIs have a high correlation (zone b). The VIs analysis results of total explained variance and variance inflation factor (VIF) values are shown in Table A1 and Table A2 (in Appendix A).

4.1. Selection of Vegetation Indexes

The best AGB models and the associated validation accuracy of the eight techniques are shown in Figure 6. After a different number of VIs (Table 5) were incorporated into the modeling, the validation accuracy of ANN, BBRT, DT, and RF (Figure 6a–c,g) flattens out. BBRT and DT performed well when using the top five VIs as the input, although the validation accuracy decreases slightly after using lower-correlation VIs as the input for modeling (Figure 6b,c). The performance of ANN, PCR, and SVM becomes complex (Figure 6a,f,h) after using lower-correlation VIs as the input for modeling.
The results, given in Table 5, are the optimum numbers of input VIs determined by the change of validation accuracy (Figure 6). The numbers were confirmed when the best estimation accuracy was obtained. For example, with the input modeling of the top 30 VIs, the highest accuracy was acquired (Figure 6a). Accordingly, for RF, the parameters ntree = 520 and mtry = 8, the ANN hidden layer is 10 and 2; c = 10 and g = −2.5 for LIBSVM. In PCR modeling, we use 85% as the threshold to determine the principal component, and when we use the top five VIs as the input, the cumulative variance of the top two principal components is 89.136%.
Considering that the correlation coefficients r among VIs were high (Figure 5), multi-collinearity may be a problem when using so many VIs for modeling. The VIs analysis results of total explained variance (Table A1) and VIF values (Table A2) also support this view. The results in Figure 6 show that different techniques have varying abilities in tackling multi-collinearity data. ANN, BBRT and RF show a good performance when dealing with collinear data (Figure 6a,b,g), which performed well in 30 VIs, 40 VIs and 54 VIs group modeling and validation. Greater modeling accuracy was obtained with MLR after using more VIs as the input for modeling (Figure 6d), but validation accuracy decreased, especially after more than 20 VIs were input (R2 (V), MAE (V) and RMSE (V) in Figure 6d).

4.2. Test with White Gaussian Noise

A comparative analysis of different estimation accuracies resulting from the eight selected techniques with white Gaussian noise (different SNR values) is presented in Figure 7. For each technique, the three figures (Figure 7) represent R2, RMSE and MAE with different SNR.
Figure 7 shows the noise immunity for the eight analytical techniques. The results indicate RF > SVM >DT > BBRT >ANN > PCR > PLSR > MLR. RF performs best in this test with a validation R2 near 0.2 and SNR = 5. For MLR, poor noise immunity was observed; its validation accuracy declined from SNR = 30 (Figure 7). Meanwhile, the validation of MLR (R2 (<0.2), RMSE (about 8 t/ha) and MAE (about 6 t/ha)) is the worst in eight techniques. As the noise increases (Figure 7, SNR <30), MLR, PCR and PLSR are extremely sensitive to it, whereas ANN, SVM, DT, and BBRT are more robust against noise; however, the latter techniques also show a poor performance with increased levels of noise (Figure 7, SNR = 5).

4.3. Stability Test

Figure 8 shows the results for modeling and validation using GRS1/3, GRS1/2, and GRS2/3 with the parameters given in Table 5. The calculated absolute difference of AGB − GRS modeling and the validation accuracy (we use ∇R2, ∇RMSE and ∇MAE to show the difference) are given in Table 6.
The difference in the modeling and validation results with eight techniques for global random sampling varies (Figure 8). For MLR and PLSR, a stable performance was observed with three sampling methods (Table 6). SVM and PCR have a poorer performance than MLR and PLSR in this test (Figure 8 and Table 6). For ANN, RF, DT and BBRT, the difference in modeling and validation accuracy is huge (Figure 8 and Table 6). As the sampling method changes (1/3, 1/2, 2/3), these four techniques still show a poor performance in all investigated techniques. In addition, the validation accuracy of almost all techniques is lower than modeling accuracy (Figure 8).
The modeling and validation results for GPS1/3, GPS1/2, and GPS2/3 are shown in Figure 9. The absolute difference of AGB − GRS between modeling and validation accuracy (we use ∇R2, ∇RMSE and ∇MAE to show the difference) appears in Table 7.
The results showed that all techniques perform differently in different modeling and validation sample sizes (Table 6 and Table 7, from 1/3 to 1/2 and 2/3). We observed a stable performance in MLR, PLSR and PCR (Figure 9 and Table 7) in the GPS sampling method. As modeling numbers increased (1/3, 1/2, 2/3), MLR and PLSR performed more stably [such as MLR: (∇R2: 0.07, 0.03, 0.03; ∇RMSE: 0.50, 0.35, 0.28; ∇MAE: 0.32, 0.30, 0.14) in Table 7)]; ANN, PLSR, RF and SVM with GPS sampling showed a clear growth in stability (Table 6 and Table 7). The above results indicate that the GPS sampling method is more effective than GRS for obtaining stable estimation models.

4.4. Estimation Accuracy with Leave One Sampling

The results for validation with LOS appear in Figure 10. PLSR provides the highest accuracy for AGB−LOS modeling. The LOS validation results shown in Figure 10 suggest that all of these techniques have impressive accuracy: the R2 values are at least 0.79 [PCR with RMSE = 1.63 t/ha, MAE = 1.24 t/ha, NRMSE = 0.10, CV(RMSE) = 0.25], and PLSR [R2 = 0.89, RMSE = 1.20 t/ha, MAE = 0.90 t/ha, NRMSE = 0.07, CV(RMSE) = 0.18] has the highest accuracy. The leave-one-out cross validation indicates that the prediction performance of these techniques can be ranked as (Figure 10i) PLSR > MLR > RF > SVM > BBRT > ANN > DT > PCR.

5. Analysis and Discussion

5.1. Analysis and Selection of Vegetation Indexes

Our correlation analysis shows that VIs are correlated with AGB to varying degrees. The results of the correlation analysis (Table 4) demonstrate that only the correlation of the top 20 VIs exceeds 0.2. Previous studies have shown that near-infrared- and red-band VIs are effective for estimating AGB [8,9,11]. The correlation between AGB and red- and near-infrared-band VIs is low in this study (Table 4). This may be because, during the reproductive stage, photosynthesis and the near-infrared reflectance [84] both clearly decrease (Figure 2), reducing the correlation between AGB and the red- and near-infrared-based VIs. This result is consistent with the results of a previous study [17]. By contrast, the correlation between the entire growth stage of winter wheat AGB estimates and the red-, green-, and blue-band spectral indexes is more promising (Table 4). Our study also demonstrated that these vegetation indices are effective in estimating AGB (Figure 6, top five VIs input). Thus, red-, green-, and blue-band spectral indexes are useful as they can be used to estimate AGB during vegetative growth and reproductive growth stages.
A serious multi-collinearity problem arises in the investigation of 54 VIs (Figure 5, Table A1 and Table A2); our results (Figure 6) show that the verification accuracy of eight techniques differs when fed multi-collinearity data (Figure 5) as the input. Previous studies have shown that machine learning techniques (ANN [81], BBRT [24] and RF [21,85]) can make full use of the narrow hyperspectral bands (strongly collinear data) and VIs. In the current study, ANN, BBRT and RF are almost unaffected by using multi-collinearity data (Figure 6a,b,g), which may indicate that these techniques are robust against noise, which may relate to the principles of these techniques. Garg et al. [86] indicated that the machine learning technique is suitable for tackling the multi-collinearity problem; our results showed that machine learning techniques have better abilities in tackling the multi-collinearity problem than that of conventional regression techniques (Figure 6). In addition, the results show that MLR has a poor performance when using multi-collinearity data to estimate AGB, which confirms the results of a previous study [87]; PLSR performs best in three conventional regression techniques for tackling the multi-collinearity problem, which confirms the results of a previous study [87]. Thus, PLSR is a useful tool that can be used to estimate several response variables simultaneously, while accounting for multi-collinearity variables [88].

5.2. Analysis of Noise Immunity

Our results show that machine learning techniques are more immune to powerful noise than conventional regression techniques (Figure 7); RF performs best in this noise test. This may be because the RF method randomly changes the input variable and validates the importance of the input data, thus generating a large number of decision trees and reducing the impact of noise; this result corresponds to the results of previous studies [25,89,90]. Our results (Figure 7) show that MLR is more sensitive to noise than PLSR, which is consistent with the findings of Zhao et al. [91]. The results of Atzberger et al. [7] indicate that the noise immunity of PCR, PLSR, and SMLR is ranked as PCR > PLSR > MLR, which is exactly the same ranking as obtained in the present work (Figure 7). The present noise immunity results (Figure 7) are important because repeated observations by remote-sensing techniques occur at different times; so techniques with poor noise immunity may lead to low accuracy because of data errors [26,78,79] (e.g., due to atmosphere, clouds, observation times, sensor noise). Our noise immunity results may explain why different regression studies of vegetation parameters based on remote sensing obtain significantly different results.

5.3. Analysis of Stability and Prediction Performance

PLSR and MLR both perform better in stability tests than machine learning techniques (Table 6 and Table 7, from 1/3 to 1/2 and 2/3 sampling). Farifteh et al. [92] indicated that PLSR performs more stably in soil salinity estimation than ANN (PLSR: R2: 0.6~0.98, RMSE% = 11.6~48%; ANN: R2 = 0.46~0.97, RMSE% = 12.5~57%). Thus, PLSR and MLR may be suitable for works in which fewer samples are available for modeling. BBRT and DT perform poorly in stability tests, and the AGB estimation model seems to be over fit because R2 is close to unity in all tests (Figure 6b,c). However, the validation accuracy of BBRT remains high, whereas that of DT is poor. Fewer studies evaluating DT for AGB estimation are available, which may be because DT does not deliver high accuracy for AGB estimation by remote sensing. Yuan et al. [81] indicated that the accuracy of the simple random sampling method is lower than stratified sampling, and our results are in agreement with that study; our results also indicate that all GPS models are more stable than GRS with 2/3 sampling (Table 6 GRS 2/3, and Table 7 GPS 2/3). This may be because the inappropriate sample selection method affects modeling and validation accuracy, which may indicate that GPS sampling is more suitable for these techniques. A previous study showed that stratified sampling helps to generate a good calibration set [82]; this may explain why GPS performed better than GRS in this study.
PLSR has the highest accuracy in leave-one-out cross validation (Figure 10a), while PCR has the lowest precision (Figure 10h). This comparison between PLSR and PCR is consistent with the results of Atzberger et al. [7], who estimated the aboveground-canopy Chl content (PLSR: R2 = 0.85, RMSE = 51; PCR: R2 = 0.57, RMSE = 82). In addition, the comparison herein between PLSR, ANN, and PCR is consistent with that of Mirzaie et al. [22] in their estimate of the water content of vegetation (PLSR: R2 = 0.93, RRMSE = 0.23; ANN: R2 = 0.83, RRMSE = 0.41; PCR: R2 = 0.78, RRMSE = 0.41). Thus, PLSR is a useful tool that can be used to estimate AGB with high accuracy.

6. Conclusions

We have provided herein a series of machine learning and conventional regression techniques to estimate hyperspectral winter wheat AGB and select input data for the sampling methods. We have also analyzed the noise immunity and prediction accuracy. The results allow the following conclusions to be drawn:
(1)
Machine learning is the correct technique for tackling the multi-collinearity problem. ANN, BBRT and RF are almost unaffected by the multi-collinearity problem (Figure 6a,b,g), while MLR and PCR could not solve it.
(2)
Machine learning techniques are much more immune to noise than conventional regression techniques. In terms of noise immunity, the techniques are ranked as follows (Figure 7): RF > SVM >DT > BBRT >ANN > PCR > PLSR > MLR. Thus, RF may be suitable for work that requires repeated observations via remote sensing.
(3)
The growth-period random sampling method performed better in stability tests. PLSR and MLR perform well in all stability tests (Figure 8 and Figure 9 and Table 6 and Table 7); these techniques and the sampling method may be suitable for work in which only a few samples are available for high-accuracy and stability estimation modeling.
(4)
This study demonstrated the potential application of VIs, machine learning and conventional regression techniques in estimating winter wheat biomass. The experimental results indicated that PLSR, MLR, and RF may be suitable for work that requires high-accuracy estimation models.

Acknowledgments

This study was supported by the National Key Research and Development Program (2016YFD0300603-5, 2016YFD0300602) the Natural Science Foundation of China (41601346, 41601369, 61661136003, 41771370, 41471285, 41471351), U.K. Science and Technology Facilities Council through the PAFiC project: Precision Agriculture for Family-farms in China (Ref.: ST/N006801/1), and the Special Funds for Technology innovation capacity building sponsored by the Beijing Academy of Agriculture and Forestry Sciences (KJCX20170423).

Author Contributions

Haikuan Feng, Zhenhai Li and Guijun Yang designed the experiments, Zhenhai Li, Jibo Yue and Haikuan Feng collected the AGB, Chl and ASD hyperspectral data. Jibo Yue and Haikuan Feng analyzed the data and wrote the manuscript, Guijun Yang and Zhenhai Li provided comments and suggestions for the manuscript and checked the writing.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Principal component analysis of the VIs was conducted by SPSS software (Statistical Product and Service Solutions, IBM, Amon, New York, NY, USA). Total explained variance is shown in Table A1. For each group of VIs, the different values represent the total explained variance with different components.
Table A1. Total explained variance of each group of VIs (%).
Table A1. Total explained variance of each group of VIs (%).
Component54VIs40VIs30VIs20VIs15VIs10VIs5VIs
175.89170.82265.58462.49565.40260.49468.135
285.00382.51380.51181.34887.25589.71189.136
392.21191.43591.55392.15794.39194.26997.232
495.22894.58694.75794.77496.54197.19999.927
597.35097.10997.17997.23598.05698.596100.000
698.14998.02598.22698.24699.03999.514-
798.70698.68398.84198.81799.48399.887-
899.21799.13799.24299.30699.77999.956-
999.42599.39199.55199.66299.88199.985-
1099.61499.63399.69699.82499.940100.000-
1199.72699.74699.80199.88999.965--
1299.81799.82999.87099.92799.986--
1399.87299.88599.90399.95899.994--
1499.90299.91499.93199.974100.000--
1599.92799.93399.95199.988100.000--
1699.94299.95099.96799.994---
1799.95699.96399.97999.998---
1899.96699.97599.985100.000---
1999.97499.98499.989100.000---
2099.98099.98999.994100.000---
2199.98599.99299.996----
2299.98899.99499.998----
2399.99199.99699.998----
2499.99399.99799.999----
2599.99499.998100.000----
2699.99699.999100.000----
2799.99699.999100.000----
2899.99799.999100.000----
2999.998100.000100.000----
3099.998100.000100.000----
3199.999100.000-----
3299.999100.000-----
3399.999100.000-----
3499.999100.000-----
35100.000100.000-----
36100.000100.000-----
37100.000100.000-----
38100.000100.000-----
39100.000100.000-----
40100.000100.000-----
41100.000------
42100.000------
43100.000------
44100.000------
45100.000------
46100.000------
47100.000------
48100.000------
49100.000------
50100.000------
51100.000------
52100.000------
53100.000------
54100.000------
Note: 54, 40, 30, 20, 15 and 5 represent each group of data; Components 1~54 represent the first to 54th component, respectively. Symbol “-” stand for “None”.
VIF (Table A2) provides an index that measures how much the variance of an estimated regression coefficient is increased because of collinearity. VI 1~VI 54 represent different VIs in Table 4 which fed into MLR modeling and validation.
Table A2. Variance inflation factor (VIF) of VIs.
Table A2. Variance inflation factor (VIF) of VIs.
VI54VIsVI40VIsVI30VIsVI20VIsVI15VIsVI10VIsVI5VIs
1121.3184.0146.6135.0120.7 113.717.1
2793.82400.92323.2228.327.226.023.3
31118.03615.03511.43317.03177.23143.83121.2
4492.34399.14346.34203.34173.3494.143.5
52673.352015.751492.15471.65258.95215.05150.6
616.2614.0613.9611.769.166.5--
71122.871526.171190.17395.07326.8719.9--
82441.381548.681290.881023.68779.88287.1--
9172.79114.8979.4941.8932.897.7--
112138.9111271.311828.011497.911438.610408.9--
121637.2121188.4121133.612805.612616.8----
13358.1131227.313956.01386.91329.4----
14175.414156.014148.314111.31472.1----
153285.7153332.0153899.8151486.715146.8----
16191.916101.71685.51639.2------
171130.6171128.617918.11723.0------
18156.218118.51879.01818.7------
216258.6213884.621886.0201280.1------
22663.822572.122470.4--------
24996.924706.224207.5--------
253304.2268495.2262382.6--------
285845.4272114.8277627.7--------
301939.9283010.1281554.9--------
328499.5301287.1302371.2--------
343063.2326092.332177.4--------
3759.03742.0--- ------
382705.938148.6----------
414203.8404259.2----------
454723.4------------
462086.8------------
504010.9------------
511151.1------------
524666.7------------
543874.4------------
Note: 54VIs, 40VIs, 30VIs, 20VIs, 15VIs and 5VIs represent different VIF values of each VI in four groups of data. Symbol “-” stand for “None”.

References

  1. Wang, J.; Zhao, C.; Huang, W. Fundamental and Application of Quantitative Remote Sensing in Agriculture; Science China Press: Beijing, China, 2008. [Google Scholar]
  2. Jin, X.; Li, Z.; Yang, G.; Yang, H.; Feng, H.; Xu, X.; Wang, J.; Li, X.; Luo, J. Winter wheat yield estimation based on multi-source medium resolution optical and radar imaging data and the AquaCrop model using the particle swarm optimization algorithm. ISPRS J. Photogramm. Remote Sens. 2017, 126, 24–37. [Google Scholar] [CrossRef]
  3. Huang, J.; Sedano, F.; Huang, Y.; Ma, H.; Li, X.; Liang, S.; Tian, L.; Zhang, X.; Fan, J.; Wu, W. Assimilating a synthetic Kalman filter leaf area index series into the WOFOST model to improve regional winter wheat yield estimation. Agric. For. Meteorol. 2016, 216, 188–202. [Google Scholar] [CrossRef]
  4. Hensgen, F.; Bühle, L.; Wachendorf, M. The effect of harvest, mulching and low-dose fertilization of liquid digestate on above ground biomass yield and diversity of lower mountain semi-natural grasslands. Agric. Ecosyst. Environ. 2016, 216, 283–292. [Google Scholar] [CrossRef]
  5. Jin, X.; Kumar, L.; Li, Z.; Xu, X.; Yang, G.; Wang, J. Estimation of winter wheat biomass and yield by combining the aquacrop model and field hyperspectral data. Remote Sens. 2016, 8. [Google Scholar] [CrossRef]
  6. Boschetti, M.; Bocchi, S.; Brivio, P.A. Assessment of pasture production in the Italian Alps using spectrometric and remote sensing information. Agric. Ecosyst. Environ. 2007, 118, 267–272. [Google Scholar] [CrossRef]
  7. Atzberger, C.; Guérif, M.; Baret, F.; Werner, W. Comparative analysis of three chemometric techniques for the spectroradiometric assessment of canopy chlorophyll content in winter wheat. Comput. Electron. Agric. 2010, 73, 165–173. [Google Scholar] [CrossRef]
  8. Fu, Y.; Yang, G.; Wang, J.; Song, X.; Feng, H. Winter wheat biomass estimation based on spectral indices, band depth analysis and partial least squares regression using hyperspectral measurements. Comput. Electron. Agric. 2014, 100, 51–59. [Google Scholar] [CrossRef]
  9. Yue, J.; Yang, G.; Li, C.; Li, Z.; Wang, Y.; Feng, H.; Xu, B. Estimation of Winter Wheat Above-Ground Biomass Using Unmanned Aerial Vehicle-Based Snapshot Hyperspectral Sensor and Crop Height Improved Models. Remote Sens. 2017, 9, 708. [Google Scholar] [CrossRef]
  10. Atzberger, C.; Darvishzadeh, R.; Immitzer, M.; Schlerf, M.; Skidmore, A.; le Maire, G. Comparative analysis of different retrieval methods for mapping grassland leaf area index using airborne imaging spectroscopy. Int. J. Appl. Earth Obs. Geoinf. 2015, 43, 19–31. [Google Scholar] [CrossRef]
  11. Cho, M.A.; Skidmore, A.; Corsi, F.; van Wieren, S.E.; Sobhan, I. Estimation of green grass/herb biomass from airborne hyperspectral imagery using spectral indices and partial least squares regression. Int. J. Appl. Earth Obs. Geoinf. 2007, 9, 414–424. [Google Scholar] [CrossRef]
  12. Schlerf, M.; Atzberger, C.; Hill, J. Remote sensing of forest biophysical variables using HyMap imaging spectrometer data. Remote Sens. Environ. 2005, 95, 177–194. [Google Scholar] [CrossRef]
  13. Galvão, L.S.; Formaggio, A.R.; Tisot, D.A. Discrimination of sugarcane varieties in Southeastern Brazil with EO-1 Hyperion data. Remote Sens. Environ. 2005, 94, 523–534. [Google Scholar] [CrossRef]
  14. Wang, W. The Associations of Photosynthesis and Gain Filling during Grain Filling Period in Flag Leaves of Wheat Species; Ocean University of China: Qingdao, China, 2007. [Google Scholar]
  15. Datt, B. A New Reflectance Index for Remote Sensing of Chlorophyll Content in Higher Plants: Tests using Eucalyptus Leaves. J. Plant Physiol. 1999, 154, 30–36. [Google Scholar] [CrossRef]
  16. Broge, N.H.; Leblanc, E. Comparing prediction power and stability of broadband and hyperspectral vegetation indices for estimation of green leaf area index and canopy chlorophyll density. Remote Sens. Environ. 2001, 76, 156–172. [Google Scholar] [CrossRef]
  17. Sun, H.; Li, M.Z.; Zhao, Y.; Zhang, Y.E.; Wang, X.M.; Li, X.H. The Spectral Characteristics and Chlorophyll Content at Winter Wheat Growth Stages. Spectrosc. Spectr. Anal. 2010, 30, 192–196. [Google Scholar] [CrossRef]
  18. Thenkabail, P.S. Biophysical and yield information for precision farming from near-real-time and historical Landsat TM images. Int. J. Remote Sens. 2003, 24, 2879–2904. [Google Scholar] [CrossRef]
  19. Rivera, J.; Verrelst, J.; Delegido, J.; Veroustraete, F.; Moreno, J. On the Semi-Automatic Retrieval of Biophysical Parameters Based on Spectral Index Optimization. Remote Sens. 2014, 6, 4927–4951. [Google Scholar] [CrossRef]
  20. Atzberger, C. Object-based retrieval of biophysical canopy variables using artificial neural nets and radiative transfer models. Remote Sens. Environ. 2004, 93, 53–67. [Google Scholar] [CrossRef]
  21. Yue, J.; Yang, G.; Feng, H. Comparative of remote sensing estimation models of winter wheat biomass based on random forest algorithm. Nongye Gongcheng Xuebao/Trans. Chin. Soc. Agric. Eng. 2016, 32, 175–182. [Google Scholar] [CrossRef]
  22. Mirzaie, M.; Darvishzadeh, R.; Shakiba, A.; Matkan, A.A.; Atzberger, C.; Skidmore, A. Comparative analysis of different uni- and multi-variate methods for estimation of vegetation water content using hyper-spectral measurements. Int. J. Appl. Earth Obs. Geoinf. 2014, 26, 1–11. [Google Scholar] [CrossRef]
  23. Du, H.; Sun, X.; Han, N.; Mao, F. RS estimation of inventory parameters and carbon storage of Moso bamboo forest based on synergistic use of object-based image analysis and decision tree. Chin. J. Appl. Ecol. 2017, 28. [Google Scholar] [CrossRef]
  24. Friedman, J.H. Greedy function approximation: A gradient boosting machine. Ann. Stat. 2001, 29, 1189–1232. [Google Scholar] [CrossRef]
  25. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef]
  26. Ali, I.; Greifeneder, F.; Stamenkovic, J.; Neumann, M.; Notarnicola, C. Review of machine learning approaches for biomass and soil moisture retrievals from remote sensing data. Remote Sens. 2015, 7, 16398–16421. [Google Scholar] [CrossRef]
  27. Cheng, J.H.; Dai, Q.; Sun, D.W.; Zeng, X.A.; Liu, D.; Pu, H. Bin Applications of non-destructive spectroscopic techniques for fish quality and safety evaluation and inspection. Trends Food Sci. Technol. 2013, 34, 18–31. [Google Scholar] [CrossRef]
  28. Pan, L.; Li, H.C.; Deng, Y.J.; Zhang, F.; Chen, X.D.; Du, Q. Hyperspectral dimensionality reduction by tensor sparse and low-rank graph-based discriminant analysis. Remote Sens. 2017, 9. [Google Scholar] [CrossRef]
  29. Gomez, C.; Lagacherie, P.; Coulouma, G. Continuum removal versus PLSR method for clay and calcium carbonate content estimation from laboratory and airborne hyperspectral measurements. Geoderma 2008, 148, 141–148. [Google Scholar] [CrossRef]
  30. Blackburn, G.A.; Pitman, J.I. Biophysical controls on the directional spectral reflectance properties of bracken (Pteridium aquilinum) canopies: Results of a field experiment. Int. J. Remote Sens. 1999, 20, 2265–2282. [Google Scholar] [CrossRef]
  31. Jin, X.; Xu, X.; Song, X.; Li, Z.; Wang, J.; Guo, W. Estimation of leaf water content in winter wheat using grey relational analysis-partial least squares modeling with hyperspectral data. Agron. J. 2013, 105, 1385–1392. [Google Scholar] [CrossRef]
  32. Pu, R.; Gong, P. Hyperspectral Remote Sensing of Vegetation Bioparameters. In Advances in Environmental Remote Sensing; CRC Press: Boca Raton, FL, USA, 2011; Volume 4, pp. 101–142. ISBN 9781420091816. [Google Scholar]
  33. Thenkabail, P.S.; Smith, R.B.; De Pauw, E. Hyperspectral vegetation indices and their relationships with agricultural crop characteristics. Remote Sens. Environ. 2000, 71, 158–182. [Google Scholar] [CrossRef]
  34. Liang, L.; Di, L.; Zhang, L.; Deng, M.; Qin, Z.; Zhao, S.; Lin, H. Estimation of crop LAI using hyperspectral vegetation indices and a hybrid inversion method. Remote Sens. Environ. 2015, 165, 123–134. [Google Scholar] [CrossRef]
  35. Gevrey, M.; Dimopoulos, I.; Lek, S. Review and comparison of methods to study the contribution of variables in artificial neural network models. Ecol. Model. 2003, 160, 249–264. [Google Scholar] [CrossRef]
  36. Cortes, C.; Vapnik, V. Support-Vector Networks. Mach. Learn. 1995, 20, 273–297. [Google Scholar] [CrossRef]
  37. Chang, C.-C.; Lin, C.-J. Libsvm. ACM Trans. Intell. Syst. Technol. 2011, 2, 1–27. [Google Scholar] [CrossRef]
  38. Utgoff, P.E.; Corporation, C.; Street, W. Decision Tree Induction Based on Efficient Tree Restructuring. Mach. Learn. 1997, 29, 5–44. [Google Scholar] [CrossRef]
  39. Fakiola, M.; Mishra, A.; Rai, M.; Singh, S.P.; O’Leary, R.A.; Ball, S.; Francis, R.W.; Firth, M.J.; Radford, B.T.; Miller, E.N.; et al. Classification and regression tree and spatial analyses reveal geographic heterogeneity in genome wide linkage study of indian visceral leishmaniasis. PLoS ONE 2010, 5, e15807. [Google Scholar] [CrossRef] [PubMed]
  40. Wold, H. Estimation of principal components and related models by iterative least squares. In Multivariate Analysis; Academic Press: New York, NY, USA, 1966; pp. 1391–1420. ISBN 0471411256. [Google Scholar]
  41. Baret, F.; Guyot, G. Potentials and limits of vegetation indices for LAI and APAR assessment. Remote Sens. Environ. 1991, 35, 161–173. [Google Scholar] [CrossRef]
  42. Sims, D.A.; Gamon, J.A. Relationships between leaf pigment content and spectral reflectance across a wide range of species, leaf structures and developmental stages. Remote Sens. Environ. 2002, 81, 337–354. [Google Scholar] [CrossRef]
  43. Huete, A.; Didan, K.; Miura, T.; Rodriguez, E.P.; Gao, X.; Ferreira, L.G. Overview of the radiometric and biophysical performance of the MODIS vegetation indices. Remote Sens. Environ. 2002, 83, 195–213. [Google Scholar] [CrossRef]
  44. Jiang, Z.; Huete, A.R.; Didan, K.; Miura, T. Development of a two-band enhanced vegetation index without a blue band. Remote Sens. Environ. 2008, 112, 3833–3845. [Google Scholar] [CrossRef]
  45. Zarco-Tejada, P.J.; Berjón, A.; López-Lozano, R.; Miller, J.R.; Martín, P.; Cachorro, V.; González, M.R.; De Frutos, A. Assessing vineyard condition with hyperspectral indices: Leaf and canopy reflectance simulation in a row-structured discontinuous canopy. Remote Sens. Environ. 2005, 99, 271–287. [Google Scholar] [CrossRef]
  46. Peñuelas, J.; Gamon, J.A.; Fredeen, A.L.; Merino, J.; Field, C.B. Reflectance indices associated with physiological changes in nitrogen- and water-limited sunflower leaves. Remote Sens. Environ. 1994, 48, 135–146. [Google Scholar] [CrossRef]
  47. Delalieux, S.; Somers, B.; Hereijgers, S.; Verstraeten, W.W.; Keulemans, W.; Coppin, P. A near-infrared narrow-waveband ratio to determine Leaf Area Index in orchards. Remote Sens. Environ. 2008, 112, 3762–3772. [Google Scholar] [CrossRef]
  48. Barnes, J.D.; Balaguer, L.; Manrique, E.; Elvira, S.; Davison, A.W. A reappraisal of the use of DMSO for the extraction and determination of chlorophylls a and b in lichens and higher plants. Environ. Exp. Bot. 1992, 32, 85–100. [Google Scholar] [CrossRef]
  49. Qi, J.; Chehbouni, A.; Huete, A.R.; Kerr, Y.H.; Sorooshian, S. A modified soil adjusted vegetation index. Remote Sens. Environ. 1994, 48, 119–126. [Google Scholar] [CrossRef]
  50. Rama Rao, N.; Garg, P.K.; Ghosh, S.K.; Dadhwal, V.K. Estimation of leaf total chlorophyll and nitrogen concentrations using hyperspectral satellite imagery. J. Agric. Sci. 2008, 146, 65–75. [Google Scholar] [CrossRef]
  51. Chen, J.M. Evaluation of vegetation indices and a modified simple ratio for boreal applications. Can. J. Remote Sens. 1996, 22, 229–242. [Google Scholar] [CrossRef]
  52. Gamon, J.A.; Peñuelas, J.; Field, C.B. A narrow-waveband spectral index that tracks diurnal changes in photosynthetic efficiency. Remote Sens. Environ. 1992, 41, 35–44. [Google Scholar] [CrossRef]
  53. Haboudane, D.; Miller, J.R.; Pattey, E.; Zarco-Tejada, P.J.; Strachan, I.B. Hyperspectral vegetation indices and novel algorithms for predicting green LAI of crop canopies: Modeling and validation in the context of precision agriculture. Remote Sens. Environ. 2004, 90, 337–352. [Google Scholar] [CrossRef]
  54. Blackburn, G.A. Quantifying chlorophylls and carotenoids at leaf and canopy scales: An evaluation of some hyperspectral approaches. Remote Sens. Environ. 1998, 66, 273–285. [Google Scholar] [CrossRef]
  55. Chappelle, E.W.; Kim, M.S.; McMurtrey, J.E. Ratio analysis of reflectance spectra (RARS): An algorithm for the remote estimation of the concentrations of chlorophyll A, chlorophyll B, and carotenoids in soybean leaves. Remote Sens. Environ. 1992, 39, 239–247. [Google Scholar] [CrossRef]
  56. Rouse, J.W.; Hass, R.H.; Schell, J.A.; Deering, D.W. Monitoring vegetation systems in the great plains with ERTS. In Proceedings of the Third Earth Resources Technology Satellite Symposium, Washington, DC, USA, 10–14 December 1973; pp. 309–317. [Google Scholar]
  57. Gamon, J.A.; Surfus, J.S. Assessing leaf pigment content and activity with a reflectometer. New Phytol. 1999, 143, 105–117. [Google Scholar] [CrossRef]
  58. Rondeaux, G.; Steven, M.; Baret, F. Optimization of soil-adjusted vegetation indices. Remote Sens. Environ. 1996, 55, 95–107. [Google Scholar] [CrossRef]
  59. Penuelas, J.; Baret, F.; Filella, I. Semi-empirical indices to assess carotenoids/chlorophyll a ratio from leaf spectral reflectance. Photosynthetica 1995, 31, 221–230. [Google Scholar]
  60. Nagler, P.L.; Daughtry, C.S.T.; Goward, S.N. Plant litter and soil reflectance. Remote Sens. Environ. 2000, 71, 207–215. [Google Scholar] [CrossRef]
  61. Roujean, J.L.; Breon, F.M. Estimating PAR absorbed by vegetation from bidirectional reflectance measurements. Remote Sens. Environ. 1995, 51, 375–384. [Google Scholar] [CrossRef]
  62. Serrano, L.; Peñuelas, J.; Ustin, S.L. Remote sensing of nitrogen and lignin in Mediterranean vegetation from AVIRIS data: Decomposing biochemical from structural signals. Remote Sens. Environ. 2002, 81, 355–364. [Google Scholar] [CrossRef]
  63. Vincini, M.; Frazzi, E.; D’Alessio, P. Angular Dependence of Maize and Sugar Beet VIs from Directional CHRIS/Proba Data. Available online: https://www.researchgate.net/profile/Ermes_Frazzi/publication/228413259_Angular_dependence_of_maize_and_sugar_beet_VIs_from_directional_CHRISProba_data/links/0046352d50c18b3fe6000000.pdf (accessed on 3 January 2018).
  64. Haboudane, D.; Miller, J.R.; Tremblay, N.; Zarco-Tejada, P.J.; Dextraze, L. Integrated narrow-band vegetation indices for prediction of crop chlorophyll content for application to precision agriculture. Remote Sens. Environ. 2002, 81, 416–426. [Google Scholar] [CrossRef]
  65. Jordan, C.F. Derivation of Leaf-Area Index from Quality of Light on the Forest Floor. Source Ecol. 1969, 50, 663–666. [Google Scholar] [CrossRef]
  66. Gitelson, A.A.; Kaufman, Y.J.; Stark, R.; Rundquist, D. Novel algorithms for remote estimation of vegetation fraction. Remote Sens. Environ. 2002, 80, 76–87. [Google Scholar] [CrossRef]
  67. Hunt, E.R.; Rock, B.N. Detection of changes in leaf water content using Near- and Middle-Infrared reflectances. Remote Sens. Environ. 1989, 30, 43–54. [Google Scholar] [CrossRef]
  68. Stamatiadis, S.; Taskos, D.; Tsadilas, C.; Christofides, C.; Tsadila, E.; Schepers, J.S. Relation of ground-sensor canopy reflectance to biomass production and grape color in two merlot vineyards. Am. J. Enol. Vitic. 2006, 57, 415–422. [Google Scholar] [CrossRef]
  69. Hardisky, M.A.; Klemas, V.; Smart, R.M. The influence of Soil Salinity, Growth Form, and Leaf Moisture on the Spectral Radiance of Spartina alterniflora Canopies. Photogramm. Eng. Remote Sens. 1983, 49, 77–83. [Google Scholar]
  70. Gitelson, A.A.; Merzlyak, M.N.; Chivkunova, O.B. Optical Properties and Nondestructive Estimation of Anthocyanin Content in Plant Leaves. Photochem. Photobiol. 2001, 74, 38–45. [Google Scholar] [CrossRef]
  71. Datt, B.; McVicar, T.R.; Van Niel, T.G.; Jupp, D.L.B.; Pearlman, J.S. Preprocessing EO-1 Hyperion hyperspectral data to support the application of agricultural indexes. IEEE Trans. Geosci. Remote Sens. 2003, 41, 1246–1259. [Google Scholar] [CrossRef]
  72. Fensholt, R.; Sandholt, I. Derivation of a shortwave infrared water stress index from MODIS near- and shortwave infrared data in a semiarid environment. Remote Sens. Environ. 2003, 87, 111–121. [Google Scholar] [CrossRef]
  73. Zarco-Tejada, P.J.; Rueda, C.A.; Ustin, S.L. Water content estimation in vegetation with MODIS reflectance data and model inversion methods. Remote Sens. Environ. 2003, 85, 109–124. [Google Scholar] [CrossRef]
  74. Daughtry, C.S.T.; Walthall, C.L.; Kim, M.S.; De Colstoun, E.B.; McMurtrey, J.E. Estimating corn leaf chlorophyll concentration from leaf and canopy reflectance. Remote Sens. Environ. 2000, 74, 229–239. [Google Scholar] [CrossRef]
  75. Penuelas, J.; Pinol, J.; Ogaya, R.; Filella, I. Estimation of plant water concentration by the reflectance Water Index WI (R900/R970). Int. J. Remote Sens. 1997, 18, 2869–2875. [Google Scholar] [CrossRef]
  76. Merzlyak, M.N.; Gitelson, A.A.; Chivkunova, O.B.; Rakitin, V.Y.U. Non-destructive optical detection of pigment changes during leaf senescence and fruit ripening. Physiol. Plant. 1999, 106, 135–141. [Google Scholar] [CrossRef]
  77. Merton, R.; Huntington, J. Early Simulation Results of the Aries-1 Satellite Sensor for Multi-Temporal Vegetation Research Derived from Aviris. In Proceedings of the Eighth Annual JPL Airborne Earth Science Workshop, Pasadena, CA, USA, 9–11 February 1999; pp. 1–10. [Google Scholar]
  78. Matsushita, Y.; Lin, S. Radiometric calibration from noise distributions. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Minneapolis, MN, USA, 17–22 June 2007. [Google Scholar] [CrossRef]
  79. Tsin, Y.; Ramesh, V.; Kanade, T. Statistical calibration of CCD imaging process. In Proceedings of the Eighth IEEE International Conference on Computer Vision, Vancouver, BC, Canada, 7–14 July 2001; pp. 480–487. [Google Scholar] [CrossRef]
  80. You, H.; Ma, Z.; Tang, Y.; Wang, Y.; Yan, J.; Ni, M.; Cen, K.; Huang, Q. Comparison of ANN (MLP), ANFIS, SVM, and RF models for the online classification of heating value of burning municipal solid waste in circulating fluidized bed incinerators. Waste Manag. 2016. [Google Scholar] [CrossRef] [PubMed]
  81. Yuan, H.; Yang, G.; Li, C.; Wang, Y.; Liu, J.; Yu, H.; Feng, H.; Xu, B.; Zhao, X.; Yang, X. Retrieving soybean leaf area index from unmanned aerial vehicle hyperspectral remote sensing: Analysis of RF, ANN, and SVM regression models. Remote Sens. 2017, 9. [Google Scholar] [CrossRef]
  82. Cochran, W.G. Sampling Technique; China Statistical Publishing House: Beijing, China, 1985. [Google Scholar]
  83. Volpe, V.; Manzoni, S.; Marani, M.; Katul, G. Leave-One-Out Cross-Validation; Springer: Berlin, Germany, 2011. [Google Scholar]
  84. Humbeck, K.; Quast, S.; Krupinska, K. Functional and molecular changes in the photosynthetic apparatus during senescence of flag leaves from field-grown barley plants. Plant Cell Environ. 1996, 19, 337–344. [Google Scholar] [CrossRef]
  85. Martínez-Muñoz, G.; Suárez, A. Out-of-bag estimation of the optimal sample size in bagging. Pattern Recognit. 2010, 43, 143–152. [Google Scholar] [CrossRef]
  86. Garg, A.; Kang, T.; Tai, K. Comparison of statistical and machine learning methods in modelling of data with multicollinearity. Int. J. Model. Identif. Control 2013, 18, 295–312. [Google Scholar] [CrossRef]
  87. Abdullah, S.; Ismail, M.; Fong, S.Y.; Ahmed, A.M.A.N. Evaluation for long term PM10 concentration forecasting using multi linear regression (MLR) and principal component regression (PCR) models. Environ. Asia 2016, 9, 101–110. [Google Scholar]
  88. Jin, X.; Ma, J.; Wen, Z.; Song, K. Estimation of maize residue cover using Landsat-8 OLI image spectral information and textural features. Remote Sens. 2015, 7, 14559–14575. [Google Scholar] [CrossRef]
  89. Altmann, A.; Toloşi, L.; Sander, O.; Lengauer, T. Permutation importance: A corrected feature importance measure. Bioinformatics 2010, 26, 1340–1347. [Google Scholar] [CrossRef] [PubMed]
  90. Deng, H.; Runger, G.; Tuv, E. Bias of importance measures for multi-valued attributes and solutions. In Proceedings of the 21st International Conference on Articial Neural Networks (ICANN 2011), Espoo, Finland, 14–17 June 2011; pp. 293–300. [Google Scholar]
  91. Zhao, W.; Hopke, P.K.; Qin, X.; Prather, K.A. Predicting bulk ambient aerosol compositions from ATOFMS data with ART-2a and multivariate analysis. Anal. Chim. Acta 2005, 549, 179–187. [Google Scholar] [CrossRef]
  92. Farifteh, J.; Van der Meer, F.; Atzberger, C.; Carranza, E.J.M. Quantitative analysis of salt-affected soil reflectance spectra: A comparison of two adaptive methods (PLSR and ANN). Remote Sens. Environ. 2007, 110, 59–78. [Google Scholar] [CrossRef]
Figure 1. (a) Location of study area shown in red. (b) Map showing Changping District in Beijing City. (c) Design of treatments and an unmanned-aerial-vehicle image of the experimental field (acquired on 12 May 2015). Three plant groups are present with two winter wheat varieties (J9843 and ZM175), three water treatments (W0, W1, and W2), and four nitrogen treatments (N0, N1, N2, and N3).
Figure 1. (a) Location of study area shown in red. (b) Map showing Changping District in Beijing City. (c) Design of treatments and an unmanned-aerial-vehicle image of the experimental field (acquired on 12 May 2015). Three plant groups are present with two winter wheat varieties (J9843 and ZM175), three water treatments (W0, W1, and W2), and four nitrogen treatments (N0, N1, N2, and N3).
Remotesensing 10 00066 g001
Figure 2. (a) Average hyperspectral reflectance spectrum for the four growing stages. (b) Average Chl and above-ground biomass (AGB) for the four growing stages.
Figure 2. (a) Average hyperspectral reflectance spectrum for the four growing stages. (b) Average Chl and above-ground biomass (AGB) for the four growing stages.
Remotesensing 10 00066 g002
Figure 3. Flowchart showing experiment methodology. Data selection, sampling methods, noise immunity, and prediction accuracy were analyzed.
Figure 3. Flowchart showing experiment methodology. Data selection, sampling methods, noise immunity, and prediction accuracy were analyzed.
Remotesensing 10 00066 g003
Figure 4. Remote-sensing imaging process and noise sources in a charge-coupled device (CCD) system.
Figure 4. Remote-sensing imaging process and noise sources in a charge-coupled device (CCD) system.
Remotesensing 10 00066 g004
Figure 5. Correlation coefficients r among 54 VIs. NPQI (first), BGI (second), ARI (16th) and TCARI (22nd) are weakly correlated with other VIs. Zone a (orange and light blue), low correlation zone; Zone b (red and dark blue), high correlation zone. VIs are ordered according to corr. coeff to AGB.
Figure 5. Correlation coefficients r among 54 VIs. NPQI (first), BGI (second), ARI (16th) and TCARI (22nd) are weakly correlated with other VIs. Zone a (orange and light blue), low correlation zone; Zone b (red and dark blue), high correlation zone. VIs are ordered according to corr. coeff to AGB.
Remotesensing 10 00066 g005
Figure 6. Best AGB modeling and validation accuracy for eight techniques with different input VIs: (a) ANN, (b) BBRT, (c) DT, (d) MLR, (e) PLSR, (f) PCR, (g) RF, (h) SVM. Note: R2 (M) and R2 (V) indicate R2 for modeling and validation, respectively. The same notation is used for MAE (M) and MAE (V), RMSE (M), and RMSE (V).
Figure 6. Best AGB modeling and validation accuracy for eight techniques with different input VIs: (a) ANN, (b) BBRT, (c) DT, (d) MLR, (e) PLSR, (f) PCR, (g) RF, (h) SVM. Note: R2 (M) and R2 (V) indicate R2 for modeling and validation, respectively. The same notation is used for MAE (M) and MAE (V), RMSE (M), and RMSE (V).
Remotesensing 10 00066 g006
Figure 7. Three measures of accuracy as a function of the signal-to-noise ratio. Note: SNR 100 represents no noise. All techniques used the first 30 sets of VIs (Table 4) as the input, with SNR = 5, 10, 15, 20, 30, 50, and no noise.
Figure 7. Three measures of accuracy as a function of the signal-to-noise ratio. Note: SNR 100 represents no noise. All techniques used the first 30 sets of VIs (Table 4) as the input, with SNR = 5, 10, 15, 20, 30, 50, and no noise.
Remotesensing 10 00066 g007
Figure 8. Modeling and validation results for global random sampling (GRS).
Figure 8. Modeling and validation results for global random sampling (GRS).
Remotesensing 10 00066 g008
Figure 9. Modeling and validation results for growth period sampling (GPS).
Figure 9. Modeling and validation results for growth period sampling (GPS).
Remotesensing 10 00066 g009
Figure 10. Measured and estimated AGB using leave one sampling. (a): PLSR; (b): MLR; (c): RF; (d): SVM; (e): BBRT; (f): ANN; (g): DT; (h): PCR; (i): R2, RMSE, MAE, NRMSE and CV(RMSE) of all techniques.
Figure 10. Measured and estimated AGB using leave one sampling. (a): PLSR; (b): MLR; (c): RF; (d): SVM; (e): BBRT; (f): ANN; (g): DT; (h): PCR; (i): R2, RMSE, MAE, NRMSE and CV(RMSE) of all techniques.
Remotesensing 10 00066 g010
Table 1. Statistics of AGB measurement in study area.
Table 1. Statistics of AGB measurement in study area.
PeriodSampleMin (t/ha)Max (t/ha)Mean (t/ha)Standard Deviation (t/ha)Coefficient of Variation (%)
Jointing481.2014.5262.5690.68526.664
Flag482.1948.2665.1141.46828.706
Flowering483.41912.7377.7901.96025.160
Grain filling485.45617.59910.9932.79325.407
Table 2. Summary of VIs used in this study.
Table 2. Summary of VIs used in this study.
VIsEquationReferenceVIsEquationReference
ATSAVIa (R800 − a R670 − b)/[(a R800 + R670 − ab + X(1 + a2)], where X = 0.08, a = 1.22, and b = 0.03[41]MND680(R800R680)/(R800 + R680 − 2R445)[42]
EVI2.5(RNIRRRed)/(RNIR + 6RRed − 7.5RBlue + 1)[43]MND705(R750R705)/(R750 + R705 − 2R445)[42]
EVI22.5(RNIRRRed)/(RNIR + 2.4RRed + 1)[44]MSR705(R750R445)/(R705R445)[42]
GIR554/R677[45]NPCI(R680R430)/(R680 + R430)[46]
LAIDIR1250/R1050[47]NPQI(R415R435)/(R415 + R435)[48]
MSAVI0.5[2R800 + 1 − ((2R800 + 1)2 − 8(R800R670))1/2][49]PBIR810/R560[50]
MSR(R800/R670 − 1)/(R800/R670 + 1)1/2[51]PRI(R531R570)/(R531 + R570)[52]
MTVI11.2[1.2(R800R550) − 2.5(R670R550)][53]PSSRR800/R500[54]
MTVI2{1.5[1.2(R800R550) − 2.5(R670R550)]}/
{(2R800 + 1)2 − [6R800 − 5(R670)1/2] − 0.5}1/2
[53]RARSR760/R500[55]
NDVI(RNIRRRed)/(RNIR + RRed)[56]RGRRRed/RGreen[57]
OSAVI1.16(R800R670)/(R800 + R670 + 0.16)[58]SIPI(R800R445)/(R800R680)[59]
PSND(R800R470)/(R800 + R470)[54]TVI0.5[120(R750R550) − 200(R670R550)][16]
PVIhyp(R1148 – a R807 − b)/(1 + a2)1/2,
where a = 1.17 and b = 3.37
[12]CAI0.5(R2020 + R2220) − R2100[60]
RDVI(R800R670)/(R800 + R670)1/2[61]NDLI[log(1/R1754) − log(1/R1680)]
/[log(1/R1754) + log(1/R1680)]
[62]
SLAIDIS(R1050R1250)/(R1050 + R1250), where S = 5[47]NDNI[log(1/R1510) − log(1/R1680)]
/[log(1/R1510) + log(1/R1680)]
[62]
SPVI0.4[3.7(R800R670) − 1.2|R530R670|][63]DSWI(R802 + R547)/(R1657 + R682)[13]
TCARI3[(R700R670) − 0.2(R700R550)(R700/R670)][64]LWVI1(R1094R893)/(R1094 + R893)[13]
SRRNIR /RRed[65]LWVI2(R1094R1205)/(R1094 + R1205)[13]
VARIgreen(RGreenRRed)/(RGreen + RRed)[66]MSIR1600/R819[67]
WDRVI(0.1 RNIRRRed)/(0.1 RNIR + RRed)[68]NDII(R819R1600)/(R819 + R1600)[69]
ARI(R550)1 − (R700)1[70]NDWI(R860R1240)/(R860 + R1240)[71]
BGIR450/R550[45]RVIhypR1088/R1148[12]
BRIR450/R690[45]SIWSI(R860R1640)/(R860 + R1640)[72]
LCI(R850R710)/(R850 + R680)[15]SRWIR860/R1240[73]
MCARI[(R701R671) − 0.2(R701R549)]/(R701/R671)[74]WIR900/R970[75]
MCARI11.2[2.5(R800R670) − 1.3(R800R550)][53]PSRI(R680R500)/R750[76]
MCARI2{1.5[2.5(R800R670) − 1.3(R800R550)]}
/{(2R800 + 1)2 − [6R800 − 5(R670)1/2] − 0.5}1/2
[53]RVSI[(R712 + R752)/2] − R732[77]
Note: RGreen, RRed, and RNIR represent bands at 470, 670, and 800 nm of hyperspectral reflectance, respectively. R470 and R800 represent bands at 470 and 800 nm of hyperspectral reflectance, etc.
Table 3. Maximum, minimum, and step length (SL) of ANN, RF, and SVM parameters.
Table 3. Maximum, minimum, and step length (SL) of ANN, RF, and SVM parameters.
MethodANNRFSVM
ParametersHidden layer 1Hidden layer 2SLntreeSLmtrySLcSLgSL
Min value11102011−100.5−100.5
Max value20202000101010
Note: Only parameters to be optimized appear in this table; other parameters were determined as per Refs. [34,80,81]. Number of modeling and verification: ANN, 20 × 20 × 7 = 2800; RF: 100 × 10 × 7 = 7000; SVM: 40 × 40 × 7 = 11,200.
Table 4. Correlation coefficients r between AGB and VIs (n = 192).
Table 4. Correlation coefficients r between AGB and VIs (n = 192).
NumberVIrNumberVIrNumberVIr
1NPQI0.757 **19MCARI10.226 **37NDLI0.106 n.s.
2BGI0.555 **20MTVI10.226 **38PSSR0.102 n.s.
3BRI0.519 **21MND6800.193 **39ATSAVI0.092 n.s.
4RVIhyp0.490 **22TCARI0.188 **40OSAVI0.088 n.s.
5NPCI0.474 **23EVI20.182 *41MND7050.086 n.s.
6CAI0.442 **24PSND0.180 *42RARS0.079 n.s.
7PVIhyp0.370 **25MSAVI0.180 *43SIWSI0.075 n.s.
8LWVI20.352 **26TVI0.177 *44NDWI0.075 n.s.
9LWVI10.337 **27GI0.176 *45PBI0.067 n.s.
10SLAIDI0.304 **28PRI0.176 *46MSR7050.061 n.s.
11SRWI0.300 **29VARIgreen0.165 *47NDVI0.039 n.s.
12LAIDI0.397 **30SIPI0.157 *48LCI0.034 n.s.
13RVSI0.299 **31MCARI20.149 *49MSR0.032 n.s.
14WI0.260 **32PSRI0.149 *50WDRVI0.031 n.s.
15SPVI0.251 **33RDVI0.146 *51SR0.029 n.s.
16ARI0.241 **34RGR0.142 *52MSI0.029 n.s.
17MCARI0.237 **35EVI0.149 *53DSWI0.027 n.s.
18NDNI0.228 **36MTVI20.121 *54NDII0.012 n.s.
Note: Probability levels are indicated by n.s., *, and ** for “not significant” (up to 0.119), 0.05 (greater than 0.141), and 0.01 (greater than 0.185), respectively. ** r (0.01, 192) = 0.185; * r (0.05, 192) = 0.141.
Table 5. Optimal number of input VIs for the eight analytical techniques.
Table 5. Optimal number of input VIs for the eight analytical techniques.
TechniqueANNBBRTDTMLRPLSRPCRRFSVM
Number of input VIs 3055201553030
Note: the number of input VIs (n) represents the top n VIs in Table 4.
Table 6. Absolute difference of AGB − GRS modeling and accuracy of validation.
Table 6. Absolute difference of AGB − GRS modeling and accuracy of validation.
TechniqueR2∇RMSE (t/ha)∇MAE (t/ha)
1/31/22/31/31/22/31/31/22/3
ANN0.17 *0.110.19 *0.450.590.650.250.690.52
MLR0.060.010.080.210.020.100.110.030.10
DT0.18*0.18*0.31 *1.12 *0.83 *0.660.560.700.91 *
BBRT0.23 *0.19*0.31 *1.64 *1.49 *1.68 *1.23 *1.17 *1.23 *
PLSR0.050.010.100.190.010.200.100.020.17
RF0.16*0.140.24 *0.85 *0.85*1.00 *0.640.660.76
SVM0.020.090.19 *0.140.300.430.060.380.41
PCR0.060.040.17 *0.070.140.260.100.110.26
Note: ∇R2 > 0.15, ∇RMSE and ∇MAE > 0.800 t/ha are marked with *. ∇R2, ∇RMSE, and ∇MAE show the absolute difference between modeling and validation R2, RMSE, MAE, respectively. The smaller the difference, the more stable and reliable the technique.
Table 7. Absolute difference AGB − GPS between modeling and validation accuracy.
Table 7. Absolute difference AGB − GPS between modeling and validation accuracy.
TechniqueR2∇RMSE (t/ha)∇MAE (t/ha)
1/31/22/31/31/22/31/31/22/3
ANN0.22 *0.070.031.20 *0.570.341.00 *0.600.26
MLR0.070.030.030.500.350.280.320.300.14
DT0.140.17 *0.17 *1.37 *1.20 *0.500.781.00 *0.77
BBRT0.23 *0.19 *0.20 *1.83 *1.64 *1.56 *1.36 *1.21 *1.23 *
PLSR0.040.040.000.230.460.070.220.400.02
RF0.140.120.111.04 *1.03 *0.87 *0.80 *0.770.70
SVM0.16 *0.070.050.81 *0.720.150.620.760.37
PCR0.010.060.010.110.470.060.100.360.01
Note: ∇R2 > 0.15, ∇RMSE and ∇MAE > 0.800 t/ha are marked with *. ∇R2, ∇RMSE, and ∇MAE show the absolute difference between modeling and validation R2, RMSE, MAE, respectively. The smaller the difference, the more stable and reliable the technique.

Share and Cite

MDPI and ACS Style

Yue, J.; Feng, H.; Yang, G.; Li, Z. A Comparison of Regression Techniques for Estimation of Above-Ground Winter Wheat Biomass Using Near-Surface Spectroscopy. Remote Sens. 2018, 10, 66. https://doi.org/10.3390/rs10010066

AMA Style

Yue J, Feng H, Yang G, Li Z. A Comparison of Regression Techniques for Estimation of Above-Ground Winter Wheat Biomass Using Near-Surface Spectroscopy. Remote Sensing. 2018; 10(1):66. https://doi.org/10.3390/rs10010066

Chicago/Turabian Style

Yue, Jibo, Haikuan Feng, Guijun Yang, and Zhenhai Li. 2018. "A Comparison of Regression Techniques for Estimation of Above-Ground Winter Wheat Biomass Using Near-Surface Spectroscopy" Remote Sensing 10, no. 1: 66. https://doi.org/10.3390/rs10010066

APA Style

Yue, J., Feng, H., Yang, G., & Li, Z. (2018). A Comparison of Regression Techniques for Estimation of Above-Ground Winter Wheat Biomass Using Near-Surface Spectroscopy. Remote Sensing, 10(1), 66. https://doi.org/10.3390/rs10010066

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