齊通通實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)

1、1 洼地填充。在 ArcMap 中加載 DEM 數(shù)據(jù)（dem-grid），在【ArcToolBox】中，選擇【Spatial Analyst Tools】，選擇【Hydrology】，雙擊【Fill】 2.在彈出的窗口中，【Input surface raster】中輸入需要被填充洼地的 DEM；【Output surface raster】中輸入 Fill 結(jié)果柵格數(shù)據(jù)的名稱和路徑；【Z limit (optional)】是高差閾值，如 果設(shè)置其為 10，則當(dāng)填充高差大于 10 時(shí)，該洼地不會(huì)被填充，默認(rèn)值是填充所有洼地（圖 4-1）。單擊【OK】得到無洼地 DEM（ 計(jì)算水流方向，選擇【H

2、ydrology】，雙擊【Flow Direction】，在彈出窗口中（圖 4-3）， 【Input surface raster】中輸入填充洼地后的 DEM；【Output flow direction raster】中輸入流向 結(jié)果柵格數(shù)據(jù)的名稱和路徑 單擊ok計(jì)算匯流累積量。選擇【Hydrology】，雙擊【Flow Accumulation】；在彈出窗口中（圖 4-5），【Input flow direction raster】中輸入計(jì)算得到的水流方向數(shù)據(jù)；【Output accumulation raster】中輸入?yún)R流累積量計(jì)算結(jié)果的文件名與路徑；【Input weight ras

3、ter (optional)】中可選 擇是否輸入權(quán)重?cái)?shù)據(jù)；【Output data type (optional)】中輸入輸出數(shù)據(jù)的類型。單擊【OK】 匯水網(wǎng)絡(luò)分級。選擇【Hydrology】，雙擊【Stream Order】，在彈出窗口中（圖 4-9）， 【Input stream raster】輸入?yún)R流累積量重分類結(jié)果；【Input flow direction raster】輸入水流方 向數(shù)據(jù)；【Output raster】輸入分級數(shù)據(jù)的文件名和路徑；【Method of stream ordering (optional)】 輸入分級方法。 建立匯水網(wǎng)絡(luò)連接。選擇【Hydrology】

4、，雙擊【Stream Link】。在彈出窗口中（圖 4-12）， 【Input stream raster】輸入?yún)R流累積量重分類結(jié)果；【Input flow direction raster】輸入水流方 向數(shù)據(jù)；【Output raster】輸入 Stream Link 數(shù)據(jù)的文件名和路徑； 矢量匯流網(wǎng)絡(luò)。選擇【Hydrology】，雙擊【Stream to Feature】（圖 4-13）?！綢nput stream raster】輸入?yún)R水網(wǎng)絡(luò)連接結(jié)果；【Input flow direction raster】輸入水流方向數(shù)據(jù)；【Output raster】 輸入矢量匯流網(wǎng)絡(luò)數(shù)據(jù)的文件名和路

5、徑；【Sinplify polylines (optional)】指是否簡化矢量線.基礎(chǔ)實(shí)驗(yàn)二 在 ArcMap 中加載 DEM 數(shù)據(jù)（dem-grid），對 DEM 洼地填充（Fill），并求取水流方 向數(shù)據(jù)（Direction）。 2. 選擇【Hydrology】，雙擊【Basin】。彈出窗口中，【Input flow direction raster】中輸入 流向數(shù)據(jù)；【Output raster】中輸入 Basin 計(jì)算結(jié)果的文件名與路徑將柵格數(shù)據(jù)轉(zhuǎn)化為矢量數(shù)據(jù)。選擇【ArcToolBox】，選擇【Conversion Tools】，選擇 【From Raster】，雙擊【Raster

6、 to Polygon】（圖 4-17），得到基于 Basin 的流域多邊形對 DEM 洼地填充，計(jì)算水流方向，計(jì)算匯流累積量，運(yùn)用重分類工具，以 1000 閾 值提取匯水網(wǎng)絡(luò)柵格。 2. 選擇【Hydrology】，雙擊【Stream Link】，在彈出窗口中，【Input stream raster】輸入 匯流累積量重分類結(jié)果；【Input flow direction raster】輸入水流方向數(shù)據(jù)；【Output raster】輸 入 Stream Link 數(shù)據(jù)的文件名和路徑雙擊【W(wǎng)atershed】，彈出窗口中，【Input flow direction raster】輸入水流方向

7、數(shù)據(jù)；【Input raster or feature pour point data】輸入 Stream Link 柵格數(shù)據(jù)；【Output raster】輸入 watershed 流域分割結(jié)果的文件名與路徑（ 將柵格數(shù)據(jù)轉(zhuǎn)化為矢量數(shù)據(jù)（Raster to Feature）（圖 4-23）。其中黑色部分是無數(shù)據(jù) 區(qū)，即在本方法下，運(yùn)用 1000 閾值無法分割的區(qū)域。 基礎(chǔ)三 進(jìn)行鄰域分析，求取大值。在 ArcMap 中加載 DEM 數(shù)據(jù)（dem-grid），選擇【 Spatial Analyst】工具，單擊【Neighborhood Statistics】。彈出窗口中，【Neighborho

8、od】選擇“Maximum”， 分析窗口設(shè)置為“17*17”。具體的分析范圍根據(jù)地貌類型有所變化，地形破碎的區(qū)域，分析 窗口應(yīng)適當(dāng)縮小，地形平坦區(qū)域，分析窗口應(yīng)適當(dāng)擴(kuò)大。選擇【Spatial Analyst】工具，單擊【Raster Calculator】，輸入公式“NbrMax of dem-grid - dem-grid”計(jì)算選擇【Spatial Analyst】工具，單擊【Reclassify】，重分類成兩類， 0 即 0，其余為 NoData. 由于有平頂山的存在，使得部分區(qū)域呈現(xiàn)山頂點(diǎn)聚集狀態(tài)。將柵格數(shù)據(jù)轉(zhuǎn)成矢量面。 選擇【ArcToolBox】，選擇【Conversion Tool

9、s】，選擇【From Raster】，雙擊【Raster to Polygon】。 彈出窗口中，【Input raster】輸入局部高點(diǎn)柵格；【Output polygon features】輸入生成多邊形 的文件名與路徑 選擇【ArcToolBox】，選擇【Data Management Tools】，選擇【Features】，雙擊【Feature To Point】。彈出窗口中，【Input Features】輸入生成的局部高點(diǎn)多邊形；【Output Feature Class】 輸入山頂點(diǎn)的文件名與路徑；. 在 ArcMap 中加載 DEM 數(shù)據(jù)（dem-grid），提取流域邊界，作為局

10、部高程點(diǎn)候選。提 取方法見前面流域分割。 2. 將流域邊界轉(zhuǎn)化為柵格數(shù)據(jù)。選擇【ArcToolBox】，選擇【Conversion】，選擇【To Raster】，雙擊【Polyline To Raster】。彈出窗口中，【Input Features】輸入流域邊界數(shù)據(jù)，【Output Raster Dataset】輸入流域邊界柵格數(shù)據(jù)的文件名和路徑；【Cellsize (optional)】選擇與原始 DEM 一致（圖 4-30）。單擊【Environments】，在彈出窗口中，單擊【General Settings】，在 【Extent】中輸入原始 DEM，得到柵格流域邊界 進(jìn)行柵格計(jì)算，使

11、得流域邊界柵格值都為 1。選擇【Spatial Analyst】工具，單擊【Raster Calculator】，輸入公式：water_bound = 0，得到邊界數(shù)據(jù) 進(jìn)行柵格計(jì)算，求得流域邊界高程。選擇【Spatial Analyst】工具，單擊【Raster Calculator】，輸入公式：Calculation * dem-grid對流域邊界高程進(jìn)行鄰域分析，求取小值。選擇【Spatial Analyst】工具，單擊 【Neighborhood Statistics】。彈出窗口中，【Statistic Type】選擇“Minimum”；分析窗口選擇 “17*17”，獲得分析結(jié)果數(shù)據(jù)

12、柵格計(jì)算，求取分析窗口下，流域邊界的局部低柵格。選擇【Spatial Analyst】工 具，單擊【Raster Calculator】，輸入公式：Calculation2 - NbrMin of Calculation2，獲得計(jì) 算結(jié)果 通過柵格計(jì)算，選擇 0 值柵格，作為鞍部點(diǎn)備選點(diǎn)。選擇【Spatial Analyst】工具， 單擊【Raster Calculator】，輸入公式：Calculation3 = 0。并將結(jié)果命名為“min_bound”。 8. 通過對原始 DEM 進(jìn)行鄰域分析，求取正負(fù)地形。選擇【Spatial Analyst】工具，單 擊【Neighborhood St

13、atistics】。彈出窗口中，【Statistic Type】選擇“Mean”；分析窗口選擇 “17*17”，獲得分析結(jié)果 對計(jì)算結(jié)果重分類。其中，大于 0 的為 1（正地形），小于 0 的為 0（負(fù)地形）。選 擇【Spatial Analyst】工具，單擊【Reclassify】，獲得柵格分類結(jié)果柵格計(jì)算，求出既是正地形上，且又是流域邊界局部低點(diǎn)的柵格。選擇【Spatial Analyst】工具，單擊【Raster Calculator】，輸入公式：min_bound * Reclass of Calculation4， 獲得計(jì)算結(jié)果 打開計(jì)算結(jié)果屬性表，選中結(jié)果為 1 的柵格。選擇【Ar

14、cToolBox】，選擇【Conversion】， 選擇【From Raster】，雙擊【Raster to Polygon】（圖 4-40）。將鞍部轉(zhuǎn)為矢量多邊形。 選擇【ArcToolBox】，選擇【Data Managerment Tools】，選擇【Features】，雙擊【Feature To Point】（圖 4-41）。將多邊形轉(zhuǎn)成點(diǎn)，消除邊界上的點(diǎn)，得到終鞍部點(diǎn)（基礎(chǔ)實(shí)驗(yàn)4 在 ArcMap 中加載 DEM 數(shù)據(jù)（dem-grid），提取 1000 閾值下的流域邊界鄰域分析，計(jì)算平均值。選擇【Spatial Analyst】工具條，單擊【Neighborhood Statistics】。 彈出窗口中，【Statistic type】選擇“Mean”，分析窗口選擇“17*17”，得到分析結(jié)果柵格計(jì)算。選擇【Spatial Analyst】工具條，單擊【Raster Calculator】。輸入公式：dem-grid - NbrMean of dem-grid，得到計(jì)算結(jié)果柵格重分類，區(qū)分正負(fù)地形。選擇【Spatial Analyst】工具條，單擊【Reclassify】，得 到重分類結(jié)果Reclass_aste1 正地形矢量多邊形生成。選擇【ArcToolBox】，選擇【