Assessing log geometry and wood quality in standing timber using terrestrial laser-scanning point clouds

Wood procurement in sawmills could be improved by resolving detailed three-dimensional stem geometry references from standing timber. This could be achieved, using the increasingly available terrestrial point clouds from various sources. Here, we collected terrestrial laser-scanning (TLS) data from 52 Scots pines (Pinus sylvestris L.) with the purpose of evaluating the accuracy of the log geometry and analysing its relationship with wood quality. For reference, the log-specific top-end diameter, volume, tapering, sweep, basic density and knottiness were measured in a sawmill. We produced stem models from the TLS data and bucked them into logs similar to those measured in the sawmill. In comparison to the sawmill data, the log-specific TLS-based top-end diameter, volume, taper and sweep estimates showed relative mean differences of 1.6, -2.4, -3.0 and 78 per cent, respectively. The correlation coefficients between increasing taper and decreasing wood density and whorl-to-whorl distances were 0.49 and -0.51, respectively. Although the stem-model geometry was resolved from the point clouds with similar accuracy to that at the sawmills, the remaining uncertainty in defining the sweep and linking the wood quality with stem geometry may currently limit the method's feasibilities. Instead of static TLS, mobile platforms would likely be more suitable for operational point cloud data acquisition.Peer reviewe

On the use of rapid-scan, low point density terrestrial laser scanning (TLS) for structural assessment of complex forest environments

Forests fulfill an important role in natural ecosystems, e.g., they provide food, fiber, habitat, and biodiversity, all of which contribute to stable ecosystems. Assessing and modeling the structure and characteristics in forests can lead to a better understanding and management of these resources. Traditional methods for collecting forest traits, known as “forest inventory”, is achieved using rough proxies, such as stem diameter, tree height, and foliar coverage; such parameters are limited in their ability to capture fine-scale structural variation in forest environments. It is in this context that terrestrial laser scanning (TLS) has come to the fore as a tool for addressing the limitations of traditional forest structure evaluation methods. However, there is a need for improving TLS data processing methods. In this work, we developed algorithms to assess the structure of complex forest environments – defined by their stem density, intricate root and stem structures, uneven-aged nature, and variable understory - using data collected by a low-cost, portable TLS system, the Compact Biomass Lidar (CBL). The objectives of this work are listed as follow: 1. Assess the utility of terrestrial lidar scanning (TLS) to accurately map elevation changes (sediment accretion rates) in mangrove forest; 2. Evaluate forest structural attributes, e.g., stems and roots, in complex forest environments toward biophysical characterization of such forests; and 3. Assess canopy-level structural traits (leaf area index; leaf area density) in complex forest environments to estimate biomass in rapidly changing environments. The low-cost system used in this research provides lower-resolution data, in terms of scan angular resolution and resulting point density, when compared to higher-cost commercial systems. As a result, the algorithms developed for evaluating the data collected by such systems should be robust to issues caused by low-resolution 3D point cloud data. The data used in various parts of this work were collected from three mangrove forests on the western Pacific island of Pohnpei in the Federated States of Micronesia, as well as tropical forests in Hawai’i, USA. Mangrove forests underscore the economy of this region, where more than half of the annual household income is derived from these forests. However, these mangrove forests are endangered by sea level rise, which necessitates an evaluation of the resilience of mangrove forests to climate change in order to better protect and manage these ecosystems. This includes the preservation of positive sediment accretion rates, and stimulating the process of root growth, sedimentation, and peat development, all of which are influenced by the forest floor elevation, relative to sea level. Currently, accretion rates are measured using surface elevation tables (SETs), which are posts permanently placed in mangrove sediments. The forest floor is measured annually with respect to the height of the SETs to evaluate changes in elevation (Cahoon et al. 2002). In this work, we evaluated the ability of the CBL system for measuring such elevation changes, to address objective #1. Digital Elevation Models (DEMs) were produced for plots, based on the point cloud resulted from co-registering eight scans, spaced 45 degree, per plot. DEMs are refined and produced using Cloth Simulation Filtering (CSF) and kriging interpolation. CSF was used because it minimizes the user input parameters, and kriging was chosen for this study due its consideration of the overall spatial arrangement of the points using semivariogram analysis, which results in a more robust model. The average consistency of the TLS-derived elevation change was 72%, with and RMSE value of 1.36 mm. However, what truly makes the TLS method more tenable, is the lower standard error (SE) values when compared to manual methods (10-70x lower). In order to achieve our second objective, we assessed structural characteristics of the above-mentioned mangrove forest and also for tropical forests in Hawaii, collected with the same CBL scanner. The same eight scans per plot (20 plots) were co-registered using pairwise registration and the Iterative Closest Point (ICP). We then removed the higher canopy using a normal change rate assessment algorithm. We used a combination of geometric classification techniques, based on the angular orientation of the planes fitted to points (facets), and machine learning 3D segmentation algorithms to detect tree stems and above-ground roots. Mangrove forests are complex forest environments, containing above-ground root mass, which can create confusion for both ground detection and structural assessment algorithms. As a result, we needed to train a supporting classifier on the roots to detect which root lidar returns were classified as stems. The accuracy and precision values for this classifier were assessed via manual investigation of the classification results in all 20 plots. The accuracy and precision for stem classification were found to be 82% and 77%, respectively. The same values for root detection were 76% and 68%, respectively. We simulated the stems using alpha shapes in order to assess their volume in the final step. The consistency of the volume evaluation was found to be 85%. This was obtained by comparing the mean stem volume (m3/ha) from field data and the TLS data in each plot. The reported accuracy is the average value for all 20 plots. Additionally, we compared the diameter-at-breast-height (DBH), recorded in the field, with the TLS-derived DBH to obtain a direct measure of the precision of our stem models. DBH evaluation resulted in an accuracy of 74% and RMSE equaled 7.52 cm. This approach can be used for automatic stem detection and structural assessment in a complex forest environment, and could contribute to biomass assessment in these rapidly changing environments. These stem and root structural assessment efforts were complemented by efforts to estimate canopy-level structural attributes of the tropical Hawai’i forest environment; we specifically estimated the leaf area index (LAI), by implementing a density-based approach. 242 scans were collected using the portable low-cost TLS (CBL), in a Hawaii Volcano National Park (HAVO) flux tower site. LAI was measured for all the plots in the site, using an AccuPAR LP-80 Instrument. The first step in this work involved detection of the higher canopy, using normal change rate assessment. After segmenting the higher canopy from the lidar point clouds, we needed to measure Leaf Area Density (LAD), using a voxel-based approach. We divided the canopy point cloud into five layers in the Z direction, after which each of these five layers were divided into voxels in the X direction. The sizes of these voxels were constrained based on interquartile analysis and the number of points in each voxel. We hypothesized that the power returned to the lidar system from woody materials, like branches, exceeds that from leaves, due to the liquid water absorption of the leaves and higher reflectivity for woody material at the 905 nm lidar wavelength. We evaluated leafy and woody materials using images from projected point clouds and determined the density of these regions to support our hypothesis. The density of points in a 3D grid size of 0.1 m, which was determined by investigating the size of the branches in the lower portion of the higher canopy, was calculated in each of the voxels. Note that “density” in this work is defined as the total number of points per grid cell, divided by the volume of that cell. Subsequently, we fitted a kernel density estimator to these values. The threshold was set based on half of the area under the curve in each of the distributions. The grid cells with a density below the threshold were labeled as leaves, while those cells with a density above the threshold were set as non-leaves. We then modeled the LAI using the point densities derived from TLS point clouds, achieving a R2 value of 0.88. We also estimated the LAI directly from lidar data by using the point densities and calculating leaf area density (LAD), which is defined as the total one-sided leaf area per unit volume. LAI can be obtained as the sum of the LAD values in all the voxels. The accuracy of LAI estimation was found to be 90%. Since the LAI values cannot be considered spatially independent throughout all the plots in this site, we performed a semivariogram analysis on the field-measured LAI data. This analysis showed that the LAI values can be assumed to be independent in plots that are at least 30 m apart. As a result, we divided the data into six subsets, where each of the plots were 30 meter spaced for each subset. LAI model R2 values for these subsets ranged between 0.84 - 0.96. The results bode well for using this method for automatic estimation of LAI values in complex forest environments, using a low-cost, low point density, rapid-scan TLS

Assessing wood properties in standing timber with laser scanning

Managed forests play crucial roles in ongoing climatic and environmental changes. Among other things, wood is capable of sinking and storing carbon in both standing timber and wood products. To promote these positive effects, more precise planning is required that will ensure sustainable forest management and maximal deposition of harvested wood for long-term applications. Information on wood properties plays a key role; i.e. the wood properties can impact the carbon stocks in forests and the suitability of wood for structural timber. With respect to the theoretical background of wood formation, stem, crown, and branching constitute potential inputs (i.e. wood quality indicators) to allometric wood property, tree biomass, and wood quality models. Due to the complex nature of wood formation, measurements of wood quality indicators that could predict wood properties along the relevant directions of variation have previously been elusive in forest inventories. However, developments in laser scanning from aerial and terrestrial platforms support more complex mapping and modeling regimes based on dense three-dimensional point clouds. The aim here was to determine how wood properties could be estimated in remote-sensing-aided forest inventories. For this purpose, methods for characterizing select wood quality indicators in standing timber, using airborne and terrestrial laser scanning (ALS and TLS, respectively) were developed and evaluated in managed boreal Scots pine (Pinus sylvestris L.) forests. Firstly, the accuracies of wood quality indicators resolved from TLS point clouds were assessed. Secondly, the results were compared with x-ray tomographic references from sawmills. Thirdly, the accuracies of tree-specific crown features delineated from the ALS data in predictive modeling of the wood quality indicators were evaluated. The results showed that the quality and density of point clouds significantly impacted the accuracies of the extracted wood quality indicators. In the assessment of wood properties, TLS should be considered as a tool for retrieving as dense stem and branching data as possible from carefully selected sample trees. Accurately retrieved morphological data could be applied to allometric wood property models. The models should use tree traits predictable with aerial remote sensing (e.g. tree height, crown dimensions) to enable extrapolations. As an outlook, terrestrial and aerial remote sensing can play an important role in filling in the knowledge gaps regarding the behavior of wood properties over different spatial and temporal extents. Further interdisciplinary cooperation will be needed to fully facilitate the use of remote sensing and spatially transferable wood property models that could become useful in tackling the challenges associated with changing climate, silviculture, and demand for wood.Hoidetuilla metsillä on useita tärkeitä rooleja muuttuvassa ilmastossa ja ympäristössä. Puu sitoo ja varastoi hiiltä niin kasvaessaan, kuin pitkäikäisiksi puutuotteiksi jalostettuna. Näiden vaikutusten huomioiminen metsänhoidossa vaatii tarkkaa suunnittelua, jolla varmistetaan metsänhoidon ja puunkäytön kestävyys. Tieto puuaineen ominaisuuksista on keskeisessä osassa, sillä ne vaikuttavat hiilivarastojen suuruuteen metsissä, sekä puun käytettävyyteen pitkäikäisenä rakennesahatavarana. Puunmuodostuksen teoreettisen taustan mukaisesti, runko, latvus ja oksarakenne ovat potentiaalisia selittäviä muuttujia (eli puun laatuindikaattoreita), kun mallinnetaan puuaineen ominaisuuksia, puubiomassaa ja puun laatua. Puunmuodostuksen monimutkaisuudesta ja moniulotteisesta vaihtelusta johtuen, tarvittavien laatuidikaattorien mittaaminen osana metsävarojen inventointia ja riittävällä yksityiskohtaisuudella on ollut aiemmin mahdotonta. Monialustaisen laserkeilauksen kehittyminen kuitenkin tukee aiempaa monipuolisempien kartoitus- ja mallinnusjärjestelmien rakentamista, jotka perustuvat tiheisiin kolmiulotteisiin pistepilviin. Tämän työn tavoitteena oli määritellä, kuinka puuaineen ominaisuuksia voidaan arvioida kaukokartoitusta hyödyntävässä metsävarojen inventoinnissa. Tätä tarkoitusta varten kehitettiin menetelmiä puun laatuindikaattorien mittaamiseksi hoidetuissa männiköissä (Pinus sylvestris L.) lento- ja maastolaserkeilauksen avulla, ja arvioitiin niiden toimivuutta. Ensin arvioitiin laatuindikaattorien mittatarkkuus pistepilvissä. Toiseksi verrattiin pistepilvimittauksia röntgentomografiamittauksiin teollisilla sahoilla. Kolmanneksi arvioitiin lentolaserkeilauksella tuotettujen latvuspiirteiden tarkkuutta laatuindikaattorien ennustamisessa. Tuloksien perusteella pistepilvien laatu ja pistetiheys vaikuttivat merkittävästi mitattujen laatuindikaattorien tarkkuuteen. Puuaineen ominaisuuksien arvioimisessa, maastolaserkeilausta tulisi käyttää työkaluna mahdollisimman yksityiskohtaisten runko- ja oksikkuustietojen keräämiseen tarkkaan valikoiduista näytepuista. Tarkasti mitatut laatuindikaattorit voivat selittää puuaineen ominaisuuksia mallinnuksessa. Käytettyjen mallien tulisi perustua laatuindikaattoreille, jotka voidaan ennustaa lentolaserkeilausaineistosta (esim. puun pituus ja latvuksen mittasuhteet), jotta ennusteet ovat yleistettävissä laajoille alueille. Tulevaisuudessa, maasta ja ilmasta tehtävällä kaukokartoituksella voi olla tärkeä rooli puuaineen ominaisuuksien aikaan ja paikkaan sidotun vaihtelun tutkimuksessa. Lisää poikkitieteellistä työtä tarvitaan, jotta kaukokartoitusta ja puuaineen ominaisuuksia ennustavia spatiaalisia malleja voidaan täysimittaisesti hyödyntää kiihtyvän ilmastonmuutoksen, muuttuvan metsänhoidon ja lisääntyvän puunkäytön tuomien haasteiden kohtaamisessa

Novel insights into Mediterranean forest structure using high-resolution remote sensing.

PhD Theses.Tree crown morphology and arrangement in three-dimensional space is a key driver of forest dynamics, determining not only the competitiveness of an individual but also the competitive effect exerted on neighbouring trees. Many theoretical frameworks aim to predict crown morphology from first principles and assumptions of Euclidean form and ultimately infer whole forest stand structure and dynamics but paucity in data has limited vigorous testing. Tree crowns are also not rigid in form and due to their sessile nature, must morphologically adapt to immediate abiotic and biotic surroundings to enhance survival. The characterisation of tree structure has been limited by the simplicity and associated error of traditional crown measurements. This project uses Terrestrial Laser Scanning data collected from a water limited Mediterranean forest community in Spain to highlight methodological opportunities presented by TLS in understanding forest structure and also the various developments required to extract ecologically meaningful metrics from these data. It then applies these novel metrics to answer questions about how tree crowns scale with size, the effects of competition and how plasticity in shape and arrangement interacts with light capture at the individual and plot scales. Modification to existing code as well as bespoke development were required to segment and calculate individual metrics from trees in this forest type. Accurate measures of crown morphology highlighted allometric scaling deviations from theoretical predictions and intra-specific differences in response to competition, calculated using more representative neighbourhood metrics. Inter-specific differences in crown plasticity and significant effects of size (height) were also evident, along with trade-offs between morphological plasticity and crown size. Light capture was positively affected by plasticity with inter-specific differences highlighting various biomass allocations strategies species undertake to acquire light. At the plot scale, mixed-genus plots intercepted less direct light and were structurally more complex rather than more volume filling

Improving TLS-based stem volume estimates by field measurements

The prediction of tree stem volumes has conventionally been based on simple field measurements and applicable allometric functions, but terrestrial laser scanning (TLS) has enabled new opportunities for extracting stem volumes of single trees. TLS-based tree dimensions are commonly estimated by automatized cylinder- or circle-based fitting approaches which, given that the stem cross-sections are relatively round and the whole stem is sufficiently covered by TLS points, enable an accurate prediction of the stem volume. The results are, however, often deteriorated by co-registration errors and occlusions, i.e., incompletely visible parts of the stem, which easily lead to poorly fitted features and problems in locating the actual treetop. As these defects are difficult to be controlled or totally avoided when collecting data at a plot level, taking advantage of additional field measurements is proposed to improve the fitting process and mitigate gross errors in the prediction of stem volumes. In this paper, this is demonstrated by modelling the stems first as cylinders by only using TLS data, after which the results are refined with the assistance of field data. The applied data consists of various field-measured stem dimensions which are used to define the acceptable diameter estimation limits and set the correct vertical extents for the analyzed tree. This approach is tested using two data sets, differing in the scanning setup, location, and the measured field variables. Adding field data improves the results and, at best, enables almost unbiased volumetric predictions with an RMSE of less than 5%. According to these results, combining TLS point clouds and simple field measurements has the potential to produce stem volume information at a considerably higher accuracy than TLS data alone