[go: up one dir, main page]
Next Article in Journal
Promising Uses of the iPad Pro Point Clouds: The Case of the Trunk Flare Diameter Estimation in the Urban Forest
Previous Article in Journal
Monitoring Asian Dust Storms from NOAA-20 CrIS Double CO2 Band Observations
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Remote Sensing
Volume 14
Issue 18
10.3390/rs14184660
remotesensing-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Winter Wheat SPAD Value Inversion Based on Multiple Pretreatment Methods

by
Lanzhi Shen
1,2,3,
Maofang Gao
1,*,
Jingwen Yan
3,
Qizhi Wang
1 and
Hua Shen
4
1
Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences/Key Laboratory of Agricultural Remote Sensing, Ministry of Agriculture and Rural Affairs, Beijing 100081, China
2
21st Century Space Technology Application Co., Ltd., Beijing 100096, China
3
Key Laboratory of Digital Signal and Image Processing of Guangdong Province, Shantou University, Shantou 515063, China
4
School of Emergency Management, Institute of Disaster Prevention, Sanhe 065201, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(18), 4660; https://doi.org/10.3390/rs14184660
Submission received: 6 August 2022 / Revised: 11 September 2022 / Accepted: 14 September 2022 / Published: 18 September 2022
Browse Figures
Graphical abstract
"> Figure 1
<p>Overview of the study area.</p> "> Figure 2
<p>The association between spectral indices and SPAD value at different growth stages.</p> "> Figure 3
<p>The estimated results of SPAD value for (<b>a</b>) training set and (<b>b</b>) validation set at jointing stage.</p> "> Figure 4
<p>The estimated results of SPAD value for (<b>a</b>) training set and (<b>b</b>) validation set at heading stage.</p> "> Figure 5
<p>The estimated results of SPAD value for (<b>a</b>) training set and (<b>b</b>) validation set at flowering stage.</p> "> Figure 6
<p>The estimated results of SPAD value for (<b>a</b>) training set and (<b>b</b>) validation set at filling stage.</p> "> Figure 7
<p>The estimated results of SPAD value for (<b>a</b>) training set and (<b>b</b>) validation set at milk stage.</p> "> Figure 8
<p>The estimated results of SPAD value for (<b>a</b>) training set and (<b>b</b>) validation set at mature stage.</p> "> Figure 9
<p>The estimated results of SPAD value for (<b>a</b>) training set and (<b>b</b>) validation set at whole growth stage.</p> "> Figure 10
<p>The estimated results of SPAD value for (<b>a</b>) training set and (<b>b</b>) validation set based on pretreatment WPD-(1/R)’-PCADR at jointing stage.</p> "> Figure 11
<p>The estimated results of SPAD value for (<b>a</b>) training set and (<b>b</b>) validation set based on pretreatment WPD-1/R-PCADR at heading stage.</p> "> Figure 12
<p>The estimated results of SPAD value for (<b>a</b>) training set and (<b>b</b>) validation set based on pretreatment WPD-(1/R)’-PCADR at flowering stage.</p> "> Figure 13
<p>The estimated results of SPAD value for (<b>a</b>) training set and (<b>b</b>) validation set based on pretreatment ND-(1/R)’-PCADR at filling stage.</p> "> Figure 14
<p>The estimated results of SPAD value for (<b>a</b>) training set and (<b>b</b>) validation set based on pretreatment WPD-(1/R)’-PCADR at milk stage.</p> "> Figure 15
<p>The estimated results of SPAD value for (<b>a</b>) training set and (<b>b</b>) validation set based on pretreatment WPD-R’-PCADR at mature stage.</p> "> Figure 16
<p>The estimated results of SPAD value for (<b>a</b>) training set and (<b>b</b>) validation set based on pretreatment WPD-R’-PCADR at whole growth stage.</p> "> Figure 17
<p>SPAD value distribution.</p> "> Figure 18
<p>The correlation diagram of model RPD and STD of sample SPAD values.</p> "> Figure 19
<p>The correlation diagram of model R<sup>2</sup><sub>V</sub> and standard deviation of sample SPAD values.</p> ">
Versions Notes

Abstract

:
SPAD value was measured by a portable chlorophyll instrument, which can reflect the relative chlorophyll content of vegetation well. Chlorophyll is an important organic chemical substance in plants that acquires and transmits energy during photosynthesis. The continuous spectral curve of winter wheat can be obtained rapidly in a specific band range by using hyperspectral remote sensing technology to estimate the SPAD value of winter wheat, which is of great significance to the growth monitoring and yield estimation research of winter wheat. In this study, with winter wheat as the research object, the spectral data and corresponding SPAD value in different growth stages were used as the data source, 20 kinds of data preprocessing spectra and sensitive spectral indices set the data as model input values, the partial least square regression (PLSR) model was established to estimate the SPAD value, and the model estimation results of different model input values at different growth stages were compared in detail. The results showed that the set of sensitive spectral indices selected in this study as input values can effectively improve the accuracy and stability of the PLSR model. In addition, the effects of 20 spectral data pretreatment methods on the estimation results of the SPAD value were compared and analyzed in different growth stages. It was found that the spectral data pretreated by the combination of wavelet packet denoising, first-order derivative transformation and principal component analysis can improve the accuracy and stability of PLSR model, and it is suitable for all growth stages. The results also showed that the estimation model is highly sensitive to the standard deviation of the SPAD value (STDchl) in sample sets. When the standard deviation is greater than 5.5 SPAD, the larger the STDchl is, the higher the model estimation accuracy is, and the more stable the model is. At this time, the model estimation accuracy is higher (R2V is greater than 0.5, ratio of performance to deviation is greater than 1.4), which can meet the estimation requirements of the SPAD value.

Graphical Abstract">

Graphical Abstract

1. Introduction

Chlorophyll is the main pigment in plants that aids in photosynthesis and plays a great role in the acquisition and transmission of energy in plants [1,2]. As an important index to evaluate the healthy growth of vegetation, it is important to obtain the chlorophyll content quickly and nondestructively [3,4]. A portable chlorophyll instrument is used to measure a relative chlorophyll content, which has real-time, rapid, and non-destructive advantages [5].
At present, hyperspectral remote sensing technology has been widely studied and applied in agriculture, including spectral feature analysis of crop leaves or canopy, monitoring and inversion of some physiological and biochemical parameters in crops, crop identification and classification, pest monitoring, etc. [6,7,8]. This technology can quickly obtain the continuous spectral curve of crops, and its small difference can reflect different characteristics such as crop growth. The visible range band (400–780 nm) of the spectrum is the strong absorption band of pigment [9]. The change of visible spectrum waveform is closely related to the content of chlorophyll in leaves. Jin et al. analyzed the spectrum of rice leaves and the corresponding chlorophyll content and found that the sensitive bands related to the chlorophyll content of rice leaves were 450–686 and 750–780 nm [10]. Lu et al. analyzed the correlation between the spectrum of cherry leaves and chlorophyll content after making various changes [11]. Using stepwise linear regression analysis, it was found that the variable log(1/R741) was highly correlated with the chlorophyll of cherry leaves. Zhang et al. conducted a variety of typical changes in the original winter wheat spectrum and found that the spectral change data based on the second derivative R′′ can accurately estimate the chlorophyll content under low temperature stress [12]. In the range of red light, the research of red edge parameters is also the focus of many researchers. Bonham Carter et al. defined the red edge position as the wavelength value corresponding to the first-order differential maximum in the red band and then analyzed its relationship with leaf chlorophyll content in detail [13]. Yang Jie et al. compared the reflectance of any two bands in the spectrum of rice leaves from 350 to 2500 nm and calculated its normalization index [14]. They proposed that the band ratio SR(R724, R709) and normalization index ND(R780, R709) are the two spectral indices most closely related to the chlorophyll content of rice leaves, and the determination coefficients are above 0.9. To extract sensitive bands, Wang et al. employed the elastic net approach to minimize the dimensionality of hyperspectral data and then retrieved the chlorophyll content of winter wheat [15]. Qiu et al. applied the SPAD value and spectral techniques to study the change of nitrogen during rape growth [16].
Because the spectral data obtained by the instrument has a certain noise, most researchers will preprocess the spectral data in the estimation of chlorophyll content of winter wheat [10,11,12,13,14]. However, single or two or three kinds of preprocessed spectral data or some spectral indices are generally used as a model input to estimate winter wheat chlorophyll content, and different preprocessing methods are used in different studies [11,12,14]. Few studies compare the effects of these common preprocessing or spectral indices on the results of a winter wheat leaf green content estimation model. Therefore, it is necessary to select some common preprocessed spectral data and select some spectral indices with high correlation with chlorophyll content as input values, and compare and analyze the model accuracy under various input values. In addition, for the remote sensing monitoring of winter wheat chlorophyll content, most researchers perform research based on the sample data of a certain growth period [17], which lacks universality for all growth periods. Therefore, it is necessary to consider the whole growth period of winter wheat and find a chlorophyll estimation model suitable for the whole growth period.
Therefore, this study will collect the spectral data of Winter Wheat in multiple growth stages (from jointing stage to maturity stage), take different combinations of preprocessed spectral data and spectral indices data [18] as input values, respectively, and model to estimate the SPAD value. In addition, PLSR [19,20] will be used to compare and analyze the estimation results of different models and explore the influencing factors of model accuracy and stability. It has certain application value for estimating the SPAD value in different growth stages and provides a theoretical basis and data preprocessing reference for estimating the chlorophyll content of vegetation from spaceborne hyperspectral satellite images.

2. Experiments and Methods

2.1. Study Area

The study area is located in the international agricultural high tech Industrial Park of Chinese Academy of Agricultural Sciences in Langfang City, Hebei Province. The geographical coordinates are 116°35′34″E–116°35′36″E and 39°35′50″N–39°35′53″N (Figure 1). The average annual temperature is 13.9 °C, and the average annual precipitation is 393.4 mm, belonging to the warm temperate continental monsoon climate. The total area of the test area is 1200 m2 (20 m wide and 60 m long). The soil texture is mainly loamy, and the soil types are mainly brown soil, moisture soil and cinnamon soil, with a unit weight of 1.5 g/cm3. The winter wheat variety selected in the experiment is Jingdong 12, and the fertilization level is consistent with the field fertilization level in this area (pure nitrogen (N) 256.5 kg·hm−2, phosphorus pentoxide (P2O5) 240 kg·hm−2, all input at one time during sowing).

2.2. Spectral and SPAD Value Data Acquisition

The data used in the study mainly include the spectral data of winter wheat leaves at various growth stages and the corresponding SPAD value.
The spectral data of winter wheat leaves were measured by portable ground object spectrometer ASD FieldSpec4. The wavelength range of the spectrometer is 350–2500 nm, the spectral sampling interval is 1.4 nm (350–1000 nm), 2 nm (1001–2500 nm) and the spectral resolution is 3 nm (700 nm), 8 nm (1400, 2100 nm). The spectra of winter wheat leaves were measured from 10:30 to 13:30 Beijing time when the weather was sunny, cloudless or light cumulus, windless or wind speed lower than 1 m/s. Before measurement, the optical fiber probe is vertically aligned with the standard reference plate (whiteboard) for calibration. The winter wheat spectrum measurement is re-calibrated with the whiteboard every 15 min to prevent the influence caused by the drift of the sensor response system and the change of the solar incidence angle [21]. During the measurement, the measurer was dressed in dark color, facing the vertical direction of the sun (perpendicular to the main plane direction), holding a pistol, and the optical fiber probe was vertically pointing downward at the no shadow in the front leaf of the top leaf of winter wheat. The spectral measurement position of winter wheat was selected at the place with less stem veins and more mesophyll. The vertical distance between the optical fiber probe and the blade is 1 cm. The spectral curves of each leaf sample are measured repeatedly, and 5 spectral curves are collected, and the average value of the 5 spectral curves is taken as the final spectral reflectance of the sample. The spectrum of green vegetation is mainly affected by pigment content in the range of 400–900 nm, and it will show obvious peak and valley characteristics, which is a representative band range distinguished from non-green vegetation [22]. The study only makes correlation analysis on the spectral data in the range of 400–900 nm. The SPAD value was measured by portable chlorophyll meter SPAD-502 and replaced by a SPAD value. After the leaf spectrum is measured, a portable chlorophyll meter is used to evenly measure five values at the spectral measurement points immediately and the average value is calculated as the SPAD value sample. Each growth stage was randomly sampled in the experimental area.
Spectral data and SPAD data were collected at the jointing stage, heading stage, flowering stage, grain filling stage, milk ripening stage and mature stage. The data of the whole growth period include the data of six single-growth periods. The set aside method is used to divide the training set and the verification set, in which 4/5 of the samples are used for training, and the remaining 1/5 are reserved for verification. In order to make the types of samples in the training set and the verification set as balanced as possible, the samples in each growth period are sorted from small to large according to the SPAD value, and then, one sample from every four is taken into the verification set according to the approximate ratio of 4:1 between the training set and the verification set, and the remaining samples are put into the training set. The data set size and SPAD value information of each growth period are shown in Table 1.

2.3. Spectral Index Collection

Affected by pigments, plant leaves have strong absorption characteristics in visible light band, especially in red band, and strong reflection characteristics in near-infrared band, which is the physical basis of remote sensing monitoring of green plants [23]. Through the reflectivity value of each band of visible light and near infrared, different vegetation indices can be obtained by different combinations of derivative values. 62 spectral indices related to the inversion of plant chlorophyll content in the published literature were collected and sorted out (Table 2). By analyzing the correlation between these spectral indices and the SPAD value, the sensitive spectral indices set will be selected to predict the SPAD value of winter wheat.

2.4. Spectral Preprocessing Methods

Denoising, data form transformation and dimension reduction single preprocessing were applied to the original winter wheat spectral data in turn (Table 3). Two different denoising approaches are no denoising (ND) and wavelet packet denoising (WPD) [63,64,65]. Five different mathematical forms of transformation are the original spectral reflectance data R, reciprocal 1/R, log(R), first derivative R’ and reciprocal first derivative (1/R)’ [66]. Two different dimensionality reduction approaches are no dimension reduction (NDR) and principle component analysis dimension reduction (PCADR) [67,68,69,70].
Denoising, data form transformation, and dimension reduction single preprocessing are combined to form a total of 20 combined preprocessing methods. The model’s input value is various combinations of preprocessed spectral data, and the model estimation is established by PLSR.

2.5. Partial Least Squares Regression

The link between vegetation spectrum and vegetation chlorophyll concentration can be studied using partial least squares regression (PLSR). This method can estimate the dependent variable through the linear combination of independent variables [71]. It is a multiple independent variables to multiple dependent variables linear regression modeling technique. When the sample size is small and there are multiple correlations between variables, the PLSR method has more advantages than the traditional regression model [71]. Due to the possible strong correlation between independent variables, it is easy to lead to overfitting of the problem when solving the linear regression problem of multiple independent variables to multiple dependent variables. However, the PLSR technique will identify several linear independent new variables to substitute the previous independent variables, allowing the diversity between independent variables to be maximized and the overfitting problem to be alleviated.
In this paper, the preprocessed spectral data and the sensitive vegetation indices with high correlation with winter wheat chlorophyll are used as the input values of the model, and the PLSR method is used to estimate and analyze the SPAD value of winter wheat.

2.6. Model Performance Evaluation Indices

The root mean square error of the training set (RMSET), the root mean square error of the validation set (RMSEV), the determination coefficient of the training set R2T, the determination coefficient of the verification set R2V, and the ratio of performance to deviation (RPD) were all used to carefully consider the model’s effectiveness in the research.
The determination coefficient R2 is the second power of the correlation coefficient R. The calculation formula of the correlation coefficient R is shown in Formula (1).
R i = cov ( x , y ) D ( x ) D ( y ) = ∑ n = 1 N ( x n i − x i ¯ ) ( y n − y ¯ ) ∑ n = 1 N ( x n i − x i ¯ ) 2 ∑ n = 1 N ( y n − y ¯ ) 2
In Formula (1), the Ri is the correlation coefficient between the winter wheat leaf spectral reflectance in band i and the winter wheat leaf SPAD value, N is the total number of samples, xni is the spectral reflectance of winter wheat leaves in band i of the nth sample, x i ¯ is the average spectral reflectance of winter wheat leaves in band i of all samples, yn is the SPAD value of winter wheat leaves in the nth sample, and y ¯ is the average SPAD value of winter wheat leaves in all samples. The value of determination coefficient is equal to R2.
The calculation formula of root mean square error (RMSE) is shown in Formula (2).
R M S E = ∑ n = 1 N ( y n t − y n p ) 2 N
In Formula (2), RMSE is root mean square error, N is the total number of samples, ynt is the measured value of winter wheat leaf SPAD value, ynp is the estimated value of winter wheat leaf SPAD value.
The calculation formula of RPD is shown in Formula (3).
RPD = SD/RMSEV
In Formula (3), RPD is the ratio of performance to deviation, SD is the standard deviation of SPAD value of winter wheat leaves in the validation set.
Generally speaking, the higher the value of the R2, the higher the correlation between the predicted value and the real value. The RMSE value should be as low as possible. The stronger the model’s stability, the closer the RMSET and RMSEV are. RPD can be categorized into three categories: The model is useful for accurately predicting the SPAD value of winter wheat when RPD is equal to or greater than 2.0. The model’s reliability can be enhanced by fine-tuning when the RPD is below 2.0. The model is unstable when RPD is equal or falls below 1.4 [72].

3. Result

3.1. Sensitive Spectral Indices Screening

By analyzing the correlation between these spectral indices and SPAD value, the sensitive spectral indices set was selected to predict the SPAD value of winter wheat.
Figure 2 shows the changes of 62 spectral indices and correlation coefficient of leaf SPAD value in jointing, heading, flowering, filling, milk maturity, maturity and whole growth stage of winter wheat. Each vertical line in the figure corresponds to a spectral index. The seven points on the vertical line represent six single growth periods and whole growth periods. The length of the vertical line is the maximum and minimum difference between spectral index and correlation coefficient of chlorophyll concentration in six single growth periods and whole growth periods. The association between each spectral index and SPAD value is more or less affected by different growth periods, as seen in the figure. The association between spectral index V11 and chlorophyll concentration at heading stage is −0.11594, but the association between spectral index V11 and leaf SPAD value at maturity stage is as high as −0.919, which fluctuates greatly. We found that the vertical line length corresponding to each spectral index ranges from 0.1 to 1.1. The sensitive spectral indices selected by comprehensive consideration can reflect the differences of different growth periods, and it is applicable to all growth periods. It is based on the spectral indices corresponding to the maximum value of positive correlation coefficient and the minimum value of negative correlation coefficient between each growth period and SPAD value, to screen the sensitive spectral indices.
Table 4 lists the spectral indices and correlation coefficients corresponding to the maximum positive correlation coefficient and the minimum negative correlation coefficient of the SPAD value in different growth periods. The maximum positive and negative correlation coefficients at jointing, milk maturity, maturity and whole growth stages are higher, above ±0.7, but the maximum positive and negative correlation coefficients at heading and filling are slightly smaller, basically below ±0.5. According to Table 3, 11 spectral indices V7, V8, V13, V24, V25, V26, V27, V35, V42, V49 and V50 are finally selected as sensitive spectral indices. The bands involved in these sensitive spectral indices are basically in the range of 670–755 nm, which are concentrated near the red edge, and some spectral indices also involve the 550 nm green peak band. This shows that the SPAD value has a great correlation with the red edge information and green peak information. The hyperspectral differences of winter wheat leaves at various growth stages are also mainly reflected in the vicinity around the green peak and red edge.

3.2. Estimation of SPAD Value Based on Sensitive Spectral Indices Set

The value range of each sensitive spectral index is quite different. In the first place, the normalized spectral indices are employed as the input value of the PLSR model to complete the subsequent SPAD value estimation. Table 5 shows the estimated findings of PLSR chlorophyll concentration for each growth period. Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9 demonstrate a scatter diagram of the estimation results of the SPAD value in various growth phases, in which figure (a) is the estimation results of training set and figure (b) is the estimation results of the validation set. The results showed that the estimation results of the heading stage, flowering stage and filling stage were poor, with RPD value less than 1.4 and R2V value less than 0.1. The estimation results of the jointing stage and milk ripening stage were better, RPD value was greater than 1.5 and R2V value was greater than 0.6. The estimation results of the mature stage and whole growth stage are the best, RPD value is basically above 2, and especially in the mature stage, the R2V value was up to 0.974.

3.3. Estimation of SPAD Value Based on Different Pretreatment Approaches

There are great differences between the estimated results of the preprocessed spectral data model in different growth stages. Table 6 lists the best model’s estimation findings under 20 kinds of combined pretreatments in different growth periods. Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15 and Figure 16 are scatter plots of the estimation results of the preprocessing optimal model in six single growth stages, as well as the entire growth stage. R2V growth value and RDP growth value in Table 6 refer to the growth values of R2V and RDP, respectively, compared with the model estimation results of the original spectral data.
Table 6 and Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15 and Figure 16 demonstrate that the chlorophyll estimation results of the optimal model in different growth periods are basically obtained by preprocessing WPD-R’-PCADR or WPD-(1/R)-PCADR, and the R2V value can be up to 0.970, indicating that the combination of WPD denoising, PCA dimensionality reduction and first derivative mathematical form transformation can well increase the PLSR SPAD value estimate model’s accuracy. This is because the first derivative transformation can highlight useful information, denoising can decrease the noise in the original spectral information, and dimensionality reduction can better select useful information and avoid overfitting in model training due to data redundancy. In addition, the chart also shows that the RPD value in the maturity stage and the whole growth stage is greater than 2, and the model R2V value is above 0.75; especially in the mature stage, the R2V value was up to 0.97, which can be a good estimate of the SPAD value of winter wheat. The RPD value of jointing stage is greater than 1.4, which can also better estimate the SPAD value of winter wheat. However, the estimation results of the SPAD value in the heading stage, flowering stage, filling stage and milk maturity stage were extremely poor, R2V was less than 0.4, and the corresponding RPD value was less than 1.4.

3.4. Comparison of SPAD Value Estimation Results of Different Model Inputs

Taking the sensitive spectral indices or the preprocessed spectral data as the model input value, the estimation results have some similarities, which are mainly reflected in different growth periods.
Firstly, the estimation accuracy of the mature stage and whole growth stage is very high, RPD value is greater than 2, and R2V can be up to 0.97 in the mature stage. The estimation accuracy of the jointing stage and milk maturity stage is relatively high, RPD value is basically greater than 1.4, and the accuracy R2V value is also greater than 0.5. The model’s accuracy in predicting the stages of heading, flowering, and filling are relatively poor, R2V value is less than 0.4, and the estimation results are obviously underfitted (Table 5 and Table 6). The estimation accuracy of different growth periods varies widely, and one of the influencing factors is the original sample set. The number of samples collected in every growth period is similar; however, the STDchl in the sample set is substantially different (Table 1). From the turning green to jointing and then to the heading stage, the temperature increases continuously, the winter wheat leaves grow continuously, and the top leaf size also increases continuously. During this period, the chlorophyll content increases sharply. From filling to the mature stage, the grain in the wheat ear takes shape, its reproductive growth is dominant, and the vegetative organs such as leaves gradually stop growing and begin to wither and fall off. Coupled with the damage of high temperature to leaf cells, the chlorophyll content in the leaves decreases significantly. Therefore, the sample’s SPAD value measured in these two periods has a large range, and the STDchl in the sample set is also large. The standard deviation of the SPAD value in the sample set at the jointing stage, milk maturity stage, maturity stage and whole growth stage is greater than 5.5. From heading stage to filling stage, the number and size of leaves basically did not change, and the chlorophyll content increased slowly. Therefore, the range of SPAD value in the sample set measured during the growth period is small, and the STDchl in the sample set is also small. The standard deviation of the SPAD value in the heading, flowering and filling stage is less than 4.
Secondly, using preprocessed spectral data or sensitive spectral indices set as the model input value can significantly boost the precision and stability of the PLSR-based SPAD value estimating model, and when the STDchl of all samples in the model sample set is greater than 5.5 (Table 1, Table 5 and Table 6), the estimation accuracy is higher (the estimation accuracy is greater than 0.5), and the model is stable (RPD is greater than 1.4). Moreover, when STDchl in the sample set is greater than 5.5, the larger the SPAD value in the sample set, the more accurate the model’s estimation. Most important of all, the STDchl in the sample set and the estimation accuracy of the model from high to low are the mature stage, whole growth stage, jointing stage and milk mature stage. It can be seen that the SPAD value estimation model based on PLSR is particularly responsive to the input sample set’s STDchl in addition to the high requirements for the spectral input value of the model (appropriate spectral data preprocessing or appropriate spectral indices set). In addition, there are great differences in the estimation results of single growth period data. If only the data of a single growth period are used to build a model to estimate the SPAD value, the universality of the estimation results is too weak. It is suggested to collect the data of mixed growth period (such as the whole growth period in this paper) for relevant estimation research, and the stability of the model is stronger.
The estimation results of input values of different models also have some differences. When the STDchl in the input value sample set is greater than 6, taking the sensitive spectral indices as the input value is more accurate than taking the preprocessed spectral data. The calculation formula of spectral indices involves the mathematical form transformation of spectral data. Taking spectral indices as the model’s input value also has a certain function of data compression and dimensionality reduction. Therefore, the spectral indices set selected in this study can significantly aid in improving the model’s accuracy and reliability.

3.5. Estimation of SPAD Value with Different Chlorophyll Standard Deviations

For the sake of further confirming the estimation model’s sensitivity to the STDchl in the input sample set, we estimated the SPAD value in different standard deviation sample sets. Within this investigation, the SPAD value in all samples ranges from 2.5 to 76.8. Taking 5 as the gradient and making statistics in the range of 0 to 80 (Figure 17), it is found that the SPAD value of most samples is concentrated in the range of 40 to 60. Based on the principle that the sample set has as many standard deviations (STD) of different SPAD values as possible, the sample set with different SPAD value ranges is set. The range of SPAD value in the sample set ranges from 5 to 80. The spectral data of each sample set pretreated by WPD-(1/R)-PCADR are used as the input value of the model, and the model estimation accuracy of different SPAD value ranges was analyzed using the PLSR SPAD value estimation model (Table 7).
The research shows that when the SPAD value range in the sample set is small (50~55, 50~60), the standard deviation of the sample set is also small, and the corresponding model R2V value is also poor (R2V is less than 0.2). As the range of SPAD values in the sample set gradually increases, the standard deviation also increases, and the corresponding model estimation accuracy (R2V) basically shows a gradually increasing trend, and the more stable the model estimation result is (RPD). Figure 18 shows the variation law of RPD with the STDchl in the sample set, and Figure 19 shows the variation law of R2V with the STDchl in the sample set. It can be seen from the two figures that the STDchl in the sample set is highly correlated with the PLSR SPAD value estimation model’s accuracy and reliability. When the STDchl in the sample set is larger than 5.5, RPD is greater than 1.4 and R2V is higher than 0.5, the model is stable, and its accuracy is high. At this time, the PLSR model could be used to predict the SPAD value with high accuracy.

4. Discussion

4.1. Pretreatment Methods

In this study, common spectral preprocessing methods were compared, and it was found that using both sensitive spectral indices and preprocessed spectral data as model input values could increase the model’s estimation accuracy and reliability to some extent. This also illustrates the importance of spectral data preprocessing. In particular, it can be seen from Table 6 that, no matter in which growth stages, the model inversion accuracy was highest when the pretreatment method was WPD-R’-PCADR or WPD-(1/R)’-PCADR. When WPD and PCADR are used combined, the benefits of the derivative technique become clear. The derivative approach, according to Oldham et al. and Li et al., is not only a useful tool for studying spectral data, but it also aids in the resolution of multicollinearity issues [73,74]. As a result, it can be utilized to improve the analysis’s sensitivity. It also aids in the reduction of noise to some extent. The first derivative (FD) provides the reflectance spectrum’s slope, whereas the second derivative (SD) provides the reflectance spectrum’s slope change [74]. Although SD allows more absorption peaks to be separated, it also introduces noise and may lead to errors. Through SD and FD changes, the spectrum is severely deformed and leads to sharp peaks. For the pretreatment of spectral data, the fractional derivative can reduce the range of spectral variation to a certain extent and retain the form properties of the original spectrum as far as possible, which is more useful than the integer derivative (FD, SD). Fractional derivatives can really be utilized to include further valuable data coming from remote sensing by extending the order to non-integers, which could also add more details to the spectrum than integer derivatives. Only FD was employed to measure SPAD values in this research. Other higher-order integer and fractional derivatives will be used to preprocess spectral data in the future.

4.2. Influence of Sample Set on Model Accuracy

The results of this paper show that (Figure 18 and Figure 19) the STDchl in the model input sample set is highly sensitive to the accuracy of the PLSR model in evaluating the SPAD value. This has to do with the characteristics of the PLSR method itself. PLSR is a machine learning algorithm, which has high requirements for the model input sample set. It needs the number of samples to be as large and diverse as possible. When the sample set is too small or the data range is too monotonous, it is easy to lead to the phenomenon of underfitting and overfitting, and the generalization ability of the model is weak. In this study, the STDchl in the sample set is too small, which is the performance of poor sample diversity. The SPAD values in the heading stage, flowering stage and milk ripening stage are relatively concentrated due to their own physiological characteristics. The standard deviation of the SPAD value in the corresponding sample set is less than 4, and the accuracy of the corresponding model is also low. However, there are not enough data to prove whether the accuracy tends to be stable when the standard deviation is large to a certain value. More growth period data need to be collected for further verification.
Table 7 reveals that when the STDchl in the sample set is greater than 5.5, the larger the STDchl in the sample set, the higher the estimation accuracy of the model. Due to the small number of samples collected from the research, whether the threshold of 5.5 can be applied to other scenes needs to be further verified. In addition, based on the biological characteristics of winter wheat itself, the content of chlorophyll in leaves has a certain range. The STDchl in the sample set we collected will always reach an upper limit. At this time, the relationship between the STDchl and the accuracy of the model also needs to be further explored and verified. In following research, we will conduct deeper exploration, research and verification on the relationship between the sample set and model accuracy.

4.3. Universality of Data

Taking winter wheat as the research object and taking the hyperspectral data and corresponding SPAD value in different growth periods as the data source, this paper establishes a model by using the PLSR method to estimate the SPAD value and analyzes the estimation results in detail. The research has achieved some results, but the scope of the study area is small. When the variety of winter wheat changes or the growing environment changes, the research results need to be further analyzed and verified. In addition, it is relatively difficult to obtain data at the leaf scale of winter wheat, because it is susceptible to weather (when there are clouds and winds, etc.), and data collectors need certain professional experience. At present, in the quantitative remote sensing’s field, combining ground non-imaging spectral data with UAV imaging spectral data and satellite image multi-spectral data to realize multi-level and multi-scale crop growth monitoring is an unavoidable tendency. Wang et al. combined ground hyperspectral data with GF1 satellite multispectral images to achieve an inversion of chlorophyll concentration geographic variation in the winter wheat canopy [75]. In the subsequent study, we attempted to expand the research scope and conduct satellite-scale estimation of the winter wheat SPAD value, so as to analyze and find efficient and high-precision estimation methods for winter wheat chlorophyll that are suitable for more conditions.

5. Conclusions

Depending on the spectral data and SPAD value in various growth stages, taking the sensitive spectral indices and preprocessed spectral data as the model input values, the model was constructed by the PLSR method to estimate the SPAD value, and the estimated results were compared and analyzed in detail. The main conclusions are as follows:
(1)
The 11 spectral indices (V13, V15, V16, V17, V25, V39, V41, V42, V43, V47 and V57) selected in this paper can be used as the input values of the model, which could boost the model’s precision and reliability while calculating the SPAD value.
(2)
Compared with the original spectral data and preprocessed spectral data as the model’s input value, especially spectral data preprocessed by WPD-(1/R)-PCADR or WPD-R’-PCADR, it can increase the accuracy of the SPAD value estimation model and enhance the model’s stability.
(3)
When the STDchl in the sample set is less than 4, the estimation results of the SPAD value are prone to underfitting; when the STDchl in the sample set is greater than 5.5, the greater the STDchl in the sample set, the higher the model’s estimating accuracy. In this case, the advantage of using sensitive spectral indices and preprocessing spectral dataset as model input values to increase the estimation model’s precision and stability is obvious. In addition, when the STDchl in the sample set was greater than 6, the model with the sensitive spectral indices as the model input had a greater estimation accuracy than the model with the pretreatment spectral dataset.
(4)
When using the spectral data and corresponding SPAD value data of a single growth period as the data source, the estimation results were not representative when using the PLSR method to predict the SPAD value, and the universality of data related to other growth periods was poor. Modeling using data from the whole growth period can improve the universality ability and stability of the model. It is recommended that spectral data for the whole fertility period be used as the data source and that the STDchl corresponding to the data source should be as large as possible (standard deviation of SPAD values > 5.5).

Author Contributions

Conceptualization, L.S., M.G. and J.Y.; methodology, L.S.; software, L.S.; validation, L.S.; formal analysis, L.S.; investigation, L.S.; resources, L.S., H.S., M.G. and Q.W.; data curation, L.S. and Q.W.; writing—original draft preparation, L.S.; writing—review and editing, M.G. and J.Y.; supervision, H.S., M.G. and J.Y.; project administration, M.G.; funding acquisition, M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by Central Public-interest Scientific Institution Basal Research Fund (grant number Y2022GH05), Fundamental Research Funds for Central Non-profit Scientific Institution (grant number 1610132021024) and Shantou University Team Building Project for Innovative and Strong University (grant number 2017KCXTD015).

Acknowledgments

The data of this field experiment were obtained with the assistance of students Ya Gao, Qiang Li, Shilei Li, Renjun Li, Xiaolin Zhu and Jingjing Yang from the Institute of agricultural resources and regional planning, Chinese Academy of Agricultural Sciences. We would like to express our great thanks to them.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Guo, Y.K.; Xu, M.; Zhang, X.J.; Liu, Y.L. Hyperspectral inversion of chlorophyll content combined with PRO-4SAIL and BP neural network. Bull. Surv. Mapp. 2020, 516, 24–27. [Google Scholar]
  2. Luo, Q.F. Research progress of chlorophyll and development and utilization of chlorophyll copper Na. For. Prod. Chem. Commun. 1995, 3, 32–33. [Google Scholar]
  3. Xiao, J. Research on Plant Chlorophyll Inversion Based on Geometric Spectral Integration; Wuhan University: Wuhan, China, 2018. [Google Scholar]
  4. Qian, B.; Huang, W.; Huichun, Y. Inversion of winter wheat chlorophyll contents based on improved algorithms for red edge position. Nongye Gongcheng Xuebao Trans. Chin. Soc. Agric. Eng. 2021, 36, 162–170. [Google Scholar]
  5. Jiang, L.F.; Shi, F.C.; Wang, H.T.; Zu, Y.G. Application tryout of chlorophyll meter SPAD-502. J. Ecol. 2005, 12, 1543–1548. [Google Scholar]
  6. Tong, Q.X.; Zhang, B.; Zhang, L.F. Current progress of hyperspectral remote sensing in China. J. Remote Sens. 2016, 20, 689–707. [Google Scholar]
  7. Wang, Y.Y.; Chen, Y.H.; Li, J.; Huang, W.J. Two new red border parameters indicating the severity of winter wheat stripe rust. J. Remote Sens. 2007, 6, 875–881. [Google Scholar]
  8. Li, D.; Chen, J.M.; Zhang, X.; Yan, Y.; Zhu, J.; Zheng, H.B.; Zhou, K.; Yao, X.; Tian, Y.C.; Zhu, Y.; et al. Improved estimation of leaf chlorophyll content of row crops from canopy reflectance spectra through minimizing canopy structural effects and optimizing off-noon observation time. Remote Sens. Environ. 2020, 248, 111985. [Google Scholar] [CrossRef]
  9. Wang, B.Z.; You, Y.M.; Lu, M.Z.; Zhao, Z.Y. Physical World Tour; Beijing Institute of Technology Press: Beijing, China, 2015; Volume 12, p. 101. [Google Scholar]
  10. Jin, Z.Y.; Tian, Q.J.; Hui, F.M. Study on the Relationship between chlorophyll concentration and spectral reflectance in rice. Remote Sens. Technol. Appl. 2003, 18, 134–137. [Google Scholar]
  11. Lu, X.; Peng, H. Predicting Cherry Leaf Chlorophyll Concentrations Based on Foliar Reflectance Spectra Variables. J. Indian Soc. Remote Sens. 2015, 43, 109–120. [Google Scholar] [CrossRef]
  12. Zhang, X.R.; Feng, M.C.; Yang, W.D.; Wang, C.; Guo, X.L.; Shi, C.C. Using spectral transformation processes to estimate chlorophyll content of winter wheat under low temperature stress. Chin. J. Eco-Agric. 2017, 25, 1351–1359. [Google Scholar]
  13. Bonham-Carter, G.F. Numerical procedures and computer program for fitting an inverted gaussian model to vegetation reflectance data. Comput. Geosci. 1988, 14, 339–356. [Google Scholar] [CrossRef]
  14. Yang, J.; Tian, Y.C.; Yao, X. Hyperspectral estimation model of chlorophyll content in upper leaves of rice. Acta Ecol. Sin. 2009, 29, 6561–6571. [Google Scholar]
  15. Wang, T.L.; Gao, M.F.; Cao, C.L.; You, J.; Zhang, X.W.; Shen, L.Z. Winter wheat chlorophyll content retrieval based on machine learning using in situ hyperspectral data. Comput. Electron. Agric. 2022, 193, 106728. [Google Scholar] [CrossRef]
  16. Qiu, Z.J.; Song, H.Y.; He, Y.; Fang, H. Variation rules of the nitrogen content of the rape growth stage using SPAD and spectial. Trans. CSAE 2007, 23, 150–154. [Google Scholar]
  17. Jin, Q. Extraction of Chlorophyll Hyperspectral Characteristics and Comprehensive Evaluation of Estimation Model in Winter Wheat Under Low Temperature Stress; Shanxi Agricultural University: Jinzhong, China, 2017. [Google Scholar]
  18. Srivastava, P.K.; Gupta, M.; Singh, U.; Prasad, R.; Pandey, P.C.; Raghubanshi, A.S.; Petropoulos, G.P. Sensitivity analysis of artificial neural network for chlorophyll prediction using hyperspectral data. Environ. Dev. Sustain. 2021, 23, 5504–5519. [Google Scholar] [CrossRef]
  19. Shen, L.Z.; Gao, M.F.; Yan, J.W.; Li, Z.L.; Duan, S.B. Hyperspectral estimation of soil organic matter content using different spectral preprocessing techniques and PLSR method. Remote Sens. 2020, 12, 1206. [Google Scholar] [CrossRef]
  20. Atzberger, C.; Gu’erif, 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]
  21. Zha, Y.S.; Lu, Z.L.; Zhang, X. Design of Micro Solar Sensor for Cubic Star. Chin. J. Sens. Actuators 2018, 31, 31–35. [Google Scholar]
  22. Wang, H. Study on Growth Metabolism and Canopy Spectral Characteristics of MAIZE Under Different Fertilization Systems; Shandong Agricultural University: Taian, China, 2015. [Google Scholar]
  23. He, Y.; Feng, L. Fundamentals of Geospatial Informatics; Zhejiang University Press: Hangzhou, China, 2010. [Google Scholar]
  24. Tucker, C.J. Red and photographic infrared linear combinations for monitoring vegetation. Remote Sens. Environ. 1979, 8, 127–150. [Google Scholar] [CrossRef]
  25. Gitelson, A.A.; Vina, A.; Ciganda, V. Remote estimation of canopy chlorophyll content in crops. Geophys. Res. Lett. 2005, 32, L08403. [Google Scholar] [CrossRef]
  26. Zarco-Tejada, P.J.; Berjon, A.; Lopezlozano, R.; Miller, J.; Martin, P.; Cachorro, V.; Gonzalez, M.; Defrutos, 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]
  27. Peñuelas, 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]
  28. 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]
  29. 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]
  30. Jordan, C.F. Derivation of Leaf-Area Index from Quality of Light on the Forest Floor. Ecology 1969, 50, 663–666. [Google Scholar] [CrossRef]
  31. McMurtrey, J.E.; Chappelle, E.W.; Kim, M.S.; Meisinger, J.J.; Corp, L.A. Distinguishing nitrogen fertilization levels in field corn (Zea mays L.) with actively induced fluorescence and passive reflectance measurements. Remote Sens. Environ. 1994, 47, 36–44. [Google Scholar] [CrossRef]
  32. Mutanga, O.; Skidmore, A.K. Narrow band vegetation indices overcome the saturation problem in biomass estimation. Int. J. Remote Sens. 2004, 25, 3999–4014. [Google Scholar] [CrossRef]
  33. Zarco-Tejada, P.J.; Pushnik, J.C.; Dobrowski, S.; Ustin, S.L. Steady-state chlorophyll a fluorescence detection from canopy derivative reflectance and double-peak red-edge effects. Remote Sens. Environ. 2003, 84, 283–294. [Google Scholar] [CrossRef]
  34. Gitelson, A.; Merzlyak, M.N. Quantitative estimation of chlorophyll-a using reflectance spectra: Experiments with autumn chestnut and maple leaves. J. Photochem. Photobiol. B Biol. 1994, 22, 247–252. [Google Scholar] [CrossRef]
  35. Gitelson, A.A.; Buschmann, C.; Lichtenthaler, H.K. The chlorophyll fluorescence ratio F735_F700 as an accurate measure of the chlorophyll content in plants. Remote Sens. Environ. 1999, 69, 296–303. [Google Scholar] [CrossRef]
  36. Datt, B. Visible/near infrared reflectance and chlorophyll content in Eucalyptus leaves. Int. J. Remote Sens. 2010, 20, 2741–2759. [Google Scholar] [CrossRef]
  37. Datt, B. Remote sensing of chlorophyll a, chlorophyll b, chlorophyll a + b, and total carotenoid content in eucalyptus leaves. Remote Sens. Environ. 1998, 66, 111–121. [Google Scholar] [CrossRef]
  38. Gitelson, A.A.; Merzlyak, M.N. Remote estimation of chlorophyll content in higher plant leaves. Int. J. Remote Sens. 1997, 18, 2691–2697. [Google Scholar] [CrossRef]
  39. Lichtenthaler, H.K.; Lang, M.; Sowinska, M.; Heisel, F.; Miehé, J.A. Detection of Vegetation Stress Via a New High Resolution Fluorescence Imaging System. J. Plant Physiol. 1996, 148, 599–612. [Google Scholar] [CrossRef]
  40. Vogelmann, J.E.; Rock, B.N. Moss, D.M. Red edge spectral measurements from sugar maple leaves. Int. J. Remote Sens. 1993, 14, 1563–1575. [Google Scholar] [CrossRef]
  41. Boochs, F.; Kupfer, G.; Dockter, K.; KÜHbauch, W. Shape of the red edge as vitality indicator for plants. Int. J. Remote Sens. 1990, 11, 1741–1753. [Google Scholar] [CrossRef]
  42. Elvidge, C.; Chen, Z. Comparison of broad-band and narrow-band red and near-infrared vegetation indices. Remote Sens. Environ. 1995, 54, 38–48. [Google Scholar] [CrossRef]
  43. Liu, J.G.; Moore, J.M. Hue image RGB colour composition. A simple technique to suppress shadow and enhance spectral signature. Int. J. Remote Sens. 1990, 11, 1521–1530. [Google Scholar] [CrossRef]
  44. 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]
  45. Filella, I.; Penuelas, J. The red edge position and shape as indicators of plant chlorophyll content, biomass and hydric status. Int. J. Remote Sens. 1994, 15, 1459–1470. [Google Scholar] [CrossRef]
  46. Zarco-Tejada, P.J.; Miller, J.R.; Noland, T.L.; Mohammed, G.H.; Sampson, P.H. Scaling-Up and Model Inversion Methods with Narrowband Optical Indices for Chlorophyll Content Estimation in Closed Forest Canopies with Hyperspectral Data. IEEE Trans. Geosci. Remote Sens. 2001, 39, 1491–1507. [Google Scholar] [CrossRef]
  47. Croft, H.; Chen, J.M.; Zhang, Y. The applicability of empirical vegetation indices for determining leaf chlorophyll content over different leaf and canopy structures. Ecol. Complex. 2014, 17, 119–130. [Google Scholar] [CrossRef]
  48. Gandia, S.; Fernández, G.; García, J.C.; Moreno, J. Retrieval of vegetation biophysical variables from CHRISPROBA data in the SPARC campaign. In Proceedings of the Proba Workshop, ESA/ESRIN, Frascati, Italy, 28–30 January 2004. [Google Scholar]
  49. Maccioni, A.; Agati, G.; Mazzinghi, P. New vegetation indices for remote measurement of chlorophylls based on leaf directional reflectance spectra. J. Photochem. Photobiol. B Biol. 2001, 61, 52–61. [Google Scholar] [CrossRef]
  50. Dash, J.; Curran, P.J. The MERIS terrestrial chlorophyll index. Int. J. Remote Sens. 2004, 25, 5403–5413. [Google Scholar] [CrossRef]
  51. Gitelson, A.A.; Kaufman, Y.J.; Merzlyak, M.N. Use of a green channel in remote sensing of global vegetation from EOS–MODIS. Remote Sens. Environ. 1996, 58, 289–298. [Google Scholar] [CrossRef]
  52. Chappelle, E.W.; Kim, M.S.; Iii, M.M. 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]
  53. Carter, G.A. Ratios of leaf reflectances in narrow wavebands as indicators of plant stress. Int. J. Remote Sens. 1994, 15, 697–703. [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. Gitelson, A.A.; Gritz, Y.; Merzlyak, M.N. Relationships between leaf chlorophyll content and spectral reflectance and algorithms for non-destructive chlorophyll assessment in higher plant leaves. J. Plant Physiol. 2003, 160, 271–282. [Google Scholar] [CrossRef]
  56. 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]
  57. Haboudane, D.; Miller, J.R.; Pattey, E.; 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]
  58. Daughtry, C.S.T.; Walthall, C.L.; Kim, M.S.; 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]
  59. 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]
  60. Huete, A.R.; Liu, H.Q.; Batchily, K.; Leeuwen, W. A comparison of vegetation indices global set of TM images for EOS-MODIS. Remote Sens. Environ. 1997, 59, 440–451. [Google Scholar] [CrossRef]
  61. Wu, C.; Niu, Z.; Tang, Q.; Huang, W. Estimating chlorophyll content from hyperspectral vegetation indices: Modeling and validation. Agric. For. Meteorol. 2008, 148, 1230–1241. [Google Scholar] [CrossRef]
  62. Maire, G.L.; François, C.; Dufrêne, E. Towards universal broad leaf chlorophyll indices using PROSPECT simulated database and hyperspectral reflectance measurements. Remote Sens. Environ. 2004, 89, 1–28. [Google Scholar] [CrossRef]
  63. Kong, L.J. Matlab Wavelet Analysis Super Learning Manual; The People’s Posts and Telecommunications Press: Beijing, China, 2014; pp. 300–301. [Google Scholar]
  64. Virmani, J.; Kumar, V.; Kalar, N. SVM-Based Characterization of Liver Ultrasound Images Using Wavele Packet Texture Descriptors. J. Digit. Imaging 2012, 26, 530–543. [Google Scholar] [CrossRef]
  65. Han, Z.Y.; Lun, S.X.; Wang, J. Research on Robust Feature Extraction and Visualization of Speech Signals; Northeastern University Press: Shenyang, China, 2012. [Google Scholar]
  66. Shen, L.Z.; Gao, M.F.; Yan, J.W.; Yao, Y.M. Estimation model of soil organic matter based on SVR and PLSR. Agric. Inf. China 2019, 31, 58–71. [Google Scholar]
  67. Pham, Q.T.; Liou, N.S. The development of on-line surface defect detection system for jujubes based on hyperspectral images. Comput. Electron. Agric. 2022, 194, 106743. [Google Scholar] [CrossRef]
  68. Cao, C.L.; Wang, T.L.; Gao, M.F.; Li, Y.; Li, D.D.; Zhang, H.J. Hyperspectral inversion of nitrogen content in maize leaves based on different dimensionality reduction algorithms. Comput. Electron. Agric. 2021, 190, 106461. [Google Scholar] [CrossRef]
  69. Juan, P.R.; Jochem, V.; Jordi, M. Hyperspectral dimensionality reduction for biophysical variable statistical retrieval. ISPRS J. Photogramm. Remote Sens. 2017, 132, 88–101. [Google Scholar]
  70. Pearson, K. On lines and planes of closest fit to systems of points in space. Phil. Mag. 1901, 2, 559–572. [Google Scholar] [CrossRef]
  71. Shi, L.; Li, Z.L.; Cui, G.M. A Model of Silicon Content in Hot Metal of Blast Furnace Based on Partial Least Square Regression. J. Inn. Mong. Univ. Nat. Sci. Ed. 2010, 41, 427–430. [Google Scholar]
  72. Chang, C.W.; Laird, D.A.; Mausbach, M.J.; Hurburgh, C.R. Near infrared reflectance spectroscopy-Principal components regression analysis of soil properties. Soil Sci. Soc. Am. J. 2001, 65, 480–490. [Google Scholar] [CrossRef]
  73. Oldham, K.B.; Spanier, J. The fractional calculus. Math. Gaz. 1974, 56, 396–400. [Google Scholar]
  74. Li, B.; Xie, W. Adaptive fractional differential approach and its application to medical image enhancement. Comput. Electr. Eng. 2015, 45, 324–335. [Google Scholar] [CrossRef]
  75. Wang, T.T. Remote Sensing Inversion of Winter Wheat Chlorophyll Based on Hyperspectral and GF-1 Satellite Images; Northwest A&F University: Shaanxi, China, 2019. [Google Scholar]
Remotesensing 14 04660 g001 550
Figure 1. Overview of the study area.
Figure 1. Overview of the study area.
Remotesensing 14 04660 g001
Remotesensing 14 04660 g002 550
Figure 2. The association between spectral indices and SPAD value at different growth stages.
Figure 2. The association between spectral indices and SPAD value at different growth stages.
Remotesensing 14 04660 g002
Remotesensing 14 04660 g003 550
Figure 3. The estimated results of SPAD value for (a) training set and (b) validation set at jointing stage.
Figure 3. The estimated results of SPAD value for (a) training set and (b) validation set at jointing stage.
Remotesensing 14 04660 g003
Remotesensing 14 04660 g004 550
Figure 4. The estimated results of SPAD value for (a) training set and (b) validation set at heading stage.
Figure 4. The estimated results of SPAD value for (a) training set and (b) validation set at heading stage.
Remotesensing 14 04660 g004
Remotesensing 14 04660 g005 550
Figure 5. The estimated results of SPAD value for (a) training set and (b) validation set at flowering stage.
Figure 5. The estimated results of SPAD value for (a) training set and (b) validation set at flowering stage.
Remotesensing 14 04660 g005
Remotesensing 14 04660 g006 550
Figure 6. The estimated results of SPAD value for (a) training set and (b) validation set at filling stage.
Figure 6. The estimated results of SPAD value for (a) training set and (b) validation set at filling stage.
Remotesensing 14 04660 g006
Remotesensing 14 04660 g007 550
Figure 7. The estimated results of SPAD value for (a) training set and (b) validation set at milk stage.
Figure 7. The estimated results of SPAD value for (a) training set and (b) validation set at milk stage.
Remotesensing 14 04660 g007
Remotesensing 14 04660 g008 550
Figure 8. The estimated results of SPAD value for (a) training set and (b) validation set at mature stage.
Figure 8. The estimated results of SPAD value for (a) training set and (b) validation set at mature stage.
Remotesensing 14 04660 g008
Remotesensing 14 04660 g009 550
Figure 9. The estimated results of SPAD value for (a) training set and (b) validation set at whole growth stage.
Figure 9. The estimated results of SPAD value for (a) training set and (b) validation set at whole growth stage.
Remotesensing 14 04660 g009
Remotesensing 14 04660 g010 550
Figure 10. The estimated results of SPAD value for (a) training set and (b) validation set based on pretreatment WPD-(1/R)’-PCADR at jointing stage.
Figure 10. The estimated results of SPAD value for (a) training set and (b) validation set based on pretreatment WPD-(1/R)’-PCADR at jointing stage.
Remotesensing 14 04660 g010
Remotesensing 14 04660 g011 550
Figure 11. The estimated results of SPAD value for (a) training set and (b) validation set based on pretreatment WPD-1/R-PCADR at heading stage.
Figure 11. The estimated results of SPAD value for (a) training set and (b) validation set based on pretreatment WPD-1/R-PCADR at heading stage.
Remotesensing 14 04660 g011
Remotesensing 14 04660 g012 550
Figure 12. The estimated results of SPAD value for (a) training set and (b) validation set based on pretreatment WPD-(1/R)’-PCADR at flowering stage.
Figure 12. The estimated results of SPAD value for (a) training set and (b) validation set based on pretreatment WPD-(1/R)’-PCADR at flowering stage.
Remotesensing 14 04660 g012
Remotesensing 14 04660 g013 550
Figure 13. The estimated results of SPAD value for (a) training set and (b) validation set based on pretreatment ND-(1/R)’-PCADR at filling stage.
Figure 13. The estimated results of SPAD value for (a) training set and (b) validation set based on pretreatment ND-(1/R)’-PCADR at filling stage.
Remotesensing 14 04660 g013
Remotesensing 14 04660 g014 550
Figure 14. The estimated results of SPAD value for (a) training set and (b) validation set based on pretreatment WPD-(1/R)’-PCADR at milk stage.
Figure 14. The estimated results of SPAD value for (a) training set and (b) validation set based on pretreatment WPD-(1/R)’-PCADR at milk stage.
Remotesensing 14 04660 g014
Remotesensing 14 04660 g015 550
Figure 15. The estimated results of SPAD value for (a) training set and (b) validation set based on pretreatment WPD-R’-PCADR at mature stage.
Figure 15. The estimated results of SPAD value for (a) training set and (b) validation set based on pretreatment WPD-R’-PCADR at mature stage.
Remotesensing 14 04660 g015
Remotesensing 14 04660 g016 550
Figure 16. The estimated results of SPAD value for (a) training set and (b) validation set based on pretreatment WPD-R’-PCADR at whole growth stage.
Figure 16. The estimated results of SPAD value for (a) training set and (b) validation set based on pretreatment WPD-R’-PCADR at whole growth stage.
Remotesensing 14 04660 g016
Remotesensing 14 04660 g017 550
Figure 17. SPAD value distribution.
Figure 17. SPAD value distribution.
Remotesensing 14 04660 g017
Remotesensing 14 04660 g018 550
Figure 18. The correlation diagram of model RPD and STD of sample SPAD values.
Figure 18. The correlation diagram of model RPD and STD of sample SPAD values.
Remotesensing 14 04660 g018
Remotesensing 14 04660 g019 550
Figure 19. The correlation diagram of model R2V and standard deviation of sample SPAD values.
Figure 19. The correlation diagram of model R2V and standard deviation of sample SPAD values.
Remotesensing 14 04660 g019
Table
Table 1. The introduction of winter wheat sample sets on each growth stages.
Table 1. The introduction of winter wheat sample sets on each growth stages.
Growth StagesData Collection DataSample SetNumber of SamplesSPAD Value
MaximumMinimumAverageStandard Deviation
Jointing stage26 April 2019Training set13255.525.243.6455.632
Validation set3256.835.546.3505.481
Total sample set16456.825.244.1735.722
Heading stage9 May 2019Training set15667.245.954.3413.511
Validation set4468.949.755.8233.652
Total sample set20068.945.954.6673.604
Flowering stage13 May 2019Training set7860.746.453.3042.652
Validation set2065.151.756.5853.196
Total sample set9865.146.453.9733.087
Filling stage22 May 2019Training set15670.146.655.1123.187
Validation set3961.443.854.7543.133
Total sample set19570.143.855.0413.180
Milk ripening stage28 May 2019Training set15665.327.050.7966.558
Validation set446217.349.5848.125
Total sample set20065.317.350.5296.969
Maturity stage31 May 2019Training set15663.63.243.01312.580
Validation set39582.542.23613.099
Total sample set19563.62.542.85712.722
Whole growth stage Training set84276.817.351.7026.432
Validation set21063.62.543.57812.592
Total sample set105276.82.550.0798.694
Table
Table 2. Spectral indices used to estimate the vegetation physiological and ecological parameters.
Table 2. Spectral indices used to estimate the vegetation physiological and ecological parameters.
NumberSpectral IndicesFormulaReferenceNumberSpectral IndicesFormulaReference
V1NDVI R 800 − R 670 R 800 + R 670 [24]V2CI R 750 ( R 700 − R 710 ) − 1 [25]
V3GI R 554 R 677 [26]V4SIPI (680) R 800 − R 455 R 800 − R 680 [27]
V5VI(700) R 700 − R 670 R 700 + R 670 [28]V6RDVI ( R 800 − R 670 ) ( R 800 + R 670 ) [29]
V7SR R 800 R 680 [30]V8SR2 R 700 R 670 [31]
V9MNDVI8 ( R 755 − R 730 ) ( R 755 + R 730 ) [32]V10DPI D 688 ∗ D 710 D 697 2 [33]
V11D2-C D 705 D 722 [33]V12D1-CF D 730 D 706 [33]
V13NDVI2-LS,Ca R 750 − R 705 R 750 + R 705 [34]V14Gitelson-LS,CTotal 1 R 700 [35]
V15Datt4-LS,C R 754 R 704 [36]V16Datt1-LS,Cb R 672 R 550 [37]
V17SIPI R 800 − R 445 R 800 − R 680 [27]V18SR5 R 752 R 690 [38]
V19SR3-LS,S R 440 R 690 [39]V20Vogelmann-LS,C R 740 R 720 [40]
V21Boochs1 D 703 [41]V22Boochs2 D 720 [41]
V23SOFDR (625–795) ∑ 625 795 D [42]V24CI R 800 - R 550 R 800 [43]
V25MSR R 800 - R 445 R 680 - R 445 [44]V26REP M a x ( D 680 , D 760 ) [45]
V27BGI R 450 R 550 [46]V28SIPI (705) R 800 - R 455 R 800 + R 705 [47]
V29ZM R 750 R 710 [46]V30Datt-LS,C R 850 - R 710 R 850 - R 680 [36]
V31SR6 R 750 R 550 [38]V32MNDVI1 R 755 - R 745 R 755 + R 745 [32]
V33D690_red D 690 [47]V34NDVI3 R 682 - R 553 R 682 + R 553 [48]
V35Maccioni-LS,C R 780 - R 710 R 780 - R 680 [49]V36MTCI R 754 - R 709 R 709 - R 681 [50]
V37NPCI R 680 - R 430 R 680 + R 430 [27]V38GNDVI-Ca R 800 - R 550 R 800 + R 550 [51]
V39Datt3-LS,C R 850 R 710 [36]V40Datt5-LS,Ca,total R 672 R 550 ∗ R 708 [36]
V41Datt2-LS R 860 R 550 ∗ R 708 [37]V42SR4 R 750 R 700 [38]
V43SR1-LS,Ca R 675 R 700 [52] V44SRPI R 430 R 680 [27]
V45Vogelmann2-LS,C R 715 R 705 [40]V46Vogelmann3-LS,C R 734 - R 747 R 715 + R 726 [40]
V47SOFDR (680–780) ∑ 680 780 D [45]V48Carter-I R 695 R 760 [53]
V49DVI R 800 − R 670 [30]V50PSSRC R 800 R 680 + R 800 R 635 [54]
V51Gitelson-RG R 800 R 550 - 1 [55]V52BI R 800 + R 670 + R 550 3 [43]
V53TCARI 3 ∗ ( ( R 700 − R 670 ) − 0.2 ∗ ( R 700 − R 550 ) ∗ ( R 700 R 670 ) ) [56]
V54MCARI1 1.2 ∗ [ 2.5 ∗ ( R 800 − R 670 ) − 1.3 ∗ ( R 800 − R 550 ) ] [57]
V55MCARI [ ( R 700 − R 670 ) − 0.2 ∗ ( R 700 − R 550 ) ] ( R 700 R 670 ) [58]
V56TVI 0.5 ∗ [ 120 ∗ ( R 750 − R 550 ) − 200 ∗ ( R 670 − R 550 ) ] [59]
V57MCARI2 1.5 ∗ [ 2.5 ∗ ( R 800 − R 670 ) ] − 1.3 ∗ ( R 800 − R 550 ) ( 2 ∗ R 800 + 1 ) − ( 6 ∗ R 800 − 5 ∗ R 670 ) − 0.5 [57]
V58EVI 2.5 ∗ ( R 800 - R 670 R 800 − ( 6 ∗ R 670 ) − ( 7.5 ∗ R 475 ) + 1 ) [60]
V59TCARI2-Wu 3 ∗ ( ( R 750 − R 705 ) − 0.2 ∗ ( R 750 − R 550 ) ∗ ( R 750 R 705 ) ) [61]
V60DD ( R 749 − R 720 ) − ( R 701 − R 672 ) [62]
V61MND (705) ( R 750 − R 705 ) ( R 750 + R 705 − 2 R 445 ) [44]
V62MNDVI ( R 800 − R 680 ) ( R 800 + R 680 − 2 R 445 ) [36]
Table
Table 3. Combination pretreatment.
Table 3. Combination pretreatment.
PreprocessingName
DenoisingND, WPD
Data form transformationR, R’, (1/R)’, 1/R, log(R)
dimension reductionNDR, PCADR
Table
Table 4. Sensitive spectral indices.
Table 4. Sensitive spectral indices.
Growth StagesHighest Positive Correlation Spectral IndicesHighest Positive Correlation CoefficientLowest Negative Correlation Spectral IndicesLowest Negative Correlation Coefficient
Jointing stageV70.795V26−0.742
Heading stageV240.415V42−0.400
Flowering stageV250.553V35−0.521
Filling stageV270.413V35−0.380
Milk ripening stageV490.720V26−0.717
Maturity stageV130.919V8−0.919
Whole growth stageV500.830V8−0.799
Table
Table 5. Estimation results of chlorophyll concentration in various growth stages with spectral indices set as input value.
Table 5. Estimation results of chlorophyll concentration in various growth stages with spectral indices set as input value.
Growth StagesRMSET (SPAD)RMSEV (SPAD)RPDR2TR2V
Jointing stage3.2043.4731.5780.6870.706
Heading stage2.9063.8940.9380.3190.010
Flowering stage2.2153.6090.8850.3100.002
Filling stage2.6463.2140.9750.3150.063
Milk ripening stage4.4225.1361.5820.5520.634
Maturity stage4.5173.6243.6150.8820.974
Whole growth stage3.8946.3991.9560.6340.792
Table
Table 6. The optimal estimation results under the data preprocessing of each growth stages.
Table 6. The optimal estimation results under the data preprocessing of each growth stages.
Growth StagesPretreatment MethodRMSET (SPAD)RMSEV (SPAD)RPDR2TR2VR2V Growth ValueRPD Growth Value
JointingWPD-(1/R)’-PCADR3.3773.4531.5870.6500.7340.3180.578
HeadingWPD-1/R-PCADR2.6163.5321.0340.4500.1800.1780.405
FloweringWPD-(1/R)’-PCADR0.5453.9480.8090.9830.0920.0710.220
FillingND-(1/R)’-PCADR2.7772.6161.1980.2440.3400.3300.225
Milk ripeningWPD-(1/R)’-PCADR4.8816.2661.2970.4520.5120.0050.101
MaturityWPD-R’-PCADR4.1433.7343.5080.9030.9700.1091.170
Whole growthWPD-R’-PCADR2.8446.1072.0620.8060.7850.0440.180
Table
Table 7. The estimation results of different SPAD value ranges.
Table 7. The estimation results of different SPAD value ranges.
SPAD Value RangeSTD (SPAD)RMSET (SPAD)RMSEV (SPAD)RPDR2TR2V
50–551.3931.2991.3051.0650.1330.121
50–602.4242.1592.2651.0830.2030.151
45–603.4622.6412.6531.3010.4200.422
45–653.7483.0183.0701.2230.3520.338
40–654.9813.2713.6571.3620.570.472
40–805.1733.6443.6991.3780.5090.485
35–805.8783.6643.6941.5780.6140.607
30–806.2424.0964.0651.5230.5730.576
25–806.8034.4694.3651.5440.5720.592
15–557.1974.5004.5651.5340.6160.589
15–807.6643.6894.7481.6030.7710.626
00–808.6804.1614.8411.7740.7730.697
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shen, L.; Gao, M.; Yan, J.; Wang, Q.; Shen, H. Winter Wheat SPAD Value Inversion Based on Multiple Pretreatment Methods. Remote Sens. 2022, 14, 4660. https://doi.org/10.3390/rs14184660

AMA Style

Shen L, Gao M, Yan J, Wang Q, Shen H. Winter Wheat SPAD Value Inversion Based on Multiple Pretreatment Methods. Remote Sensing. 2022; 14(18):4660. https://doi.org/10.3390/rs14184660

Chicago/Turabian Style

Shen, Lanzhi, Maofang Gao, Jingwen Yan, Qizhi Wang, and Hua Shen. 2022. "Winter Wheat SPAD Value Inversion Based on Multiple Pretreatment Methods" Remote Sensing 14, no. 18: 4660. https://doi.org/10.3390/rs14184660

APA Style

Shen, L., Gao, M., Yan, J., Wang, Q., & Shen, H. (2022). Winter Wheat SPAD Value Inversion Based on Multiple Pretreatment Methods. Remote Sensing, 14(18), 4660. https://doi.org/10.3390/rs14184660

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