现平面曲率的误差。所以,提取平面曲率的主要工作就是消除北坡的误差。这里,引入负地形的概念:负地形是和真实的地表高度相反,它通过一个计算(H(DEM),把原来高点变低,把原来)低点变高,H的取值可以给定一个远大于此区域高程值的值如10000。负地形的实现可以通过【A
alysis】下的【Calculator】命令,公式为DEMH(1),即可生成负地形。下面讨论如何通过负地形消除北坡的平面曲率误差。首先,利用正负地形分别提取平面曲率,因为负地形的平面曲率的误差出现在南坡,这样,正负地形的平面曲率的误差就对称
f分布于北坡和南坡,记这两个层面分别为A、B。当这两个层面相加时,没有误差的地方值翻倍,有误差的地方除了其本身的值以外,还使误差得到了累积。当这两个层面相减时,没有误差的地方值为零,有误差的地方只剩下了累积的误差,误差有正有负,对相减得到
图115
平面曲率在北坡存在误差
的新层面取绝对值。最后,对相加、相减得到的两个新层面再相减,即得到了没有误差的翻倍的平面曲率的值,对此值除以2,即可得到了没有误差的平面曲率的层面。平面曲率的提取方法如下:1、提取DEM层的坡向主题,再对此坡向主题提取坡度,得到的主题记为A;2、在【A
alysis】菜单下使用【Calculator】命令,公式为DEMH(1),提取DEM层的负地形;3、提取负地形的坡向的坡度,记为B;4、在【A
alysis】菜单下使用【Calculator】命令,公式为ABABabs2,即可求出没有误差的DEM的平面曲率。如图116
图116消除误差后的平面曲率图
f另外,根据求出的平面曲率,通过【A
alysis】菜单的【Mapquery】命令,选取平面曲率80的要素即可提取出沟脊线和沟谷线提取结果如图117。
图117提取的沟脊线和沟谷线
五、地形起伏度地形起伏度是指在一个特定的区域内,最高点海拔高度与最低点海拔高度的差值。它是描述一个区域地形的一个宏观性的指标。从地形起伏度的定义可以看出,求地形起伏度的值,首先要求出一定范围内海拔高度的最大值和最小值,然后,对其求差值即可。求一定范围内的最大值和最小值,可以通过用【A
alysis】菜单的【Neighborhoodstatistics】分别设置statistic值为最大值和最小值,邻域的设置可以为圆,也可以为矩形,邻域的大小可根据自己的要求来确定。地形起伏度的具体提取方法如下:1、激活DEM数据,在【A
alysis】菜单中使用【Neighborhoodstatistics】命令,
图118地形起伏度
f设置statistic值为最大值,邻域的类型为矩形(也可以为圆),邻域的大小为100个Mr