基于DEM的水文分析

基于DEM的水文分析介绍：基于基于DEM的水文分析的主要内容是利用水纹分析工具提取水流方向、汇流累积量、水流量积量、水流长度、河流网络、河网分级以及流域分割。（一）无洼地DEM生成DEM被认为是比较光滑的地形表面的模拟，但是由于内插的原因以及一些真实地形（如采石场或喀斯特地貌）的存在，使得DEM表面存在着一些凹陷的区域。这些区域在进行地表水流模拟时，由于低高程栅格的存在，从而使得在进行水流流向计算时得到不合理的或错误的水流方向，因此，在进行水流方向的计算之前，应该首先对原始DEM数据进行洼地填充，得到无洼地的DEM。数据：DEM数据dem（1）原始DEM数据提取水流方向执行[Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]>>[FlowDirection]在[FlowDirection]对话框中，“Forcealledgecellstoflowoutward（Optional）”的复选框前打钩，则所有在DEM数据边缘的栅格的水流方向全部流出DEM数据区域（默认为不选择）。呂船吿KDD电IB寓卩：LSZjE旳Untitled-ArcMap-Ard.nfcSpatialAnalyst呂船吿KDD电IB寓卩：LSZjE旳Untitled-ArcMap-Ard.nfcSpatialAnalyst▼二]:必獅尊□X皤Cie-n^^is^'calAnalystD卑laygauK：匕JSdediDTi354643.2223894554.187MeteFsfowdirB日Layer希El回Li^eir：jflQiriirValueHigh:227fileEditViewBoolkmairkEInsertSeiectioriToolsWindowHelp“dropraster”是该栅格在其水流方向上与其临近的栅格之间的高程差与距离的比值，以百

分比的形式记录，它反映了在整个区域中最大坡降的分布情况（可选步骤）。（2）洼地计算执行[Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]>>[Sink]。沪Sink(■=■|回Inputfclowdirectichlr：=letel'1flu'i'i'iiir-罔1Outpnitr：=leterJ|C：\UeersVwqiang\DesktopX培训勘据l空间分折勤据让蚯11\tutor\5inkT-OKCancelEnvironments...ShowHelp>>1

3)洼地深度计算1、双击Hydrology工具集中的Watershed工具。|_1M耳Ge-oistatisticdlAnal/d:Untitled-/V匚Map-Ar-clnfowater&h疝|_1M耳Ge-oistatisticdlAnal/d:Untitled-/V匚Map-Ar-clnfowater&h疝n4c31■■-&H1S143a&1251.G27Meters电■■:峙.■:■■:■Low;1□□flowdirBGLayersEl回rpiE□SinkU'alueHiph:358.file£drt里iuv/bookmarksInse-rtSelection1001=5VZndowJJe-lpDi剜巧」0呃亡］SdeEti&TfSpatialAnalyst|Lftycr：|flowair□3R昌|恳瞄(E计算每个洼地所形成的贡献区域的最低高程a.打开SpatialAnalystTools工具箱中Zonal工具集，双击ZonalStatistic工具回□KmtfESt&.+isticstyp也(epti-!irL：il]Inpu七v：回□KmtfESt&.+isticstyp也(epti-!irL：il]Inpu七v：=J.uerast已工I'ieniOutpnitr：±EtfiY『；lU?已匸2认qi泣吃订i沙kt口f'培彳I徽据、空间分析数据辽验11'■■.=皿』miInputrasterorfeaturezonedatCancelShowHelp>>Environments...沪ZonalStatisticsZonefield臣IgnoreIToDataincalculatioue(optional)计算每个洼地贡献区域出口的最低高程即洼地出水口高程a打开SpatialAnalystTools工具箱中Zonal工具集，双击Zonalfill工具

计算洼地深度a在SpatialAnalyst模块的下拉箭头。选择RasterCalculator命令，弹出RasterCalculator对话框：：：：：RasterCalculatorLayers:CalculationdemflowdirSinkwBtershsinkzonalmaxzonalmin0II〔I)I删zonalminAboutBuildingEKpressionsEvaluartaCancel>>[zonalmaz]-[zonalmm]AboutBuildingEKpressionsEvaluartaCancel>>7)洼地填充洼地填充是无洼地DEM生成的最后一个步骤。经过洼地提取之后，可以确定原始DEM上是否存在洼地，若有洼地，则需要进行填充。洼地深度的计算为填充阀值的设置提供了依据。a.执行[Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]>>[Fill]。

b.在Zlimite文本框中输入阀值。在洼地填充果过程中，洼地深度大于阀值的地方不填充，作为真实地形保留。系统默认情况是不设阀值。即所有的洼地区域将被填平。Zj血JXeeir.ii±i口□hEl口l£l口ElLjEl□Zj血JXeeir.ii±i口□hEl口l£l口ElLjEl□::lna.-iis2U-ntitled-ArcMapArclniraAnalystLow:1120Dfeplay_鈕匹.訂站訥⑷senkdep⑥®芸囂■鋼■令h售起器d詰6.506芥疗二二/甘hi-—、2on3lrri3xzonslminv/atersh寸nkSinkflowdi『demSpatialAnalystTLaya：|flpifli7E曰Layers□HfilldemValueHigh!1625.-DQFile_|TitViewB&DkmarfcsInsert^electionToolsWindowH^lp注：当一个洼地被填平之后，这个区域与附近区域再进行洼地计算，可能还会形成新的洼地。因此，洼地填充是一个不断反复的过程，直到所有的洼地都被填平、新的洼地不再产生为止。二）汇流累积量

汇流累积量数值矩阵表示区域地形没点的流水累积值。在地表径流模拟过程中，汇流累积量是基于水量方向数据计算得到的。汇流累积量的基本思想：以规则格网表示的数字地面高程模型每点处有一个单位的水量，按照自然水流从高处往低处的自然规律，根据区域地形的水流方向数据计算每点处所流过的水流量数值，便得到了该区域的汇流累积量。又水流方向数据到汇流累积量计算的过程③汇流第积耀在地表径流模拟过程中，汇流累积暈是吨于水流方向数据计算而来的。时每一个栅恪来说，K-71.隆秦积艮的大小代表着雄上游有务少牛柵格的木诡方向攝终汇流经M陵栅格.汇诡案和的数値越大，询区威越易形成地表從流。执行[Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]执行[Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]>>[FlowAccumulation]。a在inputweightraster文本框中输入权重数据，权重数据一般是降水、土壤以及植被等影响径流分布不平衡因素综合而成，对每一个栅格权重能更详细模拟该区域的地表特征。如果五权重数据，系统默认所有栅格的权重为1（三）水流长度执行[Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]>>[FlowLength]。计算方向分别为Downstream（顺流计算）和Upstream（朔流计算）。在“Inputweightraster”文本框中输入权重数据，Downstream记录着其沿着水流方向到下游流域出水口中最长距离所流经的栅格数；Upstream则记录着其沿着水流方向到上游栅格的最长的距离的栅格数。

(四)河网的生成(1)提取河流网络栅格①设定阈值。阈值的设定在河网的提取过程是很重要的，并且直接影响到河网的提取结果。阈值的设定应遵循科学、合理的原则。首先应该考虑到研究的对象，研究对象中的沟谷的最小级别，不同级别的沟谷所对应的不同的阈值；其次考虑到研究区域的状况，不同的研究区域相同级别的沟谷需要的阈值也是不同的。②提取栅格河流网络方式一：打开Arctoolbox,运行工具［SpatialAnalystTools］»［MapAlgebra］〉〉［单输出地图代数］。在［地图代数表达式］中输入公式：con(flowacc>800,1)，［输出栅格］指定为：StreamNet。计算的思想：利用所设定的栅格阈值进行条件查询并将查询结果赋予新的栅格数据中。通过此操作将流水累积量栅格flowacc］中栅格单元值(流水累积量)大于800的栅格赋值为1，而小于或等于设定阈值的栅格属性值赋为无数据，从而得到河流网络栅格［StreamNet］。方式二：在ArcMap中，点击空间分析模块〉〉计算器,在计算器中输入表达式：con(flowacc〉800,1)。(2)栅格河流网络矢量化

执行[Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]>>[streamtofeature]（3）streamlink的生成Streamlink记录着河网中的一些节点之间的连接信息（河网的结构信息）。Streamlink的每条弧段连接着两个作为出水点或汇合点的结点，或者连接着作为出水点的结点和河网起始点。通过Streamlink的计算，即得到每一个河网弧段的起始点和终止点。同样，也可以得到该汇水区域（流域）的出水口。执行[Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]>>[Streamlink]。streamLink[|回OKCancelEn^ronments・・・ShowHelp»经过计算，它将栅格河网在汇合点栅格处分割成河网片段，并将片断进行记录，在属性表中除了记录该片段的ID号之外，还记录着每个片段所包含的栅格数。

[Streamlink][Streamlink]结果可利用[streamtofeature]转换为矢量数据。須AttributesofstreamlinltIT回〔■ESWRowidVSLUE*COUNT□►131122C324144515€1€fi7ss489a9L：：：13101151112121213101214£1412Record:M|4||1k|k|JShow:|AllBel亡cted|Records(0outof1???SelectEd)d（5）河网分级在ArcGIS的水文分析中，提供两种常用的河网分级方法：Strahler分级和Shreve分级。对于Strahler分级来说，它将所有河网弧段中没有支流的河网弧段分为第1级，两个1级河网弧段汇流成的河网弧段为第2级，如此下去分别为第3级，第4级，一直到河网出水口。在这种分级中，当且仅当同级别的两条河网弧段汇流成一条河网弧段时，该弧段级别才会增加，对于那些低级弧段汇入高级弧段的情况，高级弧段的级别不会改变，这也是比较常用的一种河网分级方法。对于Shreve分级而言，其第1级河网的定义与Strahler分级是相同的,所不同的是以后更高级别的河网弧段,其级别的定义是由其汇入河网弧段的级别之和执行[Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]>>[StreamOrder]。［StreamOrder］［StreamOrder］结果可利用［streamtofeature］转换为矢量数据ASLmhlcr^fe^i分了刿坡ASLmhlcr^fe^i分了刿坡JShrtic仆级鈿果「仔了135斛（五）流域分析（1）流域盆地的确定流域盆地是由分水岭分割而成的汇水区域。它通过对水流方向数据的分析确定出所有相互连接并处于同一流域盆地的栅格。流域盆地的确定首先是要确定分析窗口边缘的出水口的位置也就是说，在进行流域盆地的划分中，所有的流域盆地的出水口均处于分析窗口的边缘。当确定了出水口的位置之后，也就是找出所有流入出水口的上游栅格的位置。执行[Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]>>[Basin]。回:卫八n目亠口5国日*E3鱼F的QB25.mKd-ArcMap-回:卫八n目亠口5国日*E3鱼F的QB25.mKd-ArcMap-ArclnfoFile-lEdil-ViewBookrriarkiTnsert^electionToolsWindowHelpLi¥"r'SpatialAr^ly^t寄5nurEESdeclnriDia^Ery■Low:1S回streflmostr■Lew:1-0shr^'tValueHigh:253守|规总OJ.pHf?Geostaftisrlic^ilAnalyst-|3333^1.5115553351.965MrarsS目LAyvn-0ValueHi@h:176（2）集水流域的生成经过上一步得到的流域盆地是一个比较大的流域盆地，在很多的水文分析中，还需要基于更小的流域单元进行分析，那么就需要进行