Comparison of Satellite and UAV-Based Multispectral Imagery for Vineyard Variability Assessment

In agriculture, remotely sensed data play a crucial role in providing valuable information on crop and soil status to perform effective management. Several spectral indices have proven to be valuable tools in describing crop spatial and temporal variability. In this paper, a detailed analysis and comparison of vineyard multispectral imagery, provided by decametric resolution satellite and low altitude Unmanned Aerial Vehicle (UAV) platforms, is presented. The effectiveness of Sentinel-2 imagery and of high-resolution UAV aerial images was evaluated by considering the well-known relation between the Normalised Difference Vegetation Index (NDVI) and crop vigour. After being pre-processed, the data from UAV was compared with the satellite imagery by computing three different NDVI indices to properly analyse the unbundled spectral contribution of the different elements in the vineyard environment considering: (i) the whole cropland surface; (ii) only the vine canopies; and (iii) only the inter-row terrain. The results show that the raw s resolution satellite imagery could not be directly used to reliably describe vineyard variability. Indeed, the contribution of inter-row surfaces to the remotely sensed dataset may affect the NDVI computation, leading to biased crop descriptors. On the contrary, vigour maps computed from the UAV imagery, considering only the pixels representing crop canopies, resulted to be more related to the in-field assessment compared to the satellite imagery. The proposed method may be extended to other crop typologies grown in rows or without intensive layout, where crop canopies do not extend to the whole surface or where the presence of weeds is significant.


Introduction
Over the last two decades, precision agriculture (PA) has received significant attention in the agricultural community [1,2].In viticulture, addressing difficulties during the production cycles by defining appropriate crop management, the PA approach has the final aim to improve vineyard yield and grape quality while reducing waste, costs and environmental impact [3,4].
A proper knowledge of the spatial variability between and within crop parcels is considered as a key factor for vine growers to estimate the outcomes in terms of yield and quality [5][6][7].In this context, remote sensing (RS) has already proved its potential and effectiveness in spatiotemporal vegetation monitoring [8][9][10][11].Indeed, data provided by optical sensors of multispectral and hyperspectral imagery systems are profitably exploited to compute a wide set of indices (such as the wide dynamic range vegetation index, the normalised difference red-edge index Index, etc.) by properly describing several crop biophysical characteristics [12,13].In the PA domain, an additional effective application is found for within-field zone management, such as for sink size estimations [14,15] and soil moisture evaluations [16,17], with particular attention to automatic procedures [18][19][20].Among the wide set of defined spectral indices, the normalized difference vegetation index (NDVI) is one of the most extensively used, since it is strictly related to crop vigour and, thus, to the estimated quality and quantity of field production [21][22][23][24][25].
Satellite multispectral imageries (MSI), due to sensors features and platforms altitude, covers extensive areas.In addition, many satellite programmes (such as Landsat, Modis, Aster, SPOT, Sentinel-1 and Sentinel-2) are nowadays providing free datasets, thus promoting satellite imagery exploitation for many agricultural applications [26][27][28][29][30], even with the multi-sensor data fusion approach [31].Examples of valuable research contributions are the low-resolution Modis and high-resolution IKONOS satellite imagery exploitation for mapping vineyard leaf area [23,32].Sentinel-2 offers decametric resolution in terms of space and time, with a ground sample distance (GSD) of up to 10 m and a revisit time of six days.Misregistration of sentinel-2A imageries was addressed in the Processing Baseline (version 02.04), and deployed by European Space Agency (ESA) on 15 June 2016 [33].The effectiveness of decametric resolution satellite imagery in describing crop status and variability, particularly when applied to arable crops, forests and extensive plantations, has been proven by several relevant studies [12,34,35].
However, when considering crops with discontinuous layouts, such as vineyards and orchards, remote sensing becomes more challenging [36].Indeed, the presence of inter-row paths and weed vegetation within the cropland may deeply affect the overall spectral indices computation, leading to a biased crop status assessment.Indeed, novel approaches and algorithms using Unmanned Aerial Vehicle (UAV) or satellite based multispectral imaging have been developed for vegetation pixels classification [37][38][39].Low altitude platforms, such as UAV and airborne sensors, by providing imagery with a high spatial resolution (even a few centimetres) and a flexible flight scheduling [40], allow differentiating between pure canopy pixel and other objects in the scene [41][42][43][44][45] or even to classify different details within canopies [46][47][48].
In this paper, a detailed analysis and comparison of vineyards MSI, provided by a decametric resolution satellite and low altitude UAV platforms, is presented.The effectiveness of the MSI from Sentinel-2 and from the UAV airborne sensors, with very high resolution, was evaluated by considering the well-known relation between the NDVI and crop vigour.In particular, the paper is structured as follows: Section 2 reports information on the considered study area, on data acquisition from the satellite and the UAV platforms, and on the performed data processing to allow a comparison of the NDVI computed from different imagery sources.The results obtained by the data processing and comparison are discussed in Section 3. Section 4 reports the conclusions and future developments.

Materials and Methods
A vineyard located in Serralunga d'Alba (Piedmont, northwest of Italy), covering a surface of about 2.5 ha, was selected as a field test.The cropland, whose latitude and longitude positions range between [44.62334 • 44.62539 • ] and [7.99855 • 8.00250 • ] (World Geodetic System 1984-WGS84), includes three vineyard parcels (named "Parcel A", "Parcel B" and "Parcel C") cultivated with the cv Nebbiolo grapevine, with an area of around 0.36, 0.69 and 0.19 ha, respectively (Figure 1).

Materials and Methods
A vineyard located in Serralunga d'Alba (Piedmont, northwest of Italy), covering a surface of about 2.5 ha, was selected as a field test.The cropland, whose latitude and longitude positions range between 44.62334° 44.62539° and 7.99855° 8.00250° (World Geodetic System 1984-WGS84), includes three vineyard parcels (named "Parcel A", "Parcel B" and "Parcel C") cultivated with the cv Nebbiolo grapevine, with an area of around 0.36, 0.69 and 0.19 ha, respectively (Figure 1).
The vineyard is on a sloped land conformation, with an elevation ranging from 330 to 420 m above sea level and a predominantly southwest orientation.Due to the irregularity of the terrain morphology in terms of altitude and soil properties, the selected vineyard is characterised by a great variation in vine vigour within and between parcels.To extend the study to several vine phenological phases, the acquisition campaigns were performed from April to September 2017.Indeed, vigour varies during the phenological cycle from its minimum, after bud break (April), to its peak in the The vineyard is on a sloped land conformation, with an elevation ranging from 330 to 420 m above sea level and a predominantly southwest orientation.Due to the irregularity of the terrain morphology in terms of altitude and soil properties, the selected vineyard is characterised by a great variation in vine vigour within and between parcels.To extend the study to several vine phenological phases, the acquisition campaigns were performed from April to September 2017.Indeed, vigour varies during the phenological cycle from its minimum, after bud break (April), to its peak in the central part of the season (June-July).Finally, vigour decreases during the final phase of grape ripening.According to the crop stages defined in the BBCH-scale [49], the acquisitions were performed in Stages 57 (inflorescences fully developed), 77 (berries beginning to touch), 85 (softening of berries) and 89 (berries ripe for harvest).The meteorological course in the 2017 season was not aligned to the average temperature and rainfall trends of the region.In particular, March temperature values above the seasonal average led the vine buds to break about 15 days early.A sudden decrease in temperatures in the second half of April caused plant stress (frost damage) in some vines, slowing the vegetative growth of the plants.The phases of flowering and of fruit set, partially influenced by this phenomenon, occurred 10-15 days early.The total rainfall in the 2017 season was 480 mm, much lower than the average of the last 12 years (about 800 mm).This affected grape ripening particularly during August and September, months in which it is refined and completed.The reduced water availability partially contributed in accelerating ripening and in anticipating the veraison and the commercial harvest by about 10 days.

Satellite Time Series Images
The Sentinel-2 satellite is equipped with a multi-spectral imaging sensor that measures Earth's Top of Atmosphere (TOA) reflected radiance in 13 spectral bands ranging from 443 nm to 2190 nm.The technical details of Sentinel-2 along with its spatial resolution and spectral ranges are summarised in Table 1.The Sentinel-2 imagery database, processed at different levels, can be downloaded from [50].In this study, cloud-free level-2A Sentinel-2 Bottom of Atmosphere (BOA) reflectance images were used.The Level-2A imagery was derived from Level 1 by applying scene classification, atmospheric and BDRF correction algorithms, using SNAP toolbox (6.0) and sen2core processor (2.5.5) provided by the ESA [51][52][53][54].Additional details about the Sentinel-2 MSI products can be found in [55].
The selected satellite tiles were acquired on four dates during the 2017 growing season (Table 2) to consider different vegetative vine statuses.Only red and near infrared bands (bands 4 and 8, respectively) were used in this study.The pixels completely included within the boundaries of the three considered "Parcel A", "Parcel B" and "Parcel C" were selected, as shown in Figure 1.All relevant information regarding the satellite imagery processed in this study is organised and summarised in Tables 1 and 2.

Satellite Time Series Images
The Sentinel-2 satellite is equipped with a multi-spectral imaging sensor that measures Earth's Top of Atmosphere (TOA) reflected radiance in 13 spectral bands ranging from 443 nm to 2190 nm.The technical details of Sentinel-2 along with its spatial resolution and spectral ranges are summarised in Table 1.The Sentinel-2 imagery database, processed at different levels, can be downloaded from [50].In this study, cloud-free level-2A Sentinel-2 Bottom of Atmosphere (BOA) reflectance images were used.The Level-2A imagery was derived from Level 1 by applying scene classification, atmospheric and BDRF correction algorithms, using SNAP toolbox (6.0) and sen2core processor (2.5.5) provided by the ESA [51][52][53][54].Additional details about the Sentinel-2 MSI products can be found in [55].
The selected satellite tiles were acquired on four dates during the 2017 growing season (Table 2) to consider different vegetative vine statuses.Only red and near infrared bands (bands 4 and 8, respectively) were used in this study.The pixels completely included within the boundaries of the three considered "Parcel A", "Parcel B" and "Parcel C" were selected, as shown in Figure 1.All relevant information regarding the satellite imagery processed in this study is organised and summarised in Tables 1 and 2. central part of the season (June-July).Finally, vigour decreases during the final phase of grape ripening.According to the crop stages defined in the BBCH-scale [49], the acquisitions were performed in Stages 57 (inflorescences fully developed), 77 (berries beginning to touch), 85 (softening of berries) and 89 (berries ripe for harvest).The meteorological course in the 2017 season was not aligned to the average temperature and rainfall trends of the region.In particular, March temperature values above the seasonal average led the vine buds to break about 15 days early.A sudden decrease in temperatures in the second half of April caused plant stress (frost damage) in some vines, slowing the vegetative growth of the plants.The phases of flowering and of fruit set, partially influenced by this phenomenon, occurred 10-15 days early.The total rainfall in the 2017 season was 480 mm, much lower than the average of the last 12 years (about 800 mm).This affected grape ripening particularly during August and September, months in which it is refined and completed.The reduced water availability partially contributed in accelerating ripening and in anticipating the veraison and the commercial harvest by about 10 days.

Satellite Time Series Images
The Sentinel-2 satellite is equipped with a multi-spectral imaging sensor that measures Earth's Top of Atmosphere (TOA) reflected radiance in 13 spectral bands ranging from 443 nm to 2190 nm.The technical details of Sentinel-2 along with its spatial resolution and spectral ranges are summarised in Table 1.The Sentinel-2 imagery database, processed at different levels, can be downloaded from [50].In this study, cloud-free level-2A Sentinel-2 Bottom of Atmosphere (BOA) reflectance images were used.The Level-2A imagery was derived from Level 1 by applying scene classification, atmospheric and BDRF correction algorithms, using SNAP toolbox (6.0) and sen2core processor (2.5.5) provided by the ESA [51][52][53][54].Additional details about the Sentinel-2 MSI products can be found in [55].
The selected satellite tiles were acquired on four dates during the 2017 growing season (Table 2) to consider different vegetative vine statuses.Only red and near infrared bands (bands 4 and 8, respectively) were used in this study.The pixels completely included within the boundaries of the three considered "Parcel A", "Parcel B" and "Parcel C" were selected, as shown in Figure 1.All relevant information regarding the satellite imagery processed in this study is organised and summarised in Tables 1 and 2.

UAV-Based Imagery
The UAV-based MSI were generated with the Agisoft PhotoScan ® software (Agisoft © , 2018 [56]) processing imagery blocks of more than 1000 aerial images acquired with an airborne Parrot Sequoia ® multispectral camera (Parrot © SA, 2017 [57]).The UAV path was planned to maintain the flight height close to 35 m with respect to the terrain by properly defining waypoint sets for each mission block on the drone guidance platform on the base of the GIS cropland map.With this specification, the aerial images GSD resulted to be 5 cm (Figure 2).A camera geometric calibration procedure was performed before the image alignment task; moreover, a radiometric calibration was applied to the image blocks by using the reference images of a Micasense calibrated reflectance panel [58] acquired before and after each UAV flight.A set of 12 A camera geometric calibration procedure was performed before the image alignment task; moreover, a radiometric calibration was applied to the image blocks by using the reference images of a Micasense calibrated reflectance panel [58] acquired before and after each UAV flight.A set of 12 ground control points, whose positions were determined with a differential GNSS system (with an accuracy of 0.1 m), was placed on selected vine trellis poles within the vineyard to georeference the MSI in a geodetic coordinates frame.
The UAV flights were performed on four different dates over the 2017 crop season (15 May, 29 June, 1 August and 23 September), according to the satellite visiting dates (Table 2).

In-Field Vigour Assessment
The vigour of the vines within the three considered parcels in the study site was evaluated based on the results of a specific in-field survey performed by trained operators and on the past experience of the farmer.In the considered study site, the vigour variability is mainly related to the pedological soil conformation and to the water availability, since irrigation is not allowed by Piedmont regulation.The vigour classification was performed by defining three classes: low "L", medium "M" and high "H".A specific data processing was performed to make the in-field vigour assessment comparable to decametric resolution imagery.In particular, a 10 m × 10 m map was obtained by rastering and clustering a vector GIS map provided by expert agronomists made by a set of three vigour class layers, according to Sentinel-2 pixel location.

Data Processing
In this section, specific methods for data processing, developed to compare and investigate the imagery derived from the two platforms with different spatial resolutions, are presented and discussed in detail.
A tile S, derived from the satellite platform, can be considered as an ordered grid of pixels s(i, j), with indices i and j representing the pixel row and column locations in the raster matrix, respectively.Each pixel s(i, j) was here defined as s(i, j) = [α s (i, j), β s (i, j), n R (i, j), n N (i, j)] T ⊂ S, where α s (i, j) and β s (i, j) are the latitude and longitude coordinates (expressed in WGS84) of the upper left corner of pixel s(i, j), respectively, and n R (i, j) and n N (i, j) are the pixel digital numbers in the red and near infrared bands (12 bit representation), respectively.
Data D derived from the UAV flights were defined as an ordered grid of pixels d(u, v) = [α d (u, v), β d (u, v), m R (u, v), m N (u, v)] T ⊂ D, where pixel d(i, j) coordinates α d (i, j) and β d (i, j) (latitude and longitude in WGS84) are related to the pixels centre and m R (i, j) and m N (i, j) are the pixel digital numbers in the red and near infrared bands (16 bit representation), respectively.
A graphical representation of the defined parameters for the satellite and UAV-based datasets is shown in Figure 3.
Remote Sens. 2019, 01, x FOR PEER REVIEW 7 of 18 ground control points, whose positions were determined with a differential GNSS system (with an accuracy of 0.1 m), was placed on selected vine trellis poles within the vineyard to georeference the MSI in a geodetic coordinates frame.
The UAV flights were performed on four different dates over the 2017 crop season (15 May, 29 June, 1 August and 23 September), according to the satellite visiting dates (Table 2).

In-Field Vigour Assessment
The vigour of the vines within the three considered parcels in the study site was evaluated based on the results of a specific in-field survey performed by trained operators and on the past experience of the farmer.In the considered study site, the vigour variability is mainly related to the pedological soil conformation and to the water availability, since irrigation is not allowed by Piedmont regulation.The vigour classification was performed by defining three classes: low "L", medium "M" and high "H".A specific data processing was performed to make the in-field vigour assessment comparable to decametric resolution imagery.In particular, a 10 m × 10 m map was obtained by rastering and clustering a vector GIS map provided by expert agronomists made by a set of three vigour class layers, according to Sentinel-2 pixel location.

Data Processing
In this section, specific methods for data processing, developed to compare and investigate the imagery derived from the two platforms with different spatial resolutions, are presented and discussed in detail.
A tile , derived from the satellite platform, can be considered as an ordered grid of pixels (i, j), with indices i and j representing the pixel row and column locations in the raster matrix, respectively.Each pixel (i, j) was here defined as (i, j) =  (i, j),  (i, j),  (i, j),  (i, j) ⊂  , where  (i, j) and  (i, j) are the latitude and longitude coordinates (expressed in WGS84) of the upper left corner of pixel (i, j), respectively, and  (i, j) and  (i, j) are the pixel digital numbers in the red and near infrared bands (12 bit representation), respectively.
Data  derived from the UAV flights were defined as an ordered grid of pixels (u, v) =  (u, v),  (u, v),  (u, v),  (u, v) ⊂  , where pixel (i, j) coordinates  (i, j) and  (i, j) (latitude and longitude in WGS84) are related to the pixels centre and  (i, j) and  (i, j) are the pixel digital numbers in the red and near infrared bands (16 bit representation), respectively.
A graphical representation of the defined parameters for the satellite and UAV-based datasets is shown in Figure 3. Ordered grid of pixels s(i, j) belonging to satellite tile S, located at latitude and longitude coordinates α s (i, j) and β s (i, j); and (b) ordered grid of pixels d(u, v) belonging to satellite imagery D, located at α d (i, j) and β d (i, j).Selected UAV pixels belonging to G(i, j), used for comparison to satellite pixel s(i, j), are highlighted in light green.
The evaluation of effectiveness in describing vineyard variability by the satellite and UAV multispectral imagery was focused on the plants vigour assessment by using the NDVI.The NDVI value for satellite pixel s(i, j) can be easily computed as by using the spectral information provided by the digital numbers n N (i, j) and n R (i, j) of the red and near infrared bands.Figure 4a shows the NDV I sat map obtained by applying Equation (1) to the entire set of selected pixels representing "Parcel A", "Parcel B" and "Parcel C" of the Sentinel-2 tile of 7 July.
To allow the comparison of the UAV-based MSI and of the Satellite imagery, a preliminary downsampling procedure of the high-resolution UAV imagery was performed.A portion of UAV dataset D, made by pixel cluster G(i, j), which is related to satellite pixel s(i, j), was defined as With this approach, satellite pixel s(i, j) and the portion of UAV map G(i, j) represent the same section of vineyard cropland, with latitude and longitude coordinates ranging between [α s (i, j + 1)α s (i, j)] and [β s (i, j)β s (i + 1, j)].As an example, an enlargement of UAV map subset G(8, 20), related to satellite pixel s(8, 20) and highlighted in Figure 1 by a yellow square on the field test map, is displayed in Figure 5.
Three specific NDVI indices were defined to perform a detailed analysis of the radiometric information provided by the UAV-based MSI, and then to compare it with the satellite one.In detail, they were computed from the UAV high-resolution data by considering: (i) the whole cropland surface represented by G(i, j); (ii) only the crop canopy pixels and, for completeness; and (iii) only the pixels representing the inter-row terrain.Using all pixels in subset G(i, j), the comprehensive NDV I uav (i, j) for the UAV imagery was defined as By applying Equation (2) to raw UAV map D, an NDV I uav map congruent (properly aligned and with the same spatial resolution) to the ones derived from the satellite imagery (NDV I sat ) can be obtained, as shown in Figure 4b for the UAV imagery acquired on 29 June 2017.
Since within a vineyard UAV orthophoto, with a GSD of 5 cm, the pixels representing the vine canopies can be detected, a more accurate crop NDVI computation with respect to NDV I uav can be performed.For this task, a pixel classification procedure is thus required for each subset G(i, j) to define two different groups of pixels G vin (i, j) and G int (i, j), with G vin (i, j) ∪ G int (i, j) = G(i, j) and G vin (i, j) ∩ G int (i, j) = ∅, representing crop canopies and inter-row surfaces, respectively.The automatic classification procedure described in Comba et al. (2015) [33] was adopted.Figure 5b reports the obtained pixel classification belonging to subset G (8,20) into the two groups G vin (8, 20) and G int (8, 20).By exploiting the spatial information concerning the location and extension of the vine canopies, an enhanced NDVI computation can be defined as An example of the enhanced NDVI definition, by considering only the NDVI of the pixels representing the vine canopies, is reported in Figure 5c, while the complete NDV I vin map for the June dataset is shown in Figure 4c.For completeness, the NDVI index was computed also for the vegetation in the inter-row, such as weed or grass, as to further evaluate the contribution of no-canopy reflectance to the comprehensive NDVI computed from satellite imagery.The obtained NDVI map for the inter-row areas, obtained by processing the UAV imagery of 29 June, is shown in Figure 4d.
Remote Sens. 2019, 01, x FOR PEER REVIEW 10 of 18 For completeness, the NDVI index was computed also for the vegetation in the inter-row, such as weed or grass, as to further evaluate the contribution of no-canopy reflectance to the comprehensive NDVI computed from satellite imagery.The obtained NDVI map for the inter-row areas, obtained by processing the UAV imagery of 29 June, is shown in Figure 4d.

Results
Three vineyard parcels (named "Parcel-A", "Parcel-B" and "Parcel-C"), selected for their peculiar spatial distributions and different micro-climate conditions, were considered in this study

Results
Three vineyard parcels (named "Parcel-A", "Parcel-B" and "Parcel-C"), selected for their peculiar spatial distributions and different micro-climate conditions, were considered in this study (Figure 1).To extend the performed analysis to different vineyard phenological phases, four imagery acquisitions were performed (UAV airborne campaigns) and considered (satellite Sentinel-2 platform) during the 2017 crop growing season (Table 2).Regarding the evaluation and the comparison of the effectiveness of the satellite and UAV-based imagery in describing and assessing the variability within and between vineyard parcels, four different NDVI maps were computed.More in detail, the NDV I sat map (Equation ( 1)) was derived from satellite imagery while imagery acquired with the UAV platform was processed to obtain three different NDVI maps: (i) a comprehensive NDV I uav map (Equation (3)) by considering the spectral information provided by both pixels representing vine canopies and inter-row surfaces; (ii) an NDV I vin map for vineyard canopies (Equation ( 4)); and (iii) an NDV I int map for inter-row paths (Equation ( 5)), by considering only one group of pixels at a time.To be able to compare UAV-based imagery with the satellite imagery, with a GSD of 10 m, high-resolution imagery from the UAV airborne platform was downsampled using Equation (2).The congruence between each spatiotemporal map pair was investigated using the Pearson correlation coefficients, adopted as a map similarity measure [59], after performing a normalisation procedure to focus on the relative differences of each map pair.
A preliminary analysis investigated the coherence of the adopted dataset by comparing the NDV I sat map and the properly downsampled comprehensive NDV I uav map for all three considered parcels and for the four acquisition campaigns.The coherence between the information provided by the two platforms was confirmed by the values of the obtained Pearson correlation coefficients, named R Sat/UAV , which were higher than 0.6 for more than 75% of the performed comparison, and never lower than 0.55.All the obtained R Sat/UAV values, which showed a considerable similarity between the maps derived from the satellite and UAV raw imagery, are organised in Table 3.The correlation plots for imagery pair D 2 /S 2 , detailed for "Parcel A", "Parcel B" and "Parcel C", are shown in Figure 6a-c, respectively.
Remote Sens. 2019, 01, x FOR PEER REVIEW 11 of 18 (Figure 1).To extend the performed analysis to different vineyard phenological phases, four imagery acquisitions were performed (UAV airborne campaigns) and considered (satellite Sentinel-2 platform) during the 2017 crop growing season (Table 2).Regarding the evaluation and the comparison of the effectiveness of the satellite and UAV-based imagery in describing and assessing the variability within and between vineyard parcels, four different NDVI maps were computed.More in detail, the  map (Equation ( 1)) was derived from satellite imagery while imagery acquired with the UAV platform was processed to obtain three different NDVI maps: (i) a comprehensive  map (Equation ( 3)) by considering the spectral information provided by both pixels representing vine canopies and inter-row surfaces; (ii) an  map for vineyard canopies (Equation ( 4)); and (iii) an  map for inter-row paths (Equation ( 5)), by considering only one group of pixels at a time.To be able to compare UAV-based imagery with the satellite imagery, with a GSD of 10 m, high-resolution imagery from the UAV airborne platform was downsampled using Equation (( 2)).The congruence between each spatiotemporal map pair was investigated using the Pearson correlation coefficients, adopted as a map similarity measure [59], after performing a normalisation procedure to focus on the relative differences of each map pair.
A preliminary analysis investigated the coherence of the adopted dataset by comparing the  map and the properly downsampled comprehensive  map for all three considered parcels and for the four acquisition campaigns.The coherence between the information provided by the two platforms was confirmed by the values of the obtained Pearson correlation coefficients, named R / , which were higher than 0.6 for more than 75% of the performed comparison, and never lower than 0.55.All the obtained R / values, which showed a considerable similarity between the maps derived from the satellite and UAV raw imagery, are organised in Table 3.The correlation plots for imagery pair  / , detailed for "Parcel A", "Parcel B" and "Parcel C", are shown in Figure 6a-c   Once the coherence of the adopted dataset was verified, the quality of the information provided by the two typologies of NDVI maps was investigated, focusing on the well-known relationship between the considered index and the vegetative condition of the vineyard.This task was performed by comparing NDV I sat and NDV I uav maps to the in-field vigour assessment provided by expert operators, which classified the different regions of the considered vineyard into three different vigour classes (Figure 7).The results of the performed analysis of variance (ANOVA), which provided p-values from 0.04 to 0.26, did not prove significant difference between vigour groups means for all the considered parcels (Figure 8a).Once the coherence of the adopted dataset was verified, the quality of the information provided by the two typologies of NDVI maps was investigated, focusing on the well-known relationship between the considered index and the vegetative condition of the vineyard.This task was performed by comparing  and  maps to the in-field vigour assessment provided by expert operators, which classified the different regions of the considered vineyard into three different vigour classes (Figure 7).The results of the performed analysis of variance (ANOVA), which provided pvalues from 0.04 to 0.26, did not prove significant difference between vigour groups means for all the considered parcels (Figure 8a).values derived from UAV imagery, considering only canopy pixels, divided into three groups on the basis of the observed in-field vigour classes "L", "M" and "H".

Discussion
The decametric resolution satellite imagery showed some limitations in directly providing reliable information regarding the status of vineyards where the crop radiometric information can be altered by other sources (e.g., inter-row paths) that, in the case of crops grown by rows, could be predominant and could negatively affect the overall assessment.This effect was confirmed by the strong relation found between the satellite  map and the  map, derived by the UAV

Discussion
The decametric resolution satellite imagery showed some limitations in directly providing reliable information regarding the status of vineyards where the crop radiometric information can be altered by other sources (e.g., inter-row paths) that, in the case of crops grown by rows, could be predominant and could negatively affect the overall assessment.This effect was confirmed by the strong relation found between the satellite NDV I sat map and the NDV I int map, derived by the UAV imagery considering only the inter-row pixels, with Pearson correlation coefficients R Sat/int of the 12 performed comparisons ranging from 0.49 to 0.67 (Table 3) and with more than 65% of spatiotemporal map pairs achieving R Sat/int > 0.60.On the contrary, the relation between the satellite NDV I sat map and the enhanced NDV I vin map, obtained by considering only pixels representing vine canopies within the UAV imagery, resulted to be weak.About 75% of the Pearson correlation coefficients R Sat/vin , obtained by comparing NDV I sat and NDV I vin maps of all three considered parcels and the four acquisition campaigns, were found to be lower than 0.41.All the obtained values of R Sat/vin are organised in Table 3.This analysis proves that, in the case of crops where the inter-row surfaces and paths involve a relevant portion of the cropland, such as vineyards, the radiometric information acquired by decametric spatial resolution satellite platforms are not sufficient to properly evaluate crop status and variability.Indeed, depending on the specific adopted crop management approach, the inter-row surface can be covered by grass, by other crops for pest and disease integrated control or it can be bare soil.In all these situations, the vineyard vigour could often be in discord with the inter-row areas, leading to biased vineyard vigour assessments from decametric spatial resolution imagery.
In addition, the effectiveness of the NDV I vin and NDV I sat maps in discriminating the vigour of the considered parcels in accordance with the experts in-field assessment was investigated with the ANOVA method.For what concern the NDV I vin map, the ANOVA, obtaining p-values ranging from 2.47 × 10 −3 to 6.88 × 10 −8 (Table 4), confirmed that the observed variability of the vineyards within the test site was well described by the NDV I vin map.Boxplots of the NDV I vin values, divided in the three vigour classes, are shown in Figure 8b.On the contrary, the ANOVA results considering NDV I sat map showed that the variability of the vineyards within the test site by satellite platform was not significantly in accordance with the experts in-field assessment (Table 5).This additional verification confirmed the main result of the presented analysis, proving that, in the case of crops where the inter-row surfaces involve a relevant portion of the cropland, such as vineyards, the radiometric information acquired by satellite platforms can have difficulties properly evaluating crop status and variability.In these situations, imagery with a high spatial resolution is required to properly assess variability within and between vineyards.

Conclusions
In this paper, a detailed analysis and comparison of multispectral imagery of vineyards, provided by decametric resolution satellite and low altitude UAV platforms, is presented.The effectiveness of the Sentinel-2 imagery and the high-resolution UAV aerial images was evaluated by considering the well-known relation between the NDVI and vineyard vigour.A cropland located in Piedmont (Serralunga d'Alba, Italy) was selected as the experiment site to perform four image acquisition campaigns, which were properly scheduled according to the main vine phenological stages.
The results show how, in the case of crops where the inter-row surfaces involve a relevant portion of the cropland, such as vineyards, the radiometric information acquired by decametric resolution satellite platforms has difficulties in properly evaluating crop status and variability.In these situations, the vigour of the vineyard could often be in discord with the inter-row areas (e.g., grass, plants for pest and disease integrated control, or soil), leading to biased vineyard vigour assessments from decametric resolution imagery, such as the Sentinel-2 imagery.This was proved by a detailed analysis of the radiometric unbundled contribution of different elements within the cropland, performed by defining three different NDVI indices from the high-resolution UAV imagery, considering: (i) the whole cropland surface; (ii) only the crop canopy pixels; and (iii) only the inter-row terrain pixels.
In this study, the NDVI maps derived from the satellite imagery were found not to be in accordance with the in-field crop vigour assessment.In addition, the satellite-based NDVI maps were found to be more related to the NDVI maps computed by the high-resolution UAV imagery, considering only the pixels representing inter-row surfaces.As a validation, a new type of NDVI map from the UAV imagery, generated by considering only the pixels representing the vine canopies, was defined.The effectiveness of this last type of map in describing the observed vineyard vigour was found to be relevant.
The proposed approach can be extended to other crop typologies that are grown by rows or without intensive layouts, where the crop canopies do not extend on the whole surface or where the presence of weeds is relevant.

Figure 1 .
Figure 1.Selected test field located in Serralunga d'Alba (Piedmont, northwest of Italy).The boundaries of the three considered parcels, named "Parcel-A", "Parcel-B" and "Parcel-C", are marked with solid green polygons.The cropland region, represented by pixel  , of the Sentinel-2 tile, is highlighted by a yellow square.The map is represented in false colours (NIR, Red and Green channels).

Figure 1 .
Figure 1.Selected test field located in Serralunga d'Alba (Piedmont, northwest of Italy).The boundaries of the three considered parcels, named "Parcel-A", "Parcel-B" and "Parcel-C", are marked with solid green polygons.The cropland region, represented by pixel s 8,20 of the Sentinel-2 tile, is highlighted by a yellow square.The map is represented in false colours (NIR, Red and Green channels).

Figure 3 .Figure 3 .
Figure 3. (a) Ordered grid of pixels (i, j) belonging to satellite tile  , located at latitude and longitude coordinates  (i, j) and  (i, j); and (b) ordered grid of pixels (u, v) belonging to Figure 3. (a)Ordered grid of pixels s(i, j) belonging to satellite tile S, located at latitude and longitude coordinates α s (i, j) and β s (i, j); and (b) ordered grid of pixels d(u, v) belonging to satellite imagery D, located at α d (i, j) and β d (i, j).Selected UAV pixels belonging to G(i, j), used for comparison to satellite pixel s(i, j), are highlighted in light green.

Figure 4 .
Figure 4. Comprehensive (a)  map, computed from satellite imagery  , and (b)  derived from UAV imagery  .(c) Enhanced vineyard  map, processing UAV imagery  by considering only canopy pixels  and (d)  map considering only inter-row surface  .In all represented NDVI maps, only pixels (i, j) completely included within "Parcel A", "Parcel B" and "Parcel C" boundaries are shown.

Figure 4 .
Figure 4. Comprehensive (a) NDV I sat map, computed from satellite imagery S 2 , and (b) NDV I uav derived from UAV imagery D 2 .(c) Enhanced vineyard NDV I vin map, processing UAV imagery D 2 by considering only canopy pixels G vin and (d) NDV I int map considering only inter-row surface G int .In all represented NDVI maps, only pixels (i, j) completely included within "Parcel A", "Parcel B" and "Parcel C" boundaries are shown.

Figure 5 .
Figure 5. (a) Enlargement of subset G(8, 20) of UAV map D 2 , highlighted by a yellow square in Figure 1, is represented in false colours (NIR, Red and Green channels); (b) classification of pixels d(u, v) ⊂ G(8, 20) into two classes: G vin , representing vine canopies (green), and G int , representing inter-row surfaces (brown); (c) computed NDVI values of vine canopies pixels G vin ; and (d) inter-row surface G int .

Figure 6 .
Figure 6.Scatter plots of NDVI values from NDV I sat map (x-axis) and: (a) the comprehensive NDV I uav map (y-axis); (b) the enhanced NDVI values of map NDV I vin (y-axis); and (c) the enhanced NDVI values of NDV I int map (y-axis), using imagery pair data D 2 /S 2 .The regression model and data pair correlation coefficients are also reported.

Figure 7 .
Figure 7. Vineyard test site classification into three vigour classes on the basis of the observed in-field vigour assessment.Classes "L", "M" and "H" refer to low, medium and high vigour, respectively.

Figure 7 . 18 Figure 8 .
Figure 7. Vineyard test site classification into three vigour classes on the basis of the observed in-field vigour assessment.Classes "L", "M" and "H" refer to low, medium and high vigour, respectively.Remote Sens. 2019, 01, x FOR PEER REVIEW 13 of 18

Figure 8 .
Figure 8. Box plots representation of: (a) NDV I sat values derived from satellite imagery; and (b) enhanced NDV I vin values derived from UAV imagery, considering only canopy pixels, divided into three groups on the basis of the observed in-field vigour classes "L", "M" and "H".

Table 1 .
Technical details of the considered and adopted platforms and sensors.

Table 1 .
Technical details of the considered and adopted platforms and sensors.

Table 1 .
Technical details of the considered and adopted platforms and sensors.

Table 2 .
Information on satellite and UAV based acquired datasets.

Table 4 .
Results of the ANOVA of UAV based NDV I vin map in relation to the three vigour classes from in-field assessment.

Table 5 .
Results of the ANOVA of satellite based NDV I sat map in relation to three vigour classes from in-field assessment.