Estimating Leaf Area Index (LAI) in Vineyards Using the PocketLAI Smart-App

Estimating leaf area index (LAI) of Vitis vinifera using indirect methods involves some critical issues, related to its discontinuous and non-homogeneous canopy. This study evaluates the smart app PocketLAI and hemispherical photography in vineyards against destructive LAI measurements. Data were collected during six surveys in an experimental site characterized by a high level of heterogeneity among plants, allowing us to explore a wide range of LAI values. During the last survey, the possibility to combine remote sensing data and in-situ PocketLAI estimates (smart scouting) was evaluated. Results showed a good agreement between PocketLAI data and direct measurements, especially for LAI ranging from 0.13 to 1.41 (R2 = 0.94, RRMSE = 17.27%), whereas the accuracy decreased when an outlying value (vineyard LAI = 2.84) was included (R2 = 0.77, RRMSE = 43.00%), due to the saturation effect in case of very dense canopies arising from lack of green pruning. The hemispherical photography showed very high values of R2, even in presence of the outlying value (R2 = 0.94), although it showed a marked and quite constant overestimation error (RRMSE = 99.46%), suggesting the need to introduce a correction factor specific for vineyards. During the smart scouting, PocketLAI showed its reliability to monitor the spatial-temporal variability of vine vigor in cordon-trained systems, and showed a potential for a wide range of applications, also in combination with remote sensing.


Introduction
Leaf area index (LAI, total one-sided area of leaf tissue per unit ground surface; [1]) is a quantitative descriptor of canopy density and it is important index for monitoring vine (Vitis vinifera L.) vigor and for providing precious information to support vineyard management. Leaf area density is significantly correlated with radiation interception and light environment in the fruit zone [2], in turn affecting the canopy microclimate and the plant carbon balance, as well as yield and grape composition [3,4]. For these reasons, attempts have been made to forecast vine yields using remotely sensed vegetation indices (i.e., normalized difference vegetation index, NDVI) related with plant vigor and LAI [5]. Concerning management support, the relationship between LAI and crop coefficient (Kc) was used to quantify water requirements and define irrigation strategies [6].
distribution of infinitely small leaves) behind its simplified transmittance model [31], for a variety of broad-leaf tree canopies [32] and for seasonal monitoring purposes [33].
The availability of a smart-app-characterized by low-cost, able to provide LAI estimates in absolute terms, in real-time, without requiring time consuming post-processing or the definition of parameters/settings-could represent an useful and user-friendly tool for operational purposes in the viticulture sector, suitable for the agricultural stakeholders (e.g., producers, insurance agents, etc.).
Recently, another smart-app (VitiCanopy) [34] was developed to estimate vine canopy vigor and porosity. VitiCanopy is based on the acquisition of images from below the canopy similarly to DHP, and on their automatic processing, analogously to PocketLAI. However, the VitiCanopy user needs to specify the gap fraction threshold, meaning the percentage of image pixel that corresponds to sky, before taking the measurements. VitiCanopy showed good agreement with LAI-2000 [34], although the authors underlined the influence of the cordon included in the image on the app accuracy. The resulting overestimation depends on intrinsic factors concerning the cordon size (e.g., vine age and training system) and on aspects related to the measurement protocol (i.e., the distance between the camera and the cordon). The authors highlighted that the assessment of the contribution of the cordon is not necessary when the aim is to monitor relative differences of vigour in time and space. On the other hand, when absolute LAI values are required, the authors suggested the calibration of a case-specific correction factor using measurements carried out on bare cordons during dormancy to quantify the cordon contribution. In this respect, the automatic acquisition of images at 57.5 • implemented in PocketLAI should make the latter free from any need of parameterization or calibration, since images are captured pointing the camera toward the leaves wall and thus excluding the cordon.
The aims of this study were to assess the reliability of vineyards LAI values retrieved with PocketLAI and to develop a specific measurement protocol for grapevine. The performances of PocketLAI were assessed against destructive measurements and compared with those of DHP. Moreover, the possibility of using in-situ PocketLAI measurements and Normalized Difference Vegetation Index (NDVI) maps for obtaining spatially distributed LAI maps was also explored via smart scouting.

Survey Scheme
The study was carried out during 2015, in the rainfed experimental vineyard of Riccagioia s.c.p.a. (Pavia Province, Italy; Lat. 44.98 N, Long. 9.09 E, 138 m a.s.l.). Vines were planted in 2008 with a row and vine spacing respectively of 2.4 m and 1.0 m, cultivated using the Guyot system and unilateral cordon-pruning. The vineyard covers an area of about 1 ha and it is characterized by a pronounced heterogeneity, mainly due to the use of different cultivars (Vitis vinifera: cv. Barbera and Chardonnay) and rootstocks (hybrids of Vitis berlandieri and Vitis riparia) within each row, as well as to different vine care practices (with or without green pruning).
Sampling areas of different vigor were identified during six surveys carried out from the emergence of inflorescence to the beginning of ripening, in order to monitor the evolution of vine canopy during the crop cycle and provide a wide range of LAI values on which testing the indirect methods against the destructive measurements. In particular, the surveys were performed at the following growth stages of the BBCH phenological scale: 56, 57, 61, 74, 77, 81 [35].
During each of the first five surveys, three sampling areas of different vigor were identified, and three vines were sampled for each area. The sampling areas were characterized by cv. Barbera vines, all subjected to green pruning. The differences of canopy vigor in-space were mainly due to the use of three different rootstock hybrids (Table 1), and in-time to the different surveyed growth stages.
During the sixth survey, five sampling areas of different vigor were identified on the basis of the relative differences of canopy vigor described by a high resolution NDVI map, and-in each of these areas-two vines were sampled. The NDVI map was derived using a sensor placed on a four-wheeler quad-bike that detected the spectral response of the vines while moving through the rows. This last survey had the aim to extend the range of LAI values on which test the performance of indirect methods and to explore the possibility to combine PocketLAI with remote sensing data. In this case, the spatial variability of canopy vigor was due to the presence of different combination of cultivar × rootstocks (Table 1), and to the involvement of vineyard areas not subjected to green pruning. On the sampled vines, PocketLAI and DHP estimates were first taken, followed by the collection of destructive samples obtained by stripping all leaves from the canopies. Moreover, during the first five surveys, the number of shoots per plant was determined, and the average number of leaves per shoot was estimated by counting the leaves on a sample of five shoots. These two measures were summarized in a descriptive index (theoretical leaf number; TL) given by the shoot number multiplied by the average number of leaves per shoot.

Indirect LAI Estimates
The method implemented in PocketLAI and the app functioning were described in previous papers [28,31,32] (documentation available at the http://www.cassandralab.com/mobiles/1 and licensing at info@cassandralab.com). Using the accelerometer and the device camera, PocketLAI automatically takes images from below the canopy at a view angle of 57.5 • while the user is rotating the device along its main axis. The gap fraction is derived using a fully automatic segmentation algorithm expressly developed to detect the sky pixels according to their chromatic values in a Hue-Saturation-Brightness (HSB) color space. The LAI value was retrieved according to the light transmittance model proposed by Warren-Wilson [29] and recently discussed by Baret et al. [30].
For this study, PocketLAI was installed on an inexpensive Samsung GT-i9105 Galaxy S II Plus smatphone. In the specific case of vineyard, measurements were taken by positioning the device at about 15 cm below the canopy and at a distance of 0.4 m from the row, orienting the device toward the vines (Figure 1a). This allowed excluding the empty space between the bottom of the canopy and the soil surface (Figure 1b), as well as the sky above the canopy (Figure 1c). For each plant, four measurement replicates were acquired by processing two shots per row side while moving along the row at about 20 and 40 cm from the vineyard pole ( Figure 1d). Averaging the four replicates, the LAI of each vine was obtained. Then, averaging the LAI of the sampled vines, the mean plant LAI (LAIp) in each sampling area was calculated.
of each vine was obtained. Then, averaging the LAI of the sampled vines, the mean plant LAI (LAIp) in each sampling area was calculated. The LAIp provided by PocketLAI was converted in vineyard LAI (LAIv) using Equation (1): where D (m) is the horizontal distance of the device from the cordon, and R (m) is the row-spacing. The equation was derived from the assumption that the LAIp obtained from PocketLAI is given by the ratio between the plant leaf area and the unit ground surface considered by the smart-app (Equation (2)). The latter is the vine-spacing multiplied by the inter-row width considered by the PocketLAI on both sides of the row, corresponding to the double of the distance between device and row (0.8 m).
The value of leaf area per plant assessed by PocketLAI was retrieved (Equation (3)) and referred to the unit of ground surface provided for single vine in the vineyard, dividing by vine-spacing × row-spacing [16] (Equation (1)). The result is the conversion of LAIp in LAIv.
where LA (m 2 ) is the plant leaf area, V (m) is the vine-spacing, D (m) is the horizontal distance of the device from the cordon. The hemispherical images were collected through fisheye lens, using a protocol similar to that used for PocketLAI. The camera was positioned below the canopy along the row line, with the lens oriented upward, acquiring two images from each plant, moving horizontally at 20 and 40 cm from the pole. The set of fisheye images from each sampling area was processed using the CAN-EYE software (v 6.314; http:www.avignon.inra.fr/can_eye) [15], since it demonstrated to be one of the best packages in a comparative study performed by Liu et al. [36]. The segmentation of the images targeted the discrimination of green organs from the sky background and from the brown pixels, in order to exclude from the LAI estimate the contribution of the cordon and other minor woody components. The 'true LAI' values provided by the software were used as LAIv. The 'true LAI' is defined as the LAI value assessed by CAN-EYE taking into account the computation of average leaf inclination angle and clumping parameter. The LAIp provided by PocketLAI was converted in vineyard LAI (LAIv) using Equation (1): where D (m) is the horizontal distance of the device from the cordon, and R (m) is the row-spacing. The equation was derived from the assumption that the LAIp obtained from PocketLAI is given by the ratio between the plant leaf area and the unit ground surface considered by the smart-app (Equation (2)). The latter is the vine-spacing multiplied by the inter-row width considered by the PocketLAI on both sides of the row, corresponding to the double of the distance between device and row (0.8 m).
The value of leaf area per plant assessed by PocketLAI was retrieved (Equation (3)) and referred to the unit of ground surface provided for single vine in the vineyard, dividing by vine-spacing × row-spacing [16] (Equation (1)). The result is the conversion of LAIp in LAIv.
where LA (m 2 ) is the plant leaf area, V (m) is the vine-spacing, D (m) is the horizontal distance of the device from the cordon. The hemispherical images were collected through fisheye lens, using a protocol similar to that used for PocketLAI. The camera was positioned below the canopy along the row line, with the lens oriented upward, acquiring two images from each plant, moving horizontally at 20 and 40 cm from the pole. The set of fisheye images from each sampling area was processed using the CAN-EYE software (v 6.314; http:www.avignon.inra.fr/can_eye) [15], since it demonstrated to be one of the best packages in a comparative study performed by Liu et al. [36]. The segmentation of the images targeted the discrimination of green organs from the sky background and from the brown pixels, in order to exclude from the LAI estimate the contribution of the cordon and other minor woody components. The 'true LAI' values provided by the software were used as LAIv. The 'true LAI' is defined as the LAI value assessed by CAN-EYE taking into account the computation of average leaf inclination angle and clumping parameter.

Destructive LAI Measurement
The LAI estimates from destructive sampling followed the approach suggested by Breda et al. [37] and Johnson and Pierce [16]. Destructive LAI measures were collected immediately upon completion of indirect estimates by removing all the leaves at the petiole from each sample plant. Then, for each plant, the leaf fresh weight (W) was determined and a disc of known area (0.01 m 2 ) was obtained for each of the leaves in a 40-leaf subsample including leaves of different ages. The 40 discs per plant were weighed and the specific leaf area (SLA; m 2 ·g −1 f.w.) was computed. The obtained values of SLA, W and plant density were used to derive the vineyard LAI (LAIv) according to Equation (4): Therefore, similarly to Equation (1), in order to compute LAIv, the leaf area per plant was referred to the unit of ground surface provided for single vine in the vineyard [16].

Data Analysis
For each sampling area, the agreement between mean LAIv values determined using destructive and indirect methods was quantified using the following metrics: relative root mean square error (RRMSE; 0 to +∞, optimum 0) [38], mean absolute error (MAE; 0 to +∞, optimum 0) [39], modelling efficiency (EF; −∞ to +1, optimum +1) [40], coefficient of residual mass (CRM; from −∞ to +∞ , optimum 0; if positive means underestimation and vice versa) [41]. Moreover, parameters of the linear regression equation between LAIv values estimated with direct and indirect methods were calculated, as well as between LAIv values and TL data. Moreover, for the sixth survey, a relationship between PocketLAI estimates and NDVI data was computed and used to derive spatially distributed LAI values for the entire vineyard.

Performance of Indirect Methods
The LAIv values obtained from destructive measurements during the campaign ranged from 0.13 to 2.84, in agreement with those normally reported in commercial vineyards [16,17,25]. The results of the comparison between LAIv values from PocketLAI, DHP and destructive measurements are shown in Table 2 and Figure 2.  Considering the whole dataset of surveys ( Table 2, Dataset-1; Figure 2a), PocketLAI values were in agreement with those from destructive measurements, with a highly significant correlation (p < 0.001) (R 2 = 0.77) and practically without under-or over-estimating behaviors (CRM = 0.06). However, the comparison was largely affected (RRMSE = 43%) by a single value (LAIv = 2.84), corresponding to the sampling area characterized by the highest vigor, identified via smart scouting before the last survey using the NDVI map. The highest LAI value corresponds to the sampled area of the vineyard untreated by the green pruning, and therefore it can be considered an outlying data.
The exclusion of this single value from the dataset led to greatly improve the coherence between PocketLAI and destructive LAIv values (Table 2, Dataset-2; Figure 2c). The new values of the agreement metrics decidedly approximated their optimum values: R 2 was 0.94, RRMSE reached 17%, and EF was 0.93.
Considering the whole dataset ( Table 2, Dataset-1; Figure 2b), DHP showed a systematic overestimation (CRM = −0.87), with high RRMSE (99.46%) and negative value of EF (meaning that the average of destructive measurements is a better estimator of each LAI value compared to the DHP values). On the other hand, DHP showed high correlation with destructive LAIv values (R 2 = 0.94) and its performances were not affected by the outlier value (i.e., no improvements detected repeating the analysis with the exclusion of the outlying value; Table 2, Dataset-2; Figure 2d).
Both indirect methods showed highly significant correlations between LAIv estimates and the descriptor of plant vigour represented by TL (Table 2), with values of R 2 similar to that achieved using the destructive method. However, a slightly better agreement between TL and LAIv was obtained for PocketLAI (R 2 = 0.96, it was 0.85 for DHP). Figure 3 shows the increasing trend of the mean LAIv during the growing season, from the beginning of May to the end of June (i.e., during the first five surveys). Data confirmed what Considering the whole dataset of surveys (Table 2, Dataset-1; Figure 2a), PocketLAI values were in agreement with those from destructive measurements, with a highly significant correlation (p < 0.001) (R 2 = 0.77) and practically without under-or over-estimating behaviors (CRM = 0.06). However, the comparison was largely affected (RRMSE = 43%) by a single value (LAIv = 2.84), corresponding to the sampling area characterized by the highest vigor, identified via smart scouting before the last survey using the NDVI map. The highest LAI value corresponds to the sampled area of the vineyard untreated by the green pruning, and therefore it can be considered an outlying data.
The exclusion of this single value from the dataset led to greatly improve the coherence between PocketLAI and destructive LAIv values (Table 2, Dataset-2; Figure 2c). The new values of the agreement metrics decidedly approximated their optimum values: R 2 was 0.94, RRMSE reached 17%, and EF was 0.93.
Considering the whole dataset ( Table 2, Dataset-1; Figure 2b), DHP showed a systematic over-estimation (CRM = −0.87), with high RRMSE (99.46%) and negative value of EF (meaning that the average of destructive measurements is a better estimator of each LAI value compared to the DHP values). On the other hand, DHP showed high correlation with destructive LAIv values (R 2 = 0.94) and its performances were not affected by the outlier value (i.e., no improvements detected repeating the analysis with the exclusion of the outlying value; Table 2, Dataset-2; Figure 2d).
Both indirect methods showed highly significant correlations between LAIv estimates and the descriptor of plant vigour represented by TL (Table 2), with values of R 2 similar to that achieved using the destructive method. However, a slightly better agreement between TL and LAIv was obtained for PocketLAI (R 2 = 0.96, it was 0.85 for DHP). Figure 3 shows the increasing trend of the mean LAIv during the growing season, from the beginning of May to the end of June (i.e., during the first five surveys). Data confirmed what discussed above, with the constant over-estimation of LAIv provided by the DHP, whereas PocketLAI estimates were always very close with destructive measurements. Moreover, at the second survey, the DHP detected an average LAI decrease of 0.19 not confirmed by the other measurement methods. The inconsistency was due to the sky lighting during the acquisition of hemispherical images. During the second survey, the conditions (sunny day and survey carried out with sun at the zenith) led to high level of sunlight affecting the evenness of foliage lighting, with resulting uncertainty in image segmentation. This led to underestimation errors, as already observed by other authors [26,42]. The PocketLAI view angle of 57.5 • minimized the effect of those conditions because, contrarily to hemispherical photography, the device camera was not oriented directly toward the zenith during image acquisition.
Sensors 2016, 16,2004 8 of 12 discussed above, with the constant over-estimation of LAIv provided by the DHP, whereas PocketLAI estimates were always very close with destructive measurements. Moreover, at the second survey, the DHP detected an average LAI decrease of 0.19 not confirmed by the other measurement methods. The inconsistency was due to the sky lighting during the acquisition of hemispherical images. During the second survey, the conditions (sunny day and survey carried out with sun at the zenith) led to high level of sunlight affecting the evenness of foliage lighting, with resulting uncertainty in image segmentation. This led to underestimation errors, as already observed by other authors [26,42]. The PocketLAI view angle of 57.5° minimized the effect of those conditions because, contrarily to hemispherical photography, the device camera was not oriented directly toward the zenith during image acquisition.

PocketLAI and NDVI
The pixels of the georeferenced NDVI map of the vineyard were clustered into five classes using the ArcGIS v.10 package, and PocketLAI estimates were collected at the sixth survey in sampling area representative of each class (Figure 4a). A regression equation (R 2 = 0.90) was then derived between the LAIv obtained from PocketLAI and NDVI values in the sampling areas (indicated by black circles in Figure 4a). The regression equation was then used to convert all the other NDVI pixels of the map into LAIv values, using the spatial analysis tool of ArcGIS for raster reclassification to obtain spatially distributed LAI data (Figure 4b). Figure 4b shows the resulting high resolution vineyard LAI map, derived from the combination of the vegetation index and only five PocketLAI estimates (smart scouting).

PocketLAI and NDVI
The pixels of the georeferenced NDVI map of the vineyard were clustered into five classes using the ArcGIS v.10 package, and PocketLAI estimates were collected at the sixth survey in sampling area representative of each class (Figure 4a). A regression equation (R 2 = 0.90) was then derived between the LAIv obtained from PocketLAI and NDVI values in the sampling areas (indicated by black circles in Figure 4a). The regression equation was then used to convert all the other NDVI pixels of the map into LAIv values, using the spatial analysis tool of ArcGIS for raster reclassification to obtain spatially distributed LAI data (Figure 4b). Figure 4b shows the resulting high resolution vineyard LAI map, derived from the combination of the vegetation index and only five PocketLAI estimates (smart scouting).

Discussion
In general, the LAI values estimated using PocketLAI and those measured using destructive method were consistent, with a coefficient of determination ranging from 0.77 (Figure 2a) to 0.94 (Figure 2d), depending on the inclusion or not in the analysis of a vineyard LAI close to 3 (LAIv = 2.84 with destructive method), observed in a unusually high-vigor sampling area due to the absence of green pruning applications. This value, converted in vine LAI (Equation (2)), corresponds to an observed plant LAI close to 7 (LAIp = 6.8).
This value can be considered as outlying in vineyards, as confirmed, e.g., by the LAI ranges reported by Pergher and Petris [12]. For similar vine training systems, this study shows LAI values ranging from 0 to 4, measured across many seasons (from 1993 to 2005) from leaf development (code 13 BBCH scale) to the beginning of ripening (code 81 BBCH scale).
Excluding this value, the maximum vineyard LAI observed was below 2 (LAIv = 1.41 with direct method), corresponding to a vine LAI lower than 4 (LAIp = 3.4). These results confirmed the tendency of PocketLAI (and of most indirect methods based on the inversion of light transmittance models) to underestimate high LAI values because of a saturation effect [26]. In case of PocketLAI, the effect was already discussed by Francone et al. [31] for dense maize and giant reed canopies. In that case, the authors obtained pronounced underestimations especially in case of giant reed canopies exceeding 5. The reasons were the saturation effect and the strong deviation of industrial giant reed plants from the ideal assumption behind light transmittance models (randomly distributed infinitely small leaves). The same explanation is valid for the case of vineyard (leaves clumped in walls). However, PocketLAI can be considered a suitable tool for monitoring Vitis vinifera LAI development, since it showed inaccuracies only in an uncommon case of vine untreated with green pruning. This practice is usually adopted in viticulture, in order to decrease the shading on the grapes and create microclimatic condition inside the canopy unfavorable to pest infection [43].
Concerning the smart scouting activities, the coupling between NDVI map and in-situ LAI data, measured with PocketLAI, proved to be an interesting approach able to provide a description of the vineyard canopy LAI in absolute terms. The possibility to convert an NDVI map, able to describe

Discussion
In general, the LAI values estimated using PocketLAI and those measured using destructive method were consistent, with a coefficient of determination ranging from 0.77 (Figure 2a) to 0.94 (Figure 2d), depending on the inclusion or not in the analysis of a vineyard LAI close to 3 (LAIv = 2.84 with destructive method), observed in a unusually high-vigor sampling area due to the absence of green pruning applications. This value, converted in vine LAI (Equation (2)), corresponds to an observed plant LAI close to 7 (LAIp = 6.8).
This value can be considered as outlying in vineyards, as confirmed, e.g., by the LAI ranges reported by Pergher and Petris [12]. For similar vine training systems, this study shows LAI values ranging from 0 to 4, measured across many seasons (from 1993 to 2005) from leaf development (code 13 BBCH scale) to the beginning of ripening (code 81 BBCH scale).
Excluding this value, the maximum vineyard LAI observed was below 2 (LAIv = 1.41 with direct method), corresponding to a vine LAI lower than 4 (LAIp = 3.4). These results confirmed the tendency of PocketLAI (and of most indirect methods based on the inversion of light transmittance models) to underestimate high LAI values because of a saturation effect [26]. In case of PocketLAI, the effect was already discussed by Francone et al. [31] for dense maize and giant reed canopies. In that case, the authors obtained pronounced underestimations especially in case of giant reed canopies exceeding 5. The reasons were the saturation effect and the strong deviation of industrial giant reed plants from the ideal assumption behind light transmittance models (randomly distributed infinitely small leaves). The same explanation is valid for the case of vineyard (leaves clumped in walls). However, PocketLAI can be considered a suitable tool for monitoring Vitis vinifera LAI development, since it showed inaccuracies only in an uncommon case of vine untreated with green pruning. This practice is usually adopted in viticulture, in order to decrease the shading on the grapes and create microclimatic condition inside the canopy unfavorable to pest infection [43].
Concerning the smart scouting activities, the coupling between NDVI map and in-situ LAI data, measured with PocketLAI, proved to be an interesting approach able to provide a description of the vineyard canopy LAI in absolute terms. The possibility to convert an NDVI map, able to describe differences in space only in relative terms, to a LAI map acquires a special relevance when the optimization of farming practices requires a description of the vineyard based on absolute values, like in the case of pesticide dose adjustment.
In case of DHP, the correlation between indirect estimates and destructive measurements was fully satisfactory (R 2 = 0.94) and the method did not show saturation effect, although estimates were consistently affected by a relevant overestimation error. The overestimation was quite constant and independent from absolute LAI values, since it was observed throughout the season. Similar results were found by Lopez-Lozano et al. [25], who compared different tools for LAI estimation in vineyards, including DHP. The authors discussed that the 'true LAI' values provided by the same software we used (CAN-EYE) led to overestimating destructive LAI values by a factor around 1.6 (the value we obtained here was around 2.6). These results suggested the possibility to apply a specific correction factor for DHP improving the performance of the method for the vineyard case.
Neither PocketLAI nor DHP showed changes in reliability for different phenological stages. Excluding the outlying value discussed above, differences in mean absolute error never exceeded 0.17 and 0.44 for PocketLAI and DHP, respectively, and no trends were observed for the evaluation metrics as a function of development stage.

Conclusions
Both indirect methods based on digital image processing, i.e., the PocketLAI smart-app and hemispherical photography, proved their reliability for monitoring the spatial and temporal variability of vine vigor in cordon-trained system and to provide LAI values. However, DHP highlighted the need to assess a correction factor because of overestimation errors already discussed for vineyards [25], whereas PocketLAI showed a certain saturation effect for vine LAI approaching or exceeding 5 or 6. However, the latter is an uncommon condition with respect to the most widespread vine training systems, because the green pruning lead to lower LAI values.
Moreover, with respect to DHP, PocketLAI entails the advantages of automatic (completely unsupervised) and real-time provision of LAI data and the suitability for operational purposes, in addition to its low-cost and user-friendliness. In this context, PocketLAI shows a potential for a wide range of applications, such as the spatial-temporal optimization of the agronomic inputs, especially in combination with remote sensing products.
Author Contributions: F.O. performed the data analysis; set the equations to convert plant-LAI into vineyard-LAI; established the protocol for indirect measurements, comparing different tests. L.B. and R.C. conceived and designed the experiments. F.O. and R.C. wrote the paper. E.M. performed the hemispherical images processing, exploring the best options setting with CanEye. E.M. and D.C. carried out the indirect measurements. S.P. and T.G. carried out the destructive LAI measurements. S.P. set the figures. V.P. performed the elaboration concerning the remote sensing data.

Conflicts of Interest:
The authors declare no conflict of interest.