D-SRCAGAN : DEM Super-resolution Generative Adversarial Network

The scarcity of high-resolution urban digital elevation model (DEM) datasets, particularly in certain developing countries, has posed a challenge for many water-related applications such as flood risk management. A solution to address this is to develop effective approaches to reconstruct high-resolution DEMs from their low-resolution equivalents that are more widely available. However, the current high-resolution DEM reconstruction approaches mainly focus on natural topography. Few attempts have been made for urban topography, which is typically an integration of complex artificial and natural features. This study proposed a novel multi-scale mapping approach based on convolutional neural network (CNN) to deal with the complex features of urban topography and to reconstruct high-resolution urban DEMs. The proposed multi-scale CNN model was firstly trained using urban DEMs that contained topographic features at different resolutions, and then used to reconstruct the urban DEM at a specified (high) resolution from a low-resolution equivalent. A two-level accuracy assessment approach was also designed to evaluate the performance of the proposed urban DEM reconstruction method, in terms of numerical accuracy and morphological accuracy. The proposed DEM reconstruction approach was applied to a 121 km2 urbanized area in London, United Kingdom. Compared with other commonly used methods, the current CNN-based approach produced superior results, providing a cost-effective innovative method to acquire high-resolution DEMs in other data-scarce regions.

A rapid high-resolution multi-sensory urban flood mapping framework via DEM upscaling

High resolution digital elevation models (DEM) are becoming increasingly available for use as source data during the process of creating synthetic environments to support simulation systems. Several data sets that provide elevation points corresponding to 1 meter intervals on the earth surface are now available. However, multiple transformations must often be applied to the raw source data before it is suitable for use by any simulation system. These transformations have an impact on the fidelity of the final (simulation) synthetic environment that is difficult to quantify. Further, intuition alone now supports any claim that higher resolution source data necessarily results in generation of higher fidelity simulation data as a product of the transformation process. This paper documents an attempt to measure fidelity differences in final simulation synthetic environments that can be directly attributed to the resolution of the source data. Specifically, several lower resolution DEM are generated from a single high resolution (1 meter horizontal spacing) source DEM and all are used as source data for TIN construction. Automated planning software is applied to each and used as a metric to measure TIN quality.

Mid-ocean ridges are unique, tectonically active geographical units on Earth that profoundly control the ocean environment and dynamics at the global scale. However, high-resolution topographic data from mid-ocean ridges are rarely available due to the difficulty in detecting ocean floors, which further limits ocean research at the global scale. Here, we divide the global mid-ocean ridge system into 2805 tiles and reconstruct their high-resolution topography by using a transfer learning approach with freely available low-resolution digital elevation models (DEMs) and limited high-resolution DEMs. A high-frequency terrain feature-based deep residual network is proposed to generate high-resolution global mid-ocean ridge DEMs. In this network, topographic knowledge related to mid-ocean ridges is integrated and quantified to improve the learning efficiency and reconstruction quality of the network. A series of verifications and evaluations demonstrate the reliability of reconstructed topographies for submarine topography research. We observe that reconstructed topography can achieve good environmental understanding and information acquisition in the global mid-ocean ridge range. We find that the complexity of the previous terrain environment is underestimated by 26.63% in terms of the slope gradient and by 14.95% in terms of terrain relief, while a 101.10% information improvement can be obtained for the reconstructed topography. The reconstructed topography indicates that diverse and intricate topographical environments of mid-ocean ridges exist among different ocean regions. The proposed transfer learning method for reconstructing high-resolution mid-ocean ridge topographies is valuable and can be utilized for reconstructing information in regions that are difficult to observe directly and lack sufficient data.

Precisely measuring the Earth’s changing surface on decadal to centennial time scales is critical for many science and engineering applications, yet long-term records of quantitative landscape change are often temporally and geographically sparse. Archives of scanned historical aerial photographs provide an opportunity to augment these records with accurate elevation measurements that capture the historical state of the Earth surface. Structure from Motion (SfM) photogrammetry workflows produce high-quality digital elevation models (DEMs) and orthoimage mosaics from these historical images, but time-intensive tasks like manual image preprocessing (e.g., fiducial marker identification) and ground control point (GCP) selection impede processing at scale. We developed an automated method to process historical images and generate self-consistent time series of high-resolution (0.5–2 m) DEMs and orthomosaics, without manual GCP selection. The method relies on SfM to correct camera interior and exterior orientation and a robust multi-stage co-registration approach using modern reference terrain datasets for geolocation refinement. We demonstrate the method using scanned images from the North American Glacier Aerial Photography (NAGAP) archive collected between 1967 and 1997. We present results for two sites with variable photo acquisition geometry and overlap — Mount Baker and South Cascade Glacier in Washington State, USA. The automated method corrects initial camera position errors of several kilometers and produces accurately georeferenced, high-resolution DEMs and orthoimages, regardless of camera configuration, acquisition geometry, terrain characteristics, and reference DEM properties. The average RMS reprojection error following bundle adjustment optimization was 0.67 px (0.15 m) for the 261 images contributing to 10 final DEM mosaics between 1970 and 1992 at Mount Baker, and 0.65 px (0.13 m) for the 243 images contributing to 18 individual DEMs between 1967 and 1997 at South Cascade Glacier. The relative accuracy of elevation values in the historical time series stacks was 0.68 m at Mount Baker and 0.37 m at South Cascade Glacier. Our products have reduced systematic error and improved accuracy compared to DEM products generated using SfM with manual GCP selection. Final elevation change measurement precision was ∼0.7–1.0 m over a 30-year period, enabling the study of processes with rates as low as ∼1-3 cm/yr. Our results demonstrate the potential of this scalable method to rapidly process archives of historical imagery and deliver new quantitative insights on long-term geodetic change and Earth surface processes.

Utilizing multispectral satellite data and digital elevation models (DEMs) has emerged as the primary approach for cartographically representing landforms. By using high-resolution satellite photos that capture spatial, temporal, spectral, and radiometric data, one may get a fresh comprehension of the geomorphology of a particular area by recognizing its landforms. In addition, a synergistic method is used by using data produced from digital elevation models (DEMs) such as Slope, Aspect, Hillshade, Curvature, Contour Patterns, and 3-D Flythrough Visuals. The increasing use of UAV (drone) technology for obtaining high-resolution digital images and elevation models has become an essential element in developing complete topographic models in landslide scars that are very unstable and prone to erosion. Comparison (differences in values) of seven (7) different DEMs between two algorithms used, i.e., QGIS and White Box Tool (WBT), were successfully attempted in the present research. The TLS, UAV and Satellite data of the study area—Kshetrapal Landslide, Chamoli (District), Uttarakhand (State), India was subjected to two different algorithms (QGIS and WBT) to evaluate and differentiate seven different DEMs (CARTOSAT, ASTER, SRTM, Alos 3D, TanDEM, MERIT, and FabDEM/FATHOM) taking into consideration various parameters viz. Aspect, Hillshade, Slope, Mean Curvature, Plan Curvature, Profile Curvature and Total Curvature. The different values of aforesaid parameters of various DEMs evaluated (using algorithms QIGS and WBT) reveal that only three parameters, i.e., Aspect, Hillshade, and Slope, show results. In contrast, the remaining ones do not show any meaningful results, and therefore, the comparison was possible only with regard to these three parameters. The comparison is drawn by comparing minimum, maximum, and elevation values (by subtracting WBT values from QGIS values) regarding Aspect, Hillshade, and Slope, arranging the differences in values as per their importance. (Increasing or decreasing order), assigning merit scores individually, and then cumulatively, and ascertaining the order of application suitability of various Dems, which stand in the order of (CARTOSAT, ASTER, SRTM, Alos 3D, TanDEM, and MERIT, and FabDEM/FATHOM).

Quantifying landscape form can provide crucial insight into the interactions between tectonics and climate. River long profile morphology, quantified by metrics such as channel steepness, is the most commonly used tool to investigate topographic form, with many studies relating long profile morphology to uplift rate, precipitation, sediment properties, or lithology, for example. River long profiles record the signal of external forcing over large spatial scales (i.e. tens of kilometres). This has many advantages: for example, it is a convenient scale for analysing variations in large-scale processes, such gradients in tectonic uplift. It also means that high resolution digital elevation models (DEMs) are not required and therefore river long profiles can be extracted globally. However, analysis of river long profiles over tens of kilometres can also result in signal smoothing and subsequent loss of finer scale tectonic or climatic signatures encoded into the landscape.Tectonic and climatic processes do not only leave their fingerprint in the long profiles of rivers. Hilltops, hillslopes, and valleys make up the majority of Earth’s landscapes by area, yet their morphology has received much less attention than that of rivers. This is in part due to the difficulty in accurately extracting hilltops and valley morphology from DEMs, especially on a global scale. Here, I show that we can now extract hilltop and valley metrics from high-resolution (< 15 m) DEMs over orogenic to continental scales using new topographic analysis techniques and high-performance computing facilities. I argue that by combining hilltop, hillslope, and valley metrics, we can obtain more information about tectonic and climatic processes than from river profiles alone.

Abstract. Digital Elevation Model (DEM) is an important input for geo-spatial analysis. For various applications like flood management, ortho rectification of remote sensing images, navigation, architectural works, defence, etc., high resolution DEM is required. TanDEM-X mission was launched in 2010 to obtain high resolution global DEM with HTRI-3 standard. SAR interferometry (InSAR) technique is used for DEM generation from TanDEM-X SAR data. The accuracy of DEM depends on many parameters like height ambiguity, incidence angle, polarization, etc. In this study, time series TanDEM-X data spanning over 3 years, had processed for generating DEM at the spatial resolution of 6 m and their accuracy had studied using DGPS elevation data and SRTM 90 m DEM. The products generated during DEM generation process are DEM, precision (or height error), coherence, layover and shadow images. Using weighted average fusion technique, ascending and descending DEMs are fused for improving the quality of DEM and to reduce invalid pixels corresponding to layover and shadow areas. Results from time series data were analysed and found RMSE error of fused DEMs is in the range of 2 m to 4 m, while individual DEM has accuracy of 3 m to 6 m with respect to DGPS elevation data. Fused DEMs are having high accuracy as well as less voids. The reduction of voids by fusion, ranges from 40 to 85 % in different combinations of data.

Satellite-derived Digital Elevation Model (DEM) selection, preparation and correction for hydrodynamic modelling in large, low-gradient and data-sparse catchments

To accurately assess landslide susceptibility and model debris flow paths, acquisition of high-resolution elevation data is essential. These data enable precise topographical analysis, considering factors like slope and curvature. But most of the time such high-resolution data will have inaccuracies due to vegetation, especially in tropical region. This study proposes a comprehensive approach that integrates field data and Structure from Motion (SfM) technology to rectify such DEM inaccuracies caused by dense vegetation. The goal is to enhance the accuracy of landslide simulations and volumetric analysis for effective post-disaster management. The Western Ghats, running parallel to India's western coast, has experienced a surge in rainfall-induced landslides, resulting in significant casualties in recent years. Notably, the Pettimudi village landslide in 2020 caused 70 deaths. The unique geomorphological features of the Western Ghats, such as concave curvature, colluvium deposits, and evidence of paleo landslides, contribute to the area's susceptibility. The study emphasizes the need for detailed assessment and mitigation strategies. The proposed method aims to improve the post-event high-resolution DEM accuracy by integrating field-collected elevation values and utilizing Agisoft Metashape software with the SfM algorithm. The rectification process involves combining elevation differences measured during a field study with photogrammetric elevation data. The field-collected elevation differences are crucial for rectifying these points, and enhancing the Digital Elevation Model (DEM) accuracy. The inaccessible source region is improved using a SfM-created DEM from drone footage, resulting in a more accurate post-event DEM. Correcting the 1 m DEM using an SfM-generated DEM proves challenging but significantly improves detail and accuracy, especially over the landslide source area. The impact of the rectification on accuracy is assessed by comparing volumes calculated from the initial DEM and the newly corrected DEM. The difference in volumes of debris depletion and accumulation, computed using initial and corrected DEMs, highlights variations, with depleted volumes being significantly larger due to the extensive depth increase over the landslide source region. These volumes are then utilized in Rapid Mass Movement Simulation (RAMMS) to validate the rectification process and enhance landslide impact predictions. Debris flow simulations in RAMMS, utilizing the rectified 1 m DEM, show specific outcomes at the landslide source region and toe of the landslide. The study emphasizes the importance of integrating SfM technology with field data to improve DEM accuracy, acknowledging the significance of additional field data for further refinement. The potential adoption of highly precise aerial imagery from Unmanned Aerial Vehicle (UAV) surveys is suggested to further enhance the SfM DEM. Debris flow modeling with RAMMS serves as a vital step in validating the accuracy and reliability of the rectified elevation model. Keywords: Landslides, Western Ghats, Pettimudi, Digital elevation model (DEM), Structure-from-Motion (SfM) technology, Rapid Mass Movement Simulation (RAMMS)

Voids caused by shadow, layover, and decorrelation usually occur in digital elevation models (DEMs) of mountainous areas that are derived from interferometric synthetic aperture radar (InSAR) datasets. The presence of voids degrades the quality and usability of the DEMs. Thus, void removal is considered as an integral part of the DEM production using InSAR data. The fusion of multiple DEMs has been widely recognized as a promising way for the void removal. Because the vertical accuracy of multiple DEMs can be different, the selection of optimum weights becomes a key problem in the fusion and is studied in this article. As a showcase, two high-resolution InSAR DEMs near Mt. Qilian in northwest China are created and then merged. The two pairs of InSAR data were acquired by TerraSAR-X from an ascending orbit and COSMO-SkyMed from a descending orbit. A maximum likelihood fusion scheme with the weights optimally determined by the height of ambiguity and the variance of phase noise is adopted to syncretize the two DEMs in our study. The fused DEM has a fine spatial resolution of 10 m and depicts the landform of the study area well. The percentage of void cells in the fused DEM is only 0.13 %, while 6.9 and 5.7 % of the cells in the COSMO-SkyMed DEM and the TerraSAR-X DEM are originally voids. Using the ICESat/GLAS elevation data and the Chinese national DEM of scale 1:50,000 as references, we evaluate vertical accuracy levels of the fused DEM as well as the original InSAR DEMs. The results show that substantial improvements could be achieved by DEM fusion after atmospheric phase screen removal. The quality of fused DEM can even meet the high-resolution terrain information (HRTI) standard.

Snow depth from ICESat laser altimetry — A test study in southern Norway

Using low resolution digital elevation model (DEM) data, reconstruction of cartographic contours is focused on. With a small number of control points, we aim to draw a smooth curve as well as to compute the area and perimeter of a closed contour accurately. The input data forms an oriented contour, each control point of which consists of three values: x, y, and tangent angle. An interpolation method based on arc spline is developed. The method can reproduce G1 continuous curves,in which the exact area and perimeter can be calculated. In a numerical experiment for analytically defined curves, accurate computation of the area and perimeter was achieved with a smaller number of control points. In another experiment using a DEM data, the reconstructed contours were smoother than those by a conventional method.

Abstract. Lunar topographic information, e.g., lunar DEM (Digital Elevation Model), is very important for lunar exploration missions and scientific research. Lunar DEMs are typically generated from photogrammetric image processing or laser altimetry, of which photogrammetric methods require multiple stereo images of an area. DEMs generated from these methods are usually achieved by various interpolation techniques, leading to interpolation artifacts in the resulting DEM. On the other hand, photometric shape reconstruction, e.g., SfS (Shape from Shading), extensively studied in the field of Computer Vision has been introduced to pixel-level resolution DEM refinement. SfS methods have the ability to reconstruct pixel-wise terrain details that explain a given image of the terrain. If the terrain and its corresponding pixel-wise albedo were to be estimated simultaneously, this is a SAfS (Shape and Albedo from Shading) problem and it will be under-determined without additional information. Previous works show strong statistical regularities in albedo of natural objects, and this is even more logically valid in the case of lunar surface due to its lower surface albedo complexity than the Earth. In this paper we suggest a method that refines a lower-resolution DEM to pixel-level resolution given a monocular image of the coverage with known light source, at the same time we also estimate the corresponding pixel-wise albedo map. We regulate the behaviour of albedo and shape such that the optimized terrain and albedo are the likely solutions that explain the corresponding image. The parameters in the approach are optimized through a kernel-based relaxation framework to gain computational advantages. In this research we experimentally employ the Lunar-Lambertian model for reflectance modelling; the framework of the algorithm is expected to be independent of a specific reflectance model. Experiments are carried out using the monocular images from Lunar Reconnaissance Orbiter (LRO) Narrow Angle Camera (NAC) (0.5 m spatial resolution), constrained by the SELENE and LRO Elevation Model (SLDEM 2015) of 60 m spatial resolution. The results indicate that local details are largely recovered by the algorithm while low frequency topographic consistency is affected by the low-resolution DEM.

Interferometric synthetic aperture radar (InSAR) processing, if lacking high quality ground control points (GCPs), may produce large errors in the final DEM. Some kind of alignment or registration can reduce these errors. This paper evaluates the accuracy of InSAR processed digital elevation models (DEM), against a high resolution DEM. InSAR DEMs were aligned and least-squares registered with USGS, 1/3 arc second NED. The accuracy evaluations after registration show InSAR DEMs alignment and registration can eliminate DEM errors caused by lack of accurate GCPs. In some cases, this process improves the accuracy of the InSAR DEM even more efficiently than including GCPs into InSAR processing.