Answers:
将您的流向量转换为值为1且具有与DEM相同的范围和像元大小的栅格。在栅格计算器中,使用地图代数表达式遵循以下方式:
Con("rivers"==1, "DEM" - 1, "DEM")
如果要在多于1个高程单位的流中燃烧,请更改"DEM" - 1
为"DEM" - a bigger value
。
要在ArcGIS 10栅格计算器中通过@elrobis答案实现Whitebox GAT衰减系数算法(Z = E-(G /(G + D))^ k * H),请使用以下方法:
"DEM" - Power (G / (G + EucDistance("rivers_raster") ), k ) * H
将G替换为您的单元格大小,将k替换为大于0(衰减系数)的值,将H替换为最大数量以降低DEM。例如:
"DEM" - Power (30 / (30 + EucDistance("rivers_raster") ), 2 ) * 5
Whitebox GAT(开源水文学和遥感软件包)在其水文实用程序中使用了这种名称的方法。Whitebox的独特之处在于它通过UI公开了分析所使用的源代码和算法(请注意查看代码按钮)。即使您打算将过程隔离到ArcGIS,尝试另一种方式也可能会有一些好处。
对于Arcgis,卢克的答案似乎是完美的,但是正如在另一个答案中提到的WhiteBox GAT一样,我不禁在其他“风味”上添加一些其他信息,以防可能有所帮助。
SAGA-GIS(免费)具有一些有趣的预处理功能-其中之一可以回答您的问题:http : //www.saga-gis.org/saga_tool_doc/2.1.3/ta_preprocessor_6.html
您将需要栅格DEM和栅格格式的河网(“ 形状转换为网格”可以完成此操作-其他任何数据都应为nodata)
模块将流网络刻录到DEM
将流网络刻录到数字高程模型(DEM)中。流单元必须使用有效的数据值进行编码,所有其他单元应设置为无数据值。前两种方法减少。第三种方法可确保稳定的下游梯度。仅当下游单元格相同或更高时才应用海拔降低。您应该提供一个带有流动方向的网格以确定下游单元格。“水槽排水路径检测”模块提供了这样的流向。
作者:O.Conrad(c)2011
规格:网格
菜单:地形分析|预处理