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表示是河网,0是非河网。
5.生成河网矢量
这里主要使用ArctoolBox->Spatial Analysis Tools->Stream to Feature工具.Input Stream raster 为第四步只有0、1值的河网栅格。流向栅格使用第二步所生成的栅格数据。
6.矢量河网处理
由于Stream to Feature工具.将所有栅格象素均转为矢量线段。所以要进行处理,方法是利用属性查询的方法把所有GRID_CODE为1的全部选择出来。导出就得到了由dem所生
成的河网矢量。
最后得到的河网如下图所示。但是由于是栅格转换而来。生成的河网并不是连续的矢量。可以根据需要做简单的处理。
ArcGIS Spatial Analyst包含一些从水文和地形信息中处理和获取新信息的专门工具。
当进行水流建模时,您需要了解水流的来源和去向. ArcGIS Spatial Analyst提供了用于进行地形表面汇流计算的工具, 它为描绘汇流网络和汇水盆地、流长计算以及确定水系级别等提供了必要的基础.该类数据通常用于将地形信息综合到水文模型中。
地下水工具可以用于进行简单的地下水二维水流和要素建模。达西流量(Darcy Flow)工具用于从水文地理数据中生成地下水流域。达西流量工具可以从一个点源出发,穿过流域,跟踪汇流路径。多孔喷嘴(Porous Puff)工具可用于计算水中化学物品或成分沿流域扩散的状况。
可利用的水文工具包括:
汇水量计算工具(Flow Accumulation)
汇水方向计算工具(Flow Direction)
流长计算工具(Flow Length)
填充(Fill)
下沉接收(Sink)
水系(Stream Link)
水系级别(Stream Order)
水系至特征(Stream To Feature )
分水岭(Watershed)
达西流量(Darcy Flow)
颗粒物跟踪(Particle Track)
多孔喷嘴(Porous Puff)