利用ArcGIS水文分析工具提取河網(wǎng)的操作ArcGIS 水文分析工具 提取 河網(wǎng)

DEM包含有多種信息，ArcToolBox提供了利用DEM提取河網(wǎng)的方法，但是操作比較煩瑣（幫助可參看Hydrologic analysis sample applications），今天結(jié)合我自己的使用將心得寫出來與大家分享。提取河網(wǎng)首先要有柵格DEM，可以利用等高線數(shù)據(jù)轉(zhuǎn)換獲得。在此基礎(chǔ)上，要經(jīng)過洼地填平、水流方向計(jì)算、水流積聚計(jì)算和河網(wǎng)矢量轉(zhuǎn)化這幾個不步驟。
1．洼地填平
DEM洼地（水流積聚地）有真是洼地和數(shù)據(jù)精度不夠高所造成的洼地。洼地填平的主要作用是避免DEM的精度不夠高所產(chǎn)生的（假的）水流積聚地。洼地填平使用ArctoolBox->Spatial Analysis Tools->Hydrology->Fill工具。

2．水流方向計(jì)算
水流方向計(jì)算就可以使用上一步所生成的DEM為源數(shù)據(jù)了（如果使用未經(jīng)洼地填平處理的數(shù)據(jù)，可能會造成精度下降）。這里主要使用ArctoolBox->Spatial Analysis Tools->Flow Direction 工具。輸入的DEM采用第一步的Fill1_exam1

3．水流積聚計(jì)算
這里主要使用ArctoolBox->Spatial Analysis Tools->Flow Accumulation工具流向。柵格數(shù)據(jù)就是第二步所獲得的數(shù)據(jù)（FlowDir_fill1）?？梢钥吹剑傻乃鞣e聚柵格已經(jīng)可以看到所產(chǎn)生的河網(wǎng)了?，F(xiàn)在所需要做的就是把這些河網(wǎng)柵格提取出來?？梢园旬a(chǎn)生的河網(wǎng)的支流的象素值作為閥值來提取河網(wǎng)柵格。

4．提取河網(wǎng)柵格
使用spatial analyst中的柵格計(jì)算器，將所有大于河網(wǎng)柵格閥值的象素全部提取出來。至于這個閥值是多少因具體情況而定。通常是要大于積聚計(jì)算后得到柵格的最低河流象素值。這里采用的是500這個值。最后生成只有0、1值的柵格數(shù)據(jù)。其中1表示是河網(wǎng)，0是非河網(wǎng)。

5．生成河網(wǎng)矢量
這里主要使用ArctoolBox->Spatial Analysis Tools->Stream to Feature工具.Input Stream raster 為第四步只有0、1值的河網(wǎng)柵格。流向柵格使用第二步所生成的柵格數(shù)據(jù)。

6．矢量河網(wǎng)處理
由于Stream to Feature工具.將所有柵格象素均轉(zhuǎn)為矢量線段。所以要進(jìn)行處理，方法是利用屬性查詢的方法把所有GRID_CODE為1的全部選擇出來。導(dǎo)出就得到了由dem所生
成的河網(wǎng)矢量。


最后得到的河網(wǎng)如下圖所示。但是由于是柵格轉(zhuǎn)換而來。生成的河網(wǎng)并不是連續(xù)的矢量?？梢愿鶕?jù)需要做簡單的處理。