水文中使用比較多，首先添加數(shù)據(jù)

在分析之前先在擴(kuò)展里啟用這個空間分析功能

開始第一個操作

進(jìn)行如下設(shè)置

？？？怎么又和教程不一樣？

教程的設(shè)置

結(jié)果如圖

然后再創(chuàng)建一個sink去檢查洼地的情況

然后進(jìn)行如下設(shè)置

結(jié)果

取消顯示下面兩層，可以看到

但是這樣對流域分析是有問題的，所以需要進(jìn)行一個填挖操作

設(shè)置輸出選項

輸出結(jié)果

這里就是填挖之后的操作

然后執(zhí)行流域分析

？？？怎么又和教程不一樣？

教程里的設(shè)置

輸出結(jié)果

教程里的結(jié)果

然后進(jìn)行接下來的分析

教程里

輸出結(jié)果

去掉其他的顯示，結(jié)果如圖

修改一下顯示效果

輸出結(jié)果

可以非常清楚的看到河網(wǎng)的雛形

再疊加上流域如圖

這里小的比如500的就提取不到了，可以再改一下表現(xiàn)與分類改成100，可以

效果如圖

對小的河流的末端也提取出來了

然后進(jìn)行下一步分析

可以看到這個分析方式與flowaccumu分類為500的時候結(jié)果幾乎是完全一樣的

重新找到這個工具集

設(shè)置如下

輸出結(jié)果

這里教程崩潰了兩次，用catalog操作了，（當(dāng)然也看不到輸出結(jié)果）

接下來對整個流域進(jìn)行勾畫

教程里的設(shè)置

我的輸出結(jié)果

教程里的輸出結(jié)果

這個時候可以看到我們之前設(shè)置的河網(wǎng)net

其實(shí)可以用arcgis批處理把這些處理全部聯(lián)系起來

先刪除這些，然后調(diào)用arcgis的工具箱

打開py編程器輸入如下代碼：

#本節(jié)用到代碼如下：

import arcpy

# 導(dǎo)入環(huán)境變量

from arcpy import env

# 導(dǎo)入空間分析的模塊（這下面sa就是空間分析的后綴，所以arcpy.sa的就是空間分析的模塊），*就是導(dǎo)入這里空間分析的所有模塊

from arcpy.sa import *

#設(shè)置工作區(qū)，也就是這個路徑，建議使用反斜杠/，而不是順斜杠\，避免報錯

env.workspace = "d:/chap14"

# 檢查一下空間分析擴(kuò)展模塊有沒有問題，教程里運(yùn)行結(jié)果是 u'Available'?，說明是有問題的，如果正常使用的話不會報錯

arcpy.CheckExtension("Spatial")

#流向分析，這里emidalat是這個文件名，教程里運(yùn)行下面這一行后，catalog里多了一個名為flowdir_emid1的文件

outflowdirection =FlowDirection("emidalat")

#確定凹陷區(qū)域，這里的outflowdirection是上一行的變量outflowdirection，教程里運(yùn)行后多了一個名為sink_flowdir1的文件

outsink = Sink(outflowdirection)

#填洼地分析，這里emidalat是這個文件名，教程執(zhí)行后多了一個名為fill_emidala1的文件

outfill = Fill("emidalat")

#確定流向，這里outfill是上面一行的變量，所以不用加雙引號，教程執(zhí)行后多了一個名為flowdir_emid1的文件

outfd = FlowDirection(outfill)

#累積流域分析，這里outfd是上面一行的變量，所以不用加雙引號，教程執(zhí)行后多了一個名為flowacc_flow1的文件

outflowac = FlowAccumulation(outfd)

#條件運(yùn)算，這里outflowac是上面一行的變量，所以不用加雙引號，教程執(zhí)行后多了一個名為con_flowacc_1的文件

outnet = Con(outflowac, 1, 0,"VALUE > 500")

#提取河網(wǎng)，因為這里 outnet , outfd 都是變量，所以都不加雙引號

outstreamlink =StreamLink(outnet,outfd)

#提取子流域，因為這里 outfd , outstreamlink 都是變量，所以都不加雙引號

outwatershed = Watershed(outfd,outstreamlink)

#輸出結(jié)果，這里修改一下，填入需要保存的文件名，這里運(yùn)行后就會保存成文件 outwatershed.tif

outwatershed.save("outwatershed.tif")

可以看到結(jié)果是一樣的