Next Article in Journal
SF-CNN: Signal Filtering Convolutional Neural Network for Precipitation Intensity Estimation
Next Article in Special Issue
Web and MATLAB-Based Platform for UAV Flight Management and Multispectral Image Processing
Previous Article in Journal
Eye-Tracker Study of Influence of Affective Disruptive Content on User’s Visual Attention and Emotional State
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Winter Wheat Nitrogen Estimation Based on Ground-Level and UAV-Mounted Sensors

1
Information Technology Research Center, Beijing Academy of Agriculture and Forestry Sciences, Beijing 100097, China
2
Key Laboratory of Quantitative Remote Sensing in Agriculture of Ministry of Agriculture, Beijing 100097, China
3
Anhui Engineering Laboratory of Agro-Ecological Big Data, Anhui University, Hefei 230601, China
4
Aerial Application Technology Research Unit, USDA-Agricultural Research Service, College Station, TX 77845, USA
5
National Engineering Research Center for Information Technology in Agriculture, Beijing 100097, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(2), 549; https://doi.org/10.3390/s22020549
Submission received: 13 December 2021 / Revised: 6 January 2022 / Accepted: 7 January 2022 / Published: 11 January 2022
(This article belongs to the Special Issue Application of UAV and Sensing in Precision Agriculture)

Abstract

:
A better understanding of wheat nitrogen status is important for improving N fertilizer management in precision farming. In this study, four different sensors were evaluated for their ability to estimate winter wheat nitrogen. A Gaussian process regression (GPR) method with the sequential backward feature removal (SBBR) routine was used to identify the best combinations of vegetation indices (VIs) sensitive to wheat N indicators for different sensors. Wheat leaf N concentration (LNC), plant N concentration (PNC), and the nutrition index (NNI) were estimated by the VIs through parametric regression (PR), multivariable linear regression (MLR), and Gaussian process regression (GPR). The study results reveal that the optical fluorescence sensor provides more accurate estimates of winter wheat N status at a low-canopy coverage condition. The Dualex Nitrogen Balance Index (NBI) is the best leaf-level indicator for wheat LNC, PNC and NNI at the early wheat growth stage. At the early growth stage, Multiplex indices are the best canopy-level indicators for LNC, PNC, and NNI. At the late growth stage, ASD VIs provide accurate estimates for wheat N indicators. This study also reveals that the GPR with SBBR analysis method provides more accurate estimates of winter wheat LNC, PNC, and NNI, with the best VI combinations for these sensors across the different winter wheat growth stages, compared with the MLR and PR methods.

1. Introduction

Nitrogen (N) is a crucial nutrient required for crop growth and grain formation. Agricultural managers can regulate N management at suitable rates and opportune moments based on the crop’s N requirements. As an essential scatheless and real-time technique, remote sensing technologies have been proved valuable for crop N status evaluation [1,2,3,4,5,6].
Much meaningful progress in sensor technology for evaluation of plant N status has been achieved in recent years. Leaf sensors, such as chlorophyll meters [7,8,9] and Dualex sensors [10], have been widely used to measure crop N status. These leaf clip sensors show a stable relation with plant N due to their direct contact with the plant. However, there are also some restrictions on these sensors because plant leaf water content, leaf structure, or other nutrient deficiencies may also easily influence the sensor readings [11]. Some studies indicate that measurement of polyphenol concentrations in the leaf is a new method to overcome such barriers [12]. The ratio of chlorophyll to polyphenol is more stable than the leaf chlorophyll distribution [11,13]. At the canopy level, the spectral features of chlorophyll from the visible to near-infrared bands are used as indicators of crop N [14,15,16]. The Multiplex 3 portable sensor has also been widely used for plant N diagnosis in recent years, and it can detect plant chlorophyll and flavonol compounds simultaneously through the chlorophyll fluorescence method [17]. As a new exploration of the low-altitude remote sensing method, unmanned aerial vehicle (UAV)-based sensing has been widely used in different fields recently due to its flexibility, affordability, and applicability for large-scale monitoring compared to handheld active sensing [18,19]. Digital color images acquired by UAV or airplane-mounted sensors have also proven to be a feasible way to estimate canopy variables such as leaf chlorophyll content [20], nitrogen status [21], wheat senescence [22], and vegetation cover [23] in a large area [24]. In general, rapid, precise, and non-destructive acquisition of N information has become an essential technique for crop nutrition and growth diagnosis [5,25], which helps dynamic regulations of N fertilizer use [26].
Previous studies attempted to retrieve N via a radiative transfer model. As an example, Jacquemoud et al. attempted to incorporate N into the PROSPECT model, but this approach was abandoned because of inconsistencies in N retrieval through model inversion [27]. More efforts have been directed at estimating crop N status using empirical methods based on observed reflectance to measured vegetation characteristics [28,29]. Many practical regression techniques using hyperspectral bands, vegetation indices (VIs), and different types of sensor data have been proved influential in plant N estimation [21,28,29]. Parametric regression (PR) methods based on band information or VIs of broadband satellite sensors are probably the oldest and largest group of variable estimation approaches [30]. VIs enhance spectral features sensitive to a vegetation property while reducing disturbances by combining some spectral bands [31,32]. The major advantage of VIs is their intrinsic simplicity. However, the selection of an optimal subset of hyperspectral bands or best VIs sensitive to plant N encounters both numerical and computational difficulties [33,34]. Multivariable linear regression (MLR) methods are attractive because of their fast performance in coping with spectroscopic data and their typical reliance on the estimation of covariances [30]. Miphokasap et al. demonstrated that the model developed by MLR led to a higher correlation coefficient and lower errors than model applications based on narrowband VIs in estimating canopy nitrogen [30,35]. In recent decades, a variety of non-linear, non-parametric methods have been developed, going beyond linear regression or linear transformation techniques. These methods, also referred to as machine learning regression algorithms, apply non-linear transformations [30]. Gaussian process regression (GPR) applied to spectroscopic and hyperspectral data started more recently, such as airborne HyMap mapping of leaf chlorophyll content [36] and spaceborne CHRIS mapping of leaf chlorophyll content, LAI, and fractional vegetation content [37]. The GPR method is deemed to be one of the most exciting machine learning regression algorithms, which can provide a full conditional statistic for the predicted variable [38,39]. GPR with the sequential backward feature removal (SBBR) routine proposed by Verrelst et al. leads to the identification of the most sensitive vegetation index for the crop nitrogen variable [34].
Nevertheless, when this technique is applied in studies with spectral data from different sensors, results on the best VIs or combinations and their performance on crop nitrogen diagnosis have rarely been reported. Therefore, in this study, we propose assessing the performance of different sensors for winter wheat N status estimation through PR, MLR, and GPR analyses. To determine which spectral features could sufficiently estimate crop N status indicators, we considered one leaf N status indicator (i.e., the leaf N concentration (LNC)), one plant N status indicator (i.e., the N plant concentration (PNC)), and one relative plant N indicator (i.e., the nutrition index (NNI)), which is calculated as the ratio of the measured N concentration and the critical N concentration. In this study, we also performed sensitivity analysis on the wheat N indicators for different sensors through the PR, MLR, and GPR methods. The main objectives of this study were to: (1) evaluate the ability of four different proximal and UAV-mounted sensors for winter wheat N status estimation; (2) identify the optimal VIs and VI combinations derived from the sensors for accurate winter wheat N status estimation; and (3) determine whether GPR models based on the combination of VIs could further improve the estimation accuracy.

2. Materials and Methods

2.1. Experimental Design

This study was carried out at the National Experimental Station for Precision Agriculture in the Changping District of Beijing, China (40°10.6′ N, 116°26.3′ E; Figure 1). The area of the study field was about 0.80 ha. Winter wheat was sown on 4 October 2013, in 15 cm-wide rows. The variety of wheat was Jindong 22. A total of 8 different nitrogen fertilizer treatments were assigned to 104 plots in this experiment, each plot measuring 7 × 7.5 m (Figure 1). There was a 1 m buffer around each fertilization plot to minimize the interference between different fertilizer treatments. Base fertilizer, including 72 kg/ha N, 60 kg/ha P2O5, and 75 kg/ha K2O, was applied to the experiment plots when the wheat was seeded, except for the CK treatments. Topdressing was applied at the wheat jointing stage (16 April, corresponding to Feekes 6), and detailed fertilizer information about the treatments is listed in Table 1.

2.2. Sensor Data Collection

The leaf and canopy spectral parameters were measured in a 1 m2 rectangular region within each plot at the wheat raising stage (8 April 2014, corresponding to Feekes 5) and the filling stage (28 May 2014, corresponding to Feekes 11). In this study, we used four different optical instruments to collect wheat spectral data. The Dualex sensor was used to measure wheat leaf spectral parameters. The Multiplex sensor and the ASD field spectrometer were used to collect wheat canopy spectral parameters. A Sony DSC-QX100 digital RGB camera mounted on a DJI S1000 UAV was used to collect wheat images for the experiment fields. Table 2 lists the detailed information on the sensors used in this study.

2.2.1. Leaf- and Canopy-Level Data Collection

For each plot, 10 wheat plants were randomly selected for Dualex measurement. The first and second fully expanded leaves of the plants were clipped at the middle location, and the measurement averages represented the Duelax NBI, Chl, and FLAV values of the plot. Multiplex measurements and ASD spectral measurements were performed and repeated 20 times at 10 cm and 1.3 m above the wheat canopy, respectively. The measurement averages represented the wheat canopy fluorescence and reflectance values for each plot. A total of 9 Multiplex fluorescence VIs and 28 ASD spectral VIs sensitive to plant N are listed in Table 3.

2.2.2. UAV-Level Data Collection

The UAV digital images were obtained at a 50 m altitude. Pix4Dmapper software (Pix4D Inc., Lausanne, Switzerland) was used for automatic image mosaicking. The original digital number (DN) values of the images were calibrated by the experimental line method (Equation (1)) [63]:
DN i = DN object DN black DN white DN black × 255
where DNi represents the calibrated DN value for band i, such as R, G, and B; DNobject is the original DN value of the images; and DNwhite and DNblack are the original DN values from the white and black panels in the UAV-DN images. Then, the images were resampled to 2 cm using the nearest neighbor method. An area of interest (AOI) (2 × 2 m) located in the middle region of each sample plot was delineated on the image to exclude the border effect. A total of 208 winter wheat AOIs (104 AOIs at wheat Feekes stage 5 and 104 AOIs at wheat Feekes stage 11) were clipped for the subsequent analysis. Then, the spectral data from the mix of winter wheat and soil background were extracted, and the VIs were calculated to indicate the winter wheat growth status for each plot. Table 4 lists the RGB image bands and VIs used in this study.

2.3. Plant Sampling Data Collection

We collected the plant samples as soon as the sensors’ measurements were completed. In this study, plants from an area of 0.12 m2 (two rows and 40 cm wide) near the spectral measurement plot were collected and sent to the laboratory for analysis. The vegetative organs (leaves, stems, and ears) were then divided and dried in the oven at 80 °C for 24 h. The biomass (g/m2) for wheat leaves (leaf biomass, LB), stems (stem biomass, SB), and ears (ear biomass, EB) per unit area was calculated based on the measured planting density and the dry weight of the samples. The leaf nitrogen concentration (LNC % (g N 100 g−1 DM)), stem nitrogen concentration (SNC % (g N 100 g−1 DM)), and ear nitrogen concentration (ENC % (g N 100 g−1 DM)) were then determined using the standard Kjeldahl method [67]. Plant N accumulation (PNA, kg/ha) was the sum of leaf, stem, and ear N accumulation (Equation (2)). Wheat plant N content (PNC %) was determined by PNA (kg/ha) divided by the plant biomass (kg/ha) (Equation (3)) [68].
PNA   ( kg / ha ) = ( LB × LNC )   + ( SB × SNC ) + ( EB × ENC )
PNC   ( % ) = PNA / ( LB + SB + EB )
The NNI was used to indicate the N status for each experimental plot, independent of the growth stage and differing biomass levels. The NNI was calculated by Equation (4) [69]:
NNI   = N a c t / N c  
where N a c t is the actual N concentration, and Nc is the critical N concentration as a percentage of the aboveground dry matter expressed in % DM. The N c for winter wheat was described by Equation (5) [70]:
N c = 4.15 × W 0.38
where W is the aboveground biomass, expressed in Mg dry matter (DM) ha−1.

2.4. Data Analysis Methods

Five traditional PR methods (linear, logarithmic, power, exponential, and second-order polynomial) and the MLR and GPR methods were used to establish relationships between wheat N parameters and spectral VIs in this study. For the PR methods, 70% of the samples were used to build the models, and 30% were used to validate them. A stepwise MLR was used to estimate the winter wheat N status by recursively applying multiple regression.
Compared to other machine learning regression approaches, the parameter optimization of GPR is simpler and can be automatically completed by maximizing the marginal likelihood in the training set [34]. The core of the kernel method for GPR is the squared exponential. The use of a flexible kernel function (covariance function) generally suffices for tackling most regression problems and is beneficial if prior knowledge is weak [37]. This is an opportunity to exploit asymmetries in the feature space by including a parameter per feature, as in the very common anisotropic squared exponential (SE) kernel function:
K ( x i , x j ) = γ e x p ( b = 1 B ( x i b x j b ) 2 2 σ b 2 )
where γ is a scaling factor; σ b   is the length scale per input feature, b = 1, …, B; and x i b represents the bth feature of the input feature vector x i . A smaller value of σ b   indicates a higher level of informative content of this certain image feature. Verrelst et al. [34] integrated this property and SBBR as a new feature selection algorithm, in which the least significant feature with the highest σ b   was removed at each iteration and a new GPR model was retrained with the remaining features only. A feature combination would eventually be identified by the feature selection algorithm. To enable the automated identification of the best-performing features for any regression problem, this feature selection algorithm was automated and integrated into a user-friendly tool named GPR-BAT [34].
In this study, GPR with SBBR was used firstly to identify the best VI combinations for different N parameters, and then MLR and GPR were used to build N parameter estimation models based on those VI combinations. The 10-fold cross-validation technique was used to determine the optimal number of latent factors based on the lowest root mean square error (RMSE). The accuracy of each model was evaluated by the coefficient of determination R2, RMSE, the mean absolute error (MAE), and the Nash–Sutcliffe modeling efficiency (NSE). The equations used to calculate RMSE, MAE, and NSE are presented in Table 5.

3. Results

3.1. Variation in Winter Wheat N Indicators

Table 6 describes the variation in the winter wheat N indicators for different growth stages. PNC had the most variation (CV = 49.46% across growth stages), followed by LNC (CV = 32.34%) and NNI (CV = 18.59%). At Feekes stage 5, the greatest difference occurred in NNI, which ranged from 0.74 to 1.58. PNC had more variation at Feekes stage 5 (CV = 13.42%) than at Feekes stages 11 (CV = 8.99%). The CV for LNC increased from 11.30% to 13.42% from Feekes 5 to Feekes 11. The CV for NNI decreased from 17.17% at Feekes 5 to 11.97% at Feekes 11.

3.2. Relationships between Wheat N and VIs for Different Sensors

After the correlation analysis between the VIs and the three N parameters, the VI with the highest correlation coefficient was selected to establish the N estimation model. Table 7 lists the VIs and their R2 and RMSE for the best models for wheat LNC, PNC, and NNI through traditional PR methods. For the PR methods, 70% of the samples (70 samples) were used to build the models, and 30% (34 samples) were used to validate them. It could be observed, at wheat Feekes stage 5, that the Dualex- and Multiplex-based VIs had a stronger correlation with wheat LNC, PNC, and NNI. The R2 for the LNC and PNC models based on the Dualex NBI reached 0.80 and 0.79, respectively. At the late wheat growth stage, canopy-level VIs performed better than the leaf-level VIs. The RGB and ASD VIs performed better than the Dualex and Multiplex VIs for LNC, PNC, and NNI estimation. RGB camera-based r performed best for LNC estimation at wheat Feekes stage 11, followed by ASD mND705, and Dualex CHI. For PNC estimation, ASD NDRE performed best, followed by the RGB camera-based ExR. Multiplex performed best on the wheat N status estimation across wheat stages (Feekes 5–11), and the R2 reached 0.87 for BRR_FRF and LNC, 0.86 for BRR_FRF and PNC, and 0.56 for NBI_G and NNI.

3.3. Wheat N Estimation through UAV-Level VIs

GPR and MLR were used to estimate the relationships between the RGB camera VIs and N variables. The 10-fold cross-validation technique was used to determine the optimal number of latent factors based on the lowest root mean square error (RMSE). Table 8 lists the cross-validation results for models built by the selected VI combinations and all VIs through the GPR and MLR methods. The VI combinations for different N variables were selected by the GPR-SBBR feature selection algorithm.
It can be observed in Table 8 that the best VI combinations for LNC, PNC, and NNI estimation at the different wheat growth stages are different. At Feekes stags 5, the VI combinations selected by the GPR-SBBR algorithm for LNC, PNC, and NNI estimation were all B and b, indicating that the RGB camera VIs B and b can be used to identify the N status at the wheat growth stage.
GPR models built by the selected VIs had a higher cross-validation R2 compared to the GPR models built by all VIs. MLR models based on the selected VIs achieved a higher R2 at Feekes stage 5 and Feekes stage 11 for LNC and NNI estimation. However, MLR models built by all VIs performed better than MLR models built by the selected VIs on N status estimation across different wheat growth stages (Feekes stage 5–11).

3.4. Wheat N Estimation through Canopy-Level VIs

Table 9 and Table 10 list the cross-validation results for the N estimation models through GPR and MLR based on canopy-level VIs. Table 9 indicates that for LNC estimation, the best ASD VI combinations were TBI1, NDIopt, and TBI2 at wheat Feekes stage 5, SR(700,670), SR(418,405), and SR(740,720) at Feekes stage 11, and SR(418,405), TBI1, and PPR across Feekes 5–11. For PNC estimation, the optimal VI combinations were TBI1, NDIopt, and TBI2 at wheat Feekes stage 5, SR(700,670), SR(418,405), and SR(740,720) at Feekes stage 11, and TBI1, PPR, REV, and REFD across Feekes 5–11. The best VIs for NNI estimation included two VIs (MSAVI and PPR) at wheat Feekes stage 5, four VIs (NDWI, MCARI, MSAVI, and PVR) at wheat Feekes stage 11, and five VIs (SR(418,405), NDIopt, MSAVI, MCARI1, and PPR) at wheat Feekes stage 5–11. It can be seen in Table 9 that the LNC, PNC, and NNI estimation models built by GPR with the selected VIs performed better than those built by the other three methods. MLR with all ASD VIs had a higher cross-validation R2 than MLR with the selected ASD VIs or GPR with all ASD VIs at wheat Feekes stage 5 and Feekes stage 11.
Table 10 lists the LNC, PNC, and NNI estimation results through both the GPR method and the MLR method. The best-performing Multiplex VI combinations for LNC estimation were SFR_G, FLAV, and NBI_R at wheat Feekes stage 5, NBI_R and SFR_R at wheat Feekes stage 11, and SFR_R, BRR_FRF, and NBI_R at wheat Feekes stage 5–11. The best Multiplex VI combinations for PNC were SFR_R, FLAV, and NBI_R at wheat Feekes stage 5, SFR_R at wheat Feekes stage 11, and SFR_R, BRR_FRF, and NBI_G at wheat Feekes stage 5–11. For NNI, Multiplex FLAV and FER_RG at wheat Feekes growth stage 5, SFR_G and SFR_R at wheat Feekes growth stage 11, and SFR_G, SFR_R, and NBI_G at wheat Feekes stage 5–11 performed well. The cross-validation R2 values for LNC, PNC, and NNI estimated by GPR with the selected VIs were higher than those estimated by GPR or MLR with all VIs at wheat Feekes stage 5 and Feekes stage 11.

3.5. Wheat N Estimation through Leaf-Level VIs

Table 11 lists the cross-validation results for the three N parameters estimated through the GPR and MLR methods based on Dualex VIs. It can be seen in Table 11 that the best-performing VI for LNC, PNC, and NNI was NBI at wheat Feekes stage 5. At wheat Feekes stage 11, Chl performed best for LNC, PNC, and NNI estimation. These results indicate that the Dualex sensor can detect the leaf N status at different wheat growth stages, especially at Feekes stage 5. During Feekes stage 5–11, the combination of the Dualex NBI and Chl performed best for estimating LNC, PNC, and NNI.

4. Discussion

Ground-based spectrometers are considered capable of monitoring crop trait expressions and nutrients [74,75]; however, data collection with handheld sensors has the limitations of instability, low sensing efficiency, and high cost [76]. In this study, variable-rate N fertilization treatments in winter wheat were carried out to represent the reality, which resulted in different winter wheat N nutrition conditions (Table 1 and Table 6). Four different sensors were used to obtain crop growth information during two winter wheat critical growth stages. This allowed a preliminary assessment of the utility of the handheld leaf and canopy sensors as well as the UAV-mounted sensor in the estimation of N status indictors of winter wheat. The predictive performances of the four sensors were compared using different modeling methods (PR, MLR, and GPR).

4.1. N Estimation Comparison for Different Sensors

We investigated the wheat N estimation capability for different sensors by comparing the cross-validation R2 of the models. Figure 2 clearly shows that different sensors have different capabilities in N status detection during the wheat growth stage.
It can be observed in Figure 2a,b that the Dualex sensor had an obvious advantage in LNC and PNC estimation at wheat Feekes stage 5, and the R2 for different modeling methods varied from 0.77 to 0.80 for LNC, and from 0.78 to 0.79 for PNC. The Multiplex sensor can also detect LNC and PNC better than the ASD spectrometer or the RGB camera at the early wheat growth stage because it can eliminate erroneous signals from bare soil or distinguish between different N treatments in shadow with full sunlight. At Feekes stage 11, due to the biomass saturation, the correlations between N status and the VIs of the Dualex and Multiplex sensors decreased quickly. The modeling mean R2 decreased from 0.79 to 0.33 for Dualex, and from 0.58 to 0.32 for Multiplex, over the winter wheat growth stages. Compared with the other sensors, the ASD spectrometer performed well in LNC and PNC estimation through GPR with the selected VIs at winter wheat Feekes stage 11. The estimation R2 reached 0.51 for LNC and 0.49 for PNC. For NNI detection, the Dualex sensor worked relatively well at wheat Feekes stage 5, with a mean R2 of 0.50 for different regression methods, followed by the Multiplex sensor, with a mean R2 of 0.48 for different regression methods. All sensors did not perform well in NNI estimation at wheat Feekes stage 11 in this study. When data from the two stages were combined (Feekes 5–11), the modeling accuracy for the three N parameters all increased obviously. It can be observed in Figure 2 that the ASD and Multiplex sensors can detect LNC and PNC better than the RGB camera and the Dualex sensor. The R2 values are all greater than 0.90 for the ASD and Multiplex sensors.
Numerous studies have revealed that, at the leaf level, there is a well-established correlation between photosynthetic capacity and N content [27], while at the scale of canopies, reflectance patterns represent the integrated effects of leaf water content, biochemical constituents, and various components of the plant structure [77]. As reviewed by Tremblay et al., determining N status in plants by means of chlorophyll fluorescence can overcome some of the limitations of reflectance-based chlorophyll methods [10]. The Dualex VIs, namely, the Nitrogen Balance Index (NBI) and the chlorophyll index (CHI), seem to be good indicators to evaluate the N conditions in wheat at the early growth stage. This was also confirmed by the results of this study. The Multiplex sensor can overcome the deficiency of a chlorophyll meter because it is able to distinguish N treatments equally well in shadow or full sunlight, and at any time during the day [10]. Multiplex performs better in estimating wheat LNC, PNC, and NNI than ASD at the early wheat growth stage, while ASD can detect the wheat N indicators more accurately than Multiplex at the late wheat growth stage.

4.2. Accuracy Evaluation of GPR, MLR, and PR Methods

In this study, the GPR model combined with the SBBR algorithm was used to obtain the optimum N-sensitive VI combinations for different sensors. Then, both the MLR and GPR models based on the VI combinations and all VIs (Table 8, Table 9, Table 10 and Table 11), as well as the PR model (Table 7), were investigated for N estimation accuracy through comparing the models’ R2 and RMSE. Figure 3 and Figure 4 show the mean absolute error (MAE) and the Nash–Sutcliffe modeling efficiency (NSE) for LNC, PNC, and NNI through the PR, MLR, and GPR modeling methods.
The ideal MAE value is 0 from the range of 0 to positive infinity, while the optimum NSE value is 1 from the range of negative infinity to 1. As shown in Figure 3, the MAE values for the LNC models built by GPR methods for different sensors are generally less than those of the PR and MLR methods. Meanwhile, the NSE values in Figure 4 indicate that the GPR methods have a higher accuracy, especially when modeling with the ASD and RGB VIs across different wheat growth stages. For PNC estimation, the PR method can obtain a lower MAE at wheat Feekes stage 5, while the GPR models perform better at wheat Feekes stage 11. For NNI estimation, GPR modeling with the selected Multiplex and ASD VIs has a higher MAE than the other methods, indicating that these two sensors may not be suitable for NNI estimation at the early wheat growth stage. It can be seen in Figure 3 and Figure 4 that when the two stages are combined, GPR modeling with all VIs or the selected VIs can estimate N parameters with a lower MAE and a higher NSE than the other modeling methods.
Figure 5 shows the scatterplots of the estimated vs. measured LNC, PNC, and NNI through the GPR-SBBR process based on the best VIs for the four sensors.
Contrary to other methods, the training phase in GPR takes place in a Bayesian framework. Verrelst et al. compared GPR with parametric methods based on established and generic VIs [34]. GPR not only outperformed parametric linear regression methods but also offered additional interesting features that go far beyond what is typically available from parametric or non-parametric approaches [35]. In this study, the GPR model with the selected ASD VIs yielded an obviously higher accuracy for LNC and PNC estimation compared to the PR and MLR models. The result from this study is consistent with previous findings.
Apart from the GPR model, the R2, MAE, and NSE value variation for the MLR and PR models across different wheat growth stages indicates that these models could not produce reliable estimates of winter wheat N status indicators. RGB image-based VIs were used to estimate the N status of winter wheat in field conditions. The observed wheat N parameters were fairly well described by the image data with the GPR model, indicating that a digital camera can be used as a low-cost tool to estimate crop N status and deliver information for N management over a broader region.
The relationships between the VIs from different sensors and N status indicators at the early wheat growth stage were generally better than those at the late stage. The growth stage was found to have a significant effect on the relationships between VIs and N status indicators. This was also confirmed by the results of previous studies [3,38]. However, because the Multiplex or ASD sensor requires measurements made over the plant, covering a small area, it is difficult to apply this type of sensor in monitoring N status in larger areas. Although multispectral or hyperspectral images with near-infrared sensitivity can often convey more information than conventional RGB images, RGB cameras have been more frequently used due to their low cost and ease of use [78]. In this study, RGB image-based VIs were used to estimate the nitrogen status of winter wheat canopies under field conditions. Compared with the other sensors, the UAV-mounted RGB camera is a better selection for crop growth and nutrition monitoring due to its ability to conduct the same type of assessment over hundreds or thousands of plots with a low cost and ease of use, as shown in this study. Different types of sensors can be carried by the UAV to detect the crop growth and nutrition status. The crop information extracted from UAV images is affected not only by the sensor type, but also by the UAV imaging conditions, such as the UAV flight height and speed. In this study, only one sensor (the RGB camera) was mounted on the UAV to test its ability to estimate crop N status. More research is needed to explore the feasibility of this type of camera and other sensors for estimating the spatial distribution of the N status of winter wheat and other crops under different UAV flight heights and speeds over large areas.

5. Conclusions

Crop nitrogen nutrition diagnosis is important for practicing precision farming. Proximal and airborne sensors provide useful information for the diagnosis of crop N nutritional status. In this study, a GPR method with the 10-fold cross-validation SBBR routine was used to identify the best VIs sensitive to wheat LNC, PNC, and NNI for four different sensors. The observed wheat N parameters were well described by the GPR and traditional PR methods throughout the two growth stages. The following conclusions can be drawn from this study:
At the early wheat growth stage, the Dualex NBI is a good leaf-level indicator for wheat LNC, PNC, and NNI estimation. Meanwhile, the combination of Multiplex SFR_G, FLAV, and NBI_R, the combination of SFR_R, FLAV, and NBI_R, and the combination of FLAV and FER_RG are the best canopy-level indicators for LNC, PNC, and NNI estimation at the early wheat growth stage. The results indicate that the optical fluorescence sensor provides more accurate estimates of winter wheat N status at a low-canopy coverage condition.
At the late wheat growth stage, the best ASD VIs, including the combination of SR(700,670), SR(418,405), and SR(740,720), the combination of SR(418,405), NDIopt, RDVI, and REFD, and the combination of NDWI, MCARI, MSAVI, and PVR, provided accurate estimates for wheat LNC, PNC, and NNI. The results indicate that the ASD sensor, which collects more detailed information, is an essential tool for N estimation at the late wheat growth stage, although the reflectance patterns of ASD represent the integrated effects of leaf water content, biochemical constituents, and various components of the plant structure.
This study reveals that the GPR with SBBR method provides more accurate estimates of winter wheat LNC, PNC, and NNI with the best VI combinations of Dualex, Multiplex, ASD, and RGB camera sensors, compared with GPR with full VIs or the MLR and PR methods.

Author Contributions

Conceptualization, X.S.; methodology, X.S. and D.Z.; investigation, G.Y. and X.X.; data curation, X.S., X.X. and H.F.; writing—original draft preparation, X.S. and D.Z.; writing—review and editing, C.Y. and G.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Key-Area Research and Development Program of Guangdong Province (2019B020214002), the National Key Technologies of Research and Development Program of China “Intergovernmental Cooperation in International Science and Technology Innovation” (2019YFE0125300), and the National Natural Science Foundation of China (41701375).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are not publicly available due to the data also forms part of an ongoing study at this time.

Acknowledgments

We appreciate H. Chang and W.G. Li for their help during field data collection.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mistele, B.; Schmidhalter, U. Estimating the nitrogen nutrition index using spectral canopy reflectance measurements. Eur. J. Agron. 2008, 29, 184–190. [Google Scholar] [CrossRef]
  2. Ziadi, N.; Brassard, M.; Bélanger, G.; Claessens, A.; Tremblay, N.; Cambouris, A.N.; Nolin, M.C.; Parent, L.E. Chlorophyll measurements and nitrogen nutrition index for the evaluation of corn nitrogen status. Agron. J. 2008, 100, 1264–1273. [Google Scholar] [CrossRef]
  3. Cao, Q.; Miao, Y.; Wang, H.; Huang, S.; Cheng, S.; Khosla, R.; Jiang, R. Nondestructive estimation of rice plant nitrogen status with Crop Circle multispectral active canopy sensor. Field Crops Res. 2013, 154, 133–144. [Google Scholar] [CrossRef]
  4. Diacono, M.; Rubino, P.; Montemurro, F. Precision nitrogen management of wheat: A review. Agron. Sustain. Dev. 2013, 33, 219–241. [Google Scholar] [CrossRef]
  5. Hansen, P.M.; Schjoerring, J.K. Reflectance measurement of canopy biomass and nitrogen status in wheat crops using normalized difference vegetation indices and partial least squares regression. Remote Sens. Environ. 2003, 86, S0034–S4257. [Google Scholar] [CrossRef]
  6. Ecarnot, M.; Compan, F.; Roumet, P. Assessing leaf nitrogen content and leafmass per unit area of wheat in the field throughout plant cycle with a portable spectrometer. Field Crops Res. 2013, 140, 44–50. [Google Scholar] [CrossRef]
  7. Debaeke, P.; Rouet, P.; Justes, E. Relationship between the normalized SPAD index and the nitrogen nutrition index: Application to durum wheat. J. Plant Nutr. 2006, 29, 75–92. [Google Scholar] [CrossRef]
  8. Prost, L.; Jeuffroy, M.H. Replacing the nitrogen nutrition index by the chlorophyll meter to assess wheat N status. Agron. Sustain. Dev. 2007, 27, 321–330. [Google Scholar] [CrossRef] [Green Version]
  9. Cao, Q.; Cui, Z.; Chen, X.; Khosla, R.; Dao, T.H.; Miao, Y. Quantifying spatial variability of indigenous nitrogen supply for precision nitrogen management in small scale farming. Precis. Agric. 2012, 13, 45–61. [Google Scholar] [CrossRef]
  10. Tremblay, N.; Wang, Z.; Cerovic, Z.G. Sensing crop nitrogen status with fluorescence indicators—A review. Agron. Sust. Dev. 2012, 32, 451–464. [Google Scholar] [CrossRef] [Green Version]
  11. Gabriel, J.L.; Zarco-Tejada, P.J.; López-Herrera, P.J.; Pérez-Martín, E.; Alonso-Ayuso, M.; Quemada, M. Airborne and ground level sensors for monitoring nitrogen status in a maize crop. Biosyst. Eng. 2017, 160, 124–133. [Google Scholar] [CrossRef]
  12. Cerovic, Z.G.; Ounis, A.; Cartelat, A.; Latouche, G.; Goulas, Y.; Meyer, S. The use of chlorophyll fluorescence excitation spectra for the non-destructive in situ assessment of UV absorbing compounds in leaves. Plant Cell Environ. 2002, 25, 1663–1676. [Google Scholar] [CrossRef]
  13. Cartelat, A.; Cerovic, Z.G.; Goulas, Y.; Meyer, S.; Lelarge, C.; Prioul, J.L. Optically assessed contents of leaf polyphenolics and chlorophyll as indicators of nitrogen deficiency in wheat (Triticum aestivum L.). Field Crops Res. 2005, 91, 35–49. [Google Scholar] [CrossRef]
  14. Lamb, D.W.; Steyn-Ross, M.; Schaare, P.; Hanna, M.M. Estimating leaf nitrogen concentration in ryegrass (Lolium spp.) pasture using the chlorophyll rededge: Theoretical modeling and experimental observations. Int. J. Remote Sens. 2002, 23, 3619–3648. [Google Scholar] [CrossRef]
  15. Reyniers, M.; Vrindts, E. Measuring wheat nitrogen status from space and ground-based platform. Int. J. Remote Sens. 2006, 27, 549–567. [Google Scholar] [CrossRef]
  16. Cho, M.A.; Skidmore, A.K. A new technique for extracting the red edge position from hyperspectral data: The linear extrapolation method. Remote Sens. Environ. 2006, 101, 181–193. [Google Scholar] [CrossRef]
  17. Ben, G.N.; Cerovic, Z.G.; Germain, C.; Toutain, S.; Latouche, G. Non-destructive optical monitoring of grape maturation by proximal sensing. Sensors 2010, 10, 10040–10068. [Google Scholar] [CrossRef]
  18. Mukhamediev, R.I.; Symagulov, A.; Kuchin, Y.; Zaitseva, E.; Bekbotayeva, A.; Yakunin, K.; Assanov, I.; Levashenko, V.; Popova, Y.; Akzhalova, A.; et al. Review of Some Applications of Unmanned Aerial Vehicles Technology in the Resource-Rich Country. Appl. Sci. 2021, 11, 10171. [Google Scholar] [CrossRef]
  19. Shakhatreh, H.; Sawalmeh, A.H.; Al-Fuqaha, A.; Dou, Z.; Almaita, E.; Khalil, I.; Othman, N.S.; Khreishah, A.; Guizani, M. Unmanned aerial vehicles (UAVs): A survey on civil applications and key research challenges. IEEE Access 2019, 7, 48572–48634. [Google Scholar] [CrossRef]
  20. Hu, H.; Zhang, J.; Sun, X.; Zhang, X. Estimation of leaf chlorophyll content of rice using image color analysis. Can. J. Remote Sens. 2013, 39, 185–190. [Google Scholar] [CrossRef]
  21. Jia, B.; He, H.B.; Ma, F.Y.; Diao, M.; Jiang, G.Y.; Zheng, Z.; Cui, J.; Fan, H. Use of a digital camera to monitor the growth and nitrogen status of cotton. Sci. World J. 2014, 2014, 602647. [Google Scholar] [CrossRef]
  22. Adamsen, F.G.; Pinter, J.P.; Barnes, M.E.; LaMorte, L.R.; Wall, W.G.; Leavitt, W.S. Measuring wheat senescence with a digital camera. Crop. Sci. 1999, 39, 719–724. [Google Scholar] [CrossRef]
  23. Lukina, E.V.; Stone, M.L.; Raun, W.R. Estimating vegetation coverage in wheat using digital images. J. Plant Nutr. 1999, 22, 341–350. [Google Scholar] [CrossRef]
  24. Elsayed, S.; Barmeier, G.; Schmidhalter, U. Passive Reflectance Sensing and Digital Image Analysis Allows for Assessing the Biomass and Nitrogen Status of Wheat in Early and Late Tillering Stages. Front. Plant. Sci. 2018, 9, 1478. [Google Scholar] [CrossRef]
  25. Song, S.; Chong, H.; Chun, S.; Cheng, Y.; Dai, H. Research advancement on crop nitrogen nutrition diagnosis. Chin. J. Soil Sci. 2006, 37, 369–372, (In Chinese with English Abstract). [Google Scholar]
  26. Guo, B.B.; Qi, S.L.; Heng, Y.R.; Duan, J.; Zhao, H.Y.; Zhang, Y.P.; Wu, W.F.; Xie, Y.X.; Zhu, Y.J. Remotely assessing leaf N uptake in winter wheat based on canopy hyperspectral red-edge absorption. Eur. J. Agron. 2017, 82, 113–124. [Google Scholar] [CrossRef]
  27. Jacquemoud, S.; Baret, F. PROSPECT: A Model of Leaf Optical Properties Spectra. Remote Sens. Environ. 1990, 34, 75–91. [Google Scholar] [CrossRef]
  28. Ollinger, S.V. Sources of variability in canopy reflectance and the convergent properties of plants. New Phytol. 2011, 189, 375–394. [Google Scholar] [CrossRef] [PubMed]
  29. Verrelst, J.; Romijn, E.; Kooistra, L. Mapping vegetation density in aheterogeneous river floodplain ecosystem using pointable CHRIS/PROBA data. Remote Sens. 2012, 4, 2866–2889. [Google Scholar] [CrossRef] [Green Version]
  30. Verrelst, J.; Camps-Valls, G.; Muñoz-Marí, J.; Rivera, J.P.; Veroustraete, F.; Clevers, J.G.P.W.; Moreno, J. Optical remote sensing and the retrieval of terrestrial vegetation bio-geophysical properties—A review. ISPRS J. Photogramm. Remote Sens. 2015, 108, 273–290. [Google Scholar] [CrossRef]
  31. Clevers, J.G.P.W. Beyond NDVI: Extraction of biophysical variables from remote sensing imagery. In Land Use and Land Cover Mapping in Europe: Practices and Trends; Manakos, I., Braun, M., Eds.; Springer: Berlin/Heidelberg, Germany, 2014; pp. 363–381. [Google Scholar]
  32. Glenn, E.P.; Huete, A.R.; Nagler, P.L.; Nelson, S.G. Relationship between remotely-sensed vegetation indices, canopy attributes and plant physiological processes: What vegetation indices can and cannot tell us about the landscape. Sensors 2008, 8, 2136–2160. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Blum, A.; Langley, P. Selection of relevant features and examples in machine learning. Artif. Intell. 1997, 97, 245–271. [Google Scholar] [CrossRef] [Green Version]
  34. Verrelst, J.; Rivera, J.P.; Gitelson, A.; Delegido, J.; Moreno, J.; Camps, V.G. Spectral band selection for vegetation properties retrieval using Gaussian processes regression. Int. J. Appl. Earth Observ. Geoinf. 2016, 52, 554–567. [Google Scholar] [CrossRef]
  35. Miphokasap, P.; Honda, K.; Vaiphasa, C.; Souris, M.; Nagai, M. Estimating canopy nitrogen concentration in sugarcane using field imaging spectroscopy. Remote Sens. 2012, 4, 1651–1670. [Google Scholar] [CrossRef] [Green Version]
  36. Verrelst, J.; Alonso, L.; Rivera Caicedo, J.; Moreno, J.; Camps-Valls, G. Gaussian process retrieval of chlorophyll content from imaging spectroscopy data. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2013, 6, 867–874. [Google Scholar] [CrossRef]
  37. Verrelst, J.; Alonso, L.; Camps-Valls, G.; Delegido, J.; Moreno, J. Retrieval of vegetation biophysical parameters using gaussian process techniques. IEEE Trans. Geosci. Remote Sens. 2012, 50, 1832–1843. [Google Scholar] [CrossRef]
  38. Verrelst, J.; Malenovský, Z.; Van der Tol, C.; Camps-Valls, G.; Gastellu-Etchegorry, J.P.; Lewis, P.; North, P.; Moreno, J. Quantifying Vegetation Biophysical Variables from Imaging Spectroscopy Data: A Review on Retrieval Methods. Surv. Geophys. 2018, 40, 589–629. [Google Scholar] [CrossRef] [Green Version]
  39. Rasmussen, C.E.; Williams, C.K.I. Gaussian Processes for Machine Learning; The MIT Press: Cambridge, MA, USA, 2006; pp. 7–30. Available online: http://www.gaussianprocess.org/gpml/chapters/RW.pdf (accessed on 28 July 2021).
  40. Li, F.; Miao, Y.; Hennig, S.D.; Gnyp, M.L.; Chen, X.; Jia, L. Evaluating hyperspectral vegetation indices for estimating nitrogen concentration of winter wheat at different growth stages. Precis. Agric. 2010, 11, 335–357. [Google Scholar] [CrossRef]
  41. Ranjan, R.; Chopra, U.K.; Sahoo, R.N.; Singh, A.K.; Pradhan, S. Assessment of plant nitrogen stress in wheat (Triticum aestivum L.) through hyperspectral indices. Int. J. Remote Sens. 2012, 33, 6342–6360. [Google Scholar] [CrossRef]
  42. Tian, Y.; Gu, K.; Chu, X.; Yao, X.; Cao, W.; Zhu, Y. Comparison of different hyperspectral vegetation indices for canopy leaf nitrogen concentration estimation in rice. Plant. Soil. 2013, 376, 193–209. [Google Scholar] [CrossRef]
  43. Song, X.; Feng, W.; He, L.; Xu, D.; Zhang, H.; Li, X.; Wang, Z.; Coburn, C.A.; Wang, C.; Guo, T. Examining view angle effects on leaf N estimation in wheat using field reflectance spectroscopy. ISPRS J. Photo. Remote Sens. 2016, 122, 57–67. [Google Scholar] [CrossRef]
  44. Tian, Y.; Yao, X.; Yang, J.; Cao, W.; Hannaway, D.B.; Zhu, Y. Assessing newly developed and published vegetation indices for estimating rice leaf nitrogen concentration with ground-and space-based hyperspectral reflectance. Field Crops Res. 2011, 120, 299–310. [Google Scholar] [CrossRef]
  45. Feng, W.; Yao, X.; Zhu, Y.; Tian, Y.C.; Cao, W.X. Monitoring leaf nitrogen status with hyperspectral reflectance in wheat. Eur. J. Agron. 2008, 28, 394–404. [Google Scholar] [CrossRef]
  46. Daniela, S.; Mirco, B.; Pietroalessandro, B.; Stefano, B. Plant nitrogen concentration in paddy rice from field canopy hyperspectral radiometry. Field Crops Res. 2009, 111, 119–129. [Google Scholar] [CrossRef]
  47. Wang, W.; Yao, X.; Tian, Y.; Liu, X.; Ni, J. Estimating leaf nitrogen concentration with three-band vegetation indices in rice and wheat. Field Crops Res. 2012, 129, 90–98. [Google Scholar] [CrossRef]
  48. Rouse, J.W.; Haas, R.H.; Schell, J.A.; Deering, D.W.; Harlan, J.C. Monitoring the Vernal Advancement and Retrogradation of Natural Vegetation; NASA/GSFC, Type III, Final Report; NTRS: Greenbelt, MD, USA, 1973; pp. 7–15. Available online: https://ntrs.nasa.gov/archive/nasa/casi.ntrs.nasa.gov/19740004927.pdf (accessed on 26 July 2021).
  49. Reujean, J.L.; Breon, F.M. Estimating PAR absorbed by vegetation from bidirectional reflectance measurements. Remote Sens Enviro. 1995, 51, 375–384. [Google Scholar] [CrossRef]
  50. Gitelson, A.; Merzlyak, M.N. Spectral reflectance changes associated with autumn senescence of Aesculus hippocastanum L. and Acer platanoides L. leaves: Spectral features and relation to chlorophyll estimation. J. Plant Phys. 1994, 143, 286. [Google Scholar] [CrossRef]
  51. Penuelas, J.; Filella, I.; Serrano, L.; Save, R. The reflectance at the 950–970 nm region as an indicator of plant water status. Int. J. Remote Sens. 1992, 14, 1887–1905. [Google Scholar] [CrossRef]
  52. Gao, B.C. NDWI—A normalized difference water index for remote sensing of vegetation liquid water from space. Remote Sens. Environ. 1996, 58, 257–266. [Google Scholar] [CrossRef]
  53. Hardisky, M.A.; Klemas, V.; Smart, R.M. The influence of Soil Salinity, Growth Form, and Leaf Moisture on the Spectral Radiance of Spartina alterniflora Canopies. Photo. Eng. Remote Sens. 1983, 48, 77–83. [Google Scholar]
  54. Daughtry, C.S.T.; Walthall, C.L.; Kim, M.S.; de colstoun, E.B.; Mcmurtrey, J.E. Estimating corn leaf chlorophyll content from leaf and canopy reflectance. Remote Sens. Environ. 2000, 74, 229–239. [Google Scholar] [CrossRef]
  55. Rondeaux, G.; Steven, M.; Baret, F. Optimization of soil-adjusted vegetation indices. Remote Sens. Environ. Feb. 1996, 55, 95–107. [Google Scholar] [CrossRef]
  56. Qi, J.; Chehbouni, A.; Huete, A.R.; Kerr, Y.H. A Modified Soil Adjusted Vegetation Index (MSAVI). Remote Sens. Environ. 1994, 48, 119–126. [Google Scholar] [CrossRef]
  57. Haboudane, D.; Miller, J.R.; Pattey, E.; Zarco-Tejada, P.J.; Strachan, I.B. Hyperspectral vegetation indices and novel algorithms for predicting green LAI of crop canopies: Modeling and validation in the context of precision agriculture. Remote Sens. Environ. 2004, 90, 337–352. [Google Scholar] [CrossRef]
  58. Metternicht, G. Vegetation indices derived from high resolution airborne videography for precision crop management. Int. J. Rem. Sens. 2003, 24, 2855–2877. [Google Scholar] [CrossRef]
  59. SpecTerra, S. Presentation and Analysis of Data; SpecTerra Services Pty Ltd.: Leederville, Australia, 1999; Available online: http://www.specterra.com.au/dmsv_data_frame.html (accessed on 10 July 2021).
  60. Demmig, A.B.; Adams, W.W. The role of xanthophylls cycle carotenoids in the protection of photosynthesis. Trends Plant. Sci. 1996, 1, 21–26. [Google Scholar] [CrossRef]
  61. Dawson, T.P.; Curran, P.J.; Plummer, S.E. LIBERTY-Modeling the effects of leaf biochemical concentration on reflectance spectra. Remote Sens. Environ. 1998, 65, 50–60. [Google Scholar] [CrossRef]
  62. Miller, J.; Hare, E.; Wu, J. Quantitative characterization of the vegetation red edge reflectance I. An inverted-Gaussian reflectance model. Int. J. Remote Sens. 1990, 11, 1755–1773. [Google Scholar] [CrossRef]
  63. Yue, J.B.; Feng, H.K.; Jin, X.L.; Yuan, H.H.; Li, Z.H.; Zhou, C.Q.; Yang, G.J.; Tian, Q.J. A Comparison of Crop Parameters Estimation Using Images from UAV-Mounted Snapshot Hyperspectral Sensor and High-Definition Digital Camera. Remote Sens. 2018, 10, 1138. [Google Scholar] [CrossRef] [Green Version]
  64. Woebbecke, D.; Meyer, G.; Von Bargen, K.; Mortensen, D. Color indices for weed identification under various soil, residue, and lighting conditions. Trans. ASAE 1995, 38, 259–269. [Google Scholar] [CrossRef]
  65. Meyer, G.E.; Neto, J.C. Verification of color vegetation indices for automated crop imaging applications. Comput. Electron. Agric. 2008, 63, 282–293. [Google Scholar] [CrossRef]
  66. Peñuelas, J.; Gamon, J.A.; Fredeen, A.L.; Merino, J.; Field, C.B. Reflectance indices associated with physiological changes in nitrogen- and water-limited sunflower leaves. Remote Sens. Environ. 1994, 48, 135–146. [Google Scholar] [CrossRef]
  67. Bradstreet, R.B. The Kjeldahl Method for Organic Nitrogen; Academic Press Incorporated: New York, NY, USA, 1965. [Google Scholar]
  68. Zhao, H.T.; Song, X.Y.; Yang, G.J.; Li, Z.H.; Zhang, D.Y.; Feng, H.K. Monitoring of Nitrogen and Grain Protein Content in Winter Wheat Based on Sentinel-2A Data. Remote Sens. 2019, 11, 1724. [Google Scholar] [CrossRef] [Green Version]
  69. Lemaire, G.; Jeuffroy, M.H.; Gastal, F. Diagnosis tool for plant and crop N status in vegetative stage: Theory and practices for crop N management. Eur. J. Agron. 2008, 28, 614–624. [Google Scholar] [CrossRef]
  70. Yue, S.; Meng, Q.; Zhao, R.; Li, F.; Chen, X.; Zhang, F.; Cui, Z. Critical nitrogen dilution curve for optimizing nitrogen management of winter wheat production in the North China Plain. Agron. J. 2012, 104, 523–529. [Google Scholar] [CrossRef]
  71. Willmott, C.J.; Matsuura, K. Advantages of the mean absolute error (MAE) over the root mean square error (RMSE) in assessing average model performance. Clim. Res. 2005, 30, 79–82. [Google Scholar] [CrossRef]
  72. Schaeffer, D.L. A model evaluation methodology applicable to environmental assessment models. Ecol. Model. 1980, 8, 275–295. [Google Scholar] [CrossRef] [Green Version]
  73. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual models part I-A discussion of principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  74. Saberioon, M.; Amin, M.S.M.; Gholizadeh, A.; Ezri, M.H. A review of optical methods for assessing nitrogen contents during rice growth. Appl. Eng. Agric. 2014, 30, 657–669. [Google Scholar] [CrossRef]
  75. Karim, S.T.; Cheng, T.; Liu, X.; Tian, Y.; Zhu, Y.; Cao, W.; Cao, Q. Potential of UAV-based active sensing for monitoring rice leaf nitrogen status. Front. Plant Sci. 2018, 9, 1834. [Google Scholar] [CrossRef] [Green Version]
  76. Yang, G.; Liu, J.; Zhao, C.; Li, Z.; Huang, Y.; Yu, H. Unmanned aerial vehicle remote sensing for field-based crop phenotyping: Current status and perspectives. Front. Plant Sci. 2017, 8, 1111. [Google Scholar] [CrossRef] [PubMed]
  77. Roberts, D.A.; Ustin, S.L.; Ogunjemiyo, S.; Greenberg, J.; Dobrowski, S.Z.; Chen, J.Q.; Hinckley, T.M. Spectral and structural measures of northwest forest vegetation at leaf to landscape scales. Ecosystems 2004, 7, 545–562. [Google Scholar] [CrossRef]
  78. Jayme, G.A.B. Detection of nutrition deficiencies in plants using proximal images and machine learning: A review. Comput. Electron. Agric. 2019, 162, 482–492. [Google Scholar] [CrossRef]
Figure 1. Study area and experimental design. NM treatment: A recommended conventional, uniform nitrogen fertilization treatment; CK treatment: non-nitrogen fertilizer treatment; BH treatment: excess nitrogen fertilizer treatment (N-rich strip treatment); S treatment: variable-rate nitrogen treatment with fertilizer rate based on the SPAD value ratio between this treatment and N-rich strip treatment; A treatment: variable-rate nitrogen treatment with fertilizer rate based on the ASD spectroradiometer vegetation index OSAVI value ratio between this treatment and N-rich strip treatment; T treatment: variable-rate nitrogen treatment with fertilizer rate based on the ASD spectroradiometer vegetation index REPLI value ratio between this treatment and N-rich strip treatment; D treatment: variable-rate nitrogen treatment with fertilizer rate based on the Dualex NBI value between this treatment and N-rich strip treatment; and M treatment: variable-rate nitrogen treatment with fertilizer rate based on the Multiplex NBI_R value ratio between this treatment and N-rich strip treatment.
Figure 1. Study area and experimental design. NM treatment: A recommended conventional, uniform nitrogen fertilization treatment; CK treatment: non-nitrogen fertilizer treatment; BH treatment: excess nitrogen fertilizer treatment (N-rich strip treatment); S treatment: variable-rate nitrogen treatment with fertilizer rate based on the SPAD value ratio between this treatment and N-rich strip treatment; A treatment: variable-rate nitrogen treatment with fertilizer rate based on the ASD spectroradiometer vegetation index OSAVI value ratio between this treatment and N-rich strip treatment; T treatment: variable-rate nitrogen treatment with fertilizer rate based on the ASD spectroradiometer vegetation index REPLI value ratio between this treatment and N-rich strip treatment; D treatment: variable-rate nitrogen treatment with fertilizer rate based on the Dualex NBI value between this treatment and N-rich strip treatment; and M treatment: variable-rate nitrogen treatment with fertilizer rate based on the Multiplex NBI_R value ratio between this treatment and N-rich strip treatment.
Sensors 22 00549 g001
Figure 2. Cross-validation R2 for LNC (a), PNC (b), and NNI (c) estimation with different sensors through PR, MLR, and GPR methods.
Figure 2. Cross-validation R2 for LNC (a), PNC (b), and NNI (c) estimation with different sensors through PR, MLR, and GPR methods.
Sensors 22 00549 g002
Figure 3. The MAE for LNC (left), PNC (middle), and NNI (right) estimation with different sensors through PR, MLR, and GPR methods.
Figure 3. The MAE for LNC (left), PNC (middle), and NNI (right) estimation with different sensors through PR, MLR, and GPR methods.
Sensors 22 00549 g003
Figure 4. The NSE for LNC (left), PNC (middle), and NNI (right) estimation with different sensors through PR, MLR, and GPR methods.
Figure 4. The NSE for LNC (left), PNC (middle), and NNI (right) estimation with different sensors through PR, MLR, and GPR methods.
Sensors 22 00549 g004
Figure 5. Measured vs. estimated LNC (left), PNC (middle), and NNI (right) values along the 1:1 line through the GPR-SBBR method based on best VIs for four different sensors. (a) LNC based on RGB sensor; (b) PNC based on RGB sensor; (c) NNI based on RGB sensor; (d) LNC based on ASD sensor; (e) PNC based on ASD sensor; (f) NNI based on ASD sensor; (g) LNC based on Multiplex sensor; (h) PNC based on Multiplex sensor; (i) NNI based on Multiplex sensor; (j) LNC based on Dualex sensor; (k) PNC based on Dualex sensor; and (l) NNI based on Dualex sensor.
Figure 5. Measured vs. estimated LNC (left), PNC (middle), and NNI (right) values along the 1:1 line through the GPR-SBBR method based on best VIs for four different sensors. (a) LNC based on RGB sensor; (b) PNC based on RGB sensor; (c) NNI based on RGB sensor; (d) LNC based on ASD sensor; (e) PNC based on ASD sensor; (f) NNI based on ASD sensor; (g) LNC based on Multiplex sensor; (h) PNC based on Multiplex sensor; (i) NNI based on Multiplex sensor; (j) LNC based on Dualex sensor; (k) PNC based on Dualex sensor; and (l) NNI based on Dualex sensor.
Sensors 22 00549 g005
Table 1. Fertilizer information for different treatments.
Table 1. Fertilizer information for different treatments.
TreatmentPlotBase FertilizerTopdressing FertilizerFertilizer Treatment Rate Statistic
NumberTimeN kg/haTimeN kg/haMean kg/haMin kg/haMax kg/haCV %
BH18Seed72Feekes 2, 451,1022252252250
NM13Seed72Feekes 6781501501500
CK13Seed0Feekes 600000
A12Seed72Feekes 678150147154.11.78
M12Seed72Feekes 678150138.6162.25.15
D12Seed72Feekes 678150141.61603.04
S12Seed72Feekes 678150144.7154.81.92
T12Seed72Feekes 678150131183.110.25
Table 2. Parameters on sensors used in this study.
Table 2. Parameters on sensors used in this study.
Sensor InformationPolyphenol and Chlorophyll MeterPolyphenol and Chlorophyll MeterField SpectrometerUAV-Based Digital Camera
Sensor TypeDualexMultiplexASDRGB Camera
Sensor nameForce-A Dualex ScientificForce-A MULTIPLEX 3ASD FieldSpec 4Sony DSC–QX100
Target samplePlant leavesPlant canopyPlant canopyPlant canopy
Field of view--25°64°
Image size---3000 × 4000
Working height-10 cm1.3 m50 m
Measurement area5 mm in diameter10 cm in diameter50 cm in diameterFull field
Spectral informationExcitation channels: UV (357 nm) and red (650 nm). Detection channels: red and far-red. Excitation channels: UV (375 nm), blue (450 nm), green (510 nm), and red (630 nm). Detection channels: yellow, red, and far-red. 350–2500 nmR,G,B
Original spectral resolution--3 nm @ 700 nm;
10 nm @ 1400 nm
-
Data spectral resolution 1 nm
Image spatial resolution---2 cm
Table 3. Vegetation indices based on Multiplex 3 and ASD sensors used in this study.
Table 3. Vegetation indices based on Multiplex 3 and ASD sensors used in this study.
SensorIDVegetation IndexFormulaReference
Multiplex1SFR_GFRF_G/RF_G[17]
2SFR_RFRF_R/RF_R[17]
3BRR_FRFBGF_UV/FRF_UV[17]
4FER_RUVFRF_R/FRF_UV[17]
5FER_RGFRF_R/FRF_G[17]
6FLAVLog(FER_RUV)[17]
7ANTHLog(FER_RG)[17]
8NBI_GFRF_UV/RF_G[17]
9NBI_RFRF_UV/RF_R[17]
ASD1SR(700,670)(R700)/(R670)[40]
2SR(418,450)(R418)/(R450)[40]
3VOGa(R740)/(R720)[41]
4SR(553,537)(R553)/(R537)[42]
5NDCI(R762 − R527)/(R762 + R527)[41]
6NDRE(R790 − R720)/(R790 + R720)[43]
7TBI1 (R434)/(R496 + R401)[44]
8mND705(R750 − R705)/(R750 + R705 − 2R445)[45]
9NDIopt(R503 − R483)/(R503 + R483)[46]
10TBI2 (R924 − R703)/(R924 − R703)[47]
11NDVI(670,780(R780 − R670)/(R780 + R670)[48]
12RDVI(R800 − R670)/(R800 + R670)1/2[49]
13SR(750,700)(R750)/(R700)[50]
14WI(R900)/(R950)[51]
15NDWI(R860 − R1240)/(R860 + R1240)[52]
16NDII(R819 − R1600)/(R819 + R1600)[53]
17MCARI[(R700 − R670) − 0.2(R700 − R550)](R700/R670)[54]
18TCARI3[(R700 − R670) − 0.2(R700 − R550)(R700/R670)][54]
19OSAVI1.16(R800 − R670)/(R800 + R670 + 0.16)[55]
20MSAVI0.5[2R800 + 1 − ((2R800 + 1)2 − 8(R800 − R670))1/2][56]
21MCARI 11.2[2.5(R800 − R670) − 1.3(R800 − R550)][57]
22MCARI 21.5[2.5(R800 − R670) − 1.3(R800 − R550)]/[(2R800 + 1)2 − (6R800 − 5(R670)1/2) − 0.5]1/2[57]
23PPR(R550 − R450)/(R550 + R450)[58]
24PVR(R550 − R650)/(R550 + R650)[59]
25PRI(R531 − R570)/(R531 + R570)[60]
26REP700 + 40[(R670 + R780)/2) − R700)/(R740 − R700)][61]
27REVReflectance value at REP[62]
28REFDFirst deviation of red edge[62]
Table 4. Vegetation indices based on digital RGB images.
Table 4. Vegetation indices based on digital RGB images.
IDVegetation IndexFull NameFormulaReference
1RDN values for red bandDNR/255[64]
2GDN values for green bandDNG/255[64]
3BDN values for blue bandDNB/255[64]
4rChromatic coordinates for redR/(R + G + B)[64]
5gChromatic coordinates for greenG/(R + G + B)[64]
6bChromatic coordinates for blueB/(R + G + B)[64]
7ExRExcess red1.4 × r − b[65]
8ExGExcess green2g − (r + b)[65]
9NDIThe normalized difference vegetation index(b − g)/(b + g)[65]
10CVI1Color vegetation index 1(r − g)[65]
11CVI2Color vegetation index 2(g − b)[65]
12CVI3Color vegetation index 3(g − b)/(r − g)[65]
13GRVIGreen–red vegetation index(g − r)/(g + r)[65]
14NPCINormalized pigment chlorophyll ratio index(b − r)/(b + r)[66]
Table 5. Statistics computed to compare the results of the different VIs used for the N parameter estimation.
Table 5. Statistics computed to compare the results of the different VIs used for the N parameter estimation.
StatisticsFormulaCharacterReference
Root mean square error (RMSE) RMSE = i n ( y i ^ y i ) 2 n from 0 to +∞, optimum 0[71]
Mean absolute error (MAE) MAE = i = 1 n | y i ^ y i | / n from 0 to +∞, optimum 0[72]
Nash–Sutcliffe efficiency (NSE) NSE = 1 i = 1 n ( y i y i ^ ) 2 / i = 1 n ( y i y ¯ ) 2 from −∞ to 1, optimum 1[73]
Notes:   y i ^ is the predicted value of the ith observation, y i   is the measured value of the ith observation, y ¯ is the average of the measured values, and n is the number of observations in the calibration set.
Table 6. Descriptive statistics of winter wheat LNC, PNC, and NNI across different growth stages.
Table 6. Descriptive statistics of winter wheat LNC, PNC, and NNI across different growth stages.
Growth StageParameterMinMaxMeanRangeStdCV (%)
Feekes 5LNC (%)3.075.164.332.090.4911.30
PNC (%)2.604.673.772.070.5113.42
NNI0.741.581.140.840.1917.17
Feekes 11LNC (%)1.082.992.321.920.3113.42
PNC (%)0.771.581.330.810.128.99
NNI0.441.170.890.730.1111.97
Feekes 5–11LNC (%)1.085.163.334.081.0832.34
PNC (%)0.774.672.553.901.2649.46
NNI0.441.581.001.140.2018.59
Table 7. Relationship between wheat N variables and VIs from different sensors at different growth stages (n = 104).
Table 7. Relationship between wheat N variables and VIs from different sensors at different growth stages (n = 104).
SensorFeekes StageLNC (%)PNC (%)NNI
ModelVIR2RMSEModelVIR2RMSEModelVIR2RMSE
RGB5PolyNPCI0.43 **0.23PolyExR0.33 **0.10PolyExR0.31 **0.09
11PolyNDI0.39 **0.38PolyNDI0.39 **0.26ExpNDI0.36 **0.16
5–11Polyr0.61 **0.68Pow ExR0.62 **0.40LogExR0.37 **0.16
ASD5PolyNDIopt0.190.44PolyNDIopt0.190.30PolySR(553,537)0.30 **0.16
11PolymND7050.40 **0.24PolyNDRE0.46 **0.09PolyNDRE0.39 **0.08
5–11PolyNDIopt0.85 **0.43PolyNDIopt0.84 **0.26LogNDIopt0.53 **0.14
Multiplex5PolyFLAV0.55 **0.33PolyFLAV0.52 **0.24PolyNBI_R0.42 **0.15
11PolySFR_R0.34 **0.25PolySFR_R0.230.11PolySFR_R0.21*0.10
5–11PolyBRR_FRF0.87 **0.39PolyBRR_FRF0.86 **0.24PolyNBI_G0.56 **0.13
Dualex5LogNBI0.80 **0.22Pow NBI0.79 **0.23LinNBI0.49 **0.14
11PolyCHI0.36 **0.27LogCHI0.160.11LogCHI0.20 *0.10
5–11PolyNBI0.57 **0.72PolyNBI0.50 **0.92ExpNBI0.46 **0.15
Note: Lin: linear, Exp: exponential, Pow: power, Poly: second-order polynomial, Log: logarithmic; ** r (0.01, 70) = 0.302, indicates significance at the 0.01 probability level; * r (0.05, 70) = 0.232, indicates significance at the 0.05 probability level.
Table 8. Cross-validation results for N estimation through GPR and MLR based on UAV-mounted RGB camera VIs.
Table 8. Cross-validation results for N estimation through GPR and MLR based on UAV-mounted RGB camera VIs.
NFeekes VIGPR-SBBRMLRVIsGPRMLR
VariableStageR2RMSER2RMSER2RMSER2RMSE
LNC5B, b0.440.370.430.37All 14 VIs0.410.370.280.43
11NPCI0.430.240.360.25All 14 VIs0.400.250.330.26
5–11B, g, CVI20.820.460.620.67All 14 VIs0.810.470.680.62
PNC5B, b0.420.380.410.38All 14 VIs0.390.390.290.44
11G, NDI0.410.090.290.10All 14 VIs0.380.090.400.09
5–11B, g, CVI20.890.430.590.82All 14 VIs0.870.450.670.74
NNI5B, b0.350.160.350.16All 14 VIs0.280.160.330.16
11R, NDI, CVI20.330.090.340.09All 14 VIs0.330.090.260.10
5–11B, G, g, NDI, CVI2, CVI30.540.140.380.16All 14 VIs0.540.130.450.15
Table 9. Cross-validation results for N estimation through GPR and MLR methods based on ASD VIs.
Table 9. Cross-validation results for N estimation through GPR and MLR methods based on ASD VIs.
NFeekesVIGPR-SBBRMLRVIsGPRMLR
VariableStageR2RMSER2RMSER2RMSER2RMSE
LNC5TBI1, NDIopt, TBI20.470.360.310.41All 28 VIs0.280.420.340.42
11SR(700,670), SR(418,405), SR(740,720)0.510.220.330.25All 28 VIs0.390.260.420.26
5–11SR(418,405), TBI1, PPR,
PRI, REFD
0.930.290.880.37All 28 VIs0.920.300.910.34
PNC5TBI1, NDIopt, TBI20.460.370.320.42All 28 VIs0.270.440.330.44
11SR(418,405), NDIopt,
RDVI, REFD
0.490.090.380.09All 28 VIs0.350.100.380.10
5–11TBI1, PPR, REV, REFD0.960.270.890.42All 28 VIs0.950.280.940.31
NNI5MSAVI, PPR0.510.140.510.14All 28 VIs0.410.150.350.17
11NDWI, MCARI,
MSAVI, PVR
0.390.080.210.10All 28 VIs0.280.100.210.10
5–11SR(418,405), NDIopt, MSAVI, MCARI1, PPR0.660.120.560.13All 28 VIs0.590.120.590.13
Table 10. Cross-validation results for N estimation through GPR and MLR methods based on Multiplex VIs.
Table 10. Cross-validation results for N estimation through GPR and MLR methods based on Multiplex VIs.
NFeekesVIGPR-SBBRMLRVIsGPRMLR
VariableStageR2RMSER2RMSER2RMSER2RMSE
LNC5SFR_G, FLAV, NBI_R0.590.310.590.31All 9 VIs 0.570.310.580.32
11SFR_R0.440.280.270.27All 9 VIs0.390.300.160.29
5–11SFR_R, BRR_FRF, NBI_R0.930.290.880.38All 9 VIs0.930.300.900.34
PNC5SFR_R, FLAV, NBI_R0.580.340.570.34All 9 VIs0.530.340.530.35
11SFR_R0.340.100.150.11All 9 VIs0.250.110.240.11
5–11SFR_R, BRR_FRF, NBI_G0.960.250.900.41All 9 VIs0.960.260.930.35
NNI5FLAV, FER_RG0.520.340.520.14All 9 VIs0.520.140.440.15
11SFR_G, SFR_R0.310.090.240.09All 9 VIs0.180.100.250.10
5–11SFR_G, SFR_R, NBI_G0.590.130.580.13All 9 VIs0.580.130.580.13
Table 11. Cross-validation results through GPR and MLR methods at different wheat growth stages.
Table 11. Cross-validation results through GPR and MLR methods at different wheat growth stages.
NFeekesVIGPR-SBBRMLRVIsGPRMLR
VariableStageR2RMSER2RMSER2RMSER2RMSE
LNC5NBI0.800.220.790.23All 3 VIs0.770.220.800.22
11Chl0.550.210.120.31All 3 VIs0.520.210.100.36
5–11NBI, Chl0.830.430.600.69All 3 VIs0.830.430.740.56
PNC5NBI0.790.230.790.23All 3 VIs0.780.220.790.23
11Chl0.200.110.110.11All 3 VIs0.120.120.160.11
5–11NBI, Chl0.820.550.550.86All 3 VIs0.810.540.670.74
NNI5NBI0.510.140.490.14All 3 VIs0.510.140.470.14
11Chl0.270.100.100.10All 3 VIs0.220.100.260.09
5–11NBI, Chl0.600.130.420.15All 3 VIs0.590.130.480.15
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Song, X.; Yang, G.; Xu, X.; Zhang, D.; Yang, C.; Feng, H. Winter Wheat Nitrogen Estimation Based on Ground-Level and UAV-Mounted Sensors. Sensors 2022, 22, 549. https://doi.org/10.3390/s22020549

AMA Style

Song X, Yang G, Xu X, Zhang D, Yang C, Feng H. Winter Wheat Nitrogen Estimation Based on Ground-Level and UAV-Mounted Sensors. Sensors. 2022; 22(2):549. https://doi.org/10.3390/s22020549

Chicago/Turabian Style

Song, Xiaoyu, Guijun Yang, Xingang Xu, Dongyan Zhang, Chenghai Yang, and Haikuan Feng. 2022. "Winter Wheat Nitrogen Estimation Based on Ground-Level and UAV-Mounted Sensors" Sensors 22, no. 2: 549. https://doi.org/10.3390/s22020549

APA Style

Song, X., Yang, G., Xu, X., Zhang, D., Yang, C., & Feng, H. (2022). Winter Wheat Nitrogen Estimation Based on Ground-Level and UAV-Mounted Sensors. Sensors, 22(2), 549. https://doi.org/10.3390/s22020549

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