自-DEM实验-基于ARCGIS的等高线实验.doc

实　 验 　二实验代码：081440213-２实验名称地形指标与地形特征信息的提取实验目的学习并掌握在AｒcGIS中提取地形指标（坡度、坡向、坡度变率、坡向变率、地形起伏度、地面粗糙度)与地形特征信息（山脊线、山谷线、鞍部）的方法和步骤实验步骤1. 坡度和坡向的提取导入 ｄem1图像对ｄem1求坡度 Ｓｐaｔial　anａlysｔ→surfａcｅ anａlｙst→sｌｏpeﻩﻩ得到sｌｏｐe of dem1对ｄem1求坡向 Spaｔiaｌ analyst→ｓuｒface aｎaｌyｓt→aspecｔ ﻩ得到Ａspｅcｔ　ｏf dem１2. 坡度变率的提取坡度变率就是坡度之坡度 对Slope of dem1再求坡度变率 Spatｉal　aｎａｌysｔ→surｆａcｅ analyｓt→sｌope得到sｌope of　slｏpe ｏf　dem１3. 坡向变率的提取坡向变率就是坡向之坡度(一) 反地形求坡向之坡度Sｐatｉal aｎaｌyst→Rａster calculatｏr输入1１５4 －　[dem１]，得到Calcuｌation即,反地形二) 再对反地形求坡向之坡度先求反地形的坡向1) 对Calcｕlaｔion求坡向 Spａtial anaｌyｓｔ→sｕrfaｃe　aｎaｌyst→aspeｃt得到Aspeｃt oｆ Calcuｌation2) 再求Aspｅcｔ ｏｆ Cａlculａtｉｏn的坡度　Ｓpatiａl anaｌyst→suｒface analyｓt→slｏｐe得到Sｌope　ｏf Aspｅｃt of Caｌcｕlation3) 再对Ａspeｃｔ oｆ dｅｍ1求坡度　Spatｉaｌ　ａnalｙst→surfacｅ　analyst→sｌｏpe得到Sloｐｅ of Aｓpeｃt of dｅm14) ＳOA ＝ ([Sloｐe ｏf Ａsｐeｃｔ of　ｄem１] +　[Ｓlope of　Aｓpecｔ of Calｃuｌation　])- Abs([Slopｅ ｏｆ Aspeｃt oｆ dem1］ - [Sｌｏpe　of　Aspｅct oｆ Calcｕlaｔioｎ　］) /2即可得到坡向变率坡向变率SOA values 　　Ｈigh：1７５.435 Low:０4. 地形起伏度的提取先求邻域Spａtｉal　analysｔ→neiｇhboｒ statｉｓtiｃs 其中Sｔatisｔic tｙｐe 设置为Mａxｉmum Neiｇhｂorhooｄ Seｔting设置为11１1 得到NbｒMax of　dem１其中Statｉsｔic　tyｐe 设置为Mｉｎimum　　 Neiｇhboｒhood Setｔｉng设置为1１１１得到NｂrMiｎ ｏf dem1再在ｒａｓｔｅr cａｌcｕlａtoｒ　中［NbrＭax　of dem1］ 　－ [NbrＭin ｏf　deｍ1]得到的Caｌculation即为地形起伏度地形起伏度5. 地面粗糙度的提取对dem1求坡度 Spatiａｌ　analyst→surface aｎａlyｓt→slｏpe得到Slope　ｏf dem1再在ｒaｓｔer　ｃalｃuｌaｔｏr 中1 ／ coｓ([Sｌoｐe of deｍ1] *　 3．141５９2654　/ 180)得到Caｌculａtion2即为地面粗糙度地面粗糙度6. 山脊线的提取先导入ｄem11) 再对ｄem1求邻域ﻩSpatiaｌ anａlysｔ→neigｈbor stａtｉstｉcｓ Inpｕt data： ｄeｍ1 只修改neigｈbor seｔｔｉnｇ中的heｉｇｈt、ｗｉdth均改为1１2) 求正地形 在rasteｒ　calｃuｌａtor中 ［ｄem１]　– [ｍeａｎｄｅｍ]求得ｃalｃulａｔｉon再rａｓter cａlcｕｌatｏr中coｎ([Cａlｃｕlatｉon ] > 0 , 1 ,　0）　 命名为zｈenｇdixiｎg3) 求负地形 同理,正地形是[Calculatiｏｎ　] ＞ 0的,而负地形就是[Calｃｕｌaｔion ] < 0在raｓter　caｌｃｕlａtor中cｏｎ（[Calcuｌａtion ]　<　０ ,　1 , 0) 　命名为ｆudixｉng4) 求fill填充洼地Ｓpaｔiａl ａｎａlyｓt　tｏols→ｈydrｏlogy→fｉｌlＩｎput　:C:\deｍ02\dｅm１Ouｔｐｕt： fill_demfｉｌl＿ｄem5) 求流向fｌow　dirｅctiｏｎSpatｉal 　ａｎalyst toolｓ→ｈｙdrｏｌoｇｙ→flow dirｅｃtionInpuｔ : fｉlｌ_deｍOutput： ｆlowｄｉr＿deｍ6) 求流量 floｗ aｃcｕmuｌａtionSpａｔial　　anａlyst　tools→hydrology→flow accumｕlatｉonＩnｐｕt :　fｌowdiｒ＿demOuｔput：　ｆｌowaccFlowａcｃ7) 求facｃ在ｒasteｒ calcuｌaｔｏr中facc０　= (［flowａｃc＿2］ =＝ ０)得到ｆacc０8) 对fａｃc0求邻域 在Sｐatial anａlyｓt→neiｇｈboｒ sｔatｉstics中ﻩﻩInｐut daｔａ: facｃ0 ﻩOuｔpｕt raster：　C:\dem02＼neｉｇhｂorfaｃｃ０9) 对dem进行阴影变化ﻩ在Sｐatiａl　aｎalyst→surfaｃe anａlｙsis→ｈilｌｓhａｄｅﻩInput surface: deｍ1ﻩOutput：　hｉllshaｄehillshade10) 对neigｈborfacｃ0进行透明处理ﻩ进入nｅighborfacc０的属性→disｐｌaｙ→transpａrentｓ 改为30%ﻩ再将neｉghｂorfａｃc0放置在ｈilｌsｈaｄｅ之上 便可以看到山脊线的大体分布11) 对neiｇhbｏrfａcc0进行重分类 在Spatｉal ａnalyst→ｒeclassｉfy 再Cｌassify 先将ｃlassｅs改为2，否则若现修改meｔｈoｄ,会无法修改 在method改为maｎual手工修改 此时会在右侧看见.５ 和 1ﻩ在将valueｓ修改为0-0.5５41ﻩ得到ｒeneｉghbｏr 12) 求山脊线在ｒasｔｅr caｌcuｌator中Reneiｇhｂor*zhｅngdixing再将属性中除了ｏｌd valｕｅs为1的其余属性值的new　ｖalueｓ设为nｏdａｔa7. 山谷线的提取基本步骤与求山脊线相似需要改动的是：求反地形ｆａndixiｎg11５4 - [ｄeｍ1]　在Spatｉal aｎaｌyｓt→ｒecｌａssｉfy将valｕes修改为０-０．6５677可以发现山谷线与流向是一致的8. 地形鞍部的提取Anｂｕ=shanjｉxiａn*shanguxａiｎRanbｕ=ａnbu*zhengdiｘｉng。