Assessment of UAV-Onboard Multispectral Sensor for Non-Destructive Site-Speciﬁc Rapeseed Crop Phenotype Variable at Di ﬀ erent Phenological Stages and Resolutions

: Unmanned aerial vehicles (UAVs) equipped with spectral sensors have become useful in the fast and non-destructive assessment of crop growth, endurance and resource dynamics. This study is intended to inspect the capabilities of UAV-onboard multispectral sensors for non-destructive phenotype variables, including leaf area index (LAI), leaf mass per area (LMA) and speciﬁc leaf area (SLA) of rapeseed oil at di ﬀ erent growth stages. In addition, the raw image data with high ground resolution (20 cm) were resampled to 30, 50 and 100 cm to determine the inﬂuence of resolution on the estimation of phenotype variables by using vegetation indices (VIs). Quadratic polynomial regression was applied to the quantitative analysis at di ﬀ erent resolutions and growth stages. The coe ﬃ cient of determination (R 2 ) and root mean square error results indicated the signiﬁcant accuracy of the LAI estimation, wherein the highest R 2 values were attained by RVI = 0.93 and MTVI2 = 0.89 at the elongation stage. The noise equivalent of sensitivity and uncertainty analyses at the di ﬀ erent growth stages accounted for the sensitivity of VIs, which revealed the optimal VIs of RVI, MTVI2 and MSAVI in the LAI estimation. LMA and SLA, which showed signiﬁcant accuracies at (R 2 = 0.85, 0.81) and (R 2 = 0.85, 0.71), were estimated on the basis of the predicted leaf dry weight and LAI at the elongation and ﬂowering stages, respectively. No signiﬁcant variations were observed in the measured regression coe ﬃ cients using di ﬀ erent resolution images. Results demonstrated the signiﬁcant potential of UAV-onboard multispectral sensor and empirical method for the non-destructive retrieval of crop canopy variables.


Introduction
Applications of the unmanned aerial vehicle (UAV) platform in precision agriculture (PA) offered a precise and reliable solution for the optimisation of crop monitoring and management [1]. UAV in PA is a new but reliable remote sensing tool that can acquire high spectral, spatial and temporal resolution data and has the advantages of low cost, flexible platform and bird's eye view for rapid data collection [2,3]. used to address the relationship between crop traits and canopy reflectance. The empirical relationship of spectral data is effective in predicting the crop biophysical traits because of its simplicity and ability to capture various canopy variations, and this approach is widely used because of its ease in computation [24].
VIs have been developed as a mathematical combination of various spectral bands of the electromagnetic spectrum and are related to various canopy parameters, which can enhance the vegetation signals by reducing the atmospheric and soil effects. The spectral VIs measure vegetation activity and demonstrate the spatial variations of the different seasons through space [7,25]. Substantial developments have been made in the interpretation and analysis of the spectral VIs, which have been applied from the field to the global level [25]. However, these indices are species-specific and are therefore not robust in different species with various leaf and canopy architectures [2,7,25].
Several studies regarding the application of UAV in rice crop yield and biomass estimation have been published [13]. The UAV application in sunflower crop has shown a significant correlation between normalized difference vegetation index (NDVI) and grain yield biomass and applied nitrogen content [2], whereas the UAV platform ('VIPtero') for onboard multispectral camera in vineyard management demonstrated a clear crop heterogeneity, which was consistent with the ground observations [9]. The multispectral sensor mounted on a UAV over potato crops demonstrated a significant correlation between NDVI and green normalized vegetation index (GNDVI) with LAI, plant cover and chlorophyll content [12].
Although the application of UAV has progressed considerably, it still has many limitations, such as payload, cost and operation consistency in agriculture. Thus, a substantial amount of work is still required. UAVs capture numerous high-resolution images that require lengthy processing times, high storage devices and high labour intensity even for professionals because of its low-altitude image acquisition capability [26,27]. Fine-resolution VIs demonstrated a significant influence on the accuracy of prediction models [28]. Although studies have determined that high-resolution images might not be the optimal choice for environmental variable and concluded that fine-resolution VIs obtain a similar prediction model accuracy [2,29,30]. The measuring scale is dependent on the association between biophysical traits and VIs from images with different resolutions, which can help select the appropriate resolution for site-specific management.
Remotely sensed VIs can offer significant information derivatives to agronomic parameters at the field scale. Therefore, the sensor that can monitor crops during the growth season at high spectral, spatial and temporal resolutions will provide useful information in the efficient and sustainable crop management with site-specific basis. Amongst the phenotypical parameters derived using remotely sensed data, LMA and SLA have received minimal attention. However, the retrieval of the mentioned parameters using the radiative transfer model (RTM) inversion could be ill-advised and computationally demanding because it requires a number of leaf and canopy variables. In this study, we optimised and validated numerous VIs for LAI and leaf dry weight (DW) prediction and investigated the potential of empirically estimated LAI and DW for the estimation of the two leaf functional traits, namely, LMA and SLA. Therefore, our study mainly aims to (1) evaluate the potential of a multispectral sensor onboard a UAV in relation to crop biophysical parameters, such as LAI and DW, during the growing season over rapeseed crops at different phenological stages, (2) assess the spectrally predicted leaf DW and LAI for the calculation of LMA and SLA at different phenological stages and (3) assess the effects of resolution on the LAI prediction over the phenological stages to provide simple, rapid and useful information for PA.

Study Area
The field experiment was conducted at the Zishi experimental station of Huazhong Agricultural University in Jingzhou (30 • 11 28.62 Province, China from 18 December 2017 to 10 March 2018. This area has a subtropical monsoon climate with an average rainfall of 1100-1300 mm, average annual temperature of 15.9 • C to 16.6 • C, hydrothermal synchronisation and consistent good climatic condition suitable for agricultural activities. Winter rapeseed was used as the experimental material in this study. To imitate the high variability of growing conditions and investigate the spectral VIs for the estimation of crop phenotype parameters, the data were collected from 16 different fields and 170 point locations. The experimental fields and sampling points are presented in Figure 1.
Remote Sens. 2020, 11, x FOR PEER REVIEW 4 of 18 hydrothermal synchronisation and consistent good climatic condition suitable for agricultural activities. Winter rapeseed was used as the experimental material in this study. To imitate the high variability of growing conditions and investigate the spectral VIs for the estimation of crop phenotype parameters, the data were collected from 16 different fields and 170 point locations. The experimental fields and sampling points are presented in Figure 1.

Biophysical Parameter Measurements
The non-destructive LAI was measured using a plant canopy analyser (SunScan, Probe type SS1, Delta-T Devices, England) over all point locations ( Figure 1) after 55 days of plantation on 18 December 2017, after 79 day of plantation on 11 January 2018, 106 days after plantation on 7 February 2018 and after 138 day of plantation on 10 March 2018 according to the crop growth stages from seedling to maturity during the growing season. Measurements were obtained on clear cloudless days. At each phenological stage, 30 subsampling points were selected randomly, and four plants were sampled destructively with their roots. The plant density was also recorded in each point, location, latitude and longitude. The samples were stored in bags and moved to the laboratory. The leaves, stems and roots of the samples were separated to obtain their fresh weights and then placed in an oven at 70 °C for 48 hours to obtain the dry mass. The fresh and dry weights of each leaf sample were determined using a high-precision digital scale. LMA and SLA were calculated as follows: where Dw and LA are the leaf DW and leaf area, respectively.

= ⁄
where LA and Dw are the leaf area and the corresponding Dw, respectively [17,31].

UAV Image Acquisition
Multispectral images were also acquired on the days of LAI measurement using the UAV platform. Prior to the flight campaign, the camera was mounted on the UAV under suitable working condition. A lightweight RedEdge MicaSense multispectral camera (MicaSense RedEdge® ) onboard a fixed-wing UAV (T-EZ; Golden Wing UAS Co., Ltd.; Chengdu, China) was used to capture the field images shown in Figure S1. A stable platform was used to adjust the camera pointing towards the nadir. The spectral bands covered the wavelength intervals 450 nm centre, 20 nm bandwidth (blue), 560 nm centre, 20 nm bandwidth (green), 668 nm centre, 10 nm bandwidth (red), 717 nm centre, 10

Biophysical Parameter Measurements
The non-destructive LAI was measured using a plant canopy analyser (SunScan, Probe type SS1, Delta-T Devices, England) over all point locations ( Figure 1) after 55 days of plantation on 18 December 2017, after 79 day of plantation on 11 January 2018, 106 days after plantation on 7 February 2018 and after 138 day of plantation on 10 March 2018 according to the crop growth stages from seedling to maturity during the growing season. Measurements were obtained on clear cloudless days. At each phenological stage, 30 subsampling points were selected randomly, and four plants were sampled destructively with their roots. The plant density was also recorded in each point, location, latitude and longitude. The samples were stored in bags and moved to the laboratory. The leaves, stems and roots of the samples were separated to obtain their fresh weights and then placed in an oven at 70 • C for 48 hours to obtain the dry mass. The fresh and dry weights of each leaf sample were determined using a high-precision digital scale. LMA and SLA were calculated as follows: where Dw and LA are the leaf DW and leaf area, respectively.
where LA and Dw are the leaf area and the corresponding Dw, respectively [17,31].

UAV Image Acquisition
Multispectral images were also acquired on the days of LAI measurement using the UAV platform. Prior to the flight campaign, the camera was mounted on the UAV under suitable working condition. A lightweight RedEdge MicaSense multispectral camera (MicaSense RedEdge®) onboard a fixed-wing Remote Sens. 2020, 12, 397 5 of 19 UAV (T-EZ; Golden Wing UAS Co., Ltd.; Chengdu, China) was used to capture the field images shown in Figure S1. A stable platform was used to adjust the camera pointing towards the nadir. The spectral bands covered the wavelength intervals 450 nm centre, 20 nm bandwidth (blue), 560 nm centre, 20 nm bandwidth (green), 668 nm centre, 10 nm bandwidth (red), 717 nm centre, 10 nm bandwidth (RedEdge), 840 nm centre and 40 nm bandwidth (near-infrared). After take-off, the UAV was programmed to follow the predefined route and complete the flight campaign. The flight altitude was 300 m and resulted in a ground resolution of 20 cm.

Image Processing
The raw multispectral images acquired by the RedEdge camera were processed in terms of image mosaicking, vignetting correction and raw digital number (DNs) for the reflectance conversion using the Pix4Dmapper (Pix4D SA, Switzerland) software. The raw digital numbers of the acquired images were firstly converted into radiance and then into surface reflectance based on the linear regression models by using the calibration target data, wherein the surface reflectance is a linear function of the digital numbers [32]. In each stage, 1024 images were captured and mosaicked from the study area, and the georeferenced images of the study area were acquired at the end of the fourth crop season. Moreover, the raw images with the ground resolution of 20 cm were resampled to 30, 50 and 100 cm by using the pixel aggregate interpolation method to assess the influence of image resolution on VIs. Lastly, the VIs were calculated.

Multispectral VIs
The canopy spectral data were used to develop the VIs, which are sensitive to canopy structure, pigments and chlorophyll, for the LAI, LMA and SLA estimation in rapeseed crops. Nine VIs were calculated, several of which have been proposed as surrogates for LAI estimation [10,25,32]. The calculated VIs, their formula and description are listed in Table 1. For all VIs, the images with the same resolution were stacked, and the mean VI values centred on the 1 m 2 location were extracted for 160 point locations in ENVI 5.1 (Exelis Visual Information Solution, Inc.; Boulder, CO, USA). Table 1. Description and formulas of the investigated VIs.

Statistical Analysis
Regression models were used for the quantitative analysis between spectral VIs as independent and LAI and leaf DW as dependent variables at the different phenological stages. A total of 140 and 30 sample points for LAI and LD were used for model calibration, respectively, whereas 30 sample points each were used for validation. To assess the model performance, the coefficient of determination (R 2 ), root mean square error (RMSE) and relative RMSE (RRMSE) were used and calculated using Equations.
(3), (4) and (5), respectively [18,42,43]. The large R 2 value and low RMSE and RRMSE values indicated the high precision and accuracy of the model as follows: where R 2 and RMSE measure the relationship of the fitting function. However, for the nonlinear functions amongst LAI and VIs, the R 2 values may be misleading because sensitivity changes continuously [10]. Thus, a different accuracy matric noise equivalent (NE) was needed to evaluate and determine the LAI estimation accuracy. NE also accounts for the slope and scattering of data points in the best-fit function as follows [10,25,44]: where RMSE is the root mean square error of the best-fit function between the spectral Vis and LAI, and d(VI)/d(LAI) is the first-order derivative of the observed relationship.

Association of VIs to the Phenolgical Stages
A wide range of variations in LAI, DW, LMA and SLA due to different canopy architectures as well as variations due to leaf sizes and shapes at different phenological stages was recorded. Table 2 presents the descriptive statistics of canopy LAI, DW, LMA and SLA of rapeseed oil at the phenological stages in the Zishi experimental station.
All the VIs investigated at the different growth stages revealed close associations with the crop reflectance characteristics. The spectral VIs demonstrated significant variations in all the phenological stages from seedling to maturity, and the results are compiled in Table 3. The VIs attained its peak values at the elongation stage and started to decline afterwards. All the VIs demonstrated a significant relation with LAI at the seedling stage. Furthermore, the MSAVI2, MTVI2 and SAVI showed the most sensitive VIs with an R 2 value of 0.87, followed by MSAVI (R 2 = 0.86), whilst NDVI and OSAVI achieved R 2 = 0.85, GNDVI R 2 = 0.82 and BNDVI obtained an R 2 value of 0.81 at the seedling stage. RVI, NDVI, MSAVI and MTVI2 showed the most influenced VIs to LAI from the seedling to maturity stages and attained the maximum R 2 (0.93-0.89, 0.88-0.86, 0.87-0.86 and 0.89-0.86) from the elongation to flowering stages. The minimum R 2 values at the maturity stage were expected due to crop senescence.

Evaluation of Spectral VIs for Estimation of Rapeseed LAI
The optimal spectral VIs for the estimation of the LAI best-fit model relationship for VIs and LAI at all the phenological stages are shown in Table 4. The spectral pattern variations and close association with crop reflectance revealed the potential capability of the explored VIs in LAI estimation. All the VIs demonstrated a high significant relationship to the temporal distribution of LAI at the different growth stages, except at the maturity stage, which revealed the lowest coefficient of determination (R 2 < 0.50) and were omitted from the calculations. Amongst the investigated phenological stages, all the VIs attained their highest coefficient of determination at the elongation stage, whereas the highest R 2 value was achieved by RVI (R 2 = 0.93), followed by MTVI2 (R 2 = 0.89) at the different phenological stages, as provided in Table 3. At the next flowering phenological stage, R 2 started to decrease until maturity, and the best-fit model of LAI with spectral VIs provided RVI (R 2 = 0.89), followed by NDVI (R 2 = 0.87) and MTVI2, MSAVI (R 2 = 0.86) and OSAVI (R 2 = 0.84), whereas SAVI attained R 2 = 0.83 and GNDVI and MSAVI2 attained R 2 = 0.82. The least R 2 value was attained by BNDVI (R 2 = 0.79). MSAVI and MTVI2 have the modifying factor to deal with the saturation problem when the LAI surpasses the saturation level (LAI > 3), whilst the sensitivity of NDVI suffered and levelled out.

Sensitivity Analysis
The sensitivity analysis results demonstrated that BNDVI and GNDVI presented the maximum insensitivity, whereas OSAVI exhibited moderately higher sensitivity and was unreliable during the entire range of LAI. NDVI and SAVI exhibited higher sensitivity for LAI < 2.5 m 2 m −2 and a decreasing sensitivity trend was observed for higher LAI > 3.0 m 2 m −2 , whereas MSAVI2 performed better than NDVI and SAVI even at LAI > 3.0 m 2 m −2 [25,45]. RVI showed higher sensitivity to LAI and is the best index for the detection of numerical variations in LAI because other indices are not robust at higher LAI values. MSAVI and MTVI2 exhibited higher and consistent sensitivity for LAI < 3.0 m 2 m −2 , and a decreasing trend was observed for higher LAI > 3.0 m 2 m −2 for MSAVI. Therefore, RVI and MTVI2 demonstrated the lowest NE values to ensure the highest sensitivities to LAI and the best multispectral indices for quantitatively detecting variations in LAI (Figure 2).

Relationship of VIs and Leaf DW
To evaluate the optimal VIs in the estimation of leaf DW (gm −2 ), the best-fit model relationship was used between the leaf DW and VIs presented in Table 1 at the different phenological stages demonstrated in Table 5. The best-fit function demonstrated the nonlinear relationship between VIs and leaf DW with R 2 that range from 0.34 to 0.68 and maximum R 2 (0.68) attained by using SAVI at the seedling stage. At the early seedling stage, SAVI performed well with the maximum R 2 value, which was demonstrated at the early growth stage, wherein the canopy is not fully developed, and the soil contributed to the reflectance. The highest R 2 (0.75) value was obtained by RVI and GNDVI, followed by NDVI with R 2 = 0.70 at the elongation stage and then began to decline afterwards. At the flowering stage, the maximum R 2 (0.72) value was achieved by MSAVI, followed by NDVI, SAVI and MSAVI2 at 0.71. All the VIs showed similar asymptotic patterns with a certain degree of scattering from the nonlinear fits.

Relationship of VIs and Leaf DW
To evaluate the optimal VIs in the estimation of leaf DW (gm −2 ), the best-fit model relationship was used between the leaf DW and VIs presented in Table 1 at the different phenological stages demonstrated in Table 5. The best-fit function demonstrated the nonlinear relationship between VIs and leaf DW with R 2 that range from 0.34 to 0.68 and maximum R 2 (0.68) attained by using SAVI at the seedling stage. At the early seedling stage, SAVI performed well with the maximum R 2 value, which was demonstrated at the early growth stage, wherein the canopy is not fully developed, and the soil contributed to the reflectance. The highest R 2 (0.75) value was obtained by RVI and GNDVI, followed by NDVI with R 2 = 0.70 at the elongation stage and then began to decline afterwards. At the flowering stage, the maximum R 2 (0.72) value was achieved by MSAVI, followed by NDVI, SAVI and MSAVI2 at 0.71. All the VIs showed similar asymptotic patterns with a certain degree of scattering from the nonlinear fits.

Estimation LMA and SLA using Spectral VIs
The performance of all the indices were compared at the different phenological stages to select the optimum index in the prediction of LAI and DW. The best-fit regression model illustrated the good relationship in the prediction of LAI and DW (Tables 4 and 5), respectively, at nearly all growth stages, except the maturity stage, which was excluded from the data. The best-fit model equations of the optimal VIs were applied for the prediction of LAI and DW at each phenological stage and used in Equations (1) and (2) in the estimation of LMA and SLA, respectively. The spatial distributions of LAI, LMA and SLA are shown in Figures S2-S4. The optimal prediction equation was used in the LAI estimation, whereas LMA and SLA were utilised in Equations (1) and (2) for the generation of maps at all the explored phenological stages, respectively. Field measurement with the estimations from composited maps provided the validation results and illustrated the potential use of LAI, LMA and SLA derivation at all the growth stages.

Validation of LAI, LMA and SLA Estimates
To validate the calibration model for LAI and estimated LMA, SLA with the predicted LAI and DW, multispectral images and ground measured data were used whilst the determined optimal VIs were considered in the prediction of LAI and DW at each growth stage and used in the estimation of LMA and SLA. The observed and estimated LAI, LMA and SLA were compared in Figures 3-6. The predictions of optimal VIs were selected on the basis of the highest R 2 and lowest RMSE, which were plotted against the ground observed LAI, LMA and SLA. The optimal VIs, RVI and MTVI2 suggested consistency between the measured and estimated LAI (R 2 > 0.80) of rapeseed oil crops at all the growth stages. By contrast, the model validation demonstrated less adequate results for LMA and SLA at the seedling stage (R 2 = 0.62, 0.25) but realised consistent results at the elongation (R 2 = 0.85, 0.80) and flowering (R 2 = 0.85, 0.71) stages ( Figure 6). Plotting the measured-estimated LMA values demonstrated less deviation from the 1:1 line, whereas SLA was strong under and over the prediction that occurs at higher SLA values, thereby suggesting the saturation of the spectral signal [46].

Evaluation of Image Resolution Effect
Four images with 20 cm pixel resolution were acquired over the rapeseed crops during the growing season and resampled to 30, 50 and 100 cm pixel resolutions. To evaluate the effects of image resolutions on the remote estimation of LAI over the rapeseed crops, the quantitative relationship

Evaluation of Image Resolution Effect
Four images with 20 cm pixel resolution were acquired over the rapeseed crops during the growing season and resampled to 30, 50 and 100 cm pixel resolutions. To evaluate the effects of image resolutions on the remote estimation of LAI over the rapeseed crops, the quantitative relationship between VIs and LAI were obtained at different image resolutions (20, 30, 50 and 100 cm). At all the mentioned image resolutions, the estimated LAI using VIs had no significant differences. The precision of the best-fit regression model between the LAI and VIs determined no significant variations and was similar at all the image resolutions. Thus, a significant high correlation between spectral VIs and LAI was acquired at the elongation and flowering stages, which are independent of image resolution. The best-fit model R 2 and RMSE are presented in Figure 7.
Remote Sens. 2020, 11, x FOR PEER REVIEW 13 of 18 between VIs and LAI were obtained at different image resolutions (20, 30, 50 and 100 cm). At all the mentioned image resolutions, the estimated LAI using VIs had no significant differences. The precision of the best-fit regression model between the LAI and VIs determined no significant variations and was similar at all the image resolutions. Thus, a significant high correlation between spectral VIs and LAI was acquired at the elongation and flowering stages, which are independent of image resolution. The best-fit model R 2 and RMSE are presented in Figure 7.

Discussion
As a remote sensing platform, UAVs are ideal for mapping and monitoring in PA [9,15]. With the recent developments in UAV technologies, numerous studies have been conducted using remote sensors onboard UAVs to evaluate its application for PA and subsequently obtained considerable

Discussion
As a remote sensing platform, UAVs are ideal for mapping and monitoring in PA [9,15]. With the recent developments in UAV technologies, numerous studies have been conducted using remote sensors onboard UAVs to evaluate its application for PA and subsequently obtained considerable attention from researchers [2,9,[12][13][14][15]25,31,47].
The key canopy biophysical and biochemical parameter retrievals can be attained adequately using certain spectral bands, that is, the green spectrum is ideal for leaf chlorophyll estimation [35]. Reflectance in the near-infrared (NIR) is recognised as sensitive to variant biomass and canopy structures, whereas the red and green reflectance values respond to the background variations and senescent vegetation [48].
This study aimed to determine the optimal VIs derived from UAV multispectral images that could be used for the remote and empirical non-destructive estimation of biophysical parameters, such as LAI, LMA and SLA over rapeseed crops at the different growth stages. The empirical method for the estimation of biophysical parameters is proficient and delivers precise estimations [24,32]. Moreover, this type of method has been widely applied for the improved quantitative accuracy of spectral VIs and used for the evaluation of crop attributes, such as LAI, N status and biomass for various crops [10,49]. In this study, the potential capability of multispectral VIs mentioned in Table 1 was evaluated over rapeseed crops at different growth stages. The results revealed that the crop variations from seedling to maturity are similar to the canopy reflectance characteristics. Variations in the crop canopy reflectance attains its peak value at the elongation stage and begins to decline afterwards until the maturity stage; this behaviour is due to crop senescence, as reported previously in [10,50]. Variations in canopy reflectance at the different spectral wave bands are synchronised with LAI because of the modification of leaf chlorophyll contents at different phenological stages [51]. Spectral reflectance from the canopy is affected not only by the biophysical characteristics of foliage but also by the direction of incidence radiation, canopy architecture and soil background [52]. The small variations in spectral reflectance at the early phenological stage in the visible region are likely due to the crop nitrogen content and soil water background. However, variations in NIR reflectance with variant leaf orientations at certain growth stages are due to the overlapping leaves that decrease the active photosynthetic size as the LAI reaches a plateau [53,54]. Table 2 demonstrates the highly significant relationship between LAI with all the aforementioned spectral VIs derived from the multispectral remote sensing data. However, a close relationship was obtained by RVI with LAI because LAI had a value of approximately 3 m 2 m −2 , followed by NDVI and MTVI2 [55]. When LAI reached its saturation point, the reduced variability of red and NIR reflectance make RVI and NDVI insensitive. However, GNDVI apparently did not reach the saturation level even at LAI with moderate to high (4)(5) values. The precision of the LAI estimate was better when green and blue bands were used instead of that that with the red band even at LAI values greater than 3 m 2 m −2 [42]. The NIR band has a considerably strong impact on the relationship between LAI and spectral reflectance and must be considered under varying crop situations [7,10,56]. The nonlinear best-fit function between LAI and VIs demonstrated inconsistent sensitivity, and R 2 and RMSE values for estimating the accuracy of LAI can be misleading. Thus, the sensitivity analysis (NE) can be applied as the precision indicator to verify the performance of VIs in LAI estimation [57,58].
LMA and SLA are important in determining plant composition and eco-physiological characterisation [17,[59][60][61][62]. Compared with the general empirical approach for the estimation of SLA, as reported in [32,48], its indirect estimation using the predicted LAI and DW illustrated highly significant retrieval accuracy. The SLA results showed remarkable variations between the different growth stages, as shown in Figure 6. The estimation results of SLA at the elongation and flowering stages are consistent with those of previously reported studies [63]. Soil background has less affected the correlation results as the data were captured at mature seedling (rosette) stage. The slight change in the performance of multispectral VIs from the UAV images could be awarded to the different viewing geometries and instrument spectral response functions [32].
Remotely sensed LMA (R 2 = 0.85, nRMSE= 0.10) and (R 2 = 0.85, nRMSE = 0.12) at the elongation and flowering stages was calculated with high estimation accuracy, as previously reported in [47]. The study estimated LMA using the LUT-based PROSAIL inversion method over the rapeseed crop. The LMA estimation results were also consistent with [19] who estimated LMA across a wide range of plant species through continuous wavelet analysis, and the study demonstrated remarkable implications in the prediction of dry matter content at the canopy level.
In terms of image resolution, our study concluded that the application of very fine-resolution remote sensing images does not reflect any significant difference. Thus, the fine-resolution images would not essentially increase the prediction accuracy of the regression model between LAI and spectral VIs. The precision of regression model between LAI and VIs did not change significantly by using images with different resolutions, as shown in Figure 7. The end users of UAV images need to consider the spatial, spectral and temporal resolutions and processing time for site-specific crop monitoring and management, as previously reported in [2,29,[64][65][66].

Conclusions
In summary, herein we reported the application of multispectral sensors onboard UAV in acquiring images of rapeseed crops and retrieving their biophysical characteristics, such as LAI, LMA and SLA, throughout the phenological stages. The results revealed that the strongest relationship of spectral VIs was obtained in the elongation stage. The sensitivity analysis between LAI and VIs revealed that RVI, MSAVI and MTVI2 were the optimal VIs for LAI estimation. LMA and SLA results demonstrated the significant estimation accuracy by using the predicted leaf DW and LAI at the elongation and flowering stages, respectively. In terms of image resolutions, robust results can still be obtained when the maximum 100 cm resolution imagery is used for oil rapeseed crop characterisation. Therefore, high-altitude images will be preferred in obtaining a decreased number of images, which will significantly influence image acquisition and processing time.
Compared with earlier remote sensing platforms, the sensor onboard UAVs had the basic advantage of reaching the targeted site and acquiring very comprehensive information throughout the growing season. Unlike satellite-and aircraft-based remote sensing platforms, the UAV remote sensing platform demonstrated better operational advantages, such as the provision of high-spatial resolution imagery with low operational cost and near-real-time data acquisition. Although satellite data were limited because of temporal resolution, cloud cover, availability at the ideal time and the operational cost for high-resolution imagery, this technology has become an important tool for large area mapping and monitoring because of its distinctive capabilities in satellite remote sensing imagery that offers extensive information in a synoptic and frequent manner.