
1)将各个图层综合到同一个图层下,生成土地利用图
2)然后增加一个属性列F(土壤有机质的含量计算)
将道路、居民地、水系、未利用地、滩涂的F值设为2.5,
水浇地、旱地、菜地设为2.75
灌溉水田为3.25.
天然草地为3.75
有林地、灌木林、疏林地设为4.5
3)将图层栅格化,这样每个栅格都会有一个F值
4)利用栅格计算器,计算K值
公式:K=exp((0.4681+0.7663*F))*ln(E)+2.884-8.1209*F
E=916+331*log10 *i
i为将30min最大雨强,这个数据老师会给。
●L值的计算:
L=(λ/72.6)m
λ表示坡长(用Feet或m表示的坡的水平投影长度),m的值与用用百分数计算的坡度有关,依据相应的坡度制定m的值:
| 坡度 | 大于5% | 3.5%-4.5% | 1%-3% | 小于1% |
| m | 0.5 | 0.4 | 0.3 | 0.2 |
S=10.8*sinθ+0.03 坡度大于9%
S=16.8* sinθ-0.5 坡度大于等于9%小于18%
S=21.91* sinθ-0.96 坡度大于18%
●P值代表认为因子
P值等于“1”表示没有人为因子影响,
| 坡度(precent)100% | P值 |
| 0-2 | 0.6 |
| 3-8 | 0.5 |
| 9-12 | 0.6 |
| 13-16 | 0.7 |
| 17-20 | 0.8 |
| 21-25 | 0.9 |
老师:论文中的方法不适合我们用,找到这样一个方法,但是NDV I不好求,没找到方法,不知道老师能不能找到方法。
先提取NDVI再用下面公式计算
C= ( NDV I - NDV I soil) / ( NDV I veg - NDV I soil)
(论文:基于ERDAS的矿区植被覆盖度遥感信息提取研究_以陕西省榆林市神府煤矿区为例.)
●R值的计算(降雨影响因子)
1、E=916+331*log10 *i,(i为将30min最大雨强)
R=(∑Ei)/N
∑Ei表示所有的降雨值,n为年数
(论文:MIKE-SHE与MUSLE耦合模拟小流域侵蚀产沙空间分布特征)
最后的A = R×K×L×S×C×P
