Python | 由高程计算坡度和坡向

写在前面

之前参加一个比赛，提供了中国的高程数据，可以基于该数据进一步计算坡度和坡向进行相关分析。

对于坡度和坡向，这里分享一个找到的库，可以方便快捷的计算。这个库为：RichDEM，官网地址如下

• https://richdem.readthedocs.io/en/latest/

安装

支持windows环境的conda安装，命令很简单

conda install conda-forge::richdem

示例

安装成功后，读取高程数据。

对于坡度和坡向的计算仅仅需要简单的两行代码：

import richdem as rd
dem_path = r'I:/chinadem_geo.nc'
#### 坡度
slope = rd.TerrainAttribute(rd.LoadGDAL(dem_path), attrib='slope_riserun')
# ## 坡向
aspect = rd.TerrainAttribute(rd.LoadGDAL(dem_path), attrib='aspect')

绘图

下面展示了一下相关绘图结果

总结

相关代码和测试数据放到了GitHub中，感兴趣的可以试试

• https://github.com/Blissful-Jasper/jianpu_record