Remotely Estimating Aerial N Uptake in Winter Wheat Using Red-Edge Area Index From Multi-Angular Hyperspectral Data

Remote sensing techniques can be efficient for non-destructive, rapid detection of wheat nitrogen (N) nutrient status. In the paper, we examined the relationships of canopy multi-angular data with aerial N uptake of winter wheat (Triticum aestivum L.) across different growing seasons, locations, years, wheat varieties, and N application rates. Seventeen vegetation indices (VIs) selected from the literature were measured for the stability in estimating aerial N uptake of wheat under 13 view zenith angles (VZAs) in the solar principal plane (SPP). In total, the back-scatter angles showed better VI behavior than the forward-scatter angles. The correlation coefficient of VIs with aerial N uptake increased with decreasing VZAs. The best linear relationship was integrated with the optimized common indices DIDA and DDn to examine dynamic changes in aerial N uptake; this led to coefficients of determination (R²) of 0.769 and 0.760 at the −10° viewing angle. Our novel area index, designed the modified right-side peak area index (mRPA), was developed in accordance with exploration of the spectral area calculation and red-edge feature using the equation: mRPA = (R₇₆₀/R₆₀₀)^1/2 × (R₇₆₀-R₇₁₈). Investigating the predictive accuracy of mRPA for aerial N uptake across VZAs demonstrated that the best performance was at −10° [R² = 0.804, p < 0.001, root mean square error (RMSE) = 3.615] and that the effect was relatively similar between −20° to +10° (R² = 0.782, p < 0.001, RMSE = 3.805). This leads us to construct a simple model under wide-angle combinations so as to improve the field operation simplicity and applicability. Fitting independent datasets to the models resulted in relative error (RE, %) values of 12.6, 14.1, and 14.9% between estimated and measured aerial N uptake for mRPA, DIDA, and DDn across the range of −20° to +10°, respectively, further confirming the superior test performance of the mRPA index. These results illustrate that the novel index mRPA represents a more accurate assessment of plant N status, which is beneficial for guiding N management in winter wheat.

Introduction

Nitrogen (N) is a vital element for higher photosynthetic functioning; N resource management is a major factor that can enhance plant growth and influence the quality of plant crops (Woodard and Bly, 1998; Smil, 2002). To ensure productivity, crop producers commonly supply plants in the field with N fertilizers. N supply generally appreciably surpasses plant N uptake, leading to the loss of nitrate through soil leaching, increased greenhouse gas (N₂O) emission, and ground water pollution (Sehy et al., 2003; Ju et al., 2006). To minimize potential N losses, N fertilizer should be applied at the correct time and according to the requirements of the crops. The development, therefore, of techniques which insure higher yield and good product while reducing ecological environment pollution attributed to unsuitable N application is essential.

In currently, remote sensing technique is among the most promising approach that has been shown to rapidly predict the spatial-temporal variability of crops and monitor crop growth status (Hansen and Schjoerring, 2003; Ecarnot et al., 2013). Many indices have been constructed by extracting characteristic spectral information for evaluating biochemical properties in crop plants (Hatfield et al., 2008; Ecarnot et al., 2013). Several researchers have demonstrated a close relationship between NDVI- and RVI-like spectral indices and aboveground N uptake (Mistele and Schmidhalter, 2008; Li et al., 2013), chlorophylls and carotenoids (Blackburn, 1998), canopy leaf biomass (Le Maire et al., 2008). An additional group of vegetation indices are constructed by the forms of three band combination. For instance, Wang et al. (2012) added the 2 × R₄₂₃ band to the NDVI (R₇₀₃, R₉₂₄) and effectively improved the sensitivity of leaf nitrogen concentration (LNC) estimation in rice and wheat. Feng et al. (2015) showed that the three-band spectral index (R₇₅₉-1.8 × R₄₁₉)/(R₇₄₂-1.8 × R₄₁₉) was a good indicator of above ground N uptake in wheat. A third group vegetation indices were developed by area-based algorithm. These include the triangle vegetation index (TVI) and modified TVI for green leaf area index (LAI) (Broge and Leblanc, 2000; Haboudane et al., 2004); the red-edge reflectance curve area for green biomass (Ren et al., 2011); and the adjusted TVI for aerial N uptake (Li et al., 2013). Additionally, remote sensing technology was also applied in the field of phenotype. Rothamsted Research reported that the data collecting from Unmanned Aerial Vehicle (UAV) based remote sensing could rapidly and accurately measure the wheat plant height and growth rate (Holman et al., 2016). Andradesanchez et al. (2014) investigated that the tractor-based phenotyping system could acquire and record data for canopy temperature, height and reflectance of cotton plants at much higher rates. However, the canopy spectral reflectance was sampled only from the vertical observation angle in prior researches, and the nadir observation were difficult to extract spatial structure from the middle and lower layers of plants (Thenkabail et al., 2000; Erdle et al., 2011).

Compared with the nadir observations, multi-angle observations contain more detailed and reliable canopy structure information that permits effective monitoring of crop N status in the middle and lower layers and provides a novel approach for quantitative remote sensing (Pocewicz et al., 2007; Huang et al., 2011). To date, many studies have shown that multi-angle measurements could improve the performance of indices when estimating the structural characteristics of ground objects (Shibayama and Wiegand, 1985; Diner et al., 1999). For instance, Galvão et al. (2009) showed that the varieties of soybean could be distinguished with the best predictive ability in the backward scattering direction. Gemmell and McDonald (2000) showed that the performance of indices under off-nadir angle can effectively discriminate forest cover and LAI. Furthermore, some studies have used to multi-angular datasets to assess plant variables, particularly biochemical components (Hasegawa et al., 2010; Huang et al., 2011). The effect of indices in estimating agronomic parameters and yields changes with the LAI and VZAs (Gemmell and McDonald, 2000; Inoue et al., 2008). Stagakis et al. (2010) focused on using satellite spectral data to estimate chlorphyll a (chl a), chl b, and carotenoids of semi-deciduous shrubs by utilizing different viewing angles and narrow-band indices. He et al. (2016) developed a multi-angular VI to enhance the estimation stability and accuracy of leaf nitrogen concentrations. No matter what the understory vegetation was (green or senesced) the relationships between canopy or total LAI and NDVI or Enhanced Vegetation Index (EVI) varied little across VZAs in pine forests (Pocewicz et al., 2007). However, only a few researches systematically tested the ability of multi-angle spectral data for predicting aerialN uptake of wheat. Taken together, these studies report the construction of a range of novel indices that use multi-angle datasets to enhance the precision and robustness of prediction indices for plant biophysical traits.

The specific aims of the paper were: to (1) study the performance of ground-based spectra and common indices to detect the aerial N uptake of winter wheat under different VZAs; (2) construct an improved, novel VI for aerial N uptake prediction; (3) compare the aerial N uptake predictive ability of the novel model with published VIs; and (4) establish the best observation angle and the best estimation model for aerial N uptake. The results of this study provide technical knowledge and a theoretical basis for monitoring N status by remote sensing technology; the information obtained from these techniques can then be used to help guide appropriate N fertilization application of wheat.

Materials and Methods

Experimental Fields

The field experiments were designed over a 4-year period in Zhengzhou and Shangshui city, China. The different locations, N fertilizer rates, wheat cultivars, and growth seasons were used (Table 1). Urea as N sources was divided into two equal doses, one administered before seeding and the rest at jointing period. Before seeding, 150 kg ha⁻¹ P₂O₅ [as Ca(H₂PO₄)₂] and 90 kg ha⁻¹ K₂O (as KCl) were used to all treatments. The N treatments with triplicates were assigned as completely random blocks in the experiment. The density of the seedlings was 3.0 × 10⁶ plants ha⁻¹.

TABLE 1

Table 1. The experimental conditions, N fertilizer levels, and measured stages.

Data Acquisition

Measurement of Canopy Multi-Angular Hyperspectral Reflectance

Canopy reflectance spectra were obtained in a 1 m² area in each plot under sunny and windless conditions between11:00 to 13:00 using an ASD (Analytical Spectral Devices Inc., Boulder, CO, USA) FieldSpec Handheld spectrometer. This spectrum instrument was equipped with 25°-field-of-view optics fiber, sampling interval of 1.6 nm and spectral resolution of 3.5 nm from 325 to 1,075 nm. The multi-angle data are obtained with a Field Goniometer System, which was designed based on the system developed by Sandmeier and Itten (1999). The goniometer is a device used to position a sensor at these different angles and azimuths (Figure 1). The observed azimuth was fixed relative to the direction of the sun, and the measured plane was defined as the SPP (Myneni et al., 1995). The VZA was divided into backward direction (the observation direction same to the sun, −) and the forward direction (the observation direction against to the sun, +); the nadir position was defined as 0°. The VZA from backward to forward direction is −60, −50, −40, −30, −20, −10, 0, 10, 20, 30, 40, 50, 60°. The VZA become larger from 0° to ±60°, regardless of the backward and forward direction. The 10 sites were averaged to a single spectral sample of each plot. The black and base-line reflectance was calculated by a 40 × 40 cm BaSO₄.

FIGURE 1

Figure 1. Dimensions and design of Field Goniometer System.

Plant Measurements

The areas of 0.20 m² of plant samplings were randomly uprooted from each test area almost simultaneously with the canopy spectral acquisition. The plant samples were weighed after desiccation in an oven at 70°C to a constant weight. The ovendried plants were grinded into powder (1 mm) for N content analysis in laboratory. The aerial nitrogen concentration was determined in line with the micro-Kjeldahl method (Isaac and Johnson, 1976).

Construction of the New VI

The red edge is a region of steep variations in spectral reflectance, and this value may provide some useful information on crop growth and N status (Figure 2) (Sims and Gamon, 2002; Cho and Skidmore, 2006; Hatfield et al., 2008; Feng et al., 2015). To date, few area-based optimized indices have been constructed using red-edge information for non-destructive, rapid assessment of plant N status. Our preliminary research found that the double-peak area parameters constructed based on analysis of the red-edge double-peak characteristics could be effective for assessing Leaf N concentration (LNC) (Feng et al., 2014). There are several techniques to divide the red edge double-peak area into the right-side peak area (RSDR) and the left-side peak area (LSDR) (Figure 3). In this study, the datasets obtained using a 0° observation angle in Exp. 1–5 were used to analyze the relationship between LSDR, RSDR(from different splitting methods), and aerial N uptake. The results showed that RSDR (R₇₆₀-R₇₁₈) divided by characteristic wavelength method had the best performance (R² = 0.740, p < 0.001; Figure 4), which suggested RSDR as potential indicator for estimating aerial N uptake. The previous research showed that the ratios of two or more bands (such as R₈₀₁/R₆₇₀, R₈₀₁/R₅₅₀) could increase sensitivity to crop physiological traits and reduce variation because of external influencial factors (Daughtry et al., 2000). Haboudane et al. (2008) inserted the (R₇₀₀/R₆₇₀)^1/2 into the TVI formula to decrease the combined impacts of the soil background reflectance. We, therefore, inserted a coefficient [in the form of (λ₁/λ₂)^1/2] into RSDR (R₇₆₀-R₇₁₈) to construct a novel VI called the mRPA. The spectral region of the above λ₁ and λ₂ were located within 400 and 900 nm. Figure 5 gave a comprehensive overview of the correlation coefficients for any two band combinations (λ₁, λ₂), and this is valid for selecting the sensitive bands of aerial N uptake. The wavebands λ₁ and λ₂ ranged between 750–900 and 550–650 nm, respectively. This area had the highest precision, with R² values above 0.78, especially R² performed best when λ₁ = 760 nm, λ₂ = 600 nm. Therefore, the final formula of mRPA was:

$\begin{array}{rc}\text{mRPA}={({\text{R}}_{760}/R_{600})}^{1/2}\times (R_{760}-R_{718})& (1)\end{array}$

FIGURE 2

Figure 2. Schematic representation of red edge reflectance curve.

FIGURE 3

Figure 3. The division of the whole red-edge double-peak area into left and right single peak areas.

FIGURE 4

Figure 4. The relationships of aerial N uptake to RSDR (A) and LSDR (B) at 0°view zenith angle.

FIGURE 5

Figure 5. Contour maps of coefficients of determination between aerial N uptake and the mRPA with formula (R_λ1/R_λ2)^1/2 × (R₇₆₀-R₇₁₈) (n = 155, p_0.001 = 0.262).

Model Calibration and Validation

The datasets from Exp. 1–5 were used to construct the evaluating models, and the dataset from Exp. 6 were used to validate these aerial N uptake evaluating models. The correlation of aerial N uptake and VIs was examined by MATLAB 7.0. In this study, 18 VIs were selected and summarized in Table 2. The quantitative relationship between the optimal VI and aerial N uptake could be established based on the highest R². The model's behavior was evaluated by employing the R², RMSE, and RE (%). Among the indices, those with the highest R² and the lowest RMSE and RE were considered as the best. RMSE and RE were calculated from actual and predicted values of samples according to Equations (2) and (3), respectively:

$\begin{array}{rc}RMSE=\sqrt{\frac{1}{n}\times {\displaystyle \sum _{i=1}^n}{(Pi-Oi)}^2}& (2)\end{array}$

$\begin{array}{rc}RE(\%)=\sqrt{\frac{1}{n}\times {\displaystyle \sum _{i=1}^n}{(\frac{Pi-Oi}{Oi})}^2}\times 100\%& (3)\end{array}$

Here, P_i, and O_i represented the estimated and measured values, respectively, and n represented the sampling number. The prediction was considered to be excellent if RE < 10%, good if RE = 10–20%, fair if RE = 20–30%, and poor if RE > 30% (Feng et al., 2014).

TABLE 2

Table 2. Summary of selected vegetation indices published in the literature.

Results

Variability of Wheat Aerial N Uptake Under Different Growth Stages

The datasets from Exp. 1 are shown in Figure 6 to illustrate the general distribution of aerial N uptake. The aerial N uptake of the two wheat cultivars increased in the vegetative period because of increasing biomass. Across the different applied N rates, the aerial N uptake of Yumai 49–198 ranged from 6.0–23.5, 7.7–28.6, 8.1–32.2, and 8.9–35.1 g kg⁻¹ in the jointing, booting, heading, and anthesis stage, respectively. The aerial N uptake of Zhengmai 9694 varied from 6.9–24.3, 9.3–25.1, 10.6–27.1, and 13.4–31.2 g kg⁻¹, respectively, in these stages. With the progression of the growth stages, coefficients of variation for aerial N uptake increased. Thus, it can be seen that the aerial N uptake was significantly influenced by the different wheat cultivars and growth period.

FIGURE 6

Figure 6. Variation in plant N uptake over jointing-anthesis growth stages in wheat cultivars of Yumai 49-198 (A) and Zhengmai 9694 (B).

Relationship Between Canopy Reflectance and Aerial N Uptake at Different VZAs

We plotted the correlation between canopy reflectance and aerial N uptake under different VZAs with the data from Exp. 1–5 (Figure 7). In the 13 VZAs, a negative relationship was detected between aerial N uptake and the reflectance in the 400–720 nm. The minimum correlation coefficient was under 560–710 nm (r < −0.57), caused by red valley and chlorophyll absorption. The highest r-value was presented in the near infrared region (NIR), and increased with decreasing VZA both in backward and forward scattering. No matter in the backward and forward scattering, the r sharply changed in the red band region (690–760 nm), and it tended to 0 near 720 nm; this indicates that r from this region was not sensitive to the VZA.

FIGURE 7

Figure 7. The correlation coefficient (r) between reflectance (R_x) and aerial N uptake at 13 viewing zenith angles (A: backward scattering, B: forward scattering) (n = 155, p_0.001 = 0.262).

Relationship Between the Aerial N Uptake and Spectral Indices at Different VZAs

Canopy spectral data of different VZAs were influenced by numerous comprehensive factors including soil background, meteorological conditions, leaves, stems, and spectral noise. In this study, 18 spectral indices, including the new index and 17 published indices are presented in Table 2. Among them, the common indices were chose in accordance with a comprehensive literature investigation of red edge characteristics. To further describe the ability of the spectral indices in assessing aerial N uptake, we used the datasets from Exp. 1–5 to compare the predictive ability. As shown in Table 3 the r from backward scattering directions were higher than forward scattering directions for most two-band indices (except for DVI and RDVI). To the performance of 13 single VZAs, the relatively high correlation coefficient within −40° to +30° was present. For the well-performance VIs, the mRPA, DIDA, and DDn were advantageous at a viewing angle of −10°, with r scores of 0.896, −0.877, and −0.872, respectively;CI_red−edge3 and mSR705 were most sensitive at the −20° viewing angle, with r scores of 0.776 and 0.771, respectively; PSSRb and MSR demonstrated the best performance at the 0° viewing angle, with r scores of 0.808 and 0.789, respectively; SDr-SDb and TVI-3 had the higher correlations at the +10° viewing angle, with r scores of 0.814 and 0.809, respectively. Notably, eight of 18 indices produced the best correlations at the −10° viewing angle. These results illustrate that the around −10° VZA may be the most suitable for aerial N uptake estimation.

TABLE 3

Table 3. The correlation coefficient (r) for the relationships of vegetation indices with aerial N uptake at different viewing zenith angles (n = 155, p_0.001 = 0.262).

The four indices [mRPA, the two best-performing common indices (DIDA and DDn, and average (corresponding average value of 18 VIs at different VZAs, shown as average)] based on the R² values of the correlations between VIs and aerial N uptake were plotted in Figure 8. The results demonstrated that R² increased with decreasing VZA in both backward and forward scattering directions, and the highest R² were obtained under −20° to +10° VZAs. DIDA, DDn, and mRPA had strong correlations (R² > 0.72) to aerial N uptake in this region. Nevertheless, the average did not show any strong correlations with aerial N uptake (R² < 0.66). Compared with the average from −20° to +10° VZAs, the R² of mRPA, DIDA, and DDn was increased by 21.3–23.9, 15.3–18.8, and 14.4–17.6%, respectively.

FIGURE 8

Figure 8. Relationship between aerial N uptake and DIDA, DDn, mRPA, and average (corresponding average value of 18 VIs at different VZAs, shown as average) at different VZAs (n = 155, p_0.001 = 0.262).

Suitable Combined Angles for Aerial N Uptake Assessment Using VIs

To ascertain the suitable range of VIs to the mRPA, R² and RMSE were selected to compare different angle combinations. As showed in Figure 9, the performance of the R² and RMSE in the back scattering direction (R² = 0.426–0.804, p < 0.001) were superior to those in the forward view angles (R² = 0.337–0.775, p < 0.001). Among the 13 VZAs, −20° to +10° VZAs showed significantly higher predictive ability (high R² and small RMSE), and the most significant viewing angle was found to be −10°, with an R² and RMSE of 0.804 and 3.615, respectively. Figure 10 revealed that these combinations including large VZAs combinations generated relatively poor correlations. There was a dominant region in range of −10° to 0° (R² = 0.796), and the predictive ability of mRPA in −20° to +10° combination was also relatively high (R² = 0.782). When the VZA was out of −40° to +20°, the predictive accuracy was relatively low (R² < 0.740). A comparison of VIs (Figure 11) demonstrated that the mRPA (R² = 0.782 and 0.734) at −20° to +10° and −30° to +20° VZAs were more sensitive than the two best-performing published index DIDA (R² = 0.740 and 0.712) and DDn (R² = 0.726 and 0.701). Figure 12 showed the quantitative relationship between aerial N uptake and mRPA. The R² increased from 0.734 (from −30° to +20° combination) to 0.782 (from −20° to +10° combination). Compared to the −10° VZA having the highest R² value, mRPA in −20° to +10° combination only had a slightly decreased R² (2.7%) and an increased RMSE (5.2%). As a result, the novel mRPA model is the most forceful index for assessing aerial N uptake because of its insensitivity to VZAs of −20° to +10°, increasing the practicality of mRPA in actual production process.

FIGURE 9

Figure 9. Comparison of the prediction power of mRPA at 13 VZAs in terms of aerial N uptake (n = 155, p_0.001 = 0.262).

FIGURE 10

Figure 10. Comparison of the predictive ability (R²) of the indices within different view zenith angles combinations in terms of aerial N uptake.

FIGURE 11

Figure 11. Comparison of the predictive ability of the indices within five kinds of view zenith angles combinations in terms of aerial N uptake.

FIGURE 12

Figure 12. Comparison of the prediction power of mRPA at different VZA combinations in terms of aerial N uptake. (A) −30° to +20° and (B) −20° to +10°.

Testing Aerial N Uptake Estimation Models

The relationship between aerial N uptake and the spectral indices (across −20° to +10° VZAs) discussed above were measured utilizing data from Exp. 6 using R², RMSE, and RE to evaluate the accuracy and applicability between observed and predicted values. Data analysis was carried out on the common best-performing VIs DIDA and DDn, and on the novel index mRPA (Figure 13). The DIDA and DDn showed acceptable performance in the tests, with R² of 0.804, RE of 14.1%, and a RMSE of 2.464 for DIDA and R² of 0.792, RE of 14.9%, and a RMSE of 2.554 for DDn; this also indicates that DIDA is a better indicator than DDn. Compared with these two common indices, the mRPA demonstrated the superior predictive ability of aerial N uptake, with R² of 0.825, RE of 12.1% and RMSE of 2.190. In summary, mRPA seems to be the best index for predicting aerial N uptake of winter wheat under different management conditions.

FIGURE 13

Figure 13. Comparison between estimated and measured aerial N uptake based on DIDA (A), DDn (B), and mRPA (C) for −20° to +10° combinations (n = 252, p < 0.001).

Discussion

Remote sensing technology is widely used in agricultural production. It mainly includes the following aspects: crop growth measuring (e.g., biomass, N content, and yield), agricultural disaster monitoring (e.g., plant diseases and insect pests, droughts and floods) and crop phenotyping (e.g., crop height, leaf size, shape, and canopy longevity), and so on (Le Maire et al., 2008; Cao et al., 2015; He et al., 2016; Holman et al., 2016; Virlet et al., 2017), which could provide technical support for production management. However, spectral reflectance has previously been shown to be significantly affected by canopy structure, planting density, the wind, the angle of the sun, and various other factors (Rondeaux et al., 1996). To extract exact information for different characteristics and improve the detection precision, area calculations have been introduced to reduce the background effects (Broge and Leblanc, 2000; Ren et al., 2011; Li et al., 2013). Delegido et al. (2010) developed the normalized area over reflectance curve index (NAOC) and this index showed a linear relation with chlorophyll content. The normalized difference of the double-peak areas (NDDA) was calculated and discovered to correlate strongly with LNC (Feng et al., 2014). Broge and Leblanc (2000) found that the TVI (constructed by the area under the concave curve of red light absorption) could effectively assess the chlorophyll content and LAI. The chlorophyll absorption integral infers chlorophyll concentration through calculating the surrounding area between a connecting line of 600 and 735 nm and the red edge (Oppelt and Mauser, 2004). In this study, five of 17 VIs have previously been constructed using the area-based algorithm. Among them, DIDA, SDr-SDb, and TVI had significant correlations with aerial N uptake, with r-values of −0.877, 0.814, and 0.809, respectively, at their advantageous viewing angles (−10, 10, and 10°, respectively). The best-performing common index DIDA had strong correlations (R² > 0.74) to aerial N uptake within −20° to +10° VZAs. In addition, DIDA showed acceptable performance in the tests, with R² of 0.804, RE of 14.1% and a RMSE of 2.464. The above results showed that the VIs constructed by area algorithm could potentially be used to precisely predict aerial N uptake.

Many vegetation indices have been also constructed, including those NDVI-, RVI-, and DVI-like spectral indices, or other derived functions to enhance accuracy of estimating models (Huete, 1988; Wang et al., 2012). Baret and Guyot (1991) and Rondeaux et al. (1996) constructed the transformed soil-adjusted vegetation index (TSAVI) and optimized SAVI (OSAVI) by adding soil line parameters into NDVI to decrease the sensitivity of the soil background reflectance at low LAI. In order to reduce the combined influences of the canopy non-photosynthetic materials and increase the sensitivity of chlorophyll concentration determination, Daughtry et al. (2000) added the R₇₀₀/R₆₇₀ to the chlorophyll absorption ratio index (CARI) to obtain modified CARI. Haboudane et al. (2008) brought the (R₇₀₀/R₆₇₀)^1/2 into the TVI to construct the triangular chlorophyll index (TCI) to increase its sensitivity of chlorophyll changes. These studies inspired us to develop a novel index by adding a coefficient. We attempted to derive the coefficient by combining different bands in a square root form. Finally, R₇₆₀ and R₆₀₀ were selected from different combinations of bands and changed into the form of (R₇₆₀/R₆₀₀)^1/2. This was integrated with RSDR (760,718) to construct the mRPA, and this novel VI had a high correlation coefficient within −20° to +10° VZAs (R² = 0.782).The mRPA makes the best of area algorithms and red-edge information, and effectively improves the monitoring accuracy of aerial N uptake.

The vegetation indices displayed anisotropy depending on the canopy structural development, shadowing, the view angles of the sensors, the inherent viewing geometry of sensors, and in some respects the underlying soil (Kimes et al., 1985). Vegetation indices merging into multiple VZAs had the potency to further enhance anisotropy estimation. Pocewicz et al. (2007) took full advantage of the hotspot effect in the backscatter direction to improve quantitative estimation of LAI. Galvão et al. (2009) highlighted that the back-scattering direction was suited to predict the yield of soybean. He et al. (2016) showed that the novel VI constructed by a four-band VI from two angles (−20° and +10°) was sensitive to the change of the LNC in wheat. In this study, the back scattering direction improved indices performance for aerial N uptake compared with the forward scattering direction. The main reason may be that back scattering observations contain more signals from sunlit branches or leaves with higher reflectance values, while forward-scatter observations derive mostly from shady branches/leaves with lower reflectance values (Stagakis et al., 2010). In addition, our results showed that the R² increased with decreasing view zenith angles in back and forward-scatter direction. This was mainly because that the spectral data obtained at small angles mainly includes the total plant characteristics (lower, middle, and upper) of wheat. In general, sampling involving upper, middle, and lower wheat leaf layers could determine the aerial N uptake of the target region. In summary, further survey analysis on the variations of relationships between canopy attributes and remote sensing observations and on the availability off-nadir in fetching information is recommended.

Multi-angular remote sensing is able to obtain three-dimensional vegetation structure information, and thus it is better than vertical measurement for monitoring the canopy structural properties and the biochemical component of ground objects (Pocewicz et al., 2007). Rautiainen et al. (2008) demonstrated that high VZAs are the best fit for detecting over-story LAI values because of the quite limited influence of bottom layer on the whole signal. The value of canopy chlorophyll inversion index(CCII) at ±50 and ±60°, ±30 and ±40°, and nadir, ±20 and ±30°, VZAs were selected for inverting the chlorophyll at the upright upper, middle and bottom layer (Huang et al., 2011). Song et al. (2016) showed that −40° VZA was suitable to effectively monitor the LNC of wheat. This study showed that the R² of mRPA changed strongly with VZAs, which the highest R² value was found at −10° VZA (R² = 0.804). The mRPA within −10° could improve the monitoring accuracy of aerial N uptake. However, it is not convenient to accurately control measuring angle under off-nadir conditions in the field operation. Our research found that mRPA was relatively insensitive between −20° and +10° VZAs region for R² changes. A comparison among spectral parameters demonstrated that performance of mRPA in −20° to +10° combined dataset was superior to the two better published indices, with R² of 0.782 for mRPA, 0.750 for DIDA, and 0.736 for DDn, respectively. Compared to the most sensitive VZA (−10°), mRPA within −20° to +10° VZAs only had a slightly reduced R² (2.7%). This allowed us to construct a unified model in variable VZAs range to assess aerial N uptake in wheat, which decrease the influence of the VZAs and increase the field operation simplicity and applicability in a wide-angle region using portable monitors. Therefore, it is vital to select appropriate VI formulas and VZAs, which could reduce the variance due to soil background and crop canopy structure. In summary, the novel index plays an important role in predicting aerial N uptake of wheat and could be utilized to more precisely regulate N fertilization rate for different cultivation sites and plant types.

In this study, the novel index mRPA had the higher predictive ability with range of −20° to +10°. This can not only provide optimized parameters for the development of the portable monitor, but also offer dynamic information for guiding precise N budgeting. In order to obtain higher yield and avoid wasting resources, it is important to consider crop-N demand as well as soil-N supply to optimum N fertilizer strategy (Ju et al., 2009; Hartmann et al., 2015). The N fertilizer requirement was calculated using the following formula: N_req = (N_target – N_uptake – N_soil)/f_NUE, where N_target is the total crop-N demand for a target yield and grain protein, calculated according to Angus (2001), N_uptake is the aerial N uptake, N_soil is the potential soil-N supply for the rest of the growing period, and f_NUE is the fertilizer-N use efficiency. The mRPA model developed in this study could effectively estimate N_uptake, which will contribute to managing the N application in winter wheat. The prediction power of VIs was affected by cultivation factors. Only if the VI was seldom influenced by the factors of cultivation, the applicability of the model was strong. We synthesized dataset from the vegetative growth stages to develop a unified model that could be easily used to assess the N status. However, this research was designed only on winter wheat in Henan province, the dependability and adaptability of this novel model ought to be tested in other crops and areas.

Conclusions

Timely assessment of aerial N uptake is important to diagnose crop N status, maximizing yields and minimizing disadvantageous environmental impacts. In this study, we compared the use of 18 VIs, including 17 common VIs and a novel index constructed in this study, to estimate aerial N uptake of wheat. The results demonstrated that back scattering observation angles improved the ability to predict aerial N uptake compared with forward-scatter viewing angles VZAs. To decrease the restrictions on the environmental conditions and to further explore the superiority of spectral information, we combined the advantages of red-edge characteristics and area-based algorithms to construct a novel index mRPA to illuminate dynamic changes in aerial N uptake. The novel VI have the characteristic of simplicity and reliability and could be developed according to the formula: mRPA = (R₇₆₀/R₆₀₀)^1/2 × (R₇₆₀-R₇₁₈). Compared with the best-performing traditional indices DIDA and DDn, the predictive ability of mRPA at −10° view angle was effectively enhanced by 4.6–5.8%. Further systematic analysis of VZA combinations showed that mRPA had the best forecasting ability when compared with the traditional indices, with small difference between combinations of −20° to +10°. This has guided us in the development of a unified model for forecasting the aerial N uptake of wheat across a wide angle range; this will increase the precision of N predictions under a range of angles using portable monitors. The integrated index mRPA was shown to be practical and exact for aerial N uptake evaluation of winter wheat. This result will be also beneficial for choosing appropriate VZA and for the construction of more precise sensors for ecosystem monitoring. Nonetheless, it is also need to further validate the reliability and stability of the novel VI and to examine its effectiveness under the condition of different cultivation environment.

Author Contributions

B-BG, WF, and Y-JZ conceived the research. B-BG, Y-PW, YZ, and X-XR performed the experiments. B-BG and WF wrote the paper. LH and YM contributed to the results analysis and discussion.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Funding

This work was supported by the Thirteenth Five-year Plan of National Key Research Project of China [grant number 2016YFD0300604], the National Natural Science Foundation of China [grant number 31671624; 31571607], the Program of Science&Technology Innovation Talents in Universities of Henan Province, China [grant number 17HASTIT036].

Abbreviations

ASD, Analytical Spectral Devices; CARI, Chlorophyll absorption ratio index; Chl a, Chlorphyll a; CIred-edge3, Red-edge chlorophyll index-3; DDn, New double difference index; DIDA, Double-peak areas; DVI, Difference vegetation index; EVI, Enhanced Vegetation Index; LAI, Leaf area index; LNC, Leaf N concentration; LSDR, Left-side peak area; MCARI-1, Modified chlorophyll absorption in reflectance index; mRER, Modified red-edge ratio; mRPA, Modified right-side peak area index; MSR, Modified simple ratio; mSR705, Modified red-edge ratio; MTVI1, Modified triangular vegetation index; N, Nitrogen; NAOC, Normalized area over reflectance curve index; NDDA, Normalized difference of the double-peak areas; NDRE, Normalized difference red-edge index; NDVI, Normalized Difference Vegetation Index; NIR, Near infrared region; PSSRb, Pigment specific simple ratio chlorophyll b; R², Coefficients of determination; RDVI, Re-normalized difference vegetation index; RE, Relative error; REFCA, Red edge reflectance curve area; RMSE, Root mean square error; RSDR, Right-side peak area; SAVI, Soil-adjusted vegetation index; SPAD, Soil and Plant Analyzer Development; SPP, Solar principal plane; TCI, Triangular chlorophyll index; TVI, Triangle vegetation index; UAV, Unmanned Aerial Vehicle; VI, Vegetation index; VIopt, Optimal vegetation index; VZA, View zenith angle.

References

Andradesanchez, P., Gore, M. A., Heun, J. T., Thorp, K. R., Carmosilva, A. E., French, A. N., et al. (2014). Development and evaluation of a field-based high-throughput phenotyping platform. Funct. Plant Biol. 41, 68–79. doi: 10.1071/FP13126

CrossRef Full Text | Google Scholar

Angus, J. F. (2001). Nitrogen supply and demand in Australian agriculture. Aust. J. Exp. Agric. 41, 277–288. doi: 10.1071/EA00141

CrossRef Full Text | Google Scholar

Baret, F., and Guyot, G. (1991). Potentials and limits of vegetation indices for LAI and APAR assessment. Remote Sens. Environ. 35, 161–173. doi: 10.1016/0034-4257(91)90009-U

CrossRef Full Text | Google Scholar

Blackburn, G. A. (1998). Quantifying chlorophylls and caroteniods at leaf and canopy Scales. Remote Sens. Environ. 66, 273–285. doi: 10.1016/S0034-4257(98)00059-5

CrossRef Full Text | Google Scholar

Broge, N. H., and Leblanc, E. (2000). Comparing prediction power and stability of broadband and hyperspectral vegetation indices for estimation of green leaf area index and canopy chlorophyll density. Remote Sens. Environ. 76, 156–172. doi: 10.1016/S0034-4257(00)00197-8

CrossRef Full Text | Google Scholar

Cao, X., Luo, Y., Zhou, Y., Fan, J., Xu, X., West, J. S., et al. (2015). Detection of powdery mildew in two winter wheat plant densities and prediction of grain yield using canopy hyperspectral reflectance. PLoS ONE 10:e121462. doi: 10.1371/journal.pone.0121462

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J. M. (2014). Evaluation of vegetation indices and a modified simple ratio for boreal applications. Can. J. Remote Sens. 22, 229–242. doi: 10.1080/07038992.1996.10855178

CrossRef Full Text | Google Scholar

Cho, M. A., and Skidmore, A. K. (2006). A new technique for extracting the red edge position from hyperspectral data: the linear extrapolation method. Remote Sens. Environ. 101, 181–193. doi: 10.1016/j.rse.2005.12.011

CrossRef Full Text | Google Scholar

Daughtry, C. S. T., Walthall, C. L., Kim, M. S., Colstoun, E. B. D., and Iii, M. M. (2000). Estimating corn leaf chlorophyll concentration from leaf and canopy reflectance. Remote Sens. Environ. 74, 229–239. doi: 10.1016/S0034-4257(00)00113-9

CrossRef Full Text | Google Scholar

Delegido, J., Alonso, L., González, G., and Moreno, J. (2010). Estimating chlorophyll content of crops from hyperspectral data using a normalized area over reflectance curve (NAOC). Int. J. Appl. Earth Obs. Geoinf. 12, 165–174. doi: 10.1016/j.jag.2010.02.003

CrossRef Full Text | Google Scholar

Diner, D. J., Asner, G. P., Davies, R., Knyazikhin, Y., Muller, J. P., Nolin, A. W., et al. (1999). New directions in earth observing: scientific applications of multiangle remote sensing. Bull. Am. Meteorol. Soc. 80, 2209–2228. doi: 10.1175/1520-0477(1999)080<2209:NDIEOS>2.0.CO;2

CrossRef Full Text | Google Scholar

Ecarnot, M., Compan, F., and Roumet, P. (2013). Assessing leaf nitrogen content and leaf mass per unit area of wheat in the field throughout plant cycle with a portable spectrometer. Field Crops Res. 140, 44–50. doi: 10.1016/j.fcr.2012.10.013

CrossRef Full Text | Google Scholar

Erdle, K., Mistele, B., and Schmidhalter, U. (2011). Comparison of active and passive spectral sensors in discriminating biomass parameters and nitrogen status in wheat cultivars. Field Crops Res. 124, 74–84. doi: 10.1016/j.fcr.2011.06.007

CrossRef Full Text | Google Scholar

Feng, W., Guo, B.-B., Wang, Z. J., He, L., Song, X., Wang, Y. H., et al. (2014). Measuring leaf nitrogen concentration in winter wheat using double-peak spectral reflection remote sensing data. Field Crops Res. 159, 43–52. doi: 10.1016/j.fcr.2014.01.010

CrossRef Full Text | Google Scholar

Feng, W., Guo, B. B., Zhang, H. Y., He, L., Zhang, Y. S., Wang, Y. H., et al. (2015). Remote estimation of above ground nitrogen uptake during vegetative growth in winter wheat using hyperspectral red-edge ratio data. Field Crops Res. 180, 197–206. doi: 10.1016/j.fcr.2015.05.020

CrossRef Full Text | Google Scholar

Feng, W., Zhu, Y., Tian, Y. C., Cao, W. X., Yao, X., and Li, Y. X. (2008). Monitoring leaf nitrogen accumulation with hyper-spectral remote sensing in wheat. Acta Ecol. Sin. 28, 23–32. doi: 10.1016/S1872-2032(08)60018-9

CrossRef Full Text | Google Scholar

Fitzgerald, G. J., Rodriguez, D., Christensen, L. K., and Belford, R. (2006). Spectral and thermal sensing for nitrogen and water status in rainfed and irrigated wheat environments. Precis. Agric. 7, 233–248. doi: 10.1007/s11119-006-9011-z

CrossRef Full Text | Google Scholar

Galvão, L. S., Roberts, D. A., Formaggio, A. R., Numata, I., and Breunig, F. M. (2009). View angle effects on the discrimination of soybean varieties and on the relationships between vegetation indices and yield using off-nadir Hyperion data. Remote Sens. Environ. 113, 846–856. doi: 10.1016/j.rse.2008.12.010

CrossRef Full Text | Google Scholar

Gemmell, F., and McDonald, A. J. (2000). View zenith angle effects on the forest information content of three spectral indices. Remote Sens. Environ. 72, 139–158. doi: 10.1016/S0034-4257(99)00086-3

CrossRef Full Text | Google Scholar

Gitelson, A. A., Viña, A., Ciganda, V., Rundquist, D. C., and Arkebauer, T. J. (2005). Remote estimation of canopy chlorophyll content in crops. Geophys. Res. Lett. 32, 93–114. doi: 10.1029/2005GL022688

CrossRef Full Text | Google Scholar

Haboudane, D., Miller, J. R., Pattey, E., Zarco-Tejada, P. J., and Strachan, I. B. (2004). 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. 90, 337–352. doi: 10.1016/j.rse.2003.12.013

CrossRef Full Text | Google Scholar

Haboudane, D., Tremblay, N., Miller, J. R., and Vigneault, P. (2008). Remote estimation of crop chlorophyll content using spectral indices derived from hyperspectral data. IEEE Trans. Geosci. Remote Sens. 46, 423–437. doi: 10.1109/TGRS.2007.904836

CrossRef Full Text | Google Scholar

Hansen, P. M., and Schjoerring, J. K. (2003). Reflectance measurement of canopy biomass and nitrogen status in wheat crops using normalized difference vegetation indices and partial least squares regression. Remote Sens. Environ. 86, 542–553. doi: 10.1016/S0034-4257(03)00131-7

CrossRef Full Text | Google Scholar

Hartmann, T. E., Yue, S. C., Schulz, R., He, X. K., Chen, X. P., Zhang, F. S., et al. (2015).Yield and N use efficiency of a maize–wheat cropping system as affected by different fertilizer management strategies in a farmer's field of the North China Plain. Field Crops Res. 174, 30–39. doi: 10.1016/j.fcr.2015.01.006

CrossRef Full Text

Hasegawa, K., Matsuyama, H., Tsuzuki, H., and Sweda, T. (2010). Improving the estimation of leaf area index by using remotely sensed NDVI with BRDF signatures. Remote Sens. Environ. 114, 514–519. doi: 10.1016/j.rse.2009.10.005

CrossRef Full Text | Google Scholar

Hatfield, J. L., Gitelson, A. A., Schepers, J. S., and Walthall, C. L. (2008). Application of spectral remote sensing for agronomic decisions. Agron. J. 100, 117–131. doi: 10.2134/agronj2006.0370c

CrossRef Full Text | Google Scholar

He, L., Zhang, H. Y., Zhang, Y. S., Song, X., Feng, W., Kang, G. Z., et al. (2016). Estimating canopy leaf nitrogen concentration in winter wheat based on multi-angular hyperspectral remote sensing. Eur. J. Agron. 73, 170–185. doi: 10.1016/j.eja.2015.11.017

CrossRef Full Text | Google Scholar

Holman, F. H., Riche, A. B., Michalski, A., Castle, M., Wooster, M. J., and Hawkesford, M. J. (2016). High throughput field phenotyping of wheat plant height and growth rate in field plot trials using UAV based remote sensing. Remote Sens. 8, 1–23. doi: 10.3390/rs8121031

CrossRef Full Text | Google Scholar

Huang, W., Wang, Z., Huang, L., and Lamb, D. W. (2011). Estimation of vertical distribution of chlorophyll concentration by bi-directional canopy reflectance spectra in winter wheat. Precis. Agric. 12, 165–178. doi: 10.1007/s11119-010-9166-5

CrossRef Full Text | Google Scholar

Huete, A. R. (1988). A Soil-Adjusted Vegetation Index (SAVI). Remote Sens. Environ., 25, 295–309. doi: 10.1016/0034-4257(88)90106-X

CrossRef Full Text | Google Scholar

Inoue, Y., Penuelas, J., Miyata, A., and Mano, M. (2008). Normalized difference spectral indices for estimating photosynthetic efficiency and capacity at a canopy scale derived from hyperspectral and CO₂ flux measurements in rice. Remote Sens. Environ. 112, 156–172. doi: 10.1016/j.rse.2007.04.011

CrossRef Full Text | Google Scholar

Isaac, R. A., and Johnson, W. C. (1976). Determination of total nitrogen in plant tissue using a bloc digestor. Assoc. Offic. Analyt. Chem. 59, 98–100.

Ju, X. T., Kou, C. L., Zhang, F. S., and Christie, P. (2006). Nitrogen balance and groundwater nitrate contamination: comparison among three intensive cropping systems on the North China Plain. Environ. Pollut. 143, 117–125. doi: 10.1016/j.envpol.2005.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Ju, X. T., Xing, G. X., Chen, X. P., Zhang, S. L., Zhang, L. J., Liu, X. J., et al. (2009). Reducing environmental risk by improving N management in intensive Chinese agricultural systems. Proc. Natl. Acad. Sci. U.S.A. 106, 3046–3051. doi: 10.1073/pnas.0813417106

PubMed Abstract | CrossRef Full Text | Google Scholar

Kimes, D. S., Newcomb, W. W., Tucker, C. J., Zonneveld, I. S., Wijngaarden, W. V., Leeuw, J. D., et al. (1985). Directional reflectance factor distributions for cover types of Northern Africa. Remote Sens. Environ. 18, 1–19. doi: 10.1016/0034-4257(85)90034-3

CrossRef Full Text | Google Scholar

le Maire, G., François, C., Soudani, K., Berveiller, D., Pontailler, J. Y., Bréda, N., et al. (2008). Calibration and validation of hyperspectral indices for the estimation of broadleaved forest leaf chlorophyll content, leaf mass per area, leaf area index and leaf canopy biomass. Remote Sens. Environ. 112, 3846–3864. doi: 10.1016/j.rse.2008.06.005

CrossRef Full Text | Google Scholar

Li, F., Mistele, B., Hu, Y., Chen, X., and Schmidhalter, U. (2013). Comparing hyperspectral index optimization algorithms to estimate aerial N uptake using multi-temporal winter wheat datasets from contrasting climatic and geographic zones in China and Germany. Agric. For. Meteorol. 180, 44–57. doi: 10.1016/j.agrformet.2013.05.003

CrossRef Full Text | Google Scholar

Mistele, B., and Schmidhalter, U. (2008). Estimating the nitrogen nutrition index using spectral canopy reflectance measurements. Eur. J. Agron. 29, 184–190. doi: 10.1016/j.eja.2008.05.007

CrossRef Full Text | Google Scholar

Myneni, R. B., Maggion, S., Iaquinta, J., Privette, J. L., Gobron, N., Pinty, B., et al. (1995). Optical remote sensing of vegetation: modeling, caveats, and algorithms. Remote Sens. Environ. 51, 169–188. doi: 10.1016/0034-4257(94)00073-V

CrossRef Full Text | Google Scholar

Oppelt, N., and Mauser, W. (2004). Hyperspectral monitoring of physiological parameters of wheat during a vegetation period using AVIS data. Int. J. Remote Sens. 25, 145–159. doi: 10.1080/0143116031000115300

CrossRef Full Text | Google Scholar

Pocewicz, A., Vierling, L. A., Lentile, L. B., and Smith, R. (2007). View angle effects on relationships between MISR vegetation indices and leaf area index in a recently burned ponderosa pine forest. Remote Sens. Environ. 107, 322–333. doi: 10.1016/j.rse.2006.06.019

CrossRef Full Text | Google Scholar

Rautiainen, M., Lang, M., Mõttus, M., Kuusk, A., Nilson, T., Kuusk, J., et al. (2008). Multi-angular reflectance properties of a hemiboreal forest: an analysis using CHRIS PROBA data. Remote Sens. Environ. 112, 2627–2642. doi: 10.1016/j.rse.2007.12.005

CrossRef Full Text | Google Scholar

Ren, H., Zhou, G., and Zhang, X. (2011). Estimation of green aboveground biomass of desert steppe in Inner Mongolia based on red-edge reflectance curve area method. Biosyst. Eng. 109, 385–395. doi: 10.1016/j.biosystemseng.2011.05.004

CrossRef Full Text | Google Scholar

Reyniers, M., Walvoort, D. J. J., and Baardemaaker, J. D. (2006). A linear model to predict with a multi-spectral radiometer the amount of nitrogen in winter wheat. Int. J. Remote Sens. 27, 4159–4178. doi: 10.1080/01431160600791650

CrossRef Full Text | Google Scholar

Richardson, A. J., and Wiegand, C. L. (1977). Distinguishing Vegetation from Soil background information. Photogramm. Eng. Rem. Sens. 43, 1541–1552.

Google Scholar

Rondeaux, G., Steven, M., and Baret, F. (1996). Optimization of soil-adjusted vegetation indices. Remote Sens. Environ. 55, 95–107. doi: 10.1016/0034-4257(95)00186-7

CrossRef Full Text | Google Scholar

Roujean, J. L., and Breon, F. M. (1995). Estimating PAR absorbed by vegetation from bidirectional reflectance measurements. Remote Sens. Environ. 51, 375–384. doi: 10.1016/0034-4257(94)00114-3

CrossRef Full Text | Google Scholar

Sandmeier, S. R., and Itten, K. I. (1999). A field goniometer system (FIGOS) for acquisition of hyperspectral BRDF data. IEEE Trans. Geosci. Rem. Sens. 37, 978–986. doi: 10.1109/36.752216

CrossRef Full Text | Google Scholar

Sehy, U., Ruser, R., and Munch, J. C. (2003). Nitrous oxide fluxes from maize fields: relationship to yield, site-specific fertilization, and soil conditions. Agric. Ecosyst. Environ. 99, 97–111. doi: 10.1016/S0167-8809(03)00139-7

CrossRef Full Text | Google Scholar

Shibayama, M., and Wiegand, C. L. (1985). View azimuth and zenith, and solar angle effects on wheat canopy reflectance. Remote Sens. Environ. 18, 91–103. doi: 10.1016/0034-4257(85)90040-9

CrossRef Full Text | Google Scholar

Sims, D. A., and Gamon, J. A. (2002). Relationships between leaf pigment content and spectral reflectance across a wide range of species, leaf structures and developmental stages. Remote Sens. Environ. 81, 337–354. doi: 10.1016/S0034-4257(02)00010-X

CrossRef Full Text | Google Scholar

Smil, V. (2002). Nitrogen and food production: proteins for human diets. Ambio 31, 126–131. doi: 10.1579/0044-7447-31.2.126

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, X., Xu, D. Y., He, L., Feng, W., Wang, Y. H., Wang, Z. J., et al. (2016).Using multi-angle hyperspectral data to monitor canopy leaf nitrogen content of wheat. Precis. Agric. 17, 721–736. doi: 10.1007/s11119-016-9445-x

CrossRef Full Text

Stagakis, S., Markos, N., Sykioti, O., and Kyparissis, A. (2010). Monitoring canopy biophysical and biochemical parameters in ecosystem scale using satellite hyperspectral imagery: an application on a Phlomis fruticosa Mediterranean ecosystem using multiangular CHRIS/PROBA observations. Remote Sens. Environ. 114, 977–994. doi: 10.1016/j.rse.2009.12.006

CrossRef Full Text | Google Scholar

Thenkabail, P. S., Smith, R. B., and Pauw, E. D. (2000). Hyperspectral vegetation indices and their relationships with agricultural crop characteristics. Remote Sens. Environ. 71, 158–182. doi: 10.1016/S0034-4257(99)00067-X

CrossRef Full Text | Google Scholar

Virlet, N., Sabermanesh, K., Sadeghi-Tehran, P., and Hawkesford, M. J. (2017). Field scanalyzer: an automated robotic field phenotyping platform for detailed crop monitoring. Funct. Plant Biol. 44, 143–153. doi: 10.1071/FP16163

CrossRef Full Text | Google Scholar

Wang, W., Yao, X., Yao, X. F., Tian, Y. C., Liu, X. J., Ni, J., et al. (2012). Estimating leaf nitrogen concentration with three-band vegetation indices in rice and wheat. Field Crops Res. 129, 90–98. doi: 10.1016/j.fcr.2012.01.014

CrossRef Full Text | Google Scholar

Woodard, H. J., and Bly, A. (1998). Relationship of nitrogen management to winter wheat yield and grain protein in South Dakota. J. Plant Nutr. 21, 217–233. doi: 10.1080/01904169809365397

CrossRef Full Text | Google Scholar