顾及地形坡度的SRTM和ASTER加权融合

为了克服现有SRTM和ASTER各自缺陷,提升公共DEM精度,本文提出了一种顾及地形坡度的SRTM和ASTER加权融合方法。首先对两种DEM进行地理配准;然后计算不同坡度等级下SRTM和ASTER的高程误差,并得到DEM融合权重;最后采用加权平均法对SRTM和ASTER进行融合。高精度控制点的检验表明：融合后DEM精度有明显提高,相比于原始SRTM和ASTER高程误差,标准差分别降低了5.65 m和1.20 m。     

The Effect of DEM Raster Resolution on First Order, Second Order and Compound Terrain Derivatives

It is well known that the grid cell size of a raster digital elevation model has significant effects on derived terrain variables such as slope, aspect, plan and profile curvature or the wetness index. In this paper the quality of DEMs derived from the interpolation of photogrammetrically derived elevation points in Alberta, Canada, is tested. DEMs with grid cell sizes ranging from 100 to 5 m were interpolated from 100 m regularly spaced elevation points and numerous surface‐specific point elevations using the ANUDEM interpolation method. In order to identify the grid resolution that matches the information content of the source data, three approaches were applied: density analysis of point elevations, an analysis of cumulative frequency distributions using the Kolmogorov‐Smirnov test and the root mean square slope measure. Results reveal that the optimum grid cell size is between 5 and 20 m, depending on terrain com‐plexity and terrain derivative. Terrain variables based on 100 m regularly sampled elevation points are compared to an independent high‐resolution DEM used as a benchmark. Subsequent correlation analysis reveals that only elevation and local slope have a strong positive relationship while all other terrain derivatives are not represented realistically when derived from a coarse DEM. Calculations of root mean square errors and relative root mean square errors further quantify the quality of terrain derivatives.     

DEM Fusion using a modified k-means clustering algorithm

Digital elevation models (DEMs) are a necessary dataset for modelling the Earth’s surface; however, all DEMs contain error. Researchers can reduce this error using DEM fusion techniques since numerous DEMs can be available for a region. However, the use of a clustering algorithm in DEM fusion has not been previously reported. In this study a new DEM fusion algorithm based on a clustering approach that works on multiple DEMs to exploit consistency in the estimates as indicators of accuracy and precision is presented. The fusion approach includes slope and elevation thresholding, k-means clustering of the elevation estimates at each cell location, as well as filtering and smoothing of the fusion product. Corroboration of the input DEMs, and the products of each step of the fusion algorithm, with a higher accuracy reference DEM enabled a detailed analysis of the effectiveness of the DEM fusion algorithm. The main findings of the research were: the k-means clustering of the elevations reduced the precision which also impacted the overall accuracy of the estimates; the number of final cluster members and the standard deviation of elevations before clustering both had a strong relationship to the error in the k-means estimates.     

An Investigation into the Causes of Errors and Inconsistencies in Predicted Viewsheds

Previous evaluations of viewshed analyses have raised concerns about the accuracy and repeatability of the process. Digital elevation model (DEM) errors, the limited spatial resolution of DEMs, and differing algorithms employed by different GIS packages have all been suggested as possible sources for inaccuracy and non‐repeatability. This study compared a field surveyed viewshed to predicted viewsheds generated using a variety of software packages and DEM databases, some of which contained known amounts of error. We found that each of the factors suggested by previous authors contributes to errors in predicted viewsheds. DEM errors contribute most to the discrepancies between surveyed and predicted viewsheds, and the majority of their negative impact occurred at very low levels of DEM error. Differing algorithms used by different GIS packages also contribute significantly to surveyed/predicted viewshed discrepancies, but more importantly, result in predicted viewsheds that disagree with one other, thereby confounding comparisons of results generated with differing software systems. Finally, the spatial resolution of DEMs also has a significant effect on the degree of agreement between surveyed and predicted viewsheds, but the magnitude of this effect is not as great as are the effects produced by DEM errors.     

Effect of Contour Intervals and Grid Cell Size on the Accuracy of DEMs and Slope Derivatives

Digital Elevation Models (DEMs) are indispensable tools in many environmental and natural resource applications. DEMs are frequently derived from contour lines. The accuracy of such DEMs depends on different factors. This research investigates the effect of sampling density used to derive contours, vertical interval between contours (spacing), grid cell size of the DEM (resolution), terrain complexity, and spatial filtering on the accuracy of the DEM and the slope derivative. The study indicated different alternatives to achieve an acceptable accuracy depending on the contour interval, the DEM resolution and the complexity of the terrain. The effect of these factors on the accuracy of the DEM and the slope derivative was quantified using models that determine the level of accuracy (RMSE). The implementation of the models will guide users to select the best combination to improve the results in areas with similar topography. For areas with variable terrain complexity, the suggestion is to generate DEMs and slope at a suitable resolution for each terrain separately and then to merge the results to produce one final layer for the whole area. This will provide accurate estimates of elevation and slope, and subsequently improve the analyses that rely on these digital derivatives.     

Spheroidal equal angular DEMs: The specificity of morphometric treatment

Digital elevation models (DEMs) are commonly constructed using two main types of regular grids: plane square grids and spheroidal equal angular grids. Methods and algorithms intended for plane square‐gridded DEMs should not be directly applied to spheroidal equal angular DEMs. This is because these grids have fundamentally different geometry. However, some researchers continue to apply square‐grid algorithms to spheroidal equal angular DEMs. It seems appropriate to consider once again the specifity of morphometric treatment of spheroidal equal angular DEMs. This article, first, demonstrates possibilities of direct calculation of local, nonlocal, and combined morphometric variables from spheroidal equal angular DEMs exemplified by slope gradient, catchment area, and topographic index. Second, the article shows computational errors when algorithms for plane square‐gridded DEMs are unreasonably applied to spheroidal equal angular DEMs. The study is exemplified by two DEMs. A medium‐resolution DEM of a relatively small, high‐mountainous area (Mount Elbrus) was extracted from the SRTM1 DEM. A low‐resolution DEM of a vast region with the diverse topography (the central and western regions of Kenya) was extracted from the SRTM30_PLUS DEM. The results show that application of square‐grid methods to spheroidal equal angular DEMs leads to substantial computational errors in models of morphometric variables.     

多源DEM融合的高差拟合神经网络方法

本文侧重于介绍智能化摄影测量机器学习的高差拟合神经网络方法。观测手段和处理方式等限制导致全球高质量无缝DEM数据的缺乏,进而制约了它在水文、地质、气象及军事等领域的应用。本文提出了一种基于高差拟合神经网络的多源DEM融合方法,尝试融合全球DEM产品SRTM1、ASTER GDEM v2和激光雷达测高数据ICESat GLAS。首先,根据ICESat GLAS的相关参数及与DEM数据的高程差值,结合坡度自适应的思想设置高差阈值对ICESat GLAS进行滤波,剔除异常数据点。然后,以ICESat GLAS数据为控制点,利用神经网络模型拟合ASTER GDEM v2的误差分布。以地形坡度信息和经纬度坐标作为网络输入,ICESat GLAS和ASTER GDEM v2的高程差值作为目标输出,训练得到预测高差,将其与ASTER GDEM v2高程值相加即可获得校正结果。最后,引入TIN差分曲面的方法,利用校正后的ASTER GDEM v2高程值对SRTM1的数据空洞进行填充,融合生成空间无缝DEM。本文通过随机选取数据进行真实试验,对模型进行了精度验证,并给出了处理结果的定量评价和目视效果。结果表明,不论是空洞还是整体区域,本文方法相比其他DEM数据集和其他方法的处理结果都能够在RMSE上表现出优势,同时,本文提出的方法能够有效克服ASTER GDEM中异常值的影响,得到空间无缝DEM。     

基于星载激光雷达和雷达高度计数据的南极冰盖表面高程制图

南极数字高程模型DEMs(Digital Elevation Models)是研究极区大气环流模式,南极冰盖动态变化和南极科学考察非常重要的基础数据。目前,科学家已经发布了五种不同的南极数字表面高程模型。这些数据都是由卫星雷达高度计,激光雷达和部分地面实测数据等制作而成。尽管如此,由于海洋与冰盖交接的南极冰盖边缘区随时间的快速变化,有必要根据新的卫星数据及时更新南极冰盖表面高程数据。因此,我们利用雷达高度计数据(Envisat RA-2)和激光雷达数据(ICESat/GLAS)制作了最新的南极冰盖高程数据。为提高ICESat/GLAS数据的精度,本文采用了五种不同的质量控制指标对GLAS数据进行处理,滤除了8.36%的不合格数据。这五种质量控制指标分别针对卫星定位误差、大气前向散射、饱和度及云的影响。同时,对Envisat RA-2数据进行干湿对流层纠正、电离层纠正、固体潮汐纠正和极潮纠正。针对两种不同的测高数据,提出了一种基于Envisat RA-2和GLAS数据光斑脚印几何相交的高程相对纠正方法,即通过分析GLAS脚印点与Envisat RA-2数据中心点重叠的点对,建立这些相交点对的高度差(GLAS-RA-2)与表征地形起伏的粗糙度之间的相关关系,对具有稳定相关关系的点对进行Envisat RA-2数据的相对纠正。通过分析南极冰盖不同区域的测高点密度,确定最终DEM的分辨率为1000 m。考虑到南极普里兹湾和内陆地区的差异性,将南极冰盖分为16个区,利用半方差分析确定最佳插值模型和参数,采用克吕金插值方法生成了1000 m分辨率的南极冰盖高程数据。利用两种机载激光雷达数据和我国多次南极科考实测的GPS数据对新的南极DEM进行了验证。结果显示,新的DEM与实测数据的差值范围为3.21—27.84 m,其误差分布与坡度密切关系。与国际上发布的南极DEM数据相比,新的DEM在坡度较大地区和快速变化的冰盖边缘地区精度有较大改进。     

A case study of using external DEM in InSAR DEM generation

Synthetic aperture radar interferometry (InSAR) has been used as an innovative technique for digital elevation model (DEM) and topographic map generation. In this paper, external DEMs are used for InSAR DEM generation to reduce the errors in data processing. The DEMs generated from repeat-pass InSAR are compared. For steep slopes and severe changes in topography, phase unwrapping quality can be improved by subtracting the phase calculated from an external DEM. It is affirmative that the absolute height accuracy of the InSAR DEM is improved by using external DEM. The data processing was undertaken without the use of ground control points and other manual operation.     

机载LiDAR点云密度对DEM精度的影响

机载激光雷达技术(LiDAR)作为一项先进的遥感技术,是植被覆盖区DEM获取的重要手段之一,而不同地形坡度条件及点云密度对DEM产品质量有重要影响。本文以辽宁省某市的机载LiDAR数据为基础,选取5种不同地形坡度的点云数据,通过随机、等间距及基于曲率3种不同的点云抽稀方法,按照点云保留率为80%、60%、40%、20%和10%共5个不同梯度的抽稀倍数对原始点云进行抽稀简化处理,生成与之对应的DEM并对其进行精度评价,以此研究地形坡度、点云抽稀方法、抽稀倍数对DEM精度的影响。结果表明,DEM精度与地形坡度呈负相关关系,即RMSE随地形坡度升高不断增加;基于曲率的抽稀方法在地形坡度＞30°时,相较于其他两种方法RMSE较小,具有明显优势;40%的点云保留率是平衡DEM精度与数据存储效率的一个节点,当点云保留率<40%时,DEM的高程RMSE会迅速增大。该研究对于利用机载LiDAR进行大范围DEM生产具有一定的指导和借鉴意义。     

A practical method for SRTM DEM correction over vegetated mountain areas

Digital elevation models (DEMs) are essential to various applications in topography, geomorphology, hydrology, and ecology. The Shuttle Radar Topographic Mission (SRTM) DEM data set is one of the most complete and most widely used DEM data sets; it provides accurate information on elevations over bare land areas. However, the accuracy of SRTM data over vegetated mountain areas is relatively low as a result of the high relief and the penetration limitation of the C-band used for obtaining global DEM products. The objective of this study is to assess the performance of SRTM DEMs and correct them over vegetated mountain areas with small-footprint airborne Light Detection and Ranging (Lidar) data, which can develop elevation products and vegetation products [e.g., vegetation height, Leaf Area Index (LAI)] of high accuracy. The assessing results show that SRTM elevations are systematically higher than those of the actual land surfaces over vegetated mountain areas. The mean difference between SRTM DEM and Lidar DEM increases with vegetation height, whereas the standard deviation of the difference increases with slope. To improve the accuracy of SRTM DEM over vegetated mountain areas, a regression model between the SRTM elevation bias and vegetation height, LAI, and slope was developed based on one control site. Without changing any coefficients, this model was proved to be applicable in all the nine study sites, which have various topography and vegetation conditions. The mean bias of the corrected SRTM DEM at the nine study sites using this model (absolute value) is 89% smaller than that of the original SRTM DEM, and the standard deviation of the corrected SRTM elevation bias is 11% smaller.     

Effect of sensor modelling methods on computation of 3-D coordinates from Cartosat-1 stereo data

The orbital and the rational polynomial coefficients (RPC) models are the two most commonly used models to compute a three-dimensional coordinates from an image stereo-pair. But it is still confusing that with the identical user provided inputs, which one of these two models provides more accurate digital elevation model (DEM), especially for mountainous terrain. This study aimed to find out the answer by evaluating the impact of used models on the vertical accuracy of DEM extracted from Cartosat-1 stereo data. We used high-accuracy photogrammetric DEM as the reference DEM. Apart from general variations in statistics, surprisingly in a few instances, both the DEMs provided contrasting results, thus proving the significance of this study. The computed root mean square errors and linear error at 90% (LE90) were lower in case of RPC DEM for various classes of slope, aspect and land cover, thus suggesting its better relative accuracy.     

Assessing uncertainty due to elevation error in a landslide susceptibility model

A landslide susceptibility model, employing a digital elevation model (DEM) and geological data, was used in a GIS to predict slope stability in a region of the H J Andrews Long-Term Research Forest, located in the Western Cascade Range in Oregon, USA. To evaluate the contribution of error in elevation to the uncertainty of the final output of the model, several different, but equally probable, versions of the input DEM were created through the addition of random, spatially autocorrelated noise (error) files. The realized DEMs were then processed to produce a family of slope stability maps from which the uncertainty effects of elevation error upon landslide susceptibility could be assessed. The ability to assess this uncertainty has the potential to help us better understand the inherent strengths and weaknesses of applying digital data and spatial information systems to this application, and to facilitate improved natural resource management decisions in relation to timber harvesting and slope stability problems.     

Accuracy assessment of lidar-derived digital elevation models

Despite the relatively high cost of airborne lidar-derived digital elevation models (DEMs), such products are usually presented without a satisfactory associated estimate of accuracy. For the most part, DEM accuracy estimates are typically provided by comparing lidar heights against a finite sample of check point coordinates from an independent source of higher accuracy, supposing a normal distribution of the derived height differences or errors. This paper proposes a new methodology to assess the vertical accuracy of lidar DEMs using confidence intervals constructed from a finite sample of errors computed at check points. A non-parametric approach has been tested where no particular error distribution is assumed, making the proposed methodology especially applicable to non-normal error distributions of the type usually found in DEMs derived from lidar. The performance of the proposed model was experimentally validated using Monte Carlo simulation on 18 vertical error data-sets. Fifteen of these data-sets were computed from original lidar data provided by the International Society for Photogrammetry and Remote Sensing Working Group III/3, using their respective filtered reference data as ground truth. The three remaining data-sets were provided by the Natural Environment Research Council's Airborne Research and Survey Facility lidar system, together with check points acquired using high precision kinematic GPS. The results proved promising, the proposed models reproducing the statistical behaviour of vertical errors of lidar using a favourable number of check points, even in the cases of data-sets with non-normally distributed residuals. This research can therefore be considered as a potentially important step towards improving the quality control of lidar-derived DEMs.     

不同精度DEM流域地形特征提取分析及坡度误差控制研究

以北京密云水库牤牛河流域为试验区,对该流域4种精度DEM提取的流域地形值进行分析,并基于此对不同精度DEM所提取的坡度值进行了误差量化控制,运用曲线回归方法,建立了坡度损失拟合方程系.试验证明,该方法对坡度误差控制较好,具有较大的应用价值.     

Underlying topography extraction over forest areas from multi-baseline PolInSAR data

In this paper, the digital elevation model (DEM) for a forest area is extracted from multi-baseline (MB) polarimetric interferometric synthetic aperture radar (PolInSAR) data. On the basis of the random-volume-over-ground (RVoG) model, the weighted complex least-squares adjustment (WCLSA) method is proposed for the ground phase estimation, so that the MB PolInSAR observations can be constrained by a generalized observation function and the observation contribution to the solution can be adjusted by a weighting strategy. A baseline length weighting strategy is then adopted to syncretize the DEMs estimated with the ground phases. The results of the simulated experiment undertaken in this study demonstrate that the WCLSA method is sensitive to the number of redundant observations and can adjust the contributions of the different observations. We also applied the WCLSA method to E-SAR L- and P-band MB PolInSAR data from the Krycklan River catchment in Northern Sweden. The results show that the two extracted DEMs are in close agreement with the Light Detection and Ranging (Lidar) DEM, with root-mean-square errors of 3.54 and 3.16 m. The DEM vertical error is correlated with the terrain slope and ground-cover condition, but not with the forest height.     

Fusion of high-resolution DEMs derived from COSMO-SkyMed and TerraSAR-X InSAR datasets

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.     

DEM地形信息提取对比研究——以坡度为例

由于DEM数据本身多尺度因素,加之地形、地貌特征具有宏观性与区域分异性的特点,直接的信息提取往往很难达到预期的目的。利用DEM制作坡度图高效、省力,但其精度有很大的不确定性,同时DEM制作过程中的误差传播、转移对坡度信息的影响缺少系统的判断依据。选取位于陕北黄土高原上的两个不同地区作为实验样区,在不同DEM生产的基础上,以高精度的1∶10 000DEM为准值,通过对1∶5万和1∶1万DEM提取定量地形要素的叠合、比较与统计分析,探讨具有不同地貌类型的区域1∶5万DEM提取地形信息的精度及其统计意义上的数量百分比关系。     

Vertical accuracy evaluation of freely available latest high-resolution (30 m) global digital elevation models over Cameroon (Central Africa) with GPS/leveling ground control points.

Digital Elevation Models (DEMs) contain topographic relief data that are vital for many geoscience applications. This study relies on the vertical accuracy of publicly available latest high-resolution (30?m) global DEMs over Cameroon. These models are (1) the ALOS World 3D-30?m (AW3D30), (2) the Shuttle Radar Topography Mission 1 Arc-Second C-Band Global DEM (SRTM 1) and (3) the Advanced Spaceborne Thermal Emission and Reflection Global DEM Version 2 (ASTER GDEM 2). After matching their coordinate systems and datums, the horizontal positional accuracy evaluation was carried out and it shows that geolocation errors significantly influence the vertical accuracy of global DEMs. After this, the three models are compared among them, in order to access random and systematic effects in the elevation data each of them contains. Further, heights from 555 GPS/leveling points distributed all over Cameroon are compared to each DEM, for their vertical accuracy determination. Traditional and robust statistical measures, normality test, outlier detection and removal were used to describe the vertical quality of the DEMs. The test of the normality rejected the hypothesis of normal distribution for all tested global DEMs. Overall vertical accuracies obtained for the three models after georeferencing and gross error removal in terms of Root Mean Square (RMS) and Normalized Median Absolute Deviation (NMAD) are: AW3D30 (13.06?m and 7.75?m), SRTM 1 (13.25?m and 7.41?m) and ASTER GDEM 2 (18.87?m and 13.30?m). Other accuracy measures (MED, 68.3% quantile, 95% quantile) supply some evidence of the good quality of AW3D30 over Cameroon. Further, the effect of land cover and slope on DEM vertical accuracy was also analyzed. All models have proved to be worse in the areas dominated by forests and shrubs areas. SRTM 1 and AW3D30 are more resilient to the effects of the scattering objects respectively in forests and cultivated areas. The dependency of DEMs accuracy on the terrain roughness is evident. In all slope intervals, AW3D30 is performing better than SRTM 1 and ASTER GDEM 2 over Cameroon. AW3D30 is more representative of the external topography over Cameroon in comparison with two others datasets and SRTM 1 can be a serious alternative to AW3D30 for a range of DEM applications in Cameroon.