数字高程模型地因素提取.doc

数字高程模型地因素提取————————————————————————————————作者:————————————————————————————————日期: 等高线三维可视化原理H=f(x,y)坡面地形因子提取1坡度打开SpacialAnalysis工具依次选择表面分析、坡度、提取坡度,输出栅格命名为(坡度)2坡向打开SpacialAnalysis工具,依次选择表面分析、坡向,提取坡向,并将输出栅格命名为(坡向)3粗糙度打开栅格计算器,输入公式1/cos[DEM*),即可以得到地表粗糙度,并命名图层为地表粗糙度。4地表起伏度选中DEM数据,打开SpacialAnalysis\邻域分析、焦点统计、后选择统计类型为最大值,邻域类型为矩形记为max;同理再次打开SpacialAnalysis\邻域分析、焦点统计、后选择统计类型为最大值,邻域类型为矩形记为min;打开SpacialAnalysis,地图代数,栅格计算器,输入公式max-min,命名产生的图层为地表起伏度山脊线、山谷线提取操作步骤:,设置默认存储路径,使用空间分析模块下拉箭头中的表面分析工具,选择坡向工具(Aspect),提取DEM的坡向数据层,命名为A。该DEM的坡向数据如下图所示:,使用空间分析模块下拉箭头中的表面分析工具,选择坡度工具(slope),提取A的坡度数据层,命名为SOA1。,记为H:,公式为(H-DEM),求反地形DEM数据如下:反地形DEM数据层calculation如下(可与原始DEM相比较)::,记为SOA2,即利用SOA方法求算反地形的坡向变率。,求没有误差的DEM的坡向变率SOA,公式为SOA=(([SOA1]+[SOA2])-Abs([SOA1]-[SOA2]))/2其中,Abs为求算绝对值,可点击右下侧将其查找出来。没有误差的DEM的坡向变率SOA如下图所示:,使用空间分析工具集中的栅格邻域计算工具(NerghborhoodStatistics);设置统计类型为平均值(mean)邻域的类型为矩形(也可以为圆),邻域的大小为11×11(这个值也可以根据自己的需要进行改变),则可得到一个邻域为11×11的矩形的平均值数据层,记为B。,求算正负地形分布区域,公式为C=[DEM]-[B]。,即可求出山脊线,公式为shanji=[C]>0&SOA>70。说明:和(&),它是比较两个或两个以上栅格数据层,如果对应的栅格值均为非0值,则输出结果为真,(赋值1),否则输出结果为假(赋值0)。选中shanji数据层,将属性值为0的设为空值,透明显示:使用表面分析中的hillshape工具,创建地形晕渲图hillshade,叠加在shanji数据层之下,地图显示效果更佳。