(3)山谷线的提取:
1)选择SpatialAnalyst下的RasterCalculator,打开栅格计算对话框,在文本框中填写反地形的计算公式:fandem=Abs(dem-2000)单击evaluation按钮,得到与原始DEM地形完全相反的反地形数据。得到反地形数据后,山谷线的提取步骤和山谷线的提取相同,最终利用重分类方法将结果二值化处理。这里不需要对反地形DEM进行洼地填充。
2)利用SpatialAnalyst菜单下的RasterCalculator将重分类后的数据与负地形数据相乘,消除那些存在于正地形区域中错误的山谷线,然后计算结果进行重分类,所有属性不为1的栅格属性值赋为nodata,即获得山谷线,如图所示。