Estimation of Winter Wheat Above-Ground Biomass Using Unmanned Aerial Vehicle-Based Snapshot Hyperspectral Sensor and Crop Height Improved Models

Correct estimation of above-ground biomass (AGB) is necessary for accurate crop growth monitoring and yield prediction. We estimated AGB based on images obtained with a snapshot hyperspectral sensor (UHD 185 firefly, Cubert GmbH, Ulm, Baden-Württemberg, Germany) mounted on an unmanned aerial vehicle (UAV). The UHD 185 images were used to calculate the crop height and hyperspectral reflectance of winter wheat canopies from hyperspectral and panchromatic images. We constructed several single-parameter models for AGB estimation based on spectral parameters, such as specific bands, spectral indices (e.g., Ratio Vegetation Index (RVI), NDVI, Greenness Index (GI) and Wide Dynamic Range VI (WDRVI)) and crop height and several models combined with spectral parameters and crop height. Comparison with experimental results indicated that incorporating crop height into the models improved the accuracy of AGB estimations (the average AGB is 6.45 t/ha). The estimation accuracy of single-parameter models was low (crop height only: R2 = 0.50, RMSE = 1.62 t/ha, MAE = 1.24 t/ha; R670 only: R2 = 0.54, RMSE = 1.55 t/ha, MAE = 1.23 t/ha; NDVI only: R2 = 0.37, RMSE = 1.81 t/ha, MAE = 1.47 t/ha; partial least squares regression R2 = 0.53, RMSE = 1.69, MAE = 1.20), but accuracy increased when crop height and spectral parameters were combined (partial least squares regression modeling: R2 = 0.78, RMSE = 1.08 t/ha, MAE = 0.83 t/ha; verification: R2 = 0.74, RMSE = 1.20 t/ha, MAE = 0.96 t/ha). Our results suggest that crop height determined from the new UAV-based snapshot hyperspectral sensor can improve AGB estimation and is advantageous for mapping applications. This new method can be used to guide agricultural management.


Introduction
Above-ground biomass (AGB) [1] reflects the status of crop growth and is related to solar energy utilization, yield and grain quality [2,3].Therefore, AGB is considered one of the most important biochemical crop parameters [4][5][6], and correct estimation of AGB can help improve crop monitoring and yield prediction [7].Estimation of AGB is also important for ecological research at regional and global scales [8] and has been employed in global carbon cycle [9][10][11] and climate change [12][13][14] studies.
Manual methods for measuring AGB are time consuming and destructive.Prior to the emergence of remote sensing technology, large AGB surveys were limited by labor and resources [4][5][6].
Hyperspectral technology has shown great potential for monitoring of crop parameters [15,16].Hyperspectral sensors can obtain spectral features based on radiation from visible to near-infrared wavelengths, and advancements in technology have increased the number of features that can be collected [17][18][19].Spectral features and indices are important variables for estimating hyperspectral crop parameters [20][21][22][23][24][25][26], and many studies have estimated AGB using a combination of narrow spectral bands and vegetation indices [27][28][29].
The development and application of unmanned aerial vehicles (UAVs) as remote sensing platforms [15,[30][31][32] have increased rapidly over the past decade.Additionally, UAVs can obtain remote sensing images with higher temporal, spatial and ground resolution than satellites.Small-sized remote sensing sensors can be mounted on UAVs, making UAV platforms increasingly popular in agricultural applications [15,31].Use of UAVs in agricultural and environmental applications has numerous advantages compared with conventional aerial vehicles, including lower cost, weight, flight speed and flight altitude [15].
Previous studies have proposed various methods and sensors for agricultural monitoring based on UAV remote sensing.The results of Ribeiro et al. [33] (geometric products), Bendig et al. [34,35] (plant height, biomass) and Jannoura et al. [36] (LAI) indicated that digital cameras with visible and infrared spectra are efficient and flexible for monitoring agricultural crop parameters.Berni et al. [37] suggested that thermal infrared cameras mounted on UAVs show potential for environmental and agricultural applications (mapping canopy conductance and crop water stress).The results of Vonbueren et al. [38], Vega et al. [39] and Honkavaara et al. [40] indicated that data from UAV-based spectral cameras can be used to monitor parameters (e.g., crop height, yield, aboveground biomass and nitrogen content) of various plants, such as grass, wheat and sunflowers.However, the use of UAVs for agricultural monitoring is limited by the weight of hyperspectral imaging systems, the complexity of image processing and the cost of sensors [41,42].
Many powerful techniques/methods had been employed by previous researchers to make full use of the hyperspectral information, such as partial least squares regression (PLSR), support vector machine regression (SVR), artificial neural network (ANN) and random forest regression (RF).The results of Jin et al. [25] (leaf water content), Fu et al. [43] (biomass) and Nguyen et al. [44] (rice leaf growth and nitrogen status) indicated that PLSR is efficient and flexible for hyperspectral monitoring of agricultural crop parameters.Yang et al. [45] (LAI) and Atzberger et al. [46] (LAI, leaf chlorophyll content, equivalent leaf water thickness) suggested that ANN has great potential for biophysical canopy variables' estimation.Yang et al. [47] suggested that SVM shows great potential for rice leaf area index estimation.Han et al. [48] suggested that RF is more efficient than SVM in apple tree canopy LAI estimation.The results of Yuan et al. [49] indicated that RF can be used to monitor the soybean leaf area index over the whole growth period.
The objective of this study was to evaluate the performance of AGB estimations based on hyperspectral data (canopy surface information) and crop height (vertical information about crop growth).A snapshot hyperspectral sensor mounted on a UAV was used to obtain hyperspectral and high-resolution panchromatic images.Using image fusion and photogrammetry, we obtained the canopy hyperspectral reflectance and crop height of winter wheat.Hyperspectral data and crop height were used to construct AGB estimation models, and model accuracy was validated with field data.

Materials and Methods
Winter wheat was grown on 48 plots with different amounts of nitrogen fertilizer and irrigation levels.We measured spectral features of the crops and analyzed the relationships between crop height, UHD 185 (UHD 185 firefly, Cubert GmbH, Ulm, Baden-Württemberg, Germany) spectrometer hyperspectral data, ASD Field Spec 3 spectrometer hyperspectral data and winter wheat AGB.

Materials and Methods
Winter wheat was grown on 48 plots with different amounts of nitrogen fertilizer and irrigation levels.We measured spectral features of the crops and analyzed the relationships between crop height, UHD 185 (UHD 185 firefly, Cubert GmbH, Ulm, Baden-Württemberg, Germany) spectrometer hyperspectral data, ASD Field Spec 3 spectrometer hyperspectral data and winter wheat AGB.

Study Area
Experiments were conducted at the Xiao Tangshan National Precision Agriculture Research Center of China, which has an area of ~2 square kilometers and is located in Changping District (115°50′17″～116°29′49″E, 40°2′18″～40°23′13″N) of Beijing City, China (Figure 1).Changping District has an average altitude of 36 m, a warm temperate semi-humid continental monsoon climate, average rainfall of 450 mm, average low temperature of −10 °C-7.5 °C and an average high temperature of 35 °C-40 °C.Meteorological data were acquired from the China Meteorological Data Service [50].

AGB and Crop Height Measurement
The sampling plan is shown in Figure 2. Crops were divided into three groups, and each group included a total of 16 winter wheat samples.Samples received different levels of irrigation and nitrogen fertilizer.The irrigation levels included rainfall only (W0), rainfall plus the normal irrigation level (W1, 100 mm) and rainfall plus twice the normal irrigation level (W2, 200 mm).The nitrogen fertilizer levels included no fertilizer (N0), one-half the normal fertilization level (N1, 195 kg/ha), the normal fertilization level (N2, 390 kg/ha) and twice the normal fertilization level (N3, 780 kg/ha).
During sampling, 20 stems were collected near the center of each plot.Crop height was measured using a straightedge.Samples were taken to the laboratory for processing and stored in paper bags.Leaves and stems were dried at 80 °C to remove moisture.Until the samples became constant weight (about 24 h), the dry weight of each sample was recorded.Winter wheat AGB was calculated using Equation ( 1) where m is the dry weight of the sample, n is the number of winter wheat ears per unit area and l is the row spacing.

AGB and Crop Height Measurement
The sampling plan is shown in Figure 2. Crops were divided into three groups, and each group included a total of 16 winter wheat samples.Samples received different levels of irrigation and nitrogen fertilizer.The irrigation levels included rainfall only (W0), rainfall plus the normal irrigation level (W1, 100 mm) and rainfall plus twice the normal irrigation level (W2, 200 mm).The nitrogen fertilizer levels included no fertilizer (N0), one-half the normal fertilization level (N1, 195 kg/ha), the normal fertilization level (N2, 390 kg/ha) and twice the normal fertilization level (N3, 780 kg/ha).
During sampling, 20 stems were collected near the center of each plot.Crop height was measured using a straightedge.Samples were taken to the laboratory for processing and stored in paper bags.Leaves and stems were dried at 80 • C to remove moisture.Until the samples became constant weight (about 24 h), the dry weight of each sample was recorded.Winter wheat AGB was calculated using Equation ( 1) where m is the dry weight of the sample, n is the number of winter wheat ears per unit area and l is the row spacing.

Snapshot Hyperspectral Sensor
The UHD 185 Firefly (UHD 185) is a new type of snapshot hyperspectral sensor, characterized by short exposure and integration time.The sensor weighs 0.47 kg, and its dimensions are 195 × 67 × 60 mm.It can obtain reflected radiation from the visible to near-infrared spectra.The sensor can capture 450-950-nm wavelengths, while maintaining a reasonable balance of image and good spectral resolution (4 nm).Radiation is recorded and processed as a 1000 × 1000 (1 band) panchromatic image and a 50 × 50 (125 bands) hyperspectral image.The acquired panchromatic images are rich in texture information, and stitching is relatively simple; however, images lack spectral information.The 50 × 50 (125 bands) hyperspectral images are characterized by rich spectral information, but lack texture information.Hyperspectral image fusion [19] was used to splice the UHD 185 hyperspectral image.After fusion, all images were stitched together using an image stitching process.

Platform
The UAV UHD 185 hyperspectral system was designed by the National Engineering Research Center for Information Technology in Agriculture (NERCITA) of China [51].The system was composed of three parts, including a ground control system, UAV and UHD 185 hyperspectral system and post-processing software.As shown in Figure 3, the ground control system included a DJI S1000 UAV (SZ DJI Technology Co., Ltd., Sham Chun, Guangdong Province, China) [52] (Figure 3), the DJI ground-control station (software, Figure 3a), reinforcement notebook computer (hardware, Figure 3a) and remote controller (Figure 3b).The UHD 185 hyperspectral system consisted of a cloud platform and the UHD 185 hyperspectral sensor (Figure 3c).The post-processing software consisted of software for stitching, as well as radiation and geometry correction (e.g., as Cubert [19], Agisoft PhotoScan [53] and ENVI 5.1 [54]).

Snapshot Hyperspectral Sensor
The UHD 185 Firefly (UHD 185) is a new type of snapshot hyperspectral sensor, characterized by short exposure and integration time.The sensor weighs 0.47 kg, and its dimensions are 195 × 67 × 60 mm.It can obtain reflected radiation from the visible to near-infrared spectra.The sensor can capture 450-950-nm wavelengths, while maintaining a reasonable balance of image and good spectral resolution (4 nm).Radiation is recorded and processed as a 1000 × 1000 (1 band) panchromatic image and a 50 × 50 (125 bands) hyperspectral image.The acquired panchromatic images are rich in texture information, and stitching is relatively simple; however, images lack spectral information.The 50 × 50 (125 bands) hyperspectral images are characterized by rich spectral information, but lack texture information.Hyperspectral image fusion [19] was used to splice the UHD 185 hyperspectral image.After fusion, all images were stitched together using an image stitching process.

Platform
The UAV UHD 185 hyperspectral system was designed by the National Engineering Research Center for Information Technology in Agriculture (NERCITA) of China [51].The system was composed of three parts, including a ground control system, UAV and UHD 185 hyperspectral system and post-processing software.As shown in Figure 3, the ground control system included a DJI S1000 UAV (SZ DJI Technology Co., Ltd., Sham Chun, Guangdong Province, China) [52] (Figure 3), the DJI ground-control station (software, Figure 3a), reinforcement notebook computer (hardware, Figure 3a) and remote controller (Figure 3b).The UHD 185 hyperspectral system consisted of a cloud platform and the UHD 185 hyperspectral sensor (Figure 3c).The post-processing software consisted of software for stitching, as well as radiation and geometry correction (e.g., as Cubert The DJI S1000 UAV is a UAV platform with eight propellers and a multirotor, which improve stability at slow flight speeds and low altitude.When equipped with two 18,000-mAh (25 V) batteries, the UAV can run for 30 min with a take-off weight of 6 kg, flying altitude of 50 m and flight speed of 8 m/s.Flights were conducted during the flagging (26 April 2015) and flowering (13 May 2015) stages of winter wheat with ground winter wheat ASD measurement.Each flight covered a field area of 6000 square meters and lasted ~20 min.Remote sensing images acquired from high-altitude platforms typically require atmospheric correction; conversely, the UHD 185 images were obtained under stable light conditions in this study, and atmospheric correction was not required.We compared hyperspectral data acquired with the UHD 185 and ground ASD.

ASD Measurement
Ground winter wheat canopy reflectance was measured using an ASD Field Spec 3 spectrometer (Analytical Spectral Devices, Boulder, Colorado, CO, USA; abbreviated as ASD).Measurements were collected during the flagging (26 April 2015) and flowering (13 May 2015) stages of winter wheat growth.After the UAV flight (from 11 a.m.-2 p.m.), ground measurements were carried out under windless conditions and stable light levels.Winter wheat canopy reflectance was measured 10 times at the center of each plot, and the average reflectance was recorded.The reflectance obtained by ground-based ASD and UAV-mounted UHD spectrometers was compared.

Methods
After the UHD 185 images were stitched together, a ground resolution of 1 cm was obtained for the hyperspectral images.We used ENVI ROI (ENVI image processing and analysis software, from ITT Visual Information Solutions, Boulder, CO, USA) tools to calculate the average canopy spectral reflectance of each plot and selected bands of 462 nm-882 nm to avoid unstable spectra; this band range should also be selected for this study to obtain the best combination of spectral bands and indices.The first two groups of crops were used to build the AGB models, and the remaining group was used as the validation set.With two UAV flights and ground measurements, we have 64 sets of data for AGB model building and 32 sets of data for AGB model validating.

Selection of Spectral Bands and Indices
To construct linear AGB estimation models, five spectral specific bands and fourteen spectral The DJI S1000 UAV is a UAV platform with eight propellers and a multirotor, which improve stability at slow flight speeds and low altitude.When equipped with two 18,000-mAh (25 V) batteries, the UAV can run for 30 min with a take-off weight of 6 kg, flying altitude of 50 m and flight speed of 8 m/s.Flights were conducted during the flagging (26 April 2015) and flowering (13 May 2015) stages of winter wheat with ground winter wheat ASD measurement.Each flight covered a field area of 6000 square meters and lasted ~20 min.Remote sensing images acquired from high-altitude platforms typically require atmospheric correction; conversely, the UHD 185 images were obtained under stable light conditions in this study, and atmospheric correction was not required.We compared hyperspectral data acquired with the UHD 185 and ground ASD.

ASD Measurement
Ground winter wheat canopy reflectance was measured using an ASD Field Spec 3 spectrometer (Analytical Spectral Devices, Boulder, Colorado, CO, USA; abbreviated as ASD).Measurements were collected during the flagging (26 April 2015) and flowering (13 May 2015) stages of winter wheat growth.After the UAV flight (from 11 a.m.-2 p.m.), ground measurements were carried out under windless conditions and stable light levels.Winter wheat canopy reflectance was measured 10 times at the center of each plot, and the average reflectance was recorded.The reflectance obtained by ground-based ASD and UAV-mounted UHD spectrometers was compared.

Methods
After the UHD 185 images were stitched together, a ground resolution of 1 cm was obtained for the hyperspectral images.We used ENVI ROI (ENVI image processing and analysis software, from ITT Visual Information Solutions, Boulder, CO, USA) tools to calculate the average canopy spectral reflectance of each plot and selected bands of 462 nm-882 nm to avoid unstable spectra; this band range should also be selected for this study to obtain the best combination of spectral bands and indices.The first two groups of crops were used to build the AGB models, and the remaining group was used as the validation set.With two UAV flights and ground measurements, we have 64 sets of data for AGB model building and 32 sets of data for AGB model validating.

MTVI2
Preserves sensitivity to LAI and resistance to Chl influence Note: UHD 185 data did not include bands centered in 680 or 800 because its spectral resolution was 4 nm.Therefore, UHD b680 and b800 were calculated as the average UHD reflectance at 682 and 678 nm and at 802 and 798 nm (written as Band 600 and Band 800 in the following passage), respectively.

Winter Wheat Height
It was impossible to set up stable ground elevation control points in the field due to the irrigation system used at our experimental site.Therefore, we calculated the relative crop height using digital elevation model (DEM) and digital surface model (DSM).DEM and DSM are two types of surface simulations used in surveying.A DEM conveys ground surface elevation, and a DSM is the sum of the DEM and vegetation surface.Thus, crop height can be obtained by subtracting the DEM from the DSM (Figure 4 and Equation ( 2)).
Height crop = Height DSM − Height DEM In this paper, the DSM was obtained using photogrammetric technology software (Agisoft PhotoScan, Agisoft, St. Petersburg, Russia), panchromatic images and the positioning and orientation system (POS) information, all acquired by the UAV.Some soil points' coordinate information was recorded using ARCGIS software (ARCGIS, Environmental Systems Research Institute, Inc., Redlands, CA, USA), then using ARCGIS extract tools to obtain all soil points' elevation.The DEM was obtained by interpolating the surface elevation using ARCGIS kriging tools.After DEM and DSM were calculated, we use ARCGIS raster calculator tools to calculate the crop height map.After the above calculation, the soil part in the crop height map is about 0 m (black area in Figure 4), and the crop zone will show the crop height value (gray area in Figure 4).In this paper, the DSM was obtained using photogrammetric technology software (Agisoft PhotoScan, Agisoft, St. Petersburg, Russia), panchromatic images and the positioning and orientation system (POS) information, all acquired by the UAV.Some soil points' coordinate information was recorded using ARCGIS software (ARCGIS, Environmental Systems Research Institute, Inc., Redlands, CA, USA), then using ARCGIS extract tools to obtain all soil points' elevation.The DEM was obtained by interpolating the surface elevation using ARCGIS kriging tools.After DEM and DSM were calculated, we use ARCGIS raster calculator tools to calculate the crop height map.After the above calculation, the soil part in the crop height map is about 0 m (black area in Figure 4), and the crop zone will show the crop height value (gray area in Figure 4).

Partial Least Squares Regression
Partial least squares regression is a statistical method introduced by Herman Wold [63] in the 1960s as an economic technique.PLSR provides an efficient way to make full use of the hyperspectral information and was also widely used in remote sensing for vegetation studying.Previous studies [64,65] indicated that PLSR has excellent capability in making use of the full spectral information.In this paper, we employed the partial least squares regression using MATLAB (MATLAB 2014a, MathWorks, Inc., Natick, MA, USA) software on Microsoft Windows.

Precision Evaluation
Spectral data collected from two plan groups were used for model construction, and data collected from the third group were used for validation.We used the coefficient of determination (R 2 ), root mean square error (RMSE) and mean absolute error (MAE) to evaluate the performance of each model.Higher R 2 values correspond to smaller RMSE and MAE values and higher model precision.Equations (3)-( 5) were used to calculate R 2 , RMSE and MAE.

Partial Least Squares Regression
Partial least squares regression is a statistical method introduced by Herman Wold [63] in the 1960s as an economic technique.PLSR provides an efficient way to make full use of the hyperspectral information and was also widely used in remote sensing for vegetation studying.Previous studies [64,65] indicated that PLSR has excellent capability in making use of the full spectral information.In this paper, we employed the partial least squares regression using MATLAB (MATLAB 2014a, MathWorks, Inc., Natick, MA, USA) software on Microsoft Windows.

Precision Evaluation
Spectral data collected from two plan groups were used for model construction, and data collected from the third group were used for validation.We used the coefficient of determination (R 2 ), root mean square error (RMSE) and mean absolute error (MAE) to evaluate the performance of each model.Higher R 2 values correspond to smaller RMSE and MAE values and higher model precision.Equations (3)-( 5) were used to calculate R 2 , RMSE and MAE.
where x i and y i are the estimated and measured values (AGB), x and y are the average estimated and measured values and n is the sample number.

Analysis of ASD and UHD Hyperspectral Data
The canopy reflectance of different plots (Plots 2-2, 4-1, 9-1) is shown in Figure 5.For wavelengths in the range of 462-882 nm, typical vegetation spectral features are shown; strong absorption is seen at 680 nm, and reflection peaks are seen at 550 nm and 800 nm.

Analysis of ASD and UHD Hyperspectral Data
The canopy reflectance of different plots (Plots 2-2, 4-1, 9-1) is shown in Figure 5.For wavelengths in the range of 462-882 nm, typical vegetation spectral features are shown; strong absorption is seen at 680 nm, and reflection peaks are seen at 550 nm and 800 nm.Significant correlations were found between AGB and canopy spectral reflectance measured with the UHD 185 and ASD from 462 nm-882 nm during the flag leaf and flowering stages (Figure 6).Correlation results for UHD185 and ASD reflectance were similar.It shows that bands centered at 470, 550, 670 and 680 have a high correlation with AGB; conversely, Band 800 was quite low.Significant correlations were found between AGB and canopy spectral reflectance measured with the UHD 185 and ASD from 462 nm-882 nm during the flag leaf and flowering stages (Figure 6).Correlation results for UHD185 and ASD reflectance were similar.It shows that bands centered at 470, 550, 670 and 680 have a high correlation with AGB; conversely, Band 800 was quite low.
Significant correlations were found between AGB and canopy spectral reflectance measured with the UHD 185 and ASD from 462 nm-882 nm during the flag leaf and flowering stages (Figure 6).Correlation results for UHD185 and ASD reflectance were similar.It shows that bands centered at 470, 550, 670 and 680 have a high correlation with AGB; conversely, Band 800 was quite low.

Crop Height Estimation
The height of each plant was extracted using ARCGIS software; estimated and measured crop height is shown in Figure 7a.The linear relationship between estimated and measured crop height was: y = 0.977x + 0.190 (R 2 = 0.69).Crop height on 26 April and 13 May 2015 is shown in Figure 7b,c.Accuracy verification (Figure 7a) and mapping (Figure 7b,c) demonstrate that the crop height is reliable, although the estimated crop height was smaller than the measured crop height by ~0.19 m; the coefficient of determination for the relationship between the two was high.Figure 7b,c also showed the crop height changing with different growing stage, N and water treatment.

. Crop Height Estimation
The height of each plant was extracted using ARCGIS software; estimated and measured crop height is shown in Figure 7a.The linear relationship between estimated and measured crop height was: y = 0.977x + 0.190 (R 2 = 0.69).Crop height on 26 April and 13 May 2015 is shown in Figure 7b,c.Accuracy verification (Figure 7a) and mapping (Figure 7b,c) demonstrate that the crop height is reliable, although the estimated crop height was smaller than the measured crop height by ~0.19 m; the coefficient of determination for the relationship between the two was high.Figure 7b,c

Hyperspectral Model
Winter wheat AGB was significantly correlated with both specific bands and spectral vegetation indices (Table 2).Variables in Table 2 with high correlation coefficients (B470, G550, R670, R680, RVI, NDVI, GI and WDRVI) were employed to build eight single-parameter AGB estimation models (M1, M2, M3, M4, M5, M6, M7 and M8, respectively) using IBM SPSS19 software (IBM Statistical Product and Service Solutions 19.0,International Business Machines Corporation, New York, NY, USA).The model input parameters, equations and accuracy indices are shown in Table 3. Crop height also was used to build a single-parameter model (M9), as shown in Table 3.

Hyperspectral Model
Winter wheat AGB was significantly correlated with both specific bands and spectral vegetation indices (Table 2).Variables in Table 2 with high correlation coefficients (B 470 , G 550 , R 670 , R 680 , RVI, NDVI, GI and WDRVI) were employed to build eight single-parameter AGB estimation models (M1, M2, M3, M4, M5, M6, M7 and M8, respectively) using IBM SPSS19 software (IBM Statistical Product and Service Solutions 19.0,International Business Machines Corporation, New York, NY, USA).The model input parameters, equations and accuracy indices are shown in Table 3. Crop height also was used to build a single-parameter model (M9), as shown in Table 3.  Note: Modeling accuracy were calculated using 64 sets of data for AGB models building.

PLSR Methods to Improve the Estimation Accuracy
Incorporation of several spectral indices in Tables 3 and 4, B 470 , G 550 , R 670 , R 680 , RVI, NDVI, GI, WDRVI, Height, HB 470 , HG 550 , HR 670 , HR 680 , HRVI, HNDVI, HGI and HWDRVI, was combined using the PLSR method (PLSR a in Table 5) to obtain better estimation results.For comparison, the other two groups of variables (PLSR b and PLSR c, list in Table 5) were also employed to calculate AGB using the PLSR method.These PLSR models also were validated using field data.Table 5 shows the modeling precision and verification quality.Note: PLSR a, b and c represent three different combinations for PLSR method.Variables in Table 4 were employed in PLSR a; Variables in Table 3 were employed in PLSR b; Spectral variables in Table 3 were employed in PLSR c respectively.

Mapping with Improved Models Based on Crop Height
A map of AGB over the study area was constructed using the HR 680 model and UHD 185 hyperspectral imagery.An increase in AGB was observed from 26 April-13 May (Figure 9a,b).Furthermore, samples with different levels of irrigation and nitrogen fertilizer were shown to be different in AGB maps; those N0 and W0 plots had lower AGB than N2, N3, W1 and W2.However, there is no strong evidence that the AGB of W2 plots is much higher than W1 plots.The same signs appear on N3 and N2.It is quite clear that excessive fertilize and water alone could not help to increase winter wheat AGB.

Discussion
The primary aim of this study was to evaluate the performance of AGB estimation using a snapshot hyperspectral sensor mounted on a UAV platform.Canopy hyperspectral reflectance and crop height were used to construct AGB estimation models.

Stability and Redundancy of UHD Hyperspectral Data
Aasen et al. [18] showed that the difference between ASD and UHD spectral reflectance is less than 1% in the range of 466-806 nm, but the offset increases rapidly for wavelengths below 466 nm and above 850 nm.We obtained similar hyperspectral data and correlation analysis results of UHD 185, ASD hyperspectral reflectance and AGB (Figures 5 and 6).Both results showed a significant positive approximation between 462 and 800 nm, but a negative approximation between 850 and 882 nm (Figure 5).Correlations between AGB and both ASD and UHD 185 reflectance were high for wavelengths between 462 and 720 nm, but low between 750 and 882 nm (Figure 6).Based on the results, we conclude that the spectrum obtained from the UHD 185 sensor is stable in the visible and near-infrared regions (462-800 nm).These findings about UHD 185 and ASD are consistent with Cui et al. [66] and Gao et al. [65].The ASD Field Spec 3 spectrometer is designed specifically to have full-range detection capacity (from 350 nm-2500 nm) and provides uniform Vis/NIR/SWIR data collection across the entire solar irradiance spectrum, and its capacity has been publicly recognized.When measuring reflection using the ASD spectrometer, the whiteboard calibration frequency was controllable, but the UHD 185 spectrometer can take this only once before UAV take off.Additionally, the UAV flew fast, leaving the UHD 185 spectrometer an instantaneous time at one location, so it was

Discussion
The primary aim of this study was to evaluate the performance of AGB estimation using a snapshot hyperspectral sensor mounted on a UAV platform.Canopy hyperspectral reflectance and crop height were used to construct AGB estimation models.

Stability and Redundancy of UHD Hyperspectral Data
Aasen et al. [18] showed that the difference between ASD and UHD spectral reflectance is less than 1% in the range of 466-806 nm, but the offset increases rapidly for wavelengths below 466 nm and above 850 nm.We obtained similar hyperspectral data and correlation analysis results of UHD 185, ASD hyperspectral reflectance and AGB (Figures 5 and 6).Both results showed a significant positive approximation between 462 and 800 nm, but a negative approximation between 850 and 882 nm (Figure 5).Correlations between AGB and both ASD and UHD 185 reflectance were high for wavelengths between 462 and 720 nm, but low between 750 and 882 nm (Figure 6).Based on the results, we conclude that the spectrum obtained from the UHD 185 sensor is stable in the visible and near-infrared regions (462-800 nm).These findings about UHD 185 and ASD are consistent with Cui et al. [66] and Gao et al. [65].The ASD Field Spec 3 spectrometer is designed specifically to have full-range detection capacity (from 350 nm-2500 nm) and provides uniform Vis/NIR/SWIR data collection across the entire solar irradiance spectrum, and its capacity has been publicly recognized.When measuring reflection using the ASD spectrometer, the whiteboard calibration frequency was controllable, but the UHD 185 spectrometer can take this only once before UAV take off.Additionally, the UAV flew fast, leaving the UHD 185 spectrometer an instantaneous time at one location, so it was impossible to obtain enough stable radiation as the ground ASD spectrometer, which also may lead to differences in hyperspectral data.
Redundancy is an important aspect of hyperspectral data that cannot be ignored.The high number of narrow spectral bands results in a high inter-correlation (Figure 10) between them and introduces redundancy into the regression equation [64].As shown in Figure 10, the band correlation of the two spectrometers is similar.It seems that from the band of 454 nm-882 nm, both ASD and UHD 185 band correlation results can be divided into six groups (a, b, c, d, e, f) by using the correlation value.In previous studies, many methods had been employed for dimensionality reduction, such as principal component analysis, sensitivity analysis and grey relational analysis [25].In this paper, B470 (blue), G550 (green), R670 (red), R680 (red) and Nir800 (NIR) were selected based on the spectral characteristics of vegetation (Figure 5), sensitivity analysis (Figure 6) and the spectral vegetation indices (Table 1 UHD 185 band correlation results can be divided into six groups (a, b, c, d, e, f) by using the correlation value.In previous studies, many methods had been employed for dimensionality reduction, such as principal component analysis, sensitivity analysis and grey relational analysis [25].

Saturation of Hyperspectral Data and Crop Height Improvement for AGB Estimation
Of the eight single-parameter models we constructed, those based on R670 and R680 achieved the best accuracy for AGB estimations (Table 3).The estimation accuracy of VI−based models (e.g., RVI, NDVI, GI and WDRVI) was lower than that of models constructed from single specific bands, possibly because the experiments were conducted during the flagging and flowering stages.Sun et al. [67] show that the predictive ability of spectral indices decreased during the middle and late stages of crop growth, while red and near-infrared band-based spectral indices lost sensitivity.During the middle and late stages of crop growth, the wheat reproductive growth was enhanced when the vegetative growth was weakened [68,69], and the winter wheat AGB continued to increase.However, the chlorophyll content of winter wheat leaves was decreased [67], leading to the nearinfrared reflectance decreases.In our study, the low correlation results between near-infrared reflectance and biomass correspond with these previous studies.We believe that the reason for the above problem was that the correlation between near-infrared spectral reflectance and crop biomass was decreased for the change of chlorophyll.Therefore, in these wheat growing stages, for biomass estimations based on sole specific bands, spectral indices are infeasible.
AGB estimation models became saturated in the late growth stage of winter wheat using sole spectral information.Many researchers had noticed this problem a long time ago, and they usually switched to using an improved vegetation index [27] or a new type of sensor (such as LiDAR [70] or SAR [71]).Some studies using PLSR can also improve these saturated problems [43], but these models lose physical/biological meaning.Our method is to provide three-dimensional information

Saturation of Hyperspectral Data and Crop Height Improvement for AGB Estimation
Of the eight single-parameter models we constructed, those based on R 670 and R 680 achieved the best accuracy for AGB estimations (Table 3).The estimation accuracy of VI−based models (e.g., RVI, NDVI, GI and WDRVI) was lower than that of models constructed from single specific bands, possibly because the experiments were conducted during the flagging and flowering stages.Sun et al. [67] show that the predictive ability of spectral indices decreased during the middle and late stages of crop growth, while red and near-infrared band-based spectral indices lost sensitivity.During the middle and late stages of crop growth, the wheat reproductive growth was enhanced when the vegetative growth was weakened [68,69], and the winter wheat AGB continued to increase.However, the chlorophyll content of winter wheat leaves was decreased [67], leading to the near-infrared reflectance decreases.In our study, the low correlation results between near-infrared reflectance and biomass correspond with these previous studies.We believe that the reason for the above problem was that the correlation between near-infrared spectral reflectance and crop biomass was decreased for the change of chlorophyll.Therefore, in these wheat growing stages, for biomass estimations based on sole specific bands, spectral indices are infeasible.
AGB estimation models became saturated in the late growth stage of winter wheat using sole spectral information.Many researchers had noticed this problem a long time ago, and they usually switched to using an improved vegetation index [27] or a new type of sensor (such as LiDAR [70] or SAR [71]).Some studies using PLSR can also improve these saturated problems [43], but these models lose physical/biological meaning.Our method is to provide three-dimensional information combining crop height and hyperspectral information to overcome the shortcoming caused by the use of canopy hyperspectral information only.Crop height is an important parameter of winter wheat growth, and many studies have shown that biomass models based on both crop height and spectral indices have high estimation accuracy.Bending et al. [34,35] combined UAV-based spectral vegetation indices (digital camera image data) and crop height (DSM generated by digital camera image) into a biomass estimation model to improve the estimation accuracy of winter wheat AGB.Aasen et al. [18] recently suggested that snapshot hyperspectral sensors are a reliable source for crop height data.Our results support those of these studies and demonstrate that the combing of hyperspectral canopy reflectance and crop height can be used to improve the estimation accuracy of winter wheat AGB.The estimation accuracy of the AGB model based on crop height only (R 2 = 0.50, RMSE 1.62 t/ha, MAE 1.24 t/ha) was similar to that of models based on single spectral vegetation indices or specific bands.When crop height was combined with spectral vegetation indices or specific bands, the estimation accuracy greatly improved.Models constructed from crop height and HR 670 , HR 680 or HRVI achieved the best AGB estimation accuracy (Table 6).Note: 'before' represents only spectral parameters AGB estimation models, and 'after' represents crop height improved AGB estimation models.

Two-Band Spectral Vegetation Indices and PLSR Methods
In this paper, by providing three-dimensional information combining crop height and hyperspectral information, the accuracy of spectral vegetation indices (Table 4) has been greatly improved, as well as becoming briefer, simpler and having clearer physical/biological meaning.However, it is also important to recognize that two-band spectral vegetation indices (maybe three-band) make use only of small parts of the information (such NDVI) [72], and the accuracy is lower.
However, using these powerful techniques/methods was limited by the complex models, for example being difficult to express in the form of a simple formula.More importantly, the models lose the crucial physical/biological meaning.Additionally, airborne hyperspectral sensors are more expensive than narrow band multispectral sensors, and two-band spectral vegetation indices may have advantages in these applications.

Crop Height and Ground Elevation Control Points
DSM, DEM and crop height values were extracted without using ground elevation control points; therefore, crop height results may still be impacted without the ground elevation control point.Elevation in DEM was measured based on multiple bare soil points from the DSM, and an interpolation method (kriging) was used to obtain the DEM.Without the help of ground elevation control points [74], these crop height values cannot be highly accurate, and thus, our study could be improved by using stable ground elevation control points.

Conclusions
Images from the UHD 185 hyperspectral sensor combined imagery and hyperspectral information.We constructed AGB estimation models based on crop height, specific bands and spectral indices, which are extremely sensitive to winter wheat growth.Eight models based on specific bands and vegetation indices were established (Table 3), as well as eight models based on crop height combined with specific bands and vegetation indices (Table 4).The results show that the combination of hyperspectral canopy reflectance data and crop height can be used to improve the estimation accuracy of winter wheat AGB (Table 6, Figure 8).
To make full use of hyperspectral data, suitable methods have to be used.In this study, concerning the two-band spectral vegetation indices (maybe three-band) making use only of small parts of the hyperspectral information, the existence of PLSR is particularly necessary to achieve higher accuracy.
The UHD 185 sensor can improve AGB monitoring performance by providing hyperspectral and stereo imagery.Crop height obtained from this new spectrometer can greatly improve biomass estimation accuracy and mapping applications (Figure 9), which are critical for guiding agricultural management.
Experiments were conducted at the Xiao Tangshan National Precision Agriculture Research Center of China, which has an area of ~2 square kilometers and is located in Changping District (115 • 50 17"~116 • 29 49 E, 40 • 2 18 ~40 • 23 13 N) of Beijing City, China (Figure 1).Changping District has an average altitude of 36 m, a warm temperate semi-humid continental monsoon climate, average rainfall of 450 mm, average low temperature of −10 • C-7.5 • C and an average high temperature of 35 • C-40 • C. Meteorological data were acquired from the China Meteorological Data Service [50].Remote Sens. 2017, 9, 708 3 of 19

Figure 1 .
Figure 1.Location of the study area.

Figure 1 .
Figure 1.Location of the study area.

Figure 2 .
Figure 2. Experimental details.Plants were divided into three groups, with a total of 16 winter wheat plants per group.Planting area was 6 × 8 m, and each group received different fertilization and irrigation treatments.

Figure 2 .
Figure 2. Experimental details.Plants were divided into three groups, with a total of 16 winter wheat plants per group.Planting area was 6 × 8 m, and each group received different fertilization and irrigation treatments.
Remote Sens. 2017, 9, 708 7 of 19 the DEM and vegetation surface.Thus, crop height can be obtained by subtracting the DEM from the DSM (Figure4and Equation (2

Figure 4 .
Figure 4. 3D surface view of the crop height map.

Figure 4 .
Figure 4. 3D surface view of the crop height map.

Figure 5 .
Figure 5. Canopy reflectance measured by the UHD 185 and ASD from 462-882 nm during the flagging stage.Grey solid and dashed lines represent the reflectance of Plot 4-1 measured by the ASD and UHD 185, respectively; black solid and dashed lines represent the reflectance of Plot 9-1 measured by the ASD and UHD 185, respectively; red solid and dashed lines represent the reflectance of Plot 2-2 measured by the ASD and UHD 185, respectively.

Figure 5 .
Figure 5. Canopy reflectance measured by the UHD 185 and ASD from 462-882 nm during the flagging stage.Grey solid and dashed lines represent the reflectance of Plot 4-1 measured by the ASD and UHD 185, respectively; black solid and dashed lines represent the reflectance of Plot 9-1 measured by the ASD and UHD 185, respectively; red solid and dashed lines represent the reflectance of Plot 2-2 measured by the ASD and UHD 185, respectively.
also showed the crop height changing with different growing stage, N and water treatment.

Figure 8 .
Figure 8.Comparison of measured and improved model estimated AGB.Figure 8. Comparison of measured and improved model estimated AGB.

Figure 8 .
Figure 8.Comparison of measured and improved model estimated AGB.Figure 8. Comparison of measured and improved model estimated AGB.

Figure 10 .
Figure 10.The correlation coefficients among different bands of the UHD 185 and ASD spectrometer.Mostly, Groups a, b and c represent band correlation values greater than 0.9; Group d represents correlation values at 0.7−0.9; the correlation values of Groups e and f are less than 0.4.

Figure 10 .
Figure 10.The correlation coefficients among different bands of the UHD 185 and ASD spectrometer.Mostly, Groups a, b and c represent band correlation values greater than 0.9; Group d represents correlation values at 0.7−0.9; the correlation values of Groups e and f are less than 0.4.

Table 3 .
Details of AGB estimate models.

Table 4 .
Details of the improved AGB model.

Table 5 .
Modeling variables and statistics (R 2 , MAE and RMSE) between the estimated and the measured AGB. ).

Table 6 .
Statistics of the coefficient of determination (R 2 ), root mean square error (RMSE) and mean absolute error (MAE) of models in modeling and verification.