Next Article in Journal
A Smart Sensor for Defending against Clock Glitching Attacks on the I2C Protocol in Robotic Applications
Previous Article in Journal
A Formaldehyde Sensor Based on Molecularly-Imprinted Polymer on a TiO2 Nanotube Array
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Using an Active-Optical Sensor to Develop an Optimal NDVI Dynamic Model for High-Yield Rice Production (Yangtze, China)

1
National Engineering and Technology Center for Information Agriculture, Key Laboratory for Crop System Analysis and Decision Making, Ministry of Agriculture, Jiangsu Key Laboratory for Information Agriculture, Jiangsu Collaborative Innovation Center for Modern Crop Production, Nanjing Agricultural University, Nanjing 210095, China
2
Department of Agronomy and Horticulture, University of Nebraska-Lincoln, Lincoln, NE 68583, USA
*
Author to whom correspondence should be addressed.
Submission received: 7 February 2017 / Revised: 17 March 2017 / Accepted: 20 March 2017 / Published: 24 March 2017
(This article belongs to the Section State-of-the-Art Sensors Technologies)

Abstract

:
The successful development of an optimal canopy vegetation index dynamic model for obtaining higher yield can offer a technical approach for real-time and nondestructive diagnosis of rice (Oryza sativa L) growth and nitrogen (N) nutrition status. In this study, multiple rice cultivars and N treatments of experimental plots were carried out to obtain: normalized difference vegetation index (NDVI), leaf area index (LAI), above-ground dry matter (DM), and grain yield (GY) data. The quantitative relationships between NDVI and these growth indices (e.g., LAI, DM and GY) were analyzed, showing positive correlations. Using the normalized modeling method, an appropriate NDVI simulation model of rice was established based on the normalized NDVI (RNDVI) and relative accumulative growing degree days (RAGDD). The NDVI dynamic model for high-yield production in rice can be expressed by a double logistic model: RNDVI = ( 1 + e 15.2829 × ( R A G D D i 0.1944 ) ) 1 ( 1 + e 11.6517 × ( R A G D D i 1.0267 ) ) 1 (R2 = 0.8577**), which can be used to accurately predict canopy NDVI dynamic changes during the entire growth period. Considering variation among rice cultivars, we constructed two relative NDVI (RNDVI) dynamic models for Japonica and Indica rice types, with R2 reaching 0.8764** and 0.8874**, respectively. Furthermore, independent experimental data were used to validate the RNDVI dynamic models. The results showed that during the entire growth period, the accuracy (k), precision (R2), and standard deviation of RNDVI dynamic models for the Japonica and Indica cultivars were 0.9991, 1.0170; 0.9084**, 0.8030**; and 0.0232, 0.0170, respectively. These results indicated that RNDVI dynamic models could accurately reflect crop growth and predict dynamic changes in high-yield crop populations, providing a rapid approach for monitoring rice growth status.

1. Introduction

The growth and nutrient indicators of crop types have their appropriate threshold ranges. To achieve the target grain yield and quality, the shoot biomass, leaf area index (LAI), and plant nitrogen (N) concentration must be at the correct threshold ranges [1,2,3]. Furthermore, the indicator threshold during the crop’s critical growth period can be used in real-time nutrition diagnosis, and can guide the quantitative application of fertilizer.
In recent years, with increased research on crop growth models, numerous scholars have concentrated on simulating the crop group tiller number and LAI, as well as other dynamic indicators [4,5,6,7,8,9,10,11,12]. Some LAI models of major cereal crops have established, for example, the universal LAI model and optimum LAI piecewise exponential model with the growth period as a driving factor [4,5]. Other drivers included the LAI linear model and logistic curve correction model set at accumulated temperatures [6,7]. On this basis, the universal LAI model was constructed with the relative accumulated temperature as a driving factor [8]. Furthermore, a normalization method considering the accumulated temperature and LAI was adopted to establish a normalized model [9] of summer maize on the Loess Plateau. In the dynamics of dry matter (DM) accumulation, Wang et al. [10] also established a logistic correction model that could dynamically simulate the DM accumulation. However, the applicability and reliability of the model relative to different varieties and different yield levels requires improvement. Relative research showed that differences due to the varieties and production levels can be eliminated by normalizing the experimental data before modeling. Moreover, the dynamic models should be universalized according to the crop growth index [11,12], although labor- and time-consuming destructive sampling is necessary to obtain the crop LAI, DM, N content, and other growth parameters.
In the last decade, remote sensing has been widely applied in agriculture. Crop growth indicators can be rapid, real-time, non-destructive, and accurately monitored by using vegetation indices (e.g., NDVI, radio vegetation index (RVI), difference vegetation index (DVI)). Among these vegetation indices, NDVI is a key indicator used in remote sensing that utilizes the low reflectivity of the red spectrum and the high reflectivity of the near-infrared spectrum. This vegetation index can characterize the growth status, biochemical characteristics, and coverage density evaluation, and has been widely used in crop identification and yield forecasting [13]. Numerous studies have reported that the spectral index at specific periods can be used to guide crop precision fertilization and predict grain yield [14,15]. For example, Franzen and Lukina [16,17] gathered the crop canopy NDVI and constructed the N fertilization optimization algorithm (NFOA) based on crop N uptake and yield, and improved the applications on sorghum, maize, and other crops [18,19]. The spectral index at the panicle initiation stage is also used to improve subsequent panicle fertilizer application in rice (Oryza sativa L.) [20,21]. Raun et al. [22] developed the concept of the response index (RI) and potential yield based on the NDVI value gathered by GreenSeeker® (Trimble Navigation Limited, Sunnyvale, CA, USA) and the target yield of wheat in each stage. In consideration of the wheat yield at the final stage to recover the N response, Raun et al. [18] obtained an estimation of wheat production and constructed the N requirement amount algorithm. Instead of using the RI approach, Franzen et al. [23] used surface soil texture and tillage categories to refine sensor-based N recommendations. Further research was carried out to improve sensor-based yield estimation by using crop height measured via acoustic sensor [24]. Xue et al. [25,26] calculated the sufficiency index (SI) of early rice as the NDVI ratio of N deficiency treatment vs. an adequate treatment in Jiangxi Province, China. The normalization approach was used to construct a N regulation algorithm based on the SI and NDVI in rice. Canopy spectral indices better monitored crop growth to allow improved stewardship of the N topdressing during later growth stages. These spectral indices–based models for estimating crop growth status are easily affected by the crop variety and sensor type, etc. As such, the versatility and universality of these estimation models are subject to certain restrictions. If crop parameters are set as a monitoring factor, the nutrition status should be diagnosed and regulated by simulating the dynamic changes of the canopy spectral index in crop groups. Until now, few studies have been reported to establish a dynamic model based on spectral indices [27,28]. Carrying out further research on the variation in the spectral index and establishing a spectral index–based dynamic model is necessary for late crop nutrition monitoring and dynamic diagnosis.
Active canopy spectroscopy has the advantage of operating simply and non-destructively without being affected by weather conditions or the solar elevation angle. It can be used to quickly gather the canopy NDVI and other spectral indices of crops. Among the various active spectrometers, the most widely used sensor is GreenSeeker at two wavelengths (red, 656 ± 10 nm and near-infrared (NIR), 774 ± 10 nm), which is widely used to monitor crop growth and nutrient conditions. For example, Osborne [29] uses GreenSeeker to monitor wheat growth and nitrogen conditions, finding that NDVI values significantly correlated with dry matter, N content and accumulation.
The first objective of this research was to perform a comprehensive analysis of the quantitative relationship of LAI, DM and grain yield (GY) with the canopy NDVI during the key growth period in rice. The second objective was to establish an appropriate rice canopy NDVI dynamic model for high-yield production by using the data normalization method. This NDVI dynamic model for obtaining high-yield production will provide a theoretical basis for precision agriculture.

2. Materials and Methods

2.1. Experimental Details

Experiments were conducted with varied N rates (0–375 kg·N·ha–1) in five Japonica rice and three Indica rice cultivars at Nanjing Experimental Station (31°56′ N, 118°59′ E) and Rugao Experimental Station (32°14′ N, 120°19′ E) in Jiangsu Province. The soil type at Nanjing Experimental Station is yellow soil, with soil total N content of 1.26 g·kg–1, available N content of 86.16 mg·kg–1, organic matter content of 18.4 g·kg–1, available phosphorus (P) content of 12.1 mg·kg–1, and available potassium (K) content of 79.6 mg·kg–1. The soil type at Rugao Experimental Station is loam, with a total N content of 1.66 g·kg–1, available N content of 96.23 mg·kg–1, organic matter content of 12.95 g·kg–1, available P content of 13.6 mg·kg–1, and available K content of 92.6 mg·kg–1. Detailed information is shown in Table 1.
In this study, we implemented five rice field experiments using a randomized block design and three replications. Experimental cultivars, sites, and N levels were presented in Table 1. P and K were applied in the beginning as basal fertilizer (135 kg·ha–1 P2O5, 203 kg·ha–1 K2O). Among the data sets, the data from experiments 3 and 4 were used to develop the NDVI dynamic model. According to the reference data for rice planting in Jiangsu Province [30], rice yield is divided into three production levels: low-yield group (yield ≤ 8.25 t·ha–1) with N application amounts of 0 kg·ha–1 and 75 kg·ha–1; middle-yield group (8.25 t·ha–1 < yield < 10.5 t·ha–1), with N application amounts of 150 kg·ha–1 and 225 kg ha–1; high-yield group (yield ≥ 10.5 t·ha–1) with N application amounts of 300 kg·ha–1 and 375 kg·ha–1. The data from experiments 1, 2, and 5 were used to validate these models. The corresponding N levels to the high-yield group in these experiments were 330 kg·ha–1, 360 kg·ha–1, and 330 kg·ha–1, respectively.

2.2. Sample Collection and Measurement

Rice canopy spectra were measured using a handheld GreenSeeker® optical sensor (Trimble Navigation Limited, Sunnyvale, CA, USA), which measures near-infrared (780 ± 6 nm) and red light (671 ± 6 nm) bands. The NDVI and relative NDVI (RNDVI) were calculated using the measured values of the above bands. All measurements were taken on sunny days with no wind or breeze. The carried sensor probe was passed over the crop at a height of approximately 0.9 m above the crop canopy. Each cell consisted of three rows, and each row had a measurement of five replications. The analysis revealed that the spectral data varied significantly, with NDVI measurements based on the average value.
Synchronized with spectral measurements, continuous sampling was taken for five hills in each plot. Green leaf area was obtained using an LI-3000 portable area meter (Li-Cor, Lincoln, NE, USA) and expressed as LAI. Samples were oven-dried for 30 min at 105 °C to quickly cease plant metabolic activities and then at 70 °C at constant weight to attain the plant DM (t·ha1).
Grain yield was determined in each plot by harvesting plants manually from three randomly identified areas of 1 m2. Spikelets were removed from panicles and final grain yields were adjusted to 13.5% and 14.5% moisture content for Japonica and Indica rice, respectively.

2.3. Data Processing and Model Construction

The data was normalized after the maximum conversion ratio method using 1st Opt.pro V1.5 (7D-Soft High Technology Inc., Beijing, China) and Origin.pro 7.5 fitting software (Origin Lab Corporation, Northampton, MA, USA). The curve was drawn with Microsoft Excel 2013 (Microsoft, Redmond, WA, USA). The analysis of variance (ANOVA) was done using IBM SPSS 20.0 software (IBM Corporation, Armonk, NY, USA).
RAGDD i = AGDD i / AGDD h
RNDVI i = NDVI i / NDVI max
where RAGDDi is the relative accumulative growing degree days (GDD) on the i-th day after transplanting (i is the number of testing days after transplanting). AGDDi is the accumulative GDD on the i-th day after transplanting. AGDDh is the accumulated GDD from transplanting to harvest. RNDVIi is the relative NDVI on the i-th day after transplanting. NDVIi represents the NDVI measured on the i-th day after transplanting. NDVImax is the maximum NDVI of the same treatment or yield level during the entire growth period, obtained from historical or experimental data.

2.4. Model Validation

Model validation was conducted using the dataset from experiments 1, 2 and 5. The R2 and root mean square error (RMSE) were used to estimate the predictive accuracy of the model.
NDVIs ( i ) = NDVI max × RNDVIs ( i )
RMSE = 1 n × i = 1 n ( P i Q i ) 2
where NDVIs(i) is the simulated NDVI value on the i-th day after transplanting. RNDVIs(i) is the relative NDVIs(i) on the i-th day derived from the models. RMSE is the root mean square error, where n is the number of samples, Pi is the model simulated value derived from the models, and Oi is the observed value.

3. Results

3.1. Dynamic Changes of the Rice Canopy NDVI

Experiment 4 was taken as an example to analyze the dynamic characteristics of the canopy NDVI during the entire growth period (Figure 1). Rice canopy NDVI values had a similar tendency, although they varied according to the cultivars and N levels. The NDVI values increased rapidly before reaching a plateau, then slowly when near the plateau. The values remained stable after reaching the plateau before a gradual decline. The varieties of both rice types achieved maximal NDVI values at the booting stage. At the early growth stages, because of the unclosed rice canopy, the bare soil and water background affected the canopy spectral sensor, reducing the reliability of using the NDVI for the estimation of crop growth indices. From the tillering stage, the NDVI value gradually showed a steep trend along with a rapid growth curve. With a more upright structure and weaker tillering capability, the Japonica rice has lower LAI and above-ground DM values compared with Indica rice, which resulted in the lower NDVI value in Japonica rice at the same yield level. After the booting stage, NDVI values were relatively stable below the maximum value for all cultivars and N treatments. Thus, a plateau appeared, and then a gradual downward trend occurred after the heading stage.
In addition, as shown in Table 2, the Indica rice usually has a larger maximum NDVI value than the Japonica rice under the same N level, and the Indica rice tends to reach the maximum NDVI value earlier than the Japonica rice. In Japonica rice, the maximum NDVI value increased along with the increasing N rates, and the NDVI of the higher N rate treatments reached the maximum value earlier than the lower N rate treatments. Therefore, an additional N application can increase the maximum NDVI value and hasten the attainment of the maximum NDVI value. Under the same conditions, Indica rice can easily achieve a greater maximum NDVI value in a shorter time period after transplanting because of its morphological structure.

3.2. Quantitative Relationships between the Rice Canopy NDVI and Population Growth Indices

The quantitative relationships of LAI, above-ground DM and GY with the canopy NDVI were analyzed at key growth stages, namely the jointing, booting, and heading stages, based on experiments 3 and 4.

3.2.1. Quantitative Relationship between NDVI and LAI

With the advancement of the growth process, the increasing rate of the rice LAI kept pace with the canopy NDVI. The correlation between them, however, declined as the growth process moved forward (Figure 2a). This is mainly due to the gradual emergence of the panicle at the heading stage, leading to affected canopy spectra and NDVI values. With the increasing number of leaves, the shade effect becomes evident. The leaves at lower positions on the main stem and tillers cannot be well recognized by the canopy sensor, so during a fitting analysis, the NDVI value increases slowly or occurs with increasing LAI at the heading stage.

3.2.2. Quantitative Relationship between NDVI and DM

Like LAI, DM is an important indicator of growth. Across growth stages, the canopy NDVI increased as the above-ground DM increased in rice, quickly at the jointing stage but then relatively slowly (Figure 2b). The above-ground DM and canopy NDVI values showed a good correlation, with R2 ranging from 0.71 to 0.79, and explained 68.76% of the NDVI variability during the entire growth period. The relationship between the NDVI and DM was relatively stable at different growth stages.

3.2.3. Quantitative Relationship between NDVI and GY

The quantitative relationships were analyzed between the rice canopy NDVI and grain yield at key growth stages (Figure 2c). The results showed that the rice canopy NDVI was significantly positively related with the grain yield (GY) from the jointing to heading stages, with a similar relationship between the NDVI and above-ground DM (Figure 2b. The R2 values of three different stages ranged stably from 0.56 to 0.62. At the booting stage, the yield was most reliably estimated using the NDVI value.

3.3. Selection of RNDVI Dynamic Model

The AGDD, NDVI, and other agronomic indices of high-yield groups (yield ≥ 10.5 t·ha–1) were normalized in experiment 3 and experiment 4. Ten fitting equations between the RAGDD and RNDVI were established for selecting optimal models, including a double logistic curve, a rational equation, and a cubic polynomial equation. Table 3 lists the five well-fitted model types between RAGDD and RNDVI. Among them, the best determination coefficients (R2) of the double logistic, cubic polynomial, and rational equations were 0.8577**, 0.8357**, and 0.8319**, respectively.
Figure 3 showed the results of the further analysis of the above three models. A good correlation of the three fitting equations can simulate the dynamic changes of the RNDVI for the high-yield rice population. However, for the polynomial equation, when x tends toward ∞, the simulated value of the RNDVI does not match the actual variation. The rational equation, which falls obviously from the booting to heading stages, cannot express the high plateau period of the NDVI curve well in rice during the same growth period. In contrast, the double logistic equation can simulate the RNDVI with the dynamics of RAGDD (Equation (5)).
R N D V I = ( 1 + e a × ( R A G D D i b ) ) 1 ( 1 + e c × ( R A G D D i d ) ) 1
where RAGDDi indicates the relative AGDD value on the i-th day after transplanting, and RNDVI corresponds to the relative NDVI value during the same period. Here, a and c are the crop’s two inflection points in the growth and senescence logistic curves, while b and d are the two time points, expressed as the RAGDD corresponding to the critical conversion time of a and c, respectively.

3.4. Establishment of the RNDVI Dynamic Models

The diversity of plant morphology exists due to different rice cultivars. To further improve the forecast accuracy of the RNDVI dynamic model, models of Japonica and Indica rice under different production levels were constructed using the double logistic method (see Figure 4). With the relevant parameters shown in Table 4, the coefficients of determination (R2) of the NDVI dynamics in the two rice types were all above 0.86. Furthermore, a significant test was performed between these model parameters (Table 4) and the double logistic model’s parameters (Table 3). The results demonstrated that all values were non-significant at the 5% level (i.e., the resulting t-values were >0.05). Thus, the relative dynamic model of the NDVI can be used to accurately simulate the NDVI dynamics of rice high-yield groups.
The RNDVI dynamic models of Indica (RNDVIIn) and Japonica (RNDVIJa) rice are as follows:
R N D V I I n = ( 1 + e 23.8261 × ( R A G D D i 0.1489 ) ) 1 ( 1 + e 12.0923 × ( R A G D D i 1.0361 ) ) 1
R N D V I J a = ( 1 + e 20.0313 × ( R A G D D i 0.2370 ) ) 1 ( 1 + e 10.9741 × ( R A G D D i 1.0195 ) ) 1
RNDVI dynamic models constructed at different yield levels for both Indica and Japonica rice showed a similar tendency (Figure 4). However, some differences still existed in the simulated dynamics of the NDVI. The greater the yield potential was, the higher the increasing rate of the NDVI before reaching the peak value, the longer the NDVI peak plateau, and the lower the decreasing rate of the NDVI after the peak value. This resulted in the higher value of parameter a, the lower value of parameter b and c, and the higher parameter d.
Figure 5 shows the differences in RNDVI changes for different rice types. A large difference between the two rice types occurred at the early growth stage; in particular, the growth rate of the RNDVI in Indica rice was faster than that in Japonica rice, and the Indica rice reached the double logistic curve inflection point earlier than the Japonica rice did, mainly due to the Indica’s plant structure. At the same RAGDD, Indica rice achieved higher RNDVI values because of the loose leaf morphology, and higher LAI and DM. The differences in the maximum RNDVI value between the two types of rice were eliminated by using the data normalization method, which is reflected in the growth rate of the RNDVI. After a stable plateau, rice populations entered the aging period, and the aging trend of the two types of rice was quite similar. The time inflection point (d value), the RNDVI decline rate (c value), and the maximum rate of decrease were not significantly different.

3.5. Model Validation

Independent observed NDVI values were used to validate the NDVI dynamic models for Indica and Japonica rice under the high-yield level. Simulated NDVI values from 15 July (tillering stage) to 6 September (heading stage) under the high-yield level were calculated by the models with RAGDD as a driving factor. The validation result is shown in Table 5, where the k value was between 0.9692 and 1.0330, and all were around 1, and R2 (representing the simulation accuracy, as below) was between 0.6102 and 0.9331, reaching a highly significant level. The RMSE values between the observed and simulated NDVI values were between 0.0064 and 0.0305, indicating high simulation accuracy. From the accuracy and precision of the simulation results, the pre- and post-simulation results better reflected the changes in the population dynamics. The comprehensive comparisons of simulated and observed values over entire growth stages revealed that the k values of Japonica and Indica rice were 0.9991 and 1.0170, respectively; the R2 values were 0.9084 and 0.8030; and the RMSE were 0.0232 and 0.0170 (Figure 6). As the sample number increased, the fitting effect throughout the growth period also reached a high degree of accuracy. Results indicated the RNDVI dynamic models could accurately reflect and simulate the population growth dynamics of high-yield rice.

4. Discussion

Population growth indices such as LAI and DM can effectively designate crop growth conditions [30,31]. Remote sensing provides a new approach to monitor growth indices, nutrition indices and yields of crops [32,33]. GreenSeeker active-optical sensors effectively assess crop growth conditions and facilitate post-assessment of the N requirement. This also has the advantage of operating simply and non-destructively, rendering it superior to other gathering methods or indicators [34,35,36,37]. Different crop cultivars and N levels can lead to different LAI values [38]; Indica rice normally has higher LAI values than Japonica rice due to its loose leaf morphology. This study found that rice canopy NDVI values had a similar tendency during the entire growth period, and both achieved maximal NDVI values at the booting stage. In addition, researchers conducted several experiments on maize, rice and wheat, and found that the spectral index can be used to monitor crops’ main growth indices [39,40,41,42,43,44]. This confirmed, to some extent, the stable relationships of the canopy NDVI with LAI, above-ground DM and GY performed at key growth stages (Figure 2). Furthermore, we found that the NDVI became saturated when the value was >0.7, a slight variation on this value compared to Goswami’s result [45]. This is due to different crop varieties and eco-sites.
When previous dynamic models based on crop growth indicators and spectral indices are applied for crop production, destructive sampling may be needed to obtain the LAI and other agronomic parameters. The simulation could proceed based on growing degree days, or a sub-function simulation may be required, which restricts the applications of these models [4,5,11,45] to some extent. In this study, the RNDVI models based on RAGDD were constructed for both Japonica and Indica high-yield rice types. The accuracy of the model reached 0.86, and the accuracy of the RNDVI model fitted separately on Indica and Japonica rice also reached 0.89 and 0.88, respectively. The normalization of the data can weaken the effects of cultivars and other factors on the model, allowing common trends in crop production. By taking crop growth into consideration, and fully analyzing the curve characteristics and range of each model, only the double logistic curve could properly express and simulate the canopy RNDVI dynamic characteristics for high-yield rice groups. This includes rapid growth after tillering, the stable plateau during booting, and a slow decline after heading. Furthermore, accumulative GDD was selected as a model variable parameter, and the AGDD for the similar crop growth period had a determined threshold range. These parameters promoted the use of the RNDVI dynamic model in production and practice. Therefore, the present study compensates for the weakness of previous models with a complex structure or more input parameters requiring destructive sampling. This study used an active-optical sensor to explore a rapid and convenient method for acquiring the dynamic changes of crop NDVI. The model has the advantages of a simple structure, fewer input parameters, and ease of use.
In addition, experiments using N levels were only used on five rice cultivars in the present study. Moreover, experiment sites were also limited to the Yangtze River area, where cultivar types and eco-sites are relatively homogeneous. These factors may affect model parameters, reducing the universality of the model [11]. Further study is required to solve these problems. At the same time, in order to deeply explain the NDVI dynamic characteristics under the crop high-yield level and to regulate top-dressing fertilizer, the relationships among the model parameters, biomass production, LAI, N accumulation, unit tiller number, and other growth indicators at different levels are still required for accurate quantification. Furthermore, more experimental data is required to evaluate and calibrate the model in the future.

5. Conclusions

The canopy NDVI has a relatively stable and positive relationship with rice growth indices such as the LAI, above-ground dry matter and grain yield, etc. By using the normalizing method, the canopy RNDVI dynamic model based on accumulative growing degree days was constructed for high-yield production of rice in the Yangtze River region. Furthermore, the variation among different cultivars was taken into consideration, requiring the construction of relative NDVI (RNDVI) dynamic models for Japonica and Indica rice types, with R2 reaching 0.8764** and 0.8874**, respectively. These results indicate that RNDVI dynamic models can accurately reflect crop growth and predict dynamic changes in high-yield crop populations, providing a rapid and non-destructive approach for diagnosing rice growth status and formulating N nutrition diagnosis. This paradigm can be used to model and predict the canopy NDVI for obtaining high-yield production in rice, which also provides a technical pattern for real-time diagnosis of crop growth and nutrition conditions.

Acknowledgments

The work was supported by the National Key Research & Development Program of China (Nos.: 2016YFD0300604; 2016YFD0200602), the Special Program for Agriculture Science and Technology from the Ministry of Agriculture in China (No.: 201303109), the Fundamental Research Funds for the Central Universities (No.: 262201602), the Priority Academic Program Development of Jiangsu Higher Education Institutions of China (PAPD), the 111 project (B16026) and Science and Technology Support Plan of Jiangsu Province (BE2015310).

Author Contributions

Xiaojun Liu and Yan Zhu conceived and designed the experiments, Hengbiao Zheng and Qiang Cao performed experiments with support by Yongchao Tian, Xiaojun Liu and Hengbiao Zheng analyzed the data, Xiaojun Liu and Yan Zhu wrote the paper, Weixing Cao and Richard B. Ferguson provided advice and edited the manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Ata-Ul-Karim, S.T.; Yao, X.; Liu, X.J.; Cao, W.X.; Zhu, Y. Development of critical nitrogen dilution curve of Japonica rice in Yangtze River Reaches. Field Crops Res. 2013, 149, 149–158. [Google Scholar] [CrossRef]
  2. Ling, Q.H. Crop Population Quality; Shanghai Scientific & Technical Publisher: Shanghai, China, 2000; pp. 40–45. [Google Scholar]
  3. Ling, Q.H. Rice Precise Quantitative Cultivation Theory and Technology; China Agriculture Press: Beijing, China, 2007; pp. 123–128. [Google Scholar]
  4. Yu, Q.; Fu, B.P.; Yao, K.M. A universal growth model on rice LAI. Agric. Meteorol. 1995, 16, 6–8. [Google Scholar]
  5. Cao, H.X.; Dong, Y.H.; Wang, X.Q.; Xu, J.F.; Gao, L.Z. Studies of dynamic simulation models of optimum leaf area index of wheat under different yielding levels. J. Triticeae Crops. 2006, 26, 128–131. [Google Scholar]
  6. Liu, Z.D.; Duan, A.W.; Gao, Y.; Liu, H. Study on dynamic model of leaf area index (LAI) for winter wheat in Xinxiang area. J. Triticeae Crops. 2008, 4, 680–685. [Google Scholar]
  7. Wang, N.Z.; Li, Z.; Huang, B.; Ma, S.J.; Chen, Y.H. The impact of precipitation to leaf area and dry matter accumulate and yield of maize on Loess Plateau. J. Maize Sci. 2007, 1, 100–102. [Google Scholar]
  8. Lin, Z.H.; Xiang, Y.Q.; Mo, X.G.; Li, J.; Wang, L. Normalized leaf area index model for summer maize. Chin. J. Eco-Agric. 2003, 11, 69–72. [Google Scholar]
  9. Zhang, X.D.; Cai, H.J.; Fu, Y.J.; Wang, J. Study on leaf area index of summer maize in loess areas. Agric. Res. Arid Areas. 2006, 24, 25–29. [Google Scholar]
  10. Wang, X.L. Logistic simulating model on crop dry matter accumulating dynamics. Agric. Meteorol. 1986, 7, 14–19. [Google Scholar]
  11. Zhang, B.; Zhao, M.; Dong, Z.Q.; Li, J.G.; Chen, C.Y.; Sun, R. Establishment and test of LAI dynamic simulation model for high yield population. Acta Agron. Sin. 2007, 33, 612–619. [Google Scholar]
  12. Hou, Y.H.; Chen, C.Y.; Guo, Z.Q.; Hou, L.B.; Dong, Z.Q.; Zhao, M. Establishment of dry matter accumulation dymamic simulation model and analysis of growth charateristc for high-yielding population of spring maize. J. Maize Sci. 2008, 16, 90–95. [Google Scholar]
  13. Chen, J.F.; Chen, K.M.; Xu, J.B. Research on the remote sensing monitoring of grassland productivity based on TM-NDVI. Agric. Sci. Technol. 2011, 12, 119–122. [Google Scholar]
  14. Sharma, L.K.; Bu, H.G.; Denton, A.; Franzen, D.W. Active-Optical Sensors Using Red NDVI Compared to Red Edge NDVI for Prediction of Corn Grain Yield in North Dakota, U.S.A. Sensors 2015, 15, 27832–27853. [Google Scholar] [CrossRef] [PubMed]
  15. Liu, N.S.; Cao, W.X.; Zhu, Y.; Zhang, J.C.; Pang, F.R.; Ni, J. The Node Deployment of Intelligent Sensor Networks Based on the Spatial Difference of Farmland Soil. Sensors 2015, 15, 28314–28339. [Google Scholar] [CrossRef] [PubMed]
  16. Lukina, E.V.; Freeman, K.W.; Wynn, K.J.; Thomason, W.E.; Mullen, R.W.; Stone, M.L.; Solie, J.B.; Klatt, A.R.; Johnson, G.V.; Elliott, R.L.; et al. Nitrogen fertilization optimization algorithm based on in-season estimates of yield and plant nitrogen uptake. J. Plant Nutr. 2007, 24, 885–898. [Google Scholar] [CrossRef]
  17. Franzen, D.W.; Sharma, L.K.; Bu, H. Active Optical Sensor Algorithms for Corn Yield Prediction and a Corn Side-Dress Nitrogen Rate Aid; North Dakota State University Extension Service: Fargo, ND, USA, 2014. [Google Scholar]
  18. Raun, W.R.; Solie, J.B.; Johnson, G.V.; Stone, M.L.; Mullen, R.W.; Freeman, K.W.; Thomason, W.E.; Lukina, E.V. Improving nitrogen use efficiency in cereal grain production with optical sensing and variable rate application. Agron J. 2002, 94, 815–820. [Google Scholar] [CrossRef]
  19. Raun, W.R.; Solie, J.B.; Stone, M.L.; Martin, K.L.; Freeman, K.W.; Mullen, R.W.; Zhang, H.N.; Schepers, J.S.; Johnson, G.V. Optical sensor-based algorithm for crop nitrogen fertilization. Soil Sci. Plant Anal. 2005, 36, 2759–2781. [Google Scholar] [CrossRef]
  20. Nguyen, H.T.; Lee, K.; Lee, B.W. Recommendation of nitrogen topdressing rates at panicle initiation stage of rice using canopy reflectance. J. Crop Sci. Biotechnol. 2008, 11, 141–150. [Google Scholar]
  21. Tubana, B.S.; Harrell, D.; Walker, T.; Phillips, S. Midseason nitrogen fertilization rate decision tool for rice using remote sensing technology. Better Crops. 2011, 95, 22–24. [Google Scholar]
  22. Raun, W.R.; Solie, J.B.; Johnson, G.V.; Stone, M.L.; Lukina, E.V.; Thomason, W.E.; Chepers, J.S. In-season prediction of potential grain yield in winter wheat using canopy reflectance. Agron. J. 2001, 93, 131–138. [Google Scholar] [CrossRef]
  23. Franzen, D.W.; Kitchen, N.; Holland, K.; Schepers, J.; Raun, W. Algorithms for in-season nutrient management in cereals. Agron J. 2016, 108, 1775–1781. [Google Scholar] [CrossRef]
  24. Sharma, L.K.; Bu, H.; Franzen, D.W.; Denton, A. Use of corn height measured with an acoustic sensor improves yield estimation with ground based active optical sensors. Comput. Electron. Agric. 2016, 124, 254–262. [Google Scholar] [CrossRef]
  25. Xue, L.H.; Li, G.H.; Qin, X.; Yang, L.Z.; Zhang, H.N. Topdressing nitrogen recommendation for early rice with an active sensor in south China. Precis. Agric. 2014, 15, 95–110. [Google Scholar] [CrossRef]
  26. Xue, L.H.; Yang, L.Z. Recommendations for nitrogen fertilizer topdressing rates in rice using canopy reflectance spectra. Biosyst. Eng. 2008, 100, 524–534. [Google Scholar] [CrossRef]
  27. Yao, X.F. Non-Destructive Monitoring and Precise Regulation of Plant Nitrogen Nutrition in Wheat; Nanjing Agricultural University: Nanjing, China, 2012. [Google Scholar]
  28. Zheng, H.B.; Cheng, T.; Yao, X.; Deng, X.Q.; Tian, Y.C.; Cao, W.X.; Zhu, Y. Detection of rice phenology through time series analysis of ground-based spectral index data. Field Crops Res. 2016, 198, 131–139. [Google Scholar] [CrossRef]
  29. Osborne, S.L. Utilization of Existing Technology to Evaluate Spring Wheat Growth and Nitrogen Nutrition in South Dakota. Commun. Soil Sci. Plant Anal. 2007, 38, 949–958. [Google Scholar] [CrossRef]
  30. Du, Y.L.; Wang, Q.S.; Wang, C.L.; Ding, Y.F. Yield potential of rice and technical approaches to high yield in Jiangsu province. Jiangsu J. Agric. Sci. 2011, 27, 926–932. [Google Scholar]
  31. Ata-Ul-Karim, S.T.; Zhu, Y.; Yao, X.; Cao, W.X. Determination of critical nitrogen dilution curve based on leaf area index in rice. Field Crops Res. 2014, 167, 76–85. [Google Scholar] [CrossRef]
  32. Yao, X.; Ata-Ul-Karim, S.T.; Zhu, Y.; Tian, Y.C.; Liu, X.J.; Cao, W.X. Development of critical nitrogen dilution curve in rice based on leaf dry matter. Eur. J. Agron. 2014, 55, 20–28. [Google Scholar] [CrossRef]
  33. Clevers, J.G.P.W.; Büker, C.; Leeuwen, H.J.C.V.; Bouman, B.A.M. A framework for monitoring crop growth by combining directional and spectral remote sensing information. Remote Sens. 1994, 50, 161–170. [Google Scholar] [CrossRef]
  34. Wang, R.C.; Wang, K.; Shen, Z.Q.; Jiang, H.X.; Zhu, D.F.; Cai, T.C. A study on modeling of rice yield estimation by remote sensing. J. Remote Sens. 1998, 2, 19–124. [Google Scholar]
  35. Lu, Y.L.; Bai, Y.L.; Yang, L.P.; Wang, L. Diagnosis on nitrogen status using GreenSeeker in spring maize. J. Maize Sci. 2008, 16, 111–114. [Google Scholar]
  36. Wang, L.; Bai, Y.L.; Lu, Y.L.; Wang, H.; Yang, L.P. NDVI analysis and yield estimation in winter wheat based on GreenSeeker. Acta Agron. Sin. 2012, 38, 747–753. [Google Scholar] [CrossRef]
  37. Ali, A.M.; Thind, H.S.; Singh, V.; Singh, B. A framework for refining nitrogen management in dry direct-seeded rice using GreenSeeker™ optical sensor. Comput. Electron. Agric. 2015, 110, 114–120. [Google Scholar] [CrossRef]
  38. Shaver, T.; Khosla, R.; Westfall, D. Evaluation of two crop canopy sensors for nitrogen recommendations in irrigated maize. J. Plant Nutr. 2014, 37, 406–419. [Google Scholar] [CrossRef]
  39. Zhang, H.C.; Gong, J.L. Research status and development discussion on high-yielding agronomy of mechanized planting rice in China. Acta Agron. Sin. 2014, 47, 1273–1289. [Google Scholar]
  40. Chen, Q.C.; Tian, Y.C.; Yao, X.; Zhu, Y.; Cao, W.X. Effects of nitrogen fertilizer top-dressing based on canopy reflectance spectra in rice. Sci. Agric. Sin. 2010, 43, 4149–4157. [Google Scholar]
  41. Calvo, R.N.I.; Rozas, H.S.; Echeverría, H.; Diovisalvi, N. Using canopy indices to quantify the economic optimum nitrogen rate in spring wheat. Agron. J. 2015, 107, 459–465. [Google Scholar] [CrossRef]
  42. Singh, B.; Singh, V.; Purba, J.; Sharma, R.K.; Jat, M.L.; Singh, Y.; Thind, H.S.; Gupta, R.K.; Chaudhary, O.P.; Chandna, P.; et al. Site-specific fertilizer nitrogen management in irrigated transplanted rice (Oryza sativa) using an optical sensor. Precis. Agric. 2015, 16, 1–21. [Google Scholar]
  43. Schmidt, J.; Beegle, D.; Zhu, Q.; Sripada, R. Improving in-season nitrogen recommendations for maize using an active sensor. Field Crops Res. 2011, 120, 94–101. [Google Scholar] [CrossRef]
  44. Yao, Y.X.; Miao, Y.X.; Huang, S.Y.; Gao, L.; Ma, X.B.; Zhao, G.M.; Jiang, R.F.; Chen, X.P.; Zhang, F.S.; Yu, K. Active canopy sensor-based precision N management strategy for rice. Agron. Sustain. Dev. 2012, 32, 925–933. [Google Scholar] [CrossRef]
  45. Goswami, S.; Gamon, J.; Vargas, S.; Tweedie, C. Relationships of NDVI, Biomass, and Leaf Area Index (LAI) for six key plant species in Barrow, Alaska. PeerJ PrePrints 2015, 3, e913v1. [Google Scholar]
Figure 1. Time series changes of NDVI value for two rice cultivars used in experiment 4 with N rates of 0, 75, 150, 225, 300, 375 kg N·ha−1. (a) Indica (YLY-1); (b) Japonica (WYJ-24).
Figure 1. Time series changes of NDVI value for two rice cultivars used in experiment 4 with N rates of 0, 75, 150, 225, 300, 375 kg N·ha−1. (a) Indica (YLY-1); (b) Japonica (WYJ-24).
Sensors 17 00672 g001
Figure 2. Quantitative relationships of NDVI to LAI (a), above-ground dry matter (b) and grain yield (c) in rice under varied N rates at jointing (JT), booting (BT) and heading (HD) growth stages.
Figure 2. Quantitative relationships of NDVI to LAI (a), above-ground dry matter (b) and grain yield (c) in rice under varied N rates at jointing (JT), booting (BT) and heading (HD) growth stages.
Sensors 17 00672 g002
Figure 3. Three types of model for RNDVI dynamics with RAGDD: (a) the double logistic model; (b) the rational equation; (c) the cubic polynomial equation.
Figure 3. Three types of model for RNDVI dynamics with RAGDD: (a) the double logistic model; (b) the rational equation; (c) the cubic polynomial equation.
Sensors 17 00672 g003
Figure 4. Comparisons of the double logistic simulated models of Indica (a) and Japonica (b) rice types based on different yield levels. RGDD1, RGDD2, and RGDD3 refer to the different RAGDDs with RNDVI peaks.
Figure 4. Comparisons of the double logistic simulated models of Indica (a) and Japonica (b) rice types based on different yield levels. RGDD1, RGDD2, and RGDD3 refer to the different RAGDDs with RNDVI peaks.
Sensors 17 00672 g004
Figure 5. Comparison of the RNDVI dynamics of Indica and Japonica rice types under the high-yield level (yield > 10.5 t·ha−1).
Figure 5. Comparison of the RNDVI dynamics of Indica and Japonica rice types under the high-yield level (yield > 10.5 t·ha−1).
Sensors 17 00672 g005
Figure 6. The relationships between the observed and simulated NDVI values of two rice cultivars (a) Japonica, (b) Indica, from tillering growth stage to flowering growth stage. The solid line is inclined at 45° to the axes.
Figure 6. The relationships between the observed and simulated NDVI values of two rice cultivars (a) Japonica, (b) Indica, from tillering growth stage to flowering growth stage. The solid line is inclined at 45° to the axes.
Sensors 17 00672 g006
Table 1. Basic information about five field experiments.
Table 1. Basic information about five field experiments.
Experiment NO.LocationCultivarN Rate (kg·ha–1)Plot Size (m2)Transplanting Date (Month/Day)Sampling Date (Month/Day)
EXP. 1 in 2008NanjingLYP9 (Indica)0, 110, 220, 3304.5 × 6.5 = 29.256/2307/16, 07/20, 07/26, 07/30, 08/03, 08/08, 08/13, 08/18, 08/23, 08/28, 08/31, 09/05
EXP. 2 in 2009NanjingLYP9 (Indica)0, 180, 3605.0 × 6.0 = 30.06/1707/15, 07/19, 07/25, 07/30, 08/04, 08/07, 08/13, 08/17, 08/22, 08/27, 09/02, 09/06
EXP. 3 in 2013RugaoWXJ14 (Japonica)0, 75, 150, 225, 300, 3755.0 × 6.0 = 30.06/2206/28, 07/03, 07/08, 07/11, 07/19, 08/08, 08/12, 08/15, 08/19, 08/25, 09/17, 09/21, 09/25, 10/02, 10/10, 10/11
SY63 (Indica)
EXP. 4 in 2014RugaoWXJ24 (Japonica)0, 75, 150, 225, 300, 3756.0 × 7.0 = 42.06/1706/23, 06/27, 06/30, 07/03, 07/06, 07/09, 07/17, 07/20, 07/24, 07/26, 07/29, 08/03, 08/06, 08/10, 08/16, 08/19, 08/23, 08/25, 08/30, 09/02, 09/05, 09/08, 09/16, 09/21, 09/25, 10/02, 10/06, 10/10, 10/14
YLY1 (Indica)
EXP. 5 in 2014RugaoWYJ24 LJ7 ZD110, 110, 220, 3305.0 × 6.0 = 30.06/177/18, 7/30, 08/06, 08/16, 08/26, 09/04
NJ4 (Japonica)
Rice cultivar: Wuxiangjing-14 (WXJ14), Wuyunjing-24 (WYJ24), Ningjing-4 (NJ4), Lianjing-7 (LJ7), Zhendao-11 (ZD11), Liangyoupei-9 (LYP9), Shanyou-63 (SY63), Y liangyou-1 (YLY1).
Table 2. The NDVImax value, days after transplanting (DAT) and AGDD when obtaining the NDVImax value for different rice cultivars and N treatments in experiments 3 and 4.
Table 2. The NDVImax value, days after transplanting (DAT) and AGDD when obtaining the NDVImax value for different rice cultivars and N treatments in experiments 3 and 4.
ItemsNDVImaxDAT(d)/AGDD (°C) When Obtaining the NDVImax Value
Year2013201420132014
CultivarWXJ14 (Japonica)SY63 (Indica)WYJ24 (Japonica)YLY1 (Indica)WXJ14 (Japonica)SY63 (Indica)WYJ24 (Japonica)YLY1 (Indica)
N Treatment00.489d0.689c0.663c0.688c65a/1285.5a65a/1155.5a70a/1123.5b69a/983.5a
750.632c0.759b0.737b0.736b65a/1285.5a59ab/1052b68b/1092.5c67b/956.5b
1500.645bc0.870a0.746b0.754ab59b/1170b59ab/1052b70a/1139a60d/869.5d
2250.659b0.880a0.769a0.769ab59b/1170b59ab/1052b64c/1030.5d63c/902.5c
3000.660b0.886a0.777a0.787a55c/1089.5c55b/979.5c64c/1030.5d63c/902.5c
3750.691a0.888a0.784a0.790a55c/1089.5c55b/979.5c60d/991.5e50e/736.5e
F-test statistical significance at 0.05 probability level.
Table 3. Coefficient of parameters, determination and RMSE of the RNDVI dynamic model with RAGDD.
Table 3. Coefficient of parameters, determination and RMSE of the RNDVI dynamic model with RAGDD.
Simulated ModelsParametersR2RMSE
abcd
y = ( 1 + e a × ( x b ) ) 1 ( 1 + e c × ( x d ) ) 1 15.28290.194411.65171.02670.85770.1161
y = a + bx + cx 2 + dx 3 –0.37965.7851–7.44372.70040.83570.1357
y = ( a + bx ) / ( cx + dx 2 + 1 ) –0.43195.3826–0.39486.11760.83190.1373
y = ( ab + cx d ) / ( b + x d ) –0.17410.00100.87203.54690.76710.1616
y = a / ( 1 + be cx ) 0.863597.944727.36410.75490.1674
x is the RAGDD, y is RNDVI.
Table 4. Parameters of the RNDVI dynamic model based on different cultivars and N rates for three yield levels.
Table 4. Parameters of the RNDVI dynamic model based on different cultivars and N rates for three yield levels.
Cultivar TypeYield LevelN RateCultivarYieldNDVImax ValueEntire Growing PeriodParameterR2RMSE
(t·ha–1)(kg·ha–1)(t·ha–1)(days)
JaponicaLow (yield ≤ 8.25 t·ha–1)0WXJ146.080.489150a: 16.45990.88340.1370
WYJ246.700.663156b: 0.3090
75WXJ147.750.632150c: 12.3144
WYJ247.870.737156d: 0.9851
Middle (8.25 t·ha–1 < yield < 10.5 t·ha–1)150WXJ148.780.645150a: 19.05440.90240.1224
WYJ248.980.746156b: 0.2629
225WXJ149.080.659150c: 11.4756
WYJ249.620.769156d: 1.0022
High (yield ≥ 10.5 t·ha–1)300WXJ1410.530.660150a: 20.03130.87640.1367
WYJ2410.540.777156b: 0.2370
375WXJ1410.610.691150c: 10.9741
WYJ2410.630.784156d: 1.0195
IndicaLow (yield ≤ 8.25 t·ha–1)0SY636.160.689153a: 14.36560.87130.1175
YLY17.010.688133b: 0.2196
75SY638.170.759153c: 14.0343
YLY18.200.736133d: 0.9972
Middle (8.25 t·ha–1 < yield < 10.5 t·ha–1)150SY639.040.870153a: 17.10280.86100.1128
YLY19.250.754133b: 0.1749
225SY639.860.880153c: 13.2413
YLY110.130.769133d: 1.0192
High (yield ≥10.5 t·kg·ha–1)300SY6310.610.886153a: 23.82610.88740.0981
YLY110.830.787133b: 0.1489
375SY6311.020.888153c: 12.0923
YLY111.260.790133d: 1.0361
Table 5. Coefficient of k, determination (R2) and RMSE of the linear correlation between the observed and simulated NDVI values at different growth stages.
Table 5. Coefficient of k, determination (R2) and RMSE of the linear correlation between the observed and simulated NDVI values at different growth stages.
Growth StagekR2RMSE
JaponicaIndicaJaponicaIndicaJaponicaIndica
Active tillering0.97491.01580.7044 **0.6102 **0.03050.0128
Middle tillering1.01871.03180.7689 **0.7990 **0.01640.0075
Jointing1.00451.03300.9331 **0.6656 **0.00790.0105
Booting1.01601.00980.6565 **0.7367 **0.01910.0113
Heading1.00860.98590.9167 **0.8211 **0.01740.0064
Flowering0.96921.03330.8762 **0.6357 **0.01190.0128
Active tillering to flowering0.99911.01700.9084 **0.8030 **0.02320.0170
** F-test statistical significance at 0.05 probability level.

Share and Cite

MDPI and ACS Style

Liu, X.; Ferguson, R.B.; Zheng, H.; Cao, Q.; Tian, Y.; Cao, W.; Zhu, Y. Using an Active-Optical Sensor to Develop an Optimal NDVI Dynamic Model for High-Yield Rice Production (Yangtze, China). Sensors 2017, 17, 672. https://doi.org/10.3390/s17040672

AMA Style

Liu X, Ferguson RB, Zheng H, Cao Q, Tian Y, Cao W, Zhu Y. Using an Active-Optical Sensor to Develop an Optimal NDVI Dynamic Model for High-Yield Rice Production (Yangtze, China). Sensors. 2017; 17(4):672. https://doi.org/10.3390/s17040672

Chicago/Turabian Style

Liu, Xiaojun, Richard B. Ferguson, Hengbiao Zheng, Qiang Cao, Yongchao Tian, Weixing Cao, and Yan Zhu. 2017. "Using an Active-Optical Sensor to Develop an Optimal NDVI Dynamic Model for High-Yield Rice Production (Yangtze, China)" Sensors 17, no. 4: 672. https://doi.org/10.3390/s17040672

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop