Vineyard Detection and Vine Variety Discrimination from Very High Resolution Satellite Data

In order to exploit remote sensing data operationally for precision agriculture applications, efficient and automated methods are required for the accurate detection of vegetation, crops and different crop varieties. To this end, we have designed, developed and evaluated an object-based classification framework towards the detection of vineyards, the vine canopy extraction and the vine variety discrimination from very high resolution multispectral data. A novel set of spectral, spatial and textural features, as well as rules, segmentation scales and a set of parameters are proposed based on object-based image analysis. The validation of the developed methodology was carried out on multitemporal WorldView-2 satellite data at four different viticulture regions in Greece. Concurrent in situ canopy reflectance observations were acquired from a portable spectroradiometer during the field campaigns. The performed quantitative evaluation indicated that the developed approach managed in all cases to detect vineyards with high completeness and correctness detection rates, i.e., over 89%. The vine canopy extraction methodology was validated with overall accuracy (OA) rates of above 96%. The quantitative evaluation regarding the vine variety discrimination task, including experiments with up to six different varieties, reached OA rates above 85% at the parcel level. The combined analysis of the experimental results with the spectral signatures from the in situ reflectance data indicated that certain vine varieties (e.g., Merlot) presented distinct spectral patterns across the VNIR spectrum.


Introduction
Multispectral sensors onboard satellite, aerial and UAV platforms are nowadays the main source of remote sensing information for precision agriculture and environmental monitoring applications.Due to certain key factors, like aerial flight regulations and cost, satellite remote sensing is, still, providing operationally cost-effective geospatial vegetation maps with up to thirty centimeters of spatial resolution.At the same time, open data policies both in the USA and EU and commercial Earth observation satellite missions are delivering an unprecedented volume of data with an increasing level of detail and accuracy.
In order to exploit such Earth observation big data operationally, for agriculture and environmental applications, efficient and automated processing methodologies are required that can accurately detect vegetation, crops and crop varieties in remote sensing images.To this end, recent studies have employed high and very high resolution satellite imagery towards vegetation, forest/tree mapping [1][2][3][4], biomass and structural parameters estimation [5][6][7] and crop type mapping and identification [8,9].Moreover, crop-based and variety-based data analysis can create valuable validated agricultural maps and products for the implementation of effective management decisions [10,11].In particular, for viticulture applications, the detection of vineyards from satellite remote sensing data and the remote discrimination of vine varieties, including numerous Vitis vinifera L. varieties that are cultivated worldwide, is of major importance both for farmers and government activities.
For more than a decade now [12][13][14][15][16], research efforts have been focusing on the exploitation of satellite image data at different spatial, spectral and temporal resolutions towards the efficient analysis and monitoring of vine crops through remote sensing sensors and techniques.In particular, apart from the within-field crop analysis and the estimation of the spatial variability in wine-grape composition and yield [17][18][19][20][21][22][23], there is plenty of research towards the classification, identification and delineation of crops in remote sensing data [24][25][26][27].However, despite recent research efforts towards the detection and delineation on medium [28][29][30] and high resolution [31,32] spatial scales, the development and validation of efficient classification frameworks for operational vineyard detection in high resolution data and over large agricultural regions still remain a challenge.
Concurrently, another more complex problem is the accurate discrimination of vine varieties from remote sensing data, a task of significant importance for farmers, wine producers, agronomists, geospatial engineers and public agencies for numerous reasons, such as the production of variety-specific maps, misplanting detection, certification labeling, vine-variety database creation, planning, etc.Therefore, along with the in-field analysis, the different mapping products (like canopy, vigor, maturity, etc.) must take into account and be optimized according to the variety [33][34][35][36][37][38], since the overall management, grape/wine quality potential and respective economic value is variety dependent [39,40].To this end, apart from the conventional analytical techniques employed to discriminate several vine or other plant varieties, such as ampelometry [41], wet chemistry methods or DNA analysis [42], the exploitation and validation of non-destructive remote sensing tools has not been much studied.Laboratory experiments based on leaf/fruit reflectance spectroscopy and statistical analysis indicated that certain vine [40] and strawberry [43] varieties were detected and separated among other ones.
Most other studies have mainly employed aerial hyperspectral datasets from manned aircrafts.Lacar et al. [44] experimented with the discrimination of two different vine varieties i.e., Cabernet Sauvignon and Syrah.Using CASI (Compact Airborne Spectrographic Imager) hyperspectral data over a vineyard in the Barossa Valley in South Australia and the spectral statistical analysis from image sampling results, this study indicated a certain variation in the spectral signatures of the two varieties in the visible region.Classification results employing 12 spectral bands from the VNIR resulted in high discrimination rates between the two varieties.Under a similar experimental set-up, Ferreiro-Arman et al. [45] employed a CASI hyperspectral dataset with spatial resolution at 3 m and 144 spectral bands over the Domain Château Luchey-Halde (Bordeaux, France).Different supervised classification techniques were tested for the discrimination of two vine varieties (i.e., Merlot Noir and Cabernet Sauvignon) and also the classes of soil and trees, resulting in high OA rates, i.e., around 90%.Moreover, based on a CASI-2 aerial hyperspectral dataset over the Pessac-Léognan region (Bordeaux, France) and on a supervised classification procedure, several discrimination experiments were performed under constant and varying within-scene illumination conditions [46].The overall evaluation indicated high (above 82%) OA rates for binary classification between pairs of, i.e., Merlot Noir, Cabernet Sauvignon and Cabernet Franc.
Towards the same direction, but addressing the problem for very high resolution satellite data, we have designed, developed and validated a remote sensing classification procedure towards the efficient vineyard detection and vine variety discrimination.We have employed an object-based classification framework [47,48] and developed an efficient set of features and rules that are able to effectively identify vineyards in high resolution multispectral pan-sharpened data in different viticulture regions.Moreover, the developed framework was validated for the extraction of the vine canopy and the discrimination of vine varieties through numerous experiments towards the separation of three and up to six different vine varieties/clones.

Contribution
Regarding the developed vineyard detection methodology, in contrast to similar research efforts towards vineyard and/or other crop identification tasks [8,9,15,[22][23][24]27,30], the proposed approach employs very high resolution multispectral data, along with a specific set of features, rules, segmentation scales and a set of parameters that deliver relatively higher detection rates.Instead of employing only spectral features (e.g., [8,22,23,27]), we employ, as in [9,24], textural features, as well.However, we solve the problem for very high resolution pan-sharpened data (at 50 cm) instead of the significantly coarser Landsat 8 and ASTER data (15 m/30 m).
Regarding the developed canopy extraction methodology in contrast to similar efforts that are based on (e.g., NDVI) thresholds and/or line detection/fitting [28,32,49], we formulate the problem under a supervised object-based classification framework and propose a set of spectral and spatial features that can extract the canopy (not just lines, center-lines, etc.) while addressing different in-between the rows materials, like soil, weeds, rock, etc.
Last, but not least, we propose a novel framework for vine variety discrimination based on a fine segmentation scale and a set of spectral features that can discriminate different varieties in very high resolution multispectral data.It is worth mentioning that similar efforts [44][45][46] were based on airborne hyperspectral data and focused on discriminating two vine varieties through classification.

Study Areas
Data acquisition campaigns were conducted in four different viticulture regions in Greece (from the central to north regions) all belonging to local Protected Geographical Indication (PGI) zones (Figure 1).Concurrent satellite and field campaigns were scheduled during the veraison period, from early July to early August, for the years 2012 and 2013.The first study area, Trapeza, is located in the Northern Peloponnese, inside the PGI zone Plagies Aigialeias, at an altitude of about 800 m, over the Trapeza village.Eight different vine varieties are cultivated in this particular region with Syrah, Stafida and Sauvignon Blanc among the most important for winemaking.The second region, Megaplatanos, is located in Central Greece, inside the PGI zone Koilada Atalantis, at an altitude of about 60 m above sea level.More than 20 different vine varieties are cultivated in this area.During this work, we focused mainly on the study of six of those, including Chardonnay, Robola and Cabernet Sauvignon.The third study area, Naoussa, is located in Western Macedonia, inside the PGI zone Naoussa and more specifically near the Giannakochori village.The altitude of the region is about 200 m above sea level.More than nine different varieties are cultivated in the region.Xinomavro, Merlot and Syrah are the dominating varieties for winemaking purposes.The forth study area, Amynteo, is located in Western Macedonia, inside the PGI zone Florina near the Amynteo village, at an altitude of about 600 m above sea level.Various varieties are cultivated in the region.During this study, we focused mainly on the Syrah, Tannat and Montepulciano varieties.

Satellite and Field Campaigns
Very high resolution panchromatic and multispectral satellite imagery was acquired from the WorldView-2 ( WV-2) satellite sensor, in GeoTIFF format of 16-bit color radiometric resolution.Images were acquired during the veraison period and in particular on 31 July 2012 for the Trapeza and Megaplatanos regions, on 29 July 2013 for the Naoussa study area and on 1 August 2013 for the Amynteo region.The WorldView-2 satellite carries an imaging instrument containing a high resolution panchromatic band with a reduced infrared and blue response and eight lower spatial resolution spectral bands.WorldView-2 captures panchromatic imagery of 0.5-m spatial resolution and 8-band multispectral imagery of 2-m resolution.Apart from the standard RGB and near infrared (NIR) spectral bands, four additional bands were available in all cases i.e., the coastal (centered at 425 nm), the yellow (centered at 605 nm), the red-edge (centered at 725 nm), which is on the onset of the high reflectivity portion of vegetation response, and a second NIR (centered at 950 nm).
In all study areas, concurrently with the satellite campaigns, intensive field campaigns were conducted in order to collect reference/ground truth data including the precise location, boundaries, variety and clone of each parcel.Existing maps with geographic information and varietal plantation were verified or updated during field surveys, with the aid of the farmers and agronomists in charge.In situ reflectance was also obtained using the GER 1500 (Spectra Vista Corporation) portable spectroradiometer, which provides spectra with 512 spectral bands, distributed in the spectral region from 350 nm to 1050 nm with 3.2-nm FWHM.More than a thousand measurements were acquired, while each measurement was the average of three successive ones, computed internally by the instrument.The instrument was held approximately 1 m above the vine canopy, resulting in an effective scene size of about 190 cm 2 .For each series of measurements, the irradiance was also obtained through measurements on a Lambertian reference panel.The position of each measurement was recorded using a portable GPS.

Data Pre-Processing
Radiometric and atmospheric corrections were performed on the satellite imagery towards the elimination of solar illumination, atmospheric and terrain effects.Digital numbers of relative radiance were converted to absolute radiance using the absolute radiometric calibration factors and effective bandwidths for each band.Atmospheric correction was conducted through ENVI's (ed.5.0) Fast Line-of-sight Atmospheric Analysis of Hypercubes (FLAASH).Moreover, the high resolution panchromatic data were fused with the lower resolution multispectral data through a pan-sharpening procedure.The high pass filter (HPF) [50] and the hyperspherical color sharpening (HCS) [51] algorithms were employed.HPF is based on transferring high spatial details extracted from the panchromatic to the re-sampled multispectral image [52,53], while the HCS algorithm, which can handle any number of input bands, is based on the transformation from the native color space to the hyperspherical color space.In all of our experiments the pan-sharpened data from the HPF algorithm were used since the algorithm managed to enhance, slightly more effectively, the linear pattern of the vineyards.
In addition, the spectral readings, from the portable spectroradiometer, were taken with nadir view in the principal plane, while absolute bi-directional reflectance was obtained by normalizing the readings with the reference panel.More precisely, bi-directional reflectance was calculated by: where R T is the absolute spectral reflectance of the target, L T the radiance of the target, L R the relative irradiance value of the reference panel and k the panel's reflectance spectrum, as a calibrating factor, measured in the laboratory.The reflectance values for each variety's canopy were depicted through reflectance charts, within the range 0 to 1, forming the corresponding spectral signatures.

Methodology
Motivated by the successful application of object-based image analysis (OBIA) in very high resolution satellite data [3,4,9,47,48,54,55], the developed methodology was based on OBIA principles, while it was designed to address vineyard detection, vine canopy extraction and vine variety discrimination.In particular, for every processing step, several experiments were performed based on features employed from the literature on similar crop identification/detection studies (e.g., [8,9,16,[22][23][24]27,30]), while the optimal ones for all datasets were selected from a larger pool, through several experiments, feature analysis tools (assessing how each feature contributes to the discrimination task) and a trial and error procedure for fine tuning their parameters.To this end, Trimble's eCognition Developer (ed.8),MathWorks' MATLAB (2015b) and in-house developed software were employed.

Vineyard Detection
The first goal was to design and develop a generic classification framework for the detection of vineyards in very high resolution multispectral pan-sharpened data.A multi-scale classification procedure along with fuzzy logic and knowledge-based rule sets [47,54] were employed.Firstly, image objects were created on several scales using a multi-resolution segmentation algorithm, based on a region-merging procedure.By minimizing the average heterogeneity and maximizing their respective homogeneity, this algorithm consecutively merges pixels or existing image objects.The procedure was based on certain user-defined parameters, like scale, color/shape and smoothness/compactness, that were weighted together to define a homogeneity criterion [48,56].The classification scheme was based on two segmentation levels, i.e., Level 1 and Level 2 (Figure 2), while the core classification was conducted at Level 2, and then, the final classification maps were derived from the finer scale Level 1.
The class hierarchy was formed including at the top level two classes, i.e., Vegetation and Non-Vegetation areas (Figure 2).Then, the Vegetation parent class was split into two child classes: Vineyards and Other-Vegetation.For the class description, certain features and rules were combined and associated with each class.Membership functions were used based on fuzzy logic associated with the class description.These functions defined the relationship between feature values and the degree of membership to the classes, e.g., : "larger than", defined by left and right border values; "smaller than", defined by left and right border values; "full range" of values.Both at Level 1 and Level 2, the class Vegetation was associated with the Normalized Difference Vegetation Index (NDVI) [57] and the membership function "larger than".Class Non-Vegetation was associated with a logical rule in order to contain every object that was not classified as Vegetation.At this point, our goal was to separate scene's vineyards from the other vegetation types.This was accomplished at the much coarser Level 2 based on a combination of eight features (Table 1, left), which were able under a fruitful synergy to highlight the distinctive properties of the Vineyards 2 and Other-Vegetation 2 subclasses.A large set of features was considered, while through several experiments and feature analysis, the optimum for all datasets of features, rules and parameters was formed (Table 1 left; Figure 2).In particular, certain spectral features, like mean blue, mean red-edge, band ratios and indices, e.g., mean red-edge/mean blue, and the modified Soil-Adjusted Vegetation Index (MSAVI) [58], were employed.Moreover, similar to [9,24], we employed two textural features (i.e., GLCM dissimilarity and GLCM standard deviation), but on very high resolution pan-sharpened data (at 50 cm) instead of the significantly coarser ASTER and Landsat 8 data.During the last processing step, the classification result was refined at Level 1 (at a significantly finer segmentation scale) in order to approximate the correct object boundaries that were smoothed at the much coarser Level 2. A class-related feature (i.e., existence of super objects) was used to project class Vineyards 2 to class Vineyards 1, and class Other-Vegetation 1 was also described by the logical rule not Vineyards 1.Based on the developed class hierarchy and rule-set, the various image objects of all images were classified through a hierarchical classification procedure ensuring that, e.g., (i) classes that contain applicable child classes are not taking part in the classification of image objects; (ii) parent classes pass on their attributes to their child classes; and (iii) classes are only applied to classification of image objects, if all contained classifiers are applicable.

Vine Canopy Extraction
For any further analysis on the pan-sharpened data, the detection of vine canopy was required in order to separate the canopy from the material in-between the rows, like soil, weeds or other crop coverage.Another segmentation level was computed on the Vineyards class, creating a new, finer segmentation level over Level 1.The goal was to create small compact (with a size of about 2 to 3 pixels) image objects.
These image objects were, then, classified based on a nearest neighbor (NN) classifier.A membership value (between zero and one) was assigned according to the distance of object's feature space to its nearest neighbor.During the training stage, a similar number of image objects was given as samples for both the canopy and the other in-between the rows materials.Training samples were covering approximately 0.01% to 0.03% of the total area in every dataset.Different features were explored for the feature space of the classifier, while the final combination (Table 1, middle) was derived through several experiments in all datasets and feature analysis tools, indicating which features contribute more to the discrimination.The set of features that worked optimally in all datasets included both spectral (e.g., mean band intensity values, NDVI) and spatial features, e.g., mean difference (diff.) to neighbors and mean difference to brighter/darker neighbors.In particular, the spatial features describe the relations between image objects and their neighbor image objects in terms of mean band intensity values, e.g., : • Mean diff. to neighbors (red, NIR1): the mean band intensity value difference between an image object and its neighbor image objects (within a radius of 4 pixels) for bands red and NIR1.
• Mean diff. to brighter/darker neighbors (red, NIR1): the mean band intensity value difference between an image object and its brighter/darker direct neighbor (with a common border) image objects for bands red and NIR1.
• Number of brighter/darker objects (red, NIR1): the number of direct neighboring (with a common border) objects with higher/lower mean band intensity values for bands red and NIR1.
• Relative border to brighter objects (red, NIR1): the ratio of the shared image border of an image object with direct neighbor brighter objects to the total border, for bands red and NIR1.

Vine Variety Discrimination
For the vine variety discrimination, which was the most challenging task of this study, a two-step methodology was developed based on supervised object-based classification and on a majority voting post-classification process.In particular, during the first step, a supervised classification procedure using the nearest neighbor (NN) classifier was applied on the detected canopy objects from the previous step.Training image objects were given as samples based on the reference/ground truth data obtained from the field campaigns.For each dataset, the area of samples to the ground truth area per class/variety varied from 0.5% to 1%.
Based on intensive experiments on all datasets and the evaluation from the feature space analysis, a group of features was formed (Table 1, right) that exploited optimally the distinctive spectral properties of vine varieties across the VNIR spectrum, e.g.,: • Mean (coastal blue, red, red-edge, NIR1, NIR2): mean band intensity values of an image object for coastal blue, red, red-edge, NIR1, NIR2 bands.
• Ratio: the amount that a given image band contributes to the total brightness of an image object for all 8 bands.
• Mean NIR1/mean red: The ratio of mean band intensity values of band NIR1 to band red of an image object.
• Mean NIR2/mean blue: The ratio of mean band intensity values of band NIR2 to band blue of an image object.
In particular, as can be observed (Table 1, right), the majority of the spectral features incorporate the two NIR bands (of the multispectral WV-2 sensor) covering approximately the range of 760 nm to 1050 nm.This implied that these were the crucial wavelengths where the different varieties possess detectable spectral differences.The latter was also supported by the analysis of the in situ reflectance signatures (Section 4.2), which indicated that the studied varieties present discrete differences beyond 750 nm.
The second step of the variety discrimination task consisted of a post-classification procedure based on a majority voting process [59] towards associating a single variety with a certain parcel.In particular, based on the polygons of plots, a priori known from the ground truth data, each parcel was labeled according to the dominating variety, i.e., the variety in which the majority of pixels was classified within this parcel.
More than twenty different vine varieties and clones were present in these study areas and recorded in the corresponding reference/ground truth data.The ones that were more broadly cultivated or/and appeared in more than one study area were selected to participate in the experimental validation.Overall, thirteen vine varieties/clones were finally chosen to participate in the validation procedures.Experiments with three and up to six different varieties/clones were performed.
In addition, for a more throughout validation of the developed framework, we employed exactly the same supervised classification scheme initially applied on the pan-sharpened datasets (50-cm spatial resolution), on the atmospherically-corrected multispectral (2-m spatial resolution) ones.Our goal was to evaluate the effect and contribution of the pan-sharpening process and the canopy extraction procedure for the vine variety discrimination task.Generally speaking, the discrimination directly on the vine canopy is expected to result in lower classification errors by eliminating the influence of the in-between the vine rows materials.On the other hand, the pan-sharpening procedure distorts, at a greater or lesser degree, the initial reflectance values.For these classification experiments, a fine segmentation was computed on the multispectral images (scale parameter: 2; shape and compactness: 0.5), and samples were collected for each variety/clone from the same areas as in the pan-sharpened datasets.The standard NN classifier employed exactly the same features for the discrimination as in the classification of the pan-sharpened images, described in the previous paragraphs.

Experimental Results and Validation
In order to validate the developed framework, a comprehensive quantitative assessment was performed based on reference/ground truth data, which was produced using the available geospatial maps of varietal plantations and on data collected during multiple field campaigns, along with an intensive and laborious manual annotation procedure over all vine parcels in all study areas.It should be noted that the reference data, including the vector boundaries per vine parcel and the corresponding planted variety or clone, were available and/or created for all (100%) of the vine parcels studied.Regarding the vine canopy extraction, experimental results were validated based on manually-collected testing samples, since the annotation of every single vine row in all study areas was not practically feasible.Moreover, the proposed object-based framework was evaluated against a pixel-based support vector machine (SVM) classification scheme.In particular, recent studies targeting the classification of agricultural landscapes or different crop types [8,22,23,26,56], different tree species [60] and different vine varieties [45,46] were based on SVM classifiers with quite promising results.To this end, we employed a linear SVM classifier from the LIBSVM [61] library and performed similar experiments with exactly the same training samples as in the proposed object-based NN classification process.

Vineyard Detection and Canopy Extraction
During the first processing step of vineyard detection, the classification procedure classified the image objects into three different classes: (i) Vineyards; (ii) Other-Vegetation types; and (iii) Non-Vegetation.In Figure 3, the resulting classification maps along with the reference/ground truth data are presented for all study areas.The detected Vineyards are shown with green color (both in the ground truth and the resulted map), the Non-Vegetation with gray color and the Other-Vegetation types with yellow color.After a close look at Figure 3, one can observe that in all study, areas the developed framework managed to detect successfully the vast majority of vineyards.Errors included the misclassification of image objects belonging to the Other-Vegetation class as Vineyards, like in the upper left region in the Megaplatanos study area, and to a lesser extent, the misclassification of image objects belonging to Vineyards as Non-Vegetation, either due to the existence of very young vine plants (i.e., less or sparser canopy) or due to stress (e.g., water, disease), which reduces the reflectance of the green and near-infrared bands.A quantitative evaluation was performed based on the calculation of the standard metrics of completeness, correctness and overall quality [62,63] using the ground truth data, which included 100% of the vine parcels.In order to assess the accuracy metrics, the measures of true positives (TP), false positives (FP) and false negatives (FN) were calculated: area o f correctly − detected objects as Vineyards (TP) area o f the ground truth (TP + FN) Correctness = area o f correctly − detected objects as Vineyards (TP) area o f all objects detected as Vineyards (TP + FP) The quantitative evaluation indicated that the developed methodology resulted in high levels of accuracy in all study areas.More specifically, in all cases, the detection completeness and correctness was above 89% with an average at 93% and 92%, respectively (Table 2).During the next processing step, the canopy of the detected vineyards was extracted, based on the developed object-based procedure described in Section 3.2.In Figure 4, examples of the detected canopy are shown.In particular, the detected canopy is shown with a green color and the other materials (e.g., soil) in-between the rows with brown.After a closer look, one can observe that the employed features managed to separate the vine canopy from the other parts effectively, although the materials in between the rows presented significant spectral variations, as can be observed in the corresponding natural color composites.Certain parts are covered with just soil and present relatively high brightness values, while other regions present darker intensities and colors.The quantitative validation was performed based on a sampling set of objects as test area (TTA) masks.The TTA area varied from 0.01% to 0.02% of the total area classified.The calculated producer, user and overall accuracy (Table 3) indicated high rates of above 96% for all datasets.3 along with the corresponding ones from the proposed framework.The object-based framework resulted in higher accuracy rates regarding the average of all three metrics.In particular, OA reached 98.2% (object-based NN) compared to the 95.9% from the pixel-based SVM (Table 3).Generally speaking, however, one can observe that both classification schemes delivered high accuracy rates, i.e., the calculated producer, user and overall accuracy rates were above 88.0% in all cases.Although SVM can be consider a more efficient classifier, the object-based framework managed to deliver slightly better results mainly due to (i) the process at the object level, which can to a certain extent address the within-field variability; and (ii) the employed spatial features, which express the spectral relations between a specific image object and its neighbor ones.In particular, these spatial features were crucial, as they exploited the alternation between the foliage and the in-between the rows materials.

Vine Variety Discrimination
For the validation of the developed vine variety discrimination framework, several experiments were performed with attempts to discriminate 3, 4, 5 and 6 different varieties/clones per experiment.The same experiments were performed both on the atmospherically-corrected multispectral (2-m resolution) and the pan-sharpened (50-cm resolution, at the extracted canopy) data.A quantitative evaluation was performed based on the reference/ground truth (GT) data, which contained 100% of the vine parcels studied.The quantitative evaluation was calculated at the pixel level, based on the calculated confusion matrices between the classification maps and the GT.For an easier interpretation of the results, the number of classified pixels for each class/variety on the confusion matrices was expressed as a percentage on the GT pixels of the corresponding variety.Thus, the diagonal of the matrix expressed the producer accuracy (PA) for each variety.The evaluation of the second step of the discrimination experiments, after the post-classification majority voting process, was performed at the parcel level.
Table 4 presents the overall accuracy (OA) rates after the validation of the developed vine variety discrimination framework on both multispectral and pan-sharpened datasets.The quantitative evaluation indicated that the variety discrimination was performed with high OA rates at the parcel level (i.e., above 83%).In particular, the vast majority of parcels were associated with the correct variety in all of our experiments.It is also worth mentioning that the vine variety discrimination on the pan-sharpened (50-cm resolution) and the multispectral (2-m resolution) data resulted in similar OA rates, presenting both an average at 96%.At the pixel level, five out of 10 experiments indicated higher OA rates on the multispectral datasets, i.e., average OA of 61%, in contrast to an average OA of 56% for the pan-sharpened ones.Generally speaking, the important spectral correlation among different vine varieties/clones and, thus, the difficulty that the classifier confronts, either based on the detected canopy (pan-sharpened data) or the multispectral data, are apparent, since the higher OA classification rate reached a 75% maximum (Table 4, Trapeza study area, experiment with four varieties).The lower rate at the pixel level was 38% (Megaplatanos, experiment with six varieties on the pan-sharpened image), while as expected, the higher the number of competing varieties, the lower the classification rates in most cases.In the following sub-sections, these results are further analyzed and assessed for each study area based on the quantitative evaluation, the vine variety misclassification cases and the study of the spectral signatures of the in situ reflectance from the portable spectroradiometer.Trapeza study area: In the study area of Trapeza, six different vine varieties and clones (i.e., Merlot, Riesling, Sauvignon Blanc Clone I, Sauvignon Blanc Clone II, Syrah Clone I and Syrah Clone II) took part in the discrimination procedures.Various experiments with four, five and six varieties were performed.The quantitative evaluation on both datasets (pan-sharpened at 50 cm, multispectral at 2 m) resulted in the same quantitative OA rates (Table 4).In all cases, the correct vine variety was assigned to the corresponding parcels (i.e., 100% OA at the parcel level).
In Figure 5, experimental results after the application of the developed discrimination framework in the Trapeza study area are presented.The resulting classification maps based on the multispectral data (Figure 5, top left) and the ground truth data (Figure 5 The above observations can be also verified by the corresponding confusion matrices presented in Table 5.In particular, for the experiments on the multispectral data (Table 5, top) Sauvignon Blanc Clone II, Sauvignon Blanc Clone I, Syrah Clone I and Merlot presented higher PA rates (above 65%), while Riesling and Syrah Clone II were lower, i.e., below 50%.In a similar way, classification experiments on the detected canopy (pan-sharpened data, Table 5, bottom) indicated the same OA performance, with Sauvignon Blanc (both clones) and Merlot achieving the highest PA rates.It should be noted that after the per-parcel majority voting, the discrimination procedure managed to associate all vine parcels with the correct varieties, even for the varieties that presented low classification rates, e.g., Riesling and Syrah Clone II.To sum up, based on all of the experiments performed in the Trapeza study area, Merlot and both Sauvignon Blanc clones presented a distinguishable behavior with high PA rates.It is also important to mention that in between the two different clones of Sauvignon Blanc and Syrah, relatively low classification errors (<15%) were reported.
The aforementioned results were compared to the spectral signatures of the in situ reflectance from the portable spectroradiometer.All varieties presented discrete differences beyond 750 nm (Figure 6, top left).In particular, the average reflectance indicated that the canopy of the Merlot vines differed from all others, since it presented the highest reflectance values across the spectrum.This particular variety presented high PA rates on the performed classification experiments as well (Table 5).The Sauvignon Blanc I and the Syrah I varieties presented relatively lower reflectance levels.Riesling also presented discrete differences compared to the other varieties and the lowest values at the NIR.However, on the image classification experiments, Riesling resulted in having generally low PA rates.

Megaplatanos study area:
In the Megaplatanos study area, six different vine varieties, namely Sauvignon Blanc, Syrah, Chardonnay, Merlot, Robola and Cabernet Sauvignon, took part in the experiments.Four, five and six different varieties were employed under the same discrimination procedure on both datasets.In general, the PA rates for each variety were higher for the classification on the multispectral image.As in the Trapeza study area, Merlot was among the most separable varieties, presenting the highest PA rates at the pixel level, of above 86% in the multispectral dataset and over 75% in the pan-sharpened data.
A typical example is presented in Figure 7 (left) for a Merlot parcel, where most of the image objects were correctly classified as class Merlot during an experiment with six competing varieties on the pan-sharpened image.Moreover, in this study area, Syrah and Sauvignon Blanc were, also, classified with relatively high PA rates, while both of them presented misclassifications mainly as Cabernet Sauvignon with rates from 15% to 30%.Chardonnay vines were mainly confused with the Cabernet Sauvignon, the Sauvignon Blanc and the Robola varieties.Among all, Robola presented the lowest PA rates (<45%) in all performed experiments.The discrimination rates (OA > 92%) for the parcel-based evaluation indicated that the vast majority of parcels was correctly labeled, apart from few Robola and Cabernet Sauvignon parcels.
The in situ reflectance data indicated that the Robola vine canopy possessed the lowest reflectance levels compared to the other varieties in both the visible and the NIR spectrum (Figure 6, top right).This characteristic, apparently, was not exploited by the classifier or the proximate observations over the Robola canopy could not simulate their spatial and spectral characteristics as they were recorded by the satellite sensor.Moreover, Syrah presented the highest reflectance values, while Sauvignon Blanc, Cabernet Sauvignon and Chardonnay presented similar reflectance values, which was also reported by their misclassification rates during the classification experiments.The spectral signature of Merlot followed a relatively different pattern across the spectrum and, in particular, in regions between 380 to 500 nm and 800 to 1000 nm.This was also in accordance with the results from the classification experiments, as Merlot resulted in being the most distinguishable variety during the discrimination procedures.
Naoussa study area: In the study area of Naoussa, four different vine varieties (i.e., Syrah, Merlot, Cabernet Sauvignon and Xinomavro) producing red wine took part in the experiments towards the discrimination of three or four different vine varieties each time.The quantitative evaluation (Table 4) indicated that the discrimination experiments in this study area resulted in having lower OA rates compared to the other three.The classification evaluation resulted in OA rates from 50% to 59%, while the discrimination task at the parcel level reached OA rates from 83% to 91%.In general several misclassification errors were reported in all experiments.The PA rates for the variety classification were slightly higher for the pan-sharpened data.Again, the Merlot variety achieved the highest PA rates compared to the other varieties.A typical example is presented in Figure 7 (right) with the resulting classification map calculated on the detected canopy (pan-sharpened data) during an experiment with four competing varieties.Here, most objects have been correctly classified as Merlot with a relatively small number of errors.In all experiments performed, Xinomavro reported several misclassification errors, mixing with all of the other varieties, with rates of 11% to 26% at the pixel level.Cabernet Sauvignon and Syrah appeared correlated (with mixing rates of about 20%), while both varieties presented misclassification errors to Merlot (20% to 25%).Due to these increased misclassification errors, the final result after the majority voting processing step presented some wrongly labeled parcels, concerning mainly the Cabernet Sauvignon and the Xinomavro varieties.However, the majority of parcels were correctly labeled in all cases, and the OA rates at parcel level stayed above 83%.
In the Naoussa study area, the in situ canopy reflectance was obtained for the Merlot, Xinomavro and Cabernet Sauvignon varieties.In general, the reflectance levels were relatively lower than those in the other three study areas (Figure 6, bottom left).In particular, the spectral signatures indicated only slight differences mainly in the NIR region (after 740 nm).Cabernet Sauvignon presented slightly higher reflectance levels, while Merlot and Xinomavro presented high similarity in their reflectance values.
Amynteo study area: In the Amynteo study area, four vine varieties and clones, namely Tannat, Montepulciano, Syrah Clone I and Syrah Clone II, took part in the discrimination experiments.Classification procedures were employed for the discrimination between three and four varieties per time.Although, the OA from the classification step presented relatively low rates, the final discrimination at the parcel level was 100% accurate in all cases.
In particular, between the two Syrah clones, the confused by the classifier pixels were relatively low, with about 15% being confused on the multispectral dataset and 22% on the pan-sharpened one.Syrah Clone I presented mixings, of rates between 22% and 36%, with Montepulciano.During the classification experiment between four competing varieties, Tannat was significantly confused (42%) with Syrah Clone II on the multispectral data, but still class Tannat gained the majority of the pixels with a 43% percentage.In a similar way, Tannat was confused with Syrah Clone II and Syrah Clone I, on the pan-sharpened data, but with relatively lower rates at 28% and 18%, respectively.In Figure 8, the resulting classification map is shown for an experiment on the multispectral dataset with four different varieties/clones in the Amynteo study area.After a closer look, one can observe that all parcels present many classification errors, while about half of Tannat parcel objects appear to be classified as Syrah Clone II.The analysis of the in situ reflectance indicated that all varieties presented discrete differences beyond 750 nm (Figure 6, bottom right).Tannat was the variety with the highest reflectance levels followed by Syrah and Montepulciano.Although for Tannat, the average in situ reflectance in both the visible and NIR spectrum stayed consistently at higher levels, differing from the two other, during the classification procedures class Tannat resulted in numerous misclassification errors.As in the previous case studies, the classifier and employed features did not manage to model this difference or the canopy in situ reflectance for Tannat could not be observed at the particular spatial scale of WV-2.Syrah and Montepulciano possessed similar reflectance values before 730 nm, while in the classification procedures, Syrah Clone I presented many misclassification errors with Montepulciano.

Variety Discrimination Using a Pixel-Based Linear SVM
The same experiments were also conducted using a pixel-based linear SVM classifier and exactly the same training samples.The features used for the pixel-based SVM framework included the eight WV-2 spectral bands, the NDVI and four band ratios, i.e., NIR1/red, NIR2/red, NIR1/blue, NIR2/blue, that had been also employed in the object-based NN framework.Quantitative results from the performed experiments are presented in Table 6 comparing the outcome from both the proposed object-based NN and pixel-based SVM framework.As has been already mentioned, the quantitative evaluation was performed based on reference data, including 100% of the vine parcels studied.In a similar way as in Table 4, the evaluation was performed at the pixel level (after the NN or the SVM classifier) and at the parcel level (after the majority voting process).Results indicated that the object-based framework outperformed the pixel-based SVM one, both at the pixel and the parcel level.In particular, after majority voting, the proposed object-based approach managed to score above 90% regarding the average OA rates.The difference in the calculated average OA at the pixel level was lower than 10%.These results are in accordance with similar studies that report on the successful application of object-based frameworks for crop type [9,24,26,27] or vegetation/tree species [1,3,4] classification.Classifying (even with a relatively small size of foliage segments) canopy objects instead of pixels can provide a more compact representation regarding the distinctive properties between different varieties/classes that a classifier searches for decision making.Thus, misclassification errors at the pixel level potentially caused due to the within-field variability can be tackled to a certain extent through object-based approaches.

Discussion
Regarding the vineyard detection task, results reached high completeness (average at 93%) and correctness (average at 92%) rates, based on the experiments from all datasets in different viticultural zones and terroirs (covering more than 1200 ha) with more than 20 vine varieties/clones and reference/ground truth data containing 100% of all vineyards in the study areas.These quantitative results are similar or better than those that have been reported in the literature [8,9,23,24,28,30] for vineyards or similar crop detection studies.It should be also noted that in the EU (in particular, in Greece), the terroirs are significantly smaller in size, usually covering non-flat viticulture regions with highly complex terrain, and thus, their detection is significantly more challenging than, e.g., in Australia, California, Chile.Moreover, the vine planting density worldwide is pretty much standard based mainly on the commercial agriculture machinery/tractors, and therefore, the textural features that map their pattern at this particular spatial resolution are expected to work equally effectively in different terroirs.
Regarding the spatial resolution of the employed satellite data, the quantitative results previously mentioned in Section 4 highlighted certain points for further discussion.Firstly, for the vineyard detection task, the pan-sharpening processing step (delivering a 50-cm spatial resolution) was absolutely required since, along with the spectral, the textural features managed to adequately describe the linear planting pattern, resulting in high detection rates.Therefore, the spatial resolution of the data during the vineyard detection step is critical.However, the pan-sharpening procedure affects the reflectance values of the initial lower resolution (2 m) multispectral bands due mainly to the applied bilinear interpolation and employed edge maps.
The performed quantitative evaluation (Table 4) showed that the pan-sharpening procedure did not impede the variety discrimination OA rates, implying that the radiometric shift (due to the pansharpening) did not alter or smooth the distinctive properties of the initial 2D signal (reflectance values of the multispectral bands).In particular, although the evaluation at the pixel level (after the supervised classification) resulted in slightly higher OA rates at the coarser (2-m spatial resolution) multispectral datasets (61% the average OA; Table 4), the final discrimination results at the parcel level (after the majority voting) delivered similar quantitative rates on both datasets (96% the average OA; Table 4).Therefore, one can claim that since the pan-sharpening significantly contributed to the vineyard detection task and did not have a negative impact on the discrimination procedure, then it should be included in the processing chain of the same or similar applications.
Moreover, the comparison of the proposed object-based framework with the pixel-based SVM one indicated that for both tasks of canopy extraction and variety discrimination, the classification of objects (even with a relatively small size) delivered higher accuracy rates during experiments using similar features and the same training and testing samples.This is, also, in accordance with similar studies [3,4,9,24,26,27] suggesting that one would expect lower misclassification errors when classifying foliage segments than single pixels with high and very high resolution satellite data.
Furthermore, regarding the vine variety discrimination, the developed framework delivered similar or better OA rates compared to similar studies, which employed airborne hyperspectral data and conducted experiments towards the separation of two competing varieties [44][45][46].One should, furthermore, note the low percentage of samples that were required during the canopy extraction (less than 0.03% of all objects to be classified) and the variety discrimination (less than 1% of the GT for each variety) procedures.This difference appears normal considering the difficulty of each problem.Moreover, for the samples' selection, we did experiments by randomly selecting different groups of samples during the training procedure, both for canopy extraction and variety discrimination.However, the overall accuracy rates were similar, and their difference was lower than 1% at the pixel level.
It is worth mentioning here that the experimental results are based on satellite data acquired around the veraison period, and therefore, the performance of the proposed framework in the beginning of the growth cycle or near its end has not been evaluated.However, veraison is considered the optimal period for data acquisition when addressing such vine/variety detection tasks, since the plant reaches a stable state regarding the development of the foliage/canopy, spending the next period on fruit ripening.
Last, but not least, through the analysis of the quantitative results from the variety classification, certain aspects were highlighted regarding the spectral behavior of the studied varieties.Merlot was the variety that achieved the highest PA rates (Figure 9) in the classification experiments.Sauvignon Blanc presented in most cases a distinguishable behavior and scored in average PA rates above 65% at the pixel level (Figure 9).Syrah generally resulted in medium PA rates and, in most cases, was discriminated from the other varieties.Moreover, the two studied clones (I and II) of Syrah were not much confused between each other, as the experiments in both Trapeza and Amynteo indicated.Furthermore, Cabernet Sauvignon presented differences (>10%) in PA rates at the pixel level between the coarser multispectral and the pan-sharpened data (Figure 9).Robola and Riesling were the varieties that in all our experiments resulted in lower PA rates.In general, classification errors largely depend on the specific varieties that are competing during the classification.If these varieties do not possess any distinctive features across the spectrum, then even under a binary classification problem, the errors will be significant.

Conclusions
In this study, an object-based framework was developed that is able to detect vineyards, extract vine canopy and discriminate several varieties from very high resolution satellite data.For the validation of the developed framework, intensive experiments were carried out, employing WorldView-2 data over four viticulture regions in Greece, along with concurrent in situ canopy reflectance observations from a portable spectroradiometer.The performed quantitative evaluation indicated that the developed approach managed in all cases to detect vineyards with high completeness and correctness rates (above 89%) on the pan-sharpened imagery.The vine canopy extraction methodology was validated with OA rates of above 96%.The various experiments for the discrimination of vine varieties, i.e., discriminating between three, four, five and six competing varieties/clones, indicated that at the parcel level the OA rates were in all study areas above 85% (average at 96%; Table 4).The comparison between the proposed object-based framework and a pixel-based one, applying a linear SVM classifier, indicated that both for the canopy extraction and the variety discrimination, the object-based approach resulted in lower misclassification errors.Last, but not least, the observations concerning the spectral characteristics and behavior of the different varieties during the classification experiments on the satellite data were in most cases in agreement with the analysis of the in situ reflectance data and the derived spectral signatures.In particular, this combined study indicated that certain vine varieties, like Merlot and Sauvignon Blanc, possess distinct spectral properties and distinguishable behavior.The overall validation and quite promising experimental results indicate that a thorough per-satellite sensor sensitivity analysis can set the basis for efficient operational tools for variety-based data analysis in viticulture.Among the future work are the integration of more efficient classifiers, experiments with more varieties and the development of a fully-automated processing framework.

Figure 1 .
Figure 1.Data acquisition campaigns were conducted in four different viticulture regions in Greece, in Protected Geographical Indication ( PGI) zones.The campaigns were scheduled during the veraison period, conducted in 2012 and 2013, from late July to early August.

Figure 2 .
Figure 2. The classification features and rules of the developed vineyard detection framework.

Figure 3 .
Figure 3. Experimental results after the application of the developed vineyard detection methodology.The detected vineyards (with green color) and the ground truth are shown for all study areas.The Other-Vegetation class is presented with yellow color.

Figure 4 .
Figure 4.The extracted vine canopy with green color and the detected in-between the rows materials, e.g., soil, with brown on two vine parcels.The corresponding natural color composites are also presented (right).

Figure 5 .
Figure 5. Experimental results from the supervised classification in the Trapeza study area for the discrimination of six different varieties.The overall classification accuracy at the pixel level (first processing step) was just above 60%, while after the majority voting (second processing step), the OA reached 100% at the parcel level (Table 4).(a) The resulting classification map based on the multispectral data (2-m resolution, left) and the ground truth data (right); (b) the resulting classification map based on the detected canopy (pan-sharpened data at 50-cm resolution) for vine parcels with Sauvignon Blanc Clone I (left) and Sauvignon Blanc Clone II (right).

Figure 6 .
Figure 6.The in situ reflectance spectral signatures of different varieties for the four study areas.

Figure 7 .
Figure 7.Typical examples of classification results based on the detected canopy (50-cm spatial resolution) over two Merlot vine parcels.

Figure 8 .
Figure 8.The classification result (left) in an experiment on the multispectral dataset (2-m resolution) with four varieties, in the Amynteo study area, and the corresponding ground truth data (right).

Figure 9 .
Figure 9.The average PA rates, at the pixel level, for the four most commonly-used varieties in the discriminative classification experiments.Quantitative results on the multispectral (2-m resolution) datasets are given with brown color, while the ones for the pan-sharpened (50-cm resolution) with green.

Table 1 .
The spectral, spatial and textural features that were associated with the developed rule-set hierarchy for the (i) vineyard detection; (ii) vine canopy extraction and (iii) vine variety discrimination tasks.

Table 2 .
The developed framework managed in all cases to detect vineyards with high detection rates; completeness and correctness rates stayed above 89% in all study areas.

Table 3 .
Comparing results obtained from the proposed object-based NN classification framework and the pixel-based SVM classifier for vine canopy extraction.

Table 4 .
Quantitative evaluation results after the application of the developed vine variety discrimination framework.

Table 5 .
The resulting confusion matrices for the vine variety discrimination between six different vine varieties/clones in the Trapeza study area.

Table 6 .
Comparing results from the proposed object-based NN framework with a pixel-based SVM classifier for the variety discrimination task.