ArcGIS实验操作(八)---地形特征提取

ArcGIS实验操作(八

地形特征信息提取

数据:在data/Ex8/文件下

·dem:分辨率为5米的栅格DEM数据。

·Result文件夹:

·shanji:提取的山脊线栅格数据;

·shangu:提取的山谷线栅格数据;

·hillshade:地形晕渲图。

要求:

利用所给区域DEM数据,提取该区域山脊线、山谷线栅格数据层。

操作步骤:

1.加载DEM数据,设置默认存储路径,使用空间分析模块下拉箭头中的表面分析工具,

选择坡向工具(Aspect,提取DEM的坡向数据层,命名为A。

该DEM的坡向数据如下图所示:

提取A的坡度数据层,命名为SOA1。

3.求取原始DEM数据层的最大高程值,记为H:

由此可见该最大高程值H为1153.79 使用栅格计算器,公式为(H-DEM,求反地形DEM数据如下:

反地形DEM数据层calculation如下(可与原始DEM相比较:

4.基于反地形DEM数据求算坡向值

反地形DEM数据层calculation的坡向数据如下:

5.提取反地形DEM坡向数据的坡度数据,记为SOA2,即利用SOA方法求算反地形的坡

向变率。

6.使用空间分析工具集中的栅格计算器,求没有误差的DEM的坡向变率SOA,公式为

SOA=(([SOA1]+[SOA2]-Abs([SOA1] -[SOA2]/2

其中,Abs为求算绝对值,可点击右下侧将其查找出来。

没有误差的DEM的坡向变率SOA如下图所示:

7.再次点击初始DEM数据,使用空间分析工具集中的栅格邻域计算工具(Nerghborhood

Statistics;设置统计类型为平均值(mean邻域的类型为矩形(也可以为圆,邻域的大小为11×11(这个值也可以根据自己的需要进行改变,则可得到一个邻域为11×11

的矩形的平均值数据层,记为B。

8.使用空间分析工具集中的栅格计算器,求算正负地形分布区域,公式为C = [DEM]-[B]。

9.使用空间分析工具集中的栅格计算器,即可求出山脊线,公式为shanji =

[C]>0&SOA>70。

说明:和(&,它是比较两个或两个以上栅格数据层,如果对应的栅格值均为非0值,则输

出结果为真,(赋值1,否则输出结果为假(赋值0。

选中 shanji 数据层，将属性值为 0 的设为空值，透明显示：

使用表面分析中的 hillshape 工具，创建地形晕渲图 hillshade，叠加在 shanji 数据层之下，地图显示效果更佳。 10. 同上，在栅格计算器中，键入公式 shangu＝[C]<0&SOA>70，可求得山谷线：

将计算结果进行重分类，所有属性不为 1 的栅格属性赋值为 NoData：

将山谷线数据层 re_shangu 与地形晕渲图 hillshade 叠加，地图窗口显示如下：