Search

Fine-resolution Light Detection and Ranging (LiDAR) data often exhibit excessive surface roughness that can hinder the characterization of topographic shape and the modeling of near-surface flow processes. Digital elevation model (DEM) smoothing methods, commonly low-pass filters, are sometimes applied to LiDAR data to subdue the roughness. These techniques can negatively impact the representation of topographic features, most notably drainage features, such as headwater streams. This paper presents the feature-preserving DEM smoothing (FPDEMS) method, which modifies surface normals to smooth the topographic surface in a similar manner to approaches that were originally designed for de-noising three-dimensional (3D) meshes. The FPDEMS method has been optimized for application with raster DEM data. The method was compared with several low-pass filters while using a 0.5-m resolution LiDAR DEM of an agricultural area in southwestern Ontario, Canada. The findings demonstrated that the technique was better at removing roughness, when compared with mean, median, and Gaussian filters, while also preserving sharp breaks-in-slope and retaining the topographic complexity at broader scales. Optimal smoothing occurred with kernel sizes of 11−21 grid cells, threshold angles of 10°−20°, and 3−15 elevation-update iterations. These parameter settings allowed for the effective reduction in roughness and DEM noise and the retention of terrace scarps, channel banks, gullies, and headwater streams.

Surface roughness is a terrain parameter that has been widely applied to the study of geomorphological processes. One of the main challenges in studying roughness is its highly scale-dependent nature. Determining appropriate mapping scales in topographically heterogenous landscapes can be difficult. A method is presented for estimating multiscale surface roughness based on the standard deviation of surface normals. This method utilizes scale partitioning and integral image processing to isolate scales of surface complexity. The computational efficiency of the method enables high scale sampling density and identification of maximum roughness for each grid cell in a digital elevation model (DEM). The approach was applied to a 0.5 m resolution LiDAR DEM of a 210 km2 area near Brantford, Canada. The case study demonstrated substantial heterogeneity in roughness properties. At shorter scales, tillage patterns and other micro-topography associated with ground beneath forest cover dominated roughness scale signatures. Extensive agricultural land-use resulted in 35.6% of the site exhibiting maximum roughness at micro-topographic scales. At larger spatial scales, rolling morainal topography and fluvial landforms, including incised channels and meander cut banks, were associated with maximum surface roughness. This method allowed for roughness mapping at spatial scales that are locally adapted to the topographic context of each individual grid cell within a DEM. Furthermore, the analysis revealed significant differences in roughness characteristics among soil texture categories, demonstrating the practical utility of locally adaptive, scale-optimized roughness.