DEM计算坡度slope问题(为什么坡度常在89°左右)
问题:通常我们下载的数据坐标都是地理坐标系,没有投影信息。直接根据不含投影信息的数据计算坡度,通常计算的结果都是89.0左右(单:degree),显然这是有问题的(不可能坡度全是近90°)。
分析:主要是因为,DEM数据的高程单位为米,而没有投影信息的数据的水平距离单位为度(°),肯定是有问题的,应该保持相同的单位(米)来计算坡度。
坡度计算图图示
- 表示坡度最为常用的方法,即两点的高程差与其水平距离的百分比,其计算公式如下:
坡度 = (高程差/水平距离)x100%
2. 度数法:用度数来表示坡度,利用反三角函数计算而得,其公式如下:(常用)
tanα(坡度)= 高程差/水平距离
所以α(坡度)= arctan(高程差/水平距离)
问题解决方法:
添加投影信息:因此,在根据DEM计算坡度时(无论是利用ARCGIS 还是ENVI 都需要预先为数据添加投影信息):
(ARCGIS: 数据管理工具--投影和变换--栅格--投影栅格--选择投影坐标系(设定))
生成slope: 添加完投影信息之后,在进行坡度计算:
(ARCGIS: 3D analyst--栅格表面--坡度)