Comparison of UAV Photogrammetry and 3D Modeling Techniques with Other Currently Used Methods for Estimation of the Tree Row Volume of a Super-High-Density Olive Orchard

: A comparison of three di ﬀ erent methods to evaluate the tree row volume (TRV) of a super-high-density olive orchard is presented in this article. The purpose was to validate the suitability of unmanned aerial vehicle (UAV) photogrammetry and 3D modeling techniques with respect to manual and traditional methods of TRV detection. The use of UAV photogrammetry can reduce the amount of estimated biomass and, therefore, reduce the volume of pesticides to be used in the ﬁeld by means of more accurate prescription maps. The presented comparison of methodologies was performed on an adult super-high-density olive orchard, planted with a density of 1660 trees per hectare. The ﬁrst method (TRV 1 ) was based on close-range photogrammetry from UAVs, the second (TRV 2 ) was based on manual in situ measurements, and the third (TRV 3 ) was based on a formula from the literature. The comparisons of TRV 2 -TRV 1 and TRV 3 -TRV 1 showed an average value of the di ﬀ erence equal to + 13% (max: + 65%; min: − 11%) and + 24% (max: + 58%; min: + 5%), respectively. The results show that the TRV 1 method has high accuracy in predicting TRV with minor working time expenditure, and the only limitation is that professionally skilled personnel is required.


Introduction
Super-high-density (SHD) olive cropping systems represent a very interesting proposal for olive oil orchard profitability because they considerably reduce production costs, thanks to full mechanization from planting to harvesting, and provide higher, more constant crop levels. SHD olive trees grow in hedgerows that form continuous fruit-bearing canopies, which are suitable for over-the-row harvesting machines. Many agronomic aspects are currently under study for SHD olive culture, such as soil management, harvesting times, cropping density, irrigation management, tree architecture, and precision farming [1][2][3][4]. Precision agriculture, especially for fruit tree crops, is an environmental and economic management strategy that uses information and communication technology to acquire data that leads to decisions aimed at agricultural production [5] in the context of climate change. The use of unmanned aerial vehicles (UAVs) for exact cultivation detection from 3D point cloud maps [6] and digital surface models (DSMs) [7], which are generated from multispectral imagery, has accelerated the extrapolation of valuable information for the optimal management of sustainable agriculture. This new frontier is exemplified by the use of data fusion for delineating homogeneous management zones in precision agriculture [8]. Furthermore, remote sensing can be used for precision agriculture to discern the "health status" of crops [9] in a measurement system net located in scattered fields. Advances in precision agriculture have allowed for prompt, site-specific management of the limiting factors of the production system, by which decisions to apply resources and agronomic practices are dictated in progress by the variability detected in the field among the crops [10]. Site-specific management means that inputs, namely, irrigation water, fertilizers, and pesticides, are applied only where and when necessary to maximize the desired result, which in most cases is linked to income maximization and reduced environmental impact [11]. For instance, pest management in agriculture substantially increases crop productivity, but at the same time, it could increase risks to human health and the environment. The Directive 2009/128/EC compels farmers to achieve the sustainable use of pesticides by reducing these impacts. All this has led to the spread of both integrated pest management and organic agriculture. Moreover, recently, precision agriculture has been used in other agronomic practices, such as estimating height [12], tree biomass [13][14][15][16][17], fruit tree crop, pollination, and pruning and harvesting amount, as well as for disease detection [18]. The use of these very new technologies has aroused greater environmental awareness and attention to human health for both consumers and fruit growers, along with awareness of the economic costs of treatment for farmers. The smart and integrated management of high-resolution environmental data from sensors and images represents a low-cost method for setting prescription maps and has been increasingly studied for agricultural applications [19].
The reduction of estimation errors in tree row volume (TRV) calculations could be a very interesting application of these new technologies in order to obtain the correct volume of pesticides for use in treatments with sprayer machines; this could allow for reductions in costs and environmental pesticide dispersions by as much as 50%. There are different techniques for obtaining tree volumes from a point cloud, such as an adjusted tree volume generated from a light detection and ranging (LiDAR) sensor [20,21] and UAV techniques based on airborne laser scanning and structure from motion (SfM) point clouds [22]. The high costs of very high-resolution imagery provided by satellites or LiDAR sensors often impede estimations of aboveground biomass and the determination of other vegetation parameters [23]. A comparison of three different methods to evaluate the TRV of a super-high-density olive orchard is presented in this article. The purpose was to validate the suitability of UAV photogrammetry and 3D modeling techniques with respect to manual and traditional methods of TRV detection. The first method (TRV 1 ) was used to estimate the volume of an SHD olive grove using UAV systems coupled with numerical software. TRV 2 estimated the volume of the SHD olive tree canopy measured manually in situ by means of measuring instruments, such as an optical meter and a portable GPS locator; TRV 2 was the reference method. Finally, TRV 3 was used to estimate the volume of the olive grove by the traditional method [24]. Subsequently, TRV 1 and TRV 3 were compared with TRV 2 . The workflow of the TRV tree detection process with the three different methods is shown in Figure 1.

Materials and Methods
The research was carried out at the experimental center "P. Martucci" of the University of Bari in Valenzano (Bari), Italy (41.025800 N, 16.907563 W, 140 m a.s.l.), where the TRV of an adult (13 years old) super-high-density olive orchard was estimated ( Figure 2).
The olive orchard had a density of 1660 trees per hectare. The field examined consisted of 24 rows including 15 olive cultivars: Arbequina (3 rows); Arbosana (3 rows); Carolea, Cima di Bitonto, and Coratina (2 rows); Don Carlo and Frantoio (2 rows); Fs-17 and I/77, Koroneiki (3 rows); and Leccino, Maurino, Nociara, Peranzana, and Urano (2 rows) ( Figure 3). The study was conducted in winter 2019 during olive tree rest. For each of the three adopted methods of TRV estimation, the entire volume of each of the 24 rows was taken into account. No pruning was performed.

Materials and Methods
The research was carried out at the experimental center "P. Martucci" of the University of Bari in Valenzano (Bari), Italy (41.025800 N, 16

Materials and Methods
The research was carried out at the experimental center "P. Martucci" of the University of Bari in Valenzano (Bari), Italy (41.025800 N, 16.907563 W, 140 m a.s.l.), where the TRV of an adult (13 years old) super-high-density olive orchard was estimated ( Figure 2).
The olive orchard had a density of 1660 trees per hectare. The field examined consisted of 24 rows including 15 olive cultivars: Arbequina (3 rows); Arbosana (3 rows); Carolea, Cima di Bitonto, and Coratina (2 rows); Don Carlo and Frantoio (2 rows); Fs-17 and I/77, Koroneiki (3 rows); and Leccino, Maurino, Nociara, Peranzana, and Urano (2 rows) ( Figure 3). The study was conducted in winter 2019 during olive tree rest. For each of the three adopted methods of TRV estimation, the entire volume of each of the 24 rows was taken into account. No pruning was performed.  The olive orchard had a density of 1660 trees per hectare. The field examined consisted of 24 rows including 15 olive cultivars: Arbequina (3 rows); Arbosana (3 rows); Carolea, Cima di Bitonto, and Coratina (2 rows); Don Carlo and Frantoio (2 rows); Fs-17 and I/77, Koroneiki (3 rows); and Leccino, Maurino, Nociara, Peranzana, and Urano (2 rows) ( Figure 3). The study was conducted in winter 2019 during olive tree rest. For each of the three adopted methods of TRV estimation, the entire volume of each of the 24 rows was taken into account. No pruning was performed.  This study compared three different methods for calculating TRV in order to assess their reliability for reproducing tree crown geometry and volume.

TRV1 Estimation
The first method used for the estimation of the TRV of the SHD olive orchard was the UAV closerange photogrammetry and 3D modeling technique. Figure 4 shows the 21 ground control points (GCPs) that were positioned on the ground in order to correlate the GPS positions of the images taken by the UAV with the GPS position measured in the field of the GCPs.  This study compared three different methods for calculating TRV in order to assess their reliability for reproducing tree crown geometry and volume.

TRV 1 Estimation
The first method used for the estimation of the TRV of the SHD olive orchard was the UAV close-range photogrammetry and 3D modeling technique. Figure 4 shows the 21 ground control points (GCPs) that were positioned on the ground in order to correlate the GPS positions of the images taken by the UAV with the GPS position measured in the field of the GCPs.
The UAV used for the experiment was a DJI SPARK, which is a low-cost UAV that was selected in order to obtain a 3D model of the field without using very expensive equipment that, at the same time, maintains good precision. The study area had an extension of 186 × 36.11 m 2 . Standard RGB images were acquired with a regular grid pattern, frontal overlap of 75%, side overlap of 60%, flight height of 20 m, image size of 3968 × 2976 pixels, a field of view (FoV) of 81.9 • , and ground sampling distance (GSD) of 0.68 cm/pixel. The flight plan was designed to minimize flight time. The UAV was equipped with a CMOS 1/2.3" sensor of 12 MP, an exposure time of 1/240 s, and an ISO of 100.
The photogrammetric 3D map, the dense point cloud, and the mesh of all the olive crowns were processed by Pix4Dmapper software (Version 4.3.33; Pix4D, Lausanne, Switzerland) [25].

TRV1 Estimation
The first method used for the estimation of the TRV of the SHD olive orchard was the UAV closerange photogrammetry and 3D modeling technique. Figure 4 shows the 21 ground control points (GCPs) that were positioned on the ground in order to correlate the GPS positions of the images taken by the UAV with the GPS position measured in the field of the GCPs.  After reconstructing the 3D model, the total TRV 1 of each row was calculated by the software using the volume function and selecting the upper half-space parallel to the surface of the ground and positioned below the crown of the trees. The GCP measurements were used for recovering the camera interior parameters and adjusting for any systematic error or block distortion. Figure 5 provides an analysis of the procedure used for one particular day to quantify the TRV. The point cloud was converted into a mesh by applying the Pix4Dmapper software. Finally, subtractions from these meshes gave us the tree crown volume of each row. The UAV used for the experiment was a DJI SPARK, which is a low-cost UAV that was selected in order to obtain a 3D model of the field without using very expensive equipment that, at the same time, maintains good precision. The study area had an extension of 186 × 36.11 m 2 . Standard RGB images were acquired with a regular grid pattern, frontal overlap of 75%, side overlap of 60%, flight height of 20 m, image size of 3968 × 2976 pixels, a field of view (FoV) of 81.9°, and ground sampling distance (GSD) of 0.68 cm/pixel. The flight plan was designed to minimize flight time. The UAV was equipped with a CMOS 1/2.3″ sensor of 12 MP, an exposure time of 1/240 s, and an ISO of 100.
The photogrammetric 3D map, the dense point cloud, and the mesh of all the olive crowns were processed by Pix4Dmapper software (Version 4.3.33; Pix4D, Lausanne, Switzerland) [25].
After reconstructing the 3D model, the total TRV1 of each row was calculated by the software using the volume function and selecting the upper half-space parallel to the surface of the ground and positioned below the crown of the trees. The GCP measurements were used for recovering the camera interior parameters and adjusting for any systematic error or block distortion. Figure 5 provides an analysis of the procedure used for one particular day to quantify the TRV. The point cloud was converted into a mesh by applying the Pix4Dmapper software. Finally, subtractions from these meshes gave us the tree crown volume of each row.

TRV2 Estimation
To estimate the TRV2 of the SHD olive grove, an in situ measurement, aimed at evaluating the TRV* of a single tree, was carried out by means of simple measuring instruments, such as an optical meter (Hilti PD-E, measuring range of 0-200 m, accuracy of ±0.001 m) and a portable GPS locator (GARMIN GPSMAP 64 S).
Once the TRV*s of 10 trees per row were calculated, the TRV2 of each row was evaluated by multiplying the number of trees in a single row by the mean value of the obtained TRV*s. This technique can be divided into two phases: (1) acquisition of the tree biometric parameters in two sampling days, and (2) data analysis and tree volume processing by means of MS Excel software.
To evaluate each tree's TRV*, the tree height was measured with a distance meter placed vertically on the ground. Measurements of the height of the vegetation and the height of the subcrown were then acquired, taking into proper consideration the real shape of each tree. Horizontally, the width of the tree crown was measured considering the end of the most protrusive branch.
To calculate the TRV*, the tree row geometry, considering its vertical section, was approximated after pruning to an isosceles trapezoid ( Figure 6); in this way, the TRV* was determined by multiplying the trapezoid area by the crown dimension in its orthogonal directions (Formula (1)):

TRV 2 Estimation
To estimate the TRV 2 of the SHD olive grove, an in situ measurement, aimed at evaluating the TRV* of a single tree, was carried out by means of simple measuring instruments, such as an optical meter (Hilti PD-E, measuring range of 0-200 m, accuracy of ±0.001 m) and a portable GPS locator (GARMIN GPSMAP 64 S).
Once the TRV*s of 10 trees per row were calculated, the TRV 2 of each row was evaluated by multiplying the number of trees in a single row by the mean value of the obtained TRV*s. This technique can be divided into two phases: (1) acquisition of the tree biometric parameters in two sampling days, and (2) data analysis and tree volume processing by means of MS Excel software.
To evaluate each tree's TRV*, the tree height was measured with a distance meter placed vertically on the ground. Measurements of the height of the vegetation and the height of the subcrown were then acquired, taking into proper consideration the real shape of each tree. Horizontally, the width of the tree crown was measured considering the end of the most protrusive branch.
To calculate the TRV*, the tree row geometry, considering its vertical section, was approximated after pruning to an isosceles trapezoid ( Figure 6); in this way, the TRV* was determined by multiplying the trapezoid area by the crown dimension in its orthogonal directions (Formula (1)): The trees crowns were highly linked to each other as a hedgerow; for this reason, P was assumed by the ratio of the entire length of the row per the number of trees of each row. Furthermore, b was assumed to be equal to 80% of B, since the base of a tree has lower vegetative growth than the apical zone and is similar to the geometry of an inverted truncated wedge. The value of 80% was obtained by measuring the average inclination of all sampling trees.
The advantages of using TRV2 include the possibility of sampling any type of tree, but above all, the simplicity of the necessary instrumentation. The drawback is the very long sampling time, which limits the number of samples that can be taken for each cultivar, and therefore, the statistical representativeness of the samples.

TRV3 Estimation
Finally, the third procedure, usually adopted for estimating TRV3 [24], consisted of multiplying the average values of the height and thickness of the crown for each row (Figure 7), multiplied by the area of 1 ha and divided by the inter-row width (Formula (2)). The crown was considered to be a rectangular shape. In our case study, as the surface was less than 1 ha, the general formula was The trees crowns were highly linked to each other as a hedgerow; for this reason, P was assumed by the ratio of the entire length of the row per the number of trees of each row. Furthermore, b was assumed to be equal to 80% of B, since the base of a tree has lower vegetative growth than the apical zone and is similar to the geometry of an inverted truncated wedge. The value of 80% was obtained by measuring the average inclination of all sampling trees.
The advantages of using TRV 2 include the possibility of sampling any type of tree, but above all, the simplicity of the necessary instrumentation. The drawback is the very long sampling time, which limits the number of samples that can be taken for each cultivar, and therefore, the statistical representativeness of the samples.

TRV 3 Estimation
Finally, the third procedure, usually adopted for estimating TRV 3 [24], consisted of multiplying the average values of the height and thickness of the crown for each row (Figure 7), multiplied by the area of 1 ha and divided by the inter-row width (Formula (2)). The crown was considered to be a rectangular shape. In our case study, as the surface was less than 1 ha, the general formula was revised by multiplying by the effective row ground area divided by 10,000 m 2 : where H is the average crown height in a row (m), C is the average crown thickness in a row (m), R is the inter-row width (m), and A is the effective row ground area (m 2 ).
Agriculture 2019, 9, x FOR PEER REVIEW 7 of 14 Figure 7. Tree crown geometry for the TRV3 estimation method.
where H is the average crown height in a row (m), C is the average crown thickness in a row (m), R is the inter-row width (m), and A is the effective row ground area (m 2 ). An overall graphical representation of the three used methods is reported in Figure 8. The TRV1 method allows for better reconstruction of the tree shape, unlike the TRV2 and TRV3 methods ( Figure 8). However, the TRV2 method is based on field measurements (10 samples per row) and should be the most representative of reality. TRV3, that is, the traditional method, is the one that most differs from the true TRV value, as it is based on a rough estimate. It should also be pointed out that the TRV3 calculated in this paper is very precise, considering that all the dimensions used in Formula (2) were derived from the averages of the biometric parameters of 10 trees per row; in traditional practice, it is taken as a reference for the calculation of the TRV3 of a single pair of trees per hectare, which produces an estimated value that is very far from the true value.  An overall graphical representation of the three used methods is reported in Figure 8. The TRV 1 method allows for better reconstruction of the tree shape, unlike the TRV 2 and TRV 3 methods (Figure 8). However, the TRV 2 method is based on field measurements (10 samples per row) and should be the most representative of reality. TRV 3 , that is, the traditional method, is the one that most differs from the true TRV value, as it is based on a rough estimate. It should also be pointed out that the TRV 3 calculated in this paper is very precise, considering that all the dimensions used in Formula (2) were derived from the averages of the biometric parameters of 10 trees per row; in traditional practice, it is taken as a reference for the calculation of the TRV 3 of a single pair of trees per hectare, which produces an estimated value that is very far from the true value.

R
where H is the average crown height in a row (m), C is the average crown thickness in a row (m), R is the inter-row width (m), and A is the effective row ground area (m 2 ). An overall graphical representation of the three used methods is reported in Figure 8. The TRV1 method allows for better reconstruction of the tree shape, unlike the TRV2 and TRV3 methods ( Figure 8). However, the TRV2 method is based on field measurements (10 samples per row) and should be the most representative of reality. TRV3, that is, the traditional method, is the one that most differs from the true TRV value, as it is based on a rough estimate. It should also be pointed out that the TRV3 calculated in this paper is very precise, considering that all the dimensions used in Formula (2) were derived from the averages of the biometric parameters of 10 trees per row; in traditional practice, it is taken as a reference for the calculation of the TRV3 of a single pair of trees per hectare, which produces an estimated value that is very far from the true value.

Results
All the results obtained from each of the described methods are shown in Table 1. In particular, for the TRV 1 method, the reported results were derived from the software evaluation; for the other two methods, all the parameters were used in their respective evaluation formulas (Formulas (1) and (2)). The first statistical analysis of the three methods (TRV 1 , TRV 2 , and TRV 3 ) on the 24 rows showed that, although the mean values of the three methods were very similar and the asymmetry was very low (close to zero), the standard deviation was very high and the kurtosis was very low compared with the threshold value of 3 of the Gaussian distributions for all of them (Table 2, Figure 9). This implies that the distribution of volumes cannot be considered Gaussian and the TRVs changed considerably from cultivar to cultivar, regardless of the method used.   Figures 10 and 11 show the scatter plots related to the data couples TRV1-TRV2 and TRV1-TRV3, respectively. Each figure also displays the regression line evaluated in order to assess the correlation between the datasets used, with a confidence interval (CI) of 95%. The correlation between TRV1-TRV2 showed a Pearson's coefficient r1-2 = 0.831 (p-value < 0.0001), while the correlation between TRV1-TRV3 showed a Pearson's coefficient r1-3 = 0.874 (p-value < 0.0001). Both obtained results indicate that evaluating tree volume by UAVs has a high level of correlation with manual and traditional methods, regardless of the analyzed cultivar.  Figures 10 and 11 show the scatter plots related to the data couples TRV 1 -TRV 2 and TRV 1 -TRV 3 , respectively. Each figure also displays the regression line evaluated in order to assess the correlation between the datasets used, with a confidence interval (CI) of 95%. The correlation between TRV 1 -TRV 2 showed a Pearson's coefficient r 1-2 = 0.831 (p-value < 0.0001), while the correlation between TRV 1 -TRV 3 showed a Pearson's coefficient r 1-3 = 0.874 (p-value < 0.0001). Both obtained results indicate that evaluating tree volume by UAVs has a high level of correlation with manual and traditional methods, regardless of the analyzed cultivar.  Figures 10 and 11 show the scatter plots related to the data couples TRV1-TRV2 and TRV1-TRV3, respectively. Each figure also displays the regression line evaluated in order to assess the correlation between the datasets used, with a confidence interval (CI) of 95%. The correlation between TRV1-TRV2 showed a Pearson's coefficient r1-2 = 0.831 (p-value < 0.0001), while the correlation between TRV1-TRV3 showed a Pearson's coefficient r1-3 = 0.874 (p-value < 0.0001). Both obtained results indicate that evaluating tree volume by UAVs has a high level of correlation with manual and traditional methods, regardless of the analyzed cultivar. Although the regression lines show a high correlation index between the different methods used, it is of significance to analyze the differences between the absolute values measured by the three TRVs per each row considering TRV 1 as the reference for the other two methods. In particular, the relative difference between TRV 2 and TRV 1 and between TRV 3 and TRV 1 was calculated. Although the regression lines show a high correlation index between the different methods used, it is of significance to analyze the differences between the absolute values measured by the three TRVs per each row considering TRV1 as the reference for the other two methods. In particular, the relative difference between TRV2 and TRV1 and between TRV3 and TRV1 was calculated.
The average value per row of the relative difference between TRV2 and TRV1 was equal to 13%, while the minimum value was equal to −11% for the Arbequina I cultivar and the maximum value was +65% for the replica of the Arbosana II cultivar (Figure 12). The outlier of 65% of the relative difference between TRV2 and TRV1 was mainly generated by the low density of the replica of the Arbosana II cultivar and the consequent underestimation effect of the volume generated by the 3D photogrammetry of the tree crowns. The average value per row of the relative difference between TRV3 and TRV1 was equal to +24%, while the minimum value was equal to 5% for the Arbequina I and Maurino cultivars and the maximum value was +58% for the Frantoio I cultivar.  (TRV3-TRV1) /TRV1 Figure 11. Linear regression graph between TRV 3 and TRV 1 .
The average value per row of the relative difference between TRV 2 and TRV 1 was equal to 13%, while the minimum value was equal to −11% for the Arbequina I cultivar and the maximum value was +65% for the replica of the Arbosana II cultivar (Figure 12). The outlier of 65% of the relative difference between TRV 2 and TRV 1 was mainly generated by the low density of the replica of the Arbosana II cultivar and the consequent underestimation effect of the volume generated by the 3D photogrammetry of the tree crowns. The average value per row of the relative difference between TRV 3 and TRV 1 was equal to +24%, while the minimum value was equal to 5% for the Arbequina I and Maurino cultivars and the maximum value was +58% for the Frantoio I cultivar.
Although the regression lines show a high correlation index between the different methods used, it is of significance to analyze the differences between the absolute values measured by the three TRVs per each row considering TRV1 as the reference for the other two methods. In particular, the relative difference between TRV2 and TRV1 and between TRV3 and TRV1 was calculated.
The average value per row of the relative difference between TRV2 and TRV1 was equal to 13%, while the minimum value was equal to −11% for the Arbequina I cultivar and the maximum value was +65% for the replica of the Arbosana II cultivar (Figure 12). The outlier of 65% of the relative difference between TRV2 and TRV1 was mainly generated by the low density of the replica of the Arbosana II cultivar and the consequent underestimation effect of the volume generated by the 3D photogrammetry of the tree crowns. The average value per row of the relative difference between TRV3 and TRV1 was equal to +24%, while the minimum value was equal to 5% for the Arbequina I and Maurino cultivars and the maximum value was +58% for the Frantoio I cultivar.   From the difference of the mean values between the various cultivars of (TRV 2 − TRV 1 )/(TRV 1 ) and (TRV 3 − TRV 1 )/(TRV 1 ), it can be noted that the TRV 1 method obtained average results lower than the TRV 2 and TRV 3 methods. This was due to the 3D modeling of the canopy that thickened the actual shape of the tree crown compared with the TRV 2 and TRV 3 methods.
Considering the total values of the volumes calculated with the three methods reported in Table 1, it is worth highlighting that TRV 2 had a difference of 10% more than the total volume calculated with TRV 1 , and TRV 3 had a difference of +23% more than the total volume calculated with TRV 1 . Also, in absolute terms, the TRV 1 method succeeded in providing a better estimate of the canopy biomass, unlike TRV 3 , which always overestimated the volume of the tree crown.
The geometrical measurements shown in Table 1 for the TRV 2 method, as mentioned before, correspond to the average of the measurements of 10 trees per row. This method is very laborious and requires a large number of specialized operators working in the field for hours. With TRV 2 , the surveys were carried out over two days, the data was backed up and processed in Excel in four days, and there were four operators employed. In comparison, the TRV 1 method required about 15 min of flight per half hectare, 4-5 h of calculation using the software, and a single specialized operator.

Discussion
The respective Pearson coefficients that were calculated demonstrated that the correlation between the methods was high and the analysis was statistically significant.
In general, the 3D modeling technique (TRV 1 ) underestimated crown volume with respect to the reference method (TRV 2 ), which did not take into account the density of the vegetation and considered empty spaces as full. This effect was indeed related to canopy density and the architecture of the different cultivars under study.
Where the canopies were normally dense, which was the case for most of the cultivars, TRV 1 "saw" all the effective density of vegetation. Indeed, the TRV 1 technique underestimated the crown volume for the rows by about 40% with respect to the other two methods, as Cima di Bitonto, Koroneiki II, Frantoio I, Carolea, and Arbosana II were characterized by low canopy density and defoliated rows [2,26]. On the contrary, no significant differences between the first method and the other two were found for Arbequina I, which was characterized by a canopy density as dense as Arbequina III and Maurino.
The comparisons between TRV 2 -TRV 1 and TRV 3 -TRV 1 showed an average value of the difference of +13% (max: +65%, min: −11%; depending on the cultivar) and +24% (max: +58%, min: +5%; depending on the cultivar), respectively. The obtained results demonstrate that the TRV 1 method is highly accurate at predicting the TRV and requires low labor effort and time expenditure for both olive trees and similar woody crops.
This result demonstrates that it is possible to measure the canopy volume in a super-high-density olive orchard using UAVs to improve the management of spatial variability and avoid tedious and time-consuming fieldwork in comparison with the other methods.
In this experiment, TRV 3 was very accurate compared with the values usually determined by farmers because the data of the geometric dimensions of the trees were well known, but normally, the value of TRV 3 is much higher and more inaccurate. The TRV 3 method overestimated crown volume, especially for the cultivars Cima di Bitonto, Frantoio, and Carolea, which are characterized by very high vigor and a trapezoidal shape.
However, all the methods used allowed us to obtain the same statistical results, and this reinforces the correctness of the experimental analysis carried out in this work. This is one of the first studies aimed at evaluating the tree row volume of a large number of cultivated genotypes for an SHD olive orchard using UAV photogrammetry and 3D modeling techniques and comparing them with other TRV evaluation methods.
The use of UAVs improves the management of spatial variability and avoids intensive fieldwork. The traditional method (TRV 3 ) has very rough TRV estimation errors; in the present work, it was +23%. The error was related to some genetic characteristics of the cultivars, such as tree vigor, crown architecture, and density. An estimate so coarse leads to cascading errors for the estimation of the leaf mass, the quantity of pesticides to be used in the field, and the estimate of the quantity of water for irrigation. The TRV 1 method can obtain greater definition with respect to traditional surveys and provides the possibility of reaching otherwise inaccessible areas in a time-and cost-saving manner. Furthermore, it guarantees the analysis of very large areas in a very short time and high survey precision. Finally, the main limitation of the first method is that professionally skilled personnel is required; in fact, the analysis and 3D modeling of the foliage assume that the operator has specific knowledge of flight and photogrammetry techniques and knows how to use complicated hardware and software equipment. It must also be considered however, that today technology has playing a major role in managing different needs of the agricultural industry, as the reduction of energy [27,28], agricultural biomass [29] and working time [30] and then increasingly specialized agricultural operators are required.