Remote and Proximal Sensing-Derived Spectral Indices and Biophysical Variables for Spatial Variation Determination in Vineyards

: Remote-sensing measurements are crucial for smart-farming applications, crop monitoring, and yield forecasting, especially in ﬁelds characterized by high heterogeneity. Therefore, in this study, Precision Viticulture (PV) methods using proximal- and remote-sensing technologies were exploited and compared in a table grape vineyard to monitor and evaluate the spatial variation of selected vegetation indices and biophysical variables throughout selected phenological stages (multi-seasonal data), from veraison to harvest. The Normalized Difference Vegetation Index and the Normalized Difference Red-Edge Index were calculated by utilizing satellite imagery (Sentinel-2) and proximal sensing (active crop canopy sensor Crop Circle ACS-470) to assess the correlation between the outputs of the different sensing methods. Moreover, numerous vegetation indices and vegetation biophysical variables (VBVs), such as the Modiﬁed Soil Adjusted Vegetation Index, the Normalized Difference Water Index, the Fraction of Vegetation Cover, and the Fraction of Absorbed Photosynthetically Active Radiation, were calculated, using the satellite data. The vegetation indices analysis revealed different degrees of correlation when using diverse sensing methods, various measurement dates, and different parts of the cultivation. The results revealed the usefulness of proximal- and remote-sensing-derived vegetation indices and variables and especially of Normalized Difference Vegetation Index and Fraction of Absorbed Photosynthetically Active Radiation in the monitoring of vineyard condition and yield examining, since they were demonstrated to have a very high degree of correlation (coefﬁcient of determination was 0.87). The adequate correlation of the vegetation indices with the yield during the latter part of the veraison stage provides valuable information for the future estimation of production in broader areas.


Introduction
Grapevines are one of the most important planted fruit crops, with a significant position in human sustenance. More than seven million hectares of vineyards are grown worldwide in many different geographical areas [1] and are commonly characterized by high heterogeneity, due to physical, biological, and chemical factors, including spatial variations in topography, climatic conditions, and soil characteristics, as well as nonstructural factors, mediated by crop practices [2,3].
Precision agriculture and farm management necessitate monitoring of crops at high spatial and frequent temporal resolution during the entire growing season [4]. Precision can be transformed subsequently into many potential spectral vegetation indices (VIs) and vegetation biophysical variables (VBVs) [41]. VIs and VBVs are important parameters and, therefore, are unique to assess variations in the physiological state and biophysical properties of vegetation. Moreover, they perform crop growth monitoring and yield assessment, detect vegetation stress, and improve crop management practices [19,31,42]. The optimal wavelengths and VIs vary for different crop biophysical parameters and growth stages (GSs). These indices have been widely implemented within remote-sensing applications by using different remote-and proximal-sensing platforms [13,17]. The Normalized Difference Vegetation Index (NDVI) is the most used vegetation index that measures the difference between the canopy reflectance in the near-infrared and visible bands [43,44]. Moreover, numerous innovative and intriguing VIs and VBVs such as the Normalized Difference Red-Edge Index (NDRE), the Modified Soil Adjusted Vegetation Index (MSAVI), the Normalized Difference Water Index (NDWI), the Fraction of Vegetation Cover (FVC) and the Fraction of Absorbed Photosynthetically Active Radiation (FAPAR) have been developed which are related to the vigorous, photosynthetically active radiation absorbed by the greener and healthier leaves. Therefore, are considered as a unique implement to interpret the spatial patterns of the canopy such as changes in canopy size, photosynthetic capacity and canopy chlorophyll content, structure and plant health status, crop productivity, nutrient or water stress and berry characteristics which can be related to changes in microclimate or other conditions [19,20,29,31,[45][46][47][48].
The research occurred in a trial vineyard in the Prefecture of Corinth in Peloponnese (Southern Greece). The Peloponnese wine region is the largest vineyard area in Greece, constituting around 31% of wine production in the country and covering an area of approximately 22,000 ha [49].
Numerous studies the last few decades have tried to connect the characteristics of the vegetation with the yield, in vineyards, in order to estimate their potential productivity, such as those conducted by Borgogno-Mondino et al. [46], Matese and Di Gennaro [2], Dobrowski et al. [50], Ballesteros et al. [51], etc. Nevertheless, the lack of plenty of studies considering table grapes instead of the most frequent examination of wine grapes, as well as the combination of remote and proximal sensing and the assessment of several VIs and VBVs reveals the importance of the present research.
This study aimed to (a) compare and evaluate proximal-and satellite-derived data; and (b) monitor and describe the temporal and spatial variability of spectral and canopybased vegetation characteristics within a vineyard and their relation to productivity, using different sensors, in various growth stages. In this context, a comparison of NDVI and NDRE derived both from remote-and proximal-sensing data sources, for the different data-sources comparison, and additionally, MSAVI, NDWI, FVC, and FAPAR derived only from remote-sensing data for the assessment of their correlation with the vineyard yield, were carried out at four different stages of the cultivation. The analysis considered the heterogeneity of the data regarding the spatial, temporal, and spectral resolution, to avoid accuracy problems.

Study Area
The study area is situated in the prefecture of Corinth (Peloponnese) in Southern Greece (37 • 54.532 N, 22 • 44.798 E, Figure 1a,b), which is characterized by a dry Mediterranean climate with hot summers and mild winters, and according to the Köppen-Geiger climate allocation, it is classified as Csa, with an average annual temperature of 18 • C and rainfall of 550-600 mm [52]. The experimentation was conducted during 2017, in a 1.4-hectare vineyard (Figure 1c). The vines were planted in 2006 with Vitis vinifera L. cv. Thompson seedless, and the variety was grafted onto a 1103 Paulsen rootstock system, having a double cross-arm trellis ordinance with 1.8 m and 2.6 m among the stems and the lines, respectively. Typical canopy management and configuration processes, growth regulators, irrigation water (about 2400 mm/ha), and fertilizers (sixteen foliar applications) were applied during the growing season to achieve profitable conditions for grape diameter and sugar content. Besides, soil physical and chemical properties, such as texture, electric conductivity, pH, organic matter, etc., were measured. The plot displays an average elevation difference of 15-20 m, which is a little higher in the western-northwestern part (Figure 2a) and shows a variance in soil composition, with mainly two different soil types, sandy loam (west and southwest part) and clay-clay loam (east and northeast part; Figure 2b), as well as a small difference in the content of iron ions, as it is higher at the western part.

Proximal and Remote-Sensing Data
Sentinel-2 satellite system carries a Multi-Spectral Imager (MSI) sensor which covers the Visible (VIS), Near Infrared (NIR), and Short-Wave Infrared (SWIR) portions of the electromagnetic spectrum. Five clouds free and atmospherically corrected (2A) images to provide Bottom of Atmosphere reflectance values in cartographic geometry (UTM/WGS84 projection), were acquired free of charge via the ESA portal (https://scihub.copernicus.eu/ accessed on 13 October 2020), for the monitoring of the study area.
The data preprocessing comprised image resampling to 10 m, since S2 spectral bands operate on the different spatial resolution of 10 (4 bands, B2, B3, B4, and B8), 20 m (6 bands, B5, B6, B7, B8A, B11, and B12), and 60 m (3 bands, B1, B9, and B10). Furthermore, given the small size of the study area, a subset and masking of the pixels outside the vineyard boundaries were applied. Finally, the collocation of 10 spectral bands (by excluding bands 1 (coastal aerosol), 9 (water vapor), and 10 (cirrus)) and four supplementary products of the atmospherically corrected images for each date was made [53]. The four supplementary are related to the directional information, including sun zenith (angle θ s ), sun azimuth (angle ϕ), view zenith mean, and view azimuth mean, which are necessary for the calculation of the VBVs (FVC and FAPAR).
The proximal measurements were conducted with a Crop Circle ACS-470 (Holland Scientific Inc., Lincoln, NE, USA) active sensor, which incorporates three optical measurement channels and is user-configurable from 440 to 800 nm, using 12.5 mm interference filters ( Figure 3a). This sensor is widely used for numerous Precision Agriculture applications to record a variety of vegetation indices [30]. The CC ACS-470 provides six spectral bands. In this study, three spectral bands were chosen and configured based on literature reviews and previous research, having as the main objective to be in conjunction with the corresponding bands of S2. Thus, the selected bands cover the red (670 nm), red edge (730 nm), and NIR (760 nm) parts of the electromagnetic spectrum [30,34]. The CC sensor was adapted to a special construction mounted on a quad bike, together with a Garmin GPS16X HVS, (Garmin, Olathe, KS, USA) to georeference the data (Figure 3b). The measurements were collected approximately 1.5 m above the soil surface and 1 m away from the side canopy area, at a rate of 10 readings per second and walking at a constant speed. Spectral-reflectance data were recorded and saved as a text file on an SD flashcard, using the Holland Scientific GeoSCOUT GLS-400 data logger (Figure 3a). A critical step before the analysis was the removal of artifacts within the data that may occur as extraneous data points that lie outside the general range of the dataset (termed outliers), or as data values that differ significantly from neighboring data values but lie within the general range of the data (termed inliers). A protocol on data trimming that was proposed by Taylor et al. [54] was used to ensure the robustness of the analysis and to yield accurate conclusions.
Reflectance data of CC and S2 were collected across the field at three different growth stages divided into five parts (Table 1 and Figure 4): (i) start of veraison (SV); (ii) mid of veraison, which is separated into two parts, MV-1 and MV-2; and (iii) technological maturity, which is also separated into two parts, H-1 and H-2. The processing stage that called veraison (derived from the French expression "verr-ray-zohn") is one of the most important moments in a grapevine's annual lifecycle is the onset of ripening when the grapes turn from green to red and naturally begin to sweeten [18,55,56]. The choice of these concreteness growth stages was made considering that there is a rapid change in berry composition, which is reproduced to the final yield and quality [27,38]. The proximalsensing measurements' dates occurred as close as possible to the dates of the satellite products (Table 1 and Figure 4). Table 1. Acquisition dates of satellite and proximal-sensing data.

Growth Stages
Proximal-Sensing Dates Satellite-Sensing Dates

Spectral Characteristics of Vegetation Indices and Vegetation Biophysical Variables
For the spatial variation comparison of proximal and satellite systems and their relation to the yield, four VIs, and two VBVs were examined. As far as it concerns the VIs, the most commonly used vegetation index, NDVI, and the innovative and intriguing calculations of NDRE, MSAVI, and NDWI were computed [13,33,[57][58][59].
NDVI discriminates the different response of vegetation to the visible and nearinfrared part of the electromagnetic spectrum that is closely related to crop status (Equation (1)) [2,43]. It takes advantage of the strong energy absorption by the chlorophyll in the red portion of the spectrum (RED), and on the high energy scattering by the internal structure of leaves in the near-infrared (NIR) part. It displays a sensitive response to green vegetation even for low vegetation-covered areas and is related to canopy structure and photosynthesis. Though, NDVI is sensitive to the effects of soil brightness and cloud or canopy shadow. The result of NDVI calculation is an image with a continuum of pixel values that range from −1 to 1, where, in general terms, negative values correspond to non-vegetated surfaces, while positive values correspond to vegetated ones, and where higher values are related to healthy photosynthetic vegetation, while lower values, especially those <0.2, are related to stressed vegetation or bare soil [21,60]. Considering that the NIR band 8 of S2 has a different wavelength range from the NIR band that is available from the CC sensor, the S2-NDVI was calculated by using the spectral band 6, which represents reflectance from radiation at 740 nm, very close to the recorded wavelength of the CC NIR band (730 nm).
NDRE is a narrowband greenness VI that was designed to provide a measure of the overall amount and quality of photosynthetic material in vegetation, which is essential for understanding the state of vegetation. Narrowband greenness VIs are considered ideal for use with imaging spectrometers. NDRE was measured by using the NIR and red-edge band (RE; 705 nm), where RE replaces the red band in the equation of NDVI (Equation (2)) [61]. Compared with NDVI, the NDRE index as a widely used red-edge-based VI has been proved to be more resistant to the saturation problem and is more sensitive than NDVI to chlorophyll content in vegetation [33,62,63].
Modified Soil Adjusted Vegetation Index 2 (MSAVI2) is a soil adjusted vegetation index that aims to overcome some of the constraints of NDVI when applied to crops with a medium or high exposure of soil surface (Equation (3)). Moreover, the use of MSAVI2 avoids the main problem of the original soil-adjusted vegetation index, which is required to specify the soil-brightness correction factor (L) through trial-and-error, based on the amount of vegetation in the study area. Qi et al. [64] determined the MSAVI2, which eliminates the need to find the soil line from a feature-space plot or even explicitly specify the soil brightness correction factor. Likewise, it has a simpler algorithm, and it is mainly used in the analysis of plant growth, yield estimation, and Leaf Area Index (LAI) assessment.
NDWI is delineated from the combination of NIR and SWIR reference units which are sensing similar depths through vegetation canopies (Equation (4)) [65]. It is sensitive to changes in vegetation water content and less sensitive to atmospheric aerosol scattering effects. High NDWI values demonstrate a high-water content of the vegetation.
Concerning the VBVs, the FVC and FAPAR biophysical parameters were estimated. FVC defines the ratio of the vertical projected area of vegetation canopy to the reference ground surface, expressed as a fraction [47,66]. FVC is a key variable related to many biophysical features, such as plant phenology, density, and yield [47,67]. FAPAR corresponds to the fraction of photosynthetically active radiation (i.e., within 400 to 700 nm) absorbed by the green parts of the canopy, and it is, consequently, an indicator of the status of the vegetation canopy, since the absorption of the photosynthetically active radiation is related to leaf chlorophyll content [45,68]. FVC and FAPAR calculation require the reflectance of 8 simulated bands (B 3 -B 7 , B 8A , B 11 , and B 12 ) and the cosine of three angles (solar zenith angle, view zenith angle, and relative azimuth angle between solar and view).

Yield Estimation and Statistical Analysis
For the comparison of VIs and VBVs deriving from S2 images with the CC canopy sensor-derived indices and their relation to the yield, a common methodology for upscaling proximal data was utilized [37,69]. Specifically, a grid-layer of 126 blocks of the same size was created with identical geometry and spatial resolution with the grid of satellite images ( Figure 5). This grid was used as a reference layer for upscaling the values generated by proximal and yield data.
The datasets acquired from Crop Circle ACS-470 were varying in terms of density and spatial distribution. Therefore, a geostatistical spatial prediction procedure of block kriging was applied to the dataset gathered for the present study [70]. The output files of CC were imported into a GIS platform for visualizing, editing and calculating the NDVI and NDRE maps. The raster grids of the two indices were upscaled into 10 m pixel size to be comparable with the spatial resolution of S2-derived indices, using the mean value of the pixel situated in the same area. Subsequently, the mean values of the pixels located within each block of the reference grid were computed by the Zonal Statistics tool of the GIS platform, to represent the indices values per block.
Finally, Pearson's correlation coefficient (r) and regression analysis (r 2 ) were used to compare the spatial similarity between the different datasets, by measuring the strength and direction of the linear relationship among the two variables [29,71].
The table grapes were harvested by hand during 16-20 August 2017. The actual yield was estimated at the harvest period by measuring the total number of bins in 50 sampling blocks and multiplying it with the average bin weight of the harvested table grapes. The sampling sites were selected randomly to represent 50 blocks of the reference grid, avoiding the selection of border pixels and inter-row central line, which separates the field into two parts, to reduce the influence of the soil cover.

Processing and Statistical Analysis Software
The processing of the remote-and proximal-sensing data was accomplished by using ArcGIS (Environmental Systems Research Institute, Redlands, CA, USA), ENVI (Exelis Visual Information Solutions, Boulder, CO, USA), and SNAP (STEP, ESA, European Union) software.
The geostatistical analysis of the Crop Circle ACS-470 datasets was accomplished utilizing the VESPER 1.6 (Variogram Estimation and Spatial Prediction plus Error) software of the Australian Centre of Precision Agriculture [70]. The software calculates local semivariograms for each neighborhood through the local kriging process and searches for the data points within the defined site, estimating the variogram cloud by fitting a model. Thus, highly accurate maps of the spatial interpolated values are generated [54].
The statistical analysis of yield and vegetation index correlation was made by utilizing the correlation coefficient estimations and regression analysis and was performed with the statistical software Statgraphics 16 (StatPoint Technologies Inc., Warrenton, VA, USA) [72,73].

Descriptive Statistics
Before performing the correlation analysis, basic statistics were calculated for the exploration of the data. The mean values of the proximal-based NDVI ranged from 0.60 to 0.66, showing a reduced seasonal variability (Table 2) (Table 2). Overall, FVC, FAPAR, and NDVI (proximal-and satellite-based) presented the highest mean values, while lower values, below 0.6, were observed in NDWI and NDRE. Moreover, according to the descriptive statistics of the measured actual yield, the values ranged from 4.5 to 19.8 tons/ha, and there seems to be a positive trend towards the eastern part of the vineyard, where it is characterized by lower production than its western part ( Table 3). The mean was 16.50 tons/ha, which is considered low. The coefficient of variance (CV) reached 18.75%, which is relatively high, considering the small area of the vineyard (1.4 ha).

Satellite and Proximal Based NDVI and NDRE Comparison
The regressions of the NDVI and NDRE datasets, throughout the season, are presented in Figure 6, along with the confidence (green lines) and the prediction limits (gray lines). The output shows the plots of fitting a linear model to describe the relationship between CCNDVI and SNDVI at all growth stages. For the selection of the model, we took into account the table "comparison of alternative models" of Statgraphics to find the best fit. The linear regression model was selected based on the values of R squared, adjusted R squared, and Mean Absolute Error (MAE), since these are the most widely used evaluation metrics for the regression models ( Table 4). The adjusted r 2 of the NDVI datasets is ranging from 0.50 to 0.64 (Table 4), while the lowest regression values (r 2 = 0.49) were recorded on MV-2 and the best fit was reached on MV(2) (r 2 = 0.64; Equations (5)- (9)). In the case of NDRE, the fitted models display 36% (r 2 = 0.36) and 33% (r 2 = 0.36) of the variability in proximal data at the stage of SV and H-1 (Table 4; Equations (10)- (14)). The recorded r 2 between proximal and satellite NDRE datasets were below 0.30 at the other stages of the crop (MV-1, MV-2, H-1, and H-2), indicating poor prediction capability of the models.   The NDVI and NDRE maps were created and classified into three classes (low, medium, and high), using the quantile classing method to demonstrate the spatial distribution of the VIs' values ( Figure 7). Following previous studies [56], the quantile classification method was used only for visual interpretation of VIs and yield. The method was not used to define management zones. There are many classification methods for visualization available, including quantiles, natural breaks, equal intervals, and standard deviations [74,75]. According to ArcGIS help, which provides a detailed description of this schema, the quantile classification is well suited to linearly distributed data, grouping features in equal numbers in each class. To this end, three classes were considered sufficient. Note that if the data distribution has a positive or negative skewness, the quantile or natural break distribution classifiers could be chosen [76]. In the current study, and by considering the distribution histograms of the datasets, the "quantile classifier" was chosen and applied as the most suitable for the visualization of our dataset.
The exploratory analysis of the data revealed considerable spatial variation, and the thematic maps presented some spatial patterns with a trend of higher values occurring at the center and western part of the field, especially in the case of S2 data.

Vegetation Indices and Yield Measurements
The VIs and VBVs maps were also created and classified into three classes (low, medium, and high), using the quantile classing method to demonstrate the spatial distribution of their values (Figure 8). The same classification was also made for the yield spatial distribution (Figure 9).
Thereinafter the correlation analysis of both VIs and VBVs values, from the proximal and remote sensors, at all the selected growth stages with the yield measurements, were conducted. The results display that the CCNDVI and CCNDRE have a lower correlation than SNDVI and SNDRE, respectively ( Table 5). The correlation of CCNDRE was, unexpectedly, extremely low, much lower than the SNDRE, and, in general, the CC values were lower than the satellite (S2) values. The stage of MV-2 appears to be the most suitable date for proximal and satellite measurements (NDVI and NDRE), due to the higher correlation values (0.87 and 0.78 for S, and 0.76 and 0.59 for CC, respectively).
As far as it concerns the correlation of both VIs and VBVs values that derived from S2 data, the maximum positive correlation with the yield was observed from NDVI and FAPAR, mainly at the MV-2 growth stage (0.87) and a suitable positive correlation at the SV and MV-1 stages, varying from 0.74 to 0.78 ( Figure 10). The maximum correlation values of the NDRE occurred also at the MV-2 stage (0.78) and the minimum primarily at the H-1 stage and less at the H-2 stages, 0.60 and 0.67 respectively. (Figure 9).
On the other hand, the other three indices MSAVI2, NDWI, and FVC demonstrate a very low correlation at all stages, and only at the MV-2 stage do they demonstrate a quite acceptable value (Figure 9). In particular, the NDWI and (less so) the FVC display very low correlation values mainly at the first two stages, SV and MV-1, varying from 0.37 to 0.54.

Proximal-and Remote-Sensing Multi-Seasonal Comparison
The use of multi-seasonal remote-and proximal-sensing sensors, which provide specific multispectral data, demonstrates to be an ideal tool to map the vineyard conditions, monitor their development through their growth stages, and assist in the estimation of the production from an early phase of crop growth. Several research studies have been developed, especially during the last few decades, dealing with the implementation of remote and proximal sensors in viticulture, such as those of Sozzi et al. [57], Gatti et al. [77], and Matese et al. [37], which used VIs derived from S2 and UAV sensors for vineyard monitoring. Bramley et al. [32], Taskos et al. [33], Anastasiou et al. [29], Reynolds et al. [78], Primicerio et al. [79], Magarreiro et al. [45], Borgogno-Mondino et al. [46], and Stamatiadis et al. [31] utilized VIs or VBVs obtained from Crop Circle and remotely sensed imagery for assessing vineyard's conditions and their relation to yield variability. Henry et al. [80], Arnó et al. [6], and Stamatiadis et al. [81] used several different proximal sensors, along with GPS and VRA to evaluate the characteristics of vineyards, while other scientists, such as Meggio et al. [19], Xue and Su [42], and Martínez-Beltrán et al. [82], utilized for the same reason various remote sensors (hyperspectral or thermal).
Multi-seasonal measurements during several different growth stages, like the selected in the present study, appear to be a reliable source of information to draw reliable conclusions about the development of the crop, as it is also examined and highlighted by several other previous studies, such as those of Anastasiou et al. [5,38], Kazmierski et al. [56], Lamb et al. [18], Taskos et al. [33], and Mathews [83]. The correlation analysis for each stage of development separately aimed to distinguish the most important period for the development of the crop and its correlation with production.
For this study, the extensively tested and analyzed broadband vegetation index, NDVI and the narrowband index, NDRE, were selected for the comparison of satellite and proximal-derived datasets [24,37,46,77,82,[84][85][86][87]. Based on the results, significant correlations were recorded mainly for the NDVI, where the statistical analysis indicated the robust relationships between the proximal and remote-sensing NDVI datasets, displaying satisfactory linear correlation coefficient results, especially for the period of veraison (MV-2) and harvest (H). These results indicate that the top canopy's vine leaves had higher photosynthetic rates at these stages and provide new insights into the monitoring of vineyards at a medium resolution scale. The fact that the differences between the minimum and maximum values of the VIs derived from the CC were very small probably due to the very high proximity of the CC scanner and the canopy, which affect insignificant the comparison and the relation of the indices. On the other hand, surprisingly, for the NDRE, the correlation was not accordingly suitable. In the case of NDRE, the low correlation with the yield is probably related to the fact that the red-edge band is associated with the visualization of chlorophyll content in leaves that, during the selected growth stages, has reached its maximum values and displays a significant decrease thereafter. Moreover, the narrowband greenness VIs, such as NDRE, are combinations of reflectance measurements sensitive to the combined effects of foliage chlorophyll concentration, canopy leaf area, foliage clumping, and canopy architecture. Thus, it is apparent from the results that these parameters in our case in the specific vineyard field are not ideal to reach a satisfactory correlation [88].

VIs-VBVs and Yield Correlation
As far as it concerns the relationship of the VIs and VBVs with the yield, the highest correlations were detected by NDVI and FAPAR vegetation indicators. Although the NDVI has been commonly used to focus on vineyards' conditions and, consequently, wine production [38,89], the disposal of FAPAR estimates from satellite data releases the possibility to evaluate the state of the vineyard health and production through a variable closer to photosynthesis processes.
On the other hand, NDRE provided quite acceptable results, while MSAVI2, NDWI, and FVC delivered inadequate and ineffective outcomes. MSAVI was tested because NDVI is sensitive to the effects of soil brightness and leaf-canopy shade that, in some cases, necessitate targeted processing and calibration [42]. Therefore, several types of research have mentioned that utilizing medium-resolution satellites to assess vineyard characteristics is risky due to their large heterogeneity and narrow line spacing, and, consequently, only a few studies have been conducted in the field of viticulture with such satellite systems [29,90]. In our experimental field, we have concluded that the gaps between the lines were too small, concerning the spatial resolution of S2 data, to achieve differences in the results by implementing MSAVI2 instead of NDVI, which is more accurate when inter-row soil lines are wide. Moreover, it was ineffective for the same reason as the attempt to separate these lines from the S2 images [57]. On the other hand, the lines were too large for trying to make the scanning of CC from a higher altitude to simulate the acquisition with the perspective view of the satellite sensor. Nevertheless, in the case of PS, a significant reduction of soil reflectance interference was noted considering that the sensor scanning takes place in close proximity to the vine canopy [81]. Additionally, the low correlation values of NDWI are probably related to the lack of irrigation and the low or zero rainfall values during the summer that create extremely dry conditions at the canopy. Their utilization could be more reliable for irrigated crops.
The low correlation values of all indices at the H-1 stage, especially of those derived from the CC sensor, probably may be explained partly by the fact that table grapes received numerous crop practices on the canopy, like the defoliation procedure intended for the grape ripening, which caused a restructuring effect on the canopy characteristics and consequently to the measurements occurred there. Likewise, the correlation between the vegetation index and vineyard yield was higher when the border pixels were removed from the calculation procedure in the case of S2 data, due to the edge effects of the vineyard surrounding area on the spectral response. However, this process is extremely risky, especially at small fields, where the border pixels comprise a large part of the vineyard [57]. This limitation will be eliminated in the following years with the expansion of the very high spatial resolutions satellites and UAVs [2,91].
Furthermore, the highest correlations were found at the same crop stage for each of the proximal and remote methods, during the middle of veraison (MV-2). This outcome agrees with the findings of Anastasiou et al. [29] and Garcia-Estevez et al. [92] which found also the highest correlation during the same growth stage (MV-2). Moreover, the highest correlation of S2-values in comparison to the CC-values, during all growth stages and especially in veraison, is probably related to the physiology of the vines. In particular, the vine's mature leaves are developed near the trunk at the initiation of vine growth, but as the vine grows to senescence, these leaves are found at the edge of the shoot. Consequently, the leaves that are at the side canopy, which are measured with proximal sensing, thus providing lower values in comparison with the S2 that measures the photosynthetic activity of the top canopy leaves.
The fact that the higher correlation was identified at the northwestern and western part of the vineyard almost certainly is related to the geomorphological and soil characteristics, since the experimental field is characterized by a slight slope that extends to the central part of the plot from north to south, and better soil properties at the western part that might affect production and vine conditions.

Conclusions
In this study, an assessment of different sensing methods (satellite and proximal) for monitoring crop conditions through several VIs and VBVs, and their relation to the yield, was conducted, during 2017, on a commercial table grape vineyard. The vegetation indicators were calculated based on the spectral information that was derived from Sentinel-2 and Crop Circle. The results demonstrated the potential of S2 images and CC data to characterize vineyard blocks' vigor and to monitor winegrowers' practices at a territorial (regional) scale.
The results were encouraging in the view of developing a wide-scale monitoring system, especially relying on the Sentinel-2 mission, allowing free access to data with high spatial and temporal resolution from most regions of the world. The conclusions from this study will be also beneficial for improving the applications of Crop Circle in production fields. Moreover, the fact that CC is independent of an external light-source existence, it is not affected by ambient conditions (weather, cloud cover, etc.) makes it a valuable tool for precise viticulture monitoring in any condition, day, or night.
The findings of proximal-and remote-sensing analysis reveal a satisfactory correlation of VIs and VBVs with the yield at specific vine growth stages (MV-2). Agronomists and farmers having this valuable knowledge as a pilot tool can implement precise spatial distributed intervention practices to improve the condition of the plants increasing, and hence, the production.
In general, the results of the present study are encouraging and give the impetus for further research over a longer period, a larger number of experimental plots, and by considering more parameters and/or vegetation indices from comparable sensors. Moreover, an added value of this work is that it presents a perspective for vineyard production estimation using data that are easily accessible, free, and available for numerous grapevines regions in alternative to continuous in situ measurements, which commonly are difficult to access and often associated with extra costs.
Our research attempted to create an operational tool, trying to relate the vineyard vegetation conditions with the yield, using different sensors, in various growth stages, in order to be used as a guide for the production estimation on a broader scale (regional or national). The achievement of such an effort necessitates further research and field measurements to accomplish a trustworthy methodology.