Accuracy of common stem volume formulae using terrestrial photogrammetric point clouds: a case study with savanna trees in Benin

Hospice A. Akpo · Gilbert Atindogbé ·Maxwell C. Obiakara · Arios B. Adjinanoukon ·Madaï Gbedolo · Noël H. Fonton

Abstract Recent applications of digital photogrammetry in forestry have highlighted its utility as a viable mensuration technique. However, in tropical regions little research has been done on the accuracy of this approach for stem volume calculation. In this study, the performance of Structure from Motion photogrammetry for estimating individual tree stem volume in relation to traditional approaches was evaluated. We selected 30 trees from five savanna species growing at the periphery of the W National Park in northern Benin and measured their circumferences at different heights using traditional tape and clinometer. Stem volumes of sample trees were estimated from the measured circumferences using nine volumetric formulae for solids of revolution, including cylinder, cone, paraboloid, neiloid and their respective fustrums. Each tree was photographed and stem volume determined using a taper function derived from tri-dimensional stem models. This reference volume was compared with the results of formulaic estimations.Tree stem prof iles were further decomposed into dif ferent portions, approximately corresponding to the stump, butt logs and logs, and the suitability of each solid of revolution was assessed for simulating the resulting shapes. Stem volumes calculated using the fustrums of paraboloid and neiloid formulae were the closest to reference volumes with a bias and root mean square error of 8.0% and 24.4%, respectively. Stems closely resembled fustrums of a paraboloid and a neiloid. Individual stem portions assumed dif ferent solids as follows: fustrums of paraboloid and neiloid were more prevalent from the stump to breast height, while a paraboloid closely matched stem shapes beyond this point. Therefore, a more accurate stem volumetric estimate was attained when stems were considered as a composite of at least three geometric solids.

Keywords Structure from motion photogrammetry ·Point cloud data · Stem volume · Savanna species · Benin

Introduction

Accurate estimates of the volume of tree stems is essential for characterizing forest stands and managing forest resources (West 2015). Stem volume is also the basis for commercial timber transactions because it is strongly correlated with the carbon storage potential of individual trees(Kankare et al. 2013; Yu et al. 2013). However, obtaining accurate stem volumes of standing trees is particularly diff icult and time-consuming, given that stems are not perfect geometric shapes. Hence, this metric can be easily overestimated, especially for trees with irregularly shaped stems(Dean 2003; Bauwens et al. 2017). Several factors have been attributed to inaccurate stem volume estimation, mainly measurement errors and the wrong choice of allometric equations (Dean 2003; West 2015).

Of the several approaches to determine stem volumes, the most widely used are based on formulae of geometric solids such as cylinder, paraboloid, neiloid and their fustrums.Other less common approaches include water displacement,graphical and integration methods (Kershaw et al. 2016).The most common formulae for measuring stem volume are based on the product of the length of a stem segment and the weighted average of its cross-sectional area (West 2015). These formulae vary depending on the number of cross-sections considered, their relative position along the stem segment and the weight assigned to each cross-section.For example, with Huberʼs formula, a cross-sectional area is taken only at the middle of the segment. Smalian’s formula uses two cross-sectional areas, at the bottom and top of the stem, and assigns equal weight to both measurements. Newton’s formula, on the other hand, dif fers from Smalian’s by an additional cross-sectional area measured at the middle of the stem segment with a fourfold weight assigned to it.

The choice of one of these common formulae depends mostly on practical considerations, such as the ability to measure more circumferences or diameters on a given stem segment. Because it involves three measurements on a stem section, Newtonʼs formula is generally more accurate (Kershaw et al. 2016). Moreover, it takes into account the wide range of shapes (paraboloid, neiloid, cylinder and cone) commonly assumed by tree stems. However, Smalian’s approach is often preferred because mid-sectional measurements are not always practical in the f ield (West 2015).

Exact stem volume determination is conventionally done using the water displacement approach, also known as xylometry, which involves the immersion of logs in a known amount of water and measuring the resulting increase in volume. Although this technique is quite accurate, it is very restrictive in its application (Özçelik et al. 2008;Akossou et al. 2013). For example, considerable logistical constraints must be overcome to make useful quantities of water available in typical f ield settings. Moreover, xylometry is destructive and may not always justify overestimations often associated with the abovementioned methods and more recent alternatives such as the centroid and centre of gravity approaches (Özçelik et al. 2008). However, toward overcoming the shortfalls of xylometry, many studies have explored the utility of point cloud data derived from terrestrial laser scanning, TLS (Kankare et al. 2013; Yu et al.2013; Liang et al. 2016; Saarinen et al. 2017; Hyyppä et al.2020) and Structure from Motion (SfM) photogrammetry(Morgenroth and Gómez 2014; Miller et al. 2015; Koeser et al. 2016; Mulverhill et al. 2019) for accurate stem volume determination.

In many cases, the accuracy of terrestrial image-based point cloud data has been shown to be marginally lower than that of TLS point clouds (Liang et al. 2014, 2015;Panagiotidis et al. 2016 ; Huang et al. 2018; Piermattei et al. 2019). In an extensive comparison involving several image acquisition scenarios for plot-based mensuration,Liang et al. ( 2015) found a small mean dif ference (using tape measurements as a reference) between laser and imagederived breast height diameter of Scots pines, in the order of 0.23–0.28 cm and − 1.34 to 3.07 cm, respectively. Another comparison between these methods (Piermattei et al. 2019)across heterogeneous forest plots produced similar breast height diameter errors ranging from 0.87 to 3.75 cm (TLS)and 1.21 to 5.07 cm (SfM photogrammetry). The recent use of highly mobile TLS devices has eliminated the logistical constraints hitherto associated with this technique in remote and hard-to-access locations (Hyyppä et al. 2020). However,costs still remain a major limitation as opposed to imagebased approaches (Huang et al. 2018; Iglhaut et al. 2019).

With recent developments in the f ields of computer vision and photogrammetry, smart phone and digital single-lens ref lex cameras have been increasingly used as low-cost devices for creating tri-dimensional models. Promising results have been reported from the application of SfM photogrammetry in forestry at scales ranging from forest plots (Liu et al. 2018; Mokroš et al. 2018; Piermattei et al.2019) to individual trees (Morgenroth and Gómez 2014;Miller et al. 2015; Bauwens et al. 2017; Huang et al. 2018;Mulverhill et al. 2019; Roberts et al. 2019; Marzulli et al.2020). However, many studies exploiting photogrammetric point cloud data seem to focus more on linear metrics, (i.e.,breast height diameter and tree height), than volumetric metrics. The f irst direct comparison between actual and point cloud-derived data of tree volume was carried out by Miller et al. ( 2015). They showed that photogrammetrically derived stem volumes of potted trees were reasonably close to their true volume—determined xylometrically—with a root mean square error and mean dif ference (bias) of 12.3% and − 8.2%,respectively. Mulverhill et al. ( 2019) also found no signif icant dif ference between stem volumes obtained from point cloud data and volumetric estimates derived from sectional measured stem diameters of felled trees.

In the context of precision forestry, the present study examines the application of ground-based SfM photogrammetry for tri-dimensional tree structure. The strength of this study was the ability to capture the full height of the mostly short-stemmed savanna species from the ground, as opposed to some recent studies that highlighted the limitations of this method, especially for excurrent trees, ones with a single leader (Mulverhill et al. 2019; Marzulli et al. 2020). Thus,our specif ic objectives were to compare stem volume estimations based on common geometric solids using terrestrial photogrammetric point clouds as a reference and to examine the suitability of each solid for simulating dif ferent stem shapes in relation to height.

Materials and methods

Study site

The study site, in the buf fer zone of the W National Park(11°26′ to 12°25′ N; 2°48′ to 3°05′ E; 164 m a.s.l.) is in northern Benin (Fig. 1). It has a strong seasonal tropical climate and is located within the West African Sudanian Savanna which stretches from Senegal to Nigeria. The mean temperature in this area is 40 °C, and rainfall varies between 600 and 700 mm a −1 .

Sampling and stem volume estimation using traditional techniques

Our sample was based on f ive of the most dominant tree species in the study area selected from the National Forest Inventory database (Projet 2009):Anogeissus leiocarpa(DC.) Guill. & Perr.,Bombax costatumPellegr. & Vuillet,Sclerocarya birrea(A. Rich.) Hochst.,Terminalia laxif loraEngl. andVitellaria paradoxaGaertn. f.. For each species, six trees were sampled in each of the following DBH classes (cm): 10 − 20, > 20 − 30, > 30 − 40, > 40 − 50, > 5 0 − 60 and > 60. Measurements were made using traditional tools (diameter tape and clinometers) and are presented in Table S1.Stem circumferences were measured at 0.10 m, 0.50 m,1.30 m and 1.50 m above ground level. Additional measurements were taken at 1-m intervals from 1.50 m up to bole height (i.e., at the top of the stem below the f irst branch)using a ladder in few cases. Prior to measuring, all heights were marked by a thin, horizontal red line drawn around the trunk. Stem segments def ined by two consecutive markings formed the basis for volumetric estimations. Volume of stem segments was estimated using the formulae in Table 1(Kershaw et al. 2016), and stem volume was obtained by summing individual segment volumes (West 2015). Values computed from each species are in Table S2.

Image acquisition and photogrammetric point cloud generation

Fig. 1 Map of study area

Tri-dimensional models of the sample trees were built using SfM (Morgenroth and Gómez 2014; Miller et al. 2015).Image acquisition was done on a circular path with a 2-m radius around the stem (Akpo et al. 2020). The distance between two consecutive camera positions was def ined by an angle of 15° on the circular path (Fig. 2). These positions were marked by stakes driven into the ground. A series of highly overlapping images were captured at a given position from the base of the stem to the top before moving to the next one. All images were taken with a Canon 77 D camera equipped with a Canon EF 50 mm f/1.8 STM lens. For optimal results, the camera was set to its maximum resolution(6000 × 4000 pixels) and held parallel to the stem in portrait orientation. The aperture priority mode (f/3.2) and lowest ISO value (100) were maintained during image acquisition to obtain clear and high quality images. Before photographing each tree, a 1.5 m ruler was vertically pinned to the stem for subsequent scaling.

Although stem circumference can be estimated on manually delineated point-based stem contours in a geographicalinformation system (GIS) environment (Bauwens et al. 2017;Akpo et al. 2020), in this study, we automated stem crosssection delineation for volume estimation in the R programming language (R Core Team 2017). Scaled tri-dimensional tree models were built using Agisoft PhotoScan Pro. 0.9.1(Agisoft LCC, St. Petersburg, Russia). The workf low in this program consists of three main steps: (1) loading images,(2) aligning images and (3) building a dense point cloud.PhotoScan users can control outputs at each step via a suite of settings depending on the desired accuracy and computer processing power (Fang and Strimbu 2017; Mulverhill et al.2019). The steps and settings used in generating tree models in this study are shown in Fig. 3, in line with Akpo et al.( 2020). Further processing was based on the spatial coordinates of points in the resulting tri-dimensional models.

Table 1 Common formulae for stem segment volume estimation

Fig. 2 Image acquisition schema

Fig. 3 Image processing workf low of PhotoScan settings used in this study

Stem volume determination using point cloud data

Spatial coordinates of point clouds were imported in the R programming environment, and horizontal sections (hereafter referred to as slices) of varying thicknesses (0.1 cm,1 cm, 50 cm, and 100 cm) were created from the base up to bole height. In most cases, slices had points well outside stem contours. To reduce errors in cross-section delineation,these outliers were detected and removed using thelofactor()function of theDMwRpackage (Torgo and Torgo 2013).Slice cross-sectional area was computed based on the Cartesian coordinates of each point of the slice. The centroid of each slice was determined and taken as the origin of a new Cartesian coordinate system in which coordinates of all points of a slice were derived. Then, distances and azimuths of all points were determined in relation to the centroid. The azimuth of a point (X1,Y1) was given by:

Data analysis

Stem volumes estimated using the nine formulae were reported in terms of arithmetic mean, standard error of the mean, minimum and maximum. The accuracy of the volume was evaluated using the relative bias (B) and root mean square error (RMSE) as def ined in the equations below:

whereVc(i)is stem volume estimated in the f ield using a tape for treei,VRef(i)is volume obtained using the reference method, andnthe total number of trees.

Results

Comparison of volume formulae

Stem volumes estimated from f ield measurements were highest forV. paradoxaand lowest forA. leiocarpa.Regardless of species, the neiloid formula produced the lowest volumes whereas cylinder estimates were the highest (Table S2).Mean bias and RMSE values between f ield and photogrammetric estimates of stem volume based on increasing stem slice thicknesses are presented in Table 2. The negative bias obtained from the cone, neiloid, paraboloid and paracone formulae suggests that these solids underestimate stem volume. Conversely, results based on the cylinder, paraboloid fustrum (Huber’s and Smalian’s formulae), and neiloid/cone fustrums (Newton’s formula), and neiloid (Eq. 5) indicated that these geometrical shapes tend to overestimate stem volume. Stem segment volume from fustrums of neiloid and Huber’s formulae produced the lowest bias (8.0%) andRMSE(20.4%) in relation to stem volumes computed by integrating surface area for 1 mm and 1 cm slices. However,a twofold increase in bias was obtained with 50 cm slices(Table 2).

Variation of segment stem volume in relation to height

The variation in bias andRMSEin relation to stem height is shown in Fig. 4. The change in bias was similar across the length of the stem for volumes calculated using neiloid fustrum, Huber’s, Smalian’s and Newton’s formulae. This pattern suggests that stems can be partitioned into the following portions: (1) butt logs (0.3 m ≤ height ≤ 1 m) for which bias values ranged from 17.0 to 18.5%, (2) middle logs (1 m < height ≤ 2.5 m) with bias values from 22.7 to 28.6%, and (3) top logs 2.5 m < height ≤ 4 m with a bias value of 19.4%. In contrast, the lowestRMSEvalues were obtained only for butt logs using both fustrums of paraboloid and neiloid formulae. Above breast height (1.3 m), the stems of the study species assumed a paraboloid form. These results are summarized in Table S3.

Fig. 4 Accuracy of stem volume computed from photogrammetric point clouds data using dif ferent geometric solids in relation to stem height

Discussion

A number of factors, including high tree density, occlusions and poor light conditions can limit the photogrammetric data acquisition process in a typical forest. These factors, combined with the low resolution of ground-operated passive sensors explain the inability to ef ficiently reconstruct the tridimensional structure of trees or forest plots for subsequent volume estimations. In this study, we assessed the accuracy of formulae for nine geometric solids for estimating stem volumes of f ive savanna species using terrestrial photogrammetry as a reference method. The volume of inherently irregular solids such as tree stems is most accurately measured by xylometry (Özçelik et al. 2008; Akossou et al. 2013).Because this approach was not practical for this study, we took advantage of terrestrial photogrammetric point cloud data to derive the stem volume of standing trees. The use of this technique as a reference can be justif ied by its accuracy, close to but a little lower than that of the more popular and costly laser scanning technologies (Panagiotidis et al.2016; Huang et al. 2018). In addition, Miller et al. ( 2015)tested the suitability of this method for stem volume estimation, whereas Koeser et al. ( 2016) underscored its utility for estimating the volume of geometrically complex shapes assumed by tree root systems.

Our accuracy assessment showed that four out of the nine geometric shapes tested were useful for determining stem volume. The most accurate of the four were the fustrums of paraboloid (Huber’s formula) and neiloid with a similar performance. Calculation of stem volume using a fustrum of neiloid requires the cross-sectional areas at the base and upper end of a stem segments whereas Huberʼs formula uses only the middle cross-sectional area. Thus, based on the principle of parsimony, Huber’s formula can be considered as the most accurate approach for estimating stem volume of the selected savanna trees. This f inding is in line with the results of previous studies (Filho et al. 2000; Akossou et al.2013). Using xylometry, these authors showed that Huber’s formula produced the most accurate log volume forTectona grandisL.f. andPinus elliottiiEngelm., respectively. Therefore, as expected for most deliquescent species, the general shape of the stem of the f ive most abundant ones in the Sudanian Savanna closely resemble fustrums of paraboloid and neiloid rather than their non-truncated equivalents (Kershaw et al. 2016). Our results also indicated that Newton’s formula is less prone to bias, thereby conf irming its versatility for most log shapes (West 2015).

The important increase in bias and RMSE obtained between stem volumes derived from integral calculus of increasing slice/sections and their formulaically calculated equivalents provide more insight into the ef fect of log length on stem volume accuracy. This agrees with earlier reports by Filho et al. ( 2000) and Akossou et al. ( 2013 ) in which increasing log lengths ofP. elliottiiandT. grandis, respectively, were compared with volumes determined by water displacement. Miller et al. ( 2015) reported a bias and RMSE of − 8.2% and 12.3%, respectively, for volumetric estimates of the stems of small potted plants. Stem volumes based on tape measurements using Huber’s formula showed a bias comparable to that of Miller et al. ( 2015). Contrary to their result, most formulaic values in this study had a positive bias, suggesting an overestimation of volume. To the best of our knowledge, the study by Miller et al. ( 2015) is the only one to directly compare photogrammetric volume to true stem volume determined by water displacement. However,the physical context of their study and their use of small potted plants in a controlled environment is dif ferent from that of the present study, which is based on trees in their natural habitat.

Mulverhill et al. ( 2019) found no signif icant dif ferences between point cloud-derived stem volumes of standing trees and volumes estimated on felled trees. They reported an accuracy (RMSE = 18.1%) comparable with that of this study. It is noteworthy that some strengths and weaknesses of the tested method were evident in recent studies in natural forests (Mulverhill et al. 2019; Marzulli et al. 2020). Because of the inability of ground-based photographs to capture the full height of stems, these authors resorted to a database of taper curves to generate metrics from the non-captured upper stem portions (above 5 m). Thus, more dif ficulties may be expected in volumetric estimation of stems of excurrent trees using photogrammetric point clouds as opposed to the mostly deliquescent savanna trees in this study. Moreso, the maximum stem height (4.5 m) was measured onB. costatum.

Our results show that stem portions assume dif ferent geometrical shapes in relation to height. Tree stems in this study were morphologically complex, and perhaps better described as a “stack” of geometrical solids, with some solids more prevalent at a given height than others (Kershaw et al. 2016).The stem could be easily partitioned into three major portions using breast height as the reference level. Stem portions below breast height assumed the shape of fustrums of paraboloid and neiloid, whereas those above closely resembled a paraboloid. These shapes are in consonance with the theory of tree stem forms (Burkhart and Tomé 2012; Kershaw et al. 2016).

Conclusions

This study identif ied the best geometrical solid for determining stem volume. Terrestrial photogrammetry is welladapted for savanna trees because of their short heights. The results show that Huber’s formula was the best approach for estimating stem volumes of tropical savanna species in West Africa. We also showed that stem section length has an ef fect on the tested volume formulae because greater log segments led to greater errors. The stem of the study species was a composite of three classical geometrical shapes.Our results should provide foresters with more appropriate options for estimating the volume of certain natural species.

Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source,provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.