基于数字高程模型的混合流向算法

从数字高程模型提取的汇水网络和汇水区等信息是分布式水文模型及应用分析的基础参数,基于地表汇水模拟的算法是提取该类信息的主要方法,其中,水流方向的确定对提取结果有着直接的影响。单流向算法因其易于实现、易于确定上游汇水区等特性,得到了广泛应用,然而单流向算法在坡度平缓区域会产生不自然的平行径流,能模拟地表水流分散径流特点的多流向算法可以在一定程度上避免此问题,但多流向算法使得不同区域的汇水单元可能存在交叉。本文结合两类流向算法各自的优点和适用性,设计实现了一种混合流向算法,以期在不同的地形条件下模拟得到更加合理的水流分配。首先,使用基于模板的形态检测方法,在给定阈值的基础上,对数字地形进行了分类,DEM被划分为山谷、山脊、鞍部、缓坡和陡坡5类。对陡坡、山谷和山脊区域运用单流向算法;对缓坡和鞍部区域采用多流向算法确定径流方向并进行水量分配。本文选取了黄土地貌和中低山丘陵的两个流域作为研究区,利用并采用了30 m和90 m两个分辨率的DEM。本文研究将混合流向算法与现有其他算法的结果进行比较。相比于多流向算法,该算法结果中的分散效应受到明显的抑制,相比于单流向算法,非自然的平行径流也大幅减少。同时,混合流向算法在较大分辨率DEM上（30 m）改进效果更加明显。     

Polarimetric SAR Data Correction and Terrain Topography Measurement Based on the Radar Target Orientation Angle

Topographic variations caused by the range and the azimuth terrain slopes induce polarization orientation changes which cause the polarization to rotate about the line of sight. The existence of these variations reduce the accuracy measurement of geophysical parameters from polarimetric synthetic aperture radar (PolSAR) images. For this reason most inversion studies are best done in area of flat earth. In area which has significant terrain variations require compensation for topography. In real situations, terrain slopes rotate the polarization basis of the polarimetric scattering matrices by an orientation angle shift, and induce significant cross-polarization power. In this paper, two methods have been investigated using the polarimetric orientation angle (PAO): the first one involves the rotation of the polarimetric scattering and coherency matrices to achieve maximum azimuthally asymmetry for polarimetric data compensation to ensure accurate estimation of geophysical parameters in rugged terrain areas. The second approach has been developed which measures azimuth and range terrain slopes that are related to shifts in polarization orientation angle. Terrain elevation maps relative to a plane parallel to the radar line of sight can then be generated by integrating these slopes requiring only one PolSAR flight pass by combing orientation angle estimation and a shape-from-shading technique (SFS) which is mostly used by the computer vision community. Experimental results with C-band polarimetric RADARSAT2 data are used evaluate the data compensation algorithm and DEM generation.     

Error assessment of grid‐based terrain shading algorithms for solar radiation modeling over complex terrain

In mountainous regions, solar radiation exhibits a strong spatial heterogeneity due to terrain shading effects. Terrain shading algorithms based on digital elevation models can be categorized into two types: area‐based and point‐specific. In this article, we evaluated two shading algorithms using designed mathematic surfaces. Theoretical shading effects over four Gauss synthetic surfaces were calculated and used to evaluate the terrain shading algorithms. We evaluated the area‐based terrain shading algorithm, Hillshade tool of ArcGIS, and the point‐specific shading algorithm from Solar Analyst (SA) in ArcGIS. Both algorithms showed shading overestimation, and Hillshade showed more accuracy with a mean absolute error (MAE) of 1.20%, as compared to the MAE of 1.66% of SA. The MAE of Hillshade increases exponentially as the spatial extent of the study area increases because the solar position for all locations on the surface is the same in Hillshade. Consequently, we suggest that the surface should be divided into more tiles in Hillshade when the discrepancy in the latitude of the whole surface is greater than 4°. Skyshed, which represents the horizon angle distribution in SA, is error‐prone over more complex terrain because horizon angle interpolation is problematic for such areas. We also propose a new terrain shading algorithm, with solar positions calculated using local latitude for each cell and the horizon angle calculated for every specific time interval, but without projections. The new model performs better than Hillshade and SA with an MAE of 0.55%.     

基于DEM坡度坡向算法精度的分析研究

坡度坡向是两个最基本的地形因子,目前对DEM坡度坡向计算模型和精度存在一些不同的甚至矛盾的观点,其原因在于没有区分误差来源和分析评价方法的不同.本文对DEM坡度坡向误差进行了理论分析,并通过实验数据对相关结论进行了验证.旨在澄清目前关于坡度坡向计算模型上的矛盾结论.     

基于DEM的地形曲率计算模型误差分析

地形曲率是地形表面几何形态和地学建模的基本变量之一。本文首先对曲率计算模型进行了归纳,然后通过数据独立的DEM误差分析方法和实际DEM的分析验证,对目前九种曲率计算的三类曲率计算模型进行了量化分析比较。研究结果表明,当DEM数据精度比较高时,高次曲面(四次曲面)能给出较高精度的曲率计算结果,而当DEM误差较大时,低次曲面(二次曲面)由于具有误差的平滑作用而能产生较高精度的曲率值。     

大钝角剖分与最小二乘法约束迭代算法优化TIN模型研究

TIN模型通过从不规则分布的数据点生成的连续三角面来逼近地形表面,就表达地形信息的角度而言,TIN模型的优点是它能以不同层次的分辨率来描述地形表面。但在相关软件对离散数据生成TIN模型时会出现边缘失真现象,特别是当一些离散数据出现凹区域时,凹区域地形逼近失真现象相当明显。而在这些凹区域大多生成的三角网以钝角三角形居多,本文对离散数据TIN模型生成DEM时存在凹区域地形逼近失真现象进行了研究。提出一种大钝角剖分与最小二乘法约束迭代算法优化TIN模型。首先对原始TIN模型中大于100°小于1800的钝角三角形进行遍历记录,然后利用最小二乘法约束迭代算法进行优化,通过此算法优化后的TIN模型生成DEM数据,分别从小区域面积稳定性与剖面拟合实际地形进行了试验与精度对比,试验显示这种算法是可行的。     

Quality Assessment Of Digital Surface Models Generated From IKONOS Imagery

The growing applications of digital surface models (DSMs) for object detection, segmentation and representation of terrestrial landscapes have provided impetus for further automation of 3D spatial information extraction processes. While new technologies such as lidar are available for almost instant DSM generation, the use of stereoscopic high-resolution satellite imagery (HRSI), coupled with image matching, affords cost-effective measurement of surface topography over large coverage areas. This investigation explores the potential of IKONOS Geo stereo imagery for producing DSMs using an alternative sensor orientation model, namely bias-corrected rational polynomial coefficients (RPCs), and a hybrid image-matching algorithm. To serve both as a reference surface and a basis for comparison, a lidar DSM was employed in the Hobart testfield, a region of differing terrain types and slope. In order to take topographic variation within the modelled surface into account, the lidar strip was divided into separate sub-areas representing differing land cover types. It is shown that over topographically diverse areas, heighting accuracy to better than 3 pixels can be readily achieved. Results improve markedly in feature-rich open and relatively flat terrain, with sub-pixel accuracy being achieved at check points surveyed using the global positioning system (GPS). This assessment demonstrates that the outlook for DSM generation from HRSI is very promising.     

Simulation Study of a Robust Algorithm for Soil Moisture and Surface Roughness Estimation Using L‐Band Radar Backscatter

Abstract

Various inversion algorithms have been developed to obtain estimates of soil moisture and surface roughness parameters from multifrequency, multiangle, and multipolarization radar reflectances. Since the penetration depth for radar signals increases with wavelength, an inversion algorithm using widely separated frequencies does not yield comparable probing depths. Furthermore, existing algorithms assume a linear relationship between the radar backscatter coefficient (in dB) and soil parameters, such as the volumetric soil moisture, soil surface roughness and surface slope. This assumption is valid only over a narrow range of soil parameters, thereby restricting its operational use under realistic conditions. Our research specifically explored the use of inversion algorithms based on L‐Band radar reflectances at 1 GHz and 2 GHz frequencies in order to retain relatively consistent probing depths. In order to extend the range of applicability, a non‐linear exponential‐type relationship was developed between radar reflectance at a specified frequency, polarization and incidence angle combination, and soil parameters of interest, viz., soil moisture, surface roughness, and surface slope. An over‐constrained inversion algorithm using a six‐parameter combination was found to yield relatively accurate estimates of soil parameters over a wide range of soil conditions even in the presence of system error.     

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.     

Generalization of DEM for terrain analysis using a compound method

This paper reports an investigation into the generalization of a grid-based digital elevation model (DEM) for the purpose of terrain analysis. The focus is on the method of restructuring the grid-based surface elevation data to form a triangulated irregular network (TIN) that is optimized to keep the important terrain features and slope morphology with the minimum number of sample points. The critical points of the terrain surface are extracted from the DEM based on their significance, measured not only by their local relief, but also by their importance in identifying inherent geomorphological and drainage features in the DEM. A compound method is proposed by integrating the traditional point-additive and feature-point methods to construct a drainage-constrained TIN. The outcome is then compared with those derived from other selected methods including filtering, point-additive or feature-point algorithms. The results show that the compound approach is capable of taking advantage of both point-additive and feature-point algorithms to maximally keep the terrain features and to maintain RMSE at an acceptable level, while reducing the elevation data points by over 99%. The analytical result also shows that the proposed method outperforms the compared methods with better control in retaining drainage features at the same level of RMSE.     

Spatial Resolution and Algorithm Choice as Modifiers of Downslope Flow Computed from Digital Elevation Models

This research investigated the dependence of computing downslope flow from digital elevation models (DEMs) on two basic elements of the flow calculation: the spatial resolution of the DEM and the algorithm used to compute the downslope flow, specifically how it defines topographic slope and aspect. Six algorithms for downslope flow and downslope accumulation were implemented in MATLAB using different assumptions about components of the computation from the terrain analysis literature. We analyzed the results of downslope flow statistically, thresholded the values to yield streams, and compared the different results with the actual surface streams in the National Hydrographic Database. We repeated the computations using DEMs of 3-m, 30-m, and 90-m resolutions, covering a test area near Santa Barbara, California. We conclude that downslope flow computations are overestimates, and that they are fraught with critical algorithmic assumptions and scale effects. They should be approached with great caution in GIS-based analyses.     

利用质心Voronoi图对地形自适应简化的算法

地形简化算法利用少量有效的地形信息表达整体地形,能很好地解决海量地形数据与计算机硬件之间的矛盾,同时满足多尺度地形应用需求。针对现有地形简化算法难以兼顾局部地形起伏与地形整体特征的问题,提出一种基于质心Voronoi图的地形自适应简化算法。首先,利用质心Voronoi图的特点,以地形起伏度作为密度函数生成质心Voronoi图;然后,利用分布在地形起伏较大区域的质心Voronoi图种子点及大多分布在地形特征线上的Voronoi区域顶点重构地形;最后,通过原始地形与重构地形的特征线验证地形简化的效果,并与三维道格拉斯-普克（3D DouglasPeucker,3D DP）算法进行精度对比。实验结果表明,从简化地形中提取的山脊线、山谷线、等高线等地形特征线与原始地形的重叠度均较高,算法能较好地保持地形整体特征;且在相同的简化级别下,算法的简化误差小于3D DP算法,具有较高的地形简化精度。     

Development and Application of 2D and 3D GRASS Modules for Simulation of Thermally Driven Slope Winds

The ability to manage and process fully three‐dimensional information has only recently been made available for a few Geographical Information Systems (GIS). An example of integrated and complementary use of 2D and 3D GRASS modules for the evaluation and representation of thermally induced slope winds over complex terrain is presented. The analytic solution provided by Prandtl (1942) to evaluate wind velocity and (potential) temperature anomaly induced by either diurnal heating or nocturnal cooling on a constant angle slope is adopted to evaluate wind and temperature profiles at any point over both idealised and real complex terrain. As these quantities depend on the slope angle of the ground and on the distance from the slope surface suitable procedures are introduced to determine the coordinate n of a point in the 3D volume measured along the direction locally normal to the terrain surface. A new GRASS module has been developed to evaluate this quantity and to generate a 3D raster file where each cell is assigned the value of the cell on the surface belonging to the normal vector. The application of the algorithm implemented in GRASS to an ideal valley and to a real valley close to the city of Trento in the Alps provides results in accordance with data reported in the literature. An extension of Prandtl's (1942) model to take into account humidity and evaporation processes on the soil is also proposed and implemented.     

DEM Generation Combining SAR Polarimetry and Shape-From-Shading Techniques

Estimation of the polarization orientation angle shifts induced by terrain azimuth slope variations is a recently developed application in radar polarimetry. In general, without any prior knowledge on the terrain, two polarimetric SAR (POLSAR) flight passes are required to derive terrain slopes in perpendicular directions for digital elevation model (DEM) generation. Moreover, we note that SAR intensity is a strong indicator of the range component of the terrain slopes. In this letter, we developed a method for DEM generation requiring only one POLSAR flight pass, by combining orientation angle estimation and a shape-from-shading technique. In particular, when limited POLSAR data are available, this POLSAR technique provides an alternative way for DEM generation. National Aeronautics and Space Administration Jet Propulsion Laboratory (NASA/JPL) AIRSAR L-band POLSAR data over Camp Roberts, California, is used to demonstrate the results of the method proposed in this letter, and a DEM derived from simultaneously measured C-band interferometric SAR from NASA/JPL topographic SAR instrument is selected as the comparative ground truth to validate the effectiveness of this single POLSAR method. Analyses and discussions are also included in this letter.     

Phenomenon-based Specification of the Digital Representation of Terrain Surfaces

The phenomenon known as 'terrain' is a continuous surface. However, when a digital terrain representation is based on a regular raster (i.e. a DEM) the digital surface is commonly not continuous. This is the case for the derivation of variables such as slope, aspect, and curvature values as performed in today's Geographic Information Systems (GIS). Often, there is no surface specified at all, as, for instance, when flow lines or watersheds are constructed. The discrepancy between the phenomenon to be modelled and its digital representation causes the terrain analysis results to be less accurate than they could be. Furthermore, if more than one type of terrain information is derived the results are likely to be based on different specifications of the seemingly same terrain surface. The combined application of the derivation results will likely introduce inconsistencies. This paper suggests founding the specification of digital terrain representations on a careful analysis of the properties of the phenomenon. The paper details the reasons for, and advantages of, continuous surface representations and emphasises the importance of a comprehensive documentation of the conceptual models underlying digital terrain representations. A review of suitable interpolation approaches for the specification of terrain surfaces is given. The paper discusses how the resulting digital surfaces are analysed and how measurement uncertainty may be accounted for.     

地形特征提取的一种简易算法

地形特征作为基本地理要素,从DEM数据中提取其隐含的地形结构信息是非常重要的。本方法采用距离倒数加权法建立格网DEM,提出一种直接比较高程值来提取地形特征的简易算法,并利用差分算法计算坡度和坡向,通过V isual C++软件平台,编程实现格网DEM地形特征的自动提取,进行三维立体显示。该算法易于编程实现,数据结构简单,程序运行效率高,提取结果与真实地形基本拟合。     

格网DEM坡度计算模型的相似性研究

当DEM的精度和分辨率确定后,坡度坡向计算结果取决于计算模型。本文设计了逐栅格单元的相对差等级评价方法,并结合基于整体的频率图和XY散点图评价指标,定量揭示了坡度计算模型之间的异同。主要结论为:①坡度计算模型有离散和连续之分,连续型模型能给出更加符合实际的坡度坡向计算结果;②坡度计算结果与计算点数有关,点数越多,计算结果越稳定;③算法间的差异随DEM分辨率的增大而增大。最后本文给出了实际应用中的坡度计算模型。     

规则格网DEM化简的一种改进方法

分析了从较高分辨率DEM化简到较低分辨率DEM的地形简化方法的优缺点,提出了考虑整体地形特征信息的一种改进方法,实现了改进后的方法具有较好的简化效果。     

基于Hutchinson的DEM建立及质量评价

建立高质量的数字高程模型(DEM)是正确计算坡度、坡向、提取流域地形特征、进行水文分析的前提。国外应用最广的是基于Hutchinson方法的DEM插值方法和应用该算法的软件ANUDEM,该软件采用有限微分内插技术和地形强化算法,自动去除伪下陷带点和生成输入数据错误文件。研究表明,通过等高线回放、DEM中误差、坡度、河流、光照模拟等方面的对比,ANUDEM生成的DEM表面光滑,比常规用TIN方法构建的TIN-DEM更能准确地表现地形起伏,其提取的坡度、光照图更准确,适宜进行水文分析。     

基于GeoMap的坡度坡向分析

紧密结合等高线的空间数据库和属性数据库,在对等高线进行离散的基础上,研究了不同方法在计算坡度坡向的特点,通过趋势面对任意点邻域进行拟合以及对GeoM ap的二次开发基础上,利用任意点的纵向和横向相邻单元的高程,求出该点的坡度值,并且在确定坡向时要顾及实际的地形情况。结果表明该算法不但精度高,而且计算效率相对与其他算法而言也是高的。