DEM包含有多种信息,ArcToolBox提供了利用DEM提取河网的方法,但是操作比较烦琐(帮助可参看Hydrologic analysis sample applications),今天结合我自己的使用将心得写出来与大家分享。提取河网首先要有栅格DEM,可以利用等高线数据转换获得。在此基础上,要经过洼地填平、水流方向计算、水流积聚计算和河网矢量转化这几个大步骤。
1.洼地填平
DEM洼地(水流积聚地)有真是洼地和数据精度不够高所造成的洼地。洼地填平的主要作用是避免DEM的精度不够高所产生的(假的)水流积聚地。洼地填平使用ArctoolBox->Spatial Analysis Tools->Hydrology -> Fill工具。
2.水流方向计算
水流方向计算就可以使用上一步所生成的DEM为源数据了(如果使用未经洼地填平处理的数据,可能会造成精度下降)。这里主要使用ArctoolBox->Spatial Analysis Tools->Flow Direction 工具。输入的DEM采用第一步的Fill1_exam1
3.水流积聚计算
这里主要使用ArctoolBox->Spatial Analysis Tools->Flow Accumulation工具流向。栅格数据就是第二步所获得的数据(FlowDir_fill1)。可以看到,生成的水流积聚栅格已经可以看到所产生的河网了。现在所需要做的就是把这些河网栅格提取出来。可以把产生的河网的支流的象素值作为阀值来提取河网栅格。
4.提取河网栅格
使用spatial analyst中的栅格计算器,将所有大于河网栅格阀值的象素全部提取出来。至于这个阀值是多少因具体情况而定。通常是要大于积聚计算后得到栅格的最低河流象素值。这里采用的是500这个值。最后生成只有0、1值的栅格数据。其中1表示是河网