Hyperspectral Estimation of Winter Wheat Leaf Water Content Based on Fractional Order Differentiation and Continuous Wavelet Transform

: Leaf water content (LWC) is one of the important indicators of crop health. It plays an important role in the physiological process of leaves, participates in almost all physiological processes of crops, and is of great signiﬁcance to the survival and growth of crops. Based on the hyperspectral (350–1350 nm) and LWC data (jointing, booting, ﬂowering, and ﬁlling periods) of winter wheat in 2020 and 2021, this work proposed to transform and process the hyperspectral data by adopting fractional order differential and continuous wavelet transform, and took a differential spectrum, wavelet coefﬁcients, and mixed variables (differential spectrum and wavelet coefﬁcients) as input variables of the model and adopted Gaussian process regression (GPR), classiﬁcation and regression decision tree (CART), and artiﬁcial neural network (ANN) methods to estimate the LWC of wheat in different growth periods. The results indicated that fractional differential and continuous wavelet transform could highlight the spectral characteristics of winter wheat canopy and improve its correlation with LWC. The three model variables had the best estimation effect on LWC in the ﬂowering period, and the average values of R 2 were 0.86 and 0.87 in modeling and veriﬁcation, which indicated that the ﬂowering period could be used as the best estimation period for LWC. Compared with the differential spectrum and wavelet coefﬁcients, LWC estimation based on mixed variables performed best. The average values of R 2 in modeling and veriﬁcation were 0.78 and 0.79. Among them, the ANN model had the highest estimation accuracy, and the R 2 in modeling and veriﬁcation could reach 0.92 and 0.91. This showed that fractional differential and continuous wavelet transform could effectively promote the sensitivity of spectral information to LWC and enhance the prediction ability and stability of wheat LWC. The outcomes of the present study have the potential to provide new ideas for the water monitoring of crops.


Introduction
Water is an essential component of crops, and water deficiency will affect the growth and development process and eventually affect the yield and quality.Thus, it is significant to obtain crop water content quickly and accurately for growth and development and yield improvement of crops [1,2].In the entire tissue of crops, crop leaves have the most vigorous metabolism.Leaf water content (LWC) is an important indicator to evaluate the water status of the plant [3], and the degree of crop water deficiency can be determined by the high or low leaf water content.
Traditional LWC measurement methods are time-consuming, laborious, and inefficient, making it difficult to achieve rapid and accurate detection.In recent years, with the wide application of hyperspectral remote sensing technology in crop quality and growth monitoring [4][5][6][7][8], the efficient and accurate estimation of crop LWC by hyperspectral remote sensing technology has attracted much attention from scholars.Currently, most studies construct spectral parameters based on raw spectra and then build crop LWC estimation models [9][10][11][12].
However, these methods cannot fully utilize the effective information in the spectra, thus limiting the further improvement of the accuracy of crop LWC estimation models.Fractional order differential can refine the crop detail information in the spectrum [13], which can further improve the estimation accuracy of the model.For example, Zhang Zhitao et al. [14] explored the possibility of combining the fractional order differentiation technique with the SVM-RF model to estimate the organic matter content of desert soils, and they found that the performance of the model was improved after fractional order differentiation treatment compared with integer order differentiation.Li Changchun et al. [13] used a fractional order differentiation technique to process hyperspectral data and combined it with a machine learning algorithm to establish an estimation model of winter wheat chlorophyll content.They found that the fractional order differentiation technique further improved the estimation effect of chlorophyll content.Wumuti Aishanjiang et al. [15] explored the feasibility of the fractional order differential technique and the continuous projection algorithm combined with a back-propagation neural network model to estimate wheat LWC.It was found that the fractional order differential enhanced the correlation between wheat LWC and hyperspectral data, helping to improve the estimation effect of the model.
Wavelet transform can decompose the spectral data for effective noise reduction and extract more spectral locations and feature parameters [16].For example, Blackburn et al. [17] employed wavelet transform to decompose leaf hyperspectral data to estimate leaf chlorophyll content, and the results indicated that wavelet-based methods outperform those based on untransformed spectra and spectral parameters, and the wavelet analysis method deserves further research to quantitatively extract information from spectral data.Lin et al. [18] adopted continuous wavelet transform combined with the partial least squares method to estimate copper content in chicory leaves, which effectively improved the estimation accuracy of the copper content in chicory leaves.Tan Xianming et al. [19] used continuous wavelet transform to decompose the corn canopy spectral data.Then, an estimation model of corn canopy chlorophyll density was established, and the results indicated that compared with the vegetation target and the free combination target of bands, using wavelet coefficient to estimate chlorophyll achieved a better effect.Wang Yancang et al. [20] constructed a winter wheat LWC inversion model based on the decomposition of spectral data information by continuous wavelet transform, and the inversion effect surpassed the spectral index.
At present, the research of crop LWC estimation based on spectral data is mostly based on single variable parameters such as vegetation index, wavelet coefficient, and spectral characteristics, and few of them adopt fractional order differential and continuous wavelet transform for comprehensive comparative analysis.Therefore, this work proposed to transform and process the hyperspectral data by adopting fractional order differential and continuous wavelet transform, and GPR, CART, and ANN methods were combined to build the LWC estimation model.The objective of this study was: (i) to analyze the correlation between differential spectra and wavelet coefficient of different growth periods and LWC; (ii) based on the differential spectrum and wavelet coefficient, evaluate the ability to construct LWC estimation models for different growth periods using the GPR, CART, and ANN methods; (iii) the ability of fractional differential and continuous wavelet transform to improve the monitoring accuracy of winter wheat leaf moisture information was studied.

Overview of the Study Area
The experimental area is located in the Winter Wheat Comprehensive Experimental Base, Shizhuang Village, Ningguo Town, Shanyang District, Jiaozuo City, Henan.The experimental area has a warm, temperate, sub-humid monsoon climate, with sufficient sunshine, mild climate, abundant precipitation, and an annual average temperature of 14.2−14.8• C. The soil type is Chao soil, with a pH of 7.1, organic matter content of 8.9 g/kg, total nitrogen content of 0.56 g/kg, alkali hydrolyzed nitrogen of 46.84 mg/kg, available phosphorus of 2.78 mg/kg, and available potassium of 112.00 mg/kg.Two wheat varieties, Jing 9843 and Zhong Mai 175, were set up in the study area, with four fertilizer gradients (0, 195, 390, and 585 kg/ha).Three water treatments were set, namely: rainfed (no irrigation), normal water (irrigation water volume is 192 mm), and twice normal water (irrigation water volume is 384 mm).There were three treatments, 16 plots per treatment, and 48 plots in total.The size of each plot is 8 m × 6 m.The field management conditions of each plot were the same, and the geographical location of the study area is shown in Figure 1.
Agronomy 2023, 13, 56 3 of 19 ing the GPR, CART, and ANN methods; (iii) the ability of fractional differential and continuous wavelet transform to improve the monitoring accuracy of winter wheat leaf moisture information was studied.

Overview of the Study Area
The experimental area is located in the Winter Wheat Comprehensive Experimental Base, Shizhuang Village, Ningguo Town, Shanyang District, Jiaozuo City, Henan.The experimental area has a warm, temperate, sub-humid monsoon climate, with sufficient sunshine, mild climate, abundant precipitation, and an annual average temperature of 14.2−14.8°C.The soil type is Chao soil, with a pH of 7.1, organic matter content of 8.9 g/kg, total nitrogen content of 0.56 g/kg, alkali hydrolyzed nitrogen of 46.84 mg/kg, available phosphorus of 2.78 mg/kg, and available potassium of 112.00 mg/kg.Two wheat varieties, Jing 9843 and Zhong Mai 175, were set up in the study area, with four fertilizer gradients (0, 195, 390, and 585 kg/ha).Three water treatments were set, namely: rainfed (no irrigation), normal water (irrigation water volume is 192 mm), and twice normal water (irrigation water volume is 384 mm).There were three treatments, 16 plots per treatment, and 48 plots in total.The size of each plot is 8 m×6 m.The field management conditions of each plot were the same, and the geographical location of the study area is shown in Figure 1.

Data Collection and Processing
The hyperspectral reflectance and LWC data of the wheat canopy were obtained in 48 plots in 2020 and 2021 separately, and the data of four key growth periods, namely, the jointing period, booting period, flowering period, and filling period, were collected.

Hyperspectral Data Collection and Processing
The American ASD Field Spec ® 4 Hi-Res (Analytica Spectra Devices., Inc, Langmont, USA) spectrometer was used to collect wheat canopy hyperspectral data.The wavelength range of the instrument was 350-2500 nm, and the resampling interval was 1 nm.The sampling interval included 1.4 nm (350-1000 nm) and 2 nm (1001-2500 nm).The specific time for data collection in clear and windless weather was from 10:30 to 14:30.The spectral data collected in this time period were more stable than those collected in other time periods.When collecting data, the probe was vertically downward with a field angle

Data Collection and Processing
The hyperspectral reflectance and LWC data of the wheat canopy were obtained in 48 plots in 2020 and 2021 separately, and the data of four key growth periods, namely, the jointing period, booting period, flowering period, and filling period, were collected.

Hyperspectral Data Collection and Processing
The American ASD Field Spec ®4 Hi-Res (Analytica Spectra Devices., Inc., Langmont, CO, USA) spectrometer was used to collect wheat canopy hyperspectral data.The wavelength range of the instrument was 350-2500 nm, and the resampling interval was 1 nm.The sampling interval included 1.4 nm (350-1000 nm) and 2 nm (1001-2500 nm).The specific time for data collection in clear and windless weather was from 10:30 to 14:30.The spectral data collected in this time period were more stable than those collected in other time periods.When collecting data, the probe was vertically downward with a field angle of 25 • , which was 30 to 40 cm above the winter wheat canopy to reduce the impact of the soil background.Moreover, to eliminate the impact of external light changes on spectral data acquisition, the BaSO 4 whiteboard calibration was performed every 10-20 min.Each plot randomly selected three sample points that could represent the overall growth of wheat.For each sample point, 10 spectral curves were collected, and the average reflectance of the three sample points was regarded as the spectral reflectance value of the plot.After data acquisition, the spectral data were preprocessed by the Savitzky-Golay filtering method.Moreover, to eliminate the influence of noise in the spectral data, hyperspectral reflectance data in the range of 350-1350 nm were selected for analysis and modeling.

LWC Data Measurement
The LWC of wheat was obtained by the drying method.Five representative samples of wheat leaves were randomly selected from each plot, and 20 fully expanded leaves were collected from the upper part of the wheat canopy, and the average LWC value of the 100 leaves was regarded as the LWC value of the plot.The collected leaves were sealed, and the fresh weight of the leaf samples was weighed by precision electronic scales (0.001 g) in the laboratory.Then, they were put into the oven (80 • C) and dried to constant weight, and the dry weight of the leaf samples was weighed; finally, the LWC was calculated using Equation (1): where FW represents the fresh weight of the leaf sample, and DW represents the dry weight of the leaf sample.

Fractional Order Differentiation
Fractional order differentiation is a further development of integer differentiation.It may decompose spectral data, deeply mine sensitive information, and is widely used in agricultural remote sensing [21,22].There are three brands of fractional differentiation.In this work, because the definition form of Grünwald-Letnikov is simple and understandable, it is used to transform and decompose hyperspectral data [23].The expression is as follows: where α is the order, λ is the spectral wavelength, Γ is the Gamma function, and n is the D-value of the bound limits of the differential equation.

Continuous Wavelet Transform
Continuous wavelet transform performs local refinement analysis of a signal based on time and spatial frequency, which can decompose a complex signal into a series of wavelet coefficients of different scales [24].In this work, the continuous wavelet transform is employed to decompose the spectral data.As a mathematically linear transformation method often used to process hyperspectral data [25][26][27], the continuous wavelet transform can better process and interpret spectral information.The calculation is shown in Equations ( 3) and (4): where λ is the spectral wavelength, a is the scale factor, b is the translation factor, f (λ) represents the hyperspectral reflectance, and ψ a,b represents the wavelet basis.

Gaussian Process Regression
Gaussian process regression (GPR) is a machine learning method that uses a GP prior to regression modeling with properties determined by the mean and covariance [28].GP is a collection of random variables all obeying joint gaussian distribution, which is essentially a multivariate gaussian distribution [29].GPR ensures that the prediction error falls within a certain range for the training data, and it can adjust the complexity of the model.Thus, it is suitable for dealing with nonlinear regression, small samples, high dimensionality, etc.

Classification and Regression Tree
A decision tree is a method based on classification and regression.When it is used in regression analysis, it usually refers to the classification and regression tree (CART) [30].The principle of regression trees is to divide the feature plane into several cells, with each division corresponding to a specific output.The model has a fast training speed and is good at obtaining nonlinear relationships in the dataset.

Artificial Neural Network
The artificial neural network (ANN) was proposed by Rumelhart and McClelland in 1986 [31].It is a multilayer feedforward neural network with backward propagation of errors, mimicking the human brain nervous system in processing complex information, so it has powerful self-learning and self-adaptive capabilities.ANN mainly contains an input layer, a hidden layer, and an output layer.The model variables are input by the input layer; then, the input attribute information is sent to the hidden layer by using the Sigmoid function; finally, the target data are output after the processing of the hidden layer.If the error between the target data and the actual data exceeds the set threshold, the error information is fed back to the hidden layer, and the data processing is continued until the error satisfies the requirement.Specifically, the model variables of the input layer include 10 fractional differential spectra, 10 wavelet energy coefficients, and mixed variables selected, and the output data of the output layer are LWC estimation data.
In each growth period, the data sets were divided at a ratio of 3:1, where 36 specimens were chosen to build the modeling dataset, and 12 specimens were chosen to build the validation dataset.The differential spectra and wavelet energy coefficients were used as independent variables with LWC as the dependent variable to build the winter wheat LWC estimation model.Each model regression framework was constructed and implemented in R language v4.0.2.

Correlation Analysis
The Pearson correlation coefficient (ρ) was adopted for correlation analysis.The value of the correlation coefficient is between −1 and 1; the larger the absolute value, the greater the correlation between the two variables; conversely, the closer the value is to 0, the lower the correlation between the two variables [32].

Model Performance Evaluation
The coefficient of determination (R 2 ), root mean squared error (RMSE), and standard root mean squared error (nRMSE) are chosen as the precision evaluation metrics of the model [33,34].

Statistical Characteristics of LWC in Wheat
The main statistical characteristics of wheat LWC in each period of different years are shown in Table 1.It can be seen from Table 1 that in 2020, the LWC change range at the jointing period was 71.04-79.89%,at the booting period was 71.85-79.48%,at the flowering period was 68.94-78.44%,and at the filling period was 29.85-71.80%.In 2021, the change range of LWC at the jointing period was 72.97-81.87%,at the booting period was 73.86-82.20%,at the flowering period was 60.53-75.27%,and at the filling period was 46.09-67.97%.On the whole, the average value of LWC in different years showed a trend of first increasing and then decreasing with the advance of the growth period.The average value of LWC in the booting period reached the maximum, and that in the filling period reached the minimum.In addition, the SD of LWC in different years increased gradually with the advance of fertility.The original hyperspectral data of four growth periods of winter wheat were processed by using the Grünwald-Letnikov form of the fractional differential method, where the order range was 0-2 and the step size was set to 0.1.The original spectra and differential spectra in the four periods were analyzed for correlation with the LWC, and the correlation coefficient is shown in Figures 2 and 3 (where red represents positive correlation, blue represents negative correlation, and the darker the same color is, the stronger the correlation between the spectral variables and the LWC).Meanwhile, 10 differential spectra with greater correlation with LWC were chosen, and their relevance matrix with LWC is presented in Figure 4.
46.09-67.97%.On the whole, the average value of LWC in different years showed a trend of first increasing and then decreasing with the advance of the growth period.The average value of LWC in the booting period reached the maximum, and that in the filling period reached the minimum.In addition, the SD of LWC in different years increased gradually with the advance of fertility.The original hyperspectral data of four growth periods of winter wheat were processed by using the Grünwald-Letnikov form of the fractional differential method, where the order range was 0-2 and the step size was set to 0.1.The original spectra and differential spectra in the four periods were analyzed for correlation with the LWC, and the correlation coefficient is shown in Figures 2 and 3 (where red represents positive correlation, blue represents negative correlation, and the darker the same color is, the stronger the correlation between the spectral variables and the LWC).Meanwhile, 10 differential spectra with greater correlation with LWC were chosen, and their relevance matrix with LWC is presented in Figure 4.The following observations can be obtained from the analysis of Figures 2, 3, and 4: 1.In the jointing period, the correlation between the original hyperspectral data and the LWC was analyzed.It could be seen from Figure 2a that there was a significant negative correlation with the LWC at a level of 0.01 within the wavelength from 350 to 721 nm.When the negative correlation was the strongest, the wavelength was 708 nm, and the correlation coefficient ρ was −0.54.The correlation between the differen- The following observations can be obtained from the analysis of Figures 2-4: 1.
In the jointing period, the correlation between the original hyperspectral data and the LWC was analyzed.It could be seen from Figure 2a that there was a significant negative correlation with the LWC at a level of 0.01 within the wavelength from 350 to 721 nm.When the negative correlation was the strongest, the wavelength was 708 nm, and the correlation coefficient ρ was −0.54.The correlation between the differential spectrum and the LWC was analyzed.It could be seen from Figure 3a that the absolute value of the |ρ| of each order of the differential spectrum and LWC was greater than 0.54.When the order was 1.0, the |ρ| reached a maximum of 0.69 at a wavelength of 926 nm.Except for orders 1.0 and 2.0, more than 347 differential spectral bands reached a significant level of 0.01.At the jointing period, 10 differential spectra were selected as follows: J1.0R926, J2.0R735, J1.0R1022, J2.0R736, J1.1R733, J1.0R785, J1.1R732, J1.2R720, J1.3R720, and J1.9R713, respectively.The correlation results were shown in Figure 4a.2.
In the booting period, the correlation between the original hyperspectral data and the LWC was analyzed.It could be seen from Figure 2b that there was a significant negative correlation with the LWC at a level of 0.01 within the wavelength from 353 to 735 nm and 1149 to 1350 nm.The negative correlation was the strongest at 355 nm, and the ρ value was −0.75.Then, the correlation between the differential spectrum and the LWC was analyzed.It could be seen from Figure 3b that the |ρ| value of each order of the differential spectrum and LWC were greater than 0.73.When the order was 1.1, the wavelength was at 433 nm, and the |ρ| value was up to a maximum of 0.80.Except for order 1.0 and order 2.0, more than 558 differential spectral bands reached a significance level of 0.01.At the booting period, 10 differential spectra were selected as follows: J1.1R433, J1.4R357, J1.8R432, J1.2R433, J1.5R357, J1.9R432, J1.7R432, J1.3R357, J1.6R432, and J0.9R418, respectively.The correlation results are shown in Figure 4b.

3.
In the flowering period, the correlation between the primitive hyperspectral data and the LWC was analyzed.According to Figure 2c, there was a significant negative correlation with the LWC at a level of 0.01 at wavelengths from 350 to 725 nm and 1316 to 1350 nm, and a significant positive correlation with the LWC at a level of 0.01 from 760 to 828 nm.When the correlation was strongest, the wavelength was 356 nm, and the ρ value was −0.74.Then, the correlation between the differential spectrum and the LWC was analyzed.According to Figure 3c, the |ρ| value of each order of the differential spectrum and LWC was greater than 0.72.When the order was 1.1, the |ρ| value was up to a maximum of 0.79 at a wavelength of 434 nm.Except for orders 1.0 and 2.0, more than 518 differential spectral bands reached a significance level of 0.01.At the jointing period, 10 differential spectra were selected as follows: J1.1R434, J0.9R374, J1.9R433, J1.2R409, J1.3R409, J1.6R351, J1.7R351, J0.8R374, J1.5R351, and J1.8R351, respectively.The correlation results are shown in Figure 4c.4.
In the filling period, the correlation between the primitive hyperspectral data and the LWC was analyzed.It could be seen from Figure 2d that there was a significant negative correlation with the LWC at a level of 0.01 at wavelengths from 350 to 725 nm and 1316 to 1350 nm, and a significant positive correlation with the LWC at a level of 0.01 from 736 to 1137 nm; the strongest correlation was obtained at 680 nm, with the ρ value of −0.76.Then, the correlation between the differential spectrum and the LWC was analyzed.As shown in Figure 3d, the maximum value of |ρ| was greater than 0.72.When the order was 1, the |ρ| value reached a maximum of 0.78 at a wavelength of 675 nm.Except for orders 1.0 and 2.0, more than 792 differential spectral bands achieved a significant level of 0.01.At the jointing period, 10 differential spectra were selected as follows: J1.0R675, J1.0R1299, J1.0R720, J1.1R693, J1.1R689, J1.0R561, J1.9R688, J1.2R693, J1.0R564, and J1.3R693, respectively.The correlation results are shown in Figure 4d.

Construction and Analysis of LWC Estimation Model
In each growth period, the 10 fractional differential spectrums screened were taken as independent variables and LWC as dependent variables.Then, 75% of the specimen data were randomly selected, and three methods, namely GPR, CART, and ANN, were used to estimate LWC for different growth periods of wheat.The remaining 25% of the specimen data were used for model accuracy verification.The results are presented in Table 2.The following observations can be obtained from Table 2: 1.
In the jointing period, the modeling accuracy of the GPR model was R 2 = 0.72, RMSE = 1.05%, and nRMSE = 1.38%, and the validation accuracy was R 2 = 0.75, RMSE = 0.97%, and nRMSE = 1.26%.Compared with the CART model and the ANN model, the modeling and validation R 2 of the GPR increased by 0.11 and 0.13, and 0.11 and 0.06, the RMSE decreased by 0.02% and 0.30%, and 0.01% and 0.01%, and the nRMSE decreased by 0.02% and 0.42%, and 0.01% and 0.02%, respectively.The GPR model achieved higher modeling and verification accuracy than the other two models.
A comprehensive analysis indicated that the GPR model had a better estimation effect in the jointing period.2.
In the booting period, the modeling accuracy of the GPR model was R 2 = 0.79, RMSE = 0.91%, and nRMSE = 1.19%, and the validation accuracy was R 2 = 0.84, RMSE = 0.83%, and nRMSE = 1.08%.Compared with the CART model and the ANN model, the modeling and validation R 2 of the GPR increased by 0.05 and 0.08, and 0.08 and 0.15, the RMSE decreased by 0.03% and 0.30%, and 0.09 and 0.19%, and the nRMSE decreased by 0.03% and 0.42%, and 0.12% and 0.27%, respectively.The GPR model achieved higher modeling and verification accuracy than the other two models.Through a comprehensive analysis, it was shown that the GPR model had a better estimation effect in the booting period.

3.
In the flowering period, the modeling accuracy of the ANN model was R 2 = 0.89, RMSE = 0.80%, and nRMSE = 1.09%, and the validation accuracy was R 2 = 0.88, RMSE = 1.14%, and nRMSE = 1.56%.Compared with the GPR model and the CART model, the modeling and validation R 2 of the ANN increased by 0.02 and 0.03, and 0.11 and 0.02, and the RMSE decreased by 0.23% and 0.10%, and 0.28% and 0.37%, and the nRMSE decreased by 0.31% and 0.12%, and 0.38% and 0.47%, respectively.The ANN model achieved higher modeling and verification accuracy than the other two models.Through a comprehensive analysis, it was found that the ANN model had the highest estimation accuracy in the flowering period.

4.
In the filling period, the modeling accuracy of the ANN model was R 2 = 0.85, RMSE = 3.82%, and nRMSE = 6.25%, and the validation accuracy was R 2 = 0.83, RMSE = 2.89%, and nRMSE = 4.61%.Compared with the GPR model and the CART model, the modeling and validation R 2 of the ANN increased by 0.00 and 0.05, and 0.13 and 0.02, and the RMSE decreased by 0.54% and 2.45%, and 1.32% and 1.49%, and the nRMSE decreased by 0.77% and 4.02%, and 2.13% and 2.21%, respectively.The ANN model achieved higher modeling and verification accuracy than the other two models.A comprehensive analysis indicated that the ANN model had a better estimation effect in the filling period.
When fractional order differential spectroscopy was combined with GPR, CART, and ANN to estimate LWC in different growth periods, the GPR achieved good estimation results in the jointing and panicle periods, and the ANN model had the highest estimation accuracy in the flowering and filling periods.

Correlation Analysis between Wavelet Coefficients and LWC
The original hyperspectral data of the four growth periods were decomposed at different scales through continuous wavelet transform to obtain wavelet coefficients at the decomposition scale from 1 to 10, in which the wavelet base was Mexican Hat, the second derivative of the Gaussian function.Then, the correlation between the wavelet coefficient and the LWC in the four periods was analyzed, and the correlation between the wavelet coefficient and the LWC in each growth period was obtained.The results are shown in Figure 5.Meanwhile, 10 wavelet coefficients with greater correlation were chosen to draw their correlation matrix with LWC, and the results are illustrated in Figure 6.

Correlation Analysis between Wavelet Coefficients and LWC
The original hyperspectral data of the four growth periods were decomposed ferent scales through continuous wavelet transform to obtain wavelet coefficient decomposition scale from 1 to 10, in which the wavelet base was Mexican Hat, the derivative of the Gaussian function.Then, the correlation between the wavelet coe and the LWC in the four periods was analyzed, and the correlation between the w coefficient and the LWC in each growth period was obtained.The results are sh Figure 5.Meanwhile, 10 wavelet coefficients with greater correlation were chosen t their correlation matrix with LWC, and the results are illustrated in Figure 6.The following observations can be obtained from the analysis of Figures 5 and 6: 1.In the jointing period, the correlation between the wavelet coefficient and the LWC was analyzed.It could be seen from Figure 5a that the correlation between the wavelet coefficient and the LWC was stronger at first and then weakened as the decomposition scale increased from 1 to 10. Within the decomposition scale of 1 to 8, the maximum sizes of the wavelet coefficients and the LWC at each scale |ρ| were greater than 0.60.When the scale was 3, the |ρ| value reached the maximum of 0.70 at a wavelength of 1142 nm.Then, the 10 wavelet coefficients with a strong correlation with LWC were chosen as the independent variables, and their decomposition scales  The following observations can be obtained from the analysis of Figures 5 and 6: 1.
In the jointing period, the correlation between the wavelet coefficient and the LWC was analyzed.It could be seen from Figure 5a that the correlation between the wavelet coefficient and the LWC was stronger at first and then weakened as the decomposition scale increased from 1 to 10. Within the decomposition scale of 1 to 8, the maximum sizes of the wavelet coefficients and the LWC at each scale |ρ| were greater than 0.60.When the scale was 3, the |ρ| value reached the maximum of 0.70 at a wavelength of 1142 nm.Then, the 10 wavelet coefficients with a strong correlation with LWC were chosen as the independent variables, and their decomposition scales and bands were C3R1142, C2R972, C1R734, C2R733, C3R733, C4R732, C6R732, C5R731, C3R731, and C2R918, respectively.The correlation results are shown in Figure 6a.2.
In the booting period, the correlation between the wavelet coefficient and the LWC was analyzed.According to Figure 5b, the correlation between the wavelet coefficient and the LWC first became stronger and then weakened as the decomposition scale increased from 1 to 10.The maximum values of |ρ| between the wavelet coefficient and the LWC within the decomposition scale of 1 to 8 were greater than 0.71.When the scale was 4, the |ρ| was up to 0.78 at 403 nm.The 10 wavelet coefficients with a strong correlation with LWC were chosen as the independent variables of the model, and their decomposition scales and bands were C4R403, C2R352, C4R396, C3R443, C2R353, C4R445, C6R1226, C4R464, C3R445, and C4R463, respectively.The correlation results are shown in Figure 6b.

3.
In the flowering period, the correlation between the wavelet coefficient and the LWC was analyzed.It could be observed from Figure 5c that the correlation between the wavelet coefficient and the LWC was first stronger and then weakened as the decomposition scale increased from 1 to 10. Within the decomposition scale of 1 to 8, the maximum sizes of the wavelet coefficients and LWC at each scale |ρ| were greater than 0.70.When the scale was 8, the |ρ| reached the maximum of 0.80 at the wavelength of 1321 nm.Then, the 10 wavelet coefficients with a strong correlation with LWC were selected as the independent variables of the model, and their decomposition scales and bands were C8R1321, C3R1326, C8R1328, C8R1309, C7R1288, C7R1282, C8R1307, C6R1223, C4R840, and C6R1216, respectively.The correlation results are shown in Figure 6c.4.
In the filling period, the correlation between the wavelet coefficient and the LWC was analyzed.According to Figure 5d, the correlation between the wavelet coefficient and the LWC was stronger at first and then weakened as the decomposition scale increased from 1 to 10.The maximum values of |ρ| between the wavelet coefficient and LWC at the decomposition scale within 1 to 8 were greater than 0.77.When the scale was 3, the |ρ| was up to 0.80 at a wavelength of 450 nm.Then, the 10 wavelet coefficients with a strong correlation with LWC were selected as the independent variables of the model, and their decomposition scales and bands were C3R450, C1R743, C6R764, C2R742, C3R742, C6R768, C4R746, C8R1347, C6R775, and C5R755, respectively.The correlation results are shown in Figure 6d.

Construction and Analysis of LWC Estimation Model
In each growth period, the 10 wavelet coefficients screened were taken as independent variables and LWC as dependent variables.Then, 75% of the specimen data were randomly selected to build LWC estimation models for different growth periods of wheat by using three methods, namely, GPR, CART, and ANN.The remaining 25% of the specimen data were used to verify the accuracy of the model.The results are presented in Table 3.The following observations are obtained from Table 3: 1.
In the jointing period, the modeling accuracy of the ANN model was R 2 = 0.72, RMSE = 0.86%, and nRMSE = 1.13%, and the validation accuracy was R 2 = 0.76, RMSE = 0.98%, and nRMSE = 1.28%.Compared with the GPR model and the CART model, the modeling and validation R 2 of the ANN increased by 0.10 and 0.13, and 0.11 and 0.15, and the RMSE decreased by 0.27% and 0.04%, 0.26% and 0.10%, and the nRMSE decreased by 0.36% and 0.07%, and 0.34% and 0.16%, respectively.The ANN model achieved higher modeling and verification accuracy than the other two models.Through a comprehensive analysis, it was found that the ANN model achieved a better estimation effect in the jointing period.2.
In the booting period, the modeling accuracy of the GPR model was R 2 = 0.78, RMSE = 1.00%, and nRMSE = 1.32%, and the validation accuracy was R 2 = 0.83, RMSE = 0.93%, and nRMSE = 1.23%.Compared with the CART model and the ANN model, the modeling and validation R 2 of the GPR increased by 0.08 and 0.08, and 0.09 and 0.10, and the RMSE decreased by 0.02% and 0.47%, and 0.00% and 0.03%, and the nRMSE decreased by 0.02% and 0.63%, and 0.00% and 0.03%, respectively.The modeling and validation accuracy of the GPR model was higher than that of the other two models.A comprehensive analysis indicated that the GPR model achieved a better estimation effect in the booting period.

3.
In the flowering period, the modeling accuracy of the ANN model was R 2 = 0.91, RMSE = 0.71%, nRMSE = 0.97%, and the validation accuracy was R 2 = 0.89, RMSE = 0.85%, and nRMSE = 1.16%.Compared with the GPR model and the CART model, the modeling and validation R 2 of the ANN increased by 0.03 and 0.01, and 0.16 and 0.12, and the RMSE decreased by 0.40% and 0.81%, and 0.38% and 1.02%, and the nRMSE decreased by 0.54% and 1.10%, and 0.51% and 1.36%, respectively.The modeling and validation accuracy of the ANN model was higher than that of the other two models.Through a comprehensive analysis, it was found that the ANN model had the highest estimation accuracy in the flowering period.4.
In the filling period, the modeling accuracy of the ANN model was R 2 = 0.85, RMSE = 3.57%, and nRMSE = 5.84%, and the validation accuracy was R 2 = 0.86, RMSE = 3.27%, and nRMSE = 5.27%.Compared with the GPR model and the CART model, the modeling and validation R 2 of the ANN increased by 0.01 and 0.01, and 0.11 and 0.09, and the RMSE decreased by 1.10% and 0.98%, and 1.36% and 4.34%, and the nRMSE decreased by 1.71% and 1.59%, and 2.27% and 7.65%, respectively.The modeling and validation accuracy of the ANN model was higher than that of the other two models.A comprehensive analysis indicated that the ANN model achieved a better estimation effect in the filling period.
When using wavelet coefficients combined with the GPR, CART, and ANN methods to estimate LWC in different growth periods, the ANN model at the joining, flowering, and filling periods achieved the highest estimation accuracy, and the GPR model obtained the best estimation effect in the booting period.

Estimation of Wheat LWC Based on Mixed Variables
In the four growth periods, the differential spectrum and wavelet coefficients screened previously were taken as model variables, and three methods of GPR, CART, and ANN were used to estimate the LWC.Table 4 shows the estimated results of LWC.
The following observations are obtained from Table 4: 1.
In the jointing period, the modeling accuracy of the GPR model was R 2 = 0.77, RMSE = 1.04%, and nRMSE = 1.37%, and the validation accuracy was R 2 = 0.77, RMSE = 0.93%, and nRMSE = 1.22%.Compared with the CART model and the ANN model, the modeling and validation R 2 of the GPR increased by 0.15 and 0.11, and 0.17 and 0.10, and the RMSE decreased by 0.06% and 0.29%, and 0.03% and 0.37%, and the nRMSE decreased by 0.08% and 0.40%, and 0.04% and 0.49%, respectively.The GPR model achieved higher modeling and verification accuracy than the other two models.Through a comprehensive analysis, it was found that the GPR model achieved a better estimation effect in the jointing period.2.
In the booting period, the modeling accuracy of the GPR model was R 2 = 0.83, RMSE = 0.86%, and nRMSE = 1.12%, and the validation accuracy was R 2 = 0.81, RMSE = 0.72%, and nRMSE = 0.94%.Compared with the CART model and the ANN model, the modeling and validation R 2 of the GPR increased by 0.08 and 0.06, and 0.10 and 0.07, and the RMSE decreased by 0.06% and 0.46%, and 0.07% and 0.19%, and the nRMSE decreased by 0.09% and 0.61%, and 0.09% and 0.25%, respectively.The modeling and validation accuracy of the GPR model was higher than that of the other two models.A comprehensive analysis indicated that the GPR model achieved a better estimation effect in the booting period.

3.
In the flowering period, the modeling accuracy of the ANN model was R 2 = 0.92, RMSE = 0.64%, nRMSE = 0.87%, and the validation accuracy was R 2 = 0.91, RMSE = 0.82%, and nRMSE = 1.12%.Compared with the GPR model and the CART model, the modeling and validation R 2 of the ANN increased by 0.02 and 0.02, and 0.09 and 0.05, and the RMSE decreased by 0.31% and 0.18%, and 0.30% and 0.68%, and the nRMSE decreased by 0.43% and 0.25%, and 0.41% and 0.90%, respectively.The modeling and validation accuracy of the ANN model was higher than that of the other two models.Through a comprehensive analysis, it was found that the ANN model had the highest estimation accuracy in the flowering period.4.
In the filling period, the modeling accuracy of the ANN model was R 2 = 0.86, RMSE = 3.14%, and nRMSE = 5.09%, and the validation accuracy was R 2 = 0.83, RMSE = 6.01%, and nRMSE = 9.57%.Compared with the GPR model and the CART model, the modeling and validation R 2 of the ANN increased by 0.01 and 0.04, and 0.10 and 0.08, and the RMSE decreased by 0.94% and 0.43%, and 1.68% and 1.84%, and the nRMSE decreased by 1.46% and 0.68%, and 2.81% and 4.14%, respectively.The modeling and validation accuracy of the ANN model was higher than that of the other two models.
A comprehensive analysis indicated that the ANN model achieved a better estimation effect in the filling period.When using mixed variables combined with GPR, CART, and ANN methods to estimate LWC in different growth periods, the GPR model obtained the best estimation ability in the joining and booting periods, the ANN model at the flowering and filling periods achieved the highest estimation accuracy.

Discussion
Analyzing the LWC estimation results of wheat in different periods, it could be seen that the overall LWC estimation effect was very good, the R 2 of modeling and verification reached 0.92 and 0.91, and the mean value reached 0.77 and 0.78, respectively.This is because the hyperspectral data have a high resolution and strong band continuity, which can correctly mirror the spectral characteristics and differences of crops and correctly obtain some agricultural information.Thus, the method is suitable for monitoring the physiological and biochemical parameters of crops.The difficulty of estimating crop LWC by adopting hyperspectral data lies in the influence of external interference factors when collecting data, which leads to noise in spectral data and affects the extraction of susceptible information [35].After processing with fractional differential and continuous wavelet transform, the influence of noise can be eliminated, and the sensitive messages in spectral data can be deeply mined, further improving the estimation accuracy of crop LWC.This is consistent with the research results in the literature [14,36,37].
The spectral characteristics of crops are closely related to their growth periods, health status, and external environment, and the spectral curves of crops in different growth periods will show varying change characteristics, so the estimation models constructed with data of one growth period generally have poor universality [38,39].Therefore, this paper collected the experimental data from four key periods of jointing, booting, flowering, and filling for modeling.By analyzing the LWC estimation results of the four growth periods, it could be seen that in the four periods, the estimation ability of the jointing period was relatively poor, and the mean value of R 2 in modeling and verification was only 0.66 and 0.70, respectively.In the booting period, the mean values of R 2 were 0.75 and 0.77, respectively, which was slightly better than that in the jointing period.In the flowering period, the average values of R 2 for modeling and verification were 0.86 and 0.87, and those in the filling period were 0.81 and 0.81, respectively.It showed that the LWC estimation accuracy in the flowering period was higher than that in the filling period.This may be because the LWC of wheat in the flowering period gradually decreases and changes greatly, and the sensitivity of LWC to hyperspectral data is relatively strong.This is consistent with the research results in the literature [20,40].
In the four growth periods, a comprehensive analysis of the LWC estimation results of using differential spectrum and wavelet coefficient combined with GPR, CART, and ANN models indicated that compared with the differential spectrum and wavelet coefficient, the LWC estimation using combined mixed variables contributed to better results.The maximum R 2 of modeling and testing could reach 0.92 and 0.91, and the average value could reach 0.78 and 0.79.This is because fractional differential and wavelet transform can eliminate the influence of noise in spectral information, improve the sensitivity of LWC to spectral data, and further enhance the stability and robustness of the model.This is compatible with the study results in the literature [13,16,41].
By comparing and analyzing the LWC estimation results based on different methods in different growth periods, it was found that in the jointing period, the mixed variables combined with the ANN model achieved a better effect in estimating LWC, and the R 2 of modeling and testing was 0.77 and 0.77, respectively.In the booting period, the mixed variables combined with the GPR model performed better in estimating LWC, and the R 2 of modeling and testing was 0.83 and 0.81, respectively.In the flowering period, using the mixed variables combined with the ANN model to estimate LWC obtained the highest accuracy, the R 2 of modeling and testing was 0.92 and 0.91, respectively.In the filling period, using the wavelet coefficient combined with the ANN model to estimate LWC contributed to high accuracy, and the R 2 of modeling and testing was 0.85 and 0.86, respectively.This work proposed to transform and process the hyperspectral data by adopting fractional order differential and continuous wavelet transform, which can effectively abate the noise impact in spectral data and improve the estimation effect of wheat LWC.The accurate monitoring of LWC can effectively diagnose wheat water status, provide basic information for wheat growth monitoring and drought monitoring, and is of great signifi-cance for achieving a stable and high yield of wheat [24].The research results can provide a reference for LWC based on remote sensing technology.In addition, fractional differential and continuous wavelet transform can also be applied to other crop phenotypic traits (such as biomass, chlorophyll, and nitrogen content).

Conclusions
Based on two data processing methods of fractional order differential and continuous wavelet transform, this work processed the original hyperspectral data of winter wheat in four periods.Then, the original hyperspectral data, differential spectrum, and wavelet coefficients were analyzed, and the analysis results were compared to obtain 10 differential spectra and 10 wavelet coefficients that were more closely related to LWC.Subsequently, the LWC estimation models for different growth periods of winter wheat were built by using GPR, CART, and ANN methods.Through analyzing the modeling and verification results of the model, the optimal LWC estimation model for each growth period and the growth period with the highest precision of LWC estimation were determined, providing new ideas for estimating crop LWC.Although this work has obtained some research results, several defects need to be addressed.
When a machine learning algorithm is used for modeling, a large number of samples are required.The work used a small amount of sample data, resulting in overfitting, which reduced the stability of the model.In future studies, the number of samples, sample years, wheat varieties, etc. can be increased to verify the model, thus further improving the stability and universality of the model.Meanwhile, because of the constriction of experimental factors, the work used data of a piece of farmland, without considering the influence of different areas on the experimental data.In the future, it is necessary to take experimental data from multiple areas and use more sample experimental data to train models, thus achieving better estimation effects and stronger robustness.

Figure 1 .
Figure 1.The schematic of the study area.

Figure 1 .
Figure 1.The schematic of the study area.

Figure 2 .
Figure 2. Correlation analysis of the raw hyperspectral with LWC in four growth periods.(a) Jointing period; (b) booting period; (c) flowering period; (d) filling period.Figure 2. Correlation analysis of the raw hyperspectral with LWC in four growth periods.(a) Jointing period; (b) booting period; (c) flowering period; (d) filling period.

Figure 4 .
Figure 4.The relevance matrix of the selected fractional order differential spectra with LWC in four growth periods.(a) Jointing period; (b) booting period; (c) flowering period; (d) filling period.

Figure 4 .
Figure 4.The relevance matrix of the selected fractional order differential spectra with LWC in four growth periods.(a) Jointing period; (b) booting period; (c) flowering period; (d) filling period.

Figure 5 .
Figure 5.The relevance between the wavelet coefficient of different growth periods and LW Jointing period; (b) booting period; (c) flowering period; (d) filling period.Figure 5.The relevance between the wavelet coefficient of different growth periods and LWC.(a) Jointing period; (b) booting period; (c) flowering period; (d) filling period.

Figure 5 .
Figure 5.The relevance between the wavelet coefficient of different growth periods and LW Jointing period; (b) booting period; (c) flowering period; (d) filling period.Figure 5.The relevance between the wavelet coefficient of different growth periods and LWC.(a) Jointing period; (b) booting period; (c) flowering period; (d) filling period.

Figure 6 .
Figure 6.The relevance matrix of the selected wavelet coefficient with LWC in four growth periods.(a) Jointing period; (b) booting period; (c) flowering period; (d) filling period.
and bands were C3R1142, C2R972, C1R734, C2R733, C3R733, C4R732, C6R732, C5R731, C3R731, and C2R918, respectively.The correlation results are shown in Figure 6a.2. In the booting period, the correlation between the wavelet coefficient and the LWC was analyzed.According to Figure 5b, the correlation between the wavelet coefficient and the LWC first became stronger and then weakened as the decomposition scale increased from 1 to 10.The maximum values of |ρ| between the wavelet coefficient and the LWC within the decomposition scale of 1 to 8 were greater than 0.71.When

Figure 6 .
Figure 6.The relevance matrix of the selected wavelet coefficient with LWC in four growth periods.(a) Jointing period; (b) booting period; (c) flowering period; (d) filling period.

Table 1 .
Statistical characteristics of LWC from each growth period in different years..2. Estimation of Wheat LWC Based on Fractional Order Differential 3.2.1.Analysis of Correlation between Original Hyperspectral and Fractional Order Differential Spectrum and LWC MAX means maximum value; MIN means minimum value; SD means standard deviation.3

Table 1 .
Statistical characteristics of LWC from each growth period in different years.

Table 2 .
The LWC estimation results obtained by adopting differential spectra combined with GPR, CART, and ANN for different growth periods.

Table 3 .
The LWC estimation results obtained by adopting wavelet coefficients combined with GPR, CART, and ANN for different growth periods.

Table 4 .
The LWC estimation results obtained by using differential spectrum and wavelet coefficients combined with GPR, CART, and ANN for different growth periods.