Assessing the Operation Parameters of a Low-Altitude UAV for the Collection of NDVI Values over a Paddy Rice Field

: Unmanned aerial vehicle (UAV) remote sensing platforms allow for normalized di ﬀ erence vegetation index (NDVI) values to be mapped with a relatively high resolution, therefore enabling an unforeseeable ability to evaluate the inﬂuence of the operation parameters on the quality of the thus acquired data. In order to better understand the e ﬀ ects of these parameters, we made a comprehensive evaluation on the e ﬀ ects of the solar zenith angle (SZA), the time of day (TOD), the ﬂight altitude (FA) and the growth level of paddy rice at a pixel-scale on UAV-acquired NDVI values. Our results show that: (1) there was an inverse relationship between the FA ( ≤ 100 m) and the mean NDVI values, (2) TOD and SZA had a greater impact on UAV–NDVIs than the FA and the growth level; (3) Better growth levels of rice—measured using the NDVI—could reduce the e ﬀ ects of the FA, TOD and SZA. We expect that our results could be used to better plan ﬂight campaigns that aim to collect NDVI values over paddy rice ﬁelds.


Introduction
Paddy rice is a staple crop for more than half the world's population, especially in Asia [1]. China is the largest paddy rice-producing country, but still faces challenges from a rapid increase in global food demand and food security [2]. Hence, it is of utmost importance to increase food productivity in a sustainable manner. Successful and timely estimations of paddy-rice growth are useful means of cultivation management, agricultural decision-making and tillage improvement [3]. Cultivating paddy rice scientifically can contribute to disease prevention, optimal nutrient use, yield increments and the reduction of environmental pollution arising due to excessive nitrogen use [4][5][6]. One essential aspect of effective field management practices is to monitor crop growth with higher efficiency.

UAV Data Collection
A Micasense RedEdge-M camera (MicaSense, Inc., Seattle, WA, USA) was attached to the Matrice-600 Pro UAV (DJI, Inc., Shenzhen, China), a six-rotor UAV, for the aerial imagery collection. This camera has five independent imagers, each with a custom narrowband filter that enables the imager to receive a precise spectrum. The RedEdge-M camera also includes a downwelling light sensor (DLS) and a GPS module. The DLS is a light sensor with five bands, which need to connect directly to the host for radiation calibration. The specifications of the RedEdge-M camera are given in Table 1.
The flight plans were designed using the Atlas Flight application (Micasense, Inc., Fremont, CA, USA). The flight altitudes were set to 40 m, 60 m, 80 m and 100 m and the flight speed was set to 5 m/s. In addition. the side overlap, and forward overlap of the camera were set to 80% to ensure effective orthorectification and mosaicking of the resulting images. Following the instructions of the Micasense manual (https://support.micasense.com), a calibrated Teflon panel was imaged by the RedEdge-M camera for radiometric calibration before and after each flight. Using these settings, the image acquisition was conducted above the study area from 6 am to 6 pm on June 11, 2018. A total of 48 flights were thus conducted and 31,980 images were collected by performing each flight mission at an interval of 15 minutes. Detailed flight parameters are given in Error! Reference source not found..

UAV Data Collection
A Micasense RedEdge-M camera (MicaSense, Inc., Seattle, WA, USA) was attached to the Matrice-600 Pro UAV (DJI, Inc., Shenzhen, China), a six-rotor UAV, for the aerial imagery collection. This camera has five independent imagers, each with a custom narrowband filter that enables the imager to receive a precise spectrum. The RedEdge-M camera also includes a downwelling light sensor (DLS) and a GPS module. The DLS is a light sensor with five bands, which need to connect directly to the host for radiation calibration. The specifications of the RedEdge-M camera are given in Table 1. The flight plans were designed using the Atlas Flight application (Micasense, Inc., Fremont, CA, USA). The flight altitudes were set to 40 m, 60 m, 80 m and 100 m and the flight speed was set to 5 m/s. In addition. the side overlap, and forward overlap of the camera were set to 80% to ensure effective orthorectification and mosaicking of the resulting images. Following the instructions of the Micasense manual (https://support.micasense.com), a calibrated Teflon panel was imaged by the RedEdge-M camera for radiometric calibration before and after each flight. Using these settings, the image acquisition was conducted above the study area from 6 am to 6 pm on 11 June 2018. A total of 48 flights were thus conducted and 31,980 images were collected by performing each flight mission at an interval of 15 min. Detailed flight parameters are given in Table 2.

Ground Data Collection
The ground NDVI data were obtained using the GreenSeeker handheld crop sensor (Trimble, Inc., Sunnyvale, CA, USA), which has two independent LEDs that can emit 774 nm near-infrared (NIR) light and 656 nm red light. Due to its own light source, the sensor is relatively unaffected by solar conditions. In this study, the GreenSeeker sensor was used to acquire NDVI ground data (called GS-NDVIs from Remote Sens. 2020, 12, 1850 5 of 16 here forward), to determine the different growth levels of the five zones (F1-F5) described above. The GreenSeeker sensor calculates the NDVI values using the following formula: where R 774 is the reflectance of near-infrared light (774 nm) and R 656 is the reflectance of red light (656 nm).
In practice, the GS-NDVI values were acquired by positioning the sensor at a 0.8-m height above the canopy. These data were acquired simultaneously with the UAV data to minimize potential discrepancies caused by asynchronous data acquisition. A total of 60 randomly distributed GS-NDVI values were collected from each zone, resulting in a total of 300 measurements across the study area.

Data Processing and Analysis
The 60 GS-NDVI measurements acquired from each of the five zones were first averaged to create a single representative value for each zone. A map was then compiled of these five GS-NDVI values using the ArcMap software (ESRI, Redlands, CA, USA).
The UAV images were processed in three steps as shown in Figure 3. First, the original images were stitched and geometrically and radiometrically corrected using the Pix4D mapper software (Pix 4D, Inc., Lausanne, Switzerland). Using these images, the UAV-NDVI values were computed using the red and NIR bands of the RedEdge-M camera: where R 840 is the reflectance of near-infrared light (840 nm) and R 668 is the reflectance of red light (668 nm).
Remote Sens. 2020, 12, x FOR PEER REVIEW 6 of 15 GreenSeeker instrument has its internal light source and hence, the effects of the TOD are negligent in the ground data. The significance of the effects of the FA, TOD, SZA and growth level on the UAV-NDVIs was determined using the Standardized Regression Coefficients [31] analysis which reflects the relative importance of different variables. This analysis was conducted using the SPSS ® software package (IBM, NY, USA) to evaluate the order of influence.

The Rice Growth Levels Determined by the GS-NDVIs
In Figure 4, the mean GS-NDVI values, obtained through ground measurements, are shown. As can be seen in figure, there is a trend from the highest NDVI values in zone F5 toward the lowest NDVI values in zone F1. This trend is consistent with fertilizer use, shown in Figure 2, that increases Second, five shapefiles were created to represent each zone using ArcMap. Finally, the Zonal Statistics tool of ArcMap was used to calculate the mean, minimum, maximum and standard deviation (STDEV) of the pixel-scale UAV-NDVI data. These statistics were calculated separately for the five zones, represented by the five shapefiles created in step two. Of the collected values, the mean UAV-NDVI reflects the mean productivity and biomass, and the standard deviation indicates the spatial variability in productivity. The minimum and maximum UAV-NDVIs can provide information on pixel mixing. The count of pixels (COP) of a certain NDVI can give indicators of the NDVI distribution and is determined by the ground sampling distance (GSD), which in turn is influenced by the FA. Therefore, the COPs of each NDVI value were computed to analyze the effect of the FA on the UAV-NDVIs.
In addition to the aforementioned, a total of 48 SZAs were calculated using the NOAA Solar Calculator tool (https://www.esrl.noaa.gov) that was given the location of the study site and the TOD of each flight mission. From herein, identifiers M6-M17 are used to represent each TOD. For example, M6 stands for TOD from 6 a.m. to 7 a.m., and M13 represents the TOD from 1 p.m. to 2 p.m., and so on. It should be noted that this experiment was only conducted for the UAV data because the GreenSeeker instrument has its internal light source and hence, the effects of the TOD are negligent in the ground data. The significance of the effects of the FA, TOD, SZA and growth level on the UAV-NDVIs was determined using the Standardized Regression Coefficients [31] analysis which reflects the relative importance of different variables. This analysis was conducted using the SPSS ® software package (IBM, NY, USA) to evaluate the order of influence.

The Rice Growth Levels Determined by the GS-NDVIs
In Figure 4, the mean GS-NDVI values, obtained through ground measurements, are shown. As can be seen in figure, there is a trend from the highest NDVI values in zone F5 toward the lowest NDVI values in zone F1. This trend is consistent with fertilizer use, shown in Figure 2, that increases from zone F1 toward zone F5. Therefore, it is obvious that the differences in the NDVI values are due to the fertilizer strategy and the zones thus represent different growth levels of paddy rice.
Remote Sens. 2020, 12, x FOR PEER REVIEW 6 of 15 GreenSeeker instrument has its internal light source and hence, the effects of the TOD are negligent in the ground data. The significance of the effects of the FA, TOD, SZA and growth level on the UAV-NDVIs was determined using the Standardized Regression Coefficients [31] analysis which reflects the relative importance of different variables. This analysis was conducted using the SPSS ® software package (IBM, NY, USA) to evaluate the order of influence.

The Rice Growth Levels Determined by the GS-NDVIs
In Figure 4, the mean GS-NDVI values, obtained through ground measurements, are shown. As can be seen in figure, there is a trend from the highest NDVI values in zone F5 toward the lowest NDVI values in zone F1. This trend is consistent with fertilizer use, shown in Figure 2, that increases from zone F1 toward zone F5. Therefore, it is obvious that the differences in the NDVI values are due to the fertilizer strategy and the zones thus represent different growth levels of paddy rice.

Effects of the Flight Altitude (FA) on the NDVI Values
In this experiment, we analyzed the influence of the FA on the UAV-NDVIs. The resulted, shown in Figure 5, suggest five trends: (i) the mean and maximum UAV-NDVIs increase with decreasing FA (Figure 5a

Effects of the Flight Altitude (FA) on the NDVI Values
In this experiment, we analyzed the influence of the FA on the UAV-NDVIs. The resulted, shown in  In order to further analyze the COP differences of the UAV-NDVIs in the field-scale, the COP distribution curves of the UAV-NDVIs were analyzed under different FAs for fields B1-B5 ( Figure  6). It could be seen that, on the whole, the higher the FA, the lower the COP, the smaller the NDVI value of the peak COP and the smaller the NDVI value overall. For instance, in B1 (Figure 6a    In order to further analyze the COP differences of the UAV-NDVIs in the field-scale, the COP distribution curves of the UAV-NDVIs were analyzed under different FAs for fields B1-B5 ( Figure 6). It could be seen that, on the whole, the higher the FA, the lower the COP, the smaller the NDVI value of the peak COP and the smaller the NDVI value overall. For instance, in B1 (Figure 6a  In order to further analyze the COP differences of the UAV-NDVIs in the field-scale, the COP distribution curves of the UAV-NDVIs were analyzed under different FAs for fields B1-B5 ( Figure  6). It could be seen that, on the whole, the higher the FA, the lower the COP, the smaller the NDVI value of the peak COP and the smaller the NDVI value overall. For instance, in B1 (Figure 6a

Dependence of the UAV-NDVIs on the Solar Zenith Angle (SZA)
As can be seen in Figure 8, two general trends can be observed in the SZA values of the study: i) the mean, minimum and maximum UAV-NDVI values generally increased as a function of increasing SZA and ii) the standard deviation values decreased as a function of increasing SZA. The trends were similar across all zones (F1-F5) and FAs. As the lowest SZAs occur around the solar noon this means that the UAV-NDVI values were lower around noon than in the morning or in the afternoon. The reverse was true for the standard deviation values. These results are consistent with those discussed in Section 3.3, which is logical, since the SZAs depend on the TOD.

Dependence of the UAV-NDVIs on the Solar Zenith Angle (SZA)
As can be seen in Figure 8, two general trends can be observed in the SZA values of the study: (i) the mean, minimum and maximum UAV-NDVI values generally increased as a function of increasing SZA and (ii) the standard deviation values decreased as a function of increasing SZA. The trends were similar across all zones (F1-F5) and FAs. As the lowest SZAs occur around the solar noon this means that the UAV-NDVI values were lower around noon than in the morning or in the afternoon. The reverse was true for the standard deviation values. These results are consistent with those discussed in Section 3.3, which is logical, since the SZAs depend on the TOD.

Effects of the Growth Levels of Rice on the UAV-NDVIs
The differences between the mean UAV-NDVI values of the five zones followed those acquired using the ground data (F1: lowest, F5: highest, see Figure 4). In contrast, the standard deviation values had a downward trend from zone F1 to zone F5. With some slight variation, these differences could be detected irrespective of the FA, TOD and SZA (see Figures 5, 7 and 8). Thus, better growth levels induced higher NDVI values and lower standard deviation values. Furthermore, as could be seen in Figure 6, the curves of each field (B1-B5) consistently show a right shift (toward higher UAV-NDVI values) and a compressed value distribution with increasing growth level (B1: lowest, B5: highest). Hence, the better the growth level, the denser the UAV-NDVI distribution at the pixel scale.

Effects of the Growth Levels of Rice on the UAV-NDVIs
The differences between the mean UAV-NDVI values of the five zones followed those acquired using the ground data (F1: lowest, F5: highest, see Figure 4). In contrast, the standard deviation values had a downward trend from zone F1 to zone F5. With some slight variation, these differences could be detected irrespective of the FA, TOD and SZA (see Figures 5, 7 and 8). Thus, better growth levels induced higher NDVI values and lower standard deviation values. Furthermore, as could be seen in Figure 6, the curves of each field (B1-B5) consistently show a right shift (toward higher UAV-NDVI values) and a compressed value distribution with increasing growth level (B1: lowest, B5: highest). Hence, the better the growth level, the denser the UAV-NDVI distribution at the pixel scale.
To further study the trends associated with the growth levels, the averaged range of the mean UAV-NDVIs are shown individually for the FAs, TODs and SZAs in Table 3. As can be seen, the maximum averaged range of the FAs was 0.031 and that of TODs and SZAs was 0.200 in zone F1. In zone 5, the range was from 0.019 to 0.085. Based on the resulted, shown on Table 3, indicating that the better growth level could reduce the effect of FAs, TODs and SZAs.  To further study the trends associated with the growth levels, the averaged range of the mean UAV-NDVIs are shown individually for the FAs, TODs and SZAs in Table 3. As can be seen, the maximum averaged range of the FAs was 0.031 and that of TODs and SZAs was 0.200 in zone F1. In zone 5, the range was from 0.019 to 0.085. Based on the resulted, shown on Table 3, indicating that the better growth level could reduce the effect of FAs, TODs and SZAs. In order to further analyze the contribution of the above factors, the standardized regression coefficients (SRC) of the linear model are shown in Table 4. According to the positive relationship between the contribution and absolute values of the SRC, the descending contribution rank regarding the mean and STDEV UAV-NDVIs can be listed as: SZAs/TODs, growth levels and FAs. It can be found that the FAs affect the values less than the other parameters and the SZAs/TODs had a significant effect on the mean UAV-NDVIs and STDEV UAV-NDVIs. Table 4. Standardized regression coefficients (SRC) of FAs, SZAs/TODs and growth levels.

Discussion
Remote sensing has been widely used to assess crop growth in different environments. UAV-based imaging technology has the potential to provide high spatial resolution (up to centimeter-scale) maps for this end, providing instant feedback needed for crop management and decision making. In this context, increasing the understanding of how the flight parameters affect the NDVI values acquired using UAV-systems can help improve the image quality of the thus acquired data. In this study, we used a lightweight UAV equipped with a multispectral camera to collect NDVI values over a paddy rice under different flight parameters as well as different crop growth levels. The results suggest that the flight parameters and growth levels have a significant effect on the UAV-based data, thus highlighting the importance of careful flight planning. Nevertheless, it is essential to mention that the parameters and environmental factors (FA, TOD, SZA and growth level) discussed in this study do not fully account for the quality of the UAV-acquired data: Factors such as the flight speed can have a significant impact on the quality of the remote sensing data.

Sensitivity of UAV-NDVIs to the Flight Altitude (FA)
Our resulted, discussed in Section 3.2, suggest that the mean, minimum, maximum and standard deviation values of the UAV-NDVIs were highly related to the FA. This was particularly true for the COP of the peak UAV-NDVI values, shown in Figure 6. Additionally, as can be seen in the Appendix A  Table A1 and Figure A1, under the same growth level, the higher the FA, the smaller the individual UAV-NDVIs and the narrower the range of the UAV-NDVIs. Regarding the effects of the FA on the NDVI values, previous studies have reported different results. Stow et al. [32] evaluated the effects of the FA on the NDVI values through field experiments, finding that the relationship between them was inconclusive. In studies by Rasmussen et al. [33] and Yu et al. [34], no significant associations were found between the FA and the NDVI values. Easterday et al. [35] concluded that of the NDVI values acquired under different FAs (30 m, 60 m, 100 m and 120 m), the lowest FAs were better for detecting water deficit in plants. Similar to our study, a report by Mesas-Carrascosa et al. [36] concluded that even though the mean and maximum UAV-NDVIs were negatively related to the FA, the reverse was true for the minimum and STDEV UAV-NDVIs.
Notwithstanding the relationship between the NDVIs and the FA, the selection of the FAs has important implications from the point of view of the spatial resolution of the remote sensing imagery. In the field of imaging technology, a higher GSD means that each pixel represents a larger area. However, with the resulting lower spatial resolution of the image, the spectral information included in each pixel becomes more mixed and hence, detecting objects and phenomena can become more difficult. At high FAs associated with low spatial resolutions, some pixels could be mixed, and the average NDVIs may become diluted [37].
Looking at the relationship between the FA and the number of pixels of this study, Figure A1 shows four subsets of the UAV-NDVI maps of the B2 field at different FAs at TOD M11. In this analysis, four subsets were selected to represent the FAs of 40, 60, 80 and 100 m. These FAs resulted in GSDs values of 2.88, 4.32, 5.81 and 7.22 cm×pixel −1 that had NDVI values between 0.5 and 0.86. As can be seen in Figure A1, increasing the GSD resulted in fewer pixels in the same sampled area and narrower COP distributions. Under normal circumstances, the higher the FA of the UAV platform, the larger the area represented by a single pixel, and thus, the lesser the ability of the UAV system to detect small features. In the case of a paddy rice field, despite their commonly high canopy densities, there are still many canopy gaps. These gaps have other elements such as soil and water, which, depending on the FA, can lead to different degrees of spectral mixing. This, in turn, can lead to overall smaller mean NDVI and standard deviation values as suggested by Figure 5a,b, respectively. Therefore, as far as the NDVI monitoring of rice is concerned, the average NDVI was inversely proportional to the FAs. This further supports the conclusion that the FA was an important factor affecting the image quality for the crops. In general, it can be concluded that FAs do impact the mean UAV-NDVIs, but it is not conclusive that the mean UAV-NDVIs will keep decreasing with increasing FA. A possible extreme flight altitude needs further exploration.
Another noticeable factor regarding the effect of FA on vegetation indexes was atmospheric effects. Yu et al. suggested normalized excess green index (ExG) and the normalized green-red difference index (NGRDI) were relatively susceptible to FAs than the NDVI values (FAs: 10 m, 30 m, 50 m and 100 m) considering of the atmospheric effects [34]. The atmospheric attenuation was small in the low-altitude atmospheres and the effect of path radiance in their model was negligible for red and infrared wavelengths. Thus, we did not make an analysis of the atmospheric effects upon NDVIs when the FAs were under 100 m. In fact, the atmospheric effects were more significant at higher FAs (manned aircrafts and satellites) [38].

Influence of the time of Day (TOD) and Solar Zenith Angle (SZA) on the UAV-NDVIs
Although the SZA values vary as a function of the time of day, date and location [39], in this study, SZAs and TODs were discussed independently, because there is the spatiotemporal difference. To the knowledge of the authors, our study is one of the first ones to report and discuss the effects of the TODs and SZAs on UAV-based NDVI values continuously from sunrise to sunset. Based on our results, the mean, minimum, maximum and standard deviation of the UAV-NDVIs are all highly related to the TODs and SZAs (see Sections 3.3 and 3.4). In previous studies, Ishihara et al. [28] suggests that the response of vegetation indices to the SZA was evident under a clear sky. Rahman [29] reported NDVIs from ground-based observations that decreased with decreasing SZA at a pasture site, but did not provide explanations to this finding. Their conclusions were directly based on the average assessment of ground-NDVIs, and hence, not remote sensing data. Although the results in this study show a similar phenomenon, we have further analyzed the UAV-NDVI distribution on a pixel-scale (including mean, minimum, maximum, standard deviation and COP of the UAV-NDVIs), an experimental setup that can reveal the relationships between the UAV-NDVIs and the flight parameters (SZAs/TODs) more clearly. Ishihara et al. suggested that surveys involving NDVI measurements be performed at a SZA of 60 • to be effective for the accurate assessment of the canopy structure and function. However, our results, shown in Figure 8, suggest no obvious differences between the mean UAV-NDVIs and standard deviation values of the UAV-NDVIs with different growth levels at a SZA of 60 • . Therefore, we do not follow the recommendation by Ishihara et al. to acquire UAV-NDVI values at a SZA of 60 • . In practice, the TODs and SZAs were the direct reasons for changing the angle between the view direction and incident light. An important phenomenon in this context is the widely discussed hot spot effect, also called the hot spot directional signature in the backscattering direction [40]. In recent studies, this phenomenon was also called the bidirectional reflectance distribution function (BRDF) effects. Generally, the radiance of the reflecting medium will decrease with an increasing angle between the view direction and incident light because of the decreased probability of seeing illuminated particles [41]. Therefore, because of the vertical downward view of the RedEdge camera, the strongest radiance and relatively complete reflectivity from the rice canopy was collected when the SZAs were relatively small. By definition, crop reflectance is the ratio of the amount of light leaving the canopy to the amount of incoming light. We suggest the reflectivity was relatively accurate at a smaller SZAs, in turn, resulting in reliable NDVI values.

Effect of the Growth Level on the UAV-NDVIs
It can also be seen in the data that the application rate of fertilizers can have a significant impact on the NDVI values ( Figure 4). In the early stages of rice growth or in the case of weak growth that results from inadequate fertilization, the canopy is distributed sparsely and there are more gaps with water between seedlings [42]. These potential gaps would influence the NDVI values. For instance, looking at the image acquired at M11 (11 a.m.-12 p.m.) at 40 m ( Figure A1), the image has a relatively high GSD of 2.88 cm at a FA of 40 m. Due to the resulting high spatial resolution, the image contains many details, including those related to non-foliar areas, and thus, small NDVI values that result from pixels with soil or water. Similarly, areas with lesser amounts of fertilizers and thus, less vigorous growth levels will have more gaps with smaller NDVI values and higher standard deviation values than result from spatial heterogeneity (vegetation, water and soil). The higher the flight altitude, the lower the spatial resolution and the lesser the amount of details in the image. Therefore, in the UAV image acquired at 100 m, the standard deviation is lower and the range of NDVI values is narrower than in the image that was acquired at 40 m (Figures 7 and A1). On the other hand, the main limitation in our growth level analysis was the lack of observations to cover entire growth stages. Such an analysis would have provided a more comprehensive quantitative framework to interpret the interrelationships between the growth levels and the flight parameters. In fact, in the early stages of rice growth, due to the sparse rice canopy, aerial remote sensing will collect a large area of ground water surface (mirror reflectance), which is a growth period that is not suitable for UAV remote sensing monitoring. Hence, the available growth stages for reliable UAV remote sensing are limited.

Conclusions
In this study, a paddy rice field was used as a test area to assess the effects of FA, TOD, SZA and growth level on the UAV-acquired NDVI values. Based on the SRC values of the linear regression resulted, all of the former parameters had a significant effect on the UAV-NDVIs. More specifically, our results suggest that (1) the SZA/TOD had the largest impact on the mean, minimum, maximum and STDEV of UAV-NDVI values, followed by the growth level and the FA; (2) the mean, maximum and STDEV of UAV-NDVIs were inversely proportional to the FA (≤100 m) and the minimum of UAV-NDVIs decreased with increasing FA; (3) the mean, minimum and maximum of UAV-NDVIs were proportional to the SZA, but it was contrary with the STDEV of UAV-NDVIs; (4) the effect of TOD on UAV-NDVIs could refer to the SZA;(5) Furthermore, according to our results, the UAV-NDVIs close to the smaller SZAs show the highest signal-to-noise ratios which we infer to reflect the most realistic growth status values. We expect that our resulted and the recommendations could provide a reference to the operating parameters of UAVs in the context of precision agriculture. like to thank the reviewer for their helpful comments. We thank AJE (https://secure.aje.com/) for its linguistic assistance during the preparation of this manuscript.

Conflicts of Interest:
The authors declare no conflicts of interest. Appendix A Figure A1. Subset of the B2 NDVI map at 40, 60, 80, 100 m with TOD of M11 (11 am-12 pm).