Uncertainties in Digital Elevation Models: Evaluation and Effects on Landform and Soil Type Classification

Digital elevation models (DEMs) are a widely used source for the digital representation of the Earth's surface in a wide range of scientific, industrial and military applications. Since many processes on Earth are influenced by the shape of the relief, a variety of different applications rely on accurate information about the topography. For instance, DEMs are used for the prediction of geohazards, climate modelling, or planning-relevant issues, such as the identification of suitable locations for renewable energies. Nowadays, DEMs can be acquired with a high geometric resolution and over large areas using various remote sensing techniques, such as photogrammetry, RADAR, or laser scanning (LiDAR). However, they are subject to uncertainties and may contain erroneous representations of the terrain. The quality and accuracy of the topographic representation in the DEM is crucial, as the use of an inaccurate dataset can negatively affect further results, such as the underestimation of landslide hazards due to a too flat representation of relief in the elevation model. Therefore, it is important for users to gain more knowledge about the accuracy of a terrain model to better assess the negative consequences of DEM uncertainties on further analysis results of a certain research application. A proper assessment of whether the purchase or acquisition of a highly accurate DEM is necessary or the use of an already existing and freely available DEM is sufficient to achieve accurate results is of great qualitative and economic importance. In this context, the first part of this thesis focuses on extending knowledge about the behaviour and presence of uncertainties in DEMs concerning terrain and land cover. Thus, the first two studies of this dissertation provide a comprehensive vertical accuracy analysis of twelve DEMs acquired from space with spatial resolutions ranging from 5 m to 90 m. The accuracy of these DEMs was investigated in two different regions of the world that are substantially different in terms of relief and land cover. The first study was conducted in the hyperarid Chilean Atacama Desert in northern Chile, with very sparse land cover and high elevation differences. The second case study was conducted in a mid-latitude region, the Rur catchment in the western part of Germany. This area has a predominantly flat to hilly terrain with relatively diverse and dense vegetation and land cover. The DEMs in both studies were evaluated with particular attention to the influence of relief and land cover on vertical accuracy. The change of error due to changing slope and land cover was quantified to determine an average loss of accuracy as a function of slope for each DEM. Additionally, these values were used to derive relief-adjusted error values for different land cover classes. The second part of this dissertation addresses the consequences that different spatial resolutions and accuracies in DEMs have on specific applications. These implications were examined in two exemplary case studies. In a geomorphometric case study, several DEMs were used to classify landforms by different approaches. The results were subsequently compared and the accuracy of the classification results with different DEMs was analysed. The second case study is settled within the field of digital soil mapping. Various soil types were predicted with machine learning algorithms (random forest and artificial neural networks) using numerous relief parameters derived from DEMs of different spatial resolutions. Subsequently, the influence of high and low resolution DEMs with the respectively derived land surface parameters on the prediction results was evaluated. The results on the vertical accuracy show that uncertainties in DEMs can have diverse reasons. Besides the spatial resolution, the acquisition technique and the degree of improvements made to the dataset significantly impact the occurrence of errors in a DEM. Furthermore, the relief and physical objects on the surface play a major role for uncertainties in DEMs. Overall, the results in steeper areas show that the loss of vertical accuracy is two to three times higher for a 90 m DEM than for DEMs of higher spatial resolutions. While very high resolution DEMs of 12 m spatial resolution or higher only lose about 1 m accuracy per 10° increase in slope steepness, 30 m DEMs lose about 2 m on average, and 90 m DEMs lose more than 3 m up to 6 m accuracy. However, the results also show significant differences for DEMs of identical spatial resolution depending on relief and land cover. With regard to different land cover classes, it can be stated that mid-latitude forested and water areas cause uncertainties in DEMs of about 6 m on average. Other tested land cover classes produced minor errors of about 1 – 2 m on average. The results of the second part of this contribution prove that a careful selection of an appropriate DEM is more crucial for certain applications than for others. The choice of different DEMs greatly impacted the landform classification results. Results from medium resolution DEMs (30 m) achieved up to 30 % lower overall accuracies than results from high resolution DEMs with a spatial resolution of 5 m. In contrast to the landform classification results, the predicted soil types in the second case study showed only minor accuracy differences of less than 2 % between the usage of a spatial high resolution DEM (15 m) and a low resolution 90 m DEM. Finally, the results of these two case studies were compared and discussed with other results from the literature in other application areas. A summary and assessment of the current state of knowledge about the impact of a particular chosen terrain model on the results of different applications was made. In summary, the vertical accuracy measures obtained for each DEM are a first attempt to determine individual error values for each DEM that can be interpreted independently of relief and land cover and can be better applied to other regions. This may help users in the future to better estimate the accuracy of a tested DEM in a particular landscape. The consequences of elevation model selection on further results are highly dependent on the topic of the study and the study area's level of detail. The current state of knowledge on the impact of uncertainties in DEMs on various applications could be established. However, the results of this work can be seen as a first step and more work is needed in the future to extend the knowledge of the effects of DEM uncertainties on further topics that have not been investigated to date

Extraction Landscape Elements from Remote Sensing Data

V této práci je popsán postup pro automatickou detekci krajinných prvků z dat pořízených bezkontaktními dálkovými metodami. Tato interpretace dat byla provedena v softwaru eCognition Developer prostřednictvím procesu klasifikace. Pro klasifikaci byla využita matoda obektově orientované analýzy, která dělí data takovým způsobem, že přiřazuje informaci o příslušnosti k nějaké třídě, například krajinnému typu, skupinám pixelů - objektům. Klasifikace byla provedena se současným využitím produktů dvou různých mapovacích technik - ortofot pořízených z leteckého snímkování a normalizovaného digitálního modelu povrchu, který byl určen pomocí LiDARU. Bylo identifikováno a klasikováno pět typů krajinných prvků.In this thesis, an approach to automatically derive information about land cover from the remotely sensed data is presented. The data interpretation was done with classification process and performed in software eCognition Developer. The Object-based image analysis, which assignes the classes - for example land cover types, to clusters of pixels (=objects), was used. For the classification, products of two different data sources were combined - the orthophotos generated from aerial imagery and Normalized Digital surface model derived from LiDAR data. Five types of landscape elements were identified and classified.

The Use of Sun Elevation Angle for Stereogrammetric Boreal Forest Height in Open Canopies

Stereogrammetry applied to globally available high resolution spaceborne imagery (HRSI; less than 5 m spatial resolution) yields fine-scaled digital surface models (DSMs) of elevation. These DSMs may represent elevations that range from the ground to the vegetation canopy surface, are produced from stereoscopic image pairs (stereo pairs) that have a variety of acquisition characteristics, and have been coupled with lidar data of forest structure and ground surface elevation to examine forest height. This work explores surface elevations from HRSI DSMs derived from two types of acquisitions in open canopy forests. We (1) apply an automated mass-production stereogrammetry workflow to along-track HRSI stereo pairs, (2) identify multiple spatially coincident DSMs whose stereo pairs were acquired under different solar geometry, (3) vertically co-register these DSMs using coincident spaceborne lidar footprints (from ICESat-GLAS) as reference, and(4) examine differences in surface elevations between the reference lidar and the co-registered HRSI DSMs associated with two general types of acquisitions (DSM types) from different sun elevation angles. We find that these DSM types, distinguished by sun elevation angle at the time of stereo pair acquisition, are associated with different surface elevations estimated from automated stereogrammetry in open canopy forests. For DSM values with corresponding reference ground surface elevation from spaceborne lidar footprints in open canopy northern Siberian Larix forests with slopes less than10, our results show that HRSI DSM acquired with sun elevation angles greater than 35deg and less than 25deg (during snow-free conditions) produced characteristic and consistently distinct distributions of elevation differences from reference lidar. The former include DSMs of near-ground surfaces with root mean square errors less than 0.68 m relative to lidar. The latter, particularly those with angles less than 10deg, show distributions with larger differences from lidar that are associated with open canopy forests whose vegetation surface elevations are captured. Terrain aspect did not have a strong effect on the distribution of vegetation surfaces. Using the two DSM types together, the distribution of DSM-differenced heights in forests (6.0 m, sigma = 1.4 m) was consistent with the distribution of plot-level mean tree heights (6.5m, sigma = 1.2 m). We conclude that the variation in sun elevation angle at time of stereo pair acquisition can create illumination conditions conducive for capturing elevations of surfaces either near the ground or associated with vegetation canopy. Knowledge of HRSI acquisition solar geometry and snow cover can be used to understand and combine stereogrammetric surface elevation estimates to co-register rand difference overlapping DSMs, providing a means to map forest height at fine scales, resolving the vertical structure of groups of trees from spaceborne platforms in open canopy forests

Automated Building Information Extraction and Evaluation from High-resolution Remotely Sensed Data

The two-dimensional (2D) footprints and three-dimensional (3D) structures of buildings are of great importance to city planning, natural disaster management, and virtual environmental simulation. As traditional manual methodologies for collecting 2D and 3D building information are often both time consuming and costly, automated methods are required for efficient large area mapping. It is challenging to extract building information from remotely sensed data, considering the complex nature of urban environments and their associated intricate building structures. Most 2D evaluation methods are focused on classification accuracy, while other dimensions of extraction accuracy are ignored. To assess 2D building extraction methods, a multi-criteria evaluation system has been designed. The proposed system consists of matched rate, shape similarity, and positional accuracy. Experimentation with four methods demonstrates that the proposed multi-criteria system is more comprehensive and effective, in comparison with traditional accuracy assessment metrics. Building height is critical for building 3D structure extraction. As data sources for height estimation, digital surface models (DSMs) that are derived from stereo images using existing software typically provide low accuracy results in terms of rooftop elevations. Therefore, a new image matching method is proposed by adding building footprint maps as constraints. Validation demonstrates that the proposed matching method can estimate building rooftop elevation with one third of the error encountered when using current commercial software. With an ideal input DSM, building height can be estimated by the elevation contrast inside and outside a building footprint. However, occlusions and shadows cause indistinct building edges in the DSMs generated from stereo images. Therefore, a “building-ground elevation difference model” (EDM) has been designed, which describes the trend of the elevation difference between a building and its neighbours, in order to find elevation values at bare ground. Experiments using this novel approach report that estimated building height with 1.5m residual, which out-performs conventional filtering methods. Finally, 3D buildings are digitally reconstructed and evaluated. Current 3D evaluation methods did not present the difference between 2D and 3D evaluation methods well; traditionally, wall accuracy is ignored. To address these problems, this thesis designs an evaluation system with three components: volume, surface, and point. As such, the resultant multi-criteria system provides an improved evaluation method for building reconstruction

Gamma-ray attenuation caused by rainforest dispersion compared to Vegetation Index: estimates on the effects in airborne gamma-spectrometry data – example from the State of Rondônia, Amazonia, Brazil

This work estimates and discusses the effects of the rainforest on the airborne gamma-spectrometric surveys, taking as case study an area in the center of the Rondônia State, Amazonia, northern Brazil, where wooded and deforested areas are frequently juxtaposed. The control of the wooded areas is made using Landsat satellite images, by the calculation of the Normalized Difference Vegetation Index (NDVI), which allows distinguishing between areas with low and high concentration of green leaves. The results show that, compared to the NDVI values for non-vegetated areas, there is an attenuation of the mean signal of this index in the rainforest areas corresponding up to 36% of the Total Count of Radiation (TC), 12% of potassium (K), 37% of equivalent thorium (eTh), and 29% of the equivalent uranium (eU) concentrations. In our interpretation, part of this effect might be caused by the radiation shielding of the biomass in rainforest areas. However, as a secondary effect for the gamma-spectrometry data processing, the canopy trees are considered on the Digital Elevation Model being used as a surface to calculate the distance between the source of radiation and the sensor, which causes a bias in the processing, underestimating the real flight height to the ground surface. This last effect could have been avoided if there were a Digital Terrain Model available for the area, which should be considered for the data processing. Improving the understanding about this phenomenon, an increase in the quality of the signal-to-noise ratio of airborne gamma-spectrometric data used for the quantitative land modeling may be achieved for regions where the presence of rainforests is significant such as the Amazon region northern Brazil

Twenty-Second Lunar and Planetary Science Conference

The papers in this collection were written for general presentation, avoiding jargon and unnecessarily complex terms. Some of the topics covered include: planetary evolution, planetary satellites, planetary composition, planetary surfaces, planetary geology, volcanology, meteorite impacts and composition, and cosmic dust. Particular emphasis is placed on Mars and the Moon

Effects of spatial resolution on radar-based precipitation estimation using sub-kilometer X-band radar measurements

Known for the ability to observe precipitation at spatial resolution higher than rain gauge networks and satellite products, weather radars allow us to measure precipitation at spatial resolutions of 1 kilometer (typical resolution for operational radars) and a few hundred meters (often used in research activities). In principle, we can operate a weather radar at resolution higher than 100m and the expectation is that radar data at higher spatial resolution can provide more information. However, there is no systematic research about whether the additional information is noise or useful data contributing to the quantitative precipitation estimation. In order to quantitatively investigate the changes, as either benefits or drawbacks, caused by increasing the spatial resolution of radar measurements, we set up an X-band radar field experiment from May to October in 2017 in the Stuttgart metropolitan region. The scan strategy consists of two quasi-simultaneous scans with a 75-m and a 250-m radial resolution respectively. They are named as the fine scan and the coarse scan, respectively. Both scans are compared to each other in terms of the radar data quality and their radar-based precipitation estimates. The primary results from these comparisons between the radar data of these two scans show that, in contrast to the coarse scan, the fine scan data are characterized with losses of weak echoes, are more subjected to external signals and second-trip echoes (drawback), are more effective in removing non-meteorological echoes (benefit), are more skillful in delineating convective storms (benefit), and show a better agreement with the external reference data (benefit)

Elevation and Deformation Extraction from TomoSAR

3D SAR tomography (TomoSAR) and 4D SAR differential tomography (Diff-TomoSAR) exploit multi-baseline SAR data stacks to provide an essential innovation of SAR Interferometry for many applications, sensing complex scenes with multiple scatterers mapped into the same SAR pixel cell. However, these are still influenced by DEM uncertainty, temporal decorrelation, orbital, tropospheric and ionospheric phase distortion and height blurring. In this thesis, these techniques are explored. As part of this exploration, the systematic procedures for DEM generation, DEM quality assessment, DEM quality improvement and DEM applications are first studied. Besides, this thesis focuses on the whole cycle of systematic methods for 3D & 4D TomoSAR imaging for height and deformation retrieval, from the problem formation phase, through the development of methods to testing on real SAR data. After DEM generation introduction from spaceborne bistatic InSAR (TanDEM-X) and airborne photogrammetry (Bluesky), a new DEM co-registration method with line feature validation (river network line, ridgeline, valley line, crater boundary feature and so on) is developed and demonstrated to assist the study of a wide area DEM data quality. This DEM co-registration method aligns two DEMs irrespective of the linear distortion model, which improves the quality of DEM vertical comparison accuracy significantly and is suitable and helpful for DEM quality assessment. A systematic TomoSAR algorithm and method have been established, tested, analysed and demonstrated for various applications (urban buildings, bridges, dams) to achieve better 3D & 4D tomographic SAR imaging results. These include applying Cosmo-Skymed X band single-polarisation data over the Zipingpu dam, Dujiangyan, Sichuan, China, to map topography; and using ALOS L band data in the San Francisco Bay region to map urban building and bridge. A new ionospheric correction method based on the tile method employing IGS TEC data, a split-spectrum and an ionospheric model via least squares are developed to correct ionospheric distortion to improve the accuracy of 3D & 4D tomographic SAR imaging. Meanwhile, a pixel by pixel orbit baseline estimation method is developed to address the research gaps of baseline estimation for 3D & 4D spaceborne SAR tomography imaging. Moreover, a SAR tomography imaging algorithm and a differential tomography four-dimensional SAR imaging algorithm based on compressive sensing, SAR interferometry phase (InSAR) calibration reference to DEM with DEM error correction, a new phase error calibration and compensation algorithm, based on PS, SVD, PGA, weighted least squares and minimum entropy, are developed to obtain accurate 3D & 4D tomographic SAR imaging results. The new baseline estimation method and consequent TomoSAR processing results showed that an accurate baseline estimation is essential to build up the TomoSAR model. After baseline estimation, phase calibration experiments (via FFT and Capon method) indicate that a phase calibration step is indispensable for TomoSAR imaging, which eventually influences the inversion results. A super-resolution reconstruction CS based study demonstrates X band data with the CS method does not fit for forest reconstruction but works for reconstruction of large civil engineering structures such as dams and urban buildings. Meanwhile, the L band data with FFT, Capon and the CS method are shown to work for the reconstruction of large manmade structures (such as bridges) and urban buildings

COBE's search for structure in the Big Bang

The launch of Cosmic Background Explorer (COBE) and the definition of Earth Observing System (EOS) are two of the major events at NASA-Goddard. The three experiments contained in COBE (Differential Microwave Radiometer (DMR), Far Infrared Absolute Spectrophotometer (FIRAS), and Diffuse Infrared Background Experiment (DIRBE)) are very important in measuring the big bang. DMR measures the isotropy of the cosmic background (direction of the radiation). FIRAS looks at the spectrum over the whole sky, searching for deviations, and DIRBE operates in the infrared part of the spectrum gathering evidence of the earliest galaxy formation. By special techniques, the radiation coming from the solar system will be distinguished from that of extragalactic origin. Unique graphics will be used to represent the temperature of the emitting material. A cosmic event will be modeled of such importance that it will affect cosmological theory for generations to come. EOS will monitor changes in the Earth's geophysics during a whole solar color cycle

Monitoring and predicting railway subsidence using InSAR and time series prediction techniques

Improvements in railway capabilities have resulted in heavier axle loads and higher speed operations, which increase the dynamic loads on the track. As a result, railway subsidence has become a threat to good railway performance and safe railway operation. The author of this thesis provides an approach for railway performance assessment through the monitoring and prediction of railway subsidence. The InSAR technique, which is able to monitor railway subsidence over a large area and long time period, was selected for railway subsidence monitoring. Future trends of railway subsidence should also be predicted using subsidence prediction models based on the time series deformation records obtained by InSAR. Three time series prediction models, which are the ARMA model, a neural network model and the grey model, are adopted in this thesis. Two case studies which monitor and predict the subsidence of the HS1 route were carried out to assess the performance of HS1. The case studies demonstrate that except for some areas with potential subsidence, no large scale subsidence has occurred on HS1 and the line is still stable after its 10 years' operation. In addition, the neural network model has the best performance in predicting the subsidence of HS1