Search

This study presents a new approach for predicting forest aboveground biomass (AGB) from airborne laser scanning (ALS) data: AGB is predicted from sequences of images depicting vertical cross-sections through the ALS point clouds. A 3D version of the VGG16 convolutional neural network (CNN) with initial weights transferred from pre-training on the ImageNet dataset was used. The approach was tested on datasets from Canada, Poland, and the Czech Republic. To analyse the effect of training sample size on model performance, different-sized samples ranging from 10 to 375 ground plots were used. The CNNs were compared with random forest models (RFs) trained on point cloud metrics. At the maximum number of training samples, the difference in RMSE between observed and predicted AGB of CNNs and RFs ranged from −2 t/ha to 5 t/ha, and the difference in squared Pearson correlation coefficient ranged from −0.05 to 0.06. Additional pre-training on synthetic data derived from virtual laser scanning of simulated forest stands could only improve the prediction performance of the CNNs when only a few real training samples (10–40) were available. While 3D CNNs trained on cross-section images derived from real data showed promising results, RFs remain a competitive alternative.

Fusing of terrestrial laser scanning (TLS) and airborne laser scanning (ALS) point cloud data has been recognized as an effective approach in forest studies. In this regard, co-registration of point clouds is considered one of the crucial steps in the integration process. Co-registering point clouds in forest environments faces various challenges, including unstable features, extensive occlusions, different viewpoints, and differences in point cloud densities. To address these intricate challenges, this study introduces an automated and robust method for co-registering TLS and ALS point clouds based on the correspondence of individual tree locations in forest environments. Initially, the positions of individual trees in both TLS and ALS data are extracted. Then, a filtering approach is applied to eliminate positions with low potential for corresponding matches in the TLS and ALS dataset. Since larger trees in the TLS data have a higher potential for corresponding matches in the ALS data, an iterative process is applied to identify correspondences between trees in both datasets. After estimating transformation parameters, the co-registration process is executed. The proposed method is applied on six datasets with varying forest complexities. The results demonstrate a high success rate up to 100% if the starting position of the TLS plots are located within ∼4 hectares (∼2000 trees). Additionally, the potential of the proposed method for co-registering TLS data with ALS data across different search areas and varying number of trees is evaluated in detail. The outcomes indicate that successful co-registration of TLS plot with 50 m diameter to ALS data is successful in the best case within a search radius of approximately 113 hectares (∼60,000 tree locations) and in the worst case for around 20 hectares (∼10,000 tree locations) depending on the forest complexity.

Abstract A key aspect of nature conservation is knowledge of which aspects of nature to conserve or restore to favor the characteristic diversity of plants in a given area. Here, we used a large plant dataset with >40 000 plots combined with airborne laser scanning (lidar) data to reveal the local characteristics of habitats having a high plant dark diversity—that is, absence of suitable species—at national extent (>43 000 km2). Such habitats have potential for reaching high realized diversity levels and hence are important in a conservation context. We calculated 10 different lidar based metrics (both terrain and vegetation structure) and combined these with seven different field‐based measures (soil chemistry and species indicators). We then used Integrated Nested Laplace Approximation for modelling plant dark diversity across 33 North European habitat types (open landscapes and forests) selected by the European communities to be important. In open habitat types high‐dark‐diversity habitats had relatively low pH, high nitrogen content, tall homogenous vegetation, and overall relatively homogenous terrains (high terrain openness) although with a rather high degree of local microtopographical variations. High‐dark‐diversity habitats in forests had relatively tall vegetation, few natural‐forest indicators, low potential solar radiation input and a low cover of small woody plants. Our results highlight important vegetation, terrain‐ and soil‐related factors that managers and policymakers should be aware of in conservation and restoration projects to ensure a natural plant diversity, for example low nutrient loads, natural microtopography and possibly also open forests with old‐growth elements such as dead wood and rot attacks.

Terrestrial LiDAR scanning (TLS) has the potential to revolutionize forestry by enabling the precise estimation of aboveground biomass, vital for forest carbon management. This study addresses the lack of comprehensive benchmarking for leaf-filtering algorithms used in TLS data processing and evaluates four widely recognized geometry-based leaf-filtering algorithms (LeWoS, TLSeparation, CANUPO, and a novel random forest model) across openly accessible TLS datasets from diverse global locations. Multiple evaluation dimensions are considered, including pointwise classification accuracy, volume comparisons using a quantitative structure model applied to wood points, computational efficiency, and visual validation. The random forest model outperformed the other algorithms in pointwise classification accuracy (overall accuracy = 0.95 ± 0.04), volume comparison (R-squared = 0.96, slope value of 0.98 compared to destructive volume), and resilience to reduced point cloud density. In contrast, TLSeparation exhibits the lowest pointwise classification accuracy (overall accuracy = 0.81 ± 0.10), while LeWoS struggles with volume comparisons (mean absolute percentage deviation ranging from 32.14 ± 29.45% to 49.14 ± 25.06%) and point cloud density variations. All algorithms show decreased performance as data density decreases. LeWoS is the fastest in terms of processing time. This study provides valuable insights for researchers to choose appropriate leaf-filtering algorithms based on their research objectives and forest conditions. It also hints at future possibilities for improved algorithm design, potentially combining radiometry and geometry to enhance forest parameter estimation accuracy.

Mapping and monitoring of habitats are requirements for protecting biodiversity. In this study, we investigated the benefit of combining airborne (laser scanning, image-based point clouds) and satellite-based (Sentinel 1 and 2) data for habitat classification. We used a two level random forest 10-fold leave-location-out cross-validation workflow to model Natura 2000 forest and grassland habitat types on a 10 m pixel scale at two study sites in Vienna, Austria. We showed that models using combined airborne and satellite-based remote sensing data perform significantly better for forests than airborne or satellite-based data alone. For frequently occurring classes, we reached class accuracies with F1-scores from 0.60 to 0.87. We identified clear difficulties of correctly assigning rare classes with model-based classification. Finally, we demonstrated the potential of the workflow to identify errors in reference data and point to the opportunities for integration in habitat mapping and monitoring.

Increasing river floods and infrastructure development in many parts of the world have created an urgent need for accurate high-resolution flood hazard mapping for more efficient flood risk management. Mapping accuracy hinges on the quality of the underlying Digital Terrain Model (DTM) and other spatial datasets. This article presents a processing strategy to ensure consistent adaption of countrywide spatial datasets to the requirements of hydraulic modelling. The suggested methods are automatized to a large extent and include (i) automatic fitting of river axis positions to the DTM, (ii) detection of culverts and obstacles in the river channel (iii) Smooth elimination of obstacles by interpolation along the river axes (iv) geometric detection of water-land borders and the top edge of embankments for (v) integration of the submerged river bed geometry into the DTM. The processing chain is applied to a river network (33880 km) and a DTM from Airborne Laser Scanning (ALS) with 1 m spatial resolution covering the entire territory of Austria (∼84000 km2). Thus, countrywide consistency of data and methods is achieved along with high local relevance. Semi-automatic validation and extensive manual checks demonstrate that processing significantly improves the DTM with respect to topographic and hydraulic consistency. However, some open issues of automatic processing remain, e.g. in case of long underground river reaches.

Microwave and optical imaging methods react differently to different land surface parameters and, thus, provide highly complementary information. However, the contribution of individual features from these two domains of the electromagnetic spectrum for tree species classification is still unclear. For large-scale forest assessments, it is moreover important to better understand the domain-specific limitations of the two sensor families, such as the impact of cloudiness and low signal-to-noise-ratio, respectively. In this study, seven deciduous and five coniferous tree species of the Austrian Biosphere Reserve Wienerwald (105,000 ha) were classified using Breiman’s random forest classifier, labeled with help of forest enterprise data. In nine test cases, variations of Sentinel-1 and Sentinel-2 imagery were passed to the classifier to evaluate their respective contributions. By solely using a high number of Sentinel-2 scenes well spread over the growing season, an overall accuracy of 83.2% was achieved. With ample Sentinel-2 scenes available, the additional use of Sentinel-1 data improved the results by 0.5 percentage points. This changed when only a single Sentinel-2 scene was supposedly available. In this case, the full set of Sentinel-1-derived features increased the overall accuracy on average by 4.7 percentage points. The same level of accuracy could be obtained using three Sentinel-2 scenes spread over the vegetation period. On the other hand, the sole use of Sentinel-1 including phenological indicators and additional features derived from the time series did not yield satisfactory overall classification accuracies (55.7%), as only coniferous species were well separated.

Timber assortments are some of the most important goods provided by forests worldwide. To quantify the amount and type of timber assortment is strongly important for socio-economic purposes, but also for accurate assessment of the carbon stored in the forest ecosystems, regardless of their main function. Terrestrial laser scanning (TLS) became a promising tool for timber assortment assessment compared to the traditional surveys, allowing reconstructing the tree architecture directly and rapidly. This study aims to introduce an approach for timber assortment assessment using TLS data in a mixed and multi-layered Mediterranean forest. It consists of five steps: (1) pre-processing, (2) timber-leaf discrimination, (3) stem detection, (4) stem reconstruction, and (5) timber assortment assessment. We assume that stem form drives the stem reconstruction, and therefore, it influences the timber assortment assessment. Results reveal that the timber-leaf discrimination accuracy is 0.98 through the Random Forests algorithm. The overall detection rate for all trees is 84.4%, and all trees with a diameter at breast height larger than 0.30 m are correctly identified. Results highlight that the main factors hindering stem reconstruction are the presence of defects outside the trunk, trees poorly covered by points, and the stem form. We expect that the proposed approach is a starting point for valorising the timber resources from unmanaged/managed forests, e.g., abandoned forests. Further studies to calibrate its performance under different forest stand conditions are furtherly required.

With the increasing occurrence of forest fires in the mid-latitudes and the alpine region, fire risk assessments become important in these regions. Fuel assessments involve the collection of information on forest structure as, e.g., the stand height or the stand density. The potential of airborne laser scanning (ALS) to provide accurate forest structure information has been demonstrated in several studies. Yet, flight acquisitions at the state level are carried out in intervals of typically five to ten years in Central Europe, which often makes the information outdated. The Sentinel-1 (S-1) synthetic aperture radar mission provides freely accessible earth observation (EO) data with short revisit times of 6 days. Forest structure information derived from this data source could, therefore, be used to update the respective ALS descriptors. In our study, we investigated the potential of S-1 time series to derive stand height and fractional cover, which is a measure of the stand density, over a temperate deciduous forest in Austria. A random forest (RF) model was used for this task, which was trained using ALS-derived forest structure parameters from 2018. The comparison of the estimated mean stand height from S-1 time series with the ALS derived stand height shows a root mean square error (RMSE) of 4.76 m and a bias of 0.09 m on a 100 m cell size, while fractional cover can be retrieved with an RMSE of 0.08 and a bias of 0.0. However, the predictions reveal a tendency to underestimate stand height and fractional cover for high-growing stands and dense areas, respectively. The stratified selection of the training set, which we investigated in order to achieve a more homogeneous distribution of the metrics for training, mitigates the underestimation tendency to some degree, yet, cannot fully eliminate it. We subsequently applied the trained model to S-1 time series of 2017 and 2019, respectively. The computed difference between the predictions suggests that large decreases in the forest height structure in this two-year interval become apparent from our RF-model, while inter-annual forest growth cannot be measured. The spatial patterns of the predicted forest height, however, are similar for both years (Pearson’s R = 0.89). Therefore, we consider that S-1 time series in combination with machine learning techniques can be applied for the derivation of forest structure information in an operational way.

The measurements of tree attributes required for forest monitoring and management planning, e.g., National Forest Inventories, are derived by rather time-consuming field measurements on sample plots, using calipers and measurement tapes. Therefore, forest managers and researchers are looking for alternative methods. Currently, terrestrial laser scanning (TLS) is the remote sensing method that provides the most accurate point clouds at the plot-level to derive these attributes from. However, the demand for even more efficient and effective solutions triggers further developments to lower the acquisition time, costs, and the expertise needed to acquire and process 3D point clouds, while maintaining the quality of extracted tree parameters. In this context, photogrammetry is considered a potential solution. Despite a variety of studies, much uncertainty still exists about the quality of photogrammetry-based methods for deriving plot-level forest attributes in natural forests. Therefore, the overall goal of this study is to evaluate the competitiveness of terrestrial photogrammetry based on structure from motion (SfM) and dense image matching for deriving tree positions, diameters at breast height (DBHs), and stem curves of forest plots by means of a consumer grade camera. We define an image capture method and we assess the accuracy of the photogrammetric results on four forest plots located in Austria and Slovakia, two in each country, selected to cover a wide range of conditions such as terrain slope, undergrowth vegetation, and tree density, age, and species. For each forest plot, the reference data of the forest parameters were obtained by conducting field surveys and TLS measurements almost simultaneously with the photogrammetric acquisitions. The TLS data were also used to estimate the accuracy of the photogrammetric ground height, which is a necessary product to derive DBHs and tree heights. For each plot, we automatically derived tree counts, tree positions, DBHs, and part of the stem curve from both TLS and SfM using a software developed at TU Wien (Forest Analysis and Inventory Tool, FAIT), and the results were compared. The images were oriented with errors of a few millimetres only, according to checkpoint residuals. The automatic tree detection rate for the SfM reconstruction ranges between 65% and 98%, where the missing trees have average DBHs of less than 12 cm. For each plot, the mean error of SfM and TLS DBH estimates is −1.13 cm and −0.77 cm with respect to the caliper measurements. The resulting stem curves show that the mean differences between SfM and TLS stem diameters is at maximum −2.45 cm up to 3 m above ground, which increases to almost +4 cm for higher elevations. This study shows that with the adopted image capture method, terrestrial SfM photogrammetry, is an accurate solution to support forest inventory for estimating the number of trees and their location, the DBHs and stem curve up to 3 m above ground.

Structural measures for retaining and distributing water—i.e., reservoirs, flood retention and power plants—play a key role to protect and feed a growing world population in a rapidly changing climate. In this work, we introduce an automated method to detect potential reservoir or retention area locations in digital terrain models. In this context, a potential reservoir is a larger terrain form that can be turned into an actual reservoir by constructing a dam. Based on contour lines derived from terrain models, potential reservoirs are found within a predefined range of dam lengths, and the locally optimal ones are then extracted. Our method is to be applied in the very early stages of project planning and for area-wide potential analysis. Tests in a 100 km2 study area bring promising results, but also show a certain sensitivity regarding terrain model quality and resolution. In total, 250–300 candidate polygons with a total volume of more than 6 million m3 were found. In order to facilitate further processing, these are stored as a GIS vector dataset.

We propose a flexible framework for automated forest patch delineations that exploits a set of canopy structure features computed from airborne laser scanning (ALS) point clouds. The approach is based on an iterative subdivision of the point cloud using k-means clustering followed by an iterative merging step to tackle oversegmentation. The framework can be adapted for different applications by selecting relevant input features that best measure the intended homogeneity. In our study, the performance of the segmentation framework was tested for the delineation of forest patches with a homogeneous canopy height structure on the one hand and with similar water cycle conditions on the other. For the latter delineation, canopy components that impact interception and evapotranspiration were used, and the delineation was mainly driven by leaf area, tree functional type, and foliage density. The framework was further tested on two scenes covering a variety of forest conditions and topographies. We demonstrate that the delineated patches capture well the spatial distributions of relevant canopy features that are used for defining the homogeneity. The consistencies range from R 2 = 0 . 84 to R 2 = 0 . 86 and from R 2 = 0 . 80 to R 2 = 0 . 91 for the most relevant features in the delineation of patches with similar height structure and water cycle conditions, respectively.

Extracting 3D tree models based on terrestrial laser scanning (TLS) point clouds is a challenging task as trees are complex objects. Current TLS devices acquire high-density data that allow a detailed reconstruction of the tree topology. However, in dense forests a fully automatic reconstruction of trees is often limited by occlusion, wind influences and co-registration issues. In this paper, a semi-automatic method for extracting branching and stem structure based on equirectangular projections (range and intensity maps) is presented. The digitization of branches and stems is based on 2D maps, which enables simple navigation and raster processing. The modeling is performed for each viewpoint individually instead of using a registered point cloud. Previously reconstructed 2D-skeletons are transformed between the maps. Therefore, wind influences, orientation imperfections of scans and data gaps can be overcome. The method is applied to a TLS dataset acquired in a forest in Germany. In total 34 scans were carried out within a managed forest to measure approximately 90 spruce trees with minimal occlusions. The results demonstrate the feasibility of the presented approach to extract tree models with a high completeness and correctness and provide an excellent input for further modeling applications.

This paper presents a methodology for the derivation of structural parameters and stem volume in forests based on Airborne Laser Scanning (ALS) data. We describe three different measures of horizontal and vertical canopy structure: (1) tree crown segmentation, (2) compactness of vegetation patches, and (3) vertical layering of vegetation patches and canopy cover. An empirical regression model for the derivation of stem volume from the ALS and forest inventory sample plot data is described and its results are validated with extensive reference data. Different study areas in Austria were used to illustrate the workflows. The presented study demonstrates the applicability of the proposed methods on study sites and ALS data of differing characteristics, as well as it points out the suitability of ALS as a tool for reliable wide area assessment of structural parameters and stem volume for forested areas.ResumoEsse artigo apresenta uma metodologia para derivação de parâmetros estruturais e de volume de madeira em florestas baseado em dados de Laser Scanner Aerotransportado (ALS). Nós descrevemos três diferentes medidas da estrutura horizontal e vertical da copa: (1) segmentação da copa da árvore, (2) compacidade das manchas de vegetação, (3) estratificação vertical das manchas de vegetação e cobertura do dossel. Um modelo empírico de regressão para derivar o volume de madeira fazendo uso de dados ALS e dados amostrais obtidos em inventário florestal é descrito e seus resultados são validados com extensivos dados de referência. Diferentes áreas na Áustria foram utilizadas para ilustrar o fluxo de trabalho. O estudo apresentado demonstra a aplicabilidade dos métodos propostos nas áreas de estudo e dos dados ALS de diferentes características, bem como aponta a adequação do ALS como ferramenta confiável para avaliação de parâmetros de estrutura e de volume de madeira de amplas áreas florestais.

This study compares methods to estimate stem volume, stem number and basal area from Airborne Laser Scanning (ALS) data for 68 field plots in a hemi-boreal, spruce dominated forest (Lat. 58°N, Long. 13°E). The stem volume was estimated with five different regression models: one model based on height and density metrics from the ALS data derived from the whole field plot, two models based on similar combinations derived from 0.5 m raster cells, and two models based on canopy volumes from the ALS data. The best result was achieved with a model based on height and density metrics derived from 0.5 m raster cells (Root Mean Square Error or RMSE 37.3%) and the worst with a model based on height and density metrics derived from the whole field plot (RMSE 41.9%). The stem number and the basal area were estimated with: (i) area-based regression models using height and density metrics from the ALS data; and (ii) single tree-based information derived from local maxima in a normalized digital surface model (nDSM) mean filtered with different conditions. The estimates from the regression model were more accurate (RMSE 52.7% for stem number and 21.5% for basal area) than those derived from the nDSM (RMSE 63.4%–91.9% and 57.0%–175.5%, respectively). The accuracy of the estimates from the nDSM varied depending on the filter size and the conditions of the applied filter. This suggests that conditional filtering is useful but sensitive to the conditions.

The delineation of forested areas is a critical task, because the resulting maps are a fundamental input for a broad field of applications and users. Different national and international forest definitions are available for manual or automatic delineation, but unfortunately most definitions lack precise geometrical descriptions for the different criteria. A mandatory criterion in forest definitions is the criterion of crown coverage (CC), which defines the proportion of the forest floor covered by the vertical projection of the tree crowns. For loosely stocked areas, this criterion is especially critical, because the size and shape of the reference area for calculating CC is not clearly defined in most definitions. Thus current forest delineations differ and tend to be non-comparable because of different settings for checking the criterion of CC in the delineation process. This paper evaluates a new approach for the automatic delineation of forested areas, based on airborne laser scanning (ALS) data with a clearly defined method for calculating CC. The new approach, the ‘tree triples’ method, is based on defining CC as a relation between the sum of the crown areas of three neighboring trees and the area of their convex hull. The approach is applied and analyzed for two study areas in Tyrol, Austria. The selected areas show a loosely stocked forest at the upper timberline and a fragmented forest on the hillside. The fully automatic method presented for delineating forested areas from ALS data shows promising results with an overall accuracy of 96%, and provides a beneficial tool for operational applications.

Roughness is an important input parameter for modeling of natural hazards such as floods, rock falls and avalanches, where it is basically assumed that flow velocities decrease with increasing roughness. Seeing roughness as a multi-scale level concept (i.e., ranging from fine-scale soil characteristics to description of understory and lower tree layer) various roughness raster products were derived from the original full-waveform airborne laser scanning (FWF-ALS) point cloud using two different types of roughness parameters, the surface roughness (SR) and the terrain roughness (TR). For the calculation of the SR, ALS terrain points within a defined height range to the terrain surface are considered. For the parameterization of the SR, two approaches are investigated. In the first approach, a geometric description by calculating the standard deviation of plane fitting residuals of terrain points is used. In the second one, the potential of the derived echo widths are analyzed for the parameterization of SR. The echo width is an indicator for roughness and the slope of the target. To achieve a comparable spatial resolution of both SR layers, the calculation of the standard deviation of detrended terrain points requires a higher terrain point density than the SR parameterization using the echo widths. The TR describes objects (i.e., point clusters) close but explicitly above the terrain surface, with 20 cm defined as threshold height value for delineation of the surface layer (i.e., forest floor layer). Two different empirically defined vegetation layers below the canopy layer were analyzed (TR I: 0.2 m to 1.0 m; TR II: 0.2 m to 3.0 m). A 1 m output grid cell size was chosen for all roughness parameters in order to provide consistency for further integration of high-resolution optical imagery. The derived roughness parameters were then jointly classified, together with a normalized Digital Surface Model (nDSM) showing the height of objects (i.e., trees) above ground. The presented approach enables the classification of forested areas in patches of different vegetation structure (e.g., varying soil roughness, understory, density of natural cover). For validation purposes in situ reference data were collected and cross-checked with the classification results, positively confirming the general feasibility of the proposed vertical concept of integrated roughness mapping on various vertical levels. Results can provide valuable input for forest mapping and monitoring, in particular with regard to natural hazard modeling.

In this study, a semi-empirical model that was originally developed for stem volume estimation is used for aboveground biomass (AGB) estimation of a spruce dominated alpine forest. The reference AGB of the available sample plots is calculated from forest inventory data by means of biomass expansion factors. Furthermore, the semi-empirical model is extended by three different canopy transparency parameters derived from airborne LiDAR data. These parameters have not been considered for stem volume estimation until now and are introduced in order to investigate the behavior of the model concerning AGB estimation. The developed additional input parameters are based on the assumption that transparency of vegetation can be measured by determining the penetration of the laser beams through the canopy. These parameters are calculated for every single point within the 3D point cloud in order to consider the varying properties of the vegetation in an appropriate way. Exploratory Data Analysis (EDA) is performed to evaluate the influence of the additional LiDAR derived canopy transparency parameters for AGB estimation. The study is carried out in a 560 km2 alpine area in Austria, where reference forest inventory data and LiDAR data are available. The investigations show that the introduction of the canopy transparency parameters does not change the results significantly according to R2 (R2 = 0.70 to R2 = 0.71) in comparison to the results derived from, the semi-empirical model, which was originally developed for stem volume estimation.

This work focuses on the accuracy estimation of canopy height models (CHMs) derived from image matching of Pléiades stereo imagery over forested mountain areas. To determine the height above ground and hence canopy height in forest areas, we use normalised digital surface models (nDSMs), computed as the differences between external high-resolution digital terrain models (DTMs) and digital surface models (DSMs) from Pléiades image matching. With the overall goal of testing the operational feasibility of Pléiades images for forest monitoring over mountain areas, two questions guide this work whose answers can help in identifying the optimal acquisition planning to derive CHMs. Specifically, we want to assess (1) the benefit of using tri-stereo images instead of stereo pairs, and (2) the impact of different viewing angles and topography. To answer the first question, we acquired new Pléiades data over a study site in Canton Ticino (Switzerland), and we compare the accuracies of CHMs from Pléiades tri-stereo and from each stereo pair combination. We perform the investigation on different viewing angles over a study area near Ljubljana (Slovenia), where three stereo pairs were acquired at one-day offsets. We focus the analyses on open stable and on tree covered areas. To evaluate the accuracy of Pléiades CHMs, we use CHMs from aerial image matching and airborne laser scanning as reference for the Ticino and Ljubljana study areas, respectively. For the two study areas, the statistics of the nDSMs in stable areas show median values close to the expected value of zero. The smallest standard deviation based on the median of absolute differences (σMAD) was 0.80 m for the forward-backward image pair in Ticino and 0.29 m in Ljubljana for the stereo images with the smallest absolute across-track angle (−5.3°). The differences between the highest accuracy Pléiades CHMs and their reference CHMs show a median of 0.02 m in Ticino with a σMAD of 1.90 m and in Ljubljana a median of 0.32 m with a σMAD of 3.79 m. The discrepancies between these results are most likely attributed to differences in forest structure, particularly tree height, density, and forest gaps. Furthermore, it should be taken into account that temporal vegetational changes between the Pléiades and reference data acquisitions introduce additional, spurious CHM differences. Overall, for narrow forward–backward angle of convergence (12°) and based on the used software and workflow to generate the nDSMs from Pléiades images, the results show that the differences between tri-stereo and stereo matching are rather small in terms of accuracy and completeness of the CHM/nDSMs. Therefore, a small angle of convergence does not constitute a major limiting factor. More relevant is the impact of a large across-track angle (19°), which considerably reduces the quality of Pléiades CHMs/nDSMs.

Airborne laser scanning (ALS) is a remote sensing technique well-suited for 3D vegetation mapping and structure characterization because the emitted laser pulses are able to penetrate small gaps in the vegetation canopy. The backscattered echoes from the foliage, woody vegetation, the terrain, and other objects are detected, leading to a cloud of points. Higher echo densities (> 20 echoes/m2) and additional classification variables from full-waveform (FWF) ALS data, namely echo amplitude, echo width and information on multiple echoes from one shot, offer new possibilities in classifying the ALS point cloud. Currently FWF sensor information is hardly used for classification purposes. This contribution presents an object-based point cloud analysis (OBPA) approach, combining segmentation and classification of the 3D FWF ALS points designed to detect tall vegetation in urban environments. The definition tall vegetation includes trees and shrubs, but excludes grassland and herbage. In the applied procedure FWF ALS echoes are segmented by a seeded region growing procedure. All echoes sorted descending by their surface roughness are used as seed points. Segments are grown based on echo width homogeneity. Next, segment statistics (mean, standard deviation, and coefficient of variation) are calculated by aggregating echo features such as amplitude and surface roughness. For classification a rule base is derived automatically from a training area using a statistical classification tree. To demonstrate our method we present data of three sites with around 500,000 echoes each. The accuracy of the classified vegetation segments is evaluated for two independent validation sites. In a point-wise error assessment, where the classification is compared with manually classified 3D points, completeness and correctness better than 90% are reached for the validation sites. In comparison to many other algorithms the proposed 3D point classification works on the original measurements directly, i.e. the acquired points. Gridding of the data is not necessary, a process which is inherently coupled to loss of data and precision. The 3D properties provide especially a good separability of buildings and terrain points respectively, if they are occluded by vegetation.

Abstract: Airborne laser scanning (ALS) is an active remote sensing technique that uses the time-of-flight measurement principle to capture the three-dimensional structure of the earth’s surface with pulsed lasers that transmit nanosecond-long laser pulses with a high pulse repetition frequency. Over forested areas most of the laser pulses are reflected by the leaves and branches of the trees, but a certain fraction of the laser pulses reaches the forest floor through small gaps in the canopy. Thus it is possible to reconstruct both the three-dimensional structure of the forest canopy and the terrain surface. For the retrieval of quantitative forest parameters such as stem volume or biomass it is necessary to use models that combine ALS with inventory data. One approach is to use multiplicative regression models that are trained with local inventory data. This method has been widely applied over boreal forest regions, but so far little experience exists with applying this method for mapping alpine forest. In this study the transferability of this approach to a 128 km2 large mountainous region in Vorarlberg, Austria, was evaluated. For the calibration of the model, inventory data as operationally collected by Austrian foresters were used. Despite these inventory data are based on variable sample plot sizes, they could be used for mapping stem volume for the entire alpine study area. The coefficient of determination R2 was 0.85 and the root mean square error (RMSE) 90.9 m3ha-1 (relative error of 21.4%) which is comparable to results of ALS studies conducted over topographically less complex environments. Due to the increasing availability, ALS data could become an operational part of Austrian’s forest inventories.

Radar sensors have the potential to retrieve vertical forest structure measurements thanks to their capability to penetrate into the foliage. However, studies are needed in order to understand better the interaction of radar beams with the canopy. The most commonly used radar technique for estimating forest parameters operates from spacecraft at different wavelength (X-, C-, and L-band). In order to assist in the interpretation of satellite data for forest applications, and as a possible complementary technique to Lidar (Light detection and ranging), the Finnish Geospatial Research Institute has developed the first helicopter-borne profiling radar system operating in Ku-band, called Tomoradar, which is able to provide a vertical canopy profile. The study focuses on the analyses of Ku-band profiling radar waveforms and the backscatter signal of boreal forest, supported by simultaneously acquired Lidar measurements, in order to detect ground and canopy profiles and quantify the ground echo ratio under different canopy coverage and the backscatter signal variation through the vegetation. The Tomoradar data was acquired simultaneously with a lightweight Velodyne VLP-16 Lidar system in October 2016 over a boreal forest located in Evo in southern Finland. Additionally, highly accurate Riegl VQ-480 Lidar data, acquired in 2014, was used as a ground reference for both lightweight systems. We analysed the co- and cross-polarized (HH and HV) Tomoradar backscatter signals of a 600 m long profile. It is found that the Ku-band Tomoradar penetrates the canopy to a similar extent as the Velodyne Lidar, i.e., the distribution of backscatter signals through the vegetation follows the vegetation density. Moreover, the ground backscatter signal strength and ground echo ratio depend strongly on the presence of gaps in the canopy. By comparing the elevation of the corresponding canopy and ground Tomoradar signal peaks with the Velodyne Lidar data, the Tomoradar ground elevation accuracy is on average −0.03 m and −0.20 m for the cross- and co-polarization, respectively, whereas the bias of the canopy elevation is, on average, −0.58 m and 1.35 m for the cross- and co-polarization, respectively. With respect to the ground height data derived from the Lidar measurements of 2014, the Tomoradar ground profile reveals, on average, higher accuracy (i.e., 0.00 m (σ = 0.41 m) and 0.04 m (σ = 0.37 m) for the co-and cross-polarizations, respectively) than the Velodyne system (−0.37 m with σ = 0.25 m).

Diameter at breast height (DBH) is one of the most important parameter in forestry. With increasing use of terrestrial and airborne laser scanning in forestry, new exceeding possibilities to directly derive DBH emerge. In particular, high resolution point clouds from laser scanners on board unmanned aerial systems (UAS) are becoming available over forest areas. In this case study, DBH estimation from a UAS point cloud based on modeling the relevant part of the tree stem with a cylinder, is analyzed with respect to accuracy and completeness. As reference, manually measured DBHs and DBHs from terrestrial laser scanning point clouds are used for comparison. We demonstrate that accuracy and completeness of the cylinder fit are depending on the stem diameter. Stems with DBH > 20 cm feature almost 100% successful reconstruction with relative differences to the reference DBH of 9% (DBH 20–30 cm) down to 1.8% for DBH > 40 cm.

Terrestrial laser scanning (TLS) is a promising technique for plot-wise acquisition of geometric attributes of forests. However, there still exists a need for TLS applications in mountain forests where tree stems’ growing directions are not vertical. This paper presents a novel method to model tree stems precisely in an alpine landslide-affected forest using TLS. Tree stems are automatically detected by a two-layer projection method. Stems are modeled by fitting a series of cylinders based on a 2D-3D random sample consensus (RANSAC)-based approach. Diameter at breast height (DBH) was manually measured in the field, and stem curves were measured from the point cloud as reference data. The results showed that all trees in the test area can be detected. The root mean square error (RMSE) of estimated DBH was 1.80 cm (5.5%). Stem curves were automatically generated and compared with reference data, as well as stem volumes. The results imply that the proposed method is able to map and model the stem curve precisely in complex forest conditions. The resulting stem parameters can be employed in single tree biomass estimation, tree growth quantification and other forest-related studies.