A Comparative Study of Three Non-Geostatistical Methods for Optimising Digital Elevation Model Interpolation

It is common to generate digital elevation models (DEMs) from aerial laser scanning (ALS) data. However, cost and lack of knowledge may preclude its use. In contrast, global navigation satellite systems (GNSS) are seldom used to collect and generate DEMs. These receivers have the potential to be considered as data sources for DEM interpolation, as they can be inexpensive, easy to use, and mobile. The data interpolation method and spatial resolution from this method needs to be optimised to create accurate DEMs. Moreover, the density of GNSS data is likely to affect DEM accuracy. This study investigates three different deterministic approaches, in combination with spatial resolution and data thinning, to determine their combined effects on DEM accuracy. Digital elevation models were interpolated, with resolutions ranging from 0.5 m to 10 m using natural neighbour (NaN), topo to raster (ANUDEM), and inverse distance weighted (IDW) methods. The GNSS data were thinned by 25% (0.389 points m−2), 50% (0.259 points m−2), and 75% (0.129 points m−2) and resulting DEMs were contrast against a DEM interpolated from unthinned data (0.519 points m−2). Digital elevation model accuracy was measured by root mean square error (RMSE) and mean absolute error (MAE). It was found that the highest resolution, 0.5 m, produced the lowest errors in resulting DEMs (RMSE = 0.428 m, MAE = 0.274 m). The ANUDEM method yielded the greatest DEM accuracy from a quantitative perspective (RMSE = 0.305 m and MAE = 0.197 m); however, NaN produced a more visually appealing surface. In all the assessments, IDW showed the lowest accuracy. Thinning the input data by 25% and even 50% had relatively little impact on DEM quality; however, accuracy decreased markedly at 75% thinning (0.129 points m−2). This study showed that, in a time where ALS is commonly used to generate DEMs, GNSS-surveyed data can be used to create accurate DEMs. This study confirmed the need for optimization to choose the appropriate interpolation method and spatial resolution in order to produce a reliable DEM