Original Research ARTICLE
Potential of UAV-Based Active Sensing for Monitoring Rice Leaf Nitrogen Status
- 1National Engineering and Technology Center for Information Agriculture (NETCIA), Key Laboratory for Crop System Analysis and Decision Making, Ministry of Agriculture, Jiangsu Key Laboratory for Information Agriculture, Jiangsu Collaborative Innovation Center for the Technology and Application of Internet of Things, Nanjing Agricultural University, Nanjing, China
- 2Key Laboratory of Soil Environment and Pollution Remediation, Institute of Soil Science, Chinese Academy of Sciences, Nanjing, China
Unmanned aerial vehicle (UAV) based active canopy sensors can serve as a promising sensing solution for the estimation of crop nitrogen (N) status with great applicability and flexibility. This study was endeavored to determine the feasibility of UAV-based active sensing to monitor the leaf N status of rice (Oryza sativa L.) and to examine the transferability of handheld-based predictive models to UAV-based active sensing. In this 3-year multi-locational study, varied N-rates (0–405 kg N ha−1) field experiments were conducted using five rice varieties. Plant samples and sensing data were collected at critical growth stages for growth analysis and monitoring. The portable active canopy sensor RapidSCAN CS-45 with red, red edge, and near infrared wavebands was used in handheld mode and aerial mode on a gimbal under a multi-rotor UAV. The results showed the great potential of UAV-based active sensing for monitoring rice leaf N status. The vegetation index-based regression models were built and evaluated based on Akaike information criterion and independent validation to predict rice leaf dry matter, leaf area index, and leaf N accumulation. Vegetation indices composed of near-infrared and red edge bands (NDRE or RERVI) acquired at a 1.5 m aviation height had a good performance for the practical application. Future studies are needed on the proper operation mode and means for precision N management with this system.
Nitrogen (N) plays a vital role in improving crop growth and productivity (Novoa and Loomis, 1981; Ata-Ul-Karim et al., 2016). Over 200 million tons of N fertilizers are estimated to be used in 2018 and continue to increase at 1.8% per year (FAO (Food and Agriculture Organization of the United Nations), 2015). However, over-application of N fertilizers is the alarming issue that has caused low N use efficiency, leading to N deposition and water eutrophication (Conant et al., 2013; Liu et al., 2013; Huang et al., 2017). Therefore, it is imperative to develop highly efficient, reliable and practical methods for monitoring crop N status to meet the demand for precision N management (Miao et al., 2011). Several traditional methods, such as the leaf color chart (Alam et al., 2005) or destructive chemical analysis (Asner and Martin, 2008) are limited by low efficiency, small-scale applicability, and professional experience requirements for accurate diagnosis. With the advances of optical sensors and remote sensing technology, crop N-status monitoring, and management based on the spectrum has been widely used in different crops (Saberioon et al., 2014; Padilla et al., 2018).
Many studies have been conducted on the utilization of spaceborne and airborne passive remote sensing for crop monitoring. Satellite remote sensing provides great possibility for large-scale crop growth monitoring and precision management, for example, satellites images of FORMOSAT-2 were used for rice (Oryza sativa L.) N-status monitoring (Huang et al., 2015); however, the quality of remote sensing images from passive sensor-based satellites would be affected by bad weather conditions like cloud and fog, leading to the lack of applicable in-season sensing datasets for crop monitoring. The unmanned aerial vehicle (UAV) emerges as a promising remote sensing platform owing to its flexibility (Yang et al., 2017), and it was widely investigated for crop monitoring with imaging sensors (Maresma et al., 2016). Rice grain yield and leaf area index (LAI) were predicted by multi-temporal vegetation indices (VIs) from UAV-based multispectral imagery (Zhou et al., 2017), and the red edge (720 nm) and near-infrared (800 nm) band-based VIs were found to be more effective in the prediction of yield and LAI. However, despite the development of semi-automatic procedures, image processing, and analysis have still been too specialized and challenging for ordinary consumers until now.
Non-imaging optical canopy sensors directly collect standardized spectral reflectance with great flexibility in data achievement and processing over imaging sensors. Sensitive wave bands and VIs have been previously utilized for crop N-status estimation using passive hyperspectral canopy sensors (Tian et al., 2014). In addition, passive multispectral sensors were also developed for crop monitoring (Ni et al., 2016), and some of them were mounted on ground vehicle platforms (Pei et al., 2014). Ni et al. (2017) designed a UAV-mounted crop-growth monitoring system based on a passive sensor with a red band and a near-infrared band, and it was proved to have potential for predicting wheat leaf nitrogen status with Normalized Difference Vegetation Index (NDVI) and Ratio Vegetation Index (RVI).
Active sensors were developed with an internal light source to avoid the calibration requirements for illumination and the light angle (Holland et al., 2012). One of the superiorities is their potential to solve the problems of cloud cover and time limitations for measurements, which limit the use of passive sensors under such conditions (Stamatiadis et al., 2010). Most studies focus on traditional two-band active canopy sensors (Danielw and Johne, 2010; Samborski et al., 2016), but three-band active canopy sensors like Crop Circle ACS-470 (Holland Scientific Inc., Lincoln, NE, USA) were reported to improve the estimating performance of winter wheat or rice N status as compared to two-band sensors (Cao et al., 2015; Shi et al., 2015). RapidSCAN CS-45 (Holland Scientific Inc., Lincoln, NE, USA) is a small-sized portable three-band active sensor which has been used in precision agriculture. The previous study on rice indicated that VIs calculated from RapidSCAN wavebands could diagnose the rice N nutrition index well (Lu et al., 2017). In addition, studies on wheat, maize, soybean, and potato and soybean also showed the potential of the handheld RapidSCAN sensor to monitor crop N-status, to predict grain yield, and for cultivar selection, as well as for making nitrogen fertilizing recommendation (Bonfil, 2016; Aranguren et al., 2018; Miller et al., 2018).
Consequently, UAV-based active sensing is expected to offer flexibility, affordability, and applicability for large-scale monitoring compared to handheld active sensing. In addition, researchers paid much attention to combining different sensing data for establishing universal sensing approaches which are suitable for a better sensing performance or wider scale application (Gevaert et al., 2015; Schirrmann et al., 2017). With a unique character of using the same kind of sensor with a similar sensing height, UAV-based active sensing has implied a hypothesis of transferring handheld-based predictive models to UAV-based active sensing.
However, considering previous studies with manned aircraft (Lamb et al., 2009, 2014) and unmanned aerial vehicles (Krienke et al., 2017), still, little attention has been paid to investigate the possibility of applying active sensors on ultra low-altitude aerial vehicles. A potential issue for these studies was the effective sensing distance to accurately collect crop canopy reflectance. A unique feature of the RapidSCAN CS-45 is its ability to conduct height-independent spectral reflectance measurements named Pseudo Solar Reflectance (PSR) measurements introduced by its manufacturer. Based on this sensor specialty, Krienke et al. (2017) utilized the RapidSCAN CS-45 on a UAV platform for testing its applicable mode and proving its performance when evaluating maize N variability. In spite of the technical potential, sensing distance evaluation for practical application is still necessary.
Although the active canopy sensor RapidSCAN CS-45 has been proved to have applicability for diagnosing rice N-status using handheld mode (Lu et al., 2017), the potential of UAV-based active sensing for rice N-status monitoring has not yet been tested. Therefore, the objectives of the current study were two-fold: (1) to determine whether UAV-based active sensing is feasible to monitor rice leaf N-status and (2) to examine the transferability of handheld-based predictive models to UAV-based active sensing.
Materials and Methods
Study Area and Experimental Design
Field trials were carried out over three rice growing seasons (June–October 2015–2017) in Jiangsu Province of east China, which is a traditional rice-farming area with a long rice planting history (Figure 1). Trials were established at Rugao Experimental Station (32.27°N and 120.75°E, central-eastern Jiangsu) in 2015 and 2016, Sihong Experimental Station (33.37°N and 118.26°E, northern Jiangsu) in 2016, and Lianyungang Experimental Station (34.56°N and 119.32°E, northern Jiangsu) in 2017. Detailed information is presented in Table 1.
Figure 1. Study site: rice experiments conducted at Sihong, Lianyungang, and Rugao Experimental Station in Jiangsu Province of China.
Experiment 1 was conducted at Sihong Experimental Station in 2016, which covered four N rates (0, 120, 240, 360 kg N ha−1) and three rice varieties. Experiment 2 was conducted at Lianyungang Experimental Station in 2017, which covered four N rates (0, 135, 270, 405 kg N ha−1), two transplanting ways (pot-seedling and carpet-seedling mechanical transplanting), and three rice varieties. Experiment 3 was conducted at Rugao Experimental Station in 2015, which covered four N rates (0, 60, 150, 240 kg N ha−1) and one rice variety. Experiment 4 was conducted at Rugao Experimental Station in 2016, which covered four N rates (0, 100, 250, 400 kg N ha−1) and two rice varieties.
All field experiments were arranged in a randomized complete block design with three replicates. Each plot size was 56 m2 (7 × 8 m) in Experiment 1, 120 m2 (8 × 15 m) in Experiment 2, and 35 m2 (5 × 7 m) in Experiment 3 and Experiment 4. N fertilizer in all field experiments was applied in the form of granular urea as three splits: 50% before transplanting, 30% at the tillering stage, and 20% at the booting stage. For each plot, based on soil analysis and recommendations from the local agriculture department, 127 kg P2O5 ha−1 was applied before transplanting in the form of Ca(H2PO4)2 and 225 kg K2O ha−1 was applied as two splits: 50% before transplanting and 50% at the stem elongation stage. Carpet rice seedlings (for all the experiments) and pot rice seedlings (for Experiment 2) were prepared in seedling fields and transplanted into the experimental fields.
Active Canopy Sensor Data Collection
The active optical crop canopy sensor RapidSCAN CS-45 with three wavebands, including red (R, 670 nm), red-edge (Re, 730 nm), and near-infrared (NIR, over 780 nm) regions, was used in this study. Spectral reflectance (%) of each band and GPS data can be automatically collected and recorded in the memory module of the sensor at 2.5 Hz (one reading per 0.4 s) by modifying the logging method of the sensor. Data can be exported as a .csv file by PC software. Owing to the lightweight (0.8 kg) and wide measurement range (0.3 to 3 m height above the rice canopy) of the sensor, it is theoretically feasible to mount it on the UAV for practical application.
The multi-rotor UAV Spreading Wings S1000+ (DJI-Innovations Inc., Shenzhen, China) with DJI D-RTK GNSS system was used to provide a stable flight condition with accurate centimeter-level 3D positioning (Figure 2). The RapidSCAN CS-45 sensor was mounted on a customized gimbal in a fixed sensing posture under the UAV. In Experiment 2, controlled by the ground station program, the UAV aviated automatically according to the pre-concerted flight path along the central axis in the row direction of each plot with heights of 1.5 and 2 m above the canopy (resulting in 0.33 and 0.58 m2 view area), respectively. The heading speed of the UAV was set as 2 m/s.
Figure 2. Overview of the sensing equipment used in this study. (A) Spreading Wings S1000+ used as the sensing platform for low-altitude rice monitoring; (B) Flight controller module with DJI D-RTK GNSS system; (C) RapidSCAN CS-45 sensor mounted on a customized gimbal under the UAV; (D) RapidSCAN CS-45 sensor in handheld mode for data acquirement.
Handheld sensing in all experiments was conducted by the operator on approximately the same path as that of the UAV-based sensing using RapidSCAN CS-45. Sensor readings were collected ~1–1.5 m above the rice canopy (resulting in 0.15–0.33 m2 view area).
Sensor data were processed in ArcMap 10.5 (ESRI, Redlands, CA, USA) and assigned to different sampling plots with the GPS position of each reading point. A buffer of 1 m was utilized to exclude data near the plot boundary. The average reflectance values collected by RapidSCAN CS-45 were computed to represent each plot both in UAV and handheld-based sensing. Calculated spectral vegetation indices used in this study are listed in Table 2.
Plant Sampling and Measurement
Rice plant samples were acquired right after collecting sensing data. Destructive plant samples of above-ground parts were randomly collected (three hills per plot) from the sensed plants according to the average number of tillers. The samples were separated into leaves and stems. The separated leaves were scanned by Li-3000c (Li-Cor., Lincoln, NE, USA) to determine the leaf area index (LAI). Each sampled component was put into the oven for enzyme deactivation under 105°C for half an hour, and then under 80°C for 72 h for weighing. The weight of each leaf sample was used to determine the leaf dry matter (LDM). Leaf nitrogen concentration (LNC) was determined by the micro-Kjeldahl method (Bremner and Mulvaney, 1982). Leaf nitrogen accumulation (LNA) was calculated by multiplying LDM and LNC. LDM, LAI, and LNA were selected as the nitrogen indicators in this study.
As shown in Table 1, to test the possibility of building universal predictive models for different applications, the handheld datasets from Experiment 1 and Experiment 2 were used for VI-based model calibration. The handheld datasets from Experiment 3 and Experiment 4 were used to validate the regression models. For examining the transferability of the models from handheld to UAV sensors, the models derived from the handheld data were applied to the UAV data from Experiment 2. Besides, the validation on UAV data were performed for two sensing heights (1.5 and 2 m above the canopy) to evaluate the stability of the UAV system over two sensing distances.
The mean value, standard deviation (SD), and coefficient of variation (CV, %) of rice agronomic parameters were calculated using SPSS 25 (SPSS Inc., Chicago, IL, USA). For simplicity and VI comparison purpose, N indicators were first predicted with single VI-based models. Therefore, linear and three types of non-linear models (quadratic, exponential, and power) were built and evaluated between each VI and each N indicator (LDM, LAI, and LNA). The basic function forms of linear (Equation 1), quadratic (Equation 2), exponential (Equation 3), and power (Equation 4) regressions are given by:
where the response variable y represents the predicted N indicator (LDM, LAI, or LNA), and x represents the VI used for N-status estimation (NDRE, RERVI, NDVI, or RVI). The parameters a0, a1, b0, b1, b2, c0, c1, d0, and d1 were estimated using the least square method implemented in the R language environment (R Core Team).
The coefficient of determination (R2) for each model was calculated to assess the calibration performance. R2 provides a measure of how well-observed outcomes are replicated by the model, based on the proportion of total variation of outcomes explained by the model (Draper and Smith, 2014). However, R2 is not suitable to evaluate the predictive and fitting performance among models with different forms and different numbers of parameters, due to the risk of overfitting for models with higher R2 but more parameters.
Therefore, the Akaike Information Criterion (AIC) were further used for model selection in this study. The basic viewpoint of AIC is that the selected model is intended to accurately predict future data rather than to infer the “true distribution” of the calibration data (Akaike, 1974; Shmueli, 2010). Since the true predictive performance of a fitted model depends largely on the number of free parameters of the model, the AIC can be used to select model by penalizing for a large number of parameters and discouraging overfitting (Akaike, 1974; Bozdogan, 1987). The AIC is calculated using Equation 5.
where k is the number of estimated parameters in the model. In this study, including the terms of residual, k was 3 for linear, exponential, and power regressions and 4 for quadratic regression. is the maximum likelihood of the model for the data.
AICs of all the single VI-based regression models were calculated using function AIC in R language environment, which were further verified using self-complied function. The lowest value of AIC indicates the preferable model. The scatter diagrams of the selected VI-based models were plotted using GraphPad Prism 6 (GraphPad Software Inc., San Diego, CA, USA). With R2 and AIC, stepwise multiple linear regression models based on spectral reflectance of R, Re, and NIR bands were also evaluated to estimate LDM, LAI, and LNA.
Prediction for N indicators was further conducted in the independent validation datasets using the models above. The practical validation performance of the models was estimated by comparing R2, relative root mean square error (RRMSE, %) and relative error (RE, %) between the predicted variable and the true observed variable. The higher the R2 and the lower the RRMSE and RE, the higher the precision and accuracy of the model for predicting plant N indicators. The formulas of RRMSE and RE are listed as below:
where Pi and Oi are the predicted and observed value of the N indicator (LDM, LAI or LNA), respectively. is the mean of observed value of the N indicator. n is the number of samples.
Variability of Rice Leaf N-Status Indicators
Nitrogen-status indicators (LDM, LAI, and LNA) of rice varied greatly across different N rates, management practices, varieties, growth stages, sites, and years (Table 3). For the calibration dataset, the LNA exhibited the most significant variation, with a CV of 73.92%, followed by LAI and LDM with a CV of 60.84 and 59.63%, respectively. Similar results were observed for the calibration dataset and validation datasets of handheld and UAV-based sensing. In total, 624 samples which have a wide range of LDM (51.93 kg ha−1 to 5313.40 kg ha−1), LAI (0.13 to 10.34), and LNA (0.96 kg ha−1 to 192.61 kg ha−1) were involved in the calibration experiments. The large variability of these parameters was supposed to cover the major possible situation and make it conceivable to evaluate the potential of using the RapidSCAN sensor for estimating and diagnosing rice leaf N status.
Table 3. Descriptive statistics of leaf dry matter (LDM), leaf area index (LAI), and leaf nitrogen accumulation (LNA) across different growth stages, varieties, sites, and years.
Relationship Between N-Status Indicators and VIs Derived From Handheld System
The relationships between each N-status indicator and each VI were built using handheld data acquired from Experiment 1 and Experiment 2 (Table 4). The performance of using individual VIs to estimate rice LDM, LAI, and LNA varied with the form of selected VI (NDRE, RERVI, NDVI, or RVI) and model type (linear, exponential, power, or quadratic regression model) across different growth stages, treatments, sites, and years.
Table 4. R2 and AIC of the regression models between single VI (NDRE, RERVI, NDVI, or RVI) calculated from handheld sensing data (Experiment 1 and Experiment 2) and each rice N-status indicators (LDM, LAI, or LNA) across different stages of rice growth.
Based on the AIC values of the models in Table 4, the VIs calculated from NIR and Re reflectance (NDRE and RERVI) had preferable performance across different types of models as compared to VIs calculated from NIR and R reflectance (NDVI and RVI) for estimating LDM, LAI, and LNA, respectively. By comparing the models of the same type (linear, quadratic, exponential, and power models respectively), NDRE and RERVI had higher R2 values over NDVI and RVI for estimating rice leaf N indicators.
A total of 12 best models were selected for each N indicator and each VI based on AIC. Scatter diagrams of the selected models are presented in Figure 3. Overall, RERVI had great potential for estimating LDM (R2 = 0.77) and LAI (R2 = 0.79) with quadratic regression models, and the NDRE had the lowest AIC in exponential model for LNA prediction (R2 = 0.83). On the other side, obvious saturation effect and relatively poor predictive results were shown by NDVI in the scatter diagrams. For practical application, the predictive models were further evaluated in the following validation analysis.
Figure 3. The best relationships based on AIC between each rice leaf N indicator [LDM (A–D), LAI (E–H), and LNA (I–L)] and each VI by active canopy sensor RapidSCAN CS-45 across different growth stages, sites, and N treatments from the calibration experiments.
Validation of the Relationships Between N-Status Indicators and VIs Using Handheld Data
The regression models between VIs and N-status indicators selected in this study were further evaluated with other handheld sensing data acquired from validation experiments (Table 5). The scatter diagrams for the best handheld validation results of VI-based predictions (determined by the lowest RE) are shown in Figure 4.
Table 5. Validation results of the selected single VI-based models for estimating N-status indicators with handheld sensing data.
Figure 4. Validation results of LDM (A), LAI (B), and LNA (C) for single VI-based predictions using validation datasets of handheld (HH) sensing.
In the validation results between the predicted and observed N indicators, comparing NDRE with NDVI and RERVI with RVI, the VIs calculated from NIR and Re reflectance had lower RRMSE and RE than the VIs calculated from NIR and R reflectance. This showed the great potential of NDRE and RERVI for N-status estimation.
Considering the relatively slight variation of RRMSE across different VI-based models to estimate each N-status indicator, REs were compared to evaluate the models for the best handheld data-based validation (Figure 4). For LDM, RERVI performed the best using a quadratic regression model (RRMSE = 30.0%, and RE = 32.0%) in validation. While for LAI estimation, using an exponential regression model, NDRE performed better than other selected models (RRMSE = 28.6%, and RE = 32.7%). RERVI had the best validation performance by using a power regression model to estimate LNA (RRMSE = 36.9%, and RE = 29.2%).
Validation of the Models Between N-Status Indicators and VIs Using UAV-Based Data
All the selected models were also evaluated with UAV-based data from validation experiments for testing the feasibility of applying the UAV-based RapidSCAN CS-45 for monitoring rice leaf N-status with models derived from handheld data (Table 6). The UAV-based data contain two parts which were achieved from the UAV-based sensing with a 1.5 and 2 m height above the rice canopy, respectively. Scatter diagrams for UAV-based validation results of VI-based predictions with the best models (determined by the lowest RE) are shown in Figure 5.
Table 6. Validation results of the selected single VI-based models for estimating rice N-status indicators with UAV-based sensing data.
Figure 5. Validation results of LDM (A), LAI (B), and LNA (C) for single VI-based predictions using validation datasets of UAV-based sensing.
Compared with the validation performance of models from 2 m UAV-based data, VIs calculated from the same bands of 1.5 m UAV-based data better estimated LDM, LAI, and LNA. Hence, the best validation results were all from 1.5 m UAV-based datasets. This showed that the height of 1.5 m above the rice canopy is a suitable flight height for monitoring rice leaf N-status with the UAV-based RapidSCAN CS-45 sensor.
Similar to the validation results of handheld data, the results of R2, RRMSE, and RE reveal that NDRE and RERVI had better performance than NDVI and RVI from the 1.5 m height sensing dataset. Determined by the REs, the best validation result for LDM estimation was presented in Figure 5 with NDRE-based exponential model (RRMSE = 15.9%, and RE = 18.3%). While RERVI performed the best for LAI (RRMSE = 18.1%, and RE = 19.2%) and LNA (RRMSE = 20.0%, and RE = 20.8%) monitoring in the UAV-based validation result based on the quadratic model and power model, respectively.
Stepwise Multiple Linear Regression Analysis
The handheld sensing data from the calibration experiments were pooled together to build stepwise regression models with reflectance of NIR, R, and Re bands to estimate N-status indicators in the calibration experiments (Table 7). The results indicated that 74.6% of LDM variability could be explained with NIR and Re bands. The first chosen band for LAI was also the NIR band, followed by the Re band (R2 = 0.751). In contrast to the results of LDM and LAI, the R band was also chosen in the stepwise multiple linear regression models to estimate LNA after NIR and Re bands, and 77.6% of the variability of LNA could be explained. As one of the input selection rules of stepwise multiple linear regression, the AIC of the models decreased following the steps for LDM, LAI, and LNA estimation. The calibration and validation results (Table 8) based on AIC indicated that these models did not perform better than the best VI-based models, as shown above.
Table 7. Stepwise multiple linear regression models based on RapidSCAN CS-45 bands (R, Re, and NIR, %) for estimating rice N-status indicators across growth stages.
Table 8. Validation results of the stepwise multiple linear regression models based on spectral reflectance of RapidSCAN CS-45 wavebands (R, Re, and NIR) for estimating rice N-status indicators with data acquired from handheld sensing, and UAV-based sensing of a 1.5 or 2 m height above the rice canopy.
Potential of UAV-Based Active Sensing for Rice N-Status Monitoring
UAV-based active sensing, as a new exploration of the low-altitude sensing method, is expected to be outstanding for crop N-status monitoring. But for testing the feasibility of the sensing system, it is critical to consider its superiorities compared to traditional methods and its actual performance from experimental results.
Ground-based sensing using handheld or ground vehicle-mounted spectrometers is recognized to be capable of describing crop trait expression and nutrient monitoring (Saberioon et al., 2014; Yang et al., 2017); however, the data collection instability, low sensing efficiency, and high cost are its limiting factors (Yang et al., 2017). The paddy field condition has an influence on the sensing speed stability of ground-based sensing, which leads to spectral data variability and instability, and the footprints or wheel ruts may increase the influence and cause roots injured. The low sensing efficiency was shown in an example where more than 40 h were needed to collect 20,000 plots' spectral data on single rows using a single vehicle at speed about 0.56 m per second, and it would be more time-consuming by handheld sensing (about 0.5 m per second walking in paddy fields; White et al., 2012). By contrast, UAV-based sensing has advantages regarding the sensing stability (controlled with an automatic ground station program), high sensing efficiency (e.g., with a 2 m per second heading aerial speed in this study), and relatively low cost (Yang et al., 2017).
Traditional UAV-based sensing focuses on image-based utilization, and relevant studies differ a lot over the sensor use and crop species. UAV image-based sensing for crop monitoring is related to several issues like image capturing and mosaicing, geometric correction, spectral radiation processing, and useful feature extraction (Zhang and Kovacs, 2012; Yang et al., 2017). These complex processes have critical effects on the data compatibility, are time-consuming, and are closely related to the operator's experience (Zhang and Kovacs, 2012). The standard utilization process and data output form are still lacking for the practical utilization of UAV image-based sensing for crop nitrogen estimation, which limits the eurytopic model building and high throughout plant-phenotypic data analysis using diverse datasets from different studies. By collecting standard spectral reflectance directly without data processing, UAV-based active sensing has convenience regarding data collection, and it is easy-to-use for common consumers compared to image-based sensing. In addition, a low sensing height offers the UAV-based active sensing with potential ability to collect more spectral information under canopy surface compared to UAV-based passive monitoring (Holland et al., 2012; Yang et al., 2017).
For testing the actual performance of UAV-based active sensing for rice leaf N monitoring, RapidSCAN CS-45 was chosen in this study. According to the manufacturer's instructions, the sensor can distinguish its own light signal from that of surrounding ambient light by modulating the light source (by rapidly pulsing the light source on and off many times a second). Therefore, this technology named PSR measurement ensures a degree of data stability within an effective sensing distance (0.3–3 m above the canopy as introduced by the manufacturer), and this leads to the possibility of building universal predictive models for handheld and UAV-based sensing.
From the study of Krienke et al. (2017), a highly significant linear relationship was found between NDRE from handheld and UAV-based RapidSCAN sensing on a maize canopy, yet, the handheld data showed less variation as compared to UAV-based data collected with a fluctuating height from 0.5 to 1.5 m. Therefore, in this study, AIC-based optimum predictive models were built based on handheld datasets to estimate rice LDM (R2 = 0.77), LAI (R2 = 0.79), and LNA (R2 = 0.83), and it was also validated by independent handheld and UAV-based data for testing the model performance and the transferability of handheld-based sensing models to the UAV-based active sensing. The results showed great potential both for handheld and UAV-based sensing.
Evaluation of N-Status Prediction and Analysis for Saturation Effect of NDVI
An obvious saturation effect was observed with NDVI for crop N-status prediction as shown in Figure 3. NDVI saturation is a common view in previous studies. For example, NDVI achieved by the active canopy sensor GreenSeeker (Trimble Navigation Limited, Sunnyvale, CA, USA) could explain 80% of wheat biomass variability (Cao et al., 2015); however, the saturation effect of NDVI is obviously existent. The applicability of NDVI for crop nutrition monitoring is probably due to the high transmittance of the NIR band and a degree of function of the red band (Knipling, 1970). This leads to the saturation effect of NDVI (Thenkabail et al., 2000).
Therefore, using wavelengths with similar penetration into the plant canopy may be one of the methods employed to overcome the NDVI saturation problems (Niel and McVicar, 2004). The existing research has shown that the saturation problem could be reduced by the Re band, and Re-based VIs could be better correlated with crop N status (Delegido et al., 2013). Re-radiation penetrates deeper into crop canopies due to lower absorption by chlorophyll compared to radiation at the red waveband. Therefore, the sensitivity of Re reflectance is higher than R reflectance (Kanke et al., 2012; Zhou et al., 2018). Considering the four kinds of vegetation indices selected in this study for rice N-status monitoring, the VIs calculated from NIR, and Re reflectance (NDRE or RERVI) performed better than the VIs calculated from NIR and R reflectance (NDVI or RVI) both in handheld and UAV-based datasets. With similar results, Re-based vegetation indices using Crop Circle ACS-470 sensor and satellite remote sensing images improved plant concentration and uptake estimation for maize (Li et al., 2014). Moreover, NDRE and RERVI also performed more stably in UAV-based experiments in this study.
Another method consists in using ratio vegetation indices. The saturation effect is related to the normalization effect from the formula of normalized VIs, which could be avoided by the ratio VI to some degree (Gnyp et al., 2014). According to Cao et al. (2015), RVI significantly reduced the saturation effect of NDVI for estimating aboveground biomass of wheat. Based on linear models (Table 4), when comparing RVI with NDVI and RERVI with NDRE, the ratio VIs performed better for the N indicator prediction as compared to the normalized VIs (calculated by the same wavebands) in this study.
The RVI based on NIR and Re band (RERVI) was a great choice which combined two sides above. According to the potato experiment conducted by Zhou et al. (2018), a linear relationship was found between N concentration and RERVI, and this agrees with the fact that RERVI increased linearly with canopy chlorophyll content, as indicated by radiative transfer models from Clevers and Kooistra (2012). RERVI also had a good performance in calibration results compared to other VIs with linear regression models in this study.
For considering the three wavebands together, stepwise multiple linear regression analysis (SMLR) was conducted for N-status monitoring. Nevertheless, the multiple regression models did not perform significantly better than the best VI-based models. As shown by the stepwise linear regression results from the research by Cao et al. (2015), with Crop Circle ACS-470 employed on wheat across all growth stages, 53% of aboveground biomass variability and 67% of plant N uptake could be explained with two to three bands. However, the validation results also indicated that the stepwise linear regression models did not perform better than the best VI-based models. This is probably for the reason that only two or three bands were used for SMLR. If four to ten discrete wavebands could be used for SMLR, the performance of N-status monitoring would be significantly improved (Nguyen et al., 2007; Cao et al., 2013).
Besides, not only SMLR, but also some other models or methods such as support vector machines (SVM), artificial neural networks (ANN), and random forest (RF), have also been proved to have potential for monitoring crop N-status(Noh et al., 2006; Wang et al., 2013; Liang et al., 2015). Superior methods (e.g., SVM, ANN, or RF) need to be tested in future studies of combining various spectral, spatial, and environmental information for a better crop N-status prediction.
Considerations for Practical Utilization of UAV-Based Active Sensing
This study has shown that the UAV mounted active canopy sensor is feasible for monitoring rice N-status, yet some points, including sensing distance, canopy perturbance from air movement, and the slightly unstable flight condition caused by the aerodynamic ground effect of a low-altitude flight, still need to be addressed for the practical use of this sensing system.
The first issue is the suitable sensing distance which differs a lot over sensor types. For instance, spectral data of Crop Circle ACS-210 were befittingly collected on an aircraft at an altitude of 3–5 m above ground level on corn (Lamb et al., 2009). The suitable measuring height was between 0.4 and 1.2 m above the wheat canopy with a passive sensor mounted on a sensor support on a UAV (Ni et al., 2017). The sensing distance is directly reflected in the flight height above the canopy in this study, as the sensor posture has been fixed by the gimbal. Our validation results of the UAV-based dataset, which indicates that the height of 1.5 m above the rice canopy is much more suitable for rice N-status monitoring than 2 m, were in consensus with the result of distance sensitivity study on the turf grass canopy with UAV-mounted RapidSCAN CS-45 by Krienke et al. (2017). And this shows that spectral reflectance was affected by sensing distance and the UAV-based RapidSCAN sensor operated effectively within a range of 0.5–1.5 m above the canopy.
Canopy perturbance from air movement generated by the UAV is also crucial to be considered while using UAV based data acquisition. The factors influencing canopy perturbance are numerous, which include but are not limited to flight height, speed, posture, surrounding air condition and air velocities influenced by aircraft design. In the study of a UAV-based passive sensor by Ni et al. (2017), after numerical computational fluid dynamics simulations, solar sensors, and two-band sensors were designed and fixed on the two ends of a long sensor support to avoid the down-wash flow field below the UAV and the system should be applied in a hovering state. However, long sensor support is not suitable for this study considering the dynamic flight state and the sensor character. In view of the technical restriction, computational fluid dynamics simulations were not conducted in this study. While, as shown in Figure 6, via several attempts by visual checks, the perturbed canopy area (marked by the yellow box) was at the back of the sensed area when the UAV was controlled to aviate with a height of above 1.5 m and heading speed of over 2 m/s.
Figure 6. Top view and practical sensing state of the UAV-based sensing system with RapidSCAN CS-45. The height of the sensor under the UAV is 1.5 m above the rice canopy and flight speed in the heading direction is 2 m/s. The arrow symbol shows the heading flight direction of the UAV. The yellow box with increased brightness shows the perturbed canopy area generated by the UAV via a visual check.
Besides, sometimes a slightly unstable flight condition would be caused by the aerodynamic ground effect of a low-altitude flight, even though the effect is much lighter on multi-rotor UAVs than single-rotor helicopters (Sanchez-Cuevas et al., 2017). Real-Time Kinematic GPS (RTK GPS) is a powerful technology, which can provide centimeter-level high accuracy 3D positioning of UAVs (Spockeli, 2015). Therefore, to overcome the unstable condition caused by the aerodynamic ground effect, an RTK system designed for UAV is highly recommended. As shown in the UAV commissioning before the actual experiments in this study, flight condition is much more stable with the DJI D-RTK GNSS System than non-RTK utilization.
The height of 1 m above the rice canopy was also considered in the first test of UAV-based active sensing in Experiment 2. However, obvious canopy perturbance was generated under the UAV in that sensing mode. Moreover, the flight condition was unstable with a 1 m-height setting, and the flight height fluctuated from about a 0.5 to 1.2 m height above the canopy even with the D-RTK GNSS system. For data stability and experiment security, UAV-based sensing at a 1 m-height was canceled in the following tests. This unstable flight condition was probably caused by mixed reasons of the aerodynamic ground effect and the technical lack of accurate aerial positioning for proximal flight.
The calibration and validation results showed the great potential of active canopy sensor RapidSCAN CS-45 to monitor rice leaf N-status using both handheld and UAV-mounted modes. Great transferability of handheld-based predictive models to UAV-based sensing was verified by the UAV data-based validation experiment. Based on model evaluation and selection by AIC, 77, 79, and 83% of the variability in LDM, LAI, and LNA were explained with the optimal VI-based regression models derived from the calibration datasets, respectively. Considering different data acquired from UAV-based sensing and handheld sensing, NDRE and RERVI exhibited a much better performance in estimating rice N-status than the traditional R-based vegetation indices (NDVI and RVI), which also displayed great potential in overcoming the saturation problem of NDVI.
The present study has put forward a novel way of monitoring rice leaf N-status by the application of a multi-rotor unmanned aerial vehicle with a portable active canopy sensor. The height of 1.5 m above the rice canopy with a heading speed of 2 m/s was suitable for practical use. Future investigations are still needed to consider the combined effect of flight height, speed, canopy perturbance, ground effect, and new low-altitude location technology. Additionally, the entire automation workflow of data collection, processing for N status prediction, and management need to be developed for this sensing system in the future.
QC, YT, YZ, and WC conceived and designed the experiments. SL, XD and QK performed the experiments. SL and QC analyzed the data and wrote the original manuscript. SA, TC, XL, YT, YZ, and WC reviewed and revised the manuscript. All authors read and approved the final manuscript.
This research was funded by the National Key Research and Development Program (2016YFD0300608), the Fundamental Research Funds for the Central University (KYZ201502, KJQN201725), the National Natural Science Foundation of China (31601222), the Natural Science Foundation of Jiangsu Province (BK20150663), the Academic Program Development of Jiangsu Higher Education Institutions (PAPD), and the project for Student Research Training (SRT) in the College of Agriculture of Nanjing Agricultural University (1711C02).
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.
We would like to thank Jun He, Zhonghua Li, and Jiameng Ni from the College of Agriculture in Nanjing Agricultural University for their field work and contributions to data collection. We would also like to thank Jufang Wang from the College of Foreign Studies in Nanjing Agricultural University for her contributions to English corrections.
Alam, M. M., Ladha, J. K., Khan, S. R., Foyjunnessa Harun-ur-Rashid, Khan, A. H., et al. (2005). Leaf color chart for managing nitrogen fertilizer in lowland rice in Bangladesh. Agron. J. 97, 949–959. doi: 10.2134/agronj2004.0206
Aranguren, M., Castellón, A., and Aizpurua, A. (2018). Topdressing nitrogen recommendation in wheat after applying organic manures: the use of field diagnostic tools. Nutr. Cycl. Agroecosyst. 110, 89–103. doi: 10.1007/S10705-017-9865-7
Ata-Ul-Karim, S. T., Liu, X., Lu, Z., Yuan, Z., Zhu, Y., and Cao, W. (2016). In-season estimation of rice grain yield using critical nitrogen dilution curve. Field Crops Res. 195, 1–8. doi: 10.1016/j.fcr.2016.04.027
Barnes, E. M., Clarke, T. R., Richards, S. E., Colaizzi, P. D., Haberland, J., Kostrzewski, M., et al. (2000). “Coincident detection of crop water stress, nitrogen status and canopy density using ground-based multispectral data,” in Proceedings of the Fifth International Conference on Precision Agriculture, Madison, WI, 16–19 July 2000, 1–15.
Cao, Q., Miao, Y., Feng, G., Gao, X., Li, F., Liu, B., et al. (2015). Active canopy sensing of winter wheat nitrogen status: an evaluation of two sensor systems. Comput. Electron. Agricult. 112, 54–67. doi: 10.1016/j.compag.2014.08.012
Cao, Q., Miao, Y., Wang, H., Huang, S., Cheng, S., Khosla, R., et al. (2013). Non-destructive estimation of rice plant nitrogen status with Crop Circle multispectral active canopy sensor. Field Crops Res. 154, 133–144. doi: 10.1016/j.fcr.2013.08.005
Clevers, J. G. P. W., and Kooistra, L. (2012). Using hyperspectral remote sensing data for retrieving canopy chlorophyll and nitrogen content. IEEE J. Select. Top. Appl. Earth Observ. Remote Sens. 5, 574–583. doi: 10.1109/JSTARS.2011.2176468
Conant, R. T., Berdanier, A. B., and Grace, P. R. (2013). Patterns and trends in nitrogen use and nitrogen recovery efficiency in world agriculture. Global Biogeochem. Cycles 27, 558–566. doi: 10.1002/gbc.20053
Delegido, J., Verrelst, J., Meza, C. M., Rivera, J. P., Alonso, L., and Moreno, J. (2013). A red-edge spectral index for remote sensing estimation of green LAI over agroecosystems. Eur. J. Agron. 46, 42–52. doi: 10.1016/j.eja.2012.12.001
Gevaert, C. M., Suomalainen, J., Tang, J., and Kooistra, L. (2015). Generation of spectral–temporal response surfaces by combining multispectral satellite and hyperspectral UAV imagery for precision agriculture applications. IEEE J. Select. Topics Appl. Earth Observ. Remote Sens. 8, 3140–3146. doi: 10.1109/JSTARS.2015.2406339
Gnyp, M. L., Miao, Y., Yuan, F., Ustin, S. L., Yu, K., Yao, Y., et al. (2014). Hyperspectral canopy sensing of paddy rice aboveground biomass at different growth stages. Field Crops Res. 155, 42–55. doi: 10.1016/j.fcr.2013.09.023
Holland, K. H., Lamb, D. W., and Schepers, J. S. (2012). Radiometry of proximal active optical sensors (AOS) for agricultural sensing. IEEE J. Select. Topics Appl. Earth Observ. Remote Sens. 5, 1793–1802. doi: 10.1109/JSTARS.2012.2198049
Huang, J., Xu, C., Ridoutt, B. G., Wang, X., and Ren, P. (2017). Nitrogen and phosphorus losses and eutrophication potential associated with fertilizer application to cropland in China. J. Cleaner Product. 159, 171–179. doi: 10.1016/J.JCLEPRO.2017.05.008
Huang, S., Miao, Y., Zhao, G., Yuan, F., Ma, X., Tan, C., et al. (2015). Satellite remote sensing-based in-season diagnosis of rice nitrogen status in Northeast China. Remote Sens. 7, 10646–10667. doi: 10.3390/rs70810646
Jasper, J., Reusch, S., and Link, A. (2009). “Active sensing of the N status of wheat using optimized wavelength combination: impact of seed rate, variety and growth stage,” in Precision agriculture '09. Papers Presented at the European Conference on Precision Agriculture, Wageningen, 23–30.
Kanke, Y., Raun, W., Solie, J., Stone, M., and Taylor, R. (2012). Red edge as a potential index for detecting differences in plant nitrogen status in winter wheat. J. Plant Nutr. 35, 1526–1541. doi: 10.1080/01904167.2012.689912
Knipling, E. B. (1970). Physical and physiological basis for the reflectance of visible and near-infrared radiation from vegetation. Remote Sens. Environ. 1, 155–159. doi: 10.1016/S0034-4257(70)80021-9
Krienke, B., Ferguson, R. B., Schlemmer, M., Holland, K., Marx, D., and Eskridge, K. (2017). Using an unmanned aerial vehicle to evaluate nitrogen variability and height effect with an active crop canopy sensor. Precis. Agric. 18, 900–915. doi: 10.1007/s11119-017-9534-5
Lamb, D. W., Schneider, D. A., and Stanley, J. N. (2014). Combination active optical and passive thermal infrared sensor for low-level airborne crop sensing. Precis. Agric. 15, 523–531. doi: 10.1007/s11119-014-9350-0
Lamb, D. W., Trotter, M. G., and Schneider, D. A. (2009). Ultra low-level airborne (ULLA) sensing of crop canopy reflectance: a case study using a CropCircleTM sensor. Comput. Electron. Agric. 69, 86–91. doi: 10.1016/j.compag.2009.07.004
Li, F., Miao, Y., Feng, G., Yuan, F., Yue, S., Gao, X., et al. (2014). Improving estimation of summer maize nitrogen status with red edge-based spectral vegetation indices. Field Crops Res. 157, 111–123. doi: 10.1016/j.fcr.2013.12.018
Liang, L., Di, L., Zhang, L., Deng, M., Qin, Z., Zhao, S., et al. (2015). Estimation of crop LAI using hyperspectral vegetation indices and a hybrid inversion method. Remote Sens. Environ. 165, 123–134. doi: 10.1016/j.rse.2015.04.032
Lu, J., Miao, Y., Shi, W., Li, J., and Yuan, F. (2017). Evaluating different approaches to non-destructive nitrogen status diagnosis of rice using portable RapidSCAN active canopy sensor. Sci. Rep. 7:14073. doi: 10.1038/s41598-017-14597-1
Maresma, Á., Ariza, M., Martínez, E., Lloveras, J., and Martínez-Casasnovas, J. (2016). Analysis of vegetation indices to determine nitrogen application and yield prediction in maize (Zea mays L.) from a standard UAV service. Remote Sens. 8:973. doi: 10.3390/rs8120973
Miller, J. J., Schepers, J. S., Shapiro, C. A., Arneson, N. J., Eskridge, K. M., Oliveira, M. C., et al. (2018). Characterizing soybean vigor and productivity using multiple crop canopy sensor readings. Field Crops Res. 216, 22–31. doi: 10.1016/j.fcr.2017.11.006
Nguyen, H. T., Lee, K. J., Fu, J. D., Yan, Y. F., Lee, B. W., and Stafford, J. V. (2007). “Estimating rice shoot biomass and nitrogen concentration from hyperspectral canopy reflectance data using the first derivative and multiple stepwise regression analyses,” in Precision Agriculture '07. Papers Presented at the European Conference on Precision Agriculture, Skiathos, 3–6 June 2007, 643–650.
Noh, H., Zhang, Q., Shin, B., Han, S., and Feng, L. (2006). A neural network model of maize crop nitrogen stress assessment for a multi-spectral imaging sensor. Biosyst. Eng. 94, 477–485. doi: 10.1016/j.biosystemseng.2006.04.009
Padilla, F. M., Gallardo, M., Peña-Fleitas, M. T., de Souza, R., and Thompson, R. B. (2018). Proximal optical sensors for nitrogen management of vegetable crops: a review. Sensors 18:2083. doi: 10.3390/s18072083
Pei, W., Lan, Y. B., Luo, X. W., Zhou, Z. Y., Wang, Z. G., and Wang, Y. H. (2014). Integrated sensor system for monitoring rice growth conditions based on unmanned ground vehicle system. Int. J. Agric. Biol. Eng. 7, 75–81. doi: 10.3965/j.ijabe.20140702.009
Rouse, J. W., Haas, R. H., Schell, J. A., and Deering, D. W. (1973). “Monitoring vegetation systems in the Great Plains with ERTS (Earth Resources Technology Satellite),” in Proceedings of Third Earth Resources Technology Satellite Symposium, Greenbelt, ON, Canada, 10–14 December 1973, 309–317. Available at: https://ntrs.nasa.gov/search.jsp?R=19740022614 (Accessed June 6, 2018)
Saberioon, M., Amin, M. S. M., Gholizadeh, A., and Ezri, M. H. (2014). A review of optical methods for assessing nitrogen contents during rice growth. Appl. Eng. Agric. 30, 657–669. doi: 10.13031/AEA.30.10478
Samborski, S. M, Gozdowski, D., Stepien, M., Walsh, O. S, and Leszczynska, E. (2016). On-farm evaluation of an active optical sensor performance for variable nitrogen application in winter wheat. Eur. J. Agron. 74, 56–67. doi: 10.1016/j.eja.2015.11.020
Sanchez-Cuevas, P., Heredia, G., and Ollero, A. (2017). Characterization of the aerodynamic ground effect and its influence in multirotor control. Int. J. Aerospace Eng. 2017, 1–17. doi: 10.1155/2017/1823056
Schirrmann, M., Hamdorf, A., Giebel, A., Gleiniger, F., Pflanz, M., and Dammer, K.-H. (2017). Regression kriging for improving crop height models fusing ultra-sonic sensing with UAV imagery. Remote Sens. 9:665. doi: 10.3390/RS9070665
Shi, W., Lu, J., Miao, Y., Cao, Q., Shen, J., Wang, H., et al. (2015). “Evaluating a crop circle active canopy sensor-based precision nitrogen management strategy for rice in Northeast China,” in 2015 Fourth International Conference on Agro-Geoinformatics (Agro-geoinformatics), Istanbul, 20-24 July 2015, 261–264. doi: 10.1109/Agro-Geoinformatics.2015.7248112
Stamatiadis, S., Taskos, D., Tsadila, E., Christofides, C., Tsadilas, C., and Schepers, J. S. (2010). Comparison of passive and active canopy sensors for the estimation of vine biomass production. Precision Agric. 11, 306–315. doi: 10.1007/S11119-009-9131-3
Thenkabail, P. S., Smith, R. B., and Pauw, E. D. (2000). Hyperspectral vegetation indices and their relationships with agricultural characteristics. Remote Sens. Environ. 71, 158–182. doi: 10.1016/S0034-4257(99)00067-X
Tian, Y.-C., Gu, K.-J., Chu, X., Yao, X., Cao, W.-X., and Zhu, Y. (2014). Comparison of different hyperspectral vegetation indices for canopy leaf nitrogen concentration estimation in rice. Plant Soil 376, 193–209. doi: 10.1007/s11104-013-1937-0
Wang, F., Huang, J., Wang, Y., Liu, Z., and Zhang, F. (2013). Estimating nitrogen concentration in rape from hyperspectral data at canopy level using support vector machines. Precision Agric. 14, 172–183. doi: 10.1007/s11119-012-9285-2
White, J. W., Andrade-Sanchez, P., Gore, M. A., Bronson, K. F., Coffelt, T. A., Conley, M. M., et al. (2012). Field-based phenomics for plant genetics research. Field Crops Res. 133, 101–112. doi: 10.1016/j.fcr.2012.04.003
Yang, G., Liu, J., Zhao, C., Li, Z., Huang, Y., Yu, H., et al. (2017). Unmanned aerial vehicle remote sensing for field-based crop phenotyping: current status and perspectives. Front. Plant Sci. 8:1111. doi: 10.3389/fpls.2017.01111
Zhou, X., Zheng, H. B., Xu, X. Q., He, J. Y., Ge, X. K., Yao, X., et al. (2017). Predicting grain yield in rice using multi-temporal vegetation indices from UAV-based multispectral and digital imagery. ISPRS J. Photogrammetry Remote Sens. 130, 246–255. doi: 10.1016/j.isprsjprs.2017.05.003
Keywords: active canopy sensor, RapidSCAN, red edge, ultra low-level airborne, sensing distance evaluation
Citation: Li S, Ding X, Kuang Q, Ata-UI-Karim ST, Cheng T, Liu X, Tian Y, Zhu Y, Cao W and Cao Q (2018) Potential of UAV-Based Active Sensing for Monitoring Rice Leaf Nitrogen Status. Front. Plant Sci. 9:1834. doi: 10.3389/fpls.2018.01834
Received: 16 July 2018; Accepted: 27 November 2018;
Published: 14 December 2018.
Edited by:Yann Guédon, Centre de coopération internationale en recherche agronomique pour le développement (CIRAD), France
Reviewed by:Ioannis Konstantinos Christodoulakis, National and Kapodistrian University of Athens, Greece
Francisco M. Padilla, University of Almería, Spain
Copyright © 2018 Li, Ding, Kuang, Ata-UI-Karim, Cheng, Liu, Tian, Zhu, Cao and Cao. 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(s) 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.
*Correspondence: Qiang Cao, firstname.lastname@example.org