arcgis水文分析实列

练习9
水文分析：根据DEM提取河流网络，计算流水累积量、流向、根据指定的流域面积大小自动划分流域
水文分分析工具 (1)
1. 数据基础：无洼地的DＥＭ3ﻩ
2．关键步骤：流向分析 ........................................................................................................................... ４３. 计算流水累积量.. (5)
4.
提取河流网络6ﻩ
5．流域分析9ﻩ
６. 其它应用--降水分析（选做) (11)
水文分分析工具
有两种途径使用水文分析功能:
(1)通过Arｃtoｏlbox:水文分析工具位于[Ｓpaｔial AnalystＴools]＞>［Hｙｄｒｏlogy］之下
如果Ｈｙdｒoｌogy工具
集没有出现,可以选中某个
工具箱后新建一个工具集
[Hydrology］，然后右键点新
建的工具集，在出现的菜单
中执行[添加]>＞[工具］
会出现如右图所示的对话
框,将需要的水文分析工具
添加到上面新建的工具集
中。

(2)另一种方法是添加[Ｈyｄrologｙ]工具栏到ArｃMａp中。

在AｒcMap中执行菜单命令:[工具]>>［定制] 命令
点击[从文件添加]按钮
找到ｅsrｉhydｒology_v２.ｄll 文件
注意：这个文件通常是在ArｃGIＳ的安装路径下,默认的情况是[C:\PrograｍＦileｓ\ＡrcＧIS＼DevｅloperKit＼samples\SｐａtialAnaｌysｔ＼ＨｙｄｒoｌoｇicＭodelｉng\Visual_Basic]
［Ｈyｄroｌogy Modeｌiｎg]工具条就被加载到ArcMap,在其前面的检查框上打上勾，如下图所示。

[Hｙdroloｇy Moｄeling］工具条就可以显示在ArcMａｐ中
注意：以下的练习基于HydrｏlogｙＭodeｌing工具
1.
数据基础:无洼地的ＤＥM
在ArｃMap中加载DEＭ数据，执行工具条[Ｈydrology Ｍodeling]中的菜单命令[Hydrolｏgy］>>[ Ｆｉll Sinks]，在出现的对话框中将[Iｎput Suｒｆacｅ］参数指定为“DＥM”
确定后得到无洼地的DEM数据:[Ｆiｌled Siｎk1 ]
2．关键步骤：流向分析
在上一步的基础上进行，执行工具条[Hydrｏlogy Mｏdｅling］中的菜单命令［Hydroｌog ｙ]>>[ Flow Directｉon ],在出现的对话框中将［Ｉｎput Surfaｃe]参数指定为“Filled Sｉnk１”
确定后得到流向栅格[ FloｗDirection１],了解流向栅格单元的数值表示的含义是什么
3．计算流水累积量
在上一步的基础上进行，执行工具条［Hydroloｇy Modｅling］中的菜单命令[ Hyｄrｏloｇｙ]>＞[ F ｌoｗAccumuｌａtｉon ］，在出现的对话框中将[Diｒｅｃtｉｏn Raster］参数指定为“FｌowＤir ｅcｔion１”
确定后得到流水累积量栅格[Flｏw Accｕmulａtion1]
４. 提取河流网络
(1)提取河流网络栅格:
在上一步的基础上进行，打开Arctoｏlｂox，运行工具[Spatｉal Analｙst Tooｌs]＞>[Mａp Aｌgebrａ]＞＞［单输出地图代数]
在[地图代数表达式]中输入公式：con (ＦloｗAｃｃuｍulａtｉon1>800,１)
[输出栅格]指定为:SｔreamＮet
说明：通过此操作将流水累积量栅格[Flow Aｃcumulation1]中栅格单元值(流水累积量）大于８00的栅格赋值为1，从而得到河流网络栅格[StreamNet]
１sｅtnulｌ
2 ｃoｎ
３pick uｐ
这三个命令很有用。

setｎull,将栅格图像中满足条件的栅格值设置为null，即nodata,这样很便于在栅格图像中删除不需要的部分。

其对应Griｄ模块命令为setnull,格式：outgrid ＝setnull ＜ingｒｉd > 5，iｎgrid>,表示，如果ingriｄ》5,则变为nodａta,否则返回为原值。

cｏｎ表示条件,如果满足条件，返回为指定值1,否者返回为指定值２,Grｉｄ命令：
outgrｉd＝con＜ingrid >５，１0 ，１0０＞,表示，如果ｉngrid值＞５,则对应栅格值改为１0,否者改为１0０。

若100不指定,则返回为noｄata。

这里可以指定多个条件。

pick up，我还没看到呢。

ps。

发现ｔｏｏｌbox有个特点，右击命令的话可以批处理哦。

得到的的河流网络栅格：ＳtreaｍＮet
关闭除[Ｓtreamneｔ]之外的其它图层
(2)提取河流网络矢量数据
在上一步的基础上进行,执行工具条[Ｈydrology Moｄelinｇ]中的菜单命令[Hydrologｙ]＞>[ Stｒeam Netｗｏrk ＡｓFeatuｒe ],在出现的对话框中将[DｉrecｔioｎＲaster］参数指
定为“FloｗDｉrecｔioｎ1”，[Accumuｌａtion Rasteｒ]参数指定为“ＦlowＡccumｕlation1”,[MｉnｉmuｍNumbｅｒｏｆCellｓfor a Ｓtream］参数指定为10０0
确定后得到河流网络矢量数据
(3)平滑处理河流网络
打开[编辑器]工具栏,执行工具栏中的命令[编辑器]>>[开始编辑],确保目标图层为河流网络图层[Shaｐｅ１]，
通过打开[Shape1属性表,并选择属性表的所有行选择图层[Ｓｈape1]中的所有要素,也可以通过要素
选择按钮选择图层中所有要素
执行[编辑器]工具栏中的命令[编辑器］＞>[更多的编辑工具]>>[高级编辑]打开工具条：［高级编辑],点击其上的[平滑］按钮(下图中前头所指）：
在[平滑]处理对话框中输入参数[允许最大偏移］:3
得到平滑后的河流网络矢量图层，执行命令：[编辑器]＞>［停止编辑],保存所做修改。

比较平滑处理后的数据与没有进行处理过的数据
5．流域分析
在上一步的基础上进行，执行工具条[Ｈｙdroｌoｇy Mｏdeliｎg]中的菜单命令[ Hydｒologｙ］>>［Ｗaｔeｒｓheｄ］，在出现的对话框中将［DiｒecｔioｎRａstｅｒ］参数指定为“FｌｏｗDiｒeｃtion1”,[Accｕｍulation Raster]参数指定为“Flow Accuｍｕlatｉon1”，［M ｉniｍｕm Nuｍｂer of Cells foｒａＳｔrｅａm]参数指定为2000
确定后得到流域栅格［Wａｔershed１］
打开[空间分析] 工具栏,执行命令:［空间分析］>＞[转换]＞>[栅格到要素］将流域栅格转换成为矢量图层，按下图所示指定参数:
得到矢量数据：［WaterSｈed.shp]
设置图层[Watershｅd]属性[符号］后得到类似上图的效果
重复以上操作步骤并将[MinｉｍumＮuｍber of Cｅlls for aＳtreaｍ］参数指定为1000０,看看结果是否有所不同。

６．其它应用--降水分析(选做)
执行命令[Hydroloｇｙ］>>[ＩnteｒａcｔivｅProｐertieｓ］
选中“Rａin Drop”按钮，（如下图红色前头所指),在DEM上任意位置点击(表示这里有降雨），则显示一条路径表示此处的降水的最终流向。