我們經(jīng)常會使用研究區(qū)的DEM數(shù)據(jù)進(jìn)行水文分析，在進(jìn)行水文分析之前，需要進(jìn)行洼地的填充，然后進(jìn)行流向、流量的分析，最后根據(jù)最符合研究區(qū)現(xiàn)實(shí)情況的閾值，提取更加精確的河網(wǎng)，然后可以將河網(wǎng)矢量化，最后根據(jù)實(shí)際情況，將不符合現(xiàn)實(shí)情況的河網(wǎng)刪去。
1、加載DEM數(shù)據(jù)，這里的DEM數(shù)據(jù)是NASA_DEM_n22e113，具體的現(xiàn)在可以到earthdata免費(fèi)下載，首先將研究區(qū)的DEM按掩膜提取，以中山市為例
2、然后進(jìn)行一個填洼的處理，將研究區(qū)范圍內(nèi)的洼地填充：
填洼之后我們可以看到，研究區(qū)的DEM最高點(diǎn)沒有變化，但是最低點(diǎn)變高了，說明研究區(qū)內(nèi)存在一定數(shù)量的洼地：
3、然后進(jìn)行流向分析，流向分析可選擇三種方法，流向工具支持三種流向建模算法。分別為 D8、多流向 (MFD) 和 D-Infinity (DINF)：
D8 流向法可對每個像元到其最陡下坡鄰域的流向進(jìn)行建模。
以 D8 流向類型運(yùn)行的流向工具的輸出是值范圍介于 1 到 255 之間的整型柵格。從中心出發(fā)的各個方向值為：
例如，如果最陡下降方向位于當(dāng)前處理像元的左側(cè)，則將該處理像元的流向編碼將為 16，如果最陡下降方向位于當(dāng)前處理像元的右側(cè)，則將該處理像元的流向編碼將為 1，如果最陡下降方向位于當(dāng)前處理像元的左上，則將該處理像元的流向編碼將為 32，以此類推。
這里我們選擇D8法：
然后進(jìn)行一個流量分析：
然后進(jìn)行柵格計(jì)算器，將流量大于300的提取出來（這里的300是一個閾值，閾值的設(shè)置不是固定的，需要根據(jù)研究區(qū)的現(xiàn)實(shí)狀況，進(jìn)行合理的選擇，也可以進(jìn)行多次試驗(yàn)進(jìn)行合理選擇，本例選擇300）：
得到河網(wǎng)后放大觀察，是樹狀的河流網(wǎng)絡(luò)：
最后，可以將柵格的河網(wǎng)進(jìn)行矢量化，得到矢量的河網(wǎng)圖層，當(dāng)然最后的結(jié)果實(shí)際情況可以疊加到衛(wèi)星圖上進(jìn)行比對，然后再針對流量的閾值進(jìn)行調(diào)整，最后使用更加精確的閾值得到河網(wǎng)，當(dāng)然最后的河網(wǎng)還是會和現(xiàn)實(shí)情況有出入，這個時候進(jìn)行矢量化之后手動刪除一些不符合現(xiàn)實(shí)情況的河網(wǎng)即可。
本文鏈接：http://www.hx-bj.cn/blog/151.html
本文標(biāo)簽：ArcGIS
