Potential of UAS-Based Remote Sensing for Estimating Tree Water Status and Yield in Sweet Cherry Trees

The present work aims to assess the usefulness of five vegetation indices (VI) derived from multispectral UAS imagery to capture the effects of deficit irrigation on the canopy structure of sweet cherry trees (Prunus avium L.) in southeastern Spain. Three irrigation treatments were assayed, a control treatment and two regulated deficit irrigation treatments. Four airborne flights were carried out during two consecutive seasons; to compare the results of the remote sensing VI, the conventional and continuous water status indicators commonly used to manage sweet cherry tree irrigation were measured, including midday stem water potential (Ψs) and maximum daily shrinkage (MDS). Simple regression between individual VIs and Ψs or MDS found stronger relationships in postharvest than in preharvest. Thus, the normalized difference vegetation index (NDVI), resulted in the strongest relationship with Ψs (r2 = 0.67) and MDS (r2 = 0.45), followed by the normalized difference red edge (NDRE). The sensitivity analysis identified the optimal soil adjusted vegetation index (OSAVI) as the VI with the highest coefficient of variation in postharvest and the difference vegetation index (DVI) in preharvest. A new index is proposed, the transformed red range vegetation index (TRRVI), which was the only VI able to statistically identify a slight water deficit applied in preharvest. The combination of the VIs studied was used in two machine learning models, decision tree and artificial neural networks, to estimate the extra labor needed for harvesting and the sweet cherry yield.


Introduction
Precision agriculture cannot be conceived without the recent evolution of technology. Technological tools such as soil-plant-atmosphere sensors and multispectral images obtained from unmanned aerial systems (UASs) are widely used to improve crop management. In recent years, several studies carried out in extensive crops and, to a lesser extend in woody crops, have reported promising results for UAS-based multispectral imagery for the remote sensing of crops and soil status as well as for detecting pests, diseases, and water and nutritional stress [1][2][3][4].
Multispectral aerial imagery is often used to calculate vegetation indices (VI) as combinations of the proportion of light reflected by leaves in the visible, red-edge, and near-infrared parts of the electromagnetic spectrum [5]. Among the most commonly used vegetation indices are the normalized tree and artificial neural networks methods, can be used as classification or regression models for discrete/continuous estimation values such as the quantity of irrigation water needed for the following week [31] or for fruit yield prediction/classification as we propose in this work.
The overall aim of this paper was to ascertain the strengths and weaknesses of remote sensing indices obtained by multispectral imagery (UAS-based data) to detect the changes in canopy structure caused by water stress and to estimate the yield of adult sweet cherry trees. To achieve this goal, the response of remote structural indices (NDVI, OSAVI, DVI, OSAVI, and TRRVI spectroscopy) to two RDI treatments and a full irrigation treatment is compared with the response obtained by a conventional indicator (midday stem water potential) and a continuous indicator (maximum daily branch shrinkage).

Study Site
The study was conducted in a 0.50 ha experimental field of 17-year-old sweet cherry trees (Prunus avium L.) combination 'Prime Giant'/SL64, in the commercial orchard "Finca Toli" (38 • 8 N; 1 • 22 W, Elev. 680 m) located in Jumilla, Spain for two consecutive seasons, 2017 and 2018. Trees were trained to a Spanish bush (vase) system and spaced 5 m between rows and 3 m apart within the row. The soil is moderately stony and has a sandy loam texture with 67.5% sand, 17.5% silt, and 15% clay. The irrigation system consisted of a single drip line per tree row and three pressure-compensated emitters (4 L h −1 ) per tree. The irrigation water was from a well, and had an average electrical conductivity (EC 25 • C ) of 0.8 dS m −1 . Full bloom was in April and the harvest in June (2 June 2017 and 14 and 19 June 2018). The annual pruning was carried out approximately 60 d after harvest in accordance with the local growers' normal practice.
Daily climatic data, including crop reference evapotranspiration (ET 0 ), air temperature, air relative humidity, and rainfall, were recorded by an automatic weather station located 13 km from the experimental orchard owned by the Spanish Consultancy Service for Irrigation Users (SIAR; http: //crea.uclm.es/siar/datmeteo/). From the temperature and humidity data, the vapor pressure deficit (VPD) was calculated according to Allen et al. [32].

Experimental Design
A completely randomized block experimental design with three irrigation treatments and four replicates (of seven trees each one) per treatment was established: a control treatment (CTL) and two regulated deficit irrigation treatments (RDI-1 and RDI-2). CTL was irrigated under non-limiting soil water conditions at 110% crop evapotranspiration (ET c ), while RDI-1 was irrigated at 90% ET c during preharvest (slight water stress), 100% ET c during flower differentiation (15-20 days after harvest, critical period) and 65% ET c during postharvest (non-critical period) and RDI-2 was irrigated at 100% ET c during preharvest and flower differentiation and 55% ET c during postharvest (severe water stress). The irrigation dose was calculated using the formula ET c = ET 0 * K c * K r , where K c is the crop-specific coefficient for sweet cherry trees reported by Marsal [33] and K r is a localization factor related to the percentage of ground covered by the crop [34].

Airborne Imagery
On 29 May and 1 August in 2017 and on 13 June and 21 August in 2018, four flight campaigns were carried out at noon. Aerial images were collected on sunny days with no noticeable wind by a multi-spectral camera (mod. Sequoia, Parrot Co. Ltd., Paris, France) of 1280 × 960 pixels equipped with GPS which was mounted on an unmannered aerial system (UAS mod. Phantom 3 Advance, DJI, Shenzhen, China). The multi-spectral camera captures images in four different bands of visible light: green, G: 550 nm center wavelength, 40 nm bandwidth; red, R: 660 nm, 40 nm; red edge, RE: 735 nm, 10 nm; near-infrared, NIR: 790 nm, 40 nm. The imagery was geo-processed to generate an orthomosaic Remote Sens. 2020, 12, 2359 4 of 18 with the Pix4Dmapper Pro software (Pix4D SA, Lausanne, Switzerland). The flight altitude for the multispectral images was set at 50 m above ground level, which provided images with a ground resolution of 6.2 cm per pixel. From each image, trees were manually identified and their canopies were individualized to arrange them into circles of 2.5 m diameter, ensuring that only pure-vegetation tree crown pixels were chosen. Remote sensing indices were calculated with the open-software system Quantum Geographical Information System (QGIS 3.12 Bucuresti), using data from multiple reflectance bands. Remote vegetation indices NDVI, OSAVI, DVI, and NDRE were calculated using standard formulas from each processed UAS image mosaic [6][7][8][9]: In addition, a new index is proposed in this work: NDVI, OSAVI, and DVI were computed in 2017 and 2018 (four flights); however, NDRE and TRRVI were only computed in 2018 (two flights).

Field Measurements
On the same day and at the same time as the UAS flight, midday stem water potential (Ψ s , a conventional water status indicator) was measured with a Scholander pressure chamber (Model 3000, Soil Moisture Equipment, Santa Barbara, CA, USA) in six trees per treatment following the methodology proposed by McCutchan and Shackel [35]. The mature and healthy leaves selected were enclosed in black polythene bags and covered with aluminum foil for 2 h before the measurement. The branch diameter fluctuations of the same trees were recorded by dendrometers (LVDT sensors, model DF ± 2.5 mm, accuracy ± 10 µm, Solartron Metrology, Bognor Regis, UK). Each dendrometer was placed on the main tree branch away from direct sunlight. To prevent thermal expansion, the sensors were installed on holders made of invar and aluminum. From the branch diameter fluctuations, the maximum daily shrinkage (MDS) (continuous water status indicator) was calculated as the daily difference between the maximum and minimum branch diameters.
Vegetative growth was measured as the mass of pruning wood from each tree, which was individually weighed in the field using electronic scales (Scaltec Instruments GmbH, mod SSH91, Heiligenstadt, Germany). Similarly, in 2018 80 trees were selected at harvest and the individual yield of each tree was weighed to determine yield per tree.

Machine Learning Methodologies for Estimating Tree Crop Yield
Decision tree (DT) and artificial neural network (ANN) algorithms were used for harvesting prediction These types of algorithms have been widely used in the application of artificial intelligence techniques to the field of agriculture [36][37][38][39]. DT was performed using random forest as a classification model between two classes, normal labor predicted for harvesting (below a defined threshold of 40 kg tree −1 ) and extra labor (above the threshold). ANN was implemented to provide sweet cherry tree yield estimation as a regression model, using the RStudio package (RStudio Inc., Boston, MA, USA). The software package used in the ANN analysis was 'neuralnet'.
Vegetation indices (VI) data were calculated from the UAS images taken during the preharvest flight in 2018 in 80 trees (predictor variables). Additional manual quantification data of the yield per tree was calculated for dataset development (response variable). The dataset was divided into Remote Sens. 2020, 12, 2359 5 of 18 57 trees (71%) for training and 23 trees (29%) for testing, performing a cross-validation process based on the grouping of training and test sets following 10-fold cross-validation, which is considered useful when few data are available. This method increases the number of validations of the model for the same amount of data by crossing different subsamples of the same data [40]. The data were split into 10 equal parts (subsamples) used for the 10-fold cross validation. Each time, a different part was used as test set and the remaining nine parts were used as training set.
Performance values, such as precision, recall, and true positives, were calculated from the error matrix (confusion matrix).
Data from all the variables were scaled between 0 and 1 to improve convergence. This process is reversed in the predicted output of the model. Five features were used as input layer and two hidden layers with 3 and 4 neurons besides a bias neuron in each hidden layer. The model used a back-propagation algorithm which computed the gradients by iterating backwards. Root mean square error (RMSE) and the coefficient of determination were evaluated for gradient thresholding. Furthermore, for each pattern (actual and predicted) minimum and maximum values, standard deviation, coefficient of variation, and the normality tests Skewness and Kurtosis were calculated to ascertain any differences between both data distributions.

Statistical Analysis
The sensitivity analysis of the water stress indicators was calculated according to Goldhamer and Fereres [41] where sensitivity (S) is the result of the division of the signal intensity (SI) by the coefficient of variation (CV), and is always greater than 0. Analysis of variance (ANOVA), analysis of covariance (ANCOVA) and post-hoc test (Duncan's multiple ranges) with a significance level of p < 0.05 were performed using the statistical software package IBM SPSS Statistic 24 (IBM Corp., Armonk, NY, USA) and Statgraphics centurion XVI (StatPoint Technologies Inc., The Plains, VA, USA) in order to assess the effects of the imposed irrigation treatments on the response of conventional and continuous water status indicators, remote vegetation indices, and vegetative growth. The regression analysis was graphed and calculated with Microsoft Excel (Windows 10 Home, Microsoft Corp., Redmond, WA, USA) and SigmaPlot 12.5 (Systat Software Inc., San Jose, CA, USA).

Environmental Conditions and Irrigation Water Applied
The weekly average reference evapotranspiration (ET 0 ) and vapor pressure deficit (VPD) and accumulated weekly precipitation recorded during the 2017 and 2018 seasons are summarized in Figure 1. Both VPD and ET 0 reflected the typical trend for Mediterranean climates both years of the study, with warm summers and mild winters.
Accumulated seasonal reference evapotranspiration from March to November reached close to 1000 mm both years of study, while the total rainfall recorded for the same period was 162 and 288 mm, for 2017 and 2018, respectively. In neither year did rainfall meet crop evapotranspiration levels, and irrigation was required; moreover, the rainfall pattern differed markedly between both years, with 124 mm falling during June-September (the period with the highest VPD values) in 2018, and 51 mm in 2017. The lower level of rainfall in 2017 than in 2018 led to a higher volume of irrigation water being applied, the volume of irrigation water decreasing by 10% in all irrigation treatments from 2017 to 2018 as a result. RDI-2 was the irrigation treatment needing the lowest amount of water in both years: 7% and 36% less water than RDI-1 and CTL, respectively. However, during preharvest, RDI-1 was the only irrigation treatment that applied deficit irrigation.

Environmental Conditions and Irrigation Water Applied
The weekly average reference evapotranspiration (ET0) and vapor pressure deficit (VPD) and accumulated weekly precipitation recorded during the 2017 and 2018 seasons are summarized in Figure 1. Both VPD and ET0 reflected the typical trend for Mediterranean climates both years of the study, with warm summers and mild winters.

Conventional and Continuous Plant Water Status Indicators and Remote Vegetation Indices
The plant water status indicators and VI provided different responses to the irrigation strategies assayed. Overall, conventional indicator (Ψ s ), continuous indicator (MDS) and remote VI (NDVI, OSAVI, DVI, NDRE, TRRVI) identified the water deficit applied to sweet cherry trees in the RDI-1 and RDI-2 treatments during postharvest compared with the trees irrigated to satisfy their water requirements (CTL). Figure 2 shows the NDVI maps taken from the orchard on two selected days in 2017, one day during preharvest and another day during postharvest.
Remote Sens. 2020, 12, x FOR PEER REVIEW 6 of 18 Accumulated seasonal reference evapotranspiration from March to November reached close to 1000 mm both years of study, while the total rainfall recorded for the same period was 162 and 288 mm, for 2017 and 2018, respectively. In neither year did rainfall meet crop evapotranspiration levels, and irrigation was required; moreover, the rainfall pattern differed markedly between both years, with 124 mm falling during June-September (the period with the highest VPD values) in 2018, and 51 mm in 2017. The lower level of rainfall in 2017 than in 2018 led to a higher volume of irrigation water being applied, the volume of irrigation water decreasing by 10% in all irrigation treatments from 2017 to 2018 as a result. RDI-2 was the irrigation treatment needing the lowest amount of water in both years: 7% and 36% less water than RDI-1 and CTL, respectively. However, during preharvest, RDI-1 was the only irrigation treatment that applied deficit irrigation.

Conventional and Continuous Plant Water Status Indicators and Remote Vegetation Indices
The plant water status indicators and VI provided different responses to the irrigation strategies assayed. Overall, conventional indicator (Ψs), continuous indicator (MDS) and remote VI (NDVI, OSAVI, DVI, NDRE, TRRVI) identified the water deficit applied to sweet cherry trees in the RDI-1 and RDI-2 treatments during postharvest compared with the trees irrigated to satisfy their water requirements (CTL). Figure 2 shows the NDVI maps taken from the orchard on two selected days in 2017, one day during preharvest and another day during postharvest. Neither the OSAVI, nor the DVI were sufficiently sensitive to detect differences between CTL and RDI-1 any year of study, although DVI was able to identify significant differences between CTL and RDI-2 trees both years of study ( Figure 3). In the postharvest of 2017, Ψs and NDVI were able to differentiate not only between CTL and RDI treatments, but also between RDI-1 and RDI-2 ( Figure  3B,F). In 2017 postharvest, trees of all treatments had lower water status than in 2018, Ψs reaching mean values of −1.04 and −1.33 MPa for RDI-1 and RDI-2, respectively, compared with the -0.74 and -0.81 MPa reached for the same treatments in 2018. In the postharvest of 2018, Ψs was the only plant water status indicator that could significantly distinguish between RDI-1 and RDI-2.
On the other hand, in preharvest the response of the plant water status indicators and VI differed from that exhibited in postharvest. During the preharvest, only the sweet cherry trees of RDI-1 were submitted to water stress. The slight water stress applied to RDI-1 trees during preharvest led to significant differences in conventional (Ψs) and continuous (MDS) water stress indicators ( Figure  3A,C). Among the remote indicators, only TRRVI was able to recognize the slight deficit applied. Neither NDVI, nor NDRE showed significant differences between RDI-1 and CTL trees during Neither the OSAVI, nor the DVI were sufficiently sensitive to detect differences between CTL and RDI-1 any year of study, although DVI was able to identify significant differences between CTL and RDI-2 trees both years of study ( Figure 3). In the postharvest of 2017, Ψ s and NDVI were able to differentiate not only between CTL and RDI treatments, but also between RDI-1 and RDI-2 ( Figure 3B preharvest ( Figure 3E,K), although both indicators were able to identify different levels of water stress, mid and strong, in the same trees during postharvest.  On the other hand, in preharvest the response of the plant water status indicators and VI differed from that exhibited in postharvest. During the preharvest, only the sweet cherry trees of RDI-1 were submitted to water stress. The slight water stress applied to RDI-1 trees during preharvest led to significant differences in conventional (Ψ s ) and continuous (MDS) water stress indicators ( Figure 3A,C). Among the remote indicators, only TRRVI was able to recognize the slight deficit applied. Neither NDVI, nor NDRE showed significant differences between RDI-1 and CTL trees during preharvest ( Figure 3E,K), although both indicators were able to identify different levels of water stress, mid and strong, in the same trees during postharvest.
In order to compare the plant water status indicators and the VI, a sensitivity analysis was carried out ( Table 1). The analysis showed that MDS was the plant water stress indicator with the highest SI (1.14 and 1.68, mean values for RDI treatments during preharvest and postharvest, respectively) followed by Ψ s and NDRE. However, MDS was also the indicator that exhibited the highest CV (0.12-0.18), followed by OSAVI, which reduced its S as indicator. As regards the remote VI, and regardless of the irrigation treatment evaluated, NDRE showed the lowest variability and the highest sensitivity in preharvest. Furthermore, NDRE, OSAVI, and TRRVI resulted in higher signal intensity and sensitivity in preharvest than in postharvest (Table 1). By contrast, NDVI and DVI provided a higher S in postharvest than in preharvest mainly because of the low CV exhibited in the postharvest flights compared with that shown in preharvest (Table 1). NDVI was the most sensitive remote indicator to both irrigation treatments in postharvest, followed by NDRE, TRRVI, OSAVI, and DVI showing similar sensitivity to Ψ s and being twice as sensitive as MDS.

Relationship between Plant Water Status Indicators and Vegetation Indices
The relationship between water status indicators and VI was examined by linear correlation analysis. The response of NDVI, OSAVI, DVI, and NDRE at two phenological moments (preharvest and postharvest) was compared in the CTL trees ( Figure 4).
VI were more strongly related to Ψ s in postharvest than in preharvest, when no relationship was found ( Figure 4). VI were not able to clearly identify slight differences within the same year in CTL trees, neither 2017, nor 2018; however, VI identified that CTL trees in 2017 showed lower values of Ψ s than in 2018. NDVI was the index which best fitted with Ψ s in 2017 (r 2 = 0.55). In 2018, VI were less related to Ψ s, and NDVI and NDRE resulted in similar, and poor, linear relationships to Ψ s (r 2 = 0.33 and 0.27, respectively). On the other hand, OSAVI and DVI showed no relationship at both times. In the postharvest baseline of NDVI (Figure 4), it can be observed that there is a clear difference between the measurements taken in 2017 and 2018, and when both years were considered the relationship between Ψ s and NDVI became stronger (r 2 = 0.78). Differences in environmental water demand and  Similarly, when data of all treatments were pooled ( Figure 5), NDVI in preharvest still showed no relationship with Ψs (r 2 = 0.07), and, while the relationship found was also strong (r 2 = 0.67) in postharvest, the equation obtained was statistically different (according to the ANCOVA analysis) from that obtained only for the CTL treatment (Figures 4 and 5A,B). When both regression lines were compared with the slopes obtained for the relationship between NDVI and Ψs, significant differences were observed (p-value = 0.0017). In the case of NDRE vs. Ψs, the ANCOVA analysis showed no difference between the slopes or intercepts of the regression lines, and for the measurements taken in postharvest, the relationship increased from r 2 = 0.29 to 0.54 (Figures 4 and 5E). The relationship between DVI and Ψs in postharvest remained poor despite the addition of the RDI treatments ( Figure  5D); however, OSAVI, which resulted poorly related with Ψs in postharvest when only CTL trees were considered, slightly improved the relationship with Ψs when all treatments were included, from r 2 = 0.30 to 0.41 ( Figure 5C). As regards the postharvest flight, NDVI showed a relationship with Ψs that was slightly stronger than that of NDRE and higher than those of OSAVI, TRRVI, and DVI which showed similar results (r 2 = 0.4 Figure 5C,D,F). The relationship between the remote indices and the Similarly, when data of all treatments were pooled ( Figure 5), NDVI in preharvest still showed no relationship with Ψ s (r 2 = 0.07), and, while the relationship found was also strong (r 2 = 0.67) in postharvest, the equation obtained was statistically different (according to the ANCOVA analysis) from that obtained only for the CTL treatment (Figures 4 and 5A,B). When both regression lines were compared with the slopes obtained for the relationship between NDVI and Ψ s , significant differences were observed (p-value = 0.0017). In the case of NDRE vs. Ψ s , the ANCOVA analysis showed no difference between the slopes or intercepts of the regression lines, and for the measurements taken in postharvest, the relationship increased from r 2 = 0.29 to 0.54 (Figures 4 and 5E). The relationship between DVI and Ψ s in postharvest remained poor despite the addition of the RDI treatments ( Figure 5D); however, OSAVI, which resulted poorly related with Ψ s in postharvest when only CTL trees were considered, slightly improved the relationship with Ψ s when all treatments were included, from r 2 = 0.30 to 0.41 ( Figure 5C). As regards the postharvest flight, NDVI showed a relationship with Ψ s that was slightly stronger than that of NDRE and higher than those of OSAVI, TRRVI, and DVI which showed similar results (r 2 = 0.4 Figure 5C,D,F). The relationship between the remote indices and the continuous plant-based water stress indicator, MDS, was lower than with Ψ s , while NDRE showed a similar coefficient of determination to NDVI ( Figure 5G,H).
Remote Sens. 2020, 12, x FOR PEER REVIEW 10 of 18 continuous plant-based water stress indicator, MDS, was lower than with Ψs, while NDRE showed a similar coefficient of determination to NDVI ( Figure 5G,H).   Figure 6 shows the relationship between the VI and the conventional and continuous indicators Ψ s and MDS. It was observed that, in general, VI were more strongly related to Ψ s than to MDS in both 2017 and 2018. Furthermore, it could be seen that the relationships increased in 2017 and decreased in 2018, which would be due to the greater range of values measured in 2017, as RDI trees experienced greater stress in 2017 than in 2018. According to the coefficients of determination obtained, the indices that were most strongly related to Ψ s and MDS were NDVI and NDRE (r 2 0.50) followed by TRRVI and OSAVI (r 2 0.30), and lastly by DVI (r 2 0.20) and NIR (r 2 0.15).
Remote Sens. 2020, 12, x FOR PEER REVIEW 11 of 18 red range vegetation index) (F). Linear relationships between maximum daily branch shrinkage and NDVI (G) and NDRE (H). Figure 6 shows the relationship between the VI and the conventional and continuous indicators Ψs and MDS. It was observed that, in general, VI were more strongly related to Ψs than to MDS in both 2017 and 2018. Furthermore, it could be seen that the relationships increased in 2017 and decreased in 2018, which would be due to the greater range of values measured in 2017, as RDI trees experienced greater stress in 2017 than in 2018. According to the coefficients of determination obtained, the indices that were most strongly related to Ψs and MDS were NDVI and NDRE (r 2 ≅ 0.50) followed by TRRVI and OSAVI (r 2 ≅ 0.30), and lastly by DVI (r 2 ≅ 0.20) and NIR (r 2 ≅ 0.15).

Assessment of Tree Vigor using Remote Indices
Increases in water stress in sweet cherry trees are closely linked to lower vegetative growth, so, pruning wood mass produced by each tree might be a good indicator of the vegetative growth of the trees and the tree vigor. In our case, the pruning wood mass for both years resulted in a mean value of 9.39 ± 0.78 kg tree -1 , with significant differences among treatments (p-value = 0.011), years (p-value < 0.001), and their interaction (p-value < 0.021). Trees under RDI treatments exhibited significantly lower pruning wood mass than CTL trees both years of the study, and trees showed more vegetative growth and consequently more pruning wood mass in 2018 than in 2017. Thus, the values obtained for each tree by the remote VI during postharvest were linearly related with the pruning wood mass of the same tree ( Figure 7). Since the pruning wood mass had an asymmetric distribution, the relationships were calculated with its logarithm (log (PW) as the logarithm of the pruning wood mass). Linear regression analysis pointed to a positive and strong relationship between NDVI and log (PW), r 2 = 0.75, Figure 7. In the

Assessment of Tree Vigor Using Remote Indices
Increases in water stress in sweet cherry trees are closely linked to lower vegetative growth, so, pruning wood mass produced by each tree might be a good indicator of the vegetative growth of the trees and the tree vigor. In our case, the pruning wood mass for both years resulted in a mean value of 9.39 ± 0.78 kg tree −1 , with significant differences among treatments (p-value = 0.011), years (p-value < 0.001), and their interaction (p-value < 0.021). Trees under RDI treatments exhibited significantly lower pruning wood mass than CTL trees both years of the study, and trees showed more vegetative growth and consequently more pruning wood mass in 2018 than in 2017. Thus, the values obtained for each tree by the remote VI during postharvest were linearly related with the pruning wood mass of the same tree (Figure 7). Remote Sens. 2020, 12, x FOR PEER REVIEW 11 of 18 red range vegetation index) (F). Linear relationships between maximum daily branch shrinkage and NDVI (G) and NDRE (H). Figure 6 shows the relationship between the VI and the conventional and continuous indicators Ψs and MDS. It was observed that, in general, VI were more strongly related to Ψs than to MDS in both 2017 and 2018. Furthermore, it could be seen that the relationships increased in 2017 and decreased in 2018, which would be due to the greater range of values measured in 2017, as RDI trees experienced greater stress in 2017 than in 2018. According to the coefficients of determination obtained, the indices that were most strongly related to Ψs and MDS were NDVI and NDRE (r 2 ≅ 0.50) followed by TRRVI and OSAVI (r 2 ≅ 0.30), and lastly by DVI (r 2 ≅ 0.20) and NIR (r 2 ≅ 0.15).

Assessment of Tree Vigor using Remote Indices
Increases in water stress in sweet cherry trees are closely linked to lower vegetative growth, so, pruning wood mass produced by each tree might be a good indicator of the vegetative growth of the trees and the tree vigor. In our case, the pruning wood mass for both years resulted in a mean value of 9.39 ± 0.78 kg tree -1 , with significant differences among treatments (p-value = 0.011), years (p-value < 0.001), and their interaction (p-value < 0.021). Trees under RDI treatments exhibited significantly lower pruning wood mass than CTL trees both years of the study, and trees showed more vegetative growth and consequently more pruning wood mass in 2018 than in 2017. Thus, the values obtained for each tree by the remote VI during postharvest were linearly related with the pruning wood mass of the same tree (Figure 7). Since the pruning wood mass had an asymmetric distribution, the relationships were calculated with its logarithm (log (PW) as the logarithm of the pruning wood mass). Linear regression analysis pointed to a positive and strong relationship between NDVI and log (PW), r 2 = 0.75, Figure 7. In the Since the pruning wood mass had an asymmetric distribution, the relationships were calculated with its logarithm (log (PW) as the logarithm of the pruning wood mass). Linear regression analysis pointed to a positive and strong relationship between NDVI and log (PW), r 2 = 0.75, Figure 7. In the same vein, OSAVI and DVI also showed a positive trend between measurements and vegetative growth, although it was weaker than that obtained with NDVI. For their part, TRRVI and NDRE measurements were poorly related with log (PW) despite the excellent results obtained by NDRE as water stress detector (Figure 4).

Tree Yield Estimation
There was no significant effect of the irrigation treatment on tree yield. The mean yield value at harvest was 42 kg tree −1 in 2018, a year that was considered a medium-high cropping year in 'Prime Giant'/SL64 sweet cherry trees. Thus, from the data of all trees, a decision tree (DT) model and an artificial neural network (ANN) approach were proposed in order to classify and estimate tree yield.
Both DT and ANN models used the plant-based remote indicators as predictor variables and the amount of cherries produced as the response variable. The mean values of the remote indicators in the 80 trees considered are shown in Table 2. With the aim of designing a DT to estimate whether the tree crop load is going to be higher than usual and so more personnel will be required for harvesting, a threshold value was set at 40 kg tree −1 . High cropping years necessitate more personnel to hand pick the cherries early in their commercial ripeness state and when market prices are higher. The DT aimed to identify two categories: above and below 40 kg tree −1 . Figure 8 shows the resulting DT, which identifies three constraints for the three remote indicators NIR, TRRVI, and NDVI. When the decision tree model was evaluated, it showed a training error of 0.228 and correctly identified 20 trees out of 23 (testing sample). When the confusion matrix was obtained for the testing sample, the precision was 0.870, the recall 0.778, and the true positive rate (specificity) 0.929.
Remote Sens. 2020, 12, x FOR PEER REVIEW 12 of 18 same vein, OSAVI and DVI also showed a positive trend between measurements and vegetative growth, although it was weaker than that obtained with NDVI. For their part, TRRVI and NDRE measurements were poorly related with log (PW) despite the excellent results obtained by NDRE as water stress detector (Figure 4).

Tree Yield Estimation
There was no significant effect of the irrigation treatment on tree yield. The mean yield value at harvest was 42 kg tree −1 in 2018, a year that was considered a medium-high cropping year in 'Prime Giant'/SL64 sweet cherry trees. Thus, from the data of all trees, a decision tree (DT) model and an artificial neural network (ANN) approach were proposed in order to classify and estimate tree yield.
Both DT and ANN models used the plant-based remote indicators as predictor variables and the amount of cherries produced as the response variable. The mean values of the remote indicators in the 80 trees considered are shown in Table 2. With the aim of designing a DT to estimate whether the tree crop load is going to be higher than usual and so more personnel will be required for harvesting, a threshold value was set at 40 kg tree -1 . High cropping years necessitate more personnel to hand pick the cherries early in their commercial ripeness state and when market prices are higher. The DT aimed to identify two categories: above and below 40 kg tree -1 . Figure 8 shows the resulting DT, which identifies three constraints for the three remote indicators NIR, TRRVI, and NDVI. When the decision tree model was evaluated, it showed a training error of 0.228 and correctly identified 20 trees out of 23 (testing sample). When the confusion matrix was obtained for the testing sample, the precision was 0.870, the recall 0.778, and the true positive rate (specificity) 0.929.  Similarly, an ANN model to predict the sweet cherry yield of each tree based on the interaction among the remote sensing indicators is proposed. After validation of the methodology to generate prediction models, average values of r 2 = 0.63 and RMSE = 6.05 kg tree −1 were obtained. Figure 9 presents a particular case of the relationship between actual and predicted yield obtained by the proposed methodology. As can be seen, there was a strong linear relationship between them. The regression analysis pointed to a good fit (r 2 = 0.67; RMSE = 5.38 kg tree −1 ) between actual and predicted values. Furthermore, the slope of the regression line obtained was close to a 1:1 line.
Remote Sens. 2020, 12, x FOR PEER REVIEW 13 of 18 Similarly, an ANN model to predict the sweet cherry yield of each tree based on the interaction among the remote sensing indicators is proposed. After validation of the methodology to generate prediction models, average values of r 2 = 0.63 and RMSE = 6.05 kg tree -1 were obtained. Figure 9 presents a particular case of the relationship between actual and predicted yield obtained by the proposed methodology. As can be seen, there was a strong linear relationship between them. The regression analysis pointed to a good fit (r 2 = 0.67; RMSE = 5.38 kg tree −1 ) between actual and predicted values. Furthermore, the slope of the regression line obtained was close to a 1:1 line. The results of a statistical analysis of both predicted and actual yield are shown in Table 3. Actual and predicted data had similar mean and maximum values and a CV close to 20%. Regarding the distribution of frequencies, both showed a similar trend towards a positive skew (platykurtic domeshaped trajectory) which, in the case of the predicted tended towards normal (both had Skewness higher than zero and Kurtosis lower than 3). The greatest limitation of the model concerned the minimum values predicted for those sweet cherry trees whose production was between 25 and 34 kg, meaning that the calculation overestimates the yield of those trees with low crop loading ( Figure 9).

Discussion
The results presented in this work showed that remote vegetation indices derived from UAS imagery can successfully assess sweet cherry tree water status; however, these indicators were not as sensitive as conventional and continuous indicators for showing significant differences among irrigation treatments when a slight water stress is applied. In general, remote indices led to significant differences in postharvest, when severe water deficit was applied in the RDI treatments. In preharvest, TRRVI was the only remote index which identified the slight water deficit applied to RDI-1, pointing to significant differences with CTL (Figure 3 L). As the vegetative growth period of sweet cherry trees lasts longer than fruit development, the measurements obtained during preharvest from the NDVI, DVI, OSAVI, and TRRVI resulted in lower values than in postharvest. By contrast, the NDRE index was lower during postharvest [9]. In August, when the postharvest flights took place, sweet cherry trees have already reached their maximum leaf area, while during preharvest (June), The results of a statistical analysis of both predicted and actual yield are shown in Table 3. Actual and predicted data had similar mean and maximum values and a CV close to 20%. Regarding the distribution of frequencies, both showed a similar trend towards a positive skew (platykurtic dome-shaped trajectory) which, in the case of the predicted tended towards normal (both had Skewness higher than zero and Kurtosis lower than 3). The greatest limitation of the model concerned the minimum values predicted for those sweet cherry trees whose production was between 25 and 34 kg, meaning that the calculation overestimates the yield of those trees with low crop loading ( Figure 9).

Discussion
The results presented in this work showed that remote vegetation indices derived from UAS imagery can successfully assess sweet cherry tree water status; however, these indicators were not as sensitive as conventional and continuous indicators for showing significant differences among irrigation treatments when a slight water stress is applied. In general, remote indices led to significant differences in postharvest, when severe water deficit was applied in the RDI treatments. In preharvest, TRRVI was the only remote index which identified the slight water deficit applied to RDI-1, pointing to significant differences with CTL ( Figure 3L). As the vegetative growth period of sweet cherry trees lasts longer than fruit development, the measurements obtained during preharvest from the NDVI, DVI, OSAVI, and TRRVI resulted in lower values than in postharvest. By contrast, the NDRE index was lower during postharvest [9]. In August, when the postharvest flights took place, sweet cherry trees have already reached their maximum leaf area, while during preharvest (June), although canopy development overlaps fruit growth (spur leaves) the maximum canopy development in early cultivars such as 'Prime Giant' is reached in postharvest [42,43].
NDVI was the most sensitive VI of all the remote indicators studied in sweet cherry trees during postharvest, and clearly reflected the result of cumulative water deficits. NDVI was able to distinguish between both RDI treatments and it resulted to have a high sensitivity (Table 1). Indeed, NDVI has been successfully employed in vines to identify different irrigation treatments [2,44]. Our results showed that CTL trees had higher NDVI values than RDI trees as water-stress reduces the capacity of plants to absorb the red energy and reflect NIR. These results are contrary to those reported in apricot trees by Ballester et al. [45], who found higher sensitivity for remote VI based on OSAVI measurements than in NDVI. This might be explained by the fact that the NIR spectral band has been stronger related than any other spectral band with variations in leaf turgidity and canopy structure [46]. In sweet cherry trees, one of the physical mechanisms that first appears in the face of water stress is leaf curling [47]. Such changes in leaf turgidity may lead to NDVI acting as a better water stress detention index in sweet cherry trees than in other fruit trees, whose leaves do not curl so readily. In our study in sweet cherry trees, DVI and OSAVI showed the highest noise of all the remote indicators during preharvest and postharvest, respectively, which decreased its sensitivity (Table 1).
Despite the strong signal intensity of NDVI in both pre-and postharvest periods, it had a high CV in preharvest, meaning that its sensitivity was lower than that obtained by NDRE and TRRVI. However, by postharvest, tree foliage had completely developed and so the CV decreased and NDVI sensitivity increased. The differences described between the preharvest and the postharvest measurements of NDVI are similar to those reported by Bellvert et al. [48] in almond and pistachio trees. Similarly to NDVI, the sensitivity of the NIR spectral band also increased in postharvest. The NIR response is directly related to canopy leaf temperature variations, and water-stressed canopies are warmer than those of well-watered trees, as leaf temperature rises when the stomata close [49]. Thus, the water deficit applied to RDI trees in postharvest led to a stomata closure.
When all the remote indicators were compared with conventional and continuous indicators ( Figure 6), NDVI showed a strong relationship with Ψ s , higher than that mentioned by Zhao et al. [5] in almond trees and lower than that reported by Beeri et al. [50] in citrus and almond trees. Ψ s has been described as the most reliable water stress indicator for assessing water status in sweet cherry trees [15,16]. Consequently, we considered the response of the trees and the relationship between each remote VI and Ψ s in order to assess their usefulness detecting changes caused by the water stress in the canopy structure. The relationship between NDVI and Ψ s was stronger when all treatments were considered ( Figure 5) than when only the CTL treatment was considered ( Figure 4). Furthermore, the relationship was even stronger when all data (2017 and 2018) were taken into account. Neither NDVI, nor any other VI were able to distinguish the small variation measured by Ψ s in CTL trees (Figure 4). NDRE and OSAVI also exhibited a strong relationship with Ψ s , the former showing similar coefficient of determination than that showed in pear trees by Van Beek [51]. As regards the relationship between NDVI and NDRE with MDS, it was similar for both VI and lower than that obtained with Ψ s ( Figure 5G,H and Figure 6). It was not possible to determine the threshold value above which NDVI is not linearly related with MDS.
NDVI was also the remote VI with the strongest relationship with tree biomass, measured as pruning wood mass, followed by DVI, OSAVI, NDRE, and TRRVI (Figure 7). This relationship agrees with the results of Conesa et al. [10] and Matese and Di Genaro [52] in nectarine and 'Sangiovese' grapes, respectively, who found solid relationships between the NDVI computed in a tree and the pruning wood mass produced by it.
Concerning crop yield estimates, the positive results of both DT and ANN obtained by the combination of remote sensing indices highlight the usefulness of UAS imagery for yield prediction in sweet cherry trees. Using both DT and ANN methodologies is not redundant, since each complements and nurtures the other. The DT predictive model showed that NIR, TRRVI, and NDVI indices played a determining role in the classification of tree yield. The results obtained for the confusion matrix and the precision of 0.87, show that the proposed model is a robust method and explains why it has previously been used with good success to predict fruit yield in other fruit trees, such as apple [53]. There were no significant differences among the irrigation treatments in terms of fruit yield, although the results demonstrated that trees with a high NDVI had fruit yields above 40 kg tree −1 (Figure 8). Higher NDVI values during preharvest are not necessarily related to higher postharvest vegetative growth, since CTL trees had higher vegetative growth and pruning wood than RDI trees. High values of NDVI in RDI trees during preharvest could be linked to a trend previously reported in sweet cherry trees submitted to water stress-that their leaf area is greater initially due to a higher number of spurs per meter; however, in postharvest the pruning wood mass of RDI trees was lower as a result of the lower growth of current the season's extension shoots compared to CTL trees [42].
As regards the ANN, the proposed methodology explained 67% of the variation in fruit yield in 2018, a year in which trees showed highly variable crop load values. The distribution of actual and predicted yields did not differ markedly (Table 3), but the approach applying ANN was not able to predict fruit yield in those trees with crop loads under 35 kg tree −1 (Figure 9). Although this could be considered a major restriction of the method, the greatest problem in sweet cherry trees occurs when the fruit yield exceeds 40 kg tree −1 and the needs of fruit-pickers increases. However, it could still be considered an effective model. Nevertheless, the proposed approaches are based on only one year, so further analysis is needed and the model proposed might be a good starting point for future works. Thus, the use of ANN for predicting fruit yield based on airborne hyperspectral imagery in sweet cherry trees provided promising results as has already been reported in other fruit trees, such as mango, avocado, and apple [25,46,54].

Conclusions
The remote VI studied were sensitive to the changes in the canopy of the sweet cherry trees caused by water deficit; however, they were not sensitive enough to detect small tree-to-tree variations within the same irrigation treatment. Moreover, they showed different results according to their phenological stage. In preharvest, when tree canopies were not completely developed none of the remote indicators studied resulted in strong correlation with Ψ s . When we compared the results of VI in postharvest, when leaf area is completely developed, NDVI, NDRE, and TRRVI were able to distinguish between the trees of CTL and RDI; however, NDVI was the only VI able to differentiate between trees of RDI-1 and RDI-2, resulting in the highest sensitivity (S > 22) and was the indicator which resulted in the strongest relationship with Ψ s (r 2 = 0.60). Our results underline the utility of NDVI as a VI in sweet cherry trees.
The predictive models proposed provide growers with more accurate information on sweet cherry tree yield. The decision tree obtained showed a precision of 0.87 when classifying the trees according to their yield, and the application of artificial neuronal networks explained 67% of the variation in fruit yield. The use of both models has positive implications for predicting high crop loads, planning the need for personnel and for scheduling optimal harvesting times.