地形鞍部的提取

地形鞍部的提取1.背景相鄰兩山頭之間呈馬鞍形的低谷凹部分稱為鞍部.鞍部點是重要的地形控制點,它和山頂點，山谷點以及山脊線，山谷線等構(gòu)成地形特征點，對地形具有很強的控制作用。因此，對這些地形特征點，線的分析研究在數(shù)字地形分析中具有很重要的意義。同時，由于鞍部點的特殊地貌形態(tài)，使得鞍部點的提取方法較山頂點和谷底點更難，目前還都存在一定的技術(shù)局限性。2.目的：利用水文分析的方法提取地形鞍部點，通過多種GIS空間分析方法的應(yīng)用，提高對知識的綜合運用能力。3.要求：利用水文分析模塊和空間分析模塊相應(yīng)功能提取樣區(qū)地形鞍部點。4.數(shù)據(jù)：25米分辨率的DEM數(shù)據(jù)，面積約為59平方公里。5算法思想：鞍部具有獨特的形態(tài)特征，可被認為是原始地形中的山脊和反地形中的山脊會合的地方，因此可以通過提取正反地形的山脊線并求其交點，獲取鞍部點。6.操作步驟：（1）正地形、等高線和暈渲圖的提?。和郊咕€和山谷線的提取中一樣，由于鞍部點的整體位置是處于山脊上的，需要提取出地形以過濾那些在負地形上的錯誤的點。正地形的提取過程與第一個例子完全相同，提取過程分別是：利用11*11的窗口進行平均值的鄰域分析，結(jié)果為meandem，原始DEM與meandem相減并以0為界進行重分類，大于0的屬性值賦值為1，小于0的賦值為0，結(jié)果命名為zhengdixing。利用SpatialAnalyst菜單下的SurfaceAnalysis菜單中的Contour和Hillshade工具分別提取樣區(qū)等高距為40米的等高線數(shù)據(jù)ctour和樣區(qū)暈渲圖hillshade（2）山脊線的提取山脊線的提取與練習(xí)1（見“小傻帽吧其實我也覺得這個網(wǎng)名起的好難聽啊，可是我覺得好聽的網(wǎng)名都被先注冊的網(wǎng)友們搶了，我就想沒人喜歡用這個網(wǎng)名，所以就用了，哎，其實我也是個白癡性格的人.呵呵，希望看到這里大家都能開心一下，加油啦！”文庫中的“山脊線、山谷線的提取”，就是用到Hydrology水文分析工具的那個word文檔）中山脊線的提取過程完全相同。分別是進行洼地填充水流方向提取匯流累計量計算匯流累積量等于0的提取。提取過程產(chǎn)生的各個數(shù)據(jù)分別為：filldem、flowdir、flowacc以及flowacc0好吧，為了方便大家，我還是在重復(fù)下山脊線提取的過程吧O(_)Oa洼地填充：使用Hydrology工具中的fill工具b無洼地水流方向計算：Hdrology中的flow Dirfillc：匯流累積量計算：Hdrology中的flowAccumulation工具d匯流累積量為0的提取:spaitialAnalyst模塊的下拉箭頭,單擊RasterCalculator,計算:e鄰域處理:f重新分類:g將二值化的neiborfacc00進行重分類為reneibor,將屬性值接近1的那一類的屬性值賦值為1,其余賦值為0h用spatialAnalyst菜單下的RasterCalculator將重分類后的reneibor數(shù)據(jù)與正地形數(shù)據(jù)zhengdixing相乘,消除那些存在于負地形區(qū)域中的錯誤山脊線,然后將計算結(jié)果進行重分類,所有值不為1的柵格賦值為nodata,記得到山脊線.（3）反地形山脊的提取反地形山脊的提取與練習(xí)1中山谷的提取過程完全相同，分別是基于原始DEM計算出反地形DEM數(shù)據(jù)（計算中是利用原始DEM減去常熟3000）；基于反地形DEM數(shù)據(jù)提取水流方向數(shù)據(jù)；基于水流方向數(shù)據(jù)進行匯流累積量數(shù)據(jù)；提取匯流累積量數(shù)據(jù)等于0的柵格。提取過程中產(chǎn)生的數(shù)據(jù)分別為：fandem，flowdirfan，flawccfan以及fanfacc0（4）鞍部點的提取1）利用SpatialAnalyst菜單下的 RasterCalculator的工具將山脊線數(shù)據(jù)flowcc0和山谷線數(shù)據(jù)fanfacc0相乘，結(jié)果命名為anbuqu2）利用同樣的方法將上一步中提取出的數(shù)據(jù)anbuqu和正地形數(shù)據(jù)zhengdixing相乘就得到了鞍部點的柵格數(shù)據(jù)，命名為rasteranbu3）重分類rasteranbu，所有0值賦為NO DATA數(shù)據(jù)，屬性為1的值保持不變，重分類后的數(shù)據(jù)為rasteranbu24）將柵格數(shù)據(jù)rasteranbu2轉(zhuǎn)換成矢量結(jié)構(gòu)