基于ARCGIS计算RUSLE中LS因子值
LS因子针对不同的研究区,计算方法不同,本文参考别人的文献,选择了计算方法:
另外还有很多计算方法
1、首先对DEM数据进行填洼处理:arcgis中fill工具。
之后计算坡度:坡度输出单位直接设置为degree就行。但在接下来的三角函数运算中,度要转换为弧度。
2、使用栅格计算器,计算S因子:
Con(“slope.tif” <= 5 ,10.8 * Sin(“slope.tif” * 0.01745) + 0.03 , Con(“slope.tif” <=10 ,16.8 * Sin(“slope.tif” * 0.01745) - 0.50 , Con(“slope.tif” <= 25,20.204 * Sin(“slope.tif” *0.01745 ) -1.2404,29.585 * Sin(“slope.tif” *0.01745 ) -5.6079 ) ) )
这里乘以0.01745的目的是将角度转换为弧度,进而计算sin值
3、进行流向分析:
在流向分析和流量分析时,使用arcmap速度很慢,尤其是数据量很大的时候,推荐使用正版的arcgis pro。申请一个试用版也许。
流向分析也是只设置出入和输出就行。
4、流量分析,也是只设置输入输出就好。
5、Lambda的计算:lanbda使用流量栅格分辨率,我这里是使用的带投影的。分辨率大约30m。
6、计算L
我这里把5和6合并在一起计算了。
Con(“slope.tif” <=1 ,Power(“flowacc.tif” * 29.4481461287987 /22.13 ,0.2) , Con(“flowacc.tif” <=3 , Power(“flowacc.tif” 29.4481461287987 /22.13 ,0.3) , Con(“slope.tif” <=5 , Power(“flowacc.tif”* 29.4481461287987 /22.13 ,0.4), Power(“flowacc.tif”* 29.4481461287987 /22.13 ,0.5)) ) )
。
也可以先使用con函数计算M,再使用power(flowacc,m)进行计算
7、计算LS=L*S
对结果重新设置样式
在运算过程中,尤其是数据量很大的情况下,尽量放在根目录下面,我不知道为啥目录很深,也会运算不出来结果。
对于结果,要根据不同人的研究,去进行进一步的验证与纠正。
基于ARCGIS计算RUSLE中LS因子值相关推荐
- 基于Python(sklearn)计算PLS中的VIP值
基于Python(sklearn)计算PLS中的VIP值 sklearn中PLS回归模型并没有计算VIP值的方法,但VIP又是很重要的筛选变量方法.下附代码思路与完整代码,若有错误,万望指正. 1.首 ...
- 计算生活中常见的Π值
目录 前言 一,原理 二:运行过程 三:整体代码 前言 Π值是很一个神奇的数,从我们上小学的时候就开始接触他了,而现在依靠计算机强大的算力,能够计算出∏值后面的几百万位甚至几十亿倍.Π值对于我们今天的 ...
- lambda :: 和计算list中某个字段值的总和
用lambda表达式计算一个list中的某个字段值的总和 int total = List.stream().mapToInt(entity::getDistance).sum(); 解释:List是 ...
- OpenCv案例(二):基于OpenCVSharp计算图像中的圆的周长和面积
需求:求下图中圆形的面积和周长: 原图如下: 处理后图像: 计算结果: 代码如下: public Mat GetObjectByImg(Mat src,out double retArea, out ...
- python取出数组大于某值_计算矩阵中大于某个值的所有值
有很多方法可以实现这一点,比如扁平化和过滤,或者简单地枚举,但是我认为使用Boolean/mask array是最简单的方法(iirc是一个更快的方法):>>> y = np.arr ...
- 多因子模型之因子(信号)测试平台----计算因子值
广告:本人的单因子测试视频教程https://edu.csdn.net/course/detail/25572 近一个半月疯狂的接触多因子模型,其中对于单个因子的回测,是最熟的.而对于单个因子,或者叫 ...
- TensorFlow中的Nan值的陷阱
北京站 | NVIDIA DLI深度学习培训 2018年1月26日 NVIDIA 深度学习学院 带你快速进入火热的DL领域 阅读全文 正文共1583 ...
- oracle sql为null值,在SQL Server Oracle MySQL当数据库中查出某值为NULL怎么办
在SQL Server Oracle MySQL当数据库中查出某值为NULL怎么办? 1.MSSQL: ISNULL() 语法 ISNULL ( check_expression , replacem ...
- python中计算如何实现_基于python如何实现计算两组数据P值
我们在做A/B试验评估的时候需要借助p_value,这篇文章记录如何利用python计算两组数据的显著性. 一.代码# TTest.py # -*- coding: utf-8 -*- ''' # C ...
最新文章
- SAP QM 执行事务代码QS23为检验特性分配Selected Set的时候报错 - You cannot use entries from catalogs 1 and 3-
- 深度:Leap Motion手势识别大揭秘
- 物联网的兴起与二维码的前景
- SpringBoot使用RabbitMQ消息队列
- 20179214 2017-2018-2 《密码与安全新技术》第七次作业
- android gradle自定义钉钉群提醒
- 最近有粉丝向我请教Java泛型,于是重温了一下,希望学弟面试顺利
- 例子---PHP与Form表单之一
- hdu3535 (分组背包,最少选一 + 最多选一 + 随意)
- 利用DB Link搞定Oracle两个数据库间的表同步
- windows防护之(一)屏蔽危险端口
- echarts 自定义 tooltip
- 阿里云的服务器居然泡在“水”里?| 数据中心参观有感
- iPhone 如何设置自定义铃声 (更新至iOS13 and macOS Catalina)
- Aip接口自动化测试框架pytest+allure+request+jsonpath+excle
- Generative Adversarial Nets[AAE]
- Tensorflow中令人困惑的点(1)
- vulcan 编程_Vulcan和AUSK之间的比较:如何充分利用Node,React和GraphQL
- c语言二维数组a中,a,a[0],a[0][0]的值与值的类型
- Qgis教程07:矢量数据属性编辑
热门文章
- 交互式二维图形编辑系统 课程设计 代写
- 计算机辅助诊断(Computer-aided diagnosis)小总结
- C#初探-抽学号软件
- cadence之贴片电阻封装绘制
- UKF 无迹卡尔曼滤波
- Docker 报错Error response from daemon: Get https://xx.xx.xx.xx/v2/: Service Unavailable
- 第6关:人民币美元汇率双向兑换计算
- 用centOS 7安装cadence搭建适合IC Design的科研环境
- java编程实例_几个java编程的典型例子
- awk、print使用