基于ArcGIS的土壤流失模型坡长坡度因子提取.pdf

基于ArcGIS的土壤流失模型坡长坡度因子提取 通用土壤流失模型（USLE）和修正的土壤流失模型（RUSLE）是当今研究某流域水土流失状况的通用模型，其适用性在全世界范围内得 到认可和推广。RUSLE的计算公式为：A=R×K×L×S×C×P，其中R为降雨侵蚀力因子，K为土壤可蚀性因子，L为坡长因子，S为坡度因 子，C为植被覆盖因子，P为水土保持因子。本文介绍在ArcGIS中利用水文分析（Hydrology）工具，对流域（区域）的DEM影像数据进行解 译，提取坡长和坡度因子。 准备工作：ArcGIS软件，DEM影像图。 1、  在网上下载流域（区域）的DEM影像图，推荐地理空间数据云（/）选择数据集为GDEMDEM 30M或GDEV2 30M 分辨率的数字高程数据。 注：ASTER GDEM V1原始数据局部地区存在异常，所以由ASTER GDEM V1加工的数字高程数据产品也存在个别区域的数据异常现 象。ASTER GDEM V2版则采用了一种先进的算法对V1版GDEM影像进行了改进，提高了数据的空间分辨率精度和高程精度。 举报 他的主页 他的动态 亲密度 用户/应用/动态 返回个人中心 缄默 一般下载的DEM影像需要拼接和裁剪出研究区域。 注：GDEMV2数字高程数据包含两个文件，dem为数字高程数据，num为质量控制数据。 2、  下一步就是拼接和裁剪出研究区域。在ArcGIS中利用数据管理工具进行栅格拼接和裁剪。 1）拼接：ArcToolbox—Data Management Tools—Raster—Mosaic To New Raster。 举报 他的主页 他的动态 亲密度 用户/应用/动态 返回个人中心 缄默 *注：在此不要选择默认的像素类型，相应的像素类型对应的输出结果如下： 1_BIT—1 位无 符 号整数 。值可以为 0 或 1。 2_BIT—2 位无 符 号整数 。支持的值为 0 到 3。 4_BIT—4 位无 符 号整数 。支持的值为 0 到 15。 8_BIT_UNSIGNED—8 位无 符 号数据类型 。支持的值为 0 到 255。 8_BIT_SIGNED—8 位有符 号数据类型 。支持的值为 -128 到 127。 16_BIT_UNSIGNED—16 位无 符 号数据类型 。取值范围为 0 到 65,535。 16_BIT_SIGNED—16 位有符 号数据类型 。取值范围为 -32,768 到 32,767。 32_BIT_UNSIGNED—32 位无 符 号数据类型 。取值范围为 0 到 4,294,967,295。 32_BIT_SIGNED—32 位有符 号数据类型 。取值范围为 -2,147,483,648 到 2,147,483,647。 32_BIT_FLOAT—支持小数的 32 位数据类型 。 64_BIT —支持小数的 64位数据类型 。 如果选择了默认的8_BIT_SIGNED类型，输出的像元值就会变成0~255。 2）裁剪：ArcToolbox—Spatial Analyst Tools—Extraction—Extract by mask。 研究区域裁剪完毕后，下面开始进行水文分析。 3、  洼地填平处理。DEM是比较光滑的地形表面的模拟，但是由于DEM误差以及一些真实地形（如喀斯特地貌）的存在，使得DEM表面存在一 些凹陷的区域。进行水流方向计算时，由于这些区域的存在，往往得到不合理的甚至错误的水流方向。因此需要对DEM进行洼地填平。 4、  水流方向计算。在栅格图像中，水流方向是指水流离开该格网时的指向。在ArcGIS中，水流方向采用D8算法，即通过计算中心栅格与 邻域栅格的最大距离权落差来确定（距离权落差是指中心栅格与邻域栅格的高程差除以两栅格间的距离，栅