地形平坦地区DEM生成算法的比较研究

数字高程模型(DEM)是GIS的基础数据之一,是许多涉及高程信息的空间分析的重要基础。DEM插值精度的好坏直接影响基于DEM的各种地学分析及应用,因此研究高精度的DEM插值算法具有重要意义。介绍了基于水文强化算法的DEM插值算法——ANUDEM,并与传统的TIN建立DEM的算法进行了比较。认为在地形复杂地区由于等高线密集,高程信息丰富,两种插值算法效果都比较好;而在地形平坦地区ANUDEM算法精度优于TINDEM,并且,ANUDEM得到的DEM很好地避免了平山顶、阶梯状地形,其派生的等高线和水系与原始数据更加吻合。相对于TINDEM算法,ANUDEM算法更适合地形平坦地区的DEM建立。     

ASTER GDEM与SRTM3高程差异影响因素分析

作为最新发布的全球地形数据,ASTER GDEM比目前常用的SRTM3数据有着更高的分辨率和更广的覆盖范围,对于相关地学分析具有重要意义。本文以华中地区为研究区域,对ASTER GDEM与SRTM3数据进行了比较,重点分析了坡度、坡向、地形起伏度、土地利用类型、植被覆盖度、生成ASTER GDEM栅格点高程数据所用的ASTER DEM影像数等因素对2种DEM数据高程差异的影响。结果表明,在研究区域内,ASTER GDEM高程比SRTM3高程平均低5.42 m,两种DEM数据高程差异的RMS值为16.90 m;ASTER GDEM与SRTM3之间的高程差异随着坡度、地形起伏度、植被覆盖度的增大而增大,而ASTER DEM影像数越大,高程差异越小;坡向、土地利用类型对高程差异也有影响。     

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

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

基于机载LiDAR的DEM数据处理及高程基准换算——以上海海岸带区域为例

以上海海岸带区域为例,介绍基于机载LiDAR采集DEM数据的基本处理流程,并重点论述将DEM数据从WGS-84大地高转换到1985国家高程基准正常高的二次曲面拟合模型的构建及精度分析。结果表明,高程基准转换模型的转换中误差为0.05 m,经模型转换后的DEM数据的高程中误差为0.33 m,能够满足海岛海岸带调查1∶10 000比例尺的数据精度要求。     

联合CryoSat-2测高数据和地面高程数据建立东南极拉斯曼丘陵地区DEM

拉斯曼丘陵地区位于东南极伊丽莎白公主地,中国南极中山站位于拉斯曼丘陵的东部,是中国南极科学考察的重要地区。数字高程模型（DEM）是南极冰盖变化研究的基础,卫星测高数据是南极地区构建DEM的主要数据来源。CryoSat-2是新一代用于极地冰盖和海冰监测的测高卫星,联合2013年和2014年南极冬季的CryoSat-2测高数据以及中国、澳大利亚、印度三个国家现场测量的60余个地面高程数据,利用克里金插值方法建立了拉斯曼丘陵地区200 m分辨率的DEM（简称LA-DEM）。利用未参与插值的地面高程数据对新建立的LA-DEM进行了验证,并与Bamber 1km DEM、ICESat DEM、RAMPv2 DEM以及BEDMAP 2等四种国际上常用的南极DEM进行比较,结果表明LA-DEM的高程精度约为19.7 m,优于其他4种南极DEM。     

顾及地形起伏特征的线性预测滤波方法

传统的线性预测滤波算法将目标点云划分为多个栅格,滤波在每一个栅格内进行,而滤波栅格大小,需要用户手动调整。针对此问题,该文提出一种从机载激光扫描数据中生成数字高程模型(DEM)的有效方法。引入统计学变量——面高程变异系数,刻画地形起伏特征,并建立线性预测滤波算法中栅格大小与面高程变异系数之间的函数关系。最后,利用几组点云数据为研究对象验证该方法的有效性,实验结果表明,该方法能自适应地根据地形的起伏特征调整滤波参数,得到比较理想的地面点数据,最终内插得到高精度的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在坡度较大地区和快速变化的冰盖边缘地区精度有较大改进。     

基于ICESat-2 ATLAS数据的DEM高程精度评价

为探究ASTER GDEMV3、SRTM1 DEM和AW3D30 DEM 3种开源DEM数据的高程精度,本文以高精度ICESat-2 ATLAS测高数据为参考数据,利用GIS统计分析、误差相关分析及数理统计对DEM的高程精度进行对比评价。结果表明：①AW3D30的质量最稳定;SRTM1 DEM在平原精度最高;在高原山地精度由高到低依次为AW3D30 DEM、ASTER GDEMV3、SRTM1 DEM。②DEM数据高程精度受地表覆盖影响较大,且与地形因素密切相关,在相同地表覆盖的两个研究区中DEM数据高程精度表现情况不一致,SRTM在平原地表覆盖下精度表现最好,平均误差为3.15 m,AW3D30 DEM在山地地表覆盖下精度表现最好,平均误差为7.61 m。③坡度对DEM数据的高程精度影响较大,在两个研究区3种DEM数据的高程误差均随坡度的增加而增加;坡向对DEM数据的高程精度影响较小,未发现明显的规律。     

南极冰盖DEM机载测高验证与分析——以西南极Thwaites冰川为例

利用冰桥计划(IceBridge)在西南极Thwaites冰川的机载激光测高数据,对ICESat卫星测高数据和目前国际常用的4种南极DEM,包括Bamber 1km DEM、ICESat DEM、RAMPv2DEM和JLB97DEM的精度进行了验证和分析。结果表明,ICESat卫星测高数据和ICESat DEM有着较高的高程可靠性,其与冰桥计划机载测高数据的平均高程差小于5m,标准差小于15m。Bamber 1km DEM高程可靠性相比ICESat卫星测高数据和ICESat DEM低一些。JLB97DEM、RAMPv2DEM与冰桥计划机载测高数据之间的标准差超过30m,尤其在坡度较大的区域,高程可靠性低。     

最佳DEM分辨率的确定及其验证分析

在玛尔挡地区格网DEM的数据上选择实验样区,以不同分辨率情况下DEM数据对地表模拟表达的逼近程度为研究对象,最优逼近时的栅格单元大小的临界值就是所求的最佳分辨率。在分析坡度中误差法和公式法等常见方法的基础上,借鉴坡度中误差的思想,选取区域地形粗糙度K、剖面线长度SL两个定量指标来综合分析确定该地区格网DEM的最佳分辨率。在ArcGIS平台上对方法进行了实验验证,得出分别以2m和8m作为玛尔挡地区1∶10 000和1∶50 000 DEM生产时是最佳分辨率的结论。研究表明这种解决办法不仅可以克服GIS空间分析中DEM分辨率确定的盲目性和随意性,而且能确保基于DEM的各种空间分析的精度,为相关研究提供重要的参考价值。     

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.     

坡度随水平分辨率变化及其空间格局研究

以黄土丘陵沟壑区的县南沟流域为研究区,基于1∶1万地形图,利用ANUDEM软件生成5m到200m分辨率DEM,并利用Arc/Info中计算坡度的方法提取了各种分辨率的坡度。研究表明,随着DEM分辨率的降低,单个样点坡度值表现出不确定性,但同一坡度级所有栅格点的坡度均值呈现一定的规律性,低坡度段表现为先升高后降低,中坡度段呈现微弱变化,陡坡度段呈现对数降低趋势;沟沿线上坡度值呈比较剧烈的下降趋势、分水线和流水线上坡度缓慢下降。     

ICESat-2/ATLAS全球高程控制点提取与分析

ICEsat-2/ATLAS是目前高程精度最高的星载激光数据,其数据覆盖全球,能够作为生产高精度全球地面参考高程的基础数据.基于ICESat-2/ATLAS全球激光数据产品ATL08,获取了全球ICESat-2陆地高程点,研究了基于参考高程数据和属性参数提取全球高程控制点的方法,并利用高精度参考高程数据验证了其精度.利...     

Digital close range photogrammetry for measurement of soil erosion

Many of the processes involved in soil erosion have dimensions on the millimetre scale. Modelling and quantification of such processes require information on soil surface topography with adequate resolution. The purpose of this study was to generate digital elevation models (DEMs) from soil surfaces with high spatial and temporal resolution. Digital photogrammetry was applied for measuring erosion rates on complex-shaped soil surfaces under laboratory rainfall conditions. A total of 60 DEMs were generated, covering a planimetric area of 16 m². The DEMs had a grid resolution of 3 mm. A vertical precision of approximately 1 mm was desired for DEM analysis. A consumer-grade digital camera was used for image acquisition. The camera was calibrated using BLUH software. Homologous points in overlapping images were identified with least squares matching software. Irregularly spaced object coordinates were interpolated to a regular grid in a geographic information system. The resulting DEMs represented the soil surface well. A precision of 1·26 mm in the vertical was attained. The precision of DEM production was limited to camera calibration. Improvements of the setup presented could include the use of better control points and more advanced image matching strategies for identification of homologous points. The DEMs allowed for detailed analysis of soil surface evolution.     

数字高程模型DEM精度研究

用不同地貌类型地区的实验数据对数字高程模型(规则格网DEM)的格网间距与中误差(RMSE)、平均高程、剖面曲率的关系进行研究,采用比较分析的方法,得到格网间距与中误差的数学模型,可以用通式:Y=A BX CX2表示,并且把地形因子(平均高程、剖面曲率)作为自变量,进一步分析通式中3个系数与地形因子的关系,发现两个地形因子与系数函数关系比较相似,并得到Y=A BX CX2 DX3的数学模型。由验证结果表明,所得模型效果良好。     

Evaluation and Comparison of Multi Resolution DEM Derived Through Cartosat-1 Stereo Pair – A Case Study of Damanganga Basin

The study evaluates and compares Digital Elevation Model (DEM) data of various grid spacing derived using high resolution Cartosat 1 stereo data for hydrologic applications. DEM is essential in modeling different environmental processes which depend on surface elevation. The accuracy of derived DEM varies with grid spacing and source. The CartoDEM is the photogrammetric DEM derived from stereo pairs. Damanganga basin lying in the Western Ghats was analysed using 11 Carto stereo pairs. The process of triangulation resulted in RMSE of 0.42. DEM was extracted at 10 m, 20 m, 30 m, 40 m, 50 m and 90 m grid spacing and compared with ASTER GDEM (30 m) and SRTM DEM (90 m). DEM accuracy was checked with Root Mean Square Error (RMSE) statistic for random points generated in different elevation zones. Extracted stream networks were compared based on Correctness Index and Figure of Merit index, calculated for all the Digital Elevation Models at varying cell sizes. In order to further evaluate the DEM’s, a simple flood simulation with no water movement and no consideration of real time precipitation data was carried out and relationship between heights of flood stage and inundation area for each Digital Elevation Model was also established.     

机载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生产具有一定的指导和借鉴意义。     

Comparative evaluation of horizontal accuracy of elevations of selected ground control points from ASTER and SRTM DEM with respect to CARTOSAT-1 DEM: a case study of Shahjahanpur district,Uttar Pradesh,India

Digital elevation model (DEM) data of Shuttle Radar Topography Mission (SRTM) are distributed at a horizontal resolution of 90 m (30 m only for US) for the world, Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) DEM data provide 30 m horizontal resolution, while CARTOSAT-1 (IRS-P5) gives 2.6 m horizontal resolution for global coverage. SRTM and ASTER data are available freely but 2.6 m CARTOSAT-1 data are costly. Hence, through this study, we found out a horizontal accuracy for selected ground control points (GCPs) from SRTM and ASTER with respect to CARTOSAT-1 DEM to implement this result (observed from horizontal accuracy) for those areas where the 2.6-m horizontal resolution data are not available. In addition to this, the present study helps in providing a benchmark against which the future DEM products (with horizontal resolution less than CARTOSAT-1) with respect to CARTOSAT-1 DEM can be evaluated. The original SRTM image contained voids that were represented digitally as ?140; such voids were initially filled using the measured values of elevation for obtaining accurate DEM. Horizontal accuracy analysis between SRTM- and ASTER-derived DEMs with respect to CARTOSAT-1 (IRS-P5) DEM allowed a qualitative assessment of the horizontal component of the error, and the appropriable statistical measures were used to estimate their horizontal accuracies. The horizontal accuracy for ASTER and SRTM DEM with respect to CARTOSAT-1 were evaluated using the root mean square error (RMSE) and relative root mean square error (R-RMSE). The results from this study revealed that the average RMSE of 20 selected GCPs was 2.17 for SRTM and 2.817 for ASTER, which are also validated using R-RMSE test which proves that SRTM data have good horizontal accuracy than ASTER with respect to CARTOSAT-1 because the average R-RMSE of 20 GCPs was 3.7 × 10^?4 and 5.3 × 10^?4 for SRTM and ASTER, respectively.     

机载LiDAR点云生成多尺度DEM方法

为了利用机载激光雷达点云生成高保真、多尺度的数字高程模型(DEM),提出了一种基于综合生成策略的方法:首先,利用点云数据中的地面点生成高分辨率、高保真的DEM作为基础DEM;然后,通过迭代的方式对上一层较高分辨的DEM进行综合获取较低分辨率、高保真的DEM。实验表明,本文方法不仅具有可行性,而且生成的多尺度DEM具有高保真的特性。     

Improving a digital elevation model by reducing source data errors and optimising interpolation algorithm parameters: An example in the Loess Plateau,China

A hydrologically correct digital elevation model (DEM) forms a basis for realistic environmental modelling, especially in complex terrain. We have performed a study in the Coarse Sandy Hilly Catchments (CSHC) of the Loess Plateau, China, which demonstrates pragmatic, yet effective methods for improving the quality of the DEM by: (1) identifying and correcting source topographic data errors and (2) optimising ANUDEM algorithm parameters. Improvement in the DEM based on fixing over 1100 errors in the input topographic data, and optimising key ANUDEM parameters was assessed using higher accuracy independent validation of 32 contributing areas and 1474 spot heights, and by semi-quantitative analysis of DEM derivatives produced from ANUDEM and Triangular Irregular Network (TIN) algorithms. Improvement in the ANUDEM DEM over the original TIN DEM was shown where the percentage of the total absolute difference in contributing areas reduced from 10.43 to 3.51%, and the bias between the spot heights and DEM elevations reduced from 45 to 32 m. Large improvement in DEM quality was gained by using ANUDEM instead of TIN, with smaller improvement gained by fixing source data errors, and optimising ANUDEM parameters.