写在前面
之前参加一个比赛,提供了中国的高程数据,可以基于该数据进一步计算坡度和坡向进行相关分析。
对于坡度和坡向,这里分享一个找到的库,可以方便快捷的计算。这个库为: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
(加微信交流群,请备注:姓名-行业-单位)
声明:本公众号分享的资源均来自于互联网的公开资源,仅供个人学习和研究使用,不代表本人对其版权、真实性、合法性、准确性等方面做出任何保证。任何使用本号分享的资源的个人或组织,需自行承担相应的法律责任和后果;本号不承担任何法律责任和连带责任。使用本号分享的资源即代表您已完全了解并同意本免责声明的所有内容。
上海位温信息科技有限公司
专业的商业气象服务商
(咨询数据详情,可添加微信,请备注:合规数据)
可提供气象产品 :
全国三千多个基准站
3213 个主要城市站点、每小时发 布 1 次,
每天发布 24 次,每次提供当前整点时刻的
城市站点实况天气(实况产品会在整点后
15 分钟、20 分钟、30 分钟进行数据更新),
要素包括:天气现象、气温、体感温度、风速、
风力、风向、相对湿度、气压、1小时降水量、
10分钟降水量、能见度、日出时间、日落时间、
紫外辐射等要素。
探空气象站
海洋气象站
雷达拼图
卫星云图(风云4、葵花8)
格点实况、再分析资料
数值预报:CMA EC GFS
强对流天气预报
强天气落区预报
强对流天气(短时强降水/冰雹/雷暴大风)概率预报图
台风、热带气旋
可提供国内与国外各类气象咨询气象服务
本公司提供服务稳定可靠,您值得拥有。
上海位温信息科技有限公司
专业的商业气象服务商
(咨询数据详情,可添加微信,请备注:合规数据)
可提供气象产品 :
全国三千多个基准站
3213 个主要城市站点、每小时发 布 1 次,
每天发布 24 次,每次提供当前整点时刻的
城市站点实况天气(实况产品会在整点后
15 分钟、20 分钟、30 分钟进行数据更新),
要素包括:天气现象、气温、体感温度、风速、
风力、风向、相对湿度、气压、1小时降水量、
10分钟降水量、能见度、日出时间、日落时间、
紫外辐射等要素。
探空气象站
海洋气象站
雷达拼图
卫星云图(风云4、葵花8)
格点实况、再分析资料
数值预报:CMA EC GFS
强对流天气预报
强天气落区预报
强对流天气(短时强降水/冰雹/雷暴大风)概率预报图
台风、热带气旋
可提供国内与国外各类气象咨询气象服务
本公司提供服务稳定可靠,您值得拥有。