arcgis下遥感图像计算坡长.doc

对坡长的计算采用GIS技术，应用ArcGIS软件进行计算。计算步骤如下：①对DEM数据进行填挖，再对填挖后的DEM数据提取负地形。②求负地形的水流方向并计算流水累积量。③提取负地形的河网，即正地形的山脊线。 可以为通过流量工具得到的栅格指定阈值。定义河流网络系统是操作的第一步。要完成此任务，可使用条件函数工具或地图代数。条件函数常用语法的一个例子是 newraster = con(accum 100, 1)。像元流入量超过 100 的所有像元均将成为河流网络的一部分。 ④对山脊线进行欧氏距离计算，得到近似坡长。⑤将计算的坡长图导入ENVI软件中，采用公式计算出坡长因子L。 对坡度与坡长因子的提取也是在实验小区上进行提取分析的，一般对L因子的提取采用Wischmeier和Smith在1978年经实验提出的的计算方法[13]。公式如下：L=(λ/22.1)^m (3-12)式中：L为坡长因子；λ为坡长；m为坡长指数，其计算公式为：m=β/(1+β) (3-13)式中：β为细沟侵蚀和纹沟侵蚀的比值，其计算公式为：β=(sin?θ/0.0896)/[3.0×(sin?θ )^0.8+0.56] (3-14)式中：θ为坡度；. 发现arcgis计算后发现图像的大小发生变化了…..,从新计算注意大小，并做笔记——计算水流方向hydrology-FLOW DIRECTION-input:DEMoutput: FLOWDIR—–计算水流累积量使用无洼dem重新计算水流方向，生成fdirfill使用此水流方向文件计算水流累积量，hydrology-fill accumulation-INPUT: fdirfillOUTPUT: FLOWACC—-河网的生成;这一步一定要注意输出结果时要根据实际遥感图像象元大小确定输出，不然结果行列号会发生变化;河网生成是基于汇流累积量的，利用ArcMap中的Spatial Analysis分析模块下的Raster Calculator来计算出所有大于设定阈值的栅格，这些栅格就是河网的潜在位置。将计算出来的栅格河网命名为streamnet首先Spatial analyst -Map algebra计算可以为通过流量工具得到的栅格指定阈值。定义河流网络系统是操作的第一步。要完成此任务，可使用条件函数工具或地图代数。条件函数常用语法的一个例子是 newraster = con(accum 100, 1)。像元流入量超过 100 的所有像元均将成为河流网络的一部分。然后:hydrology-STREAM TO FEATURE-INPUT STREAM RASTER:STREAMNETINPUT FLOW DIRECTION RASTER: fdirfillOUTPUT POLYLINE FEATURES: STREAMFEA生成矢量河网文件—-近似坡长计算对河网进行欧氏距离计算，得到近似坡长。SAT-Distance-Euclidean Distance 3D Analyst Tools-Functional Surface-Surface Length