单流向法地形湿度指数尺度效应的不同地形区差异分析

利用DEM提取地形湿度指数(Topographic Wetness Index,TWI),以满足既定地理模型构建,需考虑不同地形区TWI提取的尺度效应。该文根据实际应用中的DEM分辨率尺度,生成1m、2.5m、5m、10m、20m、30m分辨率尺度的6景DEM;利用单流向算法提取研究区TWI;通过定性和定量方法分析了在不同DEM分辨率尺度中,居民地、梯田、湖泊和自然地形区4种地形区域中TWI值分布特征。结果表明:不同地形区TWI的尺度效应差异较大;随着DEM尺度增大,TWI值域逐渐缩小,变异程度逐渐增强,反映地表形态变化的能力逐渐减小;当DEM尺度大于5m时,TWI基本不能反映人工地形区域;自然地形区的TWI变化规律与研究区整体TWI变化规律相当,居民地与梯田的TWI变化规律大致相同,湖泊区域的TWI变化相对稳定。     

面向DEM地形复杂度分析的分形方法研究

在基于栅格DEM数据的基础上,应用分形几何学方法,采用元分维模型理论,提出一种描述DEM地形复杂度的分形分析方法,并得出分析指标———地形分维指数(TFI)。实验证明,该指标可以有效描述栅格DEM数据反映的地形变化特征,分析窗口由小到大的规律性变化反映出地貌从微观到宏观的变化情况,因此其既可以作为对栅格DEM所描述地形进行评价的坡面因子,又可以作为基于DEM的地貌类型区自动划分的参照依据。     

地形信息对确定DEM适宜分辨率的影响

分辨率会直接影响基于栅格数字高程模型（DEM）的数字地形分析结果,因此在实际应用中,需要选择适宜的DEM分辨率。目前采取的基本方法,基于某种地形信息定量刻画尺度效应曲线,从而确定DEM适宜分辨率,但对于采用不同地形信息时所产生的影响尚缺乏研究。本文针对该方法中通常采用的坡度、剖面曲率、水平曲率等3 种地形信息,每种地形信息提取时,分别使用两种不同的常用算法,在3 个不同地形特征的研究区中,逐一计算其在不同分辨率下的局部方差均值,以刻画尺度效应曲线,确定相应的DEM适宜分辨率,并进行对比分析。结果表明：① 采用剖面曲率或水平曲率所得适宜分辨率结果基本相同,但采用坡度所得出的适宜分辨率结果则有明显差别,后者所得的适宜分辨率更粗;② 采用不同地形信息时,越是在平缓地形为主的研究区,所得的适宜分辨率结果越相近,在复合地形特征的研究区所得到的适宜分辨率区间均明显较宽;③ 地形属性计算时所用的算法对适宜分辨率结果的影响不明显。     

DEM 点位地形信息量化模型研究

针对DEM 点位, 首先应用微分几何法对其所负载的语法信息量进行测度, 其次根据地形特征点类型及地形结构特征确定其语义信息量, 然后基于信息学理论构建了DEM 点位地形信息综合量化模型。在此基础上, 以黄土丘陵沟壑区作为实验样区, 对DEM 点位地形信息量提取方法及其在地形简化中的初步实例应用进行了探讨和验证。实验结果显示, 所提出的DEM 点位地形信息量化方案可行;基于DEM 地形信息量指数的多尺度DEM 构建方案, 具有机理明确、易于实现的特点, 并通过优先保留地形骨架特征点, 可以有效减少地形失真, 从而满足不同层次的多尺度数字地形建模和表达要求。对DEM 点位地形信息进行有效量化, 为认识DEM 地形信息特征提供了一个新的切入点, 同时为多尺度数字地形建模提供理论依据与方法支持。     

DEM栅格单元地形异质性的量度指标研究

DEM栅格单元地形异质性可以理解为栅格单元内部地形的复杂程度,在数字地形分析中,其存在可能得出与实际不符的分析结果.然而,目前还没有一种综合量度指标可以衡量DEM栅格单元地形异质性的大小.该文在总结相关研究成果的基础上,分析了DEM栅格单元地形异质性产生的主要原因,并采用高程标准差、地形起伏度、地表粗糙度和平均坡度4个指标度量DEM栅格单元地形异质性的不同方面,最后通过归一化综合处理,得到了DEM栅格单元地形异质性的综合量度指标--DEM栅格单元地形异质性指数(DGTHI),并进行了实验验证与分析.     

地形复杂度对坡度坡向的影响

采用三阶不带权差分算法，研究了地形复杂度与坡度坡向的关系，澄清了目前关于坡度坡向误差空间分布的矛盾观点，并分别在凹向椭球和高斯合成曲面数学模型曲面DEM上对其进行验证。通过研究得出：①坡度、坡向误差与坡度值正相关；②坡度坡向误差主要分布在平坦地区；③坡向误差较坡度误差对DEM高程数据误差敏感，较小的DEM误差引起较大的坡向误差。     

基于DEM的地形单元多样性指数及其算法

在阐述地形信息表达研究进展的基础上,提出基于DEM地形单元多样性指数的概念和算法。地形单元多样性指数综合了高程、坡度、坡位、坡向、汇流量和水域信息等要素。其算法集成地形位置指数和地形湿度指数算法,采用图层叠加分析,设定分类、分级指标,进行重分类组合,划分地形单元类型,利用窗口分析法计算地形单元多样性指数。以四川省为试验区,利用精度为100 m的DEM数据和水域分布数据进行模拟计算,地形单元划分为13种典型类型,统计窗口半径设为900 m,计算出的多样性指数值小于0.5的区域仅占总面积的11%,大于0.75的区域达57%,符合实验区地形特征,并对算法进行了可行性分析和验证。结果表明,该算法提取的地形单元多样性指数可以有效反映地表形态的多样性特征及其变化。该研究结果为进一步探讨基于DEM地表形态信息的概念体系,以及从微观到宏观的地形信息空间分析研究创造了条件。     

TOPMODEL中地形指数ln(α/tanβ)的新算法

TOPMODEL中地形指数ln(α/tanβ)被用来近似表征流域径流源面积和地下水水位的空间分布特征,目前广泛使用的计算地形指数的方法为多流向算法(FD8算法).本文首先介绍了多流向算法的基本原理,并基于此对流动累积分配中的有效等高线长度精确计算提出了几何锥面内切圆算法,同时改进了传统的地形指数中单位等高线汇流面积α的计算方法,增强了多流向算法对DEM中异常栅格的处理能力.改进后的地形指数新算法在两个不同流域和不同分辨率DEM上与传统多流向算法进行了对比分析.结果表明新算法在原理上更符合地形指数物理意义,实际应用中其计算结果更准确.这种地形指数新算法的提出对于流域水文过程机理分析及陆面过程定量化研究具有一定的现实意义.     

基于DEM的中国地形起伏度适宜计算尺度研究

基于SRTM和ASTER DEM数据,在全国范围内选取13个实验区,在渐变尺度下计算平均起伏度变化曲线的"突变点",据此确定中国地形起伏度的适宜计算尺度;结合山地界定标准计算各实验区山地面积,并采用人工解译的山地范围对计算结果进行检验。研究结果表明:1)地形起伏度适宜计算尺度与所采用的DEM数据有关,DEM分辨率越小,地形起伏度适宜计算尺度越大;2)针对同一分辨率DEM数据,中国境内的地形起伏度适宜计算尺度随地貌特征变化而变化,但总体变化幅度不大;3)针对SRTM和ASTER DEM两种常用数据源,分别选择4.72km2和3.20km2作为地形起伏度适宜计算尺度是合理的,山地界定精度达90%以上。     

基于SPIHT小波的DEM自适应压缩方法研究

海量地形数据给其存储、分发和实时渲染带来了极大的挑战,因此迫切需要适合网络环境下地形可视化的数据压缩方法。该文探讨地形复杂度与DEM压缩方法的关系,研究可视化中地形复杂度的计算方法,并提出一种改进的SPIHT小波压缩方法:采用小波分解后的系数对DEM的地形复杂度进行评估,并针对地形复杂度对编码算法进行自适应调节。实验证明,这种改进的SPIHT小波压缩方法采用合适的压缩比进行DEM数据压缩,能够在满足地形可视化需要的同时提高压缩效率。     

DEM流径算法的相似性分析

流径算法是分布式水文模型、土壤侵蚀模拟等研究中的关键技术环节,决定着汇水面积、地形指数等许多重要的地形、水文参数的计算。本文以黄土高原两个典型样区的不同分辨率DEM为研究对象,对常用的五种流径算法(D8、Rho8、Dinf、MFD和DEMON)通过相对差系数、累积频率图、XY散点分布图等进行了定量的对比分析。结果表明:算法的差异主要集中在坡面区域,汇流区域各类算法的差别较小;算法差异在不同DEM尺度下都有所体现,但高分辨率下的差异会更明显;在地形复杂区域,多流向算法要优于单流向算法。研究也进一步指出汇水面积、地形指数等水文参数对流径算法具有强烈的依赖性。     

An adaptive approach to selecting a flow‐partition exponent for a multiple‐flow‐direction algorithm

Most multiple‐flow‐direction algorithms (MFDs) use a flow‐partition coefficient (exponent) to determine the fractions draining to all downslope neighbours. The commonly used MFD often employs a fixed exponent over an entire watershed. The fixed coefficient strategy cannot effectively model the impact of local terrain conditions on the dispersion of local flow. This paper addresses this problem based on the idea that dispersion of local flow varies over space due to the spatial variation of local terrain conditions. Thus, the flow‐partition exponent of an MFD should also vary over space. We present an adaptive approach for determining the flow‐partition exponent based on local topographic attribute which controls local flow partitioning. In our approach, the influence of local terrain on flow partition is modelled by a flow‐partition function which is based on local maximum downslope gradient (we refer to this approach as MFD based on maximum downslope gradient, MFD‐md for short). With this new approach, a steep terrain which induces a convergent flow condition can be modelled using a large value for the flow‐partition exponent. Similarly, a gentle terrain can be modelled using a small value for the flow‐partition exponent. MFD‐md is quantitatively evaluated using four types of mathematical surfaces and their theoretical ‘true’ value of Specific Catchment Area (SCA). The Root Mean Square Error (RMSE) shows that the error of SCA computed by MFD‐md is lower than that of SCA computed by the widely used SFD and MFD algorithms. Application of the new approach using a real DEM of a watershed in Northeast China shows that the flow accumulation computed by MFD‐md is better adapted to terrain conditions based on visual judgement.     

Pre-processing algorithms and landslide modelling on remotely sensed DEMs

Terrain analysis applications using remotely sensed Digital Elevation Models (DEMs), nowadays easily available, permit to quantify several river basin morphologic and hydrologic properties (e.g. slope, aspect, curvature, flow path lengths) and indirect hydrogeomorphic indices (e.g. specific upslope area, topographic wetness index) able to characterize the physical processes governing the landscape evolution (e.g. surface saturation, runoff, erosion, deposition). Such DEMs often contain artifacts and the automated hydrogeomorphic characterization of the watershed is influenced by terrain analysis procedures consisting in artificial depression (pit) and flat area treatment approaches combined with flow direction methods.In shallow landslide deterministic models, when applied using topographic dataset at medium scale (e.g. 30 m of resolution), the choice of the most suitable DEM-processing procedure is not trivial and can influence model results. This also affects the selection of most critical areas for further finer resolution studies or for the implementation of countermeasures aiming to landslide risk mitigation.In this paper such issue is investigated using as topographic input the ASTER DEMs and comparing two different combinations of DEM correction and flow routing schemes. The study areas comprise ten catchments in Italy for which hydrogeomorphic processes are significant. Aims of this paper are: 1) to introduce a parameter estimation procedure for the physically-based DEM correction method PEM4PIT (Physical Erosion Model for PIT removal); 2) to investigate the influence of different terrain analysis procedures on results of the slope stability model SHALSTAB (SHAllow Landsliding STABility) using remotely-sensed ASTER DEMs; 3) trying to assess which of terrain analysis methods is more appropriate for describing terrain instability.     

DEM地形信息量计算的不确定性研究

DEM地形信息量的准确度量是判定DEM数据应用适宜性与限制性关键。以黄土丘陵沟壑区DEM及其派生的坡度、坡向DTM为例进行实验,研究子集划分与分级数对DEM信息量估算精度的影响。结果表明,DEM地形信息量是实际地形信息和子集划分策略共同作用的结果,现有DEM信息量计算存在不确定性,DEM子集划分算法和分级数的选择直接影响DEM信息量的计算结果。文章借助最大熵定理的基本思路构建DEM地形信息量分级判定模型,可客观获得连续型栅格地形数据在信息量计算上最优分级,为有效估算DEM地形信息量提供科学依据。     

不同分辨率DEM提取地面坡度的不确定性模拟:以在黄土高原的试验为例

Slope is one of the crucial terrain variables in spatial analysis and land use planning,especially in the Loess Plateau area of China which is suffering from serious soil erosion. DEM based slope extracting method has been widely accepted and applied in practice. However slope accuracy derived from this method usually does not match with its popularity. A quantitative simulation to slope data uncertainty is important not only theoretically but also necessarily to applications. This paper focuses on how resolution and terrain complexity impact on the accuracy of mean slope extracted from DEMs of different resolutions in the Loess Plateau of China. Six typical geomorphologic areas are selected as test areas, representing different terrain types from smooth to rough. Their DEMs are produced from digitizing contours of 1:10,000 scale topographic maps. Field survey results show that 5 m should be the most suitable grid size for representing slope in the Loess Plateau area. Comparative and math-simulation methodology was employed for data processing and analysis. A linear correlativity between mean slope and DEM resolution was found at all test areas,but their regression coefficients related closely with the terrain complexity of the test areas. If taking stream channel density to represent terrain complexity, mean slope error could be regressed against DEM resolution (X) and stream channel density (S) at 8 resolution levels and expressed as (0.0015S2+0.031S-0.0325)X-0.0045S2-0.155S+0.1625, with a R2 value of over 0.98. Practical tests also show an effective result of this model in applications. The new development methodology applied in this study should be helpful to similar researches in spatial data uncertainty investigation.     

Simulation on slope uncertainty derived from DEMs at different resolution levels: a case study in the Loess Plateau

Slope is one of the crucial terrain variables in spatial analysis and land use planning, especially in the Loess Plateau area of China which is suffering from serious soil erosion. DEM based slope extracting method has been widely accepted and applied in practice. However slope accuracy derived from this method usually does not match with its popularity. A quantitative simulation to slope data uncertainty is important not only theoretically but also necessarily to applications. This paper focuses on how resolution and terrain complexity impact on the accuracy of mean slope extracted from DEMs of different resolutions in the Loess Plateau of China. Six typical geomorphologic areas are selected as test areas, representing different terrain types from smooth to rough. Their DEMs are produced from digitizing contours of 1:10,000 scale topographic maps. Field survey results show that 5 m should be the most suitable grid size for representing slope in the Loess Plateau area. Comparative and math-simulation methodology was employed for data processing and analysis. A linear correlativity between mean slope and DEM resolution was found at all test areas, but their regression coefficients related closely with the terrain complexity of the test areas. If taking stream channel density to represent terrain complexity, mean slope error could be regressed against DEM resolution (X) and stream channel density (S) at 8 resolution levels and expressed as (0.0015S2+0.031S-0.0325)X-0.0045S2-0.155S+0.1625, with a R2 value of over 0.98. Practical tests also show an effective result of this model in applications. The new development methodology applied in this study should be helpful to similar researches in spatial data uncertainty investigation.     

一种基于坡度分析的DEM数字水印算法

提出一种新的基于坡度分析的DEM数字水印算法.该算法以DEM坡度为主要研究对象,由于坡度误差主要集中在平坦地区,而经小波变换后的DEM低频系数信息可反映DEM区域的地形复杂度,因此,通过对DEM低频系数进行分析,可自适应地确定水印的嵌入位置;然后以加性法则将水印信息嵌入低频系数,再经小波逆变换,即可得到含水印信息的DEM数据.实验结果表明,该算法能完全满足水印的不可见性,对DEM数据的高程精度、坡度精度及等高线的提取精度影响都很小,能够满足一定的应用要求,且具有良好的抗噪能力.     

A study on DEM-derived primary topographic attributes for hydrologic applications: Sensitivity to elevation data resolution

Primary topographic attributes play a critical role in determining watershed hydrologic characteristics for water resources modeling with raster-based digital elevation models (DEM). The effects of DEM resolution on a set of important topographic derivatives are examined in this study, including slope, upslope contributing area, flow length and watershed area. The focus of the study is on how sensitive each of the attributes is to the resolution uncertainty by considering the effects of overall terrain gradient and bias from resampling. Two case study watersheds of different gradient patterns are used with their 10 m USGS DEMs. A series of DEMs up to 200 m grid size are produced from the base DEMs using three commonly used resampling methods. All the terrain variables tested vary with the grid size change. It is found that slope angles decrease and contributing area values increase constantly as DEMs are aggregated progressively to coarser resolutions. No systematic trend is observed for corresponding changes of flow path and watershed area. The analysis also suggests that gradient profile of the watershed presents an important factor for the examined sensitivities to DEM resolution.