Skip to main content

The utility of terrestrial photogrammetry for assessment of tree volume and taper in boreal mixedwood forests


Key Message

This study showed that digital terrestrial photogrammetry is able to produce accurate estimates of stem volume and diameter across a range of species and tree sizes that showed strong correspondence when compared with traditional inventory techniques. This paper demonstrates the utility of the technology for characterizing trees in complex habitats such as boreal mixedwood forests.


Accurate knowledge of tree stem taper and volume are key components of forest inventories to manage and study forest resources. Recent developments have seen the increasing use of ground-based point clouds, including from digital terrestrial photogrammetry (DTP), to provide accurate estimates of these key forest attributes.


In this study, we evaluated the utility of DTP based on a small set of photos (12 per tree) for estimating stem volume and taper on a set of 15 trees from 6 different species (Populus tremuloides, Picea glauca, Pinus contorta latifolia, Betula papyrifera, Picea mariana, Abies balsamea) in a boreal mixedwood forest in Alberta, Canada.


We constructed accurate photogrammetric point clouds and derived taper and volume from three point cloud–based methods, which were then compared with estimates from conventional, field-based measurements. All methods were evaluated for their accuracy based on field-measured taper and volume of felled trees.


Of the methods tested, we found that the point cloud–derived diameters in a taper curve matching approach performed the best at estimating diameters at the lowest parts of the stem (< 30% of total tree height), while using known DBH and height provided more accurate estimates for the upper parts of the stem (> 50% of total height). Using the field-measured DBH and height as inputs to calculate stem volume yielded the most accurate predictions; however, these were not significantly different from the best point cloud-based estimates.


The methodology confirmed that using a small set of photographs provided accurate estimates of individual tree DBH, taper, and volume across a range of species and size gradients (10.8–40.4 cm DBH).

1 Introduction

The key to sustainably managing the world’s 4 billion hectares of forested area (Bahamondez et al. 2010) is undertaking detailed and accurate inventories at varying scales (Gillis et al. 2005). Forest inventories are carried out at both operational and strategic levels, providing insight to short-term harvesting strategies or long-term environmental management, respectively (Wulder et al. 2008). The development of inventory techniques to provide robust and reliable information is critical for forest practitioners to understand dynamic forest ecosystems in changing and uncertain environmental conditions. This is particularly important in the boreal mixedwood, which represents a major component of forest habitat in the northern hemisphere, and whose structural and compositional diversity make it resilient to a variety of disturbances (Chen and Popadiouk 2002).

While some components of a forest inventory are distribution-based, such as species mixtures and age, the size of individual tree stems remains a critical measurement that provides the basis for the remainder of a forest inventory. One such measurement is the diameter at breast height (DBH), which is a key input to allometric equations calculating, for example, volume (Huang 1994) or biomass (Lambert et al. 2005). A relative frequency distribution of DBHs in a given area yields a stem size distribution (Taubert et al. 2013), which can be used directly to estimate forest attributes such as its structure, successional stage, or volume (Gobakken and Næsset 2004; Hetemäki et al. 2010; Nduwayezu et al. 2015). These inventory attributes provide valuable insights to inform stand- and landscape-level forest management decisions; therefore, their complete and accurate estimation is critical for maximizing a range of both economic and ecological forest values.

Many forest inventory attributes are manually measured as a part of a ground-based inventory, where all, or a sample, of the trees above a given diameter threshold in a sample area are measured. However, traditional inventory methods of measuring tree stem taper is difficult and typically requires felling the tree. Additionally, for trees with unconventional stem shapes resulting from varying growth patterns or environmental conditions, traditional ground-based inventory methods and equations based on diameter and height as inputs may fail to provide accurate estimates of tree volume or biomass.

In order to meet increasing demands of data accuracy and robustness, recent years have seen the incorporation of remote-sensing technologies to enhance forest inventories (White et al. 2016; Leckie and Gillis 1995; Wulder and Franklin 2003). One such technology is Terrestrial Laser Scanning (TLS), which uses Light Detection and Ranging, or LiDAR, from a ground-based sensor to more effectively characterize individual stems at a plot or individual tree scale, providing accurate estimates of tree DBH, height, stem volume, and stem biomass (Liang et al. 2016). However, TLS units are expensive and often unwieldy (Eitel et al. 2013). A faster and inexpensive alternative to derive similar data for use in forest inventory is digital photogrammetry, which creates point clouds using images taken at multiple locations. Photogrammetric point clouds derived from airborne imagery, called digital aerial photogrammetry (DAP), are used to create highly detailed surface models of the forest canopy at broad scales, but typically cannot return points from under the surface of the canopy (Tao et al. 2011; White et al. 2015). Photogrammetric point clouds can also be derived from ground-based imagery, covering a smaller area than that of DAP, but providing a much higher level of detail at the individual tree level. Recent research has shown the success of digital terrestrial photogrammetry (DTP) in the reconstruction of individual trees for attributes such as DBH (Forsman et al. 2016), location within a plot (Liang et al. 2014a), and stem shape (Bauwens et al. 2017).

Despite advancements in DTP technology and associated methods, there are current limits to its operational use. The need for manual intervention or trial and error in point cloud generation has limited the application of DTP within forest management and ecological modeling, with advances in automation of point cloud generation, an important way forward for making the technology useful as an operational tool (Mikita et al. 2016; Berveglieri et al. 2014; Hapca et al. 2007). Additionally, previous studies using DTP have relied on either the acquisition of hundreds to thousands of images over a given area (Mokroš et al. 2018), or tens to hundreds of images of single trees (Bauwens et al. 2017; Miller et al. 2015), thereby raising issues of time or data storage requirements in operational capacities. Point cloud–derived upper stem measurements, such as those from Fang and Strimbu (2017), could provide better estimates of attributes such as taper, volume, or biomass (Bauwens et al. 2017). However, most studies focus on relatively even-sized stands or a single primary species (Fang and Strimbu 2017), and, as a result, an additional analysis of point cloud accuracy across species and environmental gradients is needed to understand the utility of DTP in irregular stands such as those present in boreal mixedwood forests.

In this study, we evaluate a methodology for DTP estimation of DBH, upper stem diameter (> 1.3 m), and volume of individual trees in a boreal mixedwood forest. Limited sets of photographs taken at known locations were used to automatically generate photogrammetric point clouds for trees across a range of sizes and species. Estimates of diameters at varying heights were derived from the point clouds and used as inputs to estimation of taper and volume. The accuracy of DTP-derived estimates was assessed based on field-measured taper from felled trees and compared with traditional methods (based on height-diameter allometries) for estimation of taper and tree volume.

2 Materials and methods

2.1 Study area

The study area is a boreal mixedwood forest near the towns of Swan Hills and Slave Lake in Alberta, Canada. The study area ranges from approximately 54.97–55.26° N and approximately 115.08–115.58° W. The 700,000-ha forest management area is primarily used for timber harvesting and oil and gas extraction. Ten tree species are present, with the most common being white spruce (Picea glauca), black spruce (Picea mariana), trembling aspen (Populus tremuloides), and lodgepole pine (Pinus contorta). The region receives approximately 600 mm of annual precipitation and has mean summer and winter temperatures of 20 °C and − 21 °C, respectively (Natural Regions Committee 2006). Elevation in the study area ranges from 545 to 1575 m above sea level.

2.2 Field data

The trees measured in this study were located within plots established in July 2018. Field plots were part of a systematic random sample capturing the height and species ranges of stands present. Within plots, trees ≥ 7 cm DBH were measured, including species, height, height to live crown, crown class, and DBH. Across the ten sample plots, 15 individual trees were randomly selected and photographed as outlined in Section 2.3. Trees were later felled, and diameters were measured in 1-m increments up the stem. Field-measured volume was calculated as the sum of sections in between diameter measurements up each stem. Attributes for each tree are shown in Table 1.

Table 1 Characteristics of trees (n = 15) used in this study. Species codes are as follows: Aw trembling aspen (Populus tremuloides), Sw white spruce (Picea glauca), Pl lodgepole pine (Pinus contorta latifolia), Bw paper birch (Betula papyrifera), Sb black spruce (Picea mariana), Fb balsam fir (Abies balsamea)

2.3 Image acquisition

Before images were taken, five coded targets were positioned on and around the tree (four on the ground and one on a tree approximately 2 m high). In some cases, targets were placed on the selected tree but were filtered from resulting point clouds. Targets were generated by Agisoft Photoscan (Agisoft 2018) for automatic detection during image matching, and each target took up approximately the width of a standard letter-sized page (22 × 28 cm). The targets were used to enhance image alignment, both within and between image locations.

Figure 1 details the image acquisition method. Images were acquired using two RICOH Theta S (RICOH 2017) cameras mounted on a telescoping pole. Each camera is equipped with two fisheye lenses whose images are stitched together to generate a single spherical image with a 360° field of view. The mount ensured the two cameras remained at a fixed distance apart (70.0 cm), which allowed this distance to be input as a scale bar during point cloud processing, outlined in Section 2.4. Sets of simultaneous images, acquired with both cameras, were taken at each of three heights, approximately 2, 3, and 5 m above the ground, and at each of two locations around the tree. The result was a set of 12 images (two cameras, three heights above the ground, two locations). Using two adjacent cameras (placed ~ 70 cm from each other) allowed for high (close to 100%) overlap between image pairs. Sets of images taken at two locations meant that approximately half of the circumference of the tree was visible in the set of images. Preliminary testing showed that this methodology provided sufficient coverage around the tree for circle-fitting techniques to accurately estimate stem diameter. The locations of camera positions and coded targets were recorded relative to ground level at the location of the first image set, which was set as a center of a local coordinate system (with xyz coordinates of 0,0,0). Table 2 compares the methodology presented in this study to that of previous work producing photogrammetric point clouds from ground-based images.

Fig. 1
figure 1

Image acquisition methodology used in this study. a Outlines of the camera mount setup. b Details of the approximate location of photos and coded targets (GCPs) in relation to the tree. c The camera in use during field acquisition

Table 2 A comparison of selected studies detailing previous work with DTP to that of the proposed methodology. Single asterisk indicates values not explicitly reported but calculated based on values presented in study. Double asterisk indicates values not reported but estimated based on figures in study. NR not explicitly reported

2.4 Point cloud generation

Point clouds were processed using an automated Agisoft Photoscan workflow (Agisoft 2018). First, camera and target locations were entered and targets were automatically detected. Next, a scale bar between each set of photographs was set as the distance between the images taken in the field (~ 70 cm). Photos were aligned using a “high” setting, and the resulting tie points were filtered to remove those with high reconstruction uncertainty. Finally, dense point clouds were generated using a “high” setting and were then exported to be used in further processing. Processing was performed on a computer with an Intel Xenon E5-2630 (24 cores @ 2.3 GHz), 64 GB of DDR3 RAM, and an NVIDIA Quadro P4000 GPU.

2.5 Attribute extraction

Generated point clouds were analyzed in Computree (Piboule et al. 2015), a collaborative and open-source software to derive detailed tree-level estimates from ground-based point clouds. Analysis followed a standard Computree workflow, beginning with detection and removal of ground points, followed by noise removal (Belton et al. 2013). In the next step, horizontal clustering of points and vertical aggregation into logs was performed. For resulting logs, cylinders were fit at various heights up the stem using a least squares fitting technique (e.g., Berveglieri et al. 2017). Finally, smoothed diameters were calculated at breast height and 1-m intervals by averaging diameters of neighboring cylinders (e.g., 1.2–1.4 m for DBH).

2.6 Taper and volume assessment

As heights of the point cloud measurements did not reach the full height of the stem, estimates of taper were determined by matching point cloud–derived diameters to a database of all possible taper curves for the area. In Alberta, variable-exponent taper equations are used (Kozak 1988), and parameters of these curves are adapted to various ecoregions of the province (Huang 1994). Generation of the curve database and associated matching techniques are described below.

2.6.1 Curve database

Taper curves were generated for all possible tree dimensions in the study area based on parameters used throughout the province (Huang 1994). To generate each curve, taper equations required inputs of species, ecoregion, DBH, and total height, while outputting the diameter of the tree at any given height. All possible taper curves were thereby created using all possible combinations of the variables in our study area—all three ecoregions, six species, DBH values from 4 to 40 cm in 0.1-cm increments, and height values from 5 to 35 m, in 0.1-m increments. For each combination, equations output the diameter values at height increments of 10 cm up the stem. The list of curves was then filtered to remove trees whose allometries were unlikely to exist in our study area (e.g., trees that were 30 m tall and had a 4 cm DBH), by removing curves from the database whose height values were not within ± 5 m of the predicted height from the specified allometric equation. This limited the database to allometrically valid taper curves (e.g., those that could exist in the study area) and resulted in a final database of 1,652,778 taper curves.

2.6.2 Curve matching

Point cloud measurements of diameter at different vertical heights were used to match to possible taper curves. Two different curve matching approaches were evaluated in this study, outlined in Fig. 2. In the first method (1), diameters were not weighted and the curve chosen based on having the smallest residual between the point cloud–derived diameters and diameters from possible taper curves. The second method (2) applied weighting factors to the residuals of diameters closest to 3.28 m. This 3.28 m height was chosen as it was the average of the three camera heights, and all cameras are expected to have the lowest residual distance to this point on the stem, potentially making it the portion of the stem most accurately reconstructed by the point clouds.

Fig. 2
figure 2

A representation of the two different curve matching techniques (not weighted and weighted) used in this study. w in the equations indicates the weight applied to the residual closest to 3.28 m. In this simplified example, the blue and yellow lines represent two candidate curves coming from the generated taper curve database. In the curve database, however, there were approximately 90,000 candidate curves for each species in each ecoregion

2.6.3 Verification and accuracy assessment

As a comparison, estimates from the point cloud curve matching were compared with three other methods: (3) the selection of a taper curve based on measured DBH, height (from laser hypsometer), and species before trees were felled; (4) field-measured DBH and species were used as inputs to height-diameter allometric equations (Huang 1994) to predict a tree height, which was then input to select a taper curve; (5) the fifth dataset was similar to the fourth, but used DBH as estimated from the point cloud as an input to a height-diameter equation to predict height. In summary, five methods were compared, two of which used field measurements to match a taper equation (3 and 4), and three which were based on point cloud estimates of diameter (1, 2, and 5). An overview of these methods is shown in Fig. 3. The accuracies of the DBH, volume, and taper were evaluated by using the root-mean-squared error (RMSE), RMSE relative to the mean (RMSE%), bias, and bias relative to the mean (bias%). Results of the methods were tested to see if estimates differed significantly from one another using a t test. The equations for these statistics are as follows:

$$ RMSE=\sqrt{\frac{1}{N}\sum \limits_{i=1}^N{\left(\hat{y_i}-{y}_i\right)}^2}, $$
Fig. 3
figure 3

An overview of the five methods used: curve matching (b-c.; methods 1 and 2), inputting field measurements (d-e.; methods 3 and 4), and using the DBH as estimated from the point clouds (f.; method 5)

$$ RMSE\%=\frac{RMSE}{\overline{y}}\times 100, $$
$$ \mathrm{bias}=\frac{1}{N}\sum \limits_{i=1}^N\left({\hat{y}}_i-{y}_i\right), $$
$$ \mathrm{bias}\%=\frac{\mathrm{bias}}{\overline{y}}\times 1, $$

where N is the number of trees, yi is the reference measurement for tree i, \( {\hat{y}}_i \) is the predicted measurement for tree i, and \( \overline{y} \) is the mean of reference measurements on all trees

Data availability

An example point cloud showing DTP reconstruction of tree 11 can be found in the ResearchGate repository (Mulverhill et al. 2019) at Additionally, a web viewer showing this tree can be found at

3 Results

3.1 Point cloud reconstruction and diameter estimates

Once an efficient processing workflow was produced, the total amount of time taken for all steps was, on average, 8 min per tree (3 min for setup of locations and targets, 1 min for image acquisition, 3.5 min for point cloud generation, 0.25 min for deriving measurements from Computree, and 0.25 min for curve matching). The resulting point clouds contained between 10,000 and 62,000 stem points for the shortest (tree 14) and tallest (tree 1) trees, respectively. Points covered an area immediately around the stem and ranged from ground level to a maximum height of 4 to 8 m. Diameter estimates were derived from DTP point clouds, shown in Figs. 4 and 5, and showed good, unbiased correspondence with manual measurements (1.28 cm RMSE, 5.15 RMSE% for DBH). This degree of correspondence was observed at other heights along the stem, although the DTP-derived point clouds rarely allowed extraction of diameter measurements above 5 m. For example, for one of the tallest measured trees (tree 9, a large lodgepole pine, 26.4 cm DBH, 26.2 m height; Table 1), the generated point cloud allowed stem reconstruction up to 9 m above the ground. The shortest stem reconstruction was to a height of 3.5 m on a small black spruce (tree 14, 10.8 cm DBH, 9.7 m height). In general, lower stem heights (e.g., < 3 m) had approximately 50% of the circumference of the tree represented by points. However, point clouds further up the stem generally became more obscured by canopy or branches, resulting in less direct observation by some of the camera perspectives. This resulted in approximately 25% of the stem circumference being represented by points. Despite this, points coming from one camera location were often enough to derive a sufficiently accurate diameter estimate. For example, a separate analysis using only six images at a single location for point cloud reconstruction yielded an RMSE of 2.00 cm, or 8.10%, for DBH estimation.

Fig. 4
figure 4

Comparison between observed and predicted DBH measurements (n = 15) from DTP point clouds

Fig. 5
figure 5

Error (RMSE) of diameter estimates at various heights up the stem, with points indicating RMSE at each measured height and lines showing trend generated by smoothed conditional means. Different colored points and lines refer to the five different estimation approaches. The graph on the left shows the error in terms of the absolute height, while the one on the right shows the error at heights relative to the total tree height in 10% increments. Values over 100% on the y-axis in the right graph indicate an incorrect estimate of total tree height

3.2 Upper stem diameters and taper assessment

The relationship between measured and predicted stem diameters for all trees and all evaluation methods is shown in Fig. 5. For comparison across trees of different heights, the relationship is shown as both the measurement error by absolute height up the stem and the percentage of total tree height for individual trees. For diameters at lower sections of the stem (< 8 m), both curve matching techniques (methods 1 and 2) performed better than other approaches (~ 0.5 cm RMSE). Above 10 m, or approximately 30% of total tree height (across stems), method 3 (using the field-measured DBH and height) yielded the most accurate diameter estimates (< 1 cm RMSE). In some cases, either the curve matching or an allometric equation yielded inaccurate estimates of total tree height, producing larger discrepancies at upper parts of stems (> 75% of total tree height). Despite this, all methods were generally successful at characterizing stem diameter, with the most accurate measurements coming at points in the bottom 13 m or 50% of the stem (< 1.5 cm RMSE).

3.3 Volume predictions

The relationship between measured and predicted volumes for the different approaches is shown in Fig. 6. Method 3 (using the field-measured DBH and height) produced the most accurate overall predictions of volume (0.094 m3 RMSE, 15.5 RMSE%). Independent-samples t tests were conducted to compare mean diameter accuracy for all techniques. For all techniques, there was no significant difference in the mean accuracy of diameters, suggesting that no one technique performed better or worse than the others. All methods of volume calculation were slightly negatively biased compared with the ground-measured reference. Of the point cloud–based approaches, method 5 (the predicted DBH and allometrically assigned height) yielded the most accurate estimates (0.099 m3 RMSE, 16.3 RMSE%), but this was only marginally better and not statistically different than the curve matching approaches (methods 1 and 2). Method 1 (unweighted curve matching) produced more accurate volume estimates (0.110 m3 RMSE, 18.1 RMSE%) than method 2 (weighted curve matching; 0.120 m3 RMSE, 19.6 RMSE%), indicating that this may be the better of the two curve matching methods. However, for all 15 trees, no method performed significantly better or worse than the others.

Fig. 6
figure 6

The accuracies of the five different methods for volume estimation. Each point represents a tree and is colored according to the method presented

4 Discussion

4.1 Diameter estimates

Stem diameters at multiple heights were extracted from DTP point clouds. Point clouds produced accurate estimates of DBH (Fig. 4), showing a RMSE of 1.28 cm and a RMSE% of 5.15. Point cloud–based curve matching approaches (methods 1 and 2) produced the most accurate measurements for the lowest 8 m, or approximately 30% of the stem, while using a known DBH and height, method 3 was the most accurate method for the highest parts of the stem. For trees with irregular allometries, it was possible that the taper models used in this study inaccurately characterized the diameter at different parts of the stem. However, the models were determined to be generally accurate in their characterization of stem taper. In some cases, methods 1, 2, 4, and 5 produced inaccurate estimates of total tree height resulting in inaccurate predictions for volumes and upper stem diameters. This resulted in higher RMSE values for upper parts of trees (> 50% of total height) in scenarios where the total height was unknown. However, most tree height estimates were within 3 m of the true height (after falling), and all methods produced RMSEs of less than approximately 1.5 cm for the lowest 50% of the stem. Most tree heights measured in the field by a laser hypsometer were within 1 m of the true height (after falling), but deviated by about 3 m for the tallest two trees. These discrepancies between field-measured (hypsometer) and observed tree height are consistent with findings from Luoma et al. (2017), who determined that the standard deviation of field-measured height was 0.5 m (2.9%), up to a maximum of 4.2 m.

When compared with other studies estimating individual tree DBH from DTP point clouds (Miller et al. 2015; Bauwens et al. 2017; Fang and Strimbu 2017), we used fewer photographs (12) and evaluated more species (6) while achieving similar levels of accuracy (Table 2). For example, Fang and Strimbu (2017) reported DBH estimates with an RMSE% of 5 in a monospecific plantation. Our similar accuracy with six species indicates that there may not be a strong correlation between tree species and DBH accuracy; however, small sample sizes in both studies may not be sufficiently large to make this determination.

Our achieved accuracy with a relatively low number of images may have resulted from the inclusion of six scale bars (i.e., one between each set of adjacent images at three different heights), which was set to the distance between the cameras as they were mounted on the pole. This allowed the point clouds to be scaled more accurately than using the target locations alone. Forsman et al. (2016) also used a camera rig (multiple cameras mounted to a portable device) to scale the images with known distances while using an average of less than three images per tree to detect and measure stems on sample plots. Therefore, a rig-based system with known distances between cameras may be instrumental in producing accurate estimates of tree dimensions in cases where relatively few images per tree are captured.

4.2 Volume assessment

Photogrammetric point clouds were generally accurate in estimating tree volume. For smaller trees (i.e., under 0.5 m3), all methods of volume estimation produced similar accuracies. For the largest two trees, method 3 (using the field-measured DBH and height) produced inaccurate estimates, possibly due to inaccurate height measurements as taken from the ground, which has been shown to be influenced by stand conditions, crown class, and tree species (Wang et al. 2019). Although more accurate diameter measurements at upper parts of the stem came from method 3 (field-measured DBH and height), the majority of a tree’s volume is in the lowest portions of the stem, indicating that accurate diameter measurements at the bottom—possibly coming from DTP point clouds—could also result in more accurate volume estimates. For example, the lower 50% of tree stems in our study contained about 80% of the total tree’s volume.

Results from this study described the ability of 360° cameras to derive detailed tree-level measurements. The cameras’ large field of view means that resulting point clouds have a larger coverage area than traditional frame cameras, which are typically employed during field operations. This outlines the potential of the 360° cameras to characterize larger areas such as sample plots or stands, as larger coverage from individual images would result in fewer images being required for point cloud generation. As a comparison, the frame cameras used in Liang et al. (2014a) and Mokroš et al. (2018) were able to successfully characterize DBH for trees on sample plots, but used up to 973 and 1271 images for 900 and 1225-m2 plots, respectively. Spherical images could also provide the basis for a combination of terrestrial and aerial photogrammetric point clouds, such as in Mikita et al. (2016), who used terrestrial and aerial images to characterize tree DBHs and volumes in a 0.8-ha stand. For individual trees, however, only a subset of the entire 360° field of view was used for tree reconstruction, indicating that similar accuracies may be achieved using a rig-based system of wide-angle or fisheye lens cameras.

Of the related studies listed in Table 2, only one estimated volume from DTP point clouds (Mikita et al. 2016). Volume estimates from point clouds in their study (0.082 m3) were slightly more accurate than those reported here (0.099 m3). In this study, we evaluated the ability to derive individual tree characteristics based on point clouds created from fewer images and tested on more species than presented in Mikita et al. (2016). Additionally, Mikita et al. (2016) combined DTP and DAP point clouds, while our study was limited to ground-based images. More study is needed to determine the relationships, if any, between point cloud accuracy and characteristics such as tree size, branchiness, species, or stem density of the surrounding area. Based on an international benchmarking study of TLS by Liang et al. (2018), stem detection rates decreased with decreasing mean DBHs and stem density, while DBH estimates were stable across stand conditions. Smaller trees will have less surface area for point cloud reconstruction, meaning that there may be a resulting increase in error of DBH estimation with decreasing tree size (Ryding et al. 2015).

4.3 Applications

When compared with DTP, TLS provides more comprehensive point clouds that can be used for more detailed study such as wood quality and has the ability to return points from occluded areas such as in stands with a high stem density or on stems with many branches. However, the results seen in this study indicate the potential for DTP to provide similar levels of accuracy to TLS in DBH and volume estimates. Studies using a single TLS scan reported 1–3 cm RMSE for DBH (Maas et al. 2008; Liang and Hyyppä 2013) and ~ 10 RMSE% for volume (Liang et al. 2014b), similar to the results reported here. Liang et al. (2018) also reported accuracies of 10 and 20 RMSE% for “easy” (low stem density and high mean DBH) and “difficult” (high stem density and low mean DBH) plots, respectively. The accuracies attained in this study indicate the potential for DTP for supporting forest inventories. National Forest Inventories (NFIs) have typical accuracy requirements of 0–2 cm for DBH, 10–20% for volume, and 1–3 cm for upper stem diameters (Liang et al. 2016). Each of these accuracy requirements was met using DTP point clouds techniques in this study. Other potential applications of raw images and resulting point clouds include the estimation of canopy leaf area (Bréda 2003) or as inputs to centroid sampling of tree volume (Wiant et al. 1992).

Additionally, point clouds provide measurements which can be stored as an objective, 3D snapshot of a tree or forest condition at a given point in time. This indicates potential for the use of DTP in calibrating or validating models of forest growth. If images of the same tree are acquired at multiple times, a time series of point clouds could be generated and analyzed to monitor tree growth or change either at an individual tree or at stand level (Sheppard et al. 2016; Liang et al. 2012). Should allometric or taper models not exist, be inaccurate, or require update, point clouds could provide reference data to validation or parameterization of such models.

While TLS provides attributes such as branching structure and direct measurements of upper stem diameters, the cost of handheld cameras is in the hundreds of dollars, far less than TLS units, which can be two orders of magnitude higher. Although the focus of this study was on the structural qualities of the point clouds, their spectral attributes could be used, similar to aerial images, to assess species compositions (e.g., Packalen and Maltamo 2006) and tree condition or quality (e.g., Goodbody et al. 2018). Overall, the low cost and portability of the cameras, in addition to the objectivity and storability of the point clouds, show their value as a tool in forest inventory and modeling.

5 Conclusion

Changing resource demands and climatic conditions require quick and inexpensive means of deriving robust and accurate forest inventory measurements. DTP is one such tool that could be used to enhance traditional forest inventories. This study showed that DTP is able to produce sufficiently accurate estimates of volume and diameter that generally only slightly differed from traditional inventory techniques. The inclusion of fixed scale bars between known camera locations was critical in deriving accurate measurements from resulting point clouds. Using this methodology, accurate point clouds were produced for a variety of tree species and sizes, demonstrating the possibility of using such a technology at larger scales as well as the utility of DTP in characterizing trees present in complex and changing habitats such as the boreal mixedwood forests.


Download references


We thank the staff at Alberta Agriculture and Forestry for their input and assistance with the fieldwork.


This research was funded by the AWARE (Assessment of Wood Attributes using Remote Sensing) Natural Sciences and Engineering Research Council of Canada Collaborative Research and Development grant (NSERC File: CRDPJ 462973–14) to a team led by Nicholas Coops with support from West Fraser Timber Co. Ltd. Additional funding for field data collection was provided by Alberta Agriculture and Forestry and West Fraser Timber Co. Ltd.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Christopher Mulverhill.

Ethics declarations

Conflict of interest

The authors declare that they have no conflict of interest.

Additional information

Handling Editor: David Drew

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Contribution of the co-authors

Conceptualization: all; methodology: CM, NCC, PT, CWB; formal analysis: CM, NCC, and PT; resources: all; writing: all; funding acquisition: NCC and CWB.

This article is part of the Topical Collection on Frontiers in Modelling Future Forest Growth Yield and Wood Properties

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mulverhill, C., Coops, N.C., Tompalski, P. et al. The utility of terrestrial photogrammetry for assessment of tree volume and taper in boreal mixedwood forests. Annals of Forest Science 76, 83 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: