Original Research ARTICLE
Remotely Estimating Aerial N Uptake in Winter Wheat Using Red-Edge Area Index From Multi-Angular Hyperspectral Data
- State Key Laboratory of Wheat and Maize Crop Science, National Engineering Research Centre for Wheat, Henan Agricultural University, Zhengzhou, China
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 (R2) 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 = (R760/R600)1/2 × (R760-R718). Investigating the predictive accuracy of mRPA for aerial N uptake across VZAs demonstrated that the best performance was at −10° [R2 = 0.804, p < 0.001, root mean square error (RMSE) = 3.615] and that the effect was relatively similar between −20° to +10° (R2 = 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.
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 (N2O) 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 × R423 band to the NDVI (R703, R924) 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 (R759-1.8 × R419)/(R742-1.8 × R419) 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
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−1 P2O5 [as Ca(H2PO4)2] and 90 kg ha−1 K2O (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 × 106 plants ha−1.
Measurement of Canopy Multi-Angular Hyperspectral Reflectance
Canopy reflectance spectra were obtained in a 1 m2 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 BaSO4.
The areas of 0.20 m2 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 (R760-R718) divided by characteristic wavelength method had the best performance (R2 = 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 R801/R670, R801/R550) 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 (R700/R670)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)1/2] into RSDR (R760-R718) to construct a novel VI called the mRPA. The spectral region of the above λ1 and λ2 were located within 400 and 900 nm. Figure 5 gave a comprehensive overview of the correlation coefficients for any two band combinations (λ1, λ2), and this is valid for selecting the sensitive bands of aerial N uptake. The wavebands λ1 and λ2 ranged between 750–900 and 550–650 nm, respectively. This area had the highest precision, with R2 values above 0.78, especially R2 performed best when λ1 = 760 nm, λ2 = 600 nm. Therefore, the final formula of mRPA was:
Figure 3. The division of the whole red-edge double-peak area into left and right single peak areas.
Figure 5. Contour maps of coefficients of determination between aerial N uptake and the mRPA with formula (Rλ1/Rλ2)1/2 × (R760-R718) (n = 155, p0.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 R2. The model's behavior was evaluated by employing the R2, RMSE, and RE (%). Among the indices, those with the highest R2 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:
Here, Pi, and Oi 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).
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−1 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−1, 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. 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. The correlation coefficient (r) between reflectance (Rx) and aerial N uptake at 13 viewing zenith angles (A: backward scattering, B: forward scattering) (n = 155, p0.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;CIred−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. The correlation coefficient (r) for the relationships of vegetation indices with aerial N uptake at different viewing zenith angles (n = 155, p0.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 R2 values of the correlations between VIs and aerial N uptake were plotted in Figure 8. The results demonstrated that R2 increased with decreasing VZA in both backward and forward scattering directions, and the highest R2 were obtained under −20° to +10° VZAs. DIDA, DDn, and mRPA had strong correlations (R2 > 0.72) to aerial N uptake in this region. Nevertheless, the average did not show any strong correlations with aerial N uptake (R2 < 0.66). Compared with the average from −20° to +10° VZAs, the R2 of mRPA, DIDA, and DDn was increased by 21.3–23.9, 15.3–18.8, and 14.4–17.6%, respectively.
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, p0.001 = 0.262).
Suitable Combined Angles for Aerial N Uptake Assessment Using VIs
To ascertain the suitable range of VIs to the mRPA, R2 and RMSE were selected to compare different angle combinations. As showed in Figure 9, the performance of the R2 and RMSE in the back scattering direction (R2 = 0.426–0.804, p < 0.001) were superior to those in the forward view angles (R2 = 0.337–0.775, p < 0.001). Among the 13 VZAs, −20° to +10° VZAs showed significantly higher predictive ability (high R2 and small RMSE), and the most significant viewing angle was found to be −10°, with an R2 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° (R2 = 0.796), and the predictive ability of mRPA in −20° to +10° combination was also relatively high (R2 = 0.782). When the VZA was out of −40° to +20°, the predictive accuracy was relatively low (R2 < 0.740). A comparison of VIs (Figure 11) demonstrated that the mRPA (R2 = 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 (R2 = 0.740 and 0.712) and DDn (R2 = 0.726 and 0.701). Figure 12 showed the quantitative relationship between aerial N uptake and mRPA. The R2 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 R2 value, mRPA in −20° to +10° combination only had a slightly decreased R2 (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. Comparison of the prediction power of mRPA at 13 VZAs in terms of aerial N uptake (n = 155, p0.001 = 0.262).
Figure 10. Comparison of the predictive ability (R2) of the indices within different view zenith angles combinations in terms of aerial N uptake.
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. 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 R2, 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 R2 of 0.804, RE of 14.1%, and a RMSE of 2.464 for DIDA and R2 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 R2 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. 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).
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 (R2 > 0.74) to aerial N uptake within −20° to +10° VZAs. In addition, DIDA showed acceptable performance in the tests, with R2 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 R700/R670 to the chlorophyll absorption ratio index (CARI) to obtain modified CARI. Haboudane et al. (2008) brought the (R700/R670)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, R760 and R600 were selected from different combinations of bands and changed into the form of (R760/R600)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 (R2 = 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 R2 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 R2 of mRPA changed strongly with VZAs, which the highest R2 value was found at −10° VZA (R2 = 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 R2 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 R2 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 R2 (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: Nreq = (Ntarget – Nuptake – Nsoil)/fNUE, where Ntarget is the total crop-N demand for a target yield and grain protein, calculated according to Angus (2001), Nuptake is the aerial N uptake, Nsoil is the potential soil-N supply for the rest of the growing period, and fNUE is the fertilizer-N use efficiency. The mRPA model developed in this study could effectively estimate Nuptake, 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.
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 = (R760/R600)1/2 × (R760-R718). 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.
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.
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].
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; R2, 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.
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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 CO2 flux measurements in rice. Remote Sens. Environ. 112, 156–172. doi: 10.1016/j.rse.2007.04.011
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
Keywords: winter wheat, multi-angular hyperspectral, vegetation indices, aerial N uptake, monitoring model
Citation: Guo B-B, Zhu Y-J, Feng W, He L, Wu Y-P, Zhou Y, Ren X-X and Ma Y (2018) Remotely Estimating Aerial N Uptake in Winter Wheat Using Red-Edge Area Index From Multi-Angular Hyperspectral Data. Front. Plant Sci. 9:675. doi: 10.3389/fpls.2018.00675
Received: 30 November 2017; Accepted: 03 May 2018;
Published: 25 May 2018.
Edited by:Tracy Ann Porcelli, Independent Researcher, Canada
Reviewed by:Sung Kyeom Kim, Rural Development Administration, South Korea
Andrew Davidson, Government of Canada, Canada
Copyright © 2018 Guo, Zhu, Feng, He, Wu, Zhou, Ren and Ma. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.