ArcGIS试验-Ex19-黄土地貌鞍部地提取

1、第十一章水文分析练习2:黄土地貌鞍部的提取相邻两山头之间呈马鞍形的低凹局部称为鞍部，鞍部是两个山脊和两个山谷会合的地方。鞍部点是重要的地形控制点，它和山顶点、山谷点以与山脊线、山谷线等构成的地形特征点线，具有对地形具有很强的控制作用。因此， 对这些地形特征点、 线的分析研究在数字 地形分析中具有很重要的意义。 同时，由于鞍部点的特殊地貌形态， 使得鞍部点的提取方法 较山顶点和山谷的提取更难， 目前没有什么有效的方法来提取鞍部点， 利用水文分析的方法 可以来提取一些鞍部点，但是它还是具有一定局限性。、目的结合水文分析的思想来提取研究区域的地形鞍部点；开拓思维，拓展水文分析工具的应 用，学会利用水

2、文地质分析的思想来解决一些地形、地貌等方面的问题。三、要求结合水文地质分析的方法和空间分析的方法提取研究区域的地形鞍部点。四、数据一幅252辨率的黄土地貌 DE赚据，数据的区域大概有 59 km2。数据存放于 /ChP11/ Ex2中，请将其拷贝到 E: /ChP11/Ex2。结果保存在/ChP11/Ex2/Result 中。五、算法思想算法思想：由于鞍部是两个山脊和两个山谷会合的地方，那么对于鞍部点的提取，就可以采用分别提取山脊线和山谷线，然后再计算出山脊线与山谷线的交点，所求出的交点便是鞍部点的位置。六、操作步骤1、正地形、等高线和晕渲图的提取同山脊线与山谷线的提取中一样，由于鞍部点的整体

3、位置是处于山脊上的，需要提取出正地形以舍弃那些在负地形上的错误的提取结果。正地形的提取过程与完全一样。提取过程分别是：利用11 x 11窗口进展提取平均值的邻域分析，结果为 meandem图1；原始DEMf meandem目减图2，并以0为界进展重分类，大于0的属性重新赋值为1,小于0的赋值为0, 结果命名为zhengdixing ,图3、4所示。利用 spatial analyst菜单下的 surface analysis 菜单中的 contour 和 hillshade 工具分别提取研究区域的等高距为40mB勺等高线数据ctour 图5和研究区域的晕渲图hillshade图 6。图1用11

4、 x 11窗口进展提取平均值的邻域分析睡 Raster Calculatorpayers;ArithmeticAbs毗1LeilFkafFlooiIsNJT riflonametricSinASin 1CasACjsTanAT anPowers写qildern-meandem|Logarithms匚片 LogEk口2 I I m2ExpW LooWScjRowAbout Bidding EipressicrF. ScuEiJ SMadm EET! Q CiX 皿 1MLM111*如中 EM。| Ui嚏，必匚 式-=>图2原始DE% meander目减图3以0为界进展重分类£|?

5、hl h J.工1巾£pa.ti aJ *.上门>:*H图4重分类的结果Enputwifm卜SCEn M >Hncir! InpH.th: range-Zmih: L-f9X.fi£HCFi5i ZniK： ZJTf.OCCTWMTfe 比用EKt DOTC&LF：0川祺15】XZ i&EiljELiET b njJ »3jt difin UorHi hrjtnmHbc.;L弓笛la通num Aibif!Z3£lTdtt !lll»LAf 1 3 g tdltb：白LWedldb# Ltlor .5旧 金llllll

6、lL 3*1图5等高距为40m勺等高线数据ctourLajm sB- 口 k o+ 卜 * iH Lm . a “LmjE if CllErO«tLW& tirAis T 4m图6研究区域的晕渲图hillshade2、山脊的提取山脊的提取与练习1脊的提取过程是完全一样的，分别是进展洼地填充、然后在无洼地 的DEMt提取水流方向、基于水流方向计算汇流累积量数据、提取汇流累积量数据等于零的栅格。提取过程产生的各个数据分别为：filldem 图7、flowdir 图8、flowacc 图9以与 flowacc0 图 10、11。图7洼地填充图8在无洼地的DEMt提取水流方向图9基于

7、水流方向计算汇流累积量数据褪 Rast ex Calculatorpayers;detn filkisrnflov/accflwdirdefnHUkMdo of dem meandeinn meandemPeDla&s of Calculetioi vArithmeticAbslrt 1CeilFkafFlooiIsNJflcwacc - 0Loganrims口 pLogEk口2 |Log?ExpWLooWT rigomrnetricSin力 Sin 1Ca$ADjsTanAT arPapersjdGqr |PowAbout Bidding Eipresskrs£valjat

8、& Cancel (<图10提取汇流累积量数据等于零的栅格* 口 中 LMrImD-l料it i.Fn：3'CLcrij-rwTlLH MiLiaq if ql0Bec1«99 4-t LiLnla±a uQucECh" S'iIkTi口海？0汪 j-HKfcin |k # 臂/ M-l i "F tr-f 口*±” Hl皿ut I-?«La k” 通联 r«cicIjAibE I M I iLlL £a ATifalp'Vl Tib Lj v r liEcvm?Tq; t p

9、vIt Ulit TM11H. tldtH4cl皿 M O'bclEI«JLxKr Twit心!ifa£ Ar.yti Tu£ e C&IUvbd.DeMi>：F Dl.EiUIB«Ei*r««Uw向1"一3 C3.-& lltmiiii M it1 Hp*WLifv / IfeLJ *Z rui产 Tl'X iXERj E LOOi> 115 iarpfili-HiA 3FLw Ltn/k，SLrJEhiliif>巾由 Str-iw Lcaik lt>M to Aw

10、5g iv 1+! Ibl4f El»4JWi:P W! I.隼 | TH”串4 q比 Eui-S图11山脊线3、山谷的提取山谷的提取也与练习1谷的提取过程是完全一样的，分别是基于原始DEW算出反地形DE瞰据计算中是利用原始 DEIM去常数3000、基于反地形DE瞰据提取水流方向数据、基于水流方向数据进展汇流累积量数据、提取汇流累积量数据等于零的栅格。提取过程中产生的数据分别为：fandem图 12、13、flowfdirfan 图 14、flowaccfan 图 15以与fanfacc0 图 16。图12基于原始DEW算出反地形DE瞰据¥ Jf-HAntgll- rUwQ

11、D 080 PIL«W tc< Ih g44Vvt3hi £ rdcdv*T4vtf ini-M f-ewJte fgu3,''> |iT f I* mJ I!图13反地形DE瞰据图14基于反地形DE瞰据提取水流方向数据Fl-sr jt-rrvsixl *t ihh'FiqiiJ fl n iar-trii n- rnR vrQft-E liFt iJBswUliiUu. uLsib、匕11、匚11 rtLL wEimr IdfiCh IiHLJH-*加 I；|Lil7：IFLME图15基于水流方向数据进展汇流累积量数据rp.i.( Ji

12、M-jr4fi.iI PtKii fY；nB-iF-*iLTr 1i* *| 1呷01.I CiHt.tA'iJ.i.uaiL AciiLR f*I Lipir«- hX"Giw Ivrii 3皿i lki-rfcT-iiiXiI kNI Fck»k-Unx T-iIj |Xi, F2,I-hAhH'hwqo llw-HA ilwrdxrHi hU4«Ia)好1 « *wIwa*wm Au r 1露悭t 一” 1 n p -5SUrFw.dl. 一irh二Gaai.JLkJll一 A3IIw"il"r-Mhn

13、rLn民d出"曲l>hl fersu f密 Ret«E CAlculatoz？ X中：Arthm,eTngohorr倒 nehldem llauo: HoiacoO1 .Le47 jjjlIjJ _d _d_d _dOdd/L=<=-MdJr生Nor皿CelFlmLevari >节噌-H 1Powers由匚g。IAlar libtaccrarHoKidiihnHikkoisaf dem mcaiderin<,i>ErR成ExjIOLd>3L哈-LI州Sqil£甲Puwllowaccfanl = 1生cUt EEuldinj Ep

14、resacreCaricd图16山谷线4、鞍部点的提取(1)将提取出的山脊线数据 flowacc0和山谷线数据fanfacc0利用spatial analyst 菜单 下的raster calculator 的工具进展相乘运算，如图 17,运算结果命名为anbuqu图18。图17 raster calculator工具进展相乘运算图18山脊线数据和山谷线数据的相乘结果anbuqu(2)将上一步中提取出的数据anbuqu和正地形数据 zhengdixing 利用spatial analyst菜单下的raster calculator的工具进展相乘运算，就得到了鞍部点的栅格形式数据，命名为rasteranbu图 19。0 iii三 二一?£门 Q 4PIM*ZlWrfm llriM- fLwvrdi Uzi flkrirU3dd 印】。 miiUi >r z veinL* rrbrnnl'F! ur 6 i f alvSl rii wmlw3牍面 IIa 口口图19鞍部点的栅格形式数据图20重分类数据(3)将栅格数据rasteranbu进展重分类，所有0值和NODAT徽据赋为NODA