Combining Thermal and RGB Imaging Indices with Multivariate and Data-Driven Modeling to Estimate the Growth, Water Status, and Yield of Potato under Different Drip Irrigation Regimes

: Advances in proximal hyperspectral sensing tools, chemometric techniques, and data-driven modeling have enhanced precision irrigation management by facilitating the monitoring of several plant traits. This study investigated the performance of remote sensing indices derived from thermal and red-green-blue (RGB) images combined with stepwise multiple linear regression (SMLR) and an integrated adaptive neuro-fuzzy inference system with a genetic algorithm (ANFIS-GA) for monitoring the biomass fresh weight (BFW), biomass dry weight (BDW), biomass water content (BWC), and total tuber yield (TTY) of two potato varieties under 100%, 75%, and 50% of the estimated crop evapotranspiration (ETc). Results showed that the plant traits and indices varied signiﬁcantly between the three irrigation regimes. Furthermore, all of the indices exhibited strong relationships with BFW, CWC, and TTY (R 2 = 0.80–0.92) and moderate to weak relationships with BDW (R 2 = 0.25–0.65) when considered for each variety across the irrigation regimes, for each season across the varieties and irrigation regimes, and across all data combined, but none of the indices successfully assessed any of the plant traits when considered for each irrigation regime across the two varieties. The SMLR and ANFIS-GA models gave the best predictions for the four plant traits in the calibration and testing stages, with the exception of the SMLR testing model for BDW. Thus, the use of thermal and RGB imaging indices with ANFIS-GA models could be a practical tool for managing the growth and production of potato crops under deﬁcit irrigation regimes.


Introduction
Water shortage is one of the key challenges for sustainable agriculture in arid and semiarid countries, particularly under abrupt climate change, as the agricultural sector is highly vulnerable to continuously changing climatic patterns [1]. Dramatic climatic changes in the future are expected to increase water scarcity by approximately 20% globally [2]. Additionally, approximately 80% of the cultivated lands in arid and semiarid generation, a set of estimates is generated by selecting and combining individuals according to their level of accuracy in the problem statement domain, using operators derived from natural genetics [58].
To the best of our knowledge, no previous study has applied the hybrid methodology of ANFIS-GA to thermal and RBG imaging indices to construct robust and predictive models for the accurate estimation and monitoring of the growth, water content status, and production of potato under different irrigation water regimes. Therefore, the present study aimed to (i) quantify the response of growth, water status, and production of potato crop to different irrigation regimes, (ii) examine the performance of thermal and RGB imaging indices for indirect assessment of different traits relating to the growth, water status, and production of each potato variety across all irrigation regimes, for each irrigation regime across potato varieties, and across all conditions (varieties, irrigation regimes and growing seasons); (iii) identify the thermal and RGB imagery indices that explain the most variability in the measured plant traits across all conditions via stepwise multiple linear regression (SMLR); and (iv) evaluate the performance of ANFIS-GA models based on NRCT and RGB indices to predict the different plant traits across all conditions.

Experimental Site, Conditions, and Design
Two field experiments were conducted on drip-irrigated potato crops at the Research Station of the University of Sadat City, Egypt ( Figure 1) during the spring growing seasons of 2019 and 2020. The climate at the experimental site is characterized by warm summers and mild, cold winters with mean air temperatures of approximately 35.0 °C and 12.0 °C, respectively. During both growing seasons, there was very little rainfall and the rainwater did not contribute to the water requirements of the potato crops. The soil texture was sandy loam (69.8% sand, 22.3% silt, and 7.9% clay) with a bulk density of 1.45 g cm −3 , a field The climate at the experimental site is characterized by warm summers and mild, cold winters with mean air temperatures of approximately 35.0 • C and 12.0 • C, respectively. During both growing seasons, there was very little rainfall and the rainwater did not contribute to the water requirements of the potato crops. The soil texture was sandy loam (69.8% sand, 22.3% silt, and 7.9% clay) with a bulk density of 1.45 g cm −3 , a field capacity of 19.2%, and a wilting point of 10.1%. The electrical conductivity of the soil and irrigation water was approximately 1.20 dS m −1 .
In both growing seasons, a split-plot experimental design was used with four replications. Three drip irrigation regimes [100%, 75%, and 50% of the estimated crop evapo-Remote Sens. 2021, 13, 1679 5 of 28 transpiration (ETc)] were randomly assigned to the main plots, and two potato varieties (Bellini and Arizona) were randomly assigned to the sub-plots. The drip irrigation system was divided into four main sectors, each of which represented one replicate and consisted of eight polyethylene lateral drip lines (Twin-wall IV, 16 mm diameter) at the main plot level and four lateral drip lines assigned to each variety at the sub-plot level. Each lateral drip line was 30.0 m long, with a 0.30-m emitter spacing, and was laid out along each potato variety row at a distance of 0.75 m from adjacent lines. The lateral drip lines were connected to the sub-mainline (50 mm diameter) by T-shaped valves, which controlled the water flow to them. The drippers had a discharge rate of 4 L h −1 .
Tubers of each variety that were a uniform and consistent size were planted alongside the drippers (one tuber per dripper) on 1 February 2019 and 30 January 2020 and were dug out on 10 May in both growing seasons. Composted animal manure (75 m 3 ha −1 ), sulfur (300 kg ha −1 as elemental S), and phosphorus [200 kg ha −1 in the form of single furrowbanded calcium superphosphate (15.5% P 2 O 5 )] were applied basally during preparation of the soil for planting, while nitrogen [400 kg ha −1 as ammonium nitrate (33.5% N)] and potassium [100 kg ha −1 as potassium sulfate (50.0% K 2 O)] fertilizers were applied through the fertigation system in four and two equal doses, respectively, from the time of planting to 60 days after planting.

Irrigation Water Requirements
The Food and Agriculture Organization's (FAO's) CROPWAT software v.8 proposed by Smith [59] was used to calculate the amount of irrigation water and to determine the irrigation time for the three irrigation regimes. This software used the modified FAO Penman-Monteith equation reported by Allen et al. [60] and various daily climatic parameters collected from the nearest weather station to estimate the reference evapotranspiration (ETo). The water requirement for the 100% ETc regime was then calculated by multiplying ETo by the crop coefficient for potato (Kc). Because the FAO's recommended values of Kc for potato crops must be adjusted when the minimum relative humidity (RHmin) differs from 45% and the wind speed is greater or less than 2 m s −1 at a height of 2 m, the Kc values were adjusted based on the RHmin and wind speed data for the study area. Based on the final calculations, the total amounts of irrigation water applied for the 100% ETc regime in the first and second seasons were approximately 592.0 and 606.0 mm ha −1 , respectively. These amounts were then reduced by 25% and 50% for the 75% and 50% ETc regimes, respectively. To ensure full germination and good seedling establishment, the plants in each treatment group were irrigated regularly during the first 25 days after sowing, following which the irrigation water was applied according to the prescribed irrigation amount calculated for each irrigation regime.

Thermal Measurements
Thermal images of each sub-plot were captured at the bulking growth stage between 11:00 and 13:00 h using a handheld infrared thermal camera (Ti-32; Fluke Thermography, Glottertal, Germany) ( Figure S1). This camera was equipped with a standard lens with a field of view of 23 • × 17 • and a 320 × 240 pixel micro-bolometer sensor, which detected infrared wavebands in a wavelength range of 7.5-14 µm, which is the optimal range for imaging applications that use heat signatures [61]. The camera had an accuracy of ±0.2 • C in the temperature range −20 • C to 600 • C and a thermal sensitivity of ≤0.05 • C at 30 • C. The emissivity for measurements of the plant canopy was set to 0.96 for the dry reference and 0.95 for the wet reference and the target leaves [30]. The thermal images were captured at a height of 0.80 m above the plant canopy in the nadir orientation. The average leaf temperature within each thermal image was calculated in the SmartView Fluke IR imaging software (version 3.2; Fluke Corporation, Plymouth, MN, USA) as the average values of 50 Remote Sens. 2021, 13, 1679 6 of 28 randomly selected leaves and a polygon area that was fitted around each leaf. These data were then applied to the following equation to calculate the NRCT [33]: where T, T min , and T max are the actual infrared temperature of the canopy of potato plants measured from each sub-plot, the lowest temperature measured during the whole field trial, and the highest temperature measured during the whole field trial, respectively.

Digital RGB Imaging
In parallel with measurement of thermal images, the canopy of potato plants within each sub-plot was also photographed between 11:00 and 13:00 h using a 14-megapixel digital camera (Kodak D5100 reflex; Tokyo, Japan) with a resolution of 2454 × 2056 pixels and 8-bit RGB images. The camera was held manually and oriented vertically downward toward the potato plants at a distance of 1 m. In order to minimize illumination discrepancies between images as well as to ensure homogeneous lightning conditions under clear sky conditions, an umbrella was always held above the camera to shade the area under the camera view. The flash of the camera was always kept off during measurements. The digital images were saved in the jpeg format and then analyzed using the Python software (Version 3.7.3). Because the images sometimes contained non-potato crop canopy features, such as soil, weeds, and straw, the images need to be segmented and extracted to eliminate interference by non-canopy information during feature extraction ( Figure 2).
where T, Tmin, and Tmax are the actual infrared temperature of the canopy of potato plants measured from each sub-plot, the lowest temperature measured during the whole field trial, and the highest temperature measured during the whole field trial, respectively.

Digital RGB Imaging
In parallel with measurement of thermal images, the canopy of potato plants within each sub-plot was also photographed between 11:00 and 13:00 h using a 14-megapixel digital camera (Kodak D5100 reflex; Tokyo, Japan) with a resolution of 2454 × 2056 pixels and 8-bit RGB images. The camera was held manually and oriented vertically downward toward the potato plants at a distance of 1 m. In order to minimize illumination discrepancies between images as well as to ensure homogeneous lightning conditions under clear sky conditions, an umbrella was always held above the camera to shade the area under the camera view. The flash of the camera was always kept off during measurements. The digital images were saved in the jpeg format and then analyzed using the Python software (Version 3.7.3). Because the images sometimes contained non-potato crop canopy features, such as soil, weeds, and straw, the images need to be segmented and extracted to eliminate interference by non-canopy information during feature extraction ( Figure 2). The RGB images were converted to the Hue, Saturation, Intensity color space to extract only the color components, and plant pixels were discriminated from background pixels using fixed thresholds. A binary image of the plant was then obtained by setting the background pixels to 0 and the plant pixels to 1 [62,63]. Thresholding technique is a process of image segmentation to convert the image to grayscale [63] and produce a binary image [63] which has two possible pixel values, namely the intensity value of the image that is more than or equal to the threshold value of 1 (white or foreground) and when it is less than the threshold value or value 0 (black or background) that can be removed. Each pixel in the binary image had a size of 1 bit. The RGB images were converted to the Hue, Saturation, Intensity color space to extract only the color components, and plant pixels were discriminated from background pixels using fixed thresholds. A binary image of the plant was then obtained by setting the background pixels to 0 and the plant pixels to 1 [62,63]. Thresholding technique is a process of image segmentation to convert the image to grayscale [63] and produce a binary image [63] which has two possible pixel values, namely the intensity value of the image that is more than or equal to the threshold value of 1 (white or foreground) and when it is less than the threshold value or value 0 (black or background) that can be removed. Each pixel in the binary image had a size of 1 bit.
RGB imagery consists of three bands, so the color of each pixel is broken down into three values [64]. The values of these RGB channels in plant canopies are affected by several Remote Sens. 2021, 13, 1679 7 of 28 factors, such as the area of green biomass and the chlorophyll content. The RGB color space percentage values were extracted as sample features using the following equations: where R i , G i , and B i are the pixel values for the red, green, and blue bands, respectively, in the digital image; i and S num are the first pixel and maximum number of pixels, respectively; R , G, and B are the mean values of the red, green, and blue bands, respectively. The pseudocode that was used to train each algorithm to extract the RGB color space is presented in Supplementary Materials File S1. The formula and references of the twelve RGB imagery indices that were tested in this study are presented in Table 1. CIVE is mean color index of vegetation.

Plant Trait Measurements
After the thermal and RGB images of the potato plant canopies had been captured, 20 plants from within the imaged area of each sub-plot were randomly harvested from ground level and immediately weighed to determine the biomass fresh weight (BFW). The harvested plants were then cut into small parts, oven-dried to a constant weight at 70 • C, and weighed to determine the biomass dry weight (BDW). The BFW and BDW data were then applied to the following equation to calculate the percentage of biomass water content (BWC): Tubers were harvested when the potato plants in each treatment group reached the onset of senescence. The tubers in each sub-plot were harvested by hand over an area of 22.5 m 2 (two lateral drip lines, 15 m in length), and the total tuber yield (TTY) per ha was calculated and expressed as Mg ha −1 .
2.6. Modeling 2.6.1. Stepwise Multiple Linear Regression (SMLR) SMLR is used to analyze the response of a single dependent variable to two or more independent variables. This method uses forward selection and backward elimination approaches to choose the input (independent) variables (in this case the thermal and RGB Remote Sens. 2021, 13, 1679 8 of 28 imaging indices) for the regression equation in a multi-linear model based on their relative importance, which is determined by their influence, their distribution, and the significance of their effect on a target (dependent) variable (a given measured trait). At each step of the regression process, the input variables are inserted into the SMLR model and checked to distinguish indices that should be remove based on a significant probability level, which is generally set at 0.05 [72]. The process ends when no more indices can be removed or entered.
Extraction of the most influential thermal and RGB imaging indices that accounted for most of the variation in the measured traits (BFW, BDW, BWC, and TTY) using SMLR analysis was based on 75% of the datasets (n = 36). The best-fitting calibration equations for the distinct models with the lowest root mean square error (RMSE) values and highest coefficients of determination (R 2 ) were then used to predict each measured trait based on the remaining 25% of datasets.

Adaptive Neuro-Fuzzy Inference System (ANFIS) Model
ANFIS was developed by Jang [55] by joining ANN with FIS to allow the advantages of a fuzzy system to be combined with a learning algorithm [73,74]. ANFIS uses various output or input data to hypothesize the FIS method and then generally updates its membership criteria by using a backpropagation algorithm. Figure 3 shows a schematic diagram of the ANFIS methodology that was used in this study to predict the various measured plant traits (BFW, BDW, BWC, and TTY).
Remote Sens. 2021, 13, x FOR PEER REVIEW 8 of 28 approaches to choose the input (independent) variables (in this case the thermal and RGB imaging indices) for the regression equation in a multi-linear model based on their relative importance, which is determined by their influence, their distribution, and the significance of their effect on a target (dependent) variable (a given measured trait). At each step of the regression process, the input variables are inserted into the SMLR model and checked to distinguish indices that should be remove based on a significant probability level, which is generally set at 0.05 [72]. The process ends when no more indices can be removed or entered. Extraction of the most influential thermal and RGB imaging indices that accounted for most of the variation in the measured traits (BFW, BDW, BWC, and TTY) using SMLR analysis was based on 75% of the datasets (n = 36). The best-fitting calibration equations for the distinct models with the lowest root mean square error (RMSE) values and highest coefficients of determination (R 2 ) were then used to predict each measured trait based on the remaining 25% of datasets.

Adaptive Neuro-Fuzzy Inference System (ANFIS) Model
ANFIS was developed by Jang [55] by joining ANN with FIS to allow the advantages of a fuzzy system to be combined with a learning algorithm [73,74]. ANFIS uses various output or input data to hypothesize the FIS method and then generally updates its membership criteria by using a backpropagation algorithm. Figure 3 shows a schematic diagram of the ANFIS methodology that was used in this study to predict the various measured plant traits (BFW, BDW, BWC, and TTY). The ANFIS model was integrated with the thermal and RGB indices to predict the measured plant traits and was developed and trained using a GA approach based on a population algorithm. Figure 4 illustrates the ANFIS structural design, which had a multilayered feed-forward configuration that was linked to an incoherent x and y input network. It is essential to mention that the ANFIS model mainly consists of five functional blocks. The Sugeno model ( Figure 4) has a rule base that takes the form shown below [75,76]. The ANFIS model was integrated with the thermal and RGB indices to predict the measured plant traits and was developed and trained using a GA approach based on a population algorithm. Figure 4 illustrates the ANFIS structural design, which had a multilayered feed-forward configuration that was linked to an incoherent x and y input network.
It is essential to mention that the ANFIS model mainly consists of five functional blocks. The Sugeno model ( Figure 4) has a rule base that takes the form shown below [75,76].  Through putting a = A1 and b = B1: while through putting a = A2 and b = B2: where fi is the outputs within the inconsistent zone indicated by the fuzzy principle, Ai and Bi are the membership values, a and b are indirect identifying functions, and pi, qi, and ri are the consequential restrictions that are modified during the forward pass within the learning algorithm. If the membership functions for the fuzzy sets Ai and Bi are µAi and µBi respectively, the five layers that incorporate ANFIS are described in Figure 5. More details about ANFIS can be found in Khadr et al. [77].
In the present study, ANFIS was utilized to determine the relationships between the thermal and RGB indices and the measured plant traits and to characterize them as fuzzy "if-then" regulations. Thus, the premise segment of the fuzzy if-then rules included the thermal and RGB image data, while the consequent segment included the plant traits BFW, BDW, BWC, and TTY. The framework of the developed ANFIS archetype was composed of a Sugeno-type fuzzy system that had a linear output membership function and a typical bell-shaped input membership function. Through putting a = A 1 and b = B 1 : while through putting a = A 2 and b = B 2 : where f i is the outputs within the inconsistent zone indicated by the fuzzy principle, A i and B i are the membership values, a and b are indirect identifying functions, and p i , q i , and r i are the consequential restrictions that are modified during the forward pass within the learning algorithm. If the membership functions for the fuzzy sets A i and B i are µA i and µB i respectively, the five layers that incorporate ANFIS are described in Figure 5. More details about ANFIS can be found in Khadr et al. [77]. GAs are heuristic search approaches that are applicable to a wide range of practical optimization problems. The fundamental principles of GAs were proposed by John Holland in the 1970s and are based on the evolutionary principle of natural genetics. The al- In the present study, ANFIS was utilized to determine the relationships between the thermal and RGB indices and the measured plant traits and to characterize them as fuzzy "if-then" regulations. Thus, the premise segment of the fuzzy if-then rules included the thermal and RGB image data, while the consequent segment included the plant traits BFW, BDW, BWC, and TTY. The framework of the developed ANFIS archetype was composed of a Sugeno-type fuzzy system that had a linear output membership function and a typical bell-shaped input membership function.

Genetic Algorithm (GA)
GAs are heuristic search approaches that are applicable to a wide range of practical optimization problems. The fundamental principles of GAs were proposed by John Holland in the 1970s and are based on the evolutionary principle of natural genetics. The algorithm starts with an arbitrarily generated population of chromosomes and then uses genetic operators (i.e., selection, crossover, and mutation) to generate iterations of fitter chromosomes to optimize the assessment of solutions [78]. Following fitness function evaluation, selection is accomplished via roulette wheel selection, which is analogous to the essential tool of the "survival of the fittest" principle, to form the next generation, with each portion of the roulette wheel being assigned a fitness value. This process is repeated until an acceptable compliance value or unique solution is reached.
In the present study, GA was employed to optimize ANFIS parameters. A basic flowchart describing the GA process is shown in Figure 6. ANFIS includes two parameter types (premise and consequent parameters) that need to be updated and optimized. Premise parameters belong to the Gaussian membership function, which is given as {a i } in Figure 5, and the total number of premise parameters is equal to the sum of the parameters in all membership functions. Consequent parameters are used in the defuzzification layer, which is shown as {p i , q i , r i } in Figure 5. To minimize the error value between the outputs obtained from ANFIS and the measured data, the ANFIS parameters were optimized with GA. The error value of the solution was determined as the RMSE, which demonstrates the fit (absolute) of the model to the data points (i.e., how close the experimental data are to the model's projected values). Thus, a smaller RMSE value indicates a better fit. RMSE was calculated using the following Equation: where P o is the observed parameter value, P f is the predicted value, and n is the number of data points.

Data Analysis
The data of different destructively determined plant traits (BFW, BDW, BWC, and TTY) and the thermal and RGB imaging indices were subjected to analysis of variance (ANOVA) in order to examine the response of these measurements of the two varieties to All thermal and RGB indices were applied as inputs to train the ANFIS-GA models to estimate the measured plant traits. To create the models, the measured datasets were randomly divided into two unequal sets-training datasets (75%) and testing datasets (25%). The performances of the developed models were then evaluated using several evaluation criteria, including R 2 , the Nash-Sutcliffe coefficient (E), the RMSE, and mean absolute deviations (MADs).

Data Analysis
The data of different destructively determined plant traits (BFW, BDW, BWC, and TTY) and the thermal and RGB imaging indices were subjected to analysis of variance (ANOVA) in order to examine the response of these measurements of the two varieties to the different irrigation water regimes. The ANOVA analysis for each season was appropriate for a randomized complete block design in split-plot arrangement, with the irrigation regimes and varieties as the main factor and sub factor, respectively. However, the data of two seasons were tested for homogeneity and normality before analysis using Bartlett's test and the Shapiro-Wilk test, respectively. Because there is homogeneity of error variance for the different measurements, combined analysis of variance over two seasons of data was done. Variety, irrigation regime, and their interaction were considered a fixed effect, while growing season, replicate, and their interaction were considered a random effect. The significant differences between the mean values of the different plant traits and thermal and RGB imaging indices among the irrigation regimes for each variety were separated at the 0.05% probability level using Duncan's multiple range test. These statistical analyses were performed using Statistical Product and Service Solutions (SPSS) Version 15.0 (SPSS Inc., Chicago, IL, USA). The relationships between the plant traits and the thermal and RGB imaging indices for each variety (n = 24), each irrigation regime (n = 16), each growing season (n = 24), and across all conditions (n = 48) were fitted by linear regression using Sigma Plot (v. 11.0; SPSS, Chicago, IL, USA), and the significance of these relationships was tested by calculating the R 2 values. The significance level of R 2 values for all relationships was set at the 0.05, 0.01, and 0.001 probability level.

Combined Analysis of Variance for Measured Plant Traits, and Thermal and RGB Imaging Indices
The results of combined analysis of variance for different measured traits (Figure 7) and thermal and RGB imaging indices (Table 2), after homogeneity test for error variance, showed non-significant differences between seasons. Additionally, the measured traits and indices were not affected by the 2-way interaction of season and irrigation regime or season and variety, or by the 3-way interaction of season, irrigation regime, and variety. The results of combined analysis of variance revealed highly significant differences (P ≤ 0.01) between irrigation regimes for all measured traits and indices. However, the variety main effect, as well as the interaction of irrigation regimes and varieties, was non-significant for all measured plant traits, except TTY, and for all thermal and RGB imaging indices, except NRCT, which both showed significant differences between varieties.

Impact of the Irrigation Regime on the Growth, Water Status, and Production of Potato
Potato plants are very sensitive to soil moisture stress because of their shallow root systems, low root length density, and the fact that they extract most of their water requirements from the upper 30 cm of the soil [7,9]. Therefore, proper irrigation water management is highly necessary to avoid deficit water stress in potato plants, making it very important that the response of potato crops to different irrigation regimes is known with a high level of accuracy. Camargo et al. [79] reported that potato crops growing in Aguas Nuevas, Spain, exhibited the highest biomass accumulation when irrigated with 80-100% ETc and achieved a comparable TTY when irrigated with 80%, 100%, and 120% ETc. Similarly, Foti et al. [80] found that there was no significant difference in the TTY of potato plants irrigated with 100% and 66% ETc, while irrigation with 133% ETc gave the highest TTY. However, several other studies have reported a significant reduction in the dry matter accumulation and TTY of potato crops with a reduction in the amount of irrigation water applied [81][82][83][84], and Cantore et al. [85] reported that 50% of the full irrigation treatment and rain-fed conditions decreased the TTY by 25.9% and 63.6%, respectively, compared with the full irrigation treatment.
In the present study, the 75% and 50% ETc regimes reduced BFW by 24.0% and 48.0%, BDW by 6.5% and 17.9%, BWC by 4.2% and 10.4%, and TTY by 7.5% and 36.8%, respectively, compared with the 100% ETc regime (Figure 7). Thus, it appears that monitoring the growth and water status of potato plants under deficit irrigation regimes could help growers to make better decisions about irrigation scheduling to manage deficit irrigation and enhance crop productivity, particularly under arid and semiarid conditions and in soil with a low water retention capacity as was used in the present study (sandy loam soil with 69.8% sand). As can be seen, the 75% ETc regime resulted in a similar BDW and BWC of the potato plants and only decreased the TTY by 7.5% compared with the 100% ETc regime, despite there being a 25% saving of irrigation water (Figure 7). However, increasing the deficit irrigation stress by irrigating with 50% ETc significantly reduced all of the measured plant traits in both varieties ( Figure 5). This may be because increasing deficit irrigation stress, particularly during sensitive growth stages, significantly reduces the plant water status which, in turn, reduces the elongation and expansion of plant cells and photosynthesis efficiency and induces rapid leaf senescence [79,86]. All of these negative phenomena eventually result in a significant reduction in dry matter accumulation and final tuber yield.
Previous studies have reported that the optimum water requirements for achieving a high TTY in potato crops varies from 350 to 800 mm [15,18,87]. However, Ferreira and Carr [88] reported that the potato variety Desirée required 150-550 mm of irrigation water when growing under hot and dry climate conditions, while El-Abedin et al. [16] showed that the water requirement of potato reached 1505 mm in Saudi Arabia, which is a typical example of an arid country. In the present study, the optimum water requirements for achieving the highest TTY were 592.0 and 606 mm in the first and second growing seasons, respectively. Thus, it is clear that the water requirement for potato crops to reach the target yield can vary greatly depending on several factors, including the irrigation regime, irrigation method, climatic conditions, soil characteristics, and other management and environmental factors. This wide range of water requirements for potato crops between studies further confirms the importance of irrigation water management for enhancing crop productivity under deficit irrigation stress.

Performance of Remote Sensing Indices for Assessing the Measured Plant Traits
Several morpho-physiological plant traits are closely associated with available soil water, so their accurate detection could help determine irrigation thresholds and suitable scheduling, and play an important role in improving irrigation water management and saving irrigation water. For instance, Board and Kahlon [89] reported that the levels of soil moisture stress that cause significant reductions in the seed yield of soybean crops can be predicted through indirect assessment of the total dry matter accumulation in the early growth stage. Thus, the simultaneous and regular assessment of such plant traits using quick, simple, and non-destructive tools could play a vital role in saving irrigation water while achieving the target yield through the proper management of deficit irrigation, as explained in the following sections.

Performance of Remote Sensing Indices for Assessing the Measured Plant Traits
Several morpho-physiological plant traits are closely associated with available soil water, so their accurate detection could help determine irrigation thresholds and suitable scheduling, and play an important role in improving irrigation water management and saving irrigation water. For instance, Board and Kahlon [89] reported that the levels of soil moisture stress that cause significant reductions in the seed yield of soybean crops can be predicted through indirect assessment of the total dry matter accumulation in the early growth stage. Thus, the simultaneous and regular assessment of such plant traits using quick, simple, and non-destructive tools could play a vital role in saving irrigation water while achieving the target yield through the proper management of deficit irrigation, as explained in the following sections.

Performance of Thermal Index
Stomatal conductance plays an important role in regulating the plant water status and balance, particularly under deficit irrigation stress. Therefore, the first response of Figure 7. Varation in the biomass fresh weight (BFW), biomass dry weight (BDW), biomass water content (BWC) and total tuber yield (TTY) of both potato varieties under three irrigation regimes (100%, 75% & 50% ETc). The same letter are not significantly different from one another based on Duncan's multiple range test at a p ≤ 0.05 significance level. Error bars represents standard error of the mean.

Performance of Thermal Index
Stomatal conductance plays an important role in regulating the plant water status and balance, particularly under deficit irrigation stress. Therefore, the first response of most plants to deficit water stress is closure of the stomata to decrease the amount of water lost through transpiration. However, stomatal closure also reduces the growth rate, cell elongation and expansion, and photosynthesis and transpiration rates of plants, which ultimately leads to a substantial reduction in the biomass and final yield [32,90]. Since the transpiration rate is strongly related to plant cooling, stomatal closure for an extended period under deficit water stress leads to a considerable increase in the canopy temperature [29,30,32]. Therefore, a thermal imaging tool that measures changes in the canopy temperature non-invasively by detecting radiation in the infrared region (800-1300 nm) could be successfully used to monitor the growth and water status of plants under different irrigation regimes in a quick, easy, and non-destructive manner. Previous studies have proposed that indices based on the thermal canopy temperature should be used to exploit thermal imaging data for detecting plant traits. In particular, NRCT, which is based on the actual infrared temperature of the plant canopy, has been found to be a useful indicator for water stress monitoring in a range of crops, such as wheat, barley, and soybean, under arid and semiarid conditions [23,33,36].
In the present study, NRCT gradually increased with a decreasing amount of irrigation water (Table 2). Averaged over both seasons and varieties, NRCT was approximately 2.4 and 1.4 times higher under the 50% ETc regime than under the 100% and 75% ETc regimes, respectively, and approximately 1.8 times higher under the 75% ETc regime than under the 100% ETc regime (Table 2). Similarly, previous studies have found that a reduction in the amount of irrigation water applied increases NRCT values in wheat [36] and soybean [23] crops. These findings reflect the fact that the canopy temperature is closely associated with the balance between the amounts of water absorbed by the plant and lost via transpiration. Importantly, if deficit water stress causes the stomata to remain closed for an extended period to inhibit excessive water losses via transpiration, this will not only increase the temperature of the canopy but also effectively reduce the photosynthesis rate and nutrient uptake, ultimately leading to a substantial reduction in plant growth and productivity [32,91]. However, keeping the stomata closed under water deficit stress also allows the plant to avoid dehydration and enhances its water status [92,93]. This close relationship between soil moisture stress and stomata closure combined with the integrated response of several plant properties (e.g., canopy temperature, plant growth, and plant water status) to stomata closure suggest that thermal indices could be effectively used as rapid and non-destructive indirect indicators for the management of deficit irrigation by monitoring the response of different plant properties to soil moisture stress. Means followed by the same letter are not significantly different from one another based on Duncan's multiple range test at a p ≤ 0.05 significance level. NS, **, and *** indicate non-significant, significant at p ≤ 0.01, and highly significant p ≤ 0.001 probability level, respectively. The full names of the abbreviations of RGB imagery indices are listed in Table 1.
To test this hypothesis, the relationships between the thermal index NRCT and various plant traits related to plant growth, water status, and production were evaluated for each potato variety across the three irrigation regimes, for each irrigation regime across the both varieties, and for the pooled data across both varieties, all three irrigation regimes, and both seasons (Figures 8-10). These relationships showed that NRCT was strongly related to BFW, BWC, and TTY (R 2 = 0.79-0.92) and moderately related to BDW (R 2 = 0.25-0.56) when the data for the three irrigation regimes were combined for each variety and when the data for both varieties and all three irrigation regimes were combined (Figures 8 and 9) but was unrelated to BFW, BDW, and BWC (R 2 < 0.1) and only weakly related to TTY (R 2 ≥ 0.20) when the data for the two varieties were combined for each irrigation regime ( Figure 10). These findings indicate that caution is required when thermal indices are used to monitor plant growth and water status under different irrigation regimes across different varieties. This may be because the mechanism of transpirational cooling of the canopy and the ability of the roots to uptake water may differ between varieties under different irrigation regimes, so combining NRCT data across varieties may blur the effectiveness of this index in tracking changes in the traits related to plant growth and water status under different irrigation regimes. Means followed by the same letter are not significantly different from one another based on Duncan's multiple ran at a p ≤ 0.05 significance level. NS, **, and *** indicate non-significant, significant at p ≤ 0.01, and highly significant p probability level, respectively. The full names of the abbreviations of RGB imagery indices are listed in Table 1.   Table 1.
regimes across two years. R 2 values ≥ 0.48 are significant at p ≤ 0.05. The full names of the abbreviations of RGB imagery indices are listed in Table 1. Figure 9. The determination coefficients of relationships between the NRCT and twelve RGB imagery indices with biomass fresh weight (BFW), biomass dry weight (BDW), biomass water content (BWC), and total tuber yield (TTY) of two potato varieties subjected to different irrigation regimes for each season and across two seasons. R 2 values ≥ 0.48 are significant at p ≤ 0.05. The full names of the abbreviations of RGB imagery indices are listed in Table 1.  Table 1.
Remote Sens. 2021, 13, x FOR PEER REVIEW 17 of 28 Figure 10. The determination coefficients of relationships between the NRCT and twelve RGB imagery indices with biomass fresh weight (BFW), biomass dry weight (BDW), biomass water content (BWC), and total tuber yield (TTY) of two potato varieties under 100%, 75%, and 50% ETc. R 2 values ≥ 0.12 are significant at p ≤ 0.05. The full names of the abbreviations of RGB imagery indices are listed in Table 1.

Performance of the RGB Indices
The values of the RGB channels [i.e., the red (R), green (G), and blue (B) band percentages] showed significant variation between the three irrigation regimes for each Figure 10. The determination coefficients of relationships between the NRCT and twelve RGB imagery indices with biomass fresh weight (BFW), biomass dry weight (BDW), biomass water content (BWC), and total tuber yield (TTY) of two potato varieties under 100%, 75%, and 50% ETc. R 2 values ≥ 0.12 are significant at p ≤ 0.05. The full names of the abbreviations of RGB imagery indices are listed in Table 1.

Performance of the RGB Indices
The values of the RGB channels [i.e., the red (R), green (G), and blue (B) band percentages] showed significant variation between the three irrigation regimes for each variety, with the exception of the B% for the Arizona variety. In general, the values of R% and B% increased as the amount of irrigation water applied decreased, while the opposite was true for G% (Table 2). This result indicates that the different irrigation water regimes had a noticeable effect on the color features of the potato crop, highlighting the feasibility of using RGB as a viable and economical monitoring tool for managing deficit irrigation. These color changes can be attributed to the fact that the exposure of plants to water deficit causes noticeable alternations in several anatomical and biophysical plant properties, such as stomatal conductance, guard cell chloroplasts, leaf pigment concentrations, leaf photosynthesis, and leaf water status, which eventually induce substantial changes in their responses to visible light. For instance, the activity of guard cell chloroplasts has been shown to be sensitive to both blue and red light [94,95], and several studies have reported a close relationship between plant water status and the light spectrum in the visible region (400-700 nm) [40,96,97]. In addition, Schlemmer et al. [96] found an inverse linear relationship between the relative water content (RWC) and the blue and red spectral region, while Carter [97] reported that RWC is sensitive to the blue and red visible region centered at 480-690 nm. Thus, it appears that the three channels of visible light (RGB) could be effectively used to monitor several plant traits relating to growth and water status.
In the present study, investigation of the relationship between four plant traits and the R, G, and B percentages showed that the G band was superior to the R and B bands for estimating plant traits, exhibiting strong relationships with BFW, CWC, and TTY (R 2 ≥ 0.75) and a weak to moderate relationship with BDW (R 2 ≥ 0.35) when analyzed for each variety across all three irrigation regimes (Figure 8), and for each season across varieties and irrigation regimes as well as for the pooled data ( Figure 9). However, these three bands failed to exhibit a relationship with any of the plant traits when the data for the two varieties were combined for each irrigation regime ( Figure 10). These findings indicate that leaf color was responsive to water deficit stress mostly in the G band, suggesting that this image band could be useful for reflecting the growth, water content, and production status of different varieties of potato crops. Several previous studies on potato have similarly reported a trend for increased leaf greenness under water deficit stress conditions and have indicated the usefulness of this trait as an important indicator of drought tolerance in potato [98][99][100]. Fleisher et al. [98] also reported that leaf growth of potato crops is reduced by the same amount under short-term and long-term water deficit stress, which could explain the stability of the greenness trait as an indicator of drought tolerance under different water deficit stress scenarios, as well as the failure of the three bands to assess the four plant traits when the data were analyzed for each irrigation regime ( Figure 10). However, the stability of leaf greenness under deficit water stress may be genotype-dependent, which could also explain why the G band was effective for assessing the measured plant traits when the data were analyzed for each variety (Figure 9).
Indices based on RGB bands could also be used to assess plant growth and water status, as well as crop production [101][102][103]. In the present study, the various RGB indices obtained by the digital camera showed significant differences between the three irrigation regimes, with significantly lower values under deficit irrigation (50% ETc) than under full irrigation (100% ETc) ( Table 2). In addition, the majority of the RGB indices obtained exhibited a moderate to strong relationship with BFW, CWC, and TTY and a moderate relationship with BDW when the data were analyzed for each season across varieties and irrigation regimes (Figure 9), under all conditions combined (Figure 9), and for each variety across all three irrigation regimes ( Figure 8). However, all of the indices failed to exhibit a relationship with any of the plant traits when the data for the two varieties were combined for each irrigation regime ( Figure 10). These findings may partly be explained by the different canopy structure, canopy architecture, and leaf angle distribution of the two potato varieties examined, which resulted in different image-derived indices being calculated. However, these differences did not result in any differences in the response of these varieties to the different irrigation regimes in terms of growth or water content, which may explain why there were strong relationships between the measured plant traits and most of the RGB imagery indices for both varieties and across the entire dataset, whereas no such relationships existed for each irrigation regime.
Together, these results show that the RGB imagery indices provided information on the growth, water, and production statuses of the two potato varieties. Similarly, Marín et al. [103] found that the majority of digital camera vegetation indices significantly varied between different turfgrass mixtures growing under different irrigation water regimes, highlighting the value of these indices as a tool for identifying genotypic variation under different growing conditions. In addition, other studies have reported that the vegetation indices obtained from RGB images can provide important information for assessing the performance of different plant traits related to the growth, water content, and production statuses under different growing conditions and that these are somewhat superior to spectral reflectance indices for this purpose [42,101,104,105].

Integration of the Thermal and RGB Imagery Indices for the Assessment of the Measured Plant Traits Using Stepwise Multiple Linear Regression Models
The sensitivities of thermal and RGB imagery indices for estimating plant traits may be affected by several factors, such as the growth stage, growth conditions, plant characteristics, soil background conditions, and irrigation regime, making it difficult to generalize about which is best for a range of conditions-for instance, some indices may be highly sensitive to one plant trait and less sensitive to another. Therefore, in the present study, all of the thermal and RGB imagery indices were included as independent variables in SMLR models to identify which indices explained the most variability in each of the measured plant traits across all conditions (i.e., varieties, irrigation regimes, and seasons). The results showed that NRCT in combination with COM, NRCT together with the R, G, and B band percentages, NRCT alone, and VEG alone explained most of the variation in BFW (89%), BWC (92%), TTY (84%), and BDW (47%), respectively (Table 3). It is also interesting to note that the SMLR models based on the most influential indices showed the best prediction of BFW (R 2 val = 0.84), BWC (R 2 val = 0.89), and TTY (R 2 val = 0.73) but failed to predict BDW (Table 4). This result indicates that SMLR models can also be used to estimate plant traits and that the resulting estimates are better than those obtained from most individual digital camera vegetation indices (see . Furthermore, it is clear that the SMLR models have sufficient predicitivity, particularly for traits related to plant water status and production. This may be because SMLR models are always constructed based on the most influential variables. Similarly, Yue et al. [106] also reported that a new comprehensive growth index estimation model, which combined six plant traits (plant height, aboveground biomass, leaf area index, plant water content, plant chlorophyll content, and plant nitrogen content) and was constructed by SMLR, was better than a model based on a single spectral reflectance index (SRI) or a single RGB imagery index. In general, multivariate regression analyses, which include PLSR, SMLR, support vector machine regression (SVMR), ridge regression, random forest, and ANNs, have yielded better prediction models than SRIs for various plant traits in a range of crops and growing conditions [24,[107][108][109].

Performance of ANFIS-GA Models Based on All Thermal and RGB Imagery Indices for Predicting the Measured Plant Traits
ANFIS is already recognized as a prospective research field of great potential due to its ease of implementation and easy formulation. ANFIS models combine the transparent and linguistic representation of a fuzzy system with the learning ability of an ANN, allowing them to be trained to perform input/output mapping in the same way that an ANN can but with the added advantage of being able to provide the model's rules, which offers greater insight into the modeling process [110,111].  Data-driven modeling has commonly been used to estimate and predict various crop traits, and ANN and SVMR have been used to predict the plant water status, growth, and chlorophyll, nitrogen, phosphorus, and potassium contents [104,112]. However, the integration of various approaches, such as GA, ANN, and FIS, for estimating plant traits based on a combination of thermal and RGB imagery indices has not previously been tested.
In the present study, all 13 thermal and RGB imagery indices were used as inputs for training the ANFIS-GA model to allow the accurate estimation of the measured plant traits. The main purpose of using GA was to find the RMSE value for each output. Once the best performing ANFIS model had been selected through training, the expected values of the chosen output could be determined and compared with the actual observed values. The results showed that the ANFIS-GA models obtained from the thermal and digital imaging data could be used to predict all four measured potato traits, and the use of ANFIS-GA increased the accuracy of the estimates of these parameters (Table 5 and Figures 11 and 12) compared with the use of the NRCT and RGB imagery indices (Figures 8-10).
The performance of the proposed ANFIS-GA model was very high based on the calculated performance criteria, which included R 2 (0.99-1.00), RMSE (0-0.66), MAD (0-0.25), and E (0.99-1.00) ( Table 5), indicating the accurate results of this model. It is also interesting to note that the ANFIS-GA model showed the best testing for BFW (R 2 = 0.88), BDW (R 2 = 0.71), BWC (R 2 = 1.00), and TTY (R 2 = 0.80). The ANFIS-GA model exhibited a very small drop in the quality of the performance measures (R 2 , RMSE, MAD, and E) when moving from the training stage to the testing stage in the case of.
BWC also had an appropriate performance for BFW, BDW, and TTY. However, there was a significant drop in performance for BDW in terms of E, which became negative in the testing stage (   Figure 11. Comparison between training series (a,c) and testing series (b,d) for biomass fresh weight (BFW) and biomass dry weight (BDW) of potato varieties using the developed ANFIS model.  Figure 12. Comparison between training series (a,c) and testing series (b,d) for biomass water content (BWC) and total tuber yield (TTY) of potato varieties using the developed ANFIS model. Figure 12. Comparison between training series (a,c) and testing series (b,d) for biomass water content (BWC) and total tuber yield (TTY) of potato varieties using the developed ANFIS model.

Conclusions
Thermal and RGB imaging indices were combined with multivariate and data-driven modeling in this study in an attempt to non-destructively estimate the growth, water status, and yield of potato crops exposed to different drip irrigation regimes. The results supported the hypothesis that thermal and RGB imaging indices can be considered as an alternative approach to efficient management of deficit irrigation and maximizing the water productivity through indirect assessment of different plant traits in a rapid, easy, and cost effective manner. However, the efficiency of these indices for estimating the plant traits of potato crops depended on whether the data were analyzed specifically for each variety and irrigation regime or were combined across all conditions. Specifically, the thermal index and most of the RGB imagery indices were strongly to moderately related to the measured plant traits when the data were analyzed for each variety across irrigation regimes or for both varieties and all three irrigation regimes combined, but failed to assess the plant traits for each irrigation regime across the two varieties. The results also indicate that the combination of thermal and RGB imagery indices with the ANFIS-GA model can predict the measured plant traits more accurately. This model provided an accurate prediction of the measured plant traits in both the training and testing models (R 2 = 0.71-1.00). In conclusion, the results of this study provide critical insights into the importance of continuously monitoring the growth, water status, and production of potato crops under deficit irrigation regimes using the simple and low-cost approaches of thermal and/or RGB imaging. This type of study has rarely been conducted in irrigated arid and semiarid areas, so the information obtained will be important for supporting the on-going efforts to use drones carrying thermal and RGB devices for undertaking precision agriculture irrigation practices and improving water use efficiency in regions with water scarcity.

Data Availability Statement:
The data presented in this study are fully available in this article and Supplementary File S1.