山脊,指山高处像兽脊凸起的部分。是由两个坡向相反坡度不一的斜坡相遇组合而成条形脊状延伸的凸形地貌形态。山脊最高点的连线就是两个斜坡的交线,叫做山脊线。等高线,向低处凸的是山脊。山脊是连成一排的山峰,山峰之间连成一条长线,好像动物的脊骨有一条突出的线条,故名山脊线。这里,我们介绍一下如何使用ArcGIS生成山脊线。
数据来源本教程所使用的数据是水经注微图中下载的高程DEM数据,除了高程数据,常见的GIS数据均可以下载。
水经注微图
(私信回复“微图4.0”获取微图软件安装包)
提取正地形为了消除那些存在于负地形区域中的错误的山脊线,最后需要正地形和得到的结果进行乘法运算,所以需要先准备正地形栅格数据。
在ArcToolbox内点击“Spatial Analyst工具\邻域分析\焦点统计”,调用焦点统计工具。
调用焦点统计工具
在显示的焦点统计对话框内,输入下载的DEM数据,设置邻域高度和宽度为10,如下图所示。
焦点统计设置
点击确定之后得到统计栅格,如下图所示。
统计栅格
在ArcToolbox内点击“Spatial Analyst工具\地图代数\栅格计算器”,调用栅格计算器工具。
调用栅格计算器工具
在显示的栅格计算器对话框内,输入下面公式:"DEM.tif" - "FocalSt_tif9",如下图所示。
栅格计算器设置
点击确定之后得到计算结果,如下图所示。
栅格计算结果
在ArcToolbox内点击“Spatial Analyst工具\重分类\重分类”,调用重分类工具。
调用重分类工具
在显示的重分类对话框内,以值0为界,小于0赋值为0,大于0赋值为1,如下图所示。
重分类设置
点击确定之后可以看到重分类后的效果,如下图所示。
重分类效果
提取山脊线
在ArcToolbox内点击“Spatial Analyst工具\水文分析\填洼”,调用填洼工具。
调用填洼工具
在显示的填洼对话框内,输入表面栅格数据为下载的DEM数据,如下图所示。
填洼设置
点击确定后可以得到填洼后的栅格,如下图所示。
填洼后的栅格
在ArcToolbox内点击“Spatial Analyst工具\水文分析\流向”,调用流向工具。
调用流向工具
在显示的流向对话框内,输入表面栅格为上一步填洼后的栅格,如下图所示。
流向设置
点击确定之后可以得到流向的效果,如下图所示。
流向效果
在ArcToolbox内点击“Spatial Analyst工具\水文分析\流量”,调用流量工具。
调用流量工具
在显示的流量对话框内,输入流向栅格数据为上一步得到的流向结果,如下图所示。
流量设置
点击确定之后得到流量栅格,如下图所示。
流量栅格
在ArcToolbox内点击“Spatial Analyst工具\地图代数\栅格计算器”,调用栅格计算器工具。
调用栅格计算器工具
在显示的栅格计算器对话框内,输入下面的公式:"FlowAcc_Flow5" == 0,获取流量为0的区域。
栅格计算器设置
点击确定后得到计算结果,如下图所示。
栅格计算结果
通过邻域分析平滑山脊线,在ArcToolbox内点击“Spatial Analyst工具\邻域分析\焦点统计”,调用焦点统计工具。
调用焦点统计工具
在显示的焦点统计对话框内,输入栅格为上一步栅格计算后的栅格,如下图所示。
焦点统计设置
点击确定之后得到统计结果,如下图所示。
统计结果
为了判断重分类的阈值,需要先生成等高线和山体阴影帮助判断。
在ArcToolbox内点击“Spatial Analyst工具\表面分析\等值线”,调用等值线工具。
调用等值线工具
在显示的等值线对话框内,输入栅格选择下载的DEM数据,等值线间距设置为30,如下图所示。
等值线设置
点击确定之后生成等高线,效果如下图所示。
生成的等高线
在ArcToolbox内点击“Spatial Analyst工具\表面分析\山体阴影”,调用山体阴影工具。
调用山体阴影工具
在显示的山体阴影对话框内,输入栅格为下载的DEM数据,如下图所示。
山体阴影设置
点击确定之后可以看到生成的山体阴影,如下图所示。
生成的山体阴影
在焦点统计的栅格图上点击右键,选择属性,如下图所示。
选择属性
在显示的图层属性对话框内,选择符号系统选项卡,选择已分类,类别设置为2,以等值线图和晕渲图作为辅助判断,不断地调整分类的临界值,属性值越接近1的栅格越有可能是山脊线的位置,最终确定分类阈值为0.52585。
查找阈值
在ArcToolbox内点击“Spatial Analyst工具\重分类\重分类”,调用重分类工具。
调用重分类工具
在显示的重分类对话框内,将阈值(0.52585)以下赋值为0,阈值(0.52585)以上赋值为1,如下图所示。
重分类设置
点击确定之后可以得到重分类结果,如下图所示。
重分类结果
为了排除负地形的影响,需要将重分类结果和正地形相乘。
在ArcToolbox内点击“Spatial Analyst工具\地图代数\栅格计算器”,调用栅格计算器工具。
调用栅格计算器工具
在显示的栅格计算器对话框内,输入下面的公式:"Reclass_Foca4" * "Reclass_rast9",其中"Reclass_Foca4" 为上一步重分类结果,"Reclass_rast9"为最开始得到的正地形。
栅格计算器设置
点击确定之后可以得到栅格计算结果,如下图所示。
栅格计算结果
在ArcToolbox内点击“Spatial Analyst工具\重分类\重分类”,调用重分类工具。
调用重分类工具
在显示的重分类对话框内,输入栅格为上一步的计算结果,将0赋值为NoData,1赋值为1,如下图所示。
重分类设置
点击确定之后可以得到山脊线,如下图所示。
生成的山脊线
如果需要将山脊线提取成矢量,可以在ArcToolbox内点击“转换工具\由栅格转出\栅格转折线”将山脊线栅格转换为矢量。
山脊线转矢量
结语
以上就是如何使用ArcGIS生成山脊线的详细说明,主要包括了数据来源、提取正地形和提取山脊线等功能。如果在使用过程中有任何疑问,请联系我们客服或拔打24小时热线电话400-028-0050咨询!
,