Close Home
iForest - Biogeosciences and Forestry
vol. 11, pp. 181-188
Copyright © 2018 by the Italian Society of Silviculture and Forest Ecology
doi: 10.3832/ifor2520-010

Research Articles

Determining basic forest stand characteristics using airborne laser scanning in mixed forest stands of Central Europe

Róbert Smreček (1)Corresponding author, Zuzana Michnová (2), Ivan Sačkov (3), Zuzana Danihelová (4), Martina Levická (1), Ján Tuček (1)

Introduction 

Quantitative information on stand and tree characteristics, such as the number of trees, height, diameter and volume, is basic to make decisions in forest management ([35], [34]). Data on these characteristics and the forest environment can be obtained using traditional (ground) methods or by remote sensing (RS) techniques ([25], [36]). Currently, ground methods are being more and more replaced by automated data collection based on aerial and satellite images or data from airborne laser scanning (ALS - [15]), thanks to their cost and time efficiency. Photogrammetry methods have been applied to forest mapping for long time, and ground geodetic methods were used only in cases where details were not achievable from orthophotos ([36]). The main advantage of ALS lies in the 3D point cloud obtainable by laser impulses passing through tree crowns ([26]). The development of ALS started during the 1970s and 1980s, and suitable scanning methods were derived during the 1990s. The first applications were topographically oriented. Since then, the development of this technology was fast and ALS are currently used in a wide range of applications, including forestry ([2], [13]). The history of the development and use of ALS in forestry are described in Hyyppä et al. ([6]).

The number of trees and their position on the ground is the most important information achievable from ALS data ([27]). The best results are obtained for dominant trees in the canopy layer: in younger stands and trees in the understorey layer, the number of trees is likely to be underestimated ([21], [5]). A multitude of algorithms and procedures for the identification of individual trees have been proposed in the literature ([8], [10], [11], [28], [32]), which can be grouped in two basic approaches: (i) an area-based approach where forest characteristics are estimated using statistical analyses and models between ALS data and terrestrial measurements; and (ii) an individual tree-based approach whereby individual trees are identified from ALS data, visually or by segmentation processes, and dendrometric parameters are extracted for individual trees ([2], [10], [32]). Accuracy of tree detection is usually higher in conifer stands as compared to stands with deciduous trees ([27]). The correct identification of trees and their number is critical to derive further tree and stand characteristics such as height, diameter at breast height (DBH), volume, biomass, etc. ([28], [31]).

Tree height can be directly obtained from ALS data, though it is often underestimated as tree tops do not always reflect the laser impulses, particularly when a low scanning density is used ([17]). On the other hand, tree height is likely to be overestimated in terrain with great slopes and mountain areas ([2]). In order to accurately assess tree height a suitable density of the point cloud is crucial. Takahashi et al. ([30]) recommended a density of at least 8.8 points per m2 in order to achieve a deviation in height estimates lower than 1 m; Andersen et al. ([1]) reported a minimum number of 4-5 points per m2 to correctly identify the tree tops. With old trees having large crowns, a density of two points per m2 is assumed to be sufficient. Nevertheless, if quality results are to be achieved, 10 or more points per m2 are appropriate for data extraction ([10]).

While tree height and crown diameter can be measured directly from ALS data, DBH and volume have to be derived following known relations ([2], [6]). The DBH can be derived by regression analysis using data measured on the ground ([12], [34]) or formulas derived from existing relationships ([15], [16]). One of the methods commonly used to derive tree DBH from ALS data is its calculation using the derived height and crown size as predictors. However, crown sizes automatically obtained from ALS data are often inaccurate, and therefore their use may lead to a fairly significant degree of uncertainty ([10]). Moreover, the drawback of a large number of derived relationships between tree height and diameter is their strict dependence on specific stand structures and local natural conditions.

In this paper an automatic approach for tree identification and estimation of their dendrometric characteristics is presented, aimed at assessing the influence of spatial resolution of ALS data on both tree identification and the determination of tree parameters. Tree identification was done at four different spatial resolutions. The height of trees was determined automatically from the normalised digital surface model (nDSM) obtained from ALS data. Determination of the tree diameter was obtained by regression analysis based on reference (ground) data. ALS-derived tree parameters were then compared to those calculated from reference data obtained from field measurements. The aim was to find dependencies between ALS-derived tree parameters and the spatial resolution of ALS data.

Methodology 

Study area

The study area is a part of the University Forest Enterprise at the Technical University in Zvolen, Slovak Republic (48° 37′ N, 19° 04′ - Fig. 1). The whole territory is a part of an ancient eruptive region (Kremnické vrchy) characterized by a broken relief with different climatic characteristics.

Fig. 1 - Geographic location of the four permanent research plots in the University Forest Enterprise of the Technical University in Zvolen (Slovak Republic).

The study area contained four permament research plots (RP). RP 1 is located in the northern part of stand 350a. The stand is 105 years old, with an area of 0.5 ha and average slope of 23.3%. The most abundant tree species in RP1 is the European beach (Fagus sylvatica L.), representing almost 40% of standing trees. The share of broadleaved and coniferous trees was 92% and 8%, respectively. RP2 (area = 0.3 ha; average slope = 34.6%; stand age = 100 years) is located in the northern part of stand 314, where European silver fir (Abies alba Mill.) is the most abundant tree species (almost 36% of extant trees). The share between broadleaved and coniferous trees was 57% and 43%, respectively. RP3 (northern part of the stand 304b) has stand age of 90 years, an area of 0.25 ha and average slope of 7.1%. European silver fir and Norway spruce (Picea abies Karst.) were the most abundant species (31% and 30%, respectively), with a share of broadleaved trees of 39%. RP4 is located in the central part of stand 518a (150 years old; area = 0.25 ha; average slope = 23.7%). European beech represents 100% of trees here.

Terrain slope was calculated from the digital terrain model (DTM - see below), while the tree composition was obtained by field measurement. Stand age was obtained from the Forest management plan (see [24] and [3] for more details).

Data

ALS data were provided by a vendor in September 2011. The airborne laser scanner employed was a Riegl L-680i® (Riegl Laser Measurement Systems Gmbh, Horn, Austria), with a flight altitude of 700 m and a 50° field of view, PRR 320 kHz and SR 122 Hz. The resulting RMSE of the absolute data position was 0.047 m. The ALS data provided by the vendor were recorded as ground points and non-ground points (representing vegetation, buildings and other objects). ALS data were processed using the software Microstation V8 with the application TerraScan® by the vendor. The main characteristics of the ALS dataset are reported in Tab. 1. The images used for the manual vectorisation of data were also taken during the flight. The camera employed was a Vexcel UltraCamX® (Vexcel Imaging GmbH, Graz, Austria). GrafNav® (Waypoint, Novatel Inc., Calgary, Canada) and AEROoffice® (IGI mbh, Kreuztal, Germany) software were used for data processing. The image orthorectification was performed at the Department of Forest Management and Geodesy (Technical University in Zvolen, Slovak Republic) using the Inpho® software package (Trimble Geospatial Inc., Sunnyvale, CA, USA). The spatial resolution of the orthophotos was 10 cm.

Tab. 1 - Characteristics of the ALS dataset used for tree identification. (RP): research plot; (DSM): digital surface model; (DTM): digital terrain model.

Field measurement was conducted in each RP. The position of the RP centre was measured using the Global Navigation Satellite System (GNSS) and further corrected by visual analysis based on the normalized digital surface models. After manual correction, a sub-meter horizontal accuracy was expected for all RP. The position, height and diameter at breast height (DBH) of each tree were measured using FieldMap® (IFER Ltd., Jilove u Prahy, Czech Republic). The position in the storey of each tree was also identified. Depending on the covering exerted by tree crowns in the upper canopy layer, a metric was set for the potential visibility of each tree from above (i.e., its visibility by the ortophotos). Trees were classified in two categories: (0) no visibility (tree crown covered from above by other trees); (1) visible trees (tree crown not or partially covered from above). Most trees classified as category (0) were supressed trees in the understorey and not visible from orthophotos. Therefore, only trees classified in category (1) were taken into consideration for further analyses.

ALS data processing

ALS data were processed in the OPALS environment ([18]) that was developed at the Institute of Photogrammetry and Remote Sensing, Dept. of Geodesy and Geoinformation at the University of Technology in Vienna (Austria). The margin of each plot was extended by 15 m to eliminate the errors occurring at the plot margin during data processing. This step ensured the integrity of crowns partly falling outside the RP. The process of tree identification was carried out at spatial resolutions ranging from 0.5 to 2 m with a 0.5 m step. The required spatial resolution was defined using the Cell module in OPALS. The Cell module derives models through accumulating the selected attribute parameters. Laser height was used as the attribute and the highest attribute value was used as the “max” parameter. Subsequently, two digital surface models (DSM) were created with a grid equal to the spatial resolution that was set in the previous step. To enable DSM creation, non-ground points were imported into OPALS and stored in the OPALS data manager (ODM - for more information see [19] and [14]). The first DSM was created using the interpolation method of moving planes with the quadrant selection mode. Hereby a grid dataset was created representing the sigma 0 of the grid post adjustment (i.e., the standard deviation of the unit weight observation). The second DSM was created by the aggregation of the max values for laser height, which was set as an attribute. The final DSM was created by combining these two DSMs and the sigma 0 raster layer using the Algebra module. The Algebra module is designed to derive a new grid dataset by combining multiple input grid datasets ([18]). The logical condition is explained by the relation (eqn. 1):

\begin{equation} r[1] \;\;\text{&&}\;\; r[2]<d \;\;\text{?}\:\: r[1] : r[0] \end{equation}

where r[0] is the DSM derived by the moving planes interpolation method; r[1] is the DSM with an aggregated parameter; r[2] is the raster layer representing sigma 0; and d is the defined value.

The DTM was derived using the moving planes interpolation method with the quadrant selection mode. The grid has the same resolution as the DSM. Interpolation was made from the ground points stored in the ODM. The normalised DSM (nDSM) was created by subtracting the DSM and the DTM. Statistical filtering was performed on the nDSM using the StatFilter module from OPALS. The shape of the kernel environment was a circle and the max value was set as parameter, i.e., the cell with highest value inside the defined area, where the output grid corresponds to the input grid. The result was the identification of local maxima in the RP. Grids representing the nDSM and local maxima were combined using the Algebra module from OPALS. The result was a grid representing the positions of the trees in the RP. The cell values on this grid represent the height of individual trees.

Accuracy assessment

To assess the accuracy of tree identification using ALS data, two methods were used. The first one consisted of comparing the number of identified trees using ALS data with the number of trees measured in the field survey. The second method involved counting the omission and commission errors. To accomplish this, a reference layer of the crowns for all RP was created through vectorisation using the orthophotos with a spatial resolution of 10 cm. Three dimensional vectorization was carried out on the PLANAR system in the software environment Summit Evolution™ ver. 6.4 (DAT/EM Systems International, Anchorage, AL, USA). A commission error (C) occurs when the number of tree tops identified through ALS data is higher that the tree crowns in the reference layer, or when other objects (different from trees) were identified as trees. Contrastingly, an omission error (O) occurs when a tree top present in the reference layer was not identified based on ALS data. The detection ratio (DR) and accuracy index (AI) were calculated ([22], [29]) as follows (eqn. 2, eqn. 3):

\begin{equation} DR={\frac{N-O}{N}} \cdot 100 \end{equation}
\begin{equation} AI={\frac{N-O-C}{N}} \cdot 100 \end{equation}

where N is the total number of trees in the plot, O is the number of omission errors, and C is the number of commission errors. The accuracy of tree identification by counting the omission and commission errors was assessed only for the spatial resolution which gave the best results in terms of tree identification.

Calculation of dendrometric characteristics

Mean stand height and diameter were calculated separately for broadleaved and coniferous species. The reason for this was the high occurrence of beech and spruce.

Stand height

Two dendrometric parameters were calculated based on ALS data: (i) the average height, calculated as the average height of all trees identified in the RP; (ii) the maximum height found within each RP. Reference data were used to calculate four different parameters for each RP ([25]): (i) the mean height (hd), determined as the height corresponding to the diameter based on Weise’s rule; (ii) the dominant height calculated as the average height of the thickest 10% of trees in the RP (h10%); (iii) the dominant height calculated as the average height of the thickest 20% of trees in the RP (h20%); (iv) the average height (hA), calculated as the average height of all trees in the RP. These four parameters were calculated using only the trees classified in the visibility category 1 (visible trees - see above). Finally, the mean height values determined using ALS data were compared with those obtained by the reference data from the field measurement.

Stand diameter

The diameter from the ALS data was derived according to the regression models developed from the reference data. The variables in the regression analysis were height and DBH of trees measured in each RP during field measurement. Regression analysis was performed separately on broadleaves and conifers by pooling data from all the RPs.

ALS-derived diameter values were obtained using the above regression equation and ALS-derived height of individual trees as predictor, and then averaged over each RP. Additionally, the reference data was used to calculate the following three diameters in each RP ([25]): (i) the mean diameter (dw), according to the Weise’s rule; (ii) the dominant diameter based on the thickest 10% of trees in the RP (d10%); (iii) the dominant diameter based on the thickest 20% of trees in the RP (d20%). These three dendrometric parameters were calculated using only trees belonging to visibility category 1 (visible trees). Finally, the diameter values determined using ALS data were compared with those obtained from the reference data.

Results 

Tab. 2 reports the number of trees identified in the four research plots based on ALS data at four different spatial resolutions, along with the reference number of trees obtained from field measurements. The number of trees identified using ALS data decreases as the spatial resolution increases from 0.5 to 2 m. This is because at 0.5 m of spatial resolution the nDSM is extremely detailed, while a spatial resolution of 2 m yielded a “smoothed” nDSM. At the resolution of 0.5 m, the number of trees determined from the ALS data is extremely overestimated. The most suitable resolution for tree identification was 1.0 m for plots RP2 and RP3, and 1.5 m for RP1 and RP4. In these cases, the accuracy in tree identification was equal to or above 84%. The expectation that conifer trees are more easily identified than broadleaves based on ALS data was not confirmed here. Indeed, considering just conifers, there was an overestimation of 2 trees in RP1 using ALS data. In RP2 and RP3 the number of trees was underestimated by 17 and 9, respectively. Considering just broadleaved species, the estimate obtained was correct in RP4, while in RP1, RP2 and RP3 the number of trees was underestimated by 16, 1 and 11, respectively. The accuracy in tree identification by counting the omission and commission errors is shown in Tab. 3. The assessment was carried out for a spatial resolution of 1.0 m in the case of plots RP2 and RP3, and 1.5 m for plots RP1 and RP4.

Tab. 2 - Number of trees identified in the four research plots (RP) based on ALS data at different spatial resolutions and their comparison with the number of trees based on field measurements (Reference data). (‡): Most accurate results used in the accuracy analysis (see Tab. 3).
Tab. 3 - Accuracy of tree identification using ALS data at the spatial resolutions yielding the most accurate number of trees (see Tab. 2). (RP): research plot.

Comparing the height values derived from ALS data with those obtained from the reference dataset (Tab. 4), the error increases as the spatial resolution decreases (from 0.5 to 2 m). The same trend was also observed with the average height, except for conifers in RP1. The range of errors in the mean height (at resolution of 0.5 m) was -4.5 to 10.4%. For conifers the error range was -4.5 to 2.1%, and was larger for broadleaves (-3.1 to 10.4%). The largest error was in RP4 which is the oldest stand, with big crowns showing complicated surfaces. In this plot, the errors of the average height ranged from -6.1 to 18.9% (at resolution of 0.5 m). For conifers, the error in RP1 was -6.1% at a resolution of 0.5 m, and for all other resolutions the error was ± 1%. Hence, using a low resolution led to height values derived from the ALS data affected by larger errors as compared with those obtained using a high spatial resolution. However, this trend was not recorded in RP2 and RP4, where (in both cases for broadleaves) the best results were obtained using a spatial resolution of 0.5 m.

Tab. 4 - Mean tree height (in m) calculated from ALS data at different spatial resolutions (0.5 to 2.0 m) compared with tree height values obtained from reference data. (RP): research plot; (h10% ): dominant height based on the thickest 10% of trees in the RP; (h20%): dominant height based on the thickest 20% of trees in the RP.

The errors across all resolutions ranged from -19.6 to 16.2% for h10% and from -16.3 to 15.3% for h20%. Underestimation of h10% occurred at all spatial resolutions and ranged from 3.3% at the lowest resolution to 18.7% at the highest. As for height h20%, underestimation occurred in RP1 and RP2 (ranging from 1.7% at the lowest to 16.3% at the highest resolution), whereas in RP3 underestimation only occurred when a spatial resolution of 0.5 m was used. To determine h10% the most suitable spatial resolution was 2 m, resulting in an error range of -10.5 to -3.3%. Regarding h20%, the 2 m data resolution performed better in RP1 and RP2, but in RP3 the resolution of 1.5 m gave better results, with errors ranging from -8.5 to 2.2%. Determination of h10% and h20% in broadleaves did not always result in dependencies as is the case in conifers. In RP1 the error increased from the highest to the lowest spatial resolution, whereas the opposite trend was detected in RP4. In RP2 and RP3 there was no increase in errors from the highest spatial resolution to the lowest nor vice versa. Comparing the distribution of height values derived from ALS data with those obtained from reference data, we observed similar ranges in the case of conifers (Fig. 2). Contrastingly, close height intervals were observed for broadleaved trees only in RP1, while larger differences were recorded in the other plots.

Fig. 2 - Tree height values calculated from reference data and tree height derived from ALS data at spatial resolution from 0.5 to 2 m.

Fig. 3 shows the equations used to predict the tree diameter using ALS-derived individual height. Comparing the diameters obtained using ALS data with the diameters from reference data (Tab. 5), a trend similar to that observed for height estimation was observed. Indeed, errors in diameter estimation increased from the highest to the lowest spatial resolutions. Considering the parameters d10% and d20% the trend was the opposite. In the case of conifers, ALS-derived diameter values obtained at resolutions giving the most accurate identification of the tree number, closely approximated the mean diameter obtained from reference data. In the case of broadleaved species, the results were contradictory, as ALS-derived diameter was close to d20% in RP1, RP2 and RP4, while it was closer to the mean diameter from reference data in RP3. As expected, errors in diameter estimation using ALS data were greater than those observed for the height determination. For the mean diameter (dw), 29% of errors were <10% of the reference mean and 32% of errors were >20%. For the diameter d20%, 36% of errors were <10% and 35% of errors were >20%. Regarding d10% the situation was different: only 11% of the errors were lower than 10% from the reference mean, and 60% were greater than 20%. The errors in the mean diameter ranged from -12.1 to 15.3% at a resolution of 0.5 m (-12.1 to 5.3% for conifers, up to 15.3% for broadleaves). In the case of d10% the diameter was underestimated by 7-41%. These results were similar for both conifers and broadleaves in RP1, while better results were obtained for broadleaves in RP2 and for conifers in RP3. In the case of d20% of broadleaved species, overestimation was detected only in RP3 at a spatial resolution of 2 m, and close estimates were obtained in RP2 at a spatial resolution of 1.5 m. In all other cases d20% was underestimated. The range of errors was between -33.3 to 3.4%. Better results were obtained for broadleaves, where the error range at the spatial resolution of 2 m was between -4.2 to 3.4%. In the case of conifers, errors ranged from -10.4 to 21% at the spatial resolution of 2 m. A comparison of the diameter ranges derived from ALS data with reference data is shown in Fig. 4. In the case of conifers, there was an underestimation of the diameter, and there were no ALS-derived diameter values close to the reference maximum diameter. In the case of broadleaves, the situation was different in each RP.

Fig. 3 - Regression analysis between height and diameter based on data from field measurement. The equations used for deriving the diameter of individual trees (DBH) using the ALS-derived tree height (h) are shown for conifer and broadleaved species.
Tab. 5 - Mean tree diameter (in cm) calculated from ALS data at different spatial resolutions (0.5 to 2.0 m) compared with tree diameter values obtained from reference data. (RP): research plot; (d10%): dominant diameter based on the thickest 10% of trees in the RP; (d20%): dominant diameter based on the thickest 20% of trees in the RP.
Fig. 4 - Tree diameters calculated from reference data and tree height derived from ALS data for spatial resolution from 0.5 to 2 m.

Discussion 

The basic step in deriving stand characteristics from ALS-collected data is identifying the trees, counting their number and locating their position. In this study ALS data were processed using an automatic approach and a spatial resolution ranging from 0.5 to 2 m. Our results confirmed the conclusions of Mikita et al. ([16]) that the spatial resolution of the raster layers significantly influences the identification of trees as well as the estimation of tree and stand dendrometric characteristics. In this study, when using the highest spatial resolution (0.5 m), the number of trees was overestimated by several hundred percent. This is due to the fact that many local maxima are usually identified when the crown surface obtained from the 3D cloud is highly detailed. On the other hand, using a low spatial resolution involves a generalization or “smoothing” of the DSM, leading to underestimate the number of extant trees on the ground.

The accuracy of tree identification was strictly dependent on the spatial resolution used: 1.0 m was the best resolution in RP 2 and RP 3, whereas 1.5 m was best in RP 1 and RP 4. The highest accuracy in each RP was above 84%. Vauhkonen et al. ([32]) stated that the accuracy of tree identification in Scandinavia and Central Europe is more than 70%. In broadleaved stands, the identification accuracy is lower (50%-60%) due to the complex nature of stand canopy and structure. Jing et al. ([8]) achieved a tree identification accuracy of 83% in mixed stands and 66% in broadleaved stands. The detection ratio for the spatial resolutions giving the highest number of identified trees (1.0 and 1.5 m) was in the range 69-75%. These results are similar to those achieved by Kaartinen et al. ([10]) and Vauhkonen et al. ([32]) testing several algorithms, while Zawawi et al. ([35]) obtained an identification accuracy of 24-56%.

Regarding the tree detection ratio, Reitberger et al. ([23]) reported values ranging from 56-94%, while Sterenczak & Miscicki ([28]) achieved a detection ratio of almost 82% for spruce. Yu et al. ([34]) obtained the highest detection ratio of 96% (though with an average value of 69%), and Khosravipour et al. ([11]) reported a mean value of 82%.

The vertical structure of the stand also affects tree identification using ALS-derived data. In this study, we considered only dominant trees with the crown in the upper canopy layer, as the identification of trees in the understorey layers using ALS-data is often difficult. Previous studies reported a percentage of trees identified in the understorey varying from 3-21% ([23]) to 40-45% ([21], [5]). On the other hand, Strîmbu & Strîmbu ([29]) stated that 95% of the trees in the understorey layer were identified in their study.

Tree height derived from ALS data is expected be underestimated since the laser impulses do not always reflect from the highest point of the canopy ([5], [33]). On the other hand, several studies reported an overestimation of tree height using ALS data. For example, Farid et al. ([4]) mentioned an overestimation of 77% when determining the tree height of young cottonwood (Populus fremontii S. Wats.) and mesquite (Prosopis velutina Wooton). The authors assume that the overestimation was caused by the reflection of laser impulses from higher trees located close to each other. Heurich ([5]) achieved a mean deviation between the height derived from ALS data and ground measurements of -0.54 m. Wezyk et al. ([33]) obtained slightly better results, from -0.12 to -0.9 m.

Kaartinen & Hyyppä ([9]) compared several methods from various countries, obtaining a tree height deviation in the range from -1.53 to 3.88 m. Similar results can be found also in Kaartinen et al. ([10]), where the models compared in the study of Kaartinen & Hyyppä ([9]) were extended. In the same study, it is mentioned that the best models achieved a RMSE of 0.6-0.8 m. Zawawi et al. ([35]) stated that the lowest estimated value of tree height from ALS data was 5.5 ± 1.5m. Mikita et al. ([16]) achieved a RMSE of 1.4-1.6 m, after excluding the trees in the understorey layer, while Latifi et al. ([12]) achieved a relative RMSE of 5.5%. In this study, the mean height derived from the ALS data (using the spatial resolution that resulted in the most accurate identification of trees, 1 m and 1.5 m) had an average deviation from reference data of 8.7%. The average differences were lower for h20% (5.9%) and higher for h10% (9.7%). Our results showed a lower precision when compared with those obtained by the aforementioned authors, though they are within the range indicated by Kaartinen & Hyyppä ([9]) and Kaartinen et al. ([10]). Further, we found that the differences in the identified height among different spatial resolutions ranged from 0 to 36%. The smallest differences in mean height were obtained using higher resolutions, and using lower resolutions in the case of h10% and h20%.

Tree height can be straightly obtained from the ALS data. On the contrary, the diameter has to be derived using existing relationships between tree height and diameter ([5], [2]). According to Persson et al. ([20]), it is possible to achieve an estimation of the diameter from ALS data with a difference of 10%. Latifi et al. ([12]) achieved a relative RMSE of 7.2% and a bias of -0.1% when estimating tree diameter from ALS data. Yu et al. ([34]) reported a RMSE of 10.32% for two coniferous wood species (Picea abies Karst. and Pinus sylvestris L.), while Järnstedt et al. ([7]) obtained a worse RMSE (25.3%) for the same tree species using a denser point cloud. In the present study, the diameter was derived from the regression models obtained using reference data collected in the individual plots. The mean diameter derived from ALS data (at the spatial resolution giving the most accurate tree identification - 1 and 1.5 m) was overestimated by 15.8%, on average, while d20% and d10% where underestimated by 14.1% and 21.7%, respectively. The difference in the mean diameter using a resolution of 0.5 m ranged from -8.7 to 10.7%. In general, however, the average difference of diameters is higher than the 10% mentioned by Persson et al. ([20]). Absolute differences in the derived diameters across all spatial resolutions range from 0 to 48%. As already mentioned, the smallest differences was detected using higher spatial resolutions for the mean diameter and lower resolutions for d10% and d20%.

Conclusions 

Airborne laser scanning (ALS) is being more and more applied in forestry for tree detection and measurement. However, spatial resolution of ALS data can deeply affect the estimates of tree and stand characteristics. Using low spatial resolutions (0.5 m) we found that the number of trees in the stands is largely overestimated, while the opposite holds for high spatial resolutions (2.0 m). At intermediate resolutions of ALS data (1.0 m, 1.5 m) tree identification accuracy was over 84%.

The expectation that the most suitable spatial resolution for tree identification would be also the best for determining their height was not fulfilled. Compared with reference data (field measurements), ALS-derived mean height at resolution 0.5 m showed a deviation ranging from -3.1% to 10.4%. Optimal resolution for h10% and h20% estimation in conifers was 2.0 m, while contradictory results were obtained for broadleaved species. ALS-derived tree diameters were inferred using the regression analysis based on reference data. A spatial resolution of 0.5 m was the most suitable for determining mean tree diameter, with errors ranging from -8.7 to 10.7%. However, this method is not suitable for estimating d10% and d20%, as revealed by the large deviations found (sometimes >20%).

In this study, we applied an automatic approach for tree identification and estimation of tree dendrometric characteristics, which requires minimum interventions of the operator. Our results demonstrate that this approach is feasible and leads to mean tree height and diameter estimates which are consistent with similar studies. The main disadvantage is the need of a more generalized model for the relationships between tree diameter and height to be used for ALS-derived diameter estimation. Further analyses of larger datasets are needed to further test the feasibility of the adopted approach for the detection of tree dendrometric characteristics from ALS data.

Acknowledgements 

This research was supported by the projects: Centre of Excellence “Decision support in forest and country”, ITMS: 26220 120069, supported by the Research and Development Operational Programme funded by the ERDF (50%); and by the Slovak Research and Development Agency, project APVV-15-0393 “Innovations in the forest inventories based on progressive technologies of remote sensing” (50%).

References

(1)
Andersen HE, Reutebuch SE, McGaughey RJ (2006). A rigorous assessment of tree height measurements obtained using airborne lidar and conventional field methods. Canadian Journal of Remote Sensing 32 (5): 355-366.
::CrossRef::Google Scholar::
(2)
Balenović I, Alberti G, Marjanović H (2013). Airborne laser scanning - the status and perspectives for the application in the south-east European forestry. South-East European Forestry 4 (2): 59-79.
::CrossRef::Google Scholar::
(3)
Brodrechtová Y, Navrátil R, Sedmák R, Tuček J (2016). Using the politicized IAD framework to assess integrated forest management decision making in Slovakia. Land Use Policy. [In press]
::Google Scholar::
(4)
Farid A, Goodrich DC, Sorooshian S (2006). Using airborne lidar to discern age classes of cottonwood trees in a riparian area. Western Journal of Applied Forestry 21 (3): 149-158.
::Online::Google Scholar::
(5)
Heurich M (2008). Automatic recognition and measurement of single trees based on data from airborne laser scanning over the richly structured natural forests of the Bavarian Forest National Park. Forest Ecology and Management 255 (7): 2416-2433.
::CrossRef::Google Scholar::
(6)
Hyyppä J, Hyyppä H, Xiaowei Y, Kaartinen H, Kukko A, Holopainen M (2009). Forest inventory using small-footprint airborne LiDAR. In: “Topographic Laser Ranging and Scanning: Principles and Processing” (Shan J, Toth CK eds). CRC Press, Boca Raton, FL, USA, pp. 335-370.
::Google Scholar::
(7)
Järnstedt J, Pekkarinen A, Tuominen S, Ginzler C, Holopainen M, Viitala R (2012). Forest variable estimation using a high-resolution digital surface model. ISPRS Journal of Photogrammetry and Remote Sensing 74: 78-84.
::CrossRef::Google Scholar::
(8)
Jing L, Hu B, Li J, Noland T (2014). Automated individual tree crown delineation from LiDAR data using morphological techniques. IOP Conference Series: Earth and Environmental Science 17: 1.
::CrossRef::Google Scholar::
(9)
Kaartinen H, Hyyppä J (2008). Tree extraction - Report of EuroSDR project. Official publication no. 53, Frankfurt, Germany, pp. 60.
::Google Scholar::
(10)
Kaartinen H, Hyyppä J, Yu X, Vastaranta M, Hyyppä H, Kukko A, Holopainen M, Heipke C, Hirschmugl M, Mosdorf F, Pitkänen J, Popescu S, Solberg S, Wolf BM, Wu JC (2012). An international comparison of individual tree detection and extraction using airborne laser scanning. Remote Sensing 4: 950-974.
::CrossRef::Google Scholar::
(11)
Khosravipour A, Skidmore AK, Wang T, Isenburg M, Khoshelham K (2015). Effect of slope on treetop detection using a LiDAR canopy height model. ISPRS Journal of Photogrammetry and Remote Sensing 104: 44-52.
::CrossRef::Google Scholar::
(12)
Latifi H, Fassnacht EF, Müller J, Tharani A, Dech S, Heurich M (2015). Forest inventories by LiDAR data: a comparison of single tree segmentation and metric-based methods for inventories of a heterogeneous temperate forest. International Journal of Applied Earth Observation and Geoinformation 42: 162-174.
::CrossRef::Google Scholar::
(13)
Maas HG (2013). Forestry applications. In: “Airborne and Terrestrial Laser Scanning” (Vosselman G, Maas H-G eds). Whittles Publishing, Dunbeath, UK, pp. 213-235.
::Google Scholar::
(14)
Mandlburger G, Otepka J, Karel W, Wagner W, Pfeifer N (2009). Orientation and processing of airborne laser scanning data (OPALS) - Concept and first results of a comprehensive ALS software. In: Proceedings of the “Symposium of ISPRS - Laserscanning 2009” (Bretar F, Pierrot-Deseiligny M, Vosselman G eds). Paris (France) 1-2 Sept 2009. Paris, France, pp. 55-60.
::Online::Google Scholar::
(15)
Michnová Z (2015). Odvodenie vybraných stromových a porastových veličín metódami leteckého laserového skenovania [Determining selected tree and forest stand parameters using methods of airborne laser scanning]. PhD thesis, Department of Forest Management and Geodesy, Technical University in Zvolen, Zvolen, Slovakia, pp. 113. [in Slovak]
::Google Scholar::
(16)
Mikita T, Klimánek M, Cibulka M (2013). Evaluation of airborne laser scanning data for tree parameters and terrain modelling in forest environment. Acta Universitatis Agriculturae et Silviculturae Mendelianae Brunensis 61 (5): 1339-1347.
::CrossRef::Google Scholar::
(17)
Nelson R, Swift R, Krabill W (1988). Using airborne lasers to estimate forest canopy and stand characteristics. Journal of Forestry 86 (10): 31-38.
::Google Scholar::
(18)
OPALS Manual (2016). Orientation and processing of airborne laser scanning data: user documentation. Web site.
::Online::Google Scholar::
(19)
Otepka J, Briese C, Nothegger C (2006). First steps to a topographic information system of the next generation. In: Proceedings of the “Symposium of ISPRS Commission IV - Geo Spatial Databases for Sustainable Development”. Goa (India) 25-30 Sept 2006. Goa, India, pp. 419-425.
::Online::Google Scholar::
(20)
Persson A, Holmgren J, Söderman U (2002). Detecting and measuring individual trees using an airborne laser scanner. Photogrammetric Engineering and Remote Sensing 68 (9): 925-932.
::Online::Google Scholar::
(21)
Pitkänen J, Maltamo M, Hyyppä J, Yu X (2004). Adaptive methods for individual tree detection on airborne laser based canopy height model. In: Proceedings of the “Symposium of ISPRS - Laser-Scanners for Forest and Landscape Assessment” (Thies M, Koch B, Spiecker H, Weinacker H eds). Freiburg, (Germany) 3-6 Oct 2004. Freiburg, Germany, pp. 178-183.
::Online::Google Scholar::
(22)
Pouliot DA, King DJ, Bell FW, Pitt DG (2002). Automated tree crown detection and delineation in high-resolution digital camera imagery of coniferous forest regeneration. Remote Sensing of Environment 82: 322-334.
::CrossRef::Google Scholar::
(23)
Reitberger J, Heurich M, Krzystek P, Stilla U (2007). Single tree detection in forest areas with high-density LiDAR data. International Archives of Photogrammetry, Remote Sensing and Spatial Information Science 36 (3/W49B): 139-144.
::Online::Google Scholar::
(24)
Sedmák R, Fabrika M, Bahyl J, PôbiS I, Tuček J (2013). Application of simulation and optimazation tools for developing forest management plans in the Slovak natural and management conditions. In: “Implementation of DSS tools into the forestry practice” (Tuček J, Smreček R, Majlingová A, Garzio-Gonzalo J eds). Technická univerzita vo Zvolene, Zvolen, Slovakia, pp. 139-158.
::Google Scholar::
(25)
Smelko S (2000). Dendrometria [Forest mensuration]. Technical University in Zvolen, Zvolen, Slovakia, pp. 399. [in Slovak]
::Google Scholar::
(26)
Smreček R, Danihelová Z (2013). Forest stand height determination from low point density airborne laser scanning data in Rožnava Forest enterprise zone (Slovakia). iForest 6: 48-54.
::CrossRef::Google Scholar::
(27)
Sterenczak K (2013). Factors influencing individual tree crowns detection based on airborne laser scanning data. Forest Research Papers 74 (4): 323-333.
::CrossRef::Google Scholar::
(28)
Sterenczak K, Miscicki S (2012). Crown delineation influence on standing volume calculations in protected area. In: Proceedings of the “Symposium of ISPRS Commission VIII - Remote Sensing and Spatial Information Sciences (Shortis M, Shimoda H, Cho K eds). Melbourne (Australia) 25 Aug - 01 Sept 2012. Melbourne, Australia, pp. 441-445.
::Online::Google Scholar::
(29)
Strîmbu VF, Strîmbu BM (2015). A graph-based segmentation algorithm for tree crown extraction using airborne LiDAR data. ISPRS Journal of Photogrammetry and Remote Sensing 104: 30-43.
::CrossRef::Google Scholar::
(30)
Takahashi T, Yamamoto K, Senda Y, Tsuzuku M (2005). Estimating individual-tree heights of sugi (Cryptomeria japonica D. Don) plantations in mountainous areas using small-footprint airborne LiDAR. Journal of Forest Research 10 (4): 135-142.
::CrossRef::Google Scholar::
(31)
Vastaranta M, Kankare V, Holopainen M, Yu X, Hyyppä J, Hyyppä H (2012). Combination of individual tree detection and area-based approach in imputation of forest variables using airborne laser data. ISPRS Journal of Photogrammetry and Remote Sensing 67: 73-79.
::CrossRef::Google Scholar::
(32)
Vauhkonen J, Ene L, Gupta S, Heinzel J, Holmgren J, Pitkänen J, Solberg S, Wang Y, Weinacker H, Hauglin KM, Lien V, Packalén P, Gobakken T, Koch B, Naesset E, Tokola T, Maltamo M (2012). Comparative testing of single-tree detection algorithms under different types of forest. Forestry 85 (1): 27-40.
::CrossRef::Google Scholar::
(33)
Wezyk P, Tompalski P, Szostak M, Glista M, Pierzchalski M (2008). Describing the selected canopy layer parameters of the scots pine stands using ALS data. In: Proceedings of the “SilviLaser 2008 - 8th International Conference on LiDAR Applications in Forest Assessment and Inventory” (Hill R, Rosette J, Suárez J eds). Edinburgh (UK) 17-19 Sept 2008. Edinburg, UK, pp. 636-645.
::Online::Google Scholar::
(34)
Yu X, Hyyppä J, Holopainen M, Vastaranta M (2010). Comparison of area-based and individual tree-based methods for predicting plot-level forest attributes. Remote Sensing 2 (6): 1481-1495.
::CrossRef::Google Scholar::
(35)
Zawawi AA, Shiba M, Jemali NJN (2015). Accuracy of LiDAR-based height estimation and crown recognition in a subtropical evergreen broad-leaved forest in Okinawa, Japan. Forest Systems 24 (1): 2276-2287.
::Google Scholar::
(36)
Zíhlavník S, Scheer L (2000). Dialkový prieskum Zeme v lesníctve [Remote sensing of Earth in forestry]. Technical University in Zvolen, Zvolen, Slovakia, pp. 289. [in Slovak]
::Google Scholar::

Download

Paper Contents

Paper Sections

Paper Figures

Paper Tables

 
 
Close

 

Smreček R, Michnová Z, Sačkov I, Danihelová Z, Levická M, Tuček J (2018).
Determining basic forest stand characteristics using airborne laser scanning in mixed forest stands of Central Europe
iForest - Biogeosciences and Forestry 11: 181-188. - doi: 10.3832/ifor2520-010
Close
First Previous Next Last
 
Close
© iForest

Download Reference

Paper ID# ifor2520-010
Title Determining basic forest stand characteristics using airborne laser scanning in mixed forest stands of Central Europe
Authors Smreček R, Michnová Z, Sačkov I, Danihelová Z, Levická M, Tuček J
Format
Close Download