已知數(shù)據(jù):DEM文件ASTGTM_N00E118E.img(文末有數(shù)據(jù)下載鏈接)
使用軟件:ArcMap
要求:對(duì)數(shù)據(jù)進(jìn)行操作,提取數(shù)據(jù)文件中的山脈范圍
下面開(kāi)始操作:
1、 打開(kāi)ArcMap將DEM文件ASTGTM_N00E118E.img添加到數(shù)據(jù)框。
2、 接下來(lái)我們打開(kāi)spatial analyst工具——鄰域分析——焦點(diǎn)統(tǒng)計(jì),
輸入柵格選擇數(shù)據(jù)框中的dem數(shù)據(jù),輸出柵格名字改為MAXIMUM,輸出位置為默認(rèn)數(shù)據(jù)庫(kù)。鄰域設(shè)置高度及寬度一般為30,本次數(shù)據(jù)較小直接選擇默認(rèn)值,單位選擇像元,統(tǒng)計(jì)類型選擇MAXIMUM,選擇完成后點(diǎn)擊確定
3、 重復(fù)上一步驟打開(kāi)焦點(diǎn)統(tǒng)計(jì),輸出柵格名稱為MINIMUM,統(tǒng)計(jì)類型選擇MINIMUM,其他選項(xiàng)同上一步相同。
4、 打開(kāi)Arctoolbox——spatial analyst工具——地圖代數(shù)——柵格計(jì)算器,
將生成的兩個(gè)柵格拖入框內(nèi)進(jìn)行運(yùn)算,運(yùn)算規(guī)則 如下:用MAXIMUM減MINIMUM得出地形起伏度。
輸出柵格依舊為默認(rèn)數(shù)據(jù)庫(kù),名稱改為qfd,柵格計(jì)算時(shí)用于計(jì)算的柵格數(shù)據(jù)像元大小必須一致,否則會(huì)導(dǎo)致計(jì)算失敗。
Qfd生成后關(guān)閉其他柵格,我們可以看到以下圖像
為了更直觀的觀察起伏度可以將色帶修改,單擊色帶
修改后顯示如下
單擊識(shí)別工具查看起伏差值
點(diǎn)擊紅色區(qū)域發(fā)現(xiàn)值為0,即平地地區(qū)圖片
點(diǎn)擊山脈邊界處查看值,多選擇幾處山脈邊界查看值。以均值6位山脈邊界對(duì)柵格進(jìn)行重分類
打開(kāi)ArcToolbox——spatial analyst 工具——重分類——重分類。
輸入柵格為qfd,重分類字段選擇默認(rèn),將舊值分為兩個(gè)條目,刪除其他值
將舊值改為0-5,6-392,輸出柵格為bj
生成的山脈邊界如下
對(duì)比原始DEM
疊加后查看山脈邊界與dem的套合程度
本文鏈接:http://www.dzac.com.cn/blog/188.html
本文標(biāo)簽:ArcGIS