Detecting Individual Tree Attributes and Multispectral Indices Using Unmanned Aerial Vehicles: Applications in a Pine Clonal Orchard

Modern forestry poses new challenges that space technologies can solve thanks to the advent of unmanned aerial vehicles (UAVs). This study proposes a methodology to extract tree-level characteristics using UAVs in a spatially distributed area of pine trees on a regular basis. Analysis included different vegetation indices estimated with a high-resolution orthomosaic. Statistically reliable results were found through a three-phase workflow consisting of image acquisition, canopy analysis, and validation with field measurements. Of the 117 trees in the field, 112 (95%) were detected by the algorithm, while height, area, and crown diameter were underestimated by 1.78 m, 7.58 m2, and 1.21 m, respectively. Individual tree attributes obtained from the UAV, such as total height (H) and the crown diameter (CD), made it possible to generate good allometric equations to infer the basal diameter (BD) and diameter at breast height (DBH), with R2 of 0.76 and 0.79, respectively. Multispectral indices were useful as tree vigor parameters, although the normalized-difference vegetation index (NDVI) was highlighted as the best proxy to monitor the phytosanitary condition of the orchard. Spatial variation in individual tree productivity suggests the differential management of ramets. The consistency of the results allows for its application in the field, including the complementation of spectral information that can be generated; the increase in accuracy and efficiency poses a path to modern inventories. However, the limitation for its application in forests of more complex structures is identified; therefore, further research is recommended.


Introduction
The world's forests offer a great diversity of goods and services that have led the scientific community to continuously look for new modern and efficient management schemes [1,2]. In fact, the vulnerability of forest ecosystems to climate change [3] and the anthropogenic pressure to which they are exposed [4,5] makes it necessary to better understand the dynamics of forest growth to strengthen decision making for their sustainability.
Given the ecological complexity of forest landscapes, several studies have been directed to promote innovation in forest inventories. In particular, the accurate and automated estimation of tree attributes has been a research trend in recent years. According to a recent study [6], such technologies represent a more efficient means of obtaining detailed information of individual trees, such as total height, diameter at breast height, basal area, tree count, crown delimitation, and volume estimation. For example, Gollob, et al. [7] documented a methodology for diameter estimation, diameter-at-breast-height (DBH) measurement, and tree detection in forest inventory. Lefsky, et al. [8] conducted an approach for the estimation of forest standing volume, obtaining a 95% confidence Under the hypothesis that UAV use significantly improves tree-level estimates regarding direct field measurements, this study proposes a method to estimate tree counts, total height, diameter, and crown area, including basal-diameter (BD) and diameter-at-breast-height (DBH) inferences on the individual-tree level in a pilot area in the north of Mexico. As a complementary objective, an estimation of spectral indices is shown as a strategy for determining the vigor of the clonal orchard.

Study Area
The study site is in the Sierra Madre Occidental in northern Mexico (27 • 8 57 N, 107 • 6 41 W; 2400 m asl) ( Figure 1). This region has a variety of ecosystems where pine and oak forests dominate, which resulted from the complexity of physiographic units and climates, and it is located at the transition between the Holarctic and Neotropical regions [36]. The dominant soils are Regosols and Leptosols with loam textures. The climate of the region is semicold, subhumid, with long fresh summers and monsoon rains that are combined with winter precipitation to reach an annual average of 779 mm, and an annual average temperature between 5 and 12 • C [17,37]. This region is one of the most important wood reserves in Mexico, and it is important for commercial timber production and forest conservation [29]. In addition to supplying more than 25% of the timber in Mexico, these forests preserve the main source of springs and streams that feed rivers and reservoirs, which are used for agricultural and public water supplies in the western states of Mexico.
This area is an important and valuable source of ecosystem services (e.g., water, biodiversity, and forest resources) to local populations, and it is a strategically important region for conservation studies. Particularly since 2005, a portion of 1.10 ha has been considered a clonal orchard area of Pinus arizonica Engelm. (Figure 1).
The pine clonal seed orchard was established in 2005 using clones obtained by transplants previously carried out in nurseries. The orchard design is regular, with a real frame arrangement with approximately 6 m spacing between ramets.

Workflow
The general process for the estimates of tree-level characteristics of individual trees using a UAV is illustrated in Figure 2. The first step was flight planning that involved defining the height, and the overlap between images and between lines for the subsequent execution of the flight and image capture. Images were processed and analyzed with photogrammetric procedures to derive the digital surface model (DSM) and digital terrain model (DTM) [35]. The canopy height model (CHM) was calculated from the difference in elevation between DSM and DTM. At the same time, multispectral vegetation indices were calculated as a strategy to test the phytosanitary status of the orchard. Canopy analysis involved detecting and geolocating trees in the clonal orchard, estimating heights, and delimiting crowns to obtain diameters and crown areas. Lastly, the estimation of tree-level characteristics calculated with the UAV was evaluated with respect to those recorded in the field by correlation analysis.

UAV Image Acquisition and Processing
On 17 October 2020, the study area was overflown using a DJI Phantom 4 Multispectral (P4M) quadcopter ( Figure 3a). The P4M camera has a total of 6 imaging sensors-5 multispectral sensors (i.e., blue, green, red, red-edge, and near-infrared bands) and 1 RGB sensor, all with a global 2 MP shutter ( Table 1). The focal length of the P4M camera is 5.74 mm, image size is 1600 × 1300 pixels, and sensor size is 4.87 × 3.96 mm [38]. The P4M demonstrated good accuracy, yield, and consistent data generation [39].  The flight plan for automatic image acquisition was programmed with the DJI Ground Station Pro application [40]. A 109 × 107 m double grid was established, covering the entire clonal orchard with an overlap of 80% between photos and 70% between lines. Flight height was 60 m above ground level (AGL) with a camera angle of 90 • and 75 • for the first and second grid, respectively (Figure 3b). The P4M operates under the principle of direct onboard georeferencing. In other words, it is equipped with a global positioning system (GPS) that geotags the coordinates of each acquired image. This indicates that the acquired images were directly georeferenced by the GPS capability during the flight mission [41]. It was not necessary to use the real-time-kinematic (RTK) system because the P4M manufacturer reported that the georeferencing system can reach vertical and horizontal location pressure of ±0.1 and ±0.3 m, respectively [42].
From the captured RGB and multispectral images, photogrammetric and computer-vision procedures were applied using free and open-source software OpenDroneMap (ODM) [43]. Outputs from ODM and the commercial options are similar; however, a major advantage is cost and its inherent feature of being an open-source software offers multiple options to customize photogrammetric procedures, decrease processing time, and increase documentation available to users. [44].
This software implements the modern Structure from Motion (SfM) and Multi-View Stereo (MVS) algorithms, showing common characteristics in images with a high percentage of overlap. It then generates a 3D point cloud of 1000-20,000 points m −2 [9] from which it filters and interpolates a DSM and DTM, which, in turn, are used to orthorectify each image and construct an orthomosaic [45,46]. This process generated the RGB orthomosaic, multispectral orthomosaic, DSM, and DTM of the clonal orchard. The raster calculator was used with free and open-source QGIS software [47], and the following equation was applied [48] to generate a canopy-height model (CHM) to represent the potential height of each tree: where CHM is the canopy-height model, DSM is the digital surface model, and DTM is the digital terrain model.

Estimates of Tree-Level Attributes
There are many methods for the detection and delimitation of individual tree crowns in different forest scenarios [49,50]. The ForestTools package [51] from statistical software R [52] was used in this study. This package was designed to work with CHM LiDAR data or photogrammetric point cloud derived from UAV photogrammetry. In short, the package automatically detects tree crowns, obtains tree height (m), generates polygons, and calculates the crown area (m 2 ). From there, the following formula was used to find the crown diameter (m) of each tree present in the clonal orchard: where d is the crown diameter and a is the crown area. Initially, a Gaussian filter was applied to smooth the CHM to remove heterogeneity in pixel values [53]. The size of the filter was determine by setting the standard deviation parameter (0.75), which is in units of grid cells using free and open-source WhiteboxTools software [54]. Subsequently, the ForestTools package was implemented as a tool to individually geolocate trees and delimit tree crowns through the variable-window-filter (vwf) algorithm [55] and the marker-controlled segmentation algorithm proposed by Beucher and Meyer [56], respectively.
Field data were collected on the same day of the flight (17 October 2020) and correspond to the beginning of the species' dormancy stage [17]. Each individual tree in the clonal orchard was detected, manually tagged, and rigorously measured. Basal diameter (BD; cm) and diameter at breast height (DBH; cm) were recorded using a diameter tape. Tree heights (m) were recorded using a level staff. Average crown diameter (m) was also obtained by recording the width of the horizontal crown projections (east-west and north-south) for each tree using a measuring tape. Crown area (m 2 ) was obtained from the average crown diameter, and the following formula was applied: where a is the crown area and d is the average crown diameter. However, crown-area calculation assumes a circular geometric pattern, which in the field is usually irregular in shape and spatial extension. For this reason, delimitation of this attribute by means of a UAV recovers relevance, since it is considered to be a better approximation to the precise spatial conformation of the individual.
To compare the estimates of tree-level characteristics obtained with the UAV with respect to field measurements, we used correlation analysis to estimate the degree of correspondence between estimated variables and those recorded with field data. Moreover, simple regression models were used to estimate adjustment efficiency regarding the determination coefficient (R 2 ) and root-mean-square-error (RMSE) by means of statistical software R [52]. Furthermore, the fulfillment of regression assumptions [57] used in forest modeling was verified. Attributes were classified to increase the goodness of fit and explore the potential improvements of statistical parameters [58].
Given the goodness of estimation of certain variables from the UAV, using allometric relations equations were obtained to predict the BD and the DBH. From a practical point of view, this strategy allows for enhancing the use of individual tree-level attribute measurements, saving time and improving certainty.

Multispectral Image Analysis
Since vegetation indices are proxies of vegetation vigor or productive status [13,14], we also tested multispectral vegetation indices recorded from the UAV in the study area. These indices were estimated from the multispectral orthomosaic with the QGIS raster calculator (Equations (4)-(11)): where NDVI, the normalized difference vegetation index that is the most used indicator of the chlorophyll content in vegetation, provided information about vegetation growth and nutrients. GNDVI, green NDVI, uses a green band, which represents the surface, instead of the red used in NDVI. GNDVI can show whether vegetation is lacking water or nutrients or is suffering from postmaturity biomass recession. LCI, leaf chlorophyll content, is an important index for evaluating vegetation growth and yield. It is also an indicator for evaluating plant nutrient stress, disease, growth, and aging. NDRE, normalized difference red-edge index, uses a red-edge band instead of the red used in NDVI. This red-edge band is a spectral region in the transition zone from the red spectrum to the near-infrared spectrum. NDRE

Results
A total of 2340 images were captured. Of these, 390 images were from the RGB sensor. Figure 4 shows the RGB orthomosaic, DMT, DSM, and CHM derived from the photogrammetric process with ODM. The coverage area of these products is 1.10 ha and has a spatial resolution [60] or ground sampling distance (GSD) of 4 cm. In the field study, 117 trees in the clonal orchard were recorded and measured. Table 2 shows the descriptive statistics of tree-level characteristics recorded in the field. Using the variable-window-filter algorithm of the ForestTools package, 95% of the trees (i.e., 112 out of 117) were detected, and their tree crowns were delimited. Although the trees are the same age, the phenotypic conformation is variable [17], and their open and regular canopy ( Figure 5) was strategic for the purposes of the study.  Statistical evaluations of tree-level characteristics recorded with field data and those estimated with UAV are summarized in Table 3. Height, area, and crown diameter were underestimated by the UAV by 1.78 m, 7.58 m 2 , and 1.21 m, respectively. High and positive correlation was observed between height, area, and crown diameter recorded with field data and those estimated with the UAV. That is, the values of correlation coefficients (r) are close to 1.0 (0.97, 0.95, and 0.95, respectively). In general, RMSE and bias values indicated high consistency between measurements of tree-level characteristics. Although class rankings were made as suggested by Wang, Lehtomäki, Liang, Pyörälä, Kukko, Jaakkola, Liu, Feng, Chen, and Hyyppä [58], no better estimates were found than those obtained in total, as presented below.  Figure 6 shows the results of linear regression applied in order to determine the relationship between data recorded with field data and those estimated with UAV; a graph of residues against predicted values (i.e., height, area, and crown diameter) is also included. For height recorded in the field, corresponding to the 112 trees detected, the model was able to explain 95% of variability observed from that tree-level characteristic (R 2 = 0.95) with an RMSE of 0.36 m (Figure 7). Similarly, the linear model for the crown-area variable had a determination coefficient of 90% (R 2 = 0.90) and an RMSE of 3.88 m 2 (Figures 6b and 7). Lastly, the model for the crown diameter was able to explain variability observed in 91% of that tree-level characteristic (R 2 = 0.91) with an RMSE of 0.47 m (Figures 6c and 7).
For the inference of BD and DBH, all attributes calculated with the UAV (i.e., height, crown area, and crown diameter) were tested. The best-fitting predictor variables were determined by the following equations: where DBH is the diameter at breast height and H is the height estimate using a UAV.
where BD is the basal diameter and CD is the crown-diameter estimate using a UAV. The DBH and BD models showed favorable goodness of fit. The estimated height with UAV was the predictor variable that better explained DBH variability (R 2 = 79%, p value 2.2 × 10 −16 ) and an RMSE of 2.81 cm, while BD was better inferred by CD, explaining the variability of this parameter in 76% (R 2 ) with an RMSE of 2.73 cm. The p value of both models was significant (2.2 × 10 −16 ), which confirmed a statistical relationship (not random). The statistical significance of the coefficients of the independent variables is shown in Table 4.    Figures 8 and 9 show the scatter diagrams between each predictor and residues for each model. Therefore, it was possible to check the condition of the linear relationship between numerical predictors and response variable (i.e., if the relationship were linear, residuals would be randomly distributed around 0 with constant variability along the X axis).  Regarding vegetation-index estimates, 1950 images were captured with P4M's multispectral sensors (i.e., blue, green, red, red-edge, and near-infrared bands). Figure 9 shows the different vegetation indices calculated for the clonal orchard from the multispectral orthomosaic derived from photogrammetric processing with ODM. Moreover, Figure 10 shows the distribution of the average values of vegetation indices corresponding to tree crowns. The estimated indices were NDVI, TVI, LCI, GNDVI, RVI, OSAVI, NDRE, and NDGI; Table 5 shows the descriptive statistics of the values of these indices.  Tree-crown vegetation-index values were visually recognizable on the higher levels (i.e., green and yellow shades; Figure 9). Difference was evident between tree crowns and surfaces different from tree crowns (i.e., orange and red shades), corresponding to bare soil or provided with herbaceous vegetation. However, for all indices, high and low values could be recognized within the tree crowns, attributable to different levels of health or stress. Consequently, each index had a different histogram shape ( Figure 10). According to Minařík and Langhammer (2016), the divergent shape of histograms implies the different sensitivity of the indices to separate categories of forest disturbance and/or productivity, highlighting NDVI as the one with the best capacity.
Regarding  Figure 10 shows that this index is less contrasting to define tree crowns and highlights the photosynthetic activity in the surrounding areas with the presence of herbaceous vegetation.

Discussion
To our knowledge, this is the first approach to estimate tree-level attributes and multispectral indices using UAV in a pine clonal orchard; although it is an area with trees regularly distributed in space and of similar age, findings may lead to improving the context of current inventories.
In fact, the recent literature reported exponential interest of UAV applications in the estimation of structural attributes during the last five years [6], and contributes to improving those technologies with limitations in temporality and spatial resolution that complicate more realistic mapping in small areas [61]. Thus, our findings in terms of the calculation of height, crown attributes, and tree counts estimated using UAV offer a promising future in forest-resource management towards precision forestry [2,6,62].
Height estimation is a typical parameter in forest management since it is a basic input for calculating volume, and it is usually a parameter of forest productivity, including carbon stocks [63]. Although the study orchard was directly and rigorously telescoped (some cases involving tree climbing), accuracy and precision were limited to higher heights [58]. Optionally, our remote-sensing methodology offered a solid estimate (Table 3); however, it underestimated by 1.78 m, which could be attributed to the limited generation in the cloud of points at the tree apex. According to Krause, Sanders, Mund, and Greve [63], the photogrammetric dataset implicitly had this error, and it is considered acceptable and often attributed to the leaning angle. Furthermore, according to the standardization of typical forest inventories, an error of 10% is allowed, which is valid for practical purposes [64].
Height measurements by a technician in the field are usually performed with indirect methods where clear visibility is required from the base to the top of the tree; despite the technician's attempts to control this error, under natural conditions visibility can be restricted due to the complexity of the shape of the crown, the occlusion of neighboring trees, very steep slopes, etc. Therefore, the implication of errors in the field could be greater than that obtained by the UAV (Figure 7), so our approach is viable in operational terms. However, to significantly improve accuracy, Krause, Sanders, Mund, and Greve [63] recommended repeating the missions with the consequent gain in data quality and quantity, especially in new and fast-growing species. Likewise, georeferencing individuals to an acceptable level of precision gives them a strategic advantage over expensive and time-consuming fieldwork.
Comparatively, the results of height estimation can be considered acceptable in their estimates with those of Wu, et al. [65] (R 2 = 0.89, RMSE = 0.19 m, rRMSE = 5.37%); Roşca, et al. [66] (RMSE 0.11-0.63 m); Krause, Sanders, Mund and Greve [63] (RMSE values of 0.304 m (1.82%) and 0.34 m (2.07%), respectively (n = 34)); and Falkowski, et al. [67] measurements of tree height (r = 0.97) and crown diameter (r = 0.86), albeit in different species and under different environmental conditions, so further research should be conducted as suggested on variations in slopes, densities, complexity, and crown classes [58]. Although the same author reported that the results were improved by the classification of attributes, our study found no benefit when classifying them (R and RMES values were lower).
With respect to contour and crown-area measurements, these estimates are more important as tree-level characteristics that can be used to estimate other attributes such as volume, biomass, and carbon. Wu, Johansen, Phinn, Robson and Tu [65] highlighted the importance of this characteristic to design better cultural practices on the individual-tree level, detect physiological activity, and for phenological monitoring (production, cone quality, pruning). Above all, its metric input is potentially strategic for precision forestry [2,12]. Assessing the absolute accuracy of UAV-based results requires independent and alternative measures. However, for the measurements of diameter and crown area, it was not possible to have real data because of irregular geometrical tree characteristics. Real field measurements are laborious and inaccurate because they assume a circular geometric model, unlike the more accurate results offered by the aerial perspective where tree crowns are favorably visualized according to their exact shape and spatial extension.
Therefore, to achieve a more accurate estimation of crown area and crown diameter, it is advisable to take multiple measurements with the UAV until standard deviations are minimized to the desired level. Other mission parameters must be considered to find the optimal one (e.g., flight heights) [68]. Compared to the study by Johansen, Raharjo and McCabe [68], the GEOBIA algorithm and eCognition Developer 9.2 software found an R 2 = 0.96 and RMSE of 0.6 m, although they tested different heights, and the crown area was not measured in the field. In addition, in such a study, mapping was improved in the youngest (smallest) trees, which showed differential estimates according to tree dimensions; however, as mentioned in this study, such classification did not yield statistical improvements.
Despite the promising results (Table 3 and Figure 6), crown delineation is complicated in trees that occlude neighboring trees with their branches, so caution is required in applying this approach with those dense stands [58]. It is also recommended to complement with other technologies such as those of 2D spatial-wavelet analysis by Falkowski, Smith, Hudak, Gessler, Vierling, and Crookston [67].
Although analysis incorporating two or more species is beyond the objectives of this research, we expect that high tree densities would make it difficult to model, at least with these sensors, since tree-crown measurements in old growth stands are labor-intensive and imprecise because of the difficulty of observing them from the ground. In this case, sensors such as terrestrial laser scanning that can cross the canopy are recommended. According to Ritter and Nothdurft [69], these technologies have proven their efficiency in other geographical areas.
In relation to tree count, the potential usefulness of the developed algorithm is evident; the high rate of found precision (112/117), despite phenotypical differences in individuals, reflects the relevance of the followed approach. This coincides with recent findings generated in other latitudes, such as the case of Jiménez-Brenes, López-Granados, de Castro, Torres-Sánchez, Serrano, and Peña [27], which reported 80% of the trees in an olive orchard. It is recommended to research in more detail the configuration of different parameters of the flight [68], which could yield better tree-counting efficiency. On the other hand, the large amount of travel time required for individual field counts is drastically reduced by the procedure followed here. This gives a high degree of practical application to the results towards aerial inventories [9], which, combined with allometric equations, could calibrate forest-growth models, forest management, and studies of forest ecology.
Thus, our workflow is solid and important, particularly for monitoring small sites on a continuous and permanent basis, including the complementation of spectral information that could be generated. For example, tree-level attributed results were obtained only by RGB imaging (Figure 4a), while multispectral sensors (i.e., blue, green, red, red-edge, and near-infrared bands) augmented the data quality of the results. Taken together, the derived products are revolutionizing forest monitoring.
The use of UAVs has been increasing, and multiple innovation choices are projected for forest monitoring and assessing tree attributes [6,58]. The capacity for capturing high-resolution imagery, even in the presence of clouds, comparatively enhances the quality of data with respect to traditional Earth-observing satellites. Hence, research perspectives are expanded to tree-level attribute detection and forest-health monitoring. Particularly, the emergence of modern sensors and linked platforms are stand-out opportunities to improve the knowledge of forest dynamics.
By combining remotely sensed data products from UAV with other related disciplines, it is possible to more accurately detect tree-level characteristics, including tree physiological stress, due to the high sensibility of multispectral sensors. Furthermore, the cost of UAVs and their associated platforms is rapidly decreasing, while the number of sophistication and flexibilization procedures is increasing [6]. Opportunities offered by processing methods indicate significant progress to augment analysis framework capabilities. For instance, free and open-source software platforms are becoming popular and provide alternatives for the estimation of tree-level attributes, such as above-ground biomass and stem volume [44].
However, there is still the need to improve the distortion of camera calibration. Some authors suggested to increase the checkpoints using on-board differential GPS systems, such as real-time-kinematic (RTK) or postprocessing-kinematic (PPK) technologies, although in practical applications, they are costly, and most forestry applications on this geographic scale are enough, as shown here [19]. Results also cannot be extrapolated to forests with more complex structures, such as unevenly aged stands. They are different from those studied here, where uniformity and spatial arrangement allowed for the easy distinction of individuals. Therefore, the next challenge is to continue with further research in natural forests that are representative of forest ecosystems that are under some kind of management, including areas of conservation and protection.
The application of vegetation indices represents a significant contribution towards precision forestry. For example, spatial distribution in Figure 9 describes local variations in vigor that an orchard manager cannot detect on the field. This spatial variation in individual productivity suggests a tool for growers in terms of differential individual-tree management. Similar results were documented by Grulke, et al. [70].
The values of vegetation indices showed a different histogram shape ( Figure 10). According to Zhao, et al. [71], the divergent form of histograms reflects the different sensitivity of indices in terms of separating categories of forest disturbance and/or productivity, highlighting NDVI as the index with the greatest capacity. This is in line with the results of Brovkina, et al. [72], who considered NDVI to be a good vigor indicator.
GNDVI values were higher than those reported by Wahab, et al. [73], which may indicate greater photosynthetic activity, so it is recommended as a parameter for monitoring variations in leaf chlorophyll and vegetation health [74]. In addition, the OSAVI threshold for tree crowns in the clonal orchard suggested that green aerial biomass stored was greater than that of the subsoil [75]. According to Minařík and Langhammer [76], NDRE results suggested that P. arizonica crowns are consistent with mostly healthy trees in Picea abies forests. However, our RVI results were lower than those reported by Xing, et al. [77], which is probably attributed to differences in the used sensors since, to our knowledge, this is the first study that generated such values from a UAV; consequently, comparisons with others studies lack contrasting parameters.
We recognize the limitations of vegetation-index approaches, such as phytosanitary parameters. For example, although NDVI was the best vigor proxy, it is necessary to improve its discrimination values; grass or weed residues can decrease the predictive capacity of NDVI, and crown size and superposition can reduce the use of remote sensing on the tree level [78,79]. In this case, soil-management practices such as drainage, tillage, and traffic control are recommended, with implicit benefits for tree development. Another difficulty is related to the problem of estimating the differentiated temporal rates of physiological and phenological processes in the ramets. NDVI values reflect orchard vigor during the evaluation period from a static point of view. Consequently, further research is needed to disentangle the responses of multispectral index values on temporal scales.

Conclusions
The major goal of this research was to extract tree-level attributes from a UAV. Estimates of total height, diameter, and crown surface were consistent with those reported in the literature. Furthermore, allometric relationships proved to be efficient in predicting diameter at breast height and basal diameter from attributes calculated with the UAV. The diversity of the used indices described the phytosanitary status of the orchard, with NDVI as the most recommended to evaluate productivity without discarding periodic evaluations. In this way, our hypothesis was confirmed; these measurements are reliable, and contribute to making solid and better-informed decisions on the individual-tree level. The application of the proposed workflow demonstrated solid parameters in estimates with consistent time saving and effort. Consistency in the measurement of tree-level characteristics on the individual-tree level has potential applications to improve current inventories. Nevertheless, recognizing its scopes and limitations, it is advisable to be cautious with the application of this approach in areas of greater complexity, such as the one studied here. Therefore, it is desirable to test different flight configurations and repeat missions to optimize the results and enrich the usefulness of the data.