gis中提取山脊线和山谷线

cai8567 2022-05-10 21:06:44

山脊线、山谷线的提取

（1）提取SOA坡向变率

①点击Spatial Analyst Toola，使用Surface中的Aspect工具生成坡向，输入DEM数据，得到aspect数据。接着利用Slope工具，输入上一步提取的aspect数据，得到SOA1数据。

图 1得到SOA1数据操作

②点击Spatial Analyst Tools中的Map Algebra，使用Raster Calculator工具，输入公式为FDEM=2000-dem（得到反地形要求用大于该DEM高程的值减去DEM），得到与原来地形相反的DEM数据层FDEM。接着重复于①中的操作一致，处理得到其坡向及利用坡向继续求坡度得到SOA2。

图 2计算得到反地形

③利用栅格计算器，得到纠正后的坡向变率。点击Spatial Analyst Tools中的Map Algebra，使用Raster Calculator工具，输入公式(("SOA1" + "SOA2") - Abs("SOA1" - "SOA2") )/ 2，点击OK，完成坡向变率的提取。

图 3计算得到SOA

（2）点击Spatial Analyst Tools下的Reclass模块，选择Reclassify工具，输入坡向数据。在Classification界面中设置2个分类段，设置节点为70，返回Reclassify界面，设置大于70为1，其他为Nodata。

（3）提取正负地形部分

①求高程平均值，点击Spatial Analyst Tools中的Neighborhood分析工具，选择Focal Statistics,在弹出的对话框中输入DEM数据，并改变窗口大小为11*11，统计类型选择Mean，点击OK得到平均值。

图 5求平均高程操作

②利用栅格计算器，区分正地形和负地形部分，打开Raster Calculator工具，输入dem-mean(高程减去平均值，大于0则为正地形，小于0为负地形)。接着利用重分类的操作，将正负地形赋值为1和-1，操作如（图7）所示，在Classification界面中设置节点为0，返回Reclassify界面，进行赋值。