Airborne LiDAR for DEM generation: some critical issues

Airborne LiDAR is one of the most effective and reliable means of terrain data collection. Using LiDAR data for DEM generation is becoming a standard practice in spatial related areas. However, the effective processing of the raw LiDAR data and the generation of an efficient and high-quality DEM remain big challenges. This paper reviews the recent advances of airborne LiDAR systems and the use of LiDAR data for DEM generation, with special focus on LiDAR data filters, interpolation methods, DEM resolution, and LiDAR data reduction. Separating LiDAR points into ground and non-ground is the most critical and difficult step for DEM generation from LiDAR data. Commonly used and most recently developed LiDAR filtering methods are presented. Interpolation methods and choices of suitable interpolator and DEM resolution for LiDAR DEM generation are discussed in detail. In order to reduce the data redundancy and increase the efficiency in terms of storage and manipulation, LiDAR data reduction is required in the process of DEM generation. Feature specific elements such as breaklines contribute significantly to DEM quality. Therefore, data reduction should be conducted in such a way that critical elements are kept while less important elements are removed. Given the highdensity characteristic of LiDAR data, breaklines can be directly extracted from LiDAR data. Extraction of breaklines and integration of the breaklines into DEM generation are presented

Modelling the spatial distribution of DEM Error

Assessment of a DEM’s quality is usually undertaken by deriving a measure of DEM accuracy – how close the DEM’s elevation values are to the true elevation. Measures such as Root Mean Squared Error and standard deviation of the error are frequently used. These measures summarise elevation errors in a DEM as a single value. A more detailed description of DEM accuracy would allow better understanding of DEM quality and the consequent uncertainty associated with using DEMs in analytical applications. The research presented addresses the limitations of using a single root mean squared error (RMSE) value to represent the uncertainty associated with a DEM by developing a new technique for creating a spatially distributed model of DEM quality – an accuracy surface. The technique is based on the hypothesis that the distribution and scale of elevation error within a DEM are at least partly related to morphometric characteristics of the terrain. The technique involves generating a set of terrain parameters to characterise terrain morphometry and developing regression models to define the relationship between DEM error and morphometric character. The regression models form the basis for creating standard deviation surfaces to represent DEM accuracy. The hypothesis is shown to be true and reliable accuracy surfaces are successfully created. These accuracy surfaces provide more detailed information about DEM accuracy than a single global estimate of RMSE

Evaluating the Differences of Gridding Techniques for Digital Elevation Models Generation and Their Influence on the Modeling of Stony Debris Flows Routing: A Case Study From Rovina di Cancia Basin (North-Eastern Italian Alps)

Debris \ufb02ows are among the most hazardous phenomena in mountain areas. To cope with debris \ufb02ow hazard, it is common to delineate the risk-prone areas through routing models. The most important input to debris \ufb02ow routing models are the topographic data, usually in the form of Digital Elevation Models (DEMs). The quality of DEMs depends on the accuracy, density, and spatial distribution of the sampled points; on the characteristics of the surface; and on the applied gridding methodology. Therefore, the choice of the interpolation method affects the realistic representation of the channel and fan morphology, and thus potentially the debris \ufb02ow routing modeling outcomes. In this paper, we initially investigate the performance of common interpolation methods (i.e., linear triangulation, natural neighbor, nearest neighbor, Inverse Distance to a Power, ANUDEM, Radial Basis Functions, and ordinary kriging) in building DEMs with the complex topography of a debris \ufb02ow channel located in the Venetian Dolomites (North-eastern Italian Alps), by using small footprint full- waveform Light Detection And Ranging (LiDAR) data. The investigation is carried out through a combination of statistical analysis of vertical accuracy, algorithm robustness, and spatial clustering of vertical errors, and multi-criteria shape reliability assessment. After that, we examine the in\ufb02uence of the tested interpolation algorithms on the performance of a Geographic Information System (GIS)-based cell model for simulating stony debris \ufb02ows routing. In detail, we investigate both the correlation between the DEMs heights uncertainty resulting from the gridding procedure and that on the corresponding simulated erosion/deposition depths, both the effect of interpolation algorithms on simulated areas, erosion and deposition volumes, solid-liquid discharges, and channel morphology after the event. The comparison among the tested interpolation methods highlights that the ANUDEM and ordinary kriging algorithms are not suitable for building DEMs with complex topography. Conversely, the linear triangulation, the natural neighbor algorithm, and the thin-plate spline plus tension and completely regularized spline functions ensure the best trade-off among accuracy and shape reliability. Anyway, the evaluation of the effects of gridding techniques on debris \ufb02ow routing modeling reveals that the choice of the interpolation algorithm does not signi\ufb01cantly affect the model outcomes

The effects of data reduction on LiDAR-based digital elevation models

LiDAR data enables highly accurate terrain representations, however, various applications are hampered by data handling efficiency; specifically lengthy processing times. To address this, both point density reductions and the use of various resolution grids are compared as data reduction methods to test their effects on the accuracy and handling efficiency of the derived Digital Elevation Model (DEM). A series of point densities of 1%, 10%, 25%, 50% and 75% were interpolated along a range of horizontal resolutions (1-, 2-, 3-, 4-, 5-, 10-, and 30- m). Results indicate that resolution reduction provides the most efficient DEMs in terms of their data handling. DEMs generated at a 3 m resolution using all of the data points deviated less than 6% from the 1mDEM100%, while significantly only taking 10% of the processing time. Resolution reduction provided sufficient accuracies for varying terrain complexities

Sediment Dynamics and Channel Connectivity on Hillslopes

The pattern, magnitude, and frequency of hillslope erosion and deposition are spatially varied under the influence of micro-topography and channel geometry. This research investigates the interrelationships between erosion/deposition, micro-topography, and channel connectivity on a hillslope in Loudon, Tennessee using the centimeter (cm) level temporal Digital Elevation Models collected using laser scanning. This research addressed (1) the effect of spatial resolution on the erosion/deposition quantification, and rill delineation; (2) the influences of micro-topographic factors (e.g. slope, roughness, aspect) on erosion and deposition; (3) the relationship between the structural connectivity -- depressions and confluence of rills -- and the sedimentological connectivity. I conducted (1) visual and quantitative assessments for the erosion and deposition, and the revised automated proximity and conformity analysis for the rill network; (2) quantile regression for micro-topographic factors using segmented rill basins; and (3) cross-correlation analysis using erosion and deposition series along the channels.Overall, rills are sedimentologically more dynamic than the interrill areas. A larger grid size reduces the detectable changes in both areal and volumetric quantities, and also decreases the total length and number of rills. The offset between delineated rills and the reference increases with larger grid sizes. A larger rill basin has higher erosion and deposition with the magnitude of erosion greater than deposition. The slope has a positive influence on erosion and a negative one on deposition; roughness has a positive influence on deposition and a negative one on erosion. Areas that are more north-facing experience higher erosion and lower deposition. Rill length explains 46% of the variability for erosion and 24% for deposition. The depressions are associated with higher erosion in the downslope direction. The correlations between the erosion and the confluence are positive; the correlation between the deposition and the sink is positive. Overall, the influence of structural connectivity on the sedimentological connectivity is within 25 cm in both upstream and downstream directions. This research contributes to the understanding in how the sediment movement on hillslopes is governed by topographic variations and channel connectivity, and future work may explore hillslope channels at broader geographical and temporal scales

Benchmarking 2D hydraulic models for urban flood simulations

This paper describes benchmark testing of six two-dimensional (2D) hydraulic models (DIVAST, DIVASTTVD, TUFLOW, JFLOW, TRENT and LISFLOOD-FP) in terms of their ability to simulate surface flows in a densely urbanised area. The models are applied to a 1·0 km × 0·4 km urban catchment within the city of Glasgow, Scotland, UK, and are used to simulate a flood event that occurred at this site on 30 July 2002. An identical numerical grid describing the underlying topography is constructed for each model, using a combination of airborne laser altimetry (LiDAR) fused with digital map data, and used to run a benchmark simulation. Two numerical experiments were then conducted to test the response of each model to topographic error and uncertainty over friction parameterisation. While all the models tested produce plausible results, subtle differences between particular groups of codes give considerable insight into both the practice and science of urban hydraulic modelling. In particular, the results show that the terrain data available from modern LiDAR systems are sufficiently accurate and resolved for simulating urban flows, but such data need to be fused with digital map data of building topology and land use to gain maximum benefit from the information contained therein. When such terrain data are available, uncertainty in friction parameters becomes a more dominant factor than topographic error for typical problems. The simulations also show that flows in urban environments are characterised by numerous transitions to supercritical flow and numerical shocks. However, the effects of these are localised and they do not appear to affect overall wave propagation. In contrast, inertia terms are shown to be important in this particular case, but the specific characteristics of the test site may mean that this does not hold more generally

Defining optimal DEM resolutions and point densities for modelling hydrologically sensitive areas in agricultural catchments dominated by microtopography

AbstractDefining critical source areas (CSAs) of diffuse pollution in agricultural catchments depends upon the accurate delineation of hydrologically sensitive areas (HSAs) at highest risk of generating surface runoff pathways. In topographically complex landscapes, this delineation is constrained by digital elevation model (DEM) resolution and the influence of microtopographic features. To address this, optimal DEM resolutions and point densities for spatially modelling HSAs were investigated, for onward use in delineating CSAs. The surface runoff framework was modelled using the Topographic Wetness Index (TWI) and maps were derived from 0.25m LiDAR DEMs (40 bare-earth points m−2), resampled 1m and 2m LiDAR DEMs, and a radar generated 5m DEM. Furthermore, the resampled 1m and 2m LiDAR DEMs were regenerated with reduced bare-earth point densities (5, 2, 1, 0.5, 0.25 and 0.125 points m−2) to analyse effects on elevation accuracy and important microtopographic features. Results were compared to surface runoff field observations in two 10km2 agricultural catchments for evaluation. Analysis showed that the accuracy of modelled HSAs using different thresholds (5%, 10% and 15% of the catchment area with the highest TWI values) was much higher using LiDAR data compared to the 5m DEM (70–100% and 10–84%, respectively). This was attributed to the DEM capturing microtopographic features such as hedgerow banks, roads, tramlines and open agricultural drains, which acted as topographic barriers or channels that diverted runoff away from the hillslope scale flow direction. Furthermore, the identification of ‘breakthrough’ and ‘delivery’ points along runoff pathways where runoff and mobilised pollutants could be potentially transported between fields or delivered to the drainage channel network was much higher using LiDAR data compared to the 5m DEM (75–100% and 0–100%, respectively). Optimal DEM resolutions of 1–2m were identified for modelling HSAs, which balanced the need for microtopographic detail as well as surface generalisations required to model the natural hillslope scale movement of flow. Little loss of vertical accuracy was observed in 1–2m LiDAR DEMs with reduced bare-earth point densities of 2–5 points m−2, even at hedgerows. Further improvements in HSA models could be achieved if soil hydrological properties and the effects of flow sinks (filtered out in TWI models) on hydrological connectivity are also considered