High-Resolution Airborne UAV Imagery to Assess Olive Tree Crown Parameters Using 3 D Photo Reconstruction : Application in Breeding Trials

The development of reliable methods for the estimation of crown architecture parameters is a key issue for the quantitative evaluation of tree crop adaptation to environment conditions and/or growing system. In the present work, we developed and tested the performance of a method based on low-cost unmanned aerial vehicle (UAV) imagery for the estimation of olive crown parameters (tree height and crown diameter) in the framework of olive tree breeding programs, both on discontinuous and continuous canopy cropping systems. The workflow involved the image acquisition with consumer-grade cameras on board a UAV and orthomosaic and digital surface model generation using structure-from-motion image reconstruction (without ground point information). Finally, geographical information system analyses and object-based classification were used for the calculation of tree parameters. Results showed a high agreement between remote sensing estimation and field measurements of crown parameters. This was observed both at the individual tree/hedgerow level (relative RMSE from 6% to 20%, depending on the particular case) and also when average values for different genotypes were considered for OPEN ACCESS Remote Sens. 2015, 7 4214 phenotyping purposes (relative RMSE from 3% to 16%), pointing out the interest and applicability of these data and techniques in the selection scheme of breeding programs.


Introduction
Plant architecture is a key issue for the understanding of many aspects of plant ecology and is also one of the most important traits for the characterization and monitoring of fruit trees [1].This is especially important in fruit breeding programs aimed at developing new cultivars adapted to specific growing systems [2].In the case of olive tree, breeding could be conducted to produce cultivars suitable either to discontinuous canopy (open vase configuration) or continuous canopy (hedgerow) [3][4][5].For such purposes, reliable, cost-efficient and robust methods able to characterize the tree crown of the genotypes under selection are desirable.Few attempts have been reported to characterize tree architecture in olive breeding [6][7][8] aiming at a precise characterization of the tree crown.An accurate 3D model representing the shape and size of the tree crowns can provide information about critical aspects of plant growth and development [9] and, therefore, about its suitability for the mentioned training systems.
In the last few decades, a number of alternatives to the traditional manual-based measurement of morphological parameters for the characterization of 3D structure of trees and canopies has emerged.Most of these alternatives are based on the use of remote sensors, both indirectly by modelling 3D features from the spectral response recorded by passive sensors in the visible and infrared spectral regions or directly by active or passive sensors operating in different wavelengths capable of recording 3D information.Even if the former are particularly suited for monitoring crop condition, canopy density, etc., due to their relatively easy acquisition and analysis and their sensitivity to canopy cover and vegetation density, they tend to suffer from several shortcomings when used for the modelling of structural parameters, particularly regarding their reliability and accuracy [10,11].Among the latter, ultrasonic sensors, stereo images, radar and LIDAR or even X-ray [12,13] have been proposed, with LIDAR (light detection and ranging) laser pulse technology (both on terrestrial or airborne platforms) one of the most sound and accepted approaches for detailed measurement and monitoring of tree architecture [14,15], along with stereo vision techniques [13].
Despite the widespread use of LIDAR for the retrieval of structure parameters, it suffers from a series of shortcomings.Most of these shortcomings are related to its relative high costs of acquisition and processing, but also to the difficulties in the retrieval of the high density returns required for the accurate crown reconstruction in screening and/or phenotyping, particularly on small plots with high density tree rows.Traditional photogrammetric analyses also require high-cost cameras and platforms and rigorous processing chains for accurate results [11].Yet, the emerging techniques of computer vision enable the generation of reliable 2D and 3D imagery from collections of multi-angle imagery acquired from standard (non-metric) cameras by using 3D photo-reconstruction algorithms based on structure-from-motion (SfM) and multi-view-stereo analysis algorithms [16,17].One of the main advantages of these approaches lies in their potential for performing accurate 3D imagery reconstructions from datasets with a large number of images and a high overlap percentage, even if accurate internal and external camera orientation parameters are not available [18].
In fact, the new methods based on photo reconstruction enable the generation of very high resolution (VHR) digital surface models (DSM) and orthophotograph mosaics with high spatial and temporal resolution [18][19][20][21], and they have proven their applicability for environmental monitoring [22] and tree crown height quantification [11], being an attractive alternative to LIDAR technology, due to its performance, flexibility and relative low cost.
The simplicity, low weight and low cost of the sensors required for SfM image reconstruction, along with the development of miniaturized and low-cost inertial sensors, GPS devices and embedded computers, has enhanced the possibilities of the integration of these devices on light-weight unmanned aerial vehicles (UAVs) [23], currently referred to as unmanned aerial systems (UAS) and remotely-piloted aircraft systems (RPAS).In fact, the civil applications of UAV, such as high-resolution image acquisition, have emerged as an attractive option for agriculture and environment monitoring (see, e.g., [24][25][26]).Hence, they have been widely used in recent studies involving quantitative remote sensing applications due to their versatility.For example, [23,27,28] monitored several aspects of vegetation by using miniature thermal cameras, while [29][30][31] used narrow-band multispectral imagery for crop monitoring.Even more complex micro-hyperspectral sensors were successfully installed on board UAV platforms for crop monitoring purposes by [32][33][34].
Taking into account such progress, the present work aims at the development, assessment of the performance and potential uses of low-cost image acquisition systems on board UAV platforms for the estimation of olive crown parameters.More precisely, we integrated SfM image reconstruction techniques to derive digital surface models and ortho-mosaics and GIS and geographical object-based image analysis (GEOBIA) for the automatic estimation of crown diameter and tree height as two of the main variables of interest for tree crown architecture in plant phenotyping studies.We tested the performance of the above-mentioned methods against ground reference data obtained on two olive breeding trials designed to form a discontinuous canopy (open vase configuration) or continuous canopy (hedgerow).Special attention was paid to the assessment of the potential of the method for characterizing small differences in crown structure among the genotypes of those trials.

Study Area and Experimental Design
Two different field trials of the Córdoba breeding program were included in the present study (Figure 1).
The first trial (Study Area 1) was located at the Institute of Agricultural and Fishery Research and Training (IFAPA) Alameda del Obispo experimental farm in Córdoba, Spain (37°51′15′′N, 4°48′16′′W), and was planted in 2009.It consists of an intermediate breeding trial [35] arranged in a randomize design with four blocks and one tree per elementary plot.Trees were planted at 4 m × 2 m and trained as single trunk open vase at 1-m height to form a discontinuous canopy of scattered trees.Standard cultural practices, including drip irrigation, were followed in the orchard to ensure adequate tree growth.For the present study, 29 genotypes (4 trees per genotype) were randomly selected, and additional guard trees were also included for initial estimation of the parameters at the tree level for a total number of 150 individual trees.The second trial (Study Area 2) was located in commercial farm property of Todolivo S.L. cooperative in Pedro Abad, Córdoba, Spain (37°58′00′′N, 4°27′56′′W), and was planted in 2007.Three cultivars "Arbequina", "Arbosana" and "Sikitita" [36], and three selections of the cooperative breeding program carried out by the University of Cordoba (UC) and IFAPA (I) (UC-I6-9, UC-I7-8 and UC-I2-35) were arranged in a randomized design with 4 blocks.The elementary plot consists of 30 trees planted at 3.75 × 1.35-m spacing (1975 plants/ha) designed to form a continuous canopy (narrow hedgerow) to be harvested by over-the-row straddle harvesting machines.Trees were trained in a central leader up to 2.2 m.Irrigation and fertilization were performed to maximize accumulated yield.Additional guard hedgerows were also included for initial estimation of the parameters for a total number of 28 individual hedgerows.

Analysis Overview
The main steps of the research are summarized in Figure 2. The study was based on the analyses of collections of images acquired by consumer-grade cameras on board UAV platform for two study sites using the SfM image reconstruction approach.As detailed in the following sections, images were processed to obtain digital surface model and color infrared (CIR) orthomosaics.Data were analyzed to estimate tree heights and average crown diameters for the discontinuous canopy in Study Area 1 and hedgerow heights and widths in Study Area 2, respectively.These results were compared against field measures of these variables in order to assess the accuracies of the method, both at the individual tree and hedgerow level for the two study areas.Finally, ANOVA analyses were performed on the tree and hedgerow variables estimated and measured in order to assess the validity of the remote sensing method for comparison among genotypes.

Figure 2.
Flowchart illustrating the main analysis steps and outputs.Elements shaded in grey indicate field data and the processes used as a reference for validation.These datasets and processes were replicated in two olive orchards, one with discontinuous canopy and the other with continuous hedgerows.DifMin, difference from the minimum value.

Field Measurements
In Study Area 1, canopy height and width (measured east-west and north-south) were recorded for each tree and canopy surface, allowing an eventual estimation of crown volume from these measurements [37].

Object oriented classification
Tree/hedge height Tree/hedge crown outline

Genotype averaging
Tree/hedge diameter Tree/hedge height

Genotype average height
Genotype aver.crown diam.

Genotype average height
Genotype aver.crown diam.

CROWN ARQUITECTURE TREE/VARIETY ACCURACY COMPARISON AMONG GENOTYPES
In Study Area 2, canopy height and maximum width were measured in three trees per elementary plot (hedgerow).
All of the measurements were made in the same week as the aerial images were taken.

Image Acquisition and Pre-Processing
The airborne campaigns were conducted using an UAV operated by the Laboratory for Research Methods in Quantitative Remote Sensing (Quanta-Lab, IAS-CSIC, Spain).The UAV carried an RGB camera modified for color infrared (CIR) detection by removing the internal infrared filter.The camera was a Panasonic Lumix DMC-GF1 with a 4000 × 3000 pixel detector that captured images at f/3.2 and 1/2500 s with an angular FOV of 47.6° × 36.3° and provided ∼7-cm/pixel resolution at an approximate altitude of 300 m above ground level (AGL).The UAV platform was a 2-m wingspan fixed-wing platform with up to a 1-h endurance at 5.8 kg take-off weight (TOW) and 63-km/h ground speed (mX-SIGHT, UAV Services and Systems, Germany) that took 45 min to complete the entire flight plan at the indicated flight altitude.It was operated in the field by two trained operators and controlled by an autopilot for fully-autonomous flying (AP04, UAV Navigation, Madrid, Spain).The autopilot had a dual CPU controlling an integrated attitude heading reference system (AHRS) based on an L1GPS board, 3-axis accelerometers, gyros and a 3-axis magnetometer [23].The ground control station and the UAV were radio linked, transmitting position, attitude and status data at a 20-Hz frequency.The imagery was synchronized using the GPS position and the triggering time recorded for each image.The study areas were flown following a grid of parallel and perpendicular flight lines to ensure that each ground object was imaged in the along-and across-track directions of the UAV platform, maximizing overlap.A cloud of images comprising the flight lines covered the entire area homogeneously at a rate of 2 s between consecutive images.At an average of a 63-km/h ground speed, FOV = 47.6° × 36.3° and an approximate altitude of 300 m, the resulting average overlap obtained during the flights was extremely large and was estimated at 80%-90% in most cases.
In the process of SfM image reconstruction, 387 out of the original set of 390 (Study Area 1) and the whole original set of 255 images (Study Area 2) were used to generate each orthomosaic and the digital surface model for each orchard using pix4UAV software (®Ecublens, Switzerland).The processing steps included automatic aerial triangulation (AAT), bundle block adjustment (BBA) and DSM and orthomosaic creation.GPS tags were used for each image acquired, and there was no requirement for geo-location or orientation data.As the compact cameras generally used on board UAV platforms are extremely sensitive to temperature differences, vibrations and shocks, one of the processing steps was camera auto-calibration.This involved using information from each pixel of the images to estimate the optimal camera and lens calibration for each flight.
The imagery and the synchronized GPS position for each single image were used as input.The point cloud densification was set to high in order to maximize the density of 3D points of the digital elevation model, and the grid sampling distance in the digital elevation model point cloud (i.e., the distance between two 3D points in the DSM) was set to 100 cm.The orthomosaics cover a total of 138 ha with a spatial resolution of 0.083 m/pixel and 82 ha with a spatial resolution of 0.077 m/pixel for Study Areas 1 and 2, respectively.The camera parameters were optimized internally during the first step conducted.
The orthomosaic and DSM datasets obtained in this process were used in the subsequent analyses.Prior to further calculations, a 3 × 3 square kernel (approximately 0.25 m × 0.25 m) low-pass filter was applied to the original DSM in order to remove both tree crown and ground artifacts.
In addition, two additional parameters were computed for their integration in the crown classification.Hence, the Normalized Difference Vegetation Index (NDVI, Equation ( 1)) was calculated from the raw red and NIR band radiance values of the generated orthomosaic as a measure of the greenness of the vegetation cover.The difference from the minimum value (DifMin) was also derived from the original DSM using a search circle of 2 m around the target pixel according Equation 2: (2) where: hi = DSM value for the target pixel hmin = Minimum DSM value for a given neighborhood around the target pixel

Calculation of Tree Parameters
Tree height was retrieved from the DSM based on the identification of local maxima, likely to correspond to tree tops, and the measurement of their difference in height to a reference altitude of the ground in their neighborhood.The definition of the search window for the local maxima and ground reference varied according to the particular geometry and spatial pattern of the targeted trees for each case study, as described below.
For the IFAPA case study (discontinuous canopy), a square search window of 0.75 m was used for the identification of local maxima in the DSM corresponding to potential tree tops, and then, a circular neighborhood of 2 m was used for the calculation of the ground reference height as minimum DSM values.
For the Pedro Abad case study (continuous hedgerow canopy), a slightly different approach was used, due to the difficulties in identifying local maxima for trees in lines with closed crowns.Hence, and taking into account that each row was planted with 30 trees, the total length of each row was divided into 30 stripes of equal length.Then, the local maxima was computed for each stripe considering a rectangle of 1.80 m of across length, enough to capture the maxima of the tree height in the row, while avoiding being affected by the trees of the neighboring rows.The reference ground height for each tree was then computed as the local minimum in a circular neighborhood of a 3-m radius around the previously-identified tree top.
The crown diameter was also computed as a key variable in the estimation and prediction of olive tree productivity.Previous to diameter measurement, the tree crowns were delineated by means of an object-oriented classification using the software, eCognition Developer V 8.7 (®Trimble Germany Gmbh, Munich, Germany).We used this approach, as it was particularly suited to classification problems of high-resolution imagery where the target element is considerably larger than the image spatial resolution [38].In this object-oriented classification, the features to be classified are not single pixels, but spatial pixel clusters defined through an aggregation algorithm, known as "multi-resolution segmentation" [39].The segmentation process is controlled by the user by setting scale and homogeneity parameters.The classification was done on a multilayer image containing the three visible and IR multispectral layers of the orthoimage, the DSM along with the NDVI and the DifMin.
In the first stage, a simple level of segmentation was used, adjusting the parameters so that image objects would fit with real-world target elements.Segmentation parameter values were adjusted iteratively using a trial-error procedure, comparing segment edges against known reference areas (namely, different sections of tree crowns).A certain degree of over-segmentation was allowed in order to enable a good delineation of the elements in the image.All of the image data layers were equally weighted by one, and the scale parameter was set to 10 in the case of the IFAPA case study and to 15 in the Pedro Abad case study, due to the slight differences in the crown sizes and geometries.The shape parameter (i.e., shape vs. spectral importance, ranging from 0 to 1) was set to 0.2, and the compactness parameter (i.e., compactness vs. smoothness of segments, ranging from 0 to 1) was set to 0.5 for both study cases.
In the second stage, a set of training areas for the classes "tree crown", "background vegetation", "bare soil" and "shadow" was collected, and then, a supervised classification using the CaRT (classification and regression trees) algorithm was conducted.Finally, a classification refinement was done by applying a morphological opening filter, a threshold of object size and a watershed segmentation in order to delineate the individual tree crowns (Figure 3).Tree crown delineations were then exported as polygons, and the tree crown diameters were measured using the QGis V 2.2.0 (®Free Software Foundation, Inc., Boston, MA, USA) package.Measures were done following the same protocols as in the field measurements, i.e., cross measurements and averaging for the individual trees in the IFAPA case and width across the rows in the Pedro Abad case.

Tree Parameters Validation
In order to check the applicability of UAV imagery to the estimation of tree parameters, the first level of analyses were conducted at the individual tree level for the case study IFAPA, whereas analyses were done at the hedgerow level for the study case Pedro Abad.Hence, for each study case set, the residuals were calculated as the difference between actual individual tree (IFAPA) or hedgerow (Pedro Abad) parameters measured in the field and estimated from the DSM.The root mean square error (RSME) of these residuals was calculated along with the regression fit and the squared correlation coefficient between measured and estimated tree parameters.
As one of the objectives of this experiment was to assess the performance of UAV imagery products for the monitoring of tree biometrics differences between genotypes, a second level of analyses was conducted averaging the tree parameters per genotype.Therefore, the RMSE, regression fits and squared correlation coefficients (R 2 ) were also calculated at the genotype level for each study case.In addition, ANOVA analyses were conducted using the statistics package, R Studio V. 0.98.932(® 2009-2013 R Studio, Inc., Boston, MA, USA), in order to check the significance of differences between the tree parameters for the different genotypes assessed in each study case.These analyses were done separately on the measured and estimated datasets, including only the genotypes represented in the randomized complete block design (i.e., not including the guard rows and hedgerows) in order to compare if the two methods led to similar conclusions.

Tree Parameters Estimation
The application of image reconstruction by using SfM techniques on collections of imagery with high overlap allowed the generation of detailed DMS and orthomosaic for both study areas, as illustrated in Figure 4.
Tree height and crown diameter field measures and estimations by airborne imagery are summarized in Table 1 and Figure 5.The results showed a certain underestimation of both height and crown diameter in Case Study 1 (as shown by the mean and median statistics), whereas Case Study 2 showed similar values of measured and estimated tree heights and an overestimation of crown diameter.With the exception of tree height of Case Study 1, the estimations showed wider ranges of variation than the field measurements.
Linear fits (Figure 6) of the measured vs. estimation height at the individual tree/hedgerow level showed an R 2 = 0.07 in Case Study 1 (isolated trees), while in Case Study 2 (hedgerows), it reached 0.53.This situation was reversed in the case of crown diameters, were Study Case 1 showed higher R 2 (0.58) than Study Case 2 (0.22).RMSE also varied depending on the study case and variable, with relative errors (compared to the average value measured) ranging from 6.55% to 19.24% (Table 2).While the estimation of height and crown diameters maintained similar relative RMSE in the case of discontinuous canopy (Case Study 1), larger differences were recorded for the site with trees in hedgerows (Case Study 2).As expected from the R 2 values mentioned above, the results pointed out a higher agreement between the estimated and measured crown diameter than for the height in the case of discontinuous canopy, while results showed clearly the opposite in the case of trees in hedgerows.

Tree Parameters at Genotype Level
The data averaged by genotypes showed far better linear fits between field-measured and image-estimated parameters than the individual tree analyses, particularly in the case of height estimation in hedgerow plantations, with an R 2 around 0.9 (Figure 7).Crown diameter fits had lower R 2 for both discontinuous canopy and continuous hedgerows, with values of 0.66 and 0.53, respectively (the p-value of the latter being not significant), while the lowest R 2 (0.14) corresponded to the fit between measured and estimated height in discontinuous canopy genotypes.
RMSE showed a high agreement between the measures and estimations, particularly in the case of trees in hedgerows, with a value of 0.1 m corresponding to 3% of the average height measured (Table 3) and 0.16 m of RMSE and less than 8% of the average crown diameter measured.In the case of discontinuous canopies, the errors of estimation were less than 0.30 for both parameters estimated, meaning 11% and 16% of the average height and crown diameter measured, respectively.The ANOVA analyses checked whether similar conclusions could be drawn from either field measures or image estimations in the assessment of significant differences between genotypes.
In Case Study 1 (IFAPA, with discontinuous canopy) significant differences were found between genotypes for both field-measured and image-estimated crown diameters (Table 4).The distribution of sums of squares among sources of variation was almost the same in both measured and estimated parameters, with around 40%, 5% and 55%, due to differences between genotypes, blocks and residuals, respectively.However, significant differences between genotypes for heights were found only in the field measures in which differences between genotypes represent 41% of the sum of squares, whereas a value of 19.3% was obtained for the estimated height.crown diameter estimations and measurements in this case (7% and 13% of relative RMSE, respectively).Therefore, the relatively even architecture of the crowns along the hedgerows might have compensated for the difficulty in finding ground references for the height calculation.
In turn, the estimation of the variables in discontinuous canopies gave less accurate results (relative RMSE of around 19% for both height and crown diameter).In these particular cases, we also obtained low correlations between the estimated and measured values, pointing out the existence of certain limitations in the method at this level.One of the possible reasons for the errors in height estimation in this case could be due to some particularities of the experimental plantation.In fact, the short distance between tress in the rows and the presence of vegetation stripes between rows (cf. Figure 3) could have masked the real value of the ground values in addition to the abovementioned potential effect of a more even crown architecture of the trees in hedgerows.In any case, the height estimation accuracy obtained was comparable, yet in some cases, somewhat lower than in previous works on tree crops in patterns.For example, [11] reported an absolute RMSE between 0.3 and 0.4 m using a similar method; [44] obtained absolute RMSEs between 0.42 and 0.52 m using UAV imagery reconstruction and a pouring algorithm for the delineation of tree crowns and apex; and [7] reached an absolute RMSE of 0.4 m using a close-range Nadir ILRIS-3D sensor.
Results showed better accuracy for tree height estimation than for crown diameter in the case of hedgerows, both at the hedgerow and genotype level.In turn, this trend was not so clear in the case of discontinuous canopies, where only small differences in the accuracies of tree height and crown diameter were recorded, both at the isolated tree and genotype level.This difference between hedgerow height and diameter estimation could be due to the limitations of the remote sensing method.Yet, it could be argued, as we discuss below, that the inherent difficulties of the manual measurement of diameters across the hedgerows could compromise the validity of the reference data.It is also worth pointing out that such differences spans a few centimeters and that in experiments like the one conducted in the present work, the consideration or not of a single branch or shoot as indicative of the tree apex or crown diameter could explain these differences.In fact, some works consider that the tree height measurement is only possible to decimeter accuracy [44], sub-decimeter differences being difficult to attribute to the real accuracy of the estimation or to uncertainties in the reference data acquisition.
Estimations from imagery tend to show a slightly higher dispersion of data that could be linked to image artifacts due to acquisition or processing errors, but could also be linked to human factors in the field measurement of reference data that could mask out irregularities (protruding branches) in manual measurement.In fact, the nadir view of the tree crown in the height estimation from remote sensors could take into account tree irregularities ignored by the operator in field measurements [7].It is also worth pointing out that the present method allows the reconstruction of the whole tree or hedgerow crown projection, solving the difficulties of the manual measurement of this variable in the field [45], and it also provides a model of the 3D architecture of the upper part of the crown, which could be eventually used in modelling the whole crown shape.
Finally, the results also showed an overall agreement between the modelled and measured variables for both study cases in terms of finding statistically-significant differences between genotypes, pointing out a promising role of consumer-grade cameras on UAVs for field phenotyping purposes.In fact, the field and remote sensing methods disagreed only in the variable tree heights for discontinuous plantation, where field measurements showed significant differences between genotypes, while remote sensing estimations did not show such differences.For the other variables, the analysis of variance provided similar significant differences and distribution of sums of squares among sources of variation (genotype, block, residuals) in both measured and estimated parameters, which indicates that estimated results could be efficiently used for testing the genotype effect in olive breeding trials.

Conclusions
Breeding selection is in need of robust methods for the quantitative evaluation of a high number of genotypes and/or trials in an accurate way.In olive, as for other fruit trees, the characterization of the tree crown is a key issue to determine the adaptation of a genotype to a given growing system.In the present work, we tested the performance of a workflow based on low-cost UAV imagery for the estimation of tree crown parameters (namely, tree height and crown diameter) in the framework of olive tree breeding programs, both on discontinuous canopy (open vase configuration) and continuous canopy (hedgerow) orchards.The workflow involved image acquisition with consumer-grade cameras on board UAV platforms, orthomosaic and DSM generation using SfM image reconstruction and, eventually, GIS analyses and object-based classification for the calculation of tree parameters.
Comparison between reference field measurements and remote sensing estimation of crown parameters confirmed the performance of the workflow applied here as a quick and cost-effective alternative for the characterization of the olive tree crown, both in discontinuous canopy (open vase configuration) or continuous canopy (hedgerow).High accuracies were recorded at the individual hedgerow level (continuous canopy), both in height and crown diameter estimation.At the individual tree level (discontinuous canopy), accuracies were somehow lower, but in any case, satisfactory.These results were improved in the analyses at the genotype level for height and crown diameter, both on continuous and discontinuous canopy.With the exception of tree height for discontinuous canopy, the ANOVA analyses also showed the agreement between remote sensing estimation and field measurement in terms of the assessment of statistically-significant differences between genotypes in phenotyping studies, pointing out the interest in these data and techniques to be incorporated into the selection scheme in breeding programs.
It should be also noted that the image reconstruction processing was done without ground reference points, being a possibility for a potential improvement of the results.In any case, the ranges of disagreement between the remote sensing estimations and the field measurements were very close to the level of uncertainty in the field reference data acquisition, particularly in the case of analyses at the genotype level, and this fact could compromise the validation of further improvements in the methods.In this sense, further research could focus on the isolation of error sources, e.g., by adding an additional step in the validation as the manual measurement of the tree crown diameter, in order to identify and focus efforts on critical steps in the workflow (i.e., either on the orthomosaic and DSM production, on the image analysis or even in the reference data acquisition).Future research could also replicate the experiment using LIDAR data in order to compare its performance against the low-cost alternative of orthomosaics and DSM obtained by SfM image.

Figure 1 .
Figure 1.Study area and location of case study areas.Image source: Esri Inc.; PNOA, National Program of Aerial Orthophotograph-Ministerio de Fomento; IFAPA, Institute of Agricultural and Fishery Research and Training.

R
= pixel value of the red image band NIR = pixel value of the near infrared image band

Figure 3 .
Figure 3. Examples of crown delineation by object-oriented analysis in the Study Area 1 (Left) and Study Area 2 (Right).In yellow, the trees or hedgerow sections, depending on the case, measured in the field and used for the comparison against image estimation of tree parameters.

Figure 4 .
Figure 4. 3D blocks of the color infrared (CIR) orthomosaics draped on the DSM for two sections of Study Area 1 (Left) and Study Area 2 (Right).

Figure 5 .
Figure 5. Box and whisker plots of the measured and estimated tree variables for the two study sites at the individual tree level (150 trees for Case Study 1 and 28 hedgerows for Case Study 2).

Figure 7 .
Figure 7. Linear fits between measured and estimated tree variables for the two study sites at the genotype level (29 genotypes for Case Study 1 and six genotypes for Case Study 2).

Table 1 .
Statistics of the measured and estimated variables for the two study sites at the individual tree level (150 trees for Case Study 1 and 28 hedges for Case Study 2).H, tree height; CD, crown diameter.
Figure 6.Linear fits between measured and estimated variables for the two study sites at the individual tree (150 trees for Case Study 1) and hedgerow level (28 hedgerows for Case Study 2).

Table 2 .
Root mean square error (RMSE) and relative root mean square errors (RMSE %) computed as the percentage of the average of measured values of tree variables (Ave.meas.) for the two study sites at individual tree/hedgerow level (150 trees for Case Study 1 and 28 hedgerows for Case Study 2).

Table 3 .
Root mean square error (RMSE) and relative root mean square errors (RMSE %) computed as the percentage of the average of measured values of variables (Ave.meas.) for the two study sites at the genotype level (29 genotypes for Case Study 1 and 6 genotypes for Case Study 2).

Table 4 .
ANOVA table for the Study Case 1 (IFAPA) on the measured and estimation of height and crown diameter between genotypes.Df: degrees of freedom; Sum Sq: sum of squares.Significance: ** p < 0.01, *** p < 0.001.