Estimating forest stand density and structure using Bayesian individual tree detection, stochastic geometry, and distribution matching

Individual tree detection (ITD) and the area-based approach (ABA) are combined to generate tree-lists using airborne LiDAR data. ITD based on the Canopy Height Model (CHM) was applied for overstory trees, while ABA based on nearest neighbor (NN) imputation was applied for understory trees. Our approach is intended to compensate for the weakness of LiDAR data and ITD in estimating understory trees, keeping the strength of ITD in estimating overstory trees in tree-level. We investigated the effects of three parameters on the performance of our proposed approach: smoothing of CHM, resolution of CHM, and height cutoff (a specific height that classifies trees into overstory and understory). There was no single combination of those parameters that produced the best performance for estimating stems per ha, mean tree height, basal area, diameter distribution and height distribution. The trees in the lowest LiDAR height class yielded the largest relative bias and relative root mean squared error. Although ITD and ABA showed limited explanatory powers to estimate stems per hectare and basal area, there could be improvements from methods such as using LiDAR data with higher density, applying better algorithms for ITD and decreasing distortion of the structure of LiDAR data. Automating the procedure of finding optimal combinations of those parameters is essential to expedite forest management decisions across forest landscapes using remote sensing data.

Horvitz–Thompson-like stand density estimation is a method for estimating the stand density from tree crown objects extracted from airborne laser scanning data through individual tree detection. The estimator is based on stochastic geometry and mathematical morphology of the (planar) set formed by the detected tree crowns. This set is used to approximate the detection probabilities of trees. These probabilities are then used to calculate the estimate. The method includes a tuning parameter, which needs to be known to apply the method. We present a refinement of the method to allow more general detection conditions than those of previous papers. We also present and discuss the methods for estimating the tuning parameter of the estimator using a functional k-nearest neighbors method. We test the model fitting and prediction in two spatially separate data sets and examine the plot-level accuracy of estimation. The estimator produced a 13% lower RMSE (root-mean-square error) than the benchmark method in an external validation data set. We also analyze the effects of similarity and dissimilarity of training and validation data on the results.

More than 50% of the national lands in Japan have been surveyed by airborne laser scanning (ALS) data with different point densities; and developing an effective approach to take full advantage of these ALS data for forest management has thus become an urgent topic of study. This study attempted to assess the utility of ALS data for individual tree detection and species classification in a mixed forest with a high canopy density. For comparison, two types of tree tops and tree crowns in the study area were delineated by the individual tree crown (ITC) approach using the green band of the orthophoto imagery and the digital canopy height model (DCHM) derived from the ALS data, respectively. Then, the two types of tree crowns were classified into four classes—Pinus densiflora (Pd), Chamaecyparis obtusa (Co), Larix kaempferi (Lk), and broadleaved trees (Bl)—by a crown-based classification approach using different combinations of the three orthophoto bands with intensity and slope maps as follows: RGB (red, green and blue); RGB and intensity (RGBI); RGB and slope (RGBS); and RGB, intensity and slope (RGBIS). Finally, the tree tops were annotated with species attributes from the two best-classified tree crown maps, and the number of different tree species in each compartment was counted for comparison with the field data. The results of our study suggest that the combination of RGBIS yielded greater classification accuracy than the other combinations. In the tree crown classifications delineated by the green band and DCHM data, the improvements in the overall accuracy compared to the RGB ranged from 5.7% for the RGBS to 9.0% for the RGBIS and from 8.3% for the RGBS to 11.8% for the RGBIS. The laser intensity and slope derived from the ALS data may be valuable sources of information for tree species classification, and in terms of distinguishing species for the detection of individual trees, the findings of this study demonstrate the advantages of using DCHM instead of optical data to delineate tree crowns. In conclusion, the synthesis of individual tree delineation using DCHM data and species classification using the RGBIS combination is recommended for interpreting forest resources in the study area. However, the usefulness of this approach must be verified in future studies through its application to other forests.

Ecology and forestry sciences are using an increasing amount of data to address a wide variety of technical and research questions at the local, continental and global scales. However, one type of data remains rare: fine-grain descriptions of large landscapes. Yet, this type of data could help address the scaling issues in ecology and could prove useful for testing forest management strategies and accurately predicting the dynamics of ecosystem services. Here we present three datasets describing three large European landscapes in France, Poland and Slovenia down to the tree level. Tree diameter, height and species data were generated combining field data, vegetation maps and airborne laser scanning (ALS) data following an area-based approach. Together, these landscapes cover more than 100 000 ha and consist of more than 42 million trees of 51 different species. Alongside the data, we provide here a simple method to produce high-resolution descriptions of large landscapes using increasingly available data: inventory and ALS data. We carried out an in-depth evaluation of our workflow including, among other analyses, a leave-one-out cross validation. Overall, the landscapes we generated are in good agreement with the landscapes they aim to reproduce. In the most favourable conditions, the root mean square error (RMSE) of stand basal area (BA) and mean quadratic diameter (Dg) predictions were respectively 5.4 m 2.ha -1 and 3.9 cm, and the generated main species corresponded to the observed main species in 76.2% of cases.

Bridging spatio-temporal gaps in ALS data using Landsat time series and forest disturbance-recovery metrics via multi-task neural networks

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.

These maps depict the distribution of 12 tree species across the state of New York. The maps show where these trees do not occur (gray), occasionally occur (pale green), are a minor component (medium green), are a major component (dark green), or are the dominant species (black) in the forest, as determined by that species' total basal area. Basal area is the area of a cross-section of the trunk at 4.5 feet above ground (breast height). Basal area is a way of measuring how dominant a particular species is in a stand because of the way large trees contribute more to the total basal area than small trees. The map at the top of the circle (red maple) is the species with the greatest basal area in the state. All other maps are generally arranged to group species that commonly occur together. Of the 105 tree species recorded in New York in the 2005-2009 inventories, the species presented here are the top 12 in the state by total basal area, and together represent 74 percent of the total live tree basal area and 75 percent of the number of trees in New York. The center map shows where forest, nonforest, and water are present. In all maps, white is nonforest and water is blue. High resolution version of this map is available for download (113 MB PDF) Dataset: Live tree species basal area of the contiguous United States (2000-2009)

Updating stand-level forest inventories using airborne laser scanning and Landsat time series data

Demonstrating the transferability of forest inventory attribute models derived using airborne laser scanning data

Assessment of UAV photogrammetric DTM-independent variables for modelling and mapping forest structural indices in mixed temperate forests

Airborne laser scanning based forest inventories employ two major methods: individual tree detection (ITD) and the area-based statistical approach (ABSA). ITD is based on the assumption that trees are of a certain form and can be delineated using airborne laser scanning techniques, whereas ABSA is an empirical method based on the relations between area-level forest attributes and laser echo height distributions. These two methods are compared here within the same test area in terms of their usefulness for estimating mean forest stand characteristics and tree size distributions. All evaluations were performed using leave-one-out cross validation. The average errors in volume and basal area did not differ significantly between the methods. ABSA resulted in overall better accuracies when estimating the diameter and height of the basal area median tree and the number of stems, whereas ITD produced significantly biased estimates for the number of stems and the mean tree size. Tree size distributions were estimated with slightly better accuracy using ABSA. More comprehensive investigations revealed that both methods were not able to estimate forest structure (tree size distribution and spatial distribution of tree locations), which in turn, affected the estimation accuracies.

Conservation and sustainable use strategies for reliable management plans of valuable species are lacking due to the gap in basic scientific information on their structure and status. The objectives of this work were to study the structural characteristics and the spatial distribution of the abundance of <i>Spondias mombin</i> L. populations in the subequatorial zone of Benin. The data collection was based on the tessellation of study zone with ArcGIS analysis tools. The grid cell used was a square of 100 m side (1ha of surface). The occurrence of species in each grid was obtained using GBIF database. One hundred and seventy-five (175) trees of <i>S. mombin</i> were sampled through the four Phytodistricts of study zone. The main dendrometric variables measured from these individuals were diameter at breast height (1.30 m above ground level) and height. The parameters such as density, mean quadratic diameter, basal area and mean Lorey height were calculated. To study the spatial abundance distribution of <i>S. mombin,</i> we used Inverse Distance Weighting (IDW) interpolation method to realize abundance map from the occurrences. The results of structure showed a density of 3 trees/ha. The basal area was 1500 cm²/ha and the mean quadratic diameter was 25.10 cm. The mean of Lorey height obtained was 12 m<i>.</i> The diameter structure showed a predominance of bell-shape distribution with a relative abundance of young individuals. The low numbers of individuals in some diameter classes, the erratic distribution and the presence of old population of <i>S. mombin</i> indicate the strong anthropogenic pressure on the species. The abundance map showed that the species was present in all the phytodistricts and covers 93% of the study area’s municipalities. The conservation strategies, including both the species and its habitat, need to be developed to allow sustainable use of <i>S. mombin</i>.

Airborne laser scanning (ALS; LiDAR) data are an increasingly common data source for forest inventories, and approaches integrating ALS data with field plot measurements have become operational in several jurisdictions. As technology continues to evolve, different LiDAR sensors can provide new opportunities to incorporate LiDAR data into forest inventory workflows. Single photon LiDAR (SPL) enables efficient, large area data acquisition and merits further investigation for forest inventory applications. Herein, we investigated the capacity of leaf-on SPL data, combined with 269 field plots, for estimating forest inventory attributes in the Great Lakes–St. Lawrence mixedwood forests of southern Ontario, Canada. Inventory attribute estimates were validated at the stand level using independent reference data acquired for 27 intensively sampled stands. Top height, Lorey’s height, gross total volume for merchantable stems, merchantable stem volume, basal area, quadratic mean diameter, and total aboveground biomass were estimated with a relative RMSE of 13.52%, 7.24%, 14.61%, 16.27%, 14.42%, 12.25%, and 11.72%, respectively. Relative bias was &lt; 1% for all attributes except top height (10.34%), merchantable volume (3.37%), and basal area (1.68%). Accuracy and bias varied by forest type and stand-level validation was important for assessing model performance in different stand conditions. SPL data can be used to generate accurate, area-based forest inventories in mixedwood forests that have a multitude of tree species and complex forest management histories.

Trees are an integral component of the urban environment and important for human well-being, adaption measures to climate change and sustainable urban transformation. Understanding the small-scale impacts of urban trees and strategically managing the ecosystem services they provide requires high-resolution information on urban forest structure, which is still scarce. In contrast, there is an abundance of data portraying urban areas and an associated trend towards smart cities and digital twins as analysis platforms. A GIS workflow is presented in this paper that may close this data gap by classifying the urban forest from LiDAR point clouds, detecting and reconstructing individual crowns, and enabling a tree representation within semantic 3D city models. The workflow is designed to provide robust results for point clouds with a density of at least 4 pts/m2 that are widely available. Evaluation was conducted by mapping the urban forest of Dresden (Germany) using a point cloud with 4 pts/m². An object-based data fusion approach is implemented for the classification of the urban forest. A classification accuracy of 95 % for different urban settings is achieved by combining LiDAR with multispectral imagery and a 3D building model. Individual trees are detected by local maxima filtering and crowns are segmented using marker-controlled watershed segmentation. Evaluation highlights the influences of both urban and forest structure on individual tree detection. Substantial differences in detection accuracies are evident between trees along streets (72 %) and structurally more complex tree stands in green areas (31 %), as well as dependencies on tree height and crown diameter. Furthermore, an approach for parameterized reconstruction of tree crowns is presented, which enables efficient and realistic city-wide modeling. The suitability of LiDAR to measure individual tree metrics is illustrated as well as a framework for modeling individual tree crowns via geometric primitives.

We present a decision support tool for guiding the selection of marked stands based on airborne laser scanning (ALS) data. We describe three stages, namely (1) wall-to-wall mapping of the stands matured for cutting using low-density ALS data; (2) tree-level inventory of these stands using high-density ALS data and (3) theoretical bucking of the imputed tree stems to produce detailed information on their characteristics. We tested them in a Scots pine dominated boreal forest area in Eastern Finland, where 79 sample plots were measured in the field. The detection of the stands matured for cutting had a success rate of 95% and our results demonstrated a further potential to limit the result towards stands dominated by certain species by means of intensity values derived from the low-density ALS data. The applied single-tree detection and estimation chain produced detailed tree-level information and realistic diameter distributions, yet the detection was highly emphasised on the dominant tree layer. The error levels in the estimates were generally less than standard deviations of the field attributes. Finally, plot-level accumulations of saw-log volumes were found rather similar, whether the input was based on the imputed tree data or trees measured in the field. The results are considered useful for ranking the stands based on their properties, whether the aim in the wood procurement is to focus on certain species or to select stands suitable for production needs.