生态系统服务功能-土壤保持-土壤侵蚀可侵蚀因子K计算（仅供自己学习记录）

注意单位之间的换算

②土壤可侵蚀因子的计算（生态保护红线划定指南）

注释：《指南》里采用的是中国第二次土壤普查数据（sand,silt,clay依据国际制划分），而epic公式采用的是美国制划分方法，所以要多乘上一个公式转换K值。然而HWSD数据本身的sand,silt,clay就是美国制划分的，不用转换K值，所以用HWSD数据的话是不用乘第一个公式进行K值的转换 ，直接就是Kepic。

a.投影掩膜区域

b.链接数据表格

c.栅格计算；计算数据(选择的是T_XX)

新建字段，对应四个（浮点型）：

注意：空值区域需要删除，不然会存在栅格计算错误。

0.2+0.3*Exp (0.0256* [T_SAND]*(1- [T_SILT] /100)  )

如果算出来是0，用VB脚本：( [T_SILT] /( [T_CLAY] + [T_SILT] ))^0.3也可以

1-0.25* [T_OC]/( [T_OC] +Exp (  3.72-2.95* [T_OC] ))

注：还加一个括号哦1-0.7*(1- [T_SAND] /100)/((1- [T_SAND] /100)+Exp ( 22.9*(1- [T_SAND] /100)-5.51 ))

最后结果是四个数据相乘