Monitoring Hybrid Rice Phenology at Initial Heading Stage Based on Low-Altitude Remote Sensing Data

: Accurate monitoring of hybrid rice phenology (RP) is crucial for breeding rice cultivars and controlling fertilizing amount. The aim of this study is to monitor the exact date of hybrid rice initial heading stage (IHS DAS ) based on low-altitude remote sensing data and analyze the inﬂuence factors of RP. In this study, six ﬁeld experiments were conducted in Ezhou city and Lingshui city from 2016 to 2019, which involved different rice cultivars and nitrogen rates. Three low-altitude remote sensing platforms were used to collect rice canopy reﬂectance. Firstly, we compared the performance of normalized difference vegetation index (NDVI) and red edge chlorophyll index (CIred edge) for monitoring RP. Secondly, double logistic function (DLF), asymmetric gauss function (AGF), and symmetric gauss function (SGF) were used to ﬁt time-series CIred edge for acquiring phenological curves (PC), the feature: maximum curvature (MC) of PC was extracted to monitor IHS DAS . Finally, we analyzed the inﬂuence of rice cultivars, N rates, and air temperature on RP. The results indicated that CIred edge was more appropriate than NDVI for monitoring RP without saturation problem. Compared with DLF and AGF, SGF could ﬁt CIred edge without over ﬁtting problem. MC of SGF_CIred edge from all three platforms showed good performance in monitoring IHS DAS with good robustness, R 2 varied between 0.82 and 0.95, RMSE ranged from 2.31 to 3.81. In addition, the results demonstrated that high air temperature might cause a decrease of IHS DAS , and the growth process of rice was delayed when more nitrogen fertilizer was applied before IHS DAS . study illustrated that low-altitude remote sensing technology could be used for monitoring ﬁeld-scale hybrid rice IHS DAS accurately.


Introduction
With the rapid growth of the world population, food security has become an urgent issue, attracting worldwide attention. Rice is a very important crop that provides staple food for almost 50% of the world's population [1]. Hybrid vigor (heterosis) is a universal phenomenon in crops, and the yield of hybrid rice is usually higher than that of conventional rice [2,3]. Hybrid rice makes a great contribution to the world's food security, with more than half the growth of rice yield [4][5][6]. Hybrid rice is widely planted in large areas of Asian region, including China, Japan, and Southeast Asian countries. Breeding hybrid rice cultivars with high yield and good quality has been an important part of these countries' strategies for developing modernized agriculture [7].
Rice phenology (RP) is critical for farm management and yield evaluation [8]. RP is divided into germination, seeding, tillering, jointing, booting, heading and flowering, grain of 0.18. Bai et al. [49] utilized the fixed towers mounted with RGB cameras to achieve rice heading stage automatic observation. Lin et al. [50] calculated NDVI from a handheld radiometer to monitor RP, the results showed that the maximum NDVI had a good correlation ship with the heading and flowering stage. Han et al. [51] proposed a deep convolutional neural network (DCNN) approach to analyze hand-held RGB images for real-time RP detection, with a high overall accuracy of 0.913 and low mean absolute error of 0.090. In addition, Zheng et al. [9] compared the ability of two ground-based spectral instruments for monitoring RP, there were 2.3-9.3 days difference between the observed date and estimated date at different growth stages. In a word, these studies proved that low-altitude remote sensing sensors had good performance in monitoring field-scale RP.
To our knowledge, there are three key steps to monitor vegetation phenology based on remote sensing technology: (1) cleaning and flagging remote sensing data; (2) reconstructing vegetation phenological curves; (3) extracting features for monitoring vegetation phenology [52]. However, previous studies mainly focused on studying the phenological characteristics of rice and exploring appropriate methods for monitoring RP. Few studies compared the advantages and disadvantages of different low-altitude remote sensing sensors for field-scale RP detection, and analyzed the influence factors of RP. Thus, the aim of this study is to (1) analyze the change regularity of typical time-series VIs (NDVI and CIred edge) during the whole growth stage of hybrid rice; (2) evaluate the performance of three functions for reconstructing rice phenological curves based on three low-altitude remote sensing data, and extract features of VIs for monitoring IHS DAS ; (3) analyze the influential factors (rice cultivars, N rates, and air temperature) of RP.

Experimental Design
Experiments 1, 2, and 3 were conducted in three growing seasons spanning from 2016 to 2018. Rice seeds were sown in seedling beds. 30 days after sowing (DAS), seedlings were transplanted into rice field with a unified density 225,000 plants/hm 2 . The plot area of each rice cultivar was 40 m 2 (5.0 m × 8.0 m). For all treatments, N fertilizer (180 kg/hm 2 ) was applied in the form of urea: 50% as basal fertilizer before transplanting, 25% at the tillering stage and 25% at the booting stage. Monocalcium phosphate (P 2 O 5 90 kg/hm 2 ) and potassium chloride (K 2 O 180 kg/hm 2 ) were applied before transplanting. Other management decisions (pest control and herbicide application) followed the local standard practices. Experiment 4 was a randomized complete block design with three replications for each treatment, involving two rice cultivars and four N rates. Rice cultivars were Fengliangyou 4 (FLY4) and Luoyou 9348 (LY9348). Four N rates: 0 (N0), 120 (N1), 180 (N2) and 240 (N3) kg/hm 2 were applied in the form of urea: 50% as basal fertilizer before transplanting, 50% at the booting stage. The 24 individual plot area was 15 m 2 (5.0 m × 3.0 m), respectively. In addition, other factors were the same with experiments 1, 2, and 3. Experiment 5 and 6 were conducted from May to October of 2018 and 2019, involving 1014 and 289 rice cultivars, each rice cultivar was planted in an individual plot, and the area of plot was 1 m 2 (1.0 m × 1.0 m), respectively. In addition, other factors were the same with experiments 1, 2, and 3.
Six experiments involved different hybrid rice cultivars and N rates, as detailed in Figure 2 and Table 1.  Meteorological data were collected from the China Meteorological Data Service Center (http://www.weather.com.cn/), which included maximum air temperature (T max ), minimum air temperature (T min ) and weather conditions spanning from 2016 to 2019. In order to show weather conditions clearly, we used different numbers to represent different weather conditions: 1 represented sunny, 2 represented cloudy, 3 represented overcast, 4 represented shower, 5 represented light rain, 6 represented moderate rain, 7 represented heavy rain, and 8 represented rainstorm, respectively.

Field Phenological Data
In our study, the phenological time was recorded as days after sowing (DAS). Ground field observations and records were carried out according to the specifications for agrometeorological observation-Rice (http://www.cma.gov.cn/root7/auto13139/201903/t20190 329_519169.html) from the China Meteorological Administration. For example, when 10% of investigated rice grew spikes, the exact date was defined as the time of rice phenology at initial heading stage (IHS DAS ). Then, we investigated 40 hills of rice plants in experiments 1, 2, 3, and 4, and 20 hills of rice plants in experiments 5 and 6 were investigated by random method, respectively.

Daily Reflectance
Three sets of SKYE radiometers (SKR 1860, SKYE Instruments Ltd., Llandrindod Wells, UK) were employed to collect four bands reflectance every 30 min. In order to distinguish SKYE radiometers, we marked them as A_SKYE, B_SKYE and C_SKYE. Each set of SKYE radiometers consisted of two sensors: the upward radiometer sensor (URS) and the downward radiometer sensor (DRS). URS measured downwelling (DW) radiance with the aid of a cosine corrector, which collected radiance from sky in 180 • hemisphere above the sensor. DRS measured the 25 • field of view upwelling (UW) radiance from rice canopy. During the whole growth stage of rice, all SKYE radiometers were positioned above rice canopy at a height of 2 m with a 25 • field-of-view of 0.62 m 2 area ( Figure 3). The specifications of SKYE radiometers are listed in Table 2.
where UW λ and DW λ are calculated (mol m −2 s −1 ), respectively. In four growing seasons spanning from 2016 to 2019, SKYE radiometers collected continuously reflectance of three rice cultivars: LY 9348, FLY 4, and Yi Liang You Hua Zhan (YLYHZ), as detailed in Table 3. The median reflectance between 10:00 a.m. to 2:00 p.m. was calculated as daily reflectance to reduce the impact of weather drastic change.

The Hyperspectral Reflectance
The hyper-spectrometer ASD Field Spec 4 (ASD Inc., Boulder, CO, USA) was employed to collect rice canopy reflectance in experiment 1, 2, 3 and 4, with the spectrum ranging from 350-2500 nm. The measurements were conducted on sunny days from 10:00 am to 2:00 pm. A white reference panel (Chemical composition: BaSO 4 , Size: 25.4 × 25.4 cm) was used to make relative radiometric correction before each measurement. The optical fiber probe of ASD was positioned above rice canopy at a height of 1 m with a 25 • field-of-view of 0.15 m 2 area. Five fixed points' spectral measurements were averaged to reduce uncertainty error of spectral measurement in each plot. The hyperspectral reflectance involved different rice cultivars, as detailed in Table 4.

The UAV Multispectral Reflectance
In this study, a multispectral camera (MCA, Tetracam Inc., Chatsworth, CA, USA) was mounted on an eight-propeller UAV (S1000, DJI Technology Co., Ltd., Shenzhen, China) to acquire the panoramic multispectral images in experiment 5 and 6. The UAV flight was conducted on sunny days between 10:00 a.m. to 2:00 p.m. The flight duration varied between 10 and 15 min. The flight height was 120 m with a field of view 83.25 m × 66.50 m, with the ground sampling distance 6.5 cm/pixel.
MCA consisted of twelve individual digital miniature cameras equipped with different band pass filters centered at the wavelength of 490 nm, 520 nm, 550 nm, 570 nm, 670 nm, 680 nm, 700 nm, 720 nm, 800 nm, 850 nm, 900 nm, and 950 nm respectively ( Table 5). The multispectral images are recorded as digital number (DN) values. Prior to the flight, six calibration blankets were placed on the flat road as radiometric calibration targets ( Figure 4). One panoramic multispectral image could simultaneously contain the six calibration blankets and the total field plots within the field of view of MCA. The calibration blankets were made with highly durable woven polyester fabric at the size of 1.2 m × 1.2 m, having the relatively constant reflectance of 0.03, 0.12, 0.24, 0.36, 0.56, and 0.80, respectively. As a linear relationship was assumed between surface reflectance and at-sensor radiance, the empirical line method (ELM) was used to convert DN to surface reflectance [53,54]. Gain k and o f f set k were calculated using the least-square method from the relatively constant reflectance (0.03, 0.12, 0.24, 0.36, 0.56, and 0.80) and the corresponding DN values of calibration targets. Then Gain k and o f f set k were applied to calculate the canopy surface of each pixel in multispectral images. The specific equations are as follows.
where k is the center wavelengths of 12 bands, DN(i,k) are the digital number values of the calibration target i (i = 1, 2, 3, 4, 5, 6, respectively) in band k. Gain k and o f f set k represent the slope and intercept. R k and DN k are the canopy surface reflectance and the digital number values of multispectral images. The UAV system collected multispectral images from May to October of 2018 and 2019, involving 1303 rice cultivars in total (Table 6). PixelWrench2 software (Tetracam Inc., Chatsworth, CA, USA) was used to perform images preprocessing, including noise correction, vignetting correction, lens distortion correction, band-to-band registration and band stacking [55,56]. The radiometric calibration and reflectance extraction processing were performed using the ENVI software (Exelis Visual Information Solutions, Boulder, CO, USA).  After image preprocessing, we selected a region of interest (ROI) in the multispectral image, which included 100 pixels of 0.42 m 2 area in each plot. The plot-level reflectance was calculated as the average reflectance of all pixels within ROI.

Vegetation Indices (VIs)
Vegetation indices (VIs) refer to the spectral reflectance transformations of two or more bands, which are related to vegetation photosynthesis and nutrition status [57]. NDVI is the most commonly used VIs for quantitative remote sensing [58]. CIred edge has good performance in the estimation of chlorophyll content, leaf area index and biomass [59]. In this study, NDVI and CIred edge are employed to monitor RP. They are calculated as: where ρ N IR is the reflectance in the wavelength of 865 nm in SKYE (800 nm in ASD and MCA), ρ Red edge is the reflectance in the wavelength of 717 nm in SKYE (720 nm in ASD and MCA), ρ Red is the reflectance in the wavelength of 655 nm in SKYE (680 nm in ASD and MCA).
The calculation formulas are as follows: where DAS are the days after rice sowing, corresponding the time of VIs acquisition, e is the natural constant, a, b, c, d, f, g, h, i, j, k, l, m, p, q, and r are the best parameters of three fitting functions (DLF, AGF, and SGF), respectively.
where PC are the phenological curves: DLF_VIs, AGF_VIs, and SGF_VIs, PC is the first derivative of PC, and PC" is the second derivative of PC. MC is the maximum curvature of PC ( Figure 5).  Figure 6 shows the flowchart. Firstly, we calculated the reflectance of different remote sensing platforms. The statistical analysis of meteorological data and field phenological data were carried out. Secondly, daily NDVI and CIred edge were computed using the SKYE reflectance. We analyzed the phenological change of rice during the whole growth stage, and compared the characteristics of NDVI and CIred edge for monitoring RP. Thirdly, the VIs (NDVI or CIred edge) were computed using the reflectance of ASD and MCA. Three fitting functions (DLF, AGF and SGF) were applied to fit the suitable VIs for reconstructing PC: DLF_VIs, AGF_Vis, and SGF_VIs. We evaluated the performance of three fitting functions. Finally, we extracted the MC of PC for monitoring IHS DAS and evaluated the accuracy.

The Assessment of Models
The coefficient of determination (R 2 ) and root mean square error (RMSE) were used to evaluate the performances of three fitting functions and the accuracy of monitoring IHS DAS . Mathematically, a higher R 2 and smaller RMSE represent better model accuracy.
where x i and y i are the estimated and measured values, y is the average of the measured values, and n is the sample number, respectively. The mathematical models were performed using MATLAB R2018a (MathWorks, Inc., Natick, MA, USA). Graphics were prepared using Origin 9.1 software program (OriginLab Corporation, Northampton, MA, USA). Figure 7 shows that the weather conditions have obvious difference between two climatic zones, the climate of Lingshui city is more stable than that of Ezhou city. In Ezhou city, the weather conditions in experiment 1 and 2 fluctuate greatly. It often rains during the whole growth period of rice (Figure 7a,c). Sometimes, there is even a regional heavy rain or rainstorm. Further, the continuous high temperature (T max ≥ 35 • C) appears between DAS 76 to 88 and DAS 64 to 82 (Figure 7b,d). The weather conditions in experiment 5 and 6 are mainly sunny to cloudy with light rain occasionally (Figure 7g,i), however, the continuous high temperature (T max ≥ 35 • C) appears between DAS 49 to 79 and DAS 70 to 105 (Figure 7h,j).

Statistical Analysis of Meteorological Data
In Lingshui city, the weather conditions in experiments 3 and 4 are very stable. It is mainly sunny to cloudy with small air temperature change (Figure 7,f), there is not continuous high air temperature appearing during the whole growth period of rice. Good weather conditions in Lingshui city are helpful to obtain accurate remote sensing data.

Statistical Analysis of IHS DAS
We carried out statistical analysis of IHS DAS in six experiments (

Comparative Analysis of Daily NDVI and CIred Edge
The SKYE radiometers were used to acquire daily reflectance data of 9 plots in four rice growing seasons. Figure 8 shows that the change regularity of NDVI is different from that of CIred edge. NDVI increases rapidly and achieves the maximum value quickly. Before DAS 70, NDVI increases rapidly; between DAS 80 to 100, NDVI remains stable and is close to 0.9. When rice is at maturity stage (after DAS 100), NDVI declines slowly with the senescence of leaves and stems, and it is close to 0.5 before harvesting. Especially, as Figure 8b-d present, NDVI achieves the maximum value before DAS 55 corresponding to the middle stage of tillering, which indicates that rice grows rapidly after transplanting in field. There is an obvious saturation problem occurring before IHS DAS (the green lines in Figure 8). Thus, it is difficult to extract features from NDVI for monitoring IHS DAS .
The curves of CIred edge are different in four rice growing seasons, but the change regularity of CIred edge is similar in the same growing season. CIred edge gradually increases with an exponential growth at tillering and jointing stage. CIred edge achieves the maximum value before IHS DAS . The maximum value of CIred edge does not exceed 4 without saturation problem. CIred edge remains stable and is close to the maximum value before IHS DAS . As shown in Figure 8, there are obvious inflection points of CIred edge appearing at IHS DAS . It suggests that we can take advantage of the feature: the MC of time-series CIred edge for monitoring IHS DAS . After IHS DAS , CIred edge decreases linearly with the senescence of leaves and stems. In addition, there is some noise in time-series CIred edge. For example, Figure 8b-d display abnormal fluctuations: two peaks and one trough appear between DAS 60 to 95. Thus, it is desirable to use some smoothing and fitting methods to eliminate noise of time-series CIred edge.
The comparative analysis indicates that CIred edge presents no saturation problem during the whole growth period of rice. CIred edge is better than NDVI for field-scale RP detection. Meanwhile, CIred edge has the obvious feature (MC) for monitoring IHS DAS .

Fitting CIred Edge of SKYE
Three fitting functions (DLF, AGF and SGF) are applied to fit CIred edge of SKYE in four growing seasons. Figure 9 illustrates the relationship between the original CIred edge and the fitted phenological curves. The results show that DLF, AGF and SGF can eliminate noise of the original CIred edge and help to present the change regularity of time-series CIred edge accurately. However, the performance of fitting models varies at different growth stages. DLF has good performance in fitting CIred edge at the early stages (before DAS 80). AGF is suitable to fit CIred edge between DAS 80 to 100, and SGF shows the highest accuracy after rice begins to grow into maturity (after DAS 100). Table 8 summarizes the statistics (R 2 and RMSE) of three fitting models. For DLF models, R 2 varies between 0.94 and 0.98, RMSE ranges from 0.07 to 0.19. The model b shows the worst performance among DLF models, whereas the model e displays the best performance among DLF models. For AGF models, R 2 varies between 0.95 and 0.98, and RMSE ranges from 0.07 to 0.16. The model d shows the worst performance among AGF models, whereas the model f displays the best performance among AGF models. For SGF models, R 2 varies between 0.88 and 0.95, RMSE ranges from 0.13 to 0.36. The model b shows the worst performance among AGF models, whereas the model e displays the best performance among SGF models. Overall, compared with AGF and SGF, DLF proves to have good performance in fitting CIred edge of SKYE during the whole growth period of rice.

Fitting CIred Edge of ASD
The performances of fitting models in four experiments are shown in Table 9. In experiment 1, R 2 and RMSE vary from 0.91-1 and 0-0.11. In experiment 2, R 2 and RMSE vary from 0.25-1 and 0-0.44. In experiment 3, R 2 and RMSE vary from 0.92-1 and 0-0.20. In experiment 4, R 2 and RMSE vary from 0.95-1 and 0-0.17. Table 9. The performance of the models in fitting CIred edge of ASD. We acquire 4, 6, 6, and 7 batches of CIred edge in experiments 1, 2, 3, and 4 (Table 4). There are six, six, and three parameters that need to be optimized in DLF, AGF, and SGF models. When the frequency of times for acquiring CIred edge is low, the over fitting problem of DLF and AGF models appears (R 2 = 1, RMSE = 0). Compared with DLF and AGF models, the accuracy of SGF models is lower in four experiments. However, there are no obvious over fitting problem appearing in SGF models. The results suggest that the fitting functions with fewer parameters, such as SGF, maybe suitable to avoid over fitting problem. SGF can reconstruct vegetation phenological curves using low-frequency CIred edge. Table 10 summarizes the performance of DLF, AGF and SGF. For DLF models, R 2 varies between 0.18 and 0.99, RMSE ranges from 0.015 to 0.86. For AGF models, R 2 varies between 0.01 and 0.99, RMSE ranges from 0.02 to 1.27. For SGF models, R 2 varies between 0.24 and 0.99, RMSE ranges from 0.04 to 0.96. In addition, the performance evaluation implies that the fitting phenological curves (DLF_CIred edge, AGF_CIred edge and SGF_CIred edge) not only have a good correlation ship with the original CIred edge, but also can avoid over fitting problem (R 2 = 1, RMSE = 0). In summary, compared with AGF and SGF, DLF is more suitable to fit time-series CIred edge of MCA, with higher R 2 and lower RMSE. This result is consistent with the result of SKYE data.

Effects of Rice Cultivars and N Rates on IHS DAS
In order to understand the influence of rice cultivars and N rates on IHS DAS , experiment 4 was conducted in a randomized complete block designed with three replications for each treatment, involving two rice cultivars (FLY4 and LY9348) and four N rates (N0, N1, N2 and N3). The analysis of variance (ANOVA) indicates that rice cultivars and N rates have obvious effects on IHS DAS (Table 11). There is a significant difference between LY9348 and FLY4. The difference of IHS DAS reaches an extremely significant level under different N rates. The interaction of different rice cultivars and N rates on IHS DAS reaches a significant level. The ANOVA suggests that we would better take into account rice cultivars and N rates when monitoring IHS DAS .  [70]. In general, when using more nitrogen fertilizer before heading and flowering stage, the rice growth process might be delayed, leaves would remain green, and spikes appeared later, resulting in a longer life cycle of rice.

Daily CIred Edge for Monitoring RP
Previous studies showed that time-series NDVI was suitable for monitoring vegetation phenology [71,72]. However, the major problem of NDVI was that it became saturated for vegetation with moderate to high density [73]. Our study confirmed this fact: daily NDVI of SKYE became significantly saturated before flowers and spikes appeared in rice. Thus, it was difficult to extract features from NDVI for monitoring IHS DAS . In comparison with NDVI, daily CIred edge of SKYE had no saturation problem during the whole rice growth stage, and was more suitable for monitoring IHS DAS . Figure 8 showed that the change regularity of daily CIred edge was closely related to the actual growth process of rice. To our knowledge, CIred edge was originally proposed for high accuracy chlorophyll estimation [74], some studies have demonstrated that the canopy chlorophyll content was closely related to vegetation growth period [75,76]. Therefore CIred edge was better than NDVI for monitoring RP, and the result was consistent with the previous study [9].
In addition, the maximum air temperature (Tmax) has a significant effect on daily CIred edge. Figure 14 displays daily ρ N IR and ρ Red edge of C_SKYE and Tmax in Ezhou city and Lingshui city. In Lingshui city, Tmax is stable and always below 35 • C. However, the continuous high air temperature (Tmax ≥ 35 • C) often appears with drastic change in Ezhou city. The abnormal change of ρ N IR often occurs between DAS 60 to 95 when Tmax has drastic change ( Figure 14a). Meanwhile, the abnormal CIred edge also appeared between DAS 60 to 95 (Figure 8d). In general, rice canopy reflectance presents a spectral curve of typical green vegetation before IHS DAS , ρ Red edge is much lower than ρ N IR during different growth stage [77,78]. In our study, compared with ρ N IR , ρ Red edge slowly increases with the progress of rice growth and is not sensitive to high air temperature. However, the drastic change of Tmax could cause the sudden change of ρ N IR , resulting in the abnormal fluctuation of CIred edge appearing. Rice is a crop that is sensitive to high temperature [79] and water [80]. The continuous high air temperature or additional rain had great influence on the growth process of rice, IHS DAS of rice may be advanced or postponed [81]. Moreover, strong sunlight can burn rice leaves and change the rice canopy structure. Thus, it is better to take measures (smooth and fitting methods) to eliminate noise before time-series CIred edge was applied to monitor RP. VIs were indeed an effective tool for monitoring vegetation phenology [82,83]. However, if we can make better use of other information, for example, texture, threedimensional information and coverage, the accuracy of monitoring RP will be improved in the future [84].

Comparative Analysis of DLF, AGF and SGF for Fitting Several Source CIred Edge
In order to obtain high quality time-series CIred edge, we used DLF, AGF, and SGF to fit CIred edge. The results indicated that DLF and AGF proved to have good performance in fitting CIred edge of SKYE and MCA with higher R 2 and lower RMSE. However, when the frequency of times for acquiring CIred edge is low, Table 9 shows that the over fitting problem (R 2 = 1, RMSE = 0) of DLF and AGF models appears. Compared with DLF and AGF models, SGF models have a lower correlation ship with the original CIred edge of ASD, but without obvious over fitting problem.
There are six, six, and three parameters that need to be optimized in DLF, AGF, and SGF models. DLF and AGF are relatively complex, more parameters mean that more CIred edge data should be put into DLF and AGF models. Low-frequency measured data was more likely to cause over fitting problem when more parameters need to be optimized in fitting functions. As we know, the over fitting problem may misinterpret important phenological features and directly affect the accuracy of monitoring vegetation phenology [52]. Thus, DLF and AGF are not suitable to fit low-frequency CIred edge. Acquiring high-frequency CIred edge will help to highlight the phenological features for monitoring RP.
Overall, DLF and AGF could be effectively applied to fit dense CIred edge with higher R 2 and lower RMSE, and to describe more detailed changes of rice canopy at different growth stages. SGF had good performance in fitting low-frequency CIred edge.

Monitoring IHS DAS Based on Several Source CIred Edge
Time-series CIred edge from SKYE, ASD and MCA were reconstructed by DLF, AGF and SGF. We extracted the features: MC of DLF_CIred edge, AGF_CIred edge and SGF_CIred edge for monitoring IHS DAS .
MC of DLF_CIred edge and AGF_CIred edge from SKYE had good correlation with IHS DAS . However, MC of DLF_CIred edge and AGF_CIred edge from ASD and MCA did not perform well for monitoring IHS DAS . MC of SGF_CIred edge from three platforms all have good performance in monitoring IHS DAS , R 2 varies between 0.82 and 0.95, RMSE ranges from 2.31 to 3.81, which means that MC of SGF_CIred edge could be applied to monitor IHS DAS with good robustness. The rice flowers and spikes were usually bright yellow with contrasting colour to the rice leaves [49]. The chlorophyll content of flowers and spikes were also lower than that of green leaves. After flowers and spikes came out, the canopy chlorophyll of rice would decrease dramatically, resulting in the obvious inflection points of CIred edge appearing at IHS DAS . Thereby, the MC of time-series CIred edge was an essential feature for monitoring field-scale IHS DAS .
In this study, the accuracy of monitoring IHS DAS was different using MC extracted from three low-altitude remote sensing data. Three different platforms have the advantages and disadvantages. SKYE radiometers could continuously collect in situ field-scale rice canopy reflectance at very high temporal resolution, which helped analyze the change regularity of phenological curves and extract features for monitoring RP. However, they were mounted in the fixed positions for monitoring several fields, lacking mobility and flexibility [85]. ASD had wider spectrum range and high radiation resolution for extracting accurate VIs, but it was time-consuming to acquire the hyperspectral reflectance of different rice cultivars. The multispectral UAV system had the advantages of mobility, it could collect multispectral images for monitoring real-time large-scale vegetation growth [56,[86][87][88]. However, the processing flow of multispectral image was relatively complicated, including radiometric calibration, geometric correction, and image mosaic. It is necessary to acquire the high spatiotemporal and hyperspectral data for improving the accuracy of monitoring field-scale IHS DAS .

The Influence Factors on IHS DAS
There were many factors influencing vegetation phenology, such as illumination, air temperature, vegetation varieties and fertilizer application [89,90]. In our study, we mainly took into account rice cultivars, N rates and air temperature for monitoring IHS DAS . Table 11 indicates that rice cultivars and N rates have obvious effects on IHS DAS . The interaction of different rice cultivars and N rates on the IHS DAS reaches significant level. With the increase of N rates, IHS DAS of LY 9348 and FLY 4 both increase, which indicates that the growth process of rice is delayed when more nitrogen fertilizer is applied before heading stage, and leaves would remain green, flowers and spikes appear later [70]. In a word, rice cultivars and nitrogen fertilizer application have great influence on RP [91,92].
Air temperature could not only affect the near infrared reflectance, but also disturb RP. Compared with Lingshui city, the weather of Ezhou city fluctuates greatly, it often remains continuous high temperature (T max ≥ 35 • C). Table 3 shows that IHS DAS of LY 9348 and FLY 4 planted in Ezhou city are lower than that in Lingshui city. The results confirm that high air temperature might accelerate the whole growth process of rice and shorten growth duration [93][94][95]. The IHS DAS of hybrid rice might decrease when the same rice cultivars are planted in high temperature area.

Conclusions
In this study, we presented the processing steps for monitoring hybrid rice IHS DAS based on three low-altitude remote sensing data and analyzed the influential factors (rice cultivars, N rates, and air temperature) of RP.
The results indicated that CIred edge could show hybrid rice phenological change accurately without saturation problem, and was more appropriate than NDVI for monitoring IHS DAS . SGF could eliminate the noise of time series CIred edge without over fitting problem. MC of SGF_CIred edge from three low-altitude remote sensing platforms had good performance in monitoring IHS DAS with good robustness, as R 2 varied between 0.82 and 0.95 and RMSE ranged from 2.31 to 3.81.
This study demonstrated that high air temperature might cause a decrease of IHS DAS when the same rice cultivars were planted in high temperature area. In addition, the interaction of different rice cultivars and N rates on IHS DAS reached a significant level. The growth process of rice was delayed when more nitrogen fertilizer was applied before IHS DAS .
Future work should evaluate these results for diverse rice cultivars under different environmental conditions, by combing low-altitude remote sensing data and satellite data. In addition, new methods for reconstructing rice phenological curves will be explored for monitoring RP in the future. We expect that our study will provide useful guidance for breeding rice and fertilizer management with remote sensing technology.