首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 203 毫秒
1.
With the increasing sizes of digital elevation models (DEMs), there is a growing need to design parallel schemes for existing sequential algorithms that identify and fill depressions in raster DEMs. The Priority-Flood algorithm is the fastest sequential algorithm in the literature for depression identification and filling of raster DEMs, but it has had no parallel implementation since it was proposed approximately a decade ago. A parallel Priority-Flood algorithm based on the fastest sequential variant is proposed in this study. The algorithm partitions a DEM into stripes, processes each stripe using the sequential variant in many rounds, and progressively identifies more slope cells that are misidentified as depression cells in previous rounds. Both Open Multi-Processing (OpenMP)- and Message Passing Interface (MPI)-based implementations are presented. The speed-up ratios of the OpenMP-based implementation over the sequential algorithm are greater than four for all tested DEMs with eight computing threads. The mean speed-up ratio of our MPI-based implementation is greater than eight over TauDEM, which is a widely used MPI-based library for hydrologic information extraction. The speed-up ratios of our MPI-based implementation generally become larger with more computing nodes. This study shows that the Priority-Flood algorithm can be implemented in parallel, which makes it an ideal algorithm for depression identification and filling on both single computers and computer clusters.  相似文献   

2.
由于DEM数据不包括河流、湖泊(水库)及流域边界和堤坝等信息,因此常规仅依靠DEM提取水系的方法不能反映平坦区域及受人类影响强烈区域水系的真实特征。该文建立了融合DEM、河流、湖泊水库、流域边界和堤坝等多源数字信息的流域水系提取方法,对多源信息进行栅格化,采用所提出的高程-距离函数对DEM进行校正,使得提取的数字水系与实际水系精确拟合。将该方法应用于山区太湖西苕溪流域老石坎水库及其上游集水区和平原丘陵分布区的淮河史灌河蒋家集站、梅山、鲇鱼山水库区间流域,提取的河网结构与地图信息拟合较好。  相似文献   

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

4.
Third-order geomorphometric variables (based on third derivatives of the altitudinal field) have been neglected in geomorphometry, but their application to the delimitation of surface objects will lead to their increasing significance in future. New techniques of computation, presented and evaluated here, facilitate their use. This paper summarizes recent knowledge concerning definition, computation and geomorphologic interpretation of these variables. Formulae defining various third-order variables are unified based on the physical definition of slope gradient. Methods for their computation are compared from the point of view of method error and error generated by digital elevation model (DEM) inaccuracy. For exact mathematical test surfaces, the most natural and simple variant of the method of central differences (CD2) shows a method error 2–3 times smaller than the other methods used recently in geomorphometry. However, success in coping with DEM inaccuracy depends (for a given grid mesh) on the number and weighting of points from which the derivative is computed. This was tested for surfaces with varying degrees of random error. Here least squares-based methods are the most effective for mixed derivatives (especially for finer grids and less accurate DEMs), while a variant of the CD method, that repeats numerical evaluation of first derivatives (CD1), is the most successful for derivatives in cardinal directions. The CD2 method is generally the most successful for coarser grids where the method error is dominant.

Utilization of third-order variables is documented from examples of terrain feature (ridge, valley and edge) extraction and from a first statistical test of the hypothesis that real segments of the land surface have a tendency to a constant value of some morphometric variable. For detection of (sharp) ridges and valleys, it is shown that the rate of change of tangential curvature is inadequate: rate of change of normal curvature is also required. A basic confirmation of the constant-value tendency is provided.  相似文献   

5.
As increasingly large‐scale and higher‐resolution terrain data have become available, for example air‐form and space‐borne sensors, the volume of these datasets reveals scalability problems with existing GIS algorithms. To address this problem, a kind of serial algorithm was developed to generate viewshed on large grid‐based digital elevation models (DEMs). We first divided the whole DEM into rectangular blocks in row and column directions (called block partitioning), then processed these blocks with four axes followed by four sectors sequentially. When processing the particular block, we adopted the ‘reference plane’ algorithm to calculate the visibility of the target point on the block, and adjusted the calculation sequence according to the different spatial relationships between the block and the viewpoint since the viewpoint is not always inside the DEM. By adopting the ‘Reference Plane’ algorithm and using a block partitioning method to segment and load the DEM dynamically, it is possible to generate viewshed efficiently in PC‐based environments. Experiments showed that the divided block should be dynamically loaded whole into computer main memory when partitioning, and the suggested approach retains the accuracy of the reference plane algorithm and has near linear compute complexity.  相似文献   

6.
Digital elevation models (DEMs) given by spheroidal trapezoidal grids are more appropriate for large regional, sub-continental, continental and global geological and soil studies than square-spaced DEMs. Here we develop a method for derivation of topographic variables, specifically horizontal (k) and vertical h (k) landsurface curvatures, from spheroidal trapezoidal-spaced DEMs. First, we v derive equations for calculation of partial derivatives of elevation with DEMs of this sort. Second, we produce formulae for estimation of the method accuracy in terms of root mean square errors of partial derivatives of elevation, as well as k h and k (m and m respectively). We design the method for the case that the v kh k v Earth's shape can be ignored, that is, for DEM grid sizes of no more than 225 km. We test the method by the example of fault recognition using a DEM of a part of Central Eurasia. A comparative analysis of test results and factual geological data demonstrates that the method actually works in regions marked by complicated topographic and tectonic conditions. Upon increasing DEM grid size, one can produce generalised maps of k and k. Spatial distributions of m and m h v kh k v depend directly on the distribution of elevation RMSE. Areas with high values of m are marked by low values of m, and vice versa, areas with high values kh k v of m are marked by low values of m. Data on m and m should be utilised k v kh kh k v to control and improve applications of k and k to geological studies. The method h v developed opens up new avenues for carrying out some 'conventional' raster operations directly on geographical co-ordinates.  相似文献   

7.
Digital elevation model (DEM) elevation accuracy and spatial resolution are typically considered before a given DEM is used for the assessment of coastal flooding, sea-level rise or erosion risk. However, limitations of DEMs arising from their original data source can often be overlooked during DEM selection. Global elevation error statistics provided by DEM data suppliers can provide a useful indicator of actual DEM error, but these statistics can understate elevation errors occurring outside of idealised ground reference areas. The characteristic limitations of a range of DEM sources that may be used for the assessment of coastal inundation and erosion risk are tested using high-resolution photogrammetric, low- and medium-resolution global positioning system (GPS)-derived and very high-resolution terrestrial laser scanning point data sets. Errors detected in a high-resolution photogrammetric DEM are found to be substantially beyond quoted error, demonstrating the degree to which quoted DEM accuracy can understate local DEM error and highlighting the extent to which spatial resolution can fail to provide a reliable indicator of DEM accuracy. Superior accuracies and inundation prediction results are achieved based on much lower-resolution GPS points confirming conclusions drawn in the case of the photogrammetric DEM data. This suggests a scope for the use of GPS-derived DEMs in preference to the photogrammetric DEM data in large-scale risk-mapping studies. DEM accuracies and superior representation of micro-topography achieved using high-resolution terrestrial laser scan data confirm its advantages for the prediction of subtle inundation and erosion risk. However, the requirement for data fusion of GPS to remove ground-vegetation error highlighted limitations for the use of side-scan laser scan data in densely vegetated areas.  相似文献   

8.
典型黄土地貌类型区的地形复杂度分形研究   总被引:9,自引:0,他引:9  
龙毅  周侗  汤国安  刘学军 《山地学报》2007,25(4):385-392
黄土高原地貌形态与地形复杂度自南向北有序变化,构成了举世瞩目的独特的地理景观。选择陕西省南北剖面六个典型黄土地貌样区为基本实验区,以其1∶1万栅格DEM为数据源,探讨典型黄土地貌类型区的地形复杂度分形与空间分异特征。首先提出的元分维模型方法,以计算得到的DEM元分维值作为特征指标,研究样区的地形复杂度问题。这种通过滑动窗口的扩展分维分析方法,既可以用于分析该区域不同尺度下的地形复杂度变化情况,也可以探讨区域的局部单元复杂程度及其空间分布,从而不再局限于对全区域的单一分维评价。以此为基础,进一步应用元分维谱方法,研究地貌网格单元的元分维分级分布情况。实验结果表明:以绥德和延川为代表的黄土峁状丘陵沟壑区最为复杂,以宜君和甘泉为代表的梁状丘陵沟壑区居中,而以淳化为代表的黄土塬区和以神木为代表的风沙黄土过渡区最为平缓。实验进一步证明了扩展分形方法在黄土地貌研究中的可行性。  相似文献   

9.
数字高程模型(DEM)在表达地貌形态、认知地表过程、揭示地学机理等研究中发挥着基础性的作用,是重要的地理空间数据模型,广泛地应用于地学分析与建模中。但是,传统DEM具有属性单一的天然缺陷,难以支撑面向地学过程与机理挖掘的地球系统科学研究。亟待在传统DEM的基础上实现其数据模型的增值,服务于新地貌学研究范式和新对地观测技术背景下的数字地形建模与分析。立足于以上问题,本文构建了DEM增值的理论框架,主要包括DEM增值的概念、内涵、内容、类别、不同增值类别之间的相互关系,以及此理论框架的研究意义和应用范畴。提出了DEM增值的构建方法,包含:① 强调地上地下一体化、时间空间相耦合的DEM空间维度和时间维度增值方法;② 重视地下、地表和地上物质构成,形态属性耦合的物质属性和形态属性增值方法;③ 顾及自然过程、人工作用的地物对象、地貌形态的地物要素和形体要素增值方法。最后,分别以数字阶地模型、数字坡地模型和数字流域模型为例,阐释DEM在面向地貌学本源问题时的不同增值方法及应用场景。期望通过对DEM进行维度、属性和要素3个层面的增值,实现现代对地观测技术背景下数字高程模型表达方法的突破,并支撑知识驱动的数字地貌问题分析。  相似文献   

10.
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.  相似文献   

11.
As sea level is projected to rise throughout the twenty-first century due to climate change, there is a need to ensure that sea level rise (SLR) models accurately and defensibly represent future flood inundation levels to allow for effective coastal zone management. Digital elevation models (DEMs) are integral to SLR modelling, but are subject to error, including in their vertical resolution. Error in DEMs leads to uncertainty in the output of SLR inundation models, which if not considered, may result in poor coastal management decisions. However, DEM error is not usually described in detail by DEM suppliers; commonly only the RMSE is reported. This research explores the impact of stated vertical error in delineating zones of inundation in two locations along the Devon, United Kingdom, coastline (Exe and Otter Estuaries). We explore the consequences of needing to make assumptions about the distribution of error in the absence of detailed error data using a 1 m, publically available composite DEM with a maximum RMSE of 0.15 m, typical of recent LiDAR-derived DEMs. We compare uncertainty using two methods (i) the NOAA inundation uncertainty mapping method which assumes a normal distribution of error and (ii) a hydrologically correct bathtub method where the DEM is uniformly perturbed between the upper and lower bounds of a 95% linear error in 500 Monte Carlo Simulations (HBM+MCS). The NOAA method produced a broader zone of uncertainty (an increase of 134.9% on the HBM+MCS method), which is particularly evident in the flatter topography of the upper estuaries. The HBM+MCS method generates a narrower band of uncertainty for these flatter areas, but very similar extents where shorelines are steeper. The differences in inundation extents produced by the methods relate to a number of underpinning assumptions, and particularly, how the stated RMSE is interpreted and used to represent error in a practical sense. Unlike the NOAA method, the HBM+MCS model is computationally intensive, depending on the areas under consideration and the number of iterations. We therefore used the HBM+ MCS method to derive a regression relationship between elevation and inundation probability for the Exe Estuary. We then apply this to the adjacent Otter Estuary and show that it can defensibly reproduce zones of inundation uncertainty, avoiding the computationally intensive step of the HBM+MCS. The equation-derived zone of uncertainty was 112.1% larger than the HBM+MCS method, compared to the NOAA method which produced an uncertain area 423.9% larger. Each approach has advantages and disadvantages and requires value judgements to be made. Their use underscores the need for transparency in assumptions and communications of outputs. We urge DEM publishers to move beyond provision of a generalised RMSE and provide more detailed estimates of spatial error and complete metadata, including locations of ground control points and associated land cover.  相似文献   

12.
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.  相似文献   

13.
地形湿度指数算法误差的定量评价   总被引:2,自引:0,他引:2  
地形湿度指数(TWI)能够定量指示地形对土壤湿度空间分布的控制,是一种应用广泛的地形属性.目前基于栅格DEM的TWI计算方法结果各异,因此有必要对'TWI算法进行定量评价.对TWI算法通常是应用实际DEM数据进行评价.但实际DEM中存在的数据源误差会干扰对算法误差的评价.针对该问题,本文介绍了一种用不含数据源误差的人造...  相似文献   

14.
The smoothness of HASM   总被引:1,自引:0,他引:1  
To smooth noises inherent in uniformly sampled dataset, the smoothness of high accuracy surface modeling (HASM) was explored, and a smoothing method of HASM (HASM-SM) was developed based on a penalized least squares method. The optimal smoothing parameter of HASM-SM was automatically obtained by means of the generalized cross-validation (GCV) method. For an efficient smoothing computation, discrete cosine transform was employed to solve the system of HASM-SM and to estimate the minimum GCV score, simultaneously. Two examples including a numerical test and a real-world example were employed to compare the smoothing ability of HASM-SM with that of GCV thin plate smoothing spline (TPS) and kriging. The numerical test indicated that the minimum GCV HASM-SM is averagely more accurate than TPS and kriging for noisy surface smoothing. The real-world example of smoothing a lidar-derived Digital Elevation Model (DEM) showed that HASM-SM has an obvious smoothing effect, which is on a par with TPS. In conclusion, HASM-SM provides an efficient tool for filtering noises in grid-based surfaces like remote sensing–derived images and DEMs.  相似文献   

15.
This paper proposes a quantitative method to classify landforms using four morphometric parameters from DEM-derived thematic raster maps of slope and topographic openness. Because the different surficial processes and stages in the evolution of slopes create landscapes with different shapes, these parameters may lead to a genetic interpretation of topography. The raster maps of slope and topographic openness were constructed for Northeast Honshu, Japan, from 50-m DEMs. The mean and standard deviation of morphometric parameters within a 3050 m by 3050 m moving window on the raster maps were calculated. The results for some training areas show that constructional/depositional and erosional landforms with different relief have different morphometric characteristics. A supervised landform classification for Northeast Honshu using the knowledge from the training areas revealed a ladder geomorphological structure composed of high mountains, ranges and volcanoes. The close relationship between the ladder geomorphological structure and volcano distribution indicates that the structure reflects the magmatic plumbing system from the upper mantle to the crust of the Northeast Honshu arc.  相似文献   

16.
The detailed topographic information contained in light detection and ranging (LiDAR) digital elevation models (DEMs) can present significant challenges for modelling surface drainage patterns. These data frequently represent anthropogenic infrastructure, such as road embankments and drainage ditches. While LiDAR DEMs can improve estimates of catchment boundaries and surface flow paths, modelling efforts are often confounded difficulties associated with incomplete representation of infrastructure. The inability of DEMs to represent embankment underpasses (e.g. bridges, culverts) and the problems with existing automated techniques for dealing with these problems can lead to unsatisfactory results. This is often dealt with by manually modifying LiDAR DEMs to incorporate the effects of embankment underpasses. This paper presents a new DEM pre-processing algorithm for removing the artefact dams created by infrastructure in sites of embankment underpasses as well as enforcing flow along drainage ditches. The application of the new algorithm to a large LiDAR DEM of a site in Southwestern Ontario, Canada, demonstrated that the least-cost breaching method used by the algorithm could reliably enforce drainage pathways while minimizing the impact to the original DEM.  相似文献   

17.
Airborne laser scanning (ALS) is well suited for the production of digital elevation models (DEM), and can, in contrast to photographic methods, be used to acquire a DEM independently of surface texture and external light sources. ALS thus serves as a tool to generate DEMs of firn areas where photogram- metric methods often fail.
The potential of an integrated ALS system – comprising a laser scanner, precise differential global positioning system, and a gyro platform – for DEM generation of firn areas is currently being assessed. The Unteraargletscher, Bernese Alps, Switzerland, has been chosen as a test site. As part of a pilot project aimed at determining the mass balance distribution of that glacier without the use of in situ information, ALS measurements were conducted in autumn 1997. The DEM derived from laser measurements is extremely sensitive to the position and attitude of the aircraft. Currently the main work focuses on assessing and improving the system's accuracy by error modelling and by the development of error-correction algorithms.
Preliminary results from Unteraargletscher are presented, and the potential of this method for the generation of DEMs of firn areas is discussed.  相似文献   

18.
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.  相似文献   

19.
20.
基于DEM的地貌实体单元自动提取方法   总被引:17,自引:4,他引:13  
我国传统地貌基本形态类型分类强调地貌单元的完整性,界线划分沿地貌实体边界而非规则统计单元,目前尚缺乏地貌实体单元的有效自动提取方法。针对这一难点,本文提出一种基于DEM的地貌实体单元数字提取方法。利用坡度分级,并搜索相邻栅格单元、计算坡度级别内相互连通栅格的面积,建立坡度、面积阈值综合判别规则进行山地平原的自动划分;利用地形倒置、水文淹没分析,将山体划分的二维判别规则扩展到实际三维地形中,并结合地形结构线提取算法进行山体界线自动提取、确定山地地貌实体单元。结果表明,该方法符合我国传统地貌分类体系,能够较好实现山地/平原的自动划分和山体界线的数字提取。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号