Prediction and mapping of boreal forest fire fuel loads using high-resolution satellite stereo imagery

In Tunisia, the forests cover 1.1 million hectares. The assessment of their significant role in the biosphere depends on the good understanding of their spatial extension, their quality and their dynamics. To estimate the forest biomass, tree dendrometric measurements such as dominant height are needed. This information must be accurate, and up to date which induces new methodological challenges for the forest inventories programs. To solve this issue, the remote sensing techniques are used and especially the photogrammetry techniques using structure from motion (SfM) algorithms. SfM techniques allow 3D forest canopy modelling by exploiting stereoscopic high-resolution aerial or satellite images. They lead to various products such as the digital surface model (DSM) and the canopy high model (CHM). In this study, a DSM and the corresponding orthophoto were produced, for the first time, at a very high resolution (45 cm for both DSM and orthophoto), on forested areas in North Western of Tunisia.The 3D reconstruction was performed from overlapping aerial photos taken by a DMC camera. The issues are related to the complexity of forest areas for matching algorithms in SfM approaches. The objective of this paper is to propose and optimized multiscale processing framework and to assess the DSM accuracy obtained from overlapping aerial photos processed with structure from motion (SfM) techniques on a Mediterranean forest

ABSTRACTAirborne laser scanning (ALS) is commonly used in forest mapping. Full coverage of ALS is already available in some countries to provide high-detailed terrain elevation models. These kinds of data sets have been shown to offer great potential in forest mapping. However, it presents some drawbacks such as the resampling periods may be longer than recommended for forestry purposes or unexpected data updates. The recent development of digital photogrammetric algorithms makes dense image matching (DIM) point clouds an alternative to ALS in forest monitoring and management. Area-based approach estimations from ALS and DIM-based point clouds in a Pinus pinaster Ait. forest of Central Iberia were compared. Heights from image matching were normalized by an ALS-derived digital elevation model (DEM). A total of 50 sampling plots were used to fit non-parametric models for the estimation of forest structure variables. Plot-level validation revealed that DIM-based models predicted dominant height, stem number, basal area, and stem volume with root mean square error of 10.71%, 43.02%, 27.02%, and 26.80%, respectively. The corresponding results from ALS data were 11.06% for dominant height, 39.71% for stem number, 25.07% for basal area, and 25.60% for stem volume. This study demonstrates the usefulness of the combination of DIM with ALS-derived DEM to develop forest metrics and high-quality inventories in Mediterranean pine forests.

Archaeological prospection of forested areas using full-waveform airborne laser scanning

ALS point cloud filtering involves the separation of observations representing the physical terrain surface from those representing terrain details. A digital terrain model (DTM) is created from a subset of points representing the ground surface. The accuracy of the generated DTM is influenced by several factors, including the survey method used, the accuracy of the source data, the applied DTM generation algorithm, and the survey conditions.This article proposes the use of a new estimation method in the filtering of point clouds obtained from airborne laser scanning (ALS), provisionally called Msplit-estimation. The application of Msplit-estimation in ALS data filtering requires the determination of the appropriate functional model for the surface, which will be used in the filtering of the set of points. A polynomial terrain surface model was selected for this purpose. Two methods of filtering using the Msplit method are presented. The first is based on the estimated parameters of the polynomial describing the surface (called the ‘quality’ approach in the article). The second method (provisionally called the ‘quantity’ method) is carried out in two stages. The first stage is point cloud filtering, which results in two subsets being created. One of these is the subset of points intended for DTM creation, while the other contains the remaining points. The second stage of the approach is the creation of a DTM from the first subset.Since the Msplit method has an analytical character, the ATIN method was selected to verify the correct operation of the method. The ATIN method is based on computational geometry and uses repeated Delaunay triangulation and statistical evaluation of the geometric parameters. Comparison of Msplit with a method based on different principles mitigates errors arising from similarly functioning methods belonging to the same group of filters. The choice of the ATIN method was also dictated by its established position among filtering algorithms. The method is well-known, documented, and verified and this ensures that filtering by this method provides a reliable result that can serve as a reference for comparison with the proposed new filtering method.The theoretical discussion presented in this article was verified with two practical examples. The results obtained from computation by the Msplit method with appropriate terrain models encourage more detailed theoretical and empirical tests of this method for the filtering and segmentation of ALS data-sets.

The landslide, which occurred at Umyeon mountain (Mt. Umyeon) in Seoul, Korea in 2011, was a prime example that raised awareness about the landslide in the highly urbanized area. Although many studies have been done on Umyeon landslide, there is a lack of research that detects the area where the landslide occurred and quantifies the elevation changes through remote sensing data. In this regard, this paper aims to detect and assess topographic changes quantitatively over Mt. Umyeon by using digital elevation models (DEMs) derived from airborne laser scanning (ALS) data. Since Mt. Umyeon was hilly and covered with dense trees during summer, traces of the landslide were detected by estimating the spatially distributed uncertainty of ALS-derived DEMs. The probabilistic analysis with Bayes'™ theorem considering the spatially distributed DEM of difference (DoD) uncertainty enabled to detect the landslide traces efficiently and was less affected by the influence of ALS errors. The results indicated that ALS-derived DEMs have the potential to detect landslides with their uncertainty estimation, although the ALS data were acquired in hilly and densely vegetated areas. Moreover, quantifying topographic changes due to landslides with high reliability is considered to be beneficial and practically helpful for disaster recovery.

Airborne laser scanning (ALS) has demonstrated utility for forestry applications and has renewed interest in other forms of remotely sensed data, especially those that capture three-dimensional (3-D) forest characteristics. One such data source results from the advanced processing of high spatial resolution digital stereo imagery (DSI) to generate 3-D point clouds. From the derived point cloud, a digital surface model and forest vertical information with similarities to ALS can be generated. A key consideration is that when developing forestry related products such as a canopy height model (CHM), a high spatial resolution digital terrain model (DTM), typically from ALS, is required to normalize DSI elevations to heights above ground. In this paper we report on our investigations into the use of DSI-derived vertical information for capturing variations in forest structure and compare these results to those acquired using ALS. An ALS-derived DTM was used to provide the spatially detailed ground surface elevations to normalize DSI-derived heights. Similar metrics were calculated from the vertical information provided by both DSI and ALS. Comparisons revealed that ALS metrics provided a more detailed characterization of the canopy surface including canopy openings. Both DSI and ALS metrics had similar levels of correlation with forest structural attributes (e.g., height, volume, and biomass). DSI-based models predicted height, diameter, basal area, stem volume, and biomass with root mean square (RMS) accuracies of 11.2%, 21.7%, 23.6%, 24.5%, and 23.7%, respectively. The respective accuracies for the ALS-based predictions were 7.8%, 19.1%, 17.8%, 17.9%, and 17.5%. Change detection between ALS-derived CHM (time 1) and DSI-derived CHM (time 2) provided change estimates that demonstrated good agreement (r = 0.71) with two-date, ALS only, change outputs. For the single-layered, even-aged stands under investigation in this study, the DSI-derived vertical information is an appropriate and cost-effective data source for estimating and updating forest information. The accuracy of DSI information is based on a capability to measure the height of the upper canopy envelope with performance analogous to ALS. Forest attributes that are well captured and subsequently modeled from height metrics are best suited to estimation from DSI metrics, whereas ALS is more suitable for capturing stand density. Further investigation is required to better understand the performance of DSI-derived height products in more complex forest environments. Furthermore, the difference in variance captured between ALS and DSI-derived CHM also needs to be better understood in the context of change detection and inventory update considerations.

Dead wood and saproxylic beetles in set-aside and non set-aside forests in a boreal region

High-Resolution Topography (HRT) data sets are becoming increasingly available, improving our ability and opportunities to monitor geomorphic changes through multi-temporal Digital Terrain Models (DTMs). The use of repeated topographic surveys enables inferring the sediment dynamics of hazardous geomorphic processes such as floods, debris flows, and landslides, and allows us to derive important information on the risks often associated with these processes. The topographic surveying platforms, georeferencing systems, and processing tools have seen important developments in the last two decades, in particular Light Detection And Ranging (LiDAR) technology used in Airborne Laser Scanning (ALS) and Terrestrial Laser Scanning (TLS). Moreover, HRT data, produced through these techniques, changed a lot in terms of point cloud density, accuracy and precision over time. Therefore, old “legacy” data sets and recent surveys can often show comparison problems, especially when multi-temporal data are not homogeneous in terms of quality and uncertainties. In this context, data co-registration should be used to guarantee the coherence among multi-temporal surveys, minimizing, on stable areas, the distance between corresponding points acquired at different epochs. Although several studies highlight that this process is fundamental to properly compare multi-temporal DTMs, it is often not addressed in LiDAR post-processing workflows. In this paper we focus on the alignment of multi-temporal surveys in a topographically complex and rugged environment as the Moscardo debris-flow catchment (Eastern Italian Alps), testing various co-registration methods to align multi-temporal ALS point clouds (i.e. years 2003, 2009 and 2013) and the derived DTMs. In particular, we tested the pairwise registration with manual correspondences, the Iterative Closest Point (ICP) algorithm and a mathematical model that allows aligning simultaneously a generic number of point clouds, the so-called Generalized Procrustes Analysis (GPA), also in its GPA-ICP variant. Then, to correct the possible small inaccuracies generated from the gridding interpolation process, a custom-developed DTM co-registration tool (GRD-CoReg) was used to align gridded data. Both alignment phases (i.e. at point cloud and DTM level) proved to be fundamental and allowed us to obtain proper and reliable DTMs of Difference (DoDs), useful to quantify the debris mobilized and to detect the spatial and temporal patterns of catchment-scale erosion and deposition. The consistency of DoDs data was verified through the comparison between the erosion estimate of DoDs and the volumes of debris-flow events measured by the monitoring station close to the Moscardo torrent catchment outlet. The GPA-ICP algorithm followed by the GRD-CoReg tool proved to be the most effective solution for improving DoDs results with a decrease of systematic trend due to vertical and horizontal uncertainties between surveys, especially at steep slopes. The net volume difference (i.e. the sediment output from the catchment) of the 2003–2013 period changed from 3,237,896 m3 to 135,902 m3 in DoDs obtained from not co-registered and co-registered DTMs. The volume of debris flows measured at the catchment outlet during the same time interval amounts to 169,660 m3. The comparison with debris-flow volume measures at the monitoring station shows, therefore, that the DTMs obtained from the co-registration processes generate more reliable DoDs than those obtained from the raw DTMs (without the alignment).

Image matching as a data source for forest inventory – Comparison of Semi-Global Matching and Next-Generation Automatic Terrain Extraction algorithms in a typical managed boreal forest environment

Comparing the accuracies of forest attributes predicted from airborne laser scanning and digital aerial photogrammetry in operational forest inventories

The correct representation of the topography of terrain is an important requirement to generate photogrammetric products such as orthoimages and maps from high-resolution (HR) or very high-resolution (VHR) satellite datasets. The refining of the digital elevation model (DEM) for the generation of an orthoimage is a vital step with a direct effect on the final accuracy achieved in the orthoimages. The refined DEM has potential applications in various domains of earth sciences such as geomorphological analysis, flood inundation mapping, hydrological analysis, large-scale mapping in an urban environment, etc., impacting the resulting output accuracy. Manual editing is done in the presented study for the automatically generated DEM from IKONOS data consequent to the satellite triangulation with a root mean square error (RMSE) of 0.46, using the rational function model (RFM) and an optimal number of ground control points (GCPs). The RFM includes the rational polynomial coefficients (RPCs) to build the relation between image space and ground space. The automatically generated DEM initially represents the digital surface model (DSM), which is used to generate a digital terrain model (DTM) in this study for improving orthoimages for an area of approximately 100 km2. DSM frequently has errors due to mass points in hanging (floating) or digging, which need correction while generating DTM. The DTM assists in the removal of the geometric effects (errors) of ground relief present in the DEM (i.e., DSM here) while generating the orthoimages and thus improves the quality of orthoimages, especially in areas such as Dehradun that have highly undulating terrain with a large number of natural drainages. The difference image of reference, i.e., edited IKONOS DEM (now representing DTM) and automatically generated IKONOS DEM, i.e., DSM, has a mean difference of 1.421 m. The difference DEM (dDEM) for the reference IKONOS DEM and generated Cartosat-1 DEM at a 10 m posting interval (referred to as Carto10 DEM) results in a mean difference of 8.74 m.

Canopy height change (CHC) is one of the key characteristics of forest dynamics, associated with the fluctuations in forest above-ground biomass and carbon stocks. Field measurements and Airborne Laser Scanning (ALS) point clouds can be used to detect CHC; however, they have limited availability in space and time, making it challenging to map CHC over large areas. Alternatively, very high-resolution (VHR) satellite stereo imagery plays an increasingly vital role in estimating fine-scaled digital surface models (DSMs) across landscapes. However, its capability and potential to estimate canopy height model (CHM) and CHC has not been widely explored. Using ALS-derived CHM in 2011 and 2015 and the four-year CHC as references, we evaluated stereo-based CHM and CHC from WorldView-2 over five woody parks in Columbus, Ohio, USA. We also integrated stereo-based CHM with vegetation indices from Landsat 7 to improve CHM and CHC estimation with machine learning methods. Our results showed that VHR stereo imagery captured similar spatial patterns of CHM with ALS data but significantly overestimated CHM. Moreover, the ALS-derived CHC ranged from 1.6±1.9 m (mean ± standard deviation) to 3.1±1.2 m, as compared to from -1.3±0.7 m to 1.1±1.7 m for stereo-based CHC, indicating the limitation of CHC estimation by stereo imagery alone. Among six widely used machine learning methods, Gradient Boosting Regression method provided the most reliable estimates of CHM, with a correlation coefficient R of 0.64 and a root-mean-square error (RMSE) of 3.1 m (11.1%). Stereo-based CHM and vegetation indices explained more than 70% of CHM variability, substantially improving the estimation of 4-year CHC. Our results suggested that VHR stereo imagery alone has limitations in estimating CHM and CHC. The combination of remote-sensing structural (stereo-based CHM) and spectral (vegetation indices) information improves the CHM and CHC estimations.

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

This dissertation addresses the problem of automated vector extraction from Airborne Laser Scanner (ALS) data in urban areas. The recent popularity of Geographic Information Systems (GIS) has stimulated research on automated object extraction in order to simplify the data acquisition and update process. By automatically generating GIS inputs and updates from a single data source, the cost of data acquisition and processing will be kept to a minimum. Compared with other remote sensing techniques, extraction of objects from ALS data is in its infancy. ALS sensor technology has evolved rapidly and now allows the acquisition of very dense point clouds in a short period of time. ALS data is unique in that it explicitly contains 3D information and is acquired from an active sensor. As such, there are several benefits that can be immediately realised by using an ALS only data approach; data acquisition will not be limited to daylight hours as with other sensors and accurate height information is contained in the data. This means that registration of different data sources is not required and as only one data source is used, acquisition costs are minimised. Apart from these facts, ALS data has some other unique properties that have not been utilised to their full potential. An ALS sensor can record both height and intensity information from multiple returns in a single swath. To date, all the available information has rarely been used. For example, the intensity of a laser return has regularly been dismissed as it was considered under-sampled and noisy. This data still contains usable information and until this information is used, the optimum object recognition results will not be achieved; thus the use of as much of this information as possible is a major focus of this thesis. As ALS is an explicit 3D data source, the early stages of development were primarily focused on topographic mapping of terrain in forested areas in order to generate Digital Terrain Models (DTMs). As sensor technology has improved, so has the achievable resolution of point clouds from ALS data, and methods to extract objects from stand-alone ALS data have emerged. Attempts have been made to create city models and maps from ALS data that included buildings, roads, trees and powerlines. Each of these spatial object types has unique attributes, which means that automatic map creation is not an easy task. For example, buildings in general are easily detected in ALS data but the building outline is not easily delineated. Another difficulty with building detection is the separation of buildings and bridges as they have many of the same properties as observed by an ALS system. Bridges can usually be found in a road network but road extraction techniques typically produce poor detection rates and often require existing data and / or user interaction in semi-automatic techniques. These simple examples highlight the complexity of automatically generating vectorised maps from ALS data or in fact any data source. In this thesis, new methods are presented for the automatic creation of vectorised maps from ALS data. A two-step processing paradigm is adopted for this purpose, namely the classification of the ALS data and the vectorisation of the classification results. A classification strategy is introduced that creates a hierarchy for object detection with respect to the ALS data itself. This approach develops an ontology between the spatial object classes and the ALS data. The hierarchical framework highlights the fact that one object might not be discernable within the data without considering another. New classification algorithms are then presented within this framework. Each algorithm attempts to exploit the attributes of the data that are consistent within the spatial object class being considered as described by the classification framework. New algorithms for the classification of roads, trees and powerlines are all introduced whilst an extension to existing building classification methods is presented. Once classification is complete, a vectorisation process specific to the task at hand can be employed to yield vectorised results. These developed vectorisation processes are new and include an algorithm that has been generalised to allow the vectorisation of thick lines in images by detecting the centreline, direction and width. The primary goal of this thesis is to present a framework of new algorithms that will allow automatic spatial object detection and vectorisation whilst providing results of an acceptable quality. The algorithms presented rely solely on ALS data and require minimal operator knowledge. Each algorithm has been designed to exploit the way in which the object exhibits itself in the data. The new algorithms are integrated into a software package called JTD (Join The Dots) that will facilitate the effective automatic processing of ALS data. The results of the new algorithms have been evaluated over four 2 x 2 km areas that have been sampled with medium resolution ALS data. The results for each area are displayed and analysed to show the applicability of the whole process in an exemplary way.

Combining airborne and terrestrial laser scanning for quantifying erosion and deposition by a debris flow event