不同分辨率dem的图斑坡度精度分析.docxVIP

不同分辨率dem的图斑坡度精度分析 0 dem的精度 土壤坡度（又称坡度）是对土壤坡度的定量描述，也是基本地形特征的单株。目前,在实际应用中最常用的DEM(数字高程模型)有三种类型:规则网格(GRID),不规则网格(TIN)及数字等值线图。而生成DEM的主要方法是利用地形图数字化等高线和高程点信息等,这样则不同比例尺的地形图的等高距不同,生成DEM的精度也就不同。在生成DEM的过程中,由于不同软件采用的内插方法不同,生成的DEM精度也不同。此外,DEM生成后,不同软件采用的计算地面坡度的算法不同,生成图斑的平均坡度的算法也有不同。那么,具体采用何种等高距的等高线生成的DEM数据、采用哪种软件所计算出来的图斑的平均坡度的精度最高,就成了各种土地利用、地形分析、水文模型、滑坡监测与分析、地表物质运动、土壤侵蚀研究者所关心的问题。然而,目前关于对土地利用图斑坡度精度的研究不是太多,不能对上述问题作出全面的回答;再者,因研究区域的限制,目前还没有较为通用的方法。因此,本研究对提高土地利用图斑坡度计算的精度具有一定的重要意义,也对其它应用DEM的行业和部门有一定的借鉴。 1 实验领域的信息和基本原则 1.1 dem数据 实验样区位于云南省东北部,地理坐标在东经103°45′00″～103°52′30″、北纬25°45′00″～25°47′30″,行政区划属于沾益县。坡度类型满足研究要求。研究采用的等高距分别为5、10和20 m,以来自国家测绘部门编制的1∶10 000、1∶50 000地形图和地方测绘部门的1:5 000地形图数据作为数据基础,通过内差得到15 m DEM数据;以实测的图斑坡度数据作为基准,比较分析5、10、15、20 m分辨率下的DEM计算图斑坡度的精度及ARCGIS和MAPGIS软件平台计算图斑坡度的精度,得出最佳的坡度计算方法。本文还分析了坡度大小对DEM计算图斑坡度精度的影响。 1.2 dem方法的计算 地面某点的坡度是地形曲面函数在东西、南北方向上高程变化率的函数,即 Slope=arctanfx2+fy2?????????√?(1)Slope=arctanfx2+fy2?(1) 式中,fx是x方向高程变化率,fy是y方向高程变化率。由(1)式知,求解地表某点的坡度,关键是求解fx和fy。格网DEM是以离散形式表示地形曲面且曲面函数未知,因此在格网DEM上对fx和fy的求解,是在局部3×3移动窗口中通过数值微分方法或局部曲面拟合方法进行(见图1)。 虽然ARCGIS和MAPGIS软件计算坡度采用的都是三阶差分算法,但对fx和fy的计算方法不尽相同。其中ARCGIS采用的是三阶反距离平方权差分,即算法1;MAPGIS采用的是三阶反距离权差分即算法2,具体计算公式如下: 算法1∶fx=v8+2v1+v5?v7?2v3?v68d(2)fy=v7+2v4+v8?v6?2v2?v58d(3)算法1∶fx=v8+2v1+v5-v7-2v3-v68d(2)fy=v7+2v4+v8-v6-2v2-v58d(3) 算法2∶fx=v8+V2v1+v5?v7?V2v3?v68d(4)fy=v7+V2v4+v8?v6?V2v2?v58d(5)算法2∶fx=v8+V2v1+v5-v7-V2v3-v68d(4)fy=v7+V2v4+v8-v6-V2v2-v58d(5) 其中,v1-v8为栅格高程值,d为网格间隔(图1)。当DEM数据具有良好的精度时,数据误差可以忽略不计,fx和fy的结果主要受d的影响,此时DEM地形描述的尺度效应直接决定了坡度计算精度。在起伏变化比较剧烈的区域,DEM对地形描述的精确度受空间分辨率影响很大,直接导致基于不同分辨率DEM提取的坡度也各不相同。 2 基于arcgis的空间分析 本研究采用ARCGIS和MAPGIS软件平台。 其中,MAPGIS为武汉中地信息工程有限公司研制的具有自主版权的地理信息平台。是一个集先进的图形、图象、地质、遥感、测绘、人工智能、计算机科学于一体的大型智能软件系统,也是集数字制图、土地利用现状信息系统管理及空间分析为一体的空间信息系统。 ARCGIS为美国环境系统研究所(Environmental System Research Institute,ESRI)开发的新一代应用广泛的GIS软件之一,是一个全面的、完善的、可伸缩的GIS软件平台。由桌面、服务器、嵌入式和移动GIS软件组成,可以完成数据处理、地理分析、数据编辑、可视化和空间处理等功能,是实用性很强的GIS软件平台。 2.1 dem模型的建立 本研究主要从两个方面入手:对比由5、10、15、20 m等高距的高线生成的DEM计算坡度的精度;对比ARCGIS和MAPGIS软件平台利用DEM提取地面坡度算法的精度。最后将两者综合分