基于数字高程模型的中国岩溶地貌研究进展及前景分析

数字高程模型（DEM）蕴含丰富的地形地貌信息,基于DEM的数字地形分析方法为岩溶地貌研究提供了科学、有效的技术手段。文章针对前人应用DEM研究中国岩溶地貌所涉及的关键技术方法,从岩溶地貌识别的尺度效应、岩溶地貌的识别与分类、岩溶地貌的形态及格局分析、岩溶区生态环境变化等方面进行了梳理和分析,提出未来应构建科学的岩溶地貌数字分类体系,对岩溶地貌进行多尺度、深层次的地形分析和定量表达,并从地形现状研究拓展到地形演变的过程与机理研究,发掘出DEM在岩溶地貌研究中更多的应用。      

水沙物理模型CASC2D-SED的DEM尺度效应发生机制

DEM分辨率对分布式水沙过程模拟具有重要影响,然而,产生影响的内部机制尚不明确。改进水沙物理模型CASC2D-SED的结构,将坡度由DEM在模型内部直接提取改为由模块单独计算,并将坡度设计为模型的独立输入参数,通过单独改变坡度参数来研究坡度对水沙模拟DEM尺度效应的影响。基于改进的CASC2D-SED模型,以内蒙古准格尔旗沙圪堵镇附近的一个小流域为研究对象,以无人机航测的1m分辨率DEM数据、野外实测与室内实验获得的土壤特性数据、土地利用数据和降雨数据为基础,采用3种水沙模拟方案进行多象元尺度的水沙过程模拟,进而探索水沙过程模拟的DEM尺度效应及发生机制。研究表明:⑴在4~20m GRID分辨率区间模拟的径流量位于323.18m3和411.43m3之间,波动不大;⑵2~20m GRID分辨率区间内,模拟的侵蚀流量在3.43m3和65.61m3间变化,波动很大;(3)坡度和径流路径是水文过程模拟DEM尺度效应的两个对立影响因子,是水文过程模拟DEM尺度效应不明显的主要原因;⑷DEM尺度效应对侵蚀输沙具有重要影响,地形坡度是侵蚀输沙DEM尺度效应的主要控制因子;⑸地形坡度随DEM分辨率降低而发生的空间上的波动变化是侵蚀输沙量随DEM分辨率降低而波动变化的原因。     

基于钻孔数据任意边域地层分界面B-Rep建模研究

对相邻钻孔数据的类型进行了分析,对不同连续性和规则性的地层分界面进行了建模.考虑地层分界面的光滑性和整体控制性,运用边界表示法(B-Rep)和B样条曲面重构地层分界面,为地层建模提供了一种新的方法.建模研究表明:在钻孔数据前处理算法的基础上,B-Rep能够建立连续型、非连续型、空洞型等任意边域的地层分界面.     

数值流形法的T样条局部加密算法

在岩土工程分析中求解精度控制常常是必需的,在数值流形法中可以通过控制数学覆盖网格的稀疏和覆盖位移的阶数来达到精度的要求。提出了基于等几何分析的数值流形方法,定义了相应的数学覆盖的构造形式,推导了基于二次B样条的9节点数值流形方法分析格式;针对基于Lagrange插值函数的4节点数值流形方法提出了基于T样条思想的数学覆盖网格的局部加密方法。算例计算结果表明,相对于4节点的数值流形方法,基于非均匀有理B样条的9节点数值流形方法具有更高的精度;基于T样条思想的加密网格在保持计算精度的前提下降低了自由度的数量,表明T样条加密是一种自然的局部加密算法。     

基于多层DEM与QTPV的混合数据模型及其在地质建模中的应用

提出一种基于多层DEM与似三棱柱体元（QTPV）的混合三维数据模型。模型包含顶点、线段（棱边、三角形边）、三角形、侧面四边形、似三棱柱体元和DEM 6个基本元素,同时包含点对象、线对象、面对象、体对象、复杂对象、空间对象等6个对象。设计了6个基本元素和2种地质对象的数据结构和它们之间的拓扑关系。研究了基于原始采样点和基于内插点的两种建模方法。以模拟数据和地质钻探数据对开发的原型系统进行验证。研究表明,基于多层DEM与似三棱柱体元的混合三维数据模型具有同时表达空间对象的表面和内部结构的能力,适合地质勘探领域的三维建模。     

基于矩形特征的无空间参考DEM图像模糊匹配

多源、多时相与多尺度DEM的匹配是DEM(Digital Elevation Model)应用中的关键技术。对无空间参考的DEM进行匹配时,若特征点不足或匹配DEM之间的分辨率、地形特征等差异较大,已有的DEM匹配技术往往难以成功。为此,这里提出了基于矩形特征的模糊匹配方法。首先利用矩形特征构建了五种图像相似度指标BRFSI并统计分析了各个指标的分类性能;然后,使用Gentle AdaBoost算法对大量样本进行训练后,得到了用于判定"匹配"与"非匹配"的分类器;最后建立了进行DEM图像模糊匹配的匹配模型。经过四种典型数据验证,该方法很好地解决了无空间参考DEM匹配中特征点不足、DEM之间相似度较低的难点问题,为进一步的DEM精匹配提供了良好的初始条件。     

生态地理建模中的多尺度问题

本文在分析生态地理建模内涵的基础上,讨论了生态地理建模中的尺度转换问题、跨尺度相互作用问题、空间尺度与时间尺度的关联问题和多尺度数据处理问题.由于生态地理问题的非线性、生态环境的异质性和随机事件,简单的线性尺度转换方法远不能满足生态地理建模的要求.为了从根本上解决生态地理建模中的时空尺度问题,除需要运用微分几何学和等级理论等经典方法外,还需要引入格点生成法和网格计算等现代理论和技术手段.     

二维随机介质对声波传播影响

地层的形成是随机性和不确定性,用基于随机理论的方法来描述更为适合。利用随机介质建模的方法来描述地下介质中的小尺度非均匀性。为了考察小尺度非均匀性介质中声波的传播特性,建立了不同尺度的二维随机介质模型,利用有限差分法模拟随机介质中的声波场。分别对比不同自相关长度条件下地震记录及其波形与振幅谱,研究这些特征量与随机介质模型的自相关长度的关系。结果表明,自相关长度作为描述小尺度非均匀性尺度的参量,对声波在随机介质中传播影响很大。     

土壤水分运动空间变异性尺度效应的染色示踪入渗试验研究

探讨了土壤水非均匀流动特性和描述方法,通过染色示踪剂调查了三种试验尺度条件下非均匀流动模式,并采用随机层叠模型对不同实验尺度条件下非均匀流动模式进行了模拟。随机层叠模型中具有对数正态分布性质的随机层叠发生器被用来描述水流入渗过程,不同的方法被用于模型参数求解。试验观测和模拟计算结果均表明,尺度特性是非均匀流动的重要影响因素之一,准确的描述不同研究尺度下的非均匀流动特征,须同时考虑流动在水平和垂直方向的变异性。随着研究尺度的增加,流动的非均匀性变异程度更加明显。     

空间分辨率与取样方式对DEM流域特征提取的影响

随着数字水文的兴起和分布式水文模型研究的发展, 利用DEM提取水文特征, 进而进行水文模拟的方法越来越广泛地为水文学者所采用. 空间分辨率的改变与DEM重新取样方式对水文模拟都会产生重要影响. 采取不同取样方法获得多种尺度的DEM, 对不同分辨率下的流域特征值进行了统计分析与比较, 引入熵的概念度量不同分辨率的DEM包含的信息量, 以及不同取样方式对信息量的影响. 并计算了以50 m DEM所包含的信息量为基准, 在不同的信息损失下所要求的最低分辨率.     

Spatial probabilistic modeling of slope failure using an integrated GIS Monte Carlo simulation approach

Spatial probabilistic modeling of slope failure using a combined Geographic Information System (GIS), infinite-slope stability model and Monte Carlo simulation approach is proposed and applied in the landslide-prone area of Sasebo city, southern Japan. A digital elevation model (DEM) for the study area has been created at a scale of 1/2500. Calculated results of slope angle and slope aspect derived from the DEM are discussed. Through the spatial interpolation of the identified stream network, the thickness distribution of the colluvium above Tertiary strata is determined with precision. Finally, by integrating an infinite-slope stability model and Monte Carlo simulation with GIS, and applying spatial processing, a slope failure probability distribution map is obtained for the case of both low and high water levels.     

Accuracy Analysis of Digital Elevation Model Relating to Spatial Resolution and Terrain Slope by Bilinear Interpolation

This study investigates the accuracy analysis of the digital elevation model (DEM) with respect to the following two major factors that strongly affect the interpolated accuracy: (1) spatial resolution of a DEM and (2) terrain slope. Unlike existing studies based mainly on a simulation approach, this research first provides an analytical approach in order to build the relationship between the interpolated DEM accuracy and its influencing factors. The bi-linear interpolation model was adopted to produce this analytic model formalized as inequalities. Then, our analytic models were verified and further rectified by means of experimental studies in order to derive a practical formula for estimating the DEM accuracy together with an optimization model for calculating the required resolution when a prescribed upper bound to the DEM accuracy is given. Moreover, this analytic approach can cope with either a grid-based DEM or a randomly scattered scenario whose efficacies have been validated by the experiments using both synthetic and realistic data sets. In particular, these findings first establish the rules for directly correlating the horizontal resolution of DEM data with vertical accuracy.     

Development of an urban landslide cellular automata model: a case study of North Vancouver, Canada

Many GIS-based landslide models require detailed datasets that are ideally collected from field measurements, which can incur high costs for carrying out surveys. Even when the data is on hand, implementing physics-based slope stability techniques can be difficult. Common research practice uses differential equations to characterize the dynamic flow of a landslide, but it is often laborious without making substantial simplifications. A possible solution is to implement a cellular automata modeling approach, which represents both spatial and temporal components, to simulate the dynamics of the landslide propagation process. In this study, a simplified cellular automata model is developed for the effective prediction of landslide runouts, where the data requirement is a high resolution digital elevation model (DEM). Parameters, such as slope and slope curvature features, are derived from the DEM and coupled with logistic regression. The developed model is implemented on the Patrick and Dawson-Chu Slide in North Vancouver, Canada. The results from this study site were favorable, given almost 90% agreement between simulated landslides and data obtained for real landslides. In addition, sensitivity analysis was performed on the initiation sites to test the model logic and outputs of the landslide flow.     

数字高程模型分辨率对流域地形特征参数的影响

地形特征(如高程和坡度)和水文特征(如河流长度和河流坡度)是分布式流域水文水质模型的基础输入参数,用于量化描述模型模拟流域的自然特征。这些特征参数的准确性直接影响水文水质过程模拟的准确性。应用数字高程模型(Digital Elevation Model,DEM)在4个不同地形的子流域研究了10种不同分辨率DEM对平均高程、流域面积、坡度、河流坡度、最长河长等参数的影响。结果表明,随着DEM分辨率降低,流域地形变缓,流域平均坡度逐渐减小;随着DEM网格分辨率的变化,子流域划分范围和河道位置也都可能发生变化,且该变化在地形起伏较小的丘陵平原地区较明显,子流域集水面积和河长进一步随之改变;河流坡度随DEM分辨率降低则呈无规则变化。从地形和水文参数两方面揭示了DEM 分辨率在分布式流域模型中的不确定性影响。     

Application of sky view factor technique to the interpretation and reactivation assessment of landslide activity

High-resolution digital elevation models are crucial to the investigation of natural disasters, and a variety of methods based on visualization and relief map compilations have been proposed. In this study, the sky view factor (SVF) is applied to slope maps and a digital elevation model (DEM) of the Oso landslide, a deadly landslide that occurred in Washington State on March 22, 2014, to demonstrate the effectiveness of SVF-enhanced relief maps in mapping and evaluating large-scale or deep-seated landslide hazards. A procedure for combining the SVF-enhanced DEM with slope and elevation maps is also presented. Then the maps are used to extract the landslide-prone areas and perform a reactivation analysis of the post-Oso landslide using an analytic hierarchy process (AHP). By using the SVF-enhanced DEM to perform the AHP assessment on multi-period images, we accurately evaluate hazard of the landslide for both pre and post-2014 conditions. Finally, different visualization maps, limitation and recommend parameters for generating SVF relief map are presented in the paper.     

基于GIS的福建省地形分析

数字高程模型（DEM）是DTM的一种特例,是对地球表面地形地貌的一种离散的数字表达。DEM是地理信息系统进行地形分析的基础数据。采用1：10万数字等高线数据建构30m空间分辨率格网DEM,并基于DEM进行分级,提取坡度、地形起伏度及坡度聚类,坡谱等地形分析;比较准确地表达了福建的总体地形特征,为地理分析及环境评价等提供了基础的数据。     

Comparison of different algorithms for calculating the shading effects of topography on solar irradiance in a mountainous area

In mountainous areas, topographic shading effects exert significant influence on the spatial distribution of solar irradiance on the land surface, which should be considered in estimating and modeling physical, biological, and hydrological processes. Various algorithms have been proposed for quantitative assessment of the impact of topography on the arriving solar radiation on slopes in mountainous area based on digital elevation model (DEM). However, different algorithms tend to produce different results even based on the same DEM. In this study, two algorithms for determining the maximum horizontal viewing angle of terrain obstruction were compared in a catchment in the Qilian Mountains, Northwest China, which were identical except for the way in interpreting DEM. Two algorithms for calculating sky view factor were also compared. Then, different algorithm combinations were compared and their sensitivities to the horizontal resolution of DEM were investigated. Results showed that the differences in outputs of different algorithms and combinations were significant and cannot be ignored in estimating the spatial distribution of solar irradiance in mountainous areas. A simple ground validation concluded that the algorithm assuming a linear variation in elevation was a little better than that assuming a stepped variation in elevation in and among neighboring DEM cells. For calculating diffuse solar radiation, the algorithm based on a sloped coordinate system was theoretically better than that based on a horizontal coordinate system. It is worth noting that even though the estimated annual mean of solar irradiance over the study area was almost the same among some algorithm combinations, the difference in spatial distribution was significant, which may lead to consequent difference in assessing or modeling some other environmental variables relevant to surface energy exchange, e.g., evaporation, snowmelt, and ground temperature. However, the analysis in this study was mainly theoretical, and further ground-based validations should be made in future.     

The effectiveness of high-resolution LiDAR data combined with PSInSAR data in landslide study

The spatial resolution of digital elevation models (DEMs) is an important factor for reliable landslide studies. Multi-interferometric techniques such as persistent scatterer interferometric synthetic aperture radar (PSInSAR) are used to evaluate the landslide state of activity and its ground deformation velocity, which is commonly measured along the satellite line of sight (LOS). In order to compare velocities measured by different satellites in different periods, their values can be projected along the steepest slope direction, which is the most probable direction of real movement. In order to achieve this result, DEM-derived products are needed. In this paper, the effectiveness of different DEM resolutions was evaluated in order to project ground deformation velocities measured by means of PSInSAR technique in two different case studies in the Messina Province (Sicily, southern Italy): San Fratello and Giampilieri. Three DEMs were used: (i) a 20-m resolution DEM of the Italian Military Geographic Institute (IGM), (ii) a 2-m resolution DEM derived from airborne laser scanning (ALS) light detection and ranging (LiDAR) data for the San Fratello 2010 landslide, and (iii) a 1-m resolution DEM derived from ALS LiDAR data for the area of Giampilieri. The evaluation of the applied method effectiveness was performed by comparing the DEMs elevation with those of each single permanent scatterer (PS) and projecting the measured velocities along the steepest slope direction. Results highlight that the higher DEM resolution is more suitable for this type of analysis; in particular, the PS located nearby the watershed divides is affected by geometrical problems when their velocities are projected along the steepest slope.