《地形指数计算》由会员分享,可在线阅读,更多相关《地形指数计算(5页珍藏版)》请在金锄头文库上搜索。
1、关于在ArcGIS中地形指数的计算,网上能找到不少的文章,但是不容易找到具有的方法。 这是一个自己摸索的方法,不知道对不对,欢迎讨论指正。基本公式:tan p原始数据:DEM栅格过程数据:洼地填充后的洼地流向flowdirfill和汇流累积量flowacc。二者计算过程参考汤国 安地理信息系统空间分析实验教程第十一章.1 计算 SCA(specific catchment area)SCA=CA/flow widthCA= (flowacc+1) *单位栅格的面积(分辨率人2) (+1很重要,在修正这个之前,忽略加1 导致最后结果出现负值。)(1)计算汇水面积CA。通过hydrology模块,
2、计算出填充洼地后的flowacc*分辨率平方得到汇水面积CA国a昌黑昭Cl x(2)确定流向宽度flow width。对于D8算法而言,当流向为格网对角线方向时,流向宽度为sqrt (2) *格网间距L;东西 南北四方向的,流向宽度为格网间距L。经洼地填充后生成的Flowdirfill即8个流向,用con条件语句进行计算,指定东西南北方向 的格网间距为L,对角线方向的为sqrt(2)*L。得到流向宽度。聲 Faster CalculatorKTc-ols Pio *E * O M4Power:wTngon&nirfHccan(flawdufdl =- 1詁5(1郴M训=-ISSJAowdirf
3、iU=佰齧, Rowdi阿=62SH|flcwdutkll=讹砕岡叫 2H(hw血和躺 &2S歧寧雌|4flfflwJRfiH = 辽西哼举甲期佃询何=s 13SJ5*Sqr1Abodl fiuildng km&旳&raEvaluateCanotiArtshmeticExample指定分辨率为25米,条件语句如下:con (flowdirfill = 1,25,flowdirfill = 4,25,flowdirfill = 16,25,flowdirfill= 64,25,flowdirfill = 2,25*Sqrt(2),flowdirfill = 8,25*Sqrt(2),flowdi
4、rfill= 32,25*Sqrt(2),flowdirfill = 128,25*Sqrt(2)(3)计算 SCA。将汇水面积除以流向宽度得到单位汇水面积。2计算slopes:.s 0:0:0:田:田:田:田:田:田:田:.E.用arctoolbox里面的slope功能直接对原始数据进行坡度的求取。注意得到的坡度单位 是。degree- a:oerver ioo-Spat 1 al Analyst ToolConditionalDensityBist ance ract ronCeneralrzat ion.Crcundwat erHydrologyInterpolation.LocalMa
5、p Algeb匚nMathJrtultivariat eneighborhoodOverlay.Kaster Great ionlieclasESolar Radi at ionSurface护 Aspect护 Corit ou.r护 Coritou.r List护 匚口rvat u.r e护 CutZFill参 Hillshade参 Observer Points护 Slope参 Vieu/shed+ 觀 Zonal3计算地形指数在这里说明一下。生产坡度的时候单位是度。但是在栅格计算器中的tan函数默认值为弧度。 所以直接tan(slope)会出现问题,需要将角度的坡度转化成弧度进行计算。公
6、式如上图。用 spatial analysis tools 里面的 raster calculator 工具进行计算得到地形指数值,从porperties中可以得到地形指数的最小值最大值平均值等统计指标。遇到error000539的处理方法:VBScript (the default ETP/V)01 . gp. Calcul ateFiel d_rna h age me n+ (nn_fi 1 e, p_rB3 ?二_山荷幻 / LB00 T )Python:0: . EP. Calcul ateFJe d.manaEeine n 十(rin_fig, p.-rea, !=績川 / Fl oat (1W0), PYTHON)