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

Spectral filtering as a method of visualising and removing striped artefacts in digital elevation data

Spectral filtering was compared with traditional mean spatial filters to assess their ability to identify and remove striped artefacts in digital elevation data. The techniques were applied to two datasets: a 100 m contour derived digital elevation model (DEM) of southern Norway and a 2 m LiDAR DSM of the Lake District, UK. Both datasets contained diagonal data artefacts that were found to propagate into subsequent terrain analysis. Spectral filtering used fast Fourier transformation (FFT) frequency data to identify these data artefacts in both datasets. These were removed from the data by applying a cut filter, prior to the inverse transform. Spectral filtering showed considerable advantages over mean spatial filters, when both the absolute and spatial distribution of elevation changes made were examined. Elevation changes from the spectral filtering were restricted to frequencies removed by the cut filter, were small in magnitude and consequently avoided any global smoothing. Spectral filtering was found to avoid the smoothing of kernel based data editing, and provided a more informative measure of data artefacts present in the FFT frequency domain. Artefacts were found to be heterogeneous through the surfaces, a result of their strong correlations with spatially autocorrelated variables: landcover and landsurface geometry. Spectral filtering performed better on the 100 m DEM, where signal and artefact were clearly distinguishable in the frequency data. Spectrally filtered digital elevation datasets were found to provide a superior and more precise representation of the landsurface and be a more appropriate dataset for any subsequent geomorphological applications

High-resolution DEM generated from LiDAR data for water resource management

Terrain patterns play an important role in determining the nature of water resources and related hydrological modelling. Digital Elevation Models (DEMs), offering an efficient way to represent ground surface, allow automated direct extraction of hydrological features (Garbrecht and Martz, 1999), thus bringing advantages in terms of processing efficiency, cost effectiveness, and accuracy assessment, compared with traditional methods based on topographic maps, field surveys, or photographic interpretations. However, researchers have found that DEM quality and resolution affect the accuracy of any extracted hydrological features (Kenward et al., 2000). Therefore, DEM quality and resolution must be specified according to the nature and application of the hydrological features. The most commonly used DEM in Victoria, Australia is Vicmap Elevation delivered by the Land Victoria, Department of Sustainability and Environment. It was produced by using elevation data mainly derived from existing contour map at a scale of 1:25,000 and digital stereo capture, providing a state-wide terrain surface representation with a horizontal resolution of 20 metres. The claimed standard deviations, vertical and horizontal, are 5 metres and 10 metres respectively (Land- Victoria, 2002). In worst case, horizontal errors could be up to ±30m. Although high resolution stereo aerial photos provide a potential way to generate high resolution DEMs, under the limitations of currently used technologies by prevalent commercial photogrammetry software, only DSMs (Digital Surface Models) other than DEMs can be directly generated. Manual removal of the nonground data so that the DSM is transformed into a DEM is time consuming. Therefore, using stereo aerial photos to produce DEM with currently available techniques is not an accurate and costeffective method. Light Detection and Ranging (LiDAR) data covering 6900 km² of the Corangamite Catchment area of Victoria were collected over the period 19 July 2003 to 10 August 2003. It will be used to support a series of salinity and water management projects for the Corangamite Catchment Management Authority (CCMA). The DEM derived from the LiDAR data has a vertical accuracy of 0.5m and a horizontal accuracy of 1.5m. The high quality DEM leads to derive much detailed terrain and hydrological attributes with high accuracy. Available data sources of DEMs in a catchment management area were evaluated in this study, including the Vicmap DEM, a DEM generated from stereo aerial photos, and LiDAR-derived DEM. LiDAR technology and LiDAR derived DEM were described. In order to assess the capability of LiDAR-derived DEM for improving the quality of extracted hydrological features, sub-catchment boundaries and drainage networks were generated from the Vicmap DEM and the LiDAR-derived DEM. Results were compared and analysed in terms of accuracy and resolution of DEMs. Elevation differences between Vicmap and LiDAR-derived DEMs are significant, up to 65m in some areas. Subcatchment boundaries derived from these two DEMs are also quite different. In spite of using same resolution for the Vicmap DEM and the LiDARderived DEM, high accuracy LiDAR-derived DEM gave a detailed delineation of sub-catchment. Compared with results derived from LiDAR DEM, the drainage networks derived from Vicmap DEM do not give a detailed description, and even lead to discrepancies in some areas. It is demonstrated that a LiDAR-derived DEM with high accuracy and high resolution offers the capability of improving the quality of hydrological features extracted from DEMs

Techniques for augmenting the visualisation of dynamic raster surfaces

Despite their aesthetic appeal and condensed nature, dynamic raster surface representations such as a temporal series of a landform and an attribute series of a socio-economic attribute of an area, are often criticised for the lack of an effective information delivery and interactivity.In this work, we readdress some of the earlier raised reasons for these limitations -information-laden quality of surface datasets, lack of spatial and temporal continuity in the original data, and a limited scope for a real-time interactivity. We demonstrate with examples that the use of four techniques namely the re-expression of the surfaces as a framework of morphometric features, spatial generalisation, morphing, graphic lag and brushing can augment the visualisation of dynamic raster surfaces in temporal and attribute series

Developing Guidelines for Two-Dimensional Model Review and Acceptance

Two independent modelers ran two hydraulic models, SRH-2D and HEC-RAS 2D. The models were applied to the Lakina River (MP 44 McCarthy Road) and to Quartz Creek (MP 0.7 Quartz Creek Road), which approximately represent straight and bend flow conditions, respectively. We compared the results, including water depth, depth averaged velocity, and bed shear stress, from the two models for both modelers. We found that the extent and density of survey data were insufficient for Quartz Creek. Neither model was calibrated due to the lack of basic field data (i.e., discharge, water surface elevation, and sediment characteristics). Consequently, we were unable to draw any conclusion about the accuracy of the models. Concerning the time step and the equations used (simplified or full) to solve the momentum equation in the HEC-RAS 2D model, we found that the minimum time step allowed by the model must be used if the diffusion wave equation is used in the simulations. A greater time step can be used if the full momentum equation is used in the simulations. We developed a set of guidelines for reviewing model results, and developed and provided a two-day training workshop on the two models for ADOT&PF hydraulic engineers

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

Accuracy assessment in glacier change analysis

This thesis assesses the accuracy of digital elevation models (DEM) generated from contour lines and LiDAR points (Light Detection and Ranging) employing several interpolation methods at different resolutions. The study area is Jostefonn glacier that is situated in Sogn og Fjordane county, Norway. There are several ways to assess accuracy of DEMs including simple ways such as visual comparison and more sophisticated methods like relative and absolute comparison. Digital elevation models of the Jostefonn glacier were created from contour lines for years 1966 and 1993. LiDAR data from year 2011 was used as a reference data set. Of all the interpolation methods tested Natural Neighbours (NN) and Triangular Irregular Network (TIN) algorithms rendered the best results and proved to be superior to other interpolation methods. Several resolutions were tested (the cell size of 5 m, 10 m, 20 m and 50 m) and the best outcome was achieved by as small cell size as possible. The digital elevation models were compared to a reference data set outside the glacier area both on a cell-by-cell basis and extracting information at test points. Both methods rendered the same results that are presented in this thesis. Several techniques were employed to assess the accuracy of digital elevation models including visualization and statistical analysis. Visualization techniques included comparison of the original contour lines with those generated from DEMs. Root mean square error, mean absolute error and other accuracy measures were statistically analysed. The greatest elevation difference between the digital elevation model of interest and the reference data set was observed in the areas of a steep terrain. The steeper the terrain, the greater the observed error. The magnitude of the errors can be reduced by using a smaller cell size but that this is offset by a larger amount of data and increased data processing time.Popular science Glaciers are very sensitive indicators of climate change. The major cause of melting glaciers is global warming. This rapid rate of melting has serious negative impact on the earth causing flooding, leaving impact on flora and fauna, resulting in shortage of freshwater and hydroelectricity. The long-term monitoring of glaciers and the knowledge gained from it can help governments, environmental and water resource managers to make plans to cope with impacts of climate change. Results from glacier monitoring ought to be precise, showing the actual situation compared to the situation in the past as well as predicting possible glacier changes in the future. The aim of this thesis was to investigate how sensitive the results were to different methods used in glacier change detection focusing on the quality of Digital Elevation Models (DEMs). The study area of this thesis was the Jostefonn glacier situated in Sogn and Fjordane county, Norway. Digital elevation models were created from contour lines for years 1966 and 1993. LiDAR data from year 2011 was used as a reference data set. Several techniques were employed to estimate the accuracy of digital elevation models including visualization, statistical analysis, analysing the accuracy of digital elevation models for terrain on different slopes, comparison to a reference data set outside the glacier area that was considered to be stable and where no elevation change was expected. The original contour lines (1966 and 1993) were compared with the ones generated from the created terrain models (glacier area) as well as with the contour lines from the reference data set (outside the glacier area) by visualization techniques. Accuracy measures (Root Mean Square Error, Mean Absolute Error and others) were statistically analysed. Natural Neighbours and Triangular Irregular Network interpolators proved to be superior to other algorithms used to create the terrain models. The best outcome was achieved by using as small cell size as possible. 5 m resolution rendered the best results from the resolutions tested (5 m, 10 m, 20 m and 50 m). The greatest elevation differences were observed in the areas of a steep terrain. The steeper the terrain, the greater the elevation difference. The terracing effect was noticed in the digital elevation models due to the high density of elevation points on the contour lines and hardly any points between them. Useful information can be obtained by estimating accuracy of digital elevation models. The accuracy of terrain models determines the reliability of glacier change analysis and that is why the digital elevation model must represent the terrain as accurately as possible. The different methods used in this thesis rendered very similar results and that indicated that the results were reliable and the terrain models created with Natural Neighbours and Triangular Irregular Network interpolators (resolution of 5 m) can be employed in further glacier change analysis