降雨插值_ArcGIS计算土壤侵蚀模数（二）降雨侵蚀力因子R计算

本次采用中国土壤流失方程CSLE计算土壤侵蚀模数，计算公式为：

A=RKLSBET

式中，式中：A—土壤侵蚀模数。t·hm-2a-1；R—降雨侵蚀力因子，MJ·mm·hm-2·h-1·a-1；K—土壤可蚀性因子，t·hm2·h·hm-2·MJ·mm-1；L—坡长因子，无量纲；S—坡度因子，无量纲；B—植被覆盖度与生物措施因子，无量纲；E—工程措施因子，无量纲；T—耕作措施因子，无量纲。

本次分享R因子计算过程，采用ArcGIS计算R因子首先需要降雨量数据，研究区降雨量插值计算详见第二篇文章。

1.降雨侵蚀力因子R计算公式

R因子计算根据所获取的降雨资料分别采用不同计算公式：

（1）当获得逐日降雨资料时，按如下公式计算：

式中：

Rr—第r日的降雨侵蚀力因子，MJ·mm/（hm2·h）；

Pr—第r日的日降雨量，mm，需单日降雨量≥12mm，否则按0计；

α，β—计算降雨侵蚀因子统计系数，无量纲；

Pd12—日降雨量≥12mm的日降雨量，mm；

Py12—日降雨量≥12mm的年降雨量，mm；

获得逐日降雨侵蚀力因子后，可累加得多日、月、年降雨侵蚀力因子。

（2）当获得月降雨资料时，按如下公式计算：

式中：

Rm—第m月的降雨侵蚀力因子，MJ·mm/（hm2·h）；

Pm—第m月的降雨量，mm。

获得月降雨侵蚀力因子后，可累加得多月、年降雨侵蚀力因子。

（3）可获得年降雨资料时，按如下公式计算：

式中：

Rn—年降雨侵蚀力因子，MJ·mm/（hm2·h）；

Pn—年降雨量，mm。

（4）可获得多年平均降雨资料时，按如下公式计算：

式中：

Rd—多年平均降雨侵蚀力因子，MJ·mm/（hm2·h）；

Pd—多年平均降雨量，mm。

由中国气象数据网下载数据为累年月值降雨量，可分别计算12个月的R因子和多年平均R因子，本次只计算多年平均R因子做示范，按照公式（4）计算

2.降雨侵蚀力因子计算

打开插值后的降雨量数据，打开Spatial Analyst Tools—地图代数—栅格计算器工具

输入计算公式

点击确定，即可计算得出内黄县R因子图层。

好了，降雨侵蚀力因子计算完成，R因子计算很简单，重点是降雨量数据的插值计算，需要注意的是不同的降雨资料选用不同的计算公式，如果学习过程中缺少行政区边界数据，可留言需要的行政区，后续我也会分享如何根据影像数据自己制作区域边界。