ArcGIS水文分析綜合作業(yè)實(shí)習(xí)報(bào)告

./實(shí)習(xí)水文分析班級(jí)級(jí)班專(zhuān)業(yè)地理信息系統(tǒng)學(xué)號(hào)姓名得分一實(shí)習(xí)內(nèi)容和意義二數(shù)據(jù)準(zhǔn)備數(shù)據(jù)來(lái)源：國(guó)際科學(xué)數(shù)據(jù)平臺(tái)<中國(guó)科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心科學(xué)數(shù)據(jù)中心>數(shù)據(jù)類(lèi)型：SRTM90M數(shù)據(jù)的選取：登陸國(guó)際科學(xué)數(shù)據(jù)平臺(tái),在原始dem數(shù)據(jù)中選擇SRTM90M數(shù)據(jù),進(jìn)入下載頁(yè)面,圖1選擇合適的范圍,下載的ASC=2\*ROMANII數(shù)據(jù)名稱(chēng)依次為srtm_59_04.zip、srtm_59_05.zip、srtm_59_06.zip。三涉及的基本概念四具體操作步驟〔以汾河為實(shí)例1、流域提取1.1數(shù)據(jù)轉(zhuǎn)換：將解壓后的數(shù)據(jù)<.asc>添加到global_mapper中,圖2將數(shù)據(jù)的投影坐標(biāo)系統(tǒng)轉(zhuǎn)化為所需要的高斯克呂格坐標(biāo),輸出為帶有高程的柵格數(shù)據(jù)。圖31.2數(shù)據(jù)處理打開(kāi)ArcMap,加載工具箱,在toraster下拉條中選擇ASC=2\*ROMANIItoraster,圖4打開(kāi)如圖所示對(duì)話框,加入數(shù)據(jù)。圖5轉(zhuǎn)化為柵格數(shù)據(jù)后,打開(kāi)SpatialAnalystTools中的hydrology第一步：洼地填充:打開(kāi)fill對(duì)話框進(jìn)行洼地填充圖6第二步：計(jì)算水流方向數(shù)據(jù)圖7第三步：計(jì)算水流累計(jì)量數(shù)據(jù)圖8根據(jù)計(jì)算的水流方向數(shù)據(jù)和累計(jì)量數(shù)據(jù)和實(shí)驗(yàn)的流域閾值找出流域范圍。圖9選取所需要的流域區(qū)域：圖10提取流域的原始dem數(shù)據(jù)圖112、高程情況2.1Tin圖的形成使用arctoolbox中3d的convertion進(jìn)行數(shù)據(jù)格式的轉(zhuǎn)化,如圖所示：圖12得到流域tin圖圖132.2高程分級(jí)選擇spatialanalyst下的reclassify對(duì)dem數(shù)據(jù)的高程重新分類(lèi)得到新的高程分類(lèi)數(shù)據(jù)圖14將高程分級(jí)數(shù)據(jù)屬性表導(dǎo)出,在excel中進(jìn)行面積比例計(jì)算,如下圖圖152.3坡度分級(jí)在spatialanalyst下的surface中打開(kāi)slope對(duì)話框,加載數(shù)據(jù)如圖所示,得到坡度數(shù)據(jù)圖16打開(kāi)reclassify對(duì)話框,加入坡度數(shù)據(jù),進(jìn)行重分類(lèi)：圖17得到坡度分級(jí)數(shù)據(jù)：圖18將屬性表導(dǎo)出,得到不同坡度的面積及所占的比例圖192.4地表起伏度利用spatialanalyst下的鄰域分析,設(shè)置參數(shù)如下圖所示,得到起伏度數(shù)據(jù)圖202.5地面粗糙度利用上一步得到的坡度數(shù)據(jù),利用柵格計(jì)算器得到地面粗糙度數(shù)據(jù)層圖21圖222.6山脊線、山谷線1點(diǎn)擊dem數(shù)據(jù),使用表面分析中的坡向工具,提取dem的坡向數(shù)據(jù)層,命名為A2點(diǎn)擊數(shù)據(jù)層A,使用表面分析中的坡度工具,提取數(shù)據(jù)層A的坡度數(shù)據(jù),命名為soa13>求取原始dem數(shù)據(jù)層的最大高程值,記為H,使用空間分析工具集中的柵格計(jì)算器,公式為H-dem,得到與原來(lái)地形相反的數(shù)據(jù)層,即反地形dem數(shù)據(jù)圖234基于反地形dem數(shù)據(jù)求算坡向值5利用soa方法求算反地形的坡向變率,記為soa26利用空間分析工具集中的柵格計(jì)算器,公式為soa=〔〔soa1+soa2-abs〔soa1-soa2/2即可求出沒(méi)有誤差的dem坡向變率soa圖247再次點(diǎn)擊初始dem數(shù)據(jù),使用空間分析工具集中的柵格鄰域計(jì)算工具；設(shè)置統(tǒng)計(jì)類(lèi)型為平均值,鄰域類(lèi)型為矩形,鄰域大小為11*11,可得到平均值數(shù)據(jù)層B8使用空間分析工具集中的柵格計(jì)算器,公式為c=dem-B,即可求出正負(fù)地形分布區(qū)域9使用空間分析工具集中的柵格計(jì)算器公式為shanji=c>0&soa>50,即可求出山脊線圖2510同理,在柵格計(jì)算器中,鍵入公式為山谷=c<0&soa>50,即可求出山谷線圖262.7山頂點(diǎn)單擊空間分析下拉菜單,選擇鄰域分析,設(shè)置參數(shù)如圖所示,單擊ok按鈕,提取11*11分析窗口最大值圖27單擊空間分析下拉箭頭,選擇柵格計(jì)算器輸入公式sd=H-dem==0,提取山頂點(diǎn)區(qū)域圖28對(duì)sd數(shù)據(jù)進(jìn)行重分類(lèi)得到山頂點(diǎn)數(shù)據(jù)圖29選擇convert子菜單并單擊rastertofeature,輸出山頂點(diǎn)數(shù)據(jù)圖302.8河網(wǎng)利用原始的dem數(shù)據(jù)得到填充洼地的dem、水流方向數(shù)據(jù)、水流累積量數(shù)據(jù)。利用空間分析工具下的柵格計(jì)算器計(jì)算柵格河網(wǎng),得到柵格河網(wǎng)圖31利用工具箱的轉(zhuǎn)換工具將柵格河網(wǎng)轉(zhuǎn)為矢量河網(wǎng)圖322.9溝壑密度打開(kāi)矢量河網(wǎng)數(shù)據(jù)屬性表,添加length字段,計(jì)算長(zhǎng)度,在字段上右鍵選擇statistic,進(jìn)行統(tǒng)