1准备数据
某流域DEM数据
2准备matlab工具包
在官网下载最新TopoToolbox,需要邮箱注册
3数据预处理
在ArcGIS中进行DEM "定义投影",不要地理坐标系,要投影坐标系(UTM),具体带号可参考WGS 1984地理坐标系的墨卡托投影分度带(UTM ZONE)选择方法-CSDN博客
4参考官网ksn的实例代码(如下图),在matlab脚本中进行工具包文件读取和DEM数据替换,再根据实际情况微调
注意1:不同版本的工具包对于函数的命名会有修改,需根据实例代码和工具包函数适当改写,保持一致。
注意2:王一舟老师的B站课程“RiverProAnalysis河流高程剖面分析”,另一种计算ksn指数的途径,可以参考。