生态系统服务功能之水源涵养（五）

上一期关于流速系数的计算，我采用的是GEE实现的，我发现用arcgis直接对土地利用数据进行重分类更方便。

1、土地利用类型重分类

 确定，结果即为流速系数V，这时我们发现流速系数为整数型，需要将其转换为浮点型，方便后期进行水源涵养的计算。

2、整型转浮点型

 确定，运行结束后即转换成功。

我们再来看一遍水源涵养的计算公式：

 𝑊𝑅=𝑀𝑖𝑛(249/𝑉,1)∗𝑀𝑖𝑛(1,0.9∗𝑇1/3)∗𝑀𝑖𝑛(1,𝐾/300)∗𝑌

其中：WR为水源涵养量（mm）；V为流速系数；K为土壤饱和导水率；Y为产水量；T1为地形指数。

3、计算𝑀𝑖𝑛(249/𝑉,1)

我们先计算249/V：

再计算 𝑀𝑖𝑛(249/𝑉,1)，用con函数，输入公式：Con("V_249_2020.tif" >= 1,1,"V_249_2020.tif")

公式Con("V_249_2020.tif" >= 1,1,"V_249_2020.tif")表示 V_249_2020.tif的值大于等于1时，就取1，即得到了 𝑀𝑖𝑛(249/𝑉,1)。

4、计算𝑀𝑖𝑛(1,0.9∗𝑇1/3)

我们先来计算0.9*T1/3:

 结果为：

 再计算𝑀𝑖𝑛(1,0.9∗𝑇1/3)，栅格计算器：Con("TWI_03.tif" >= 1,1,"TWI_03.tif")

 运行结果：

 5、再计算𝑀𝑖𝑛(1,𝐾/300)

先计算K/300，再计算𝑀𝑖𝑛(1,𝐾/300)

 运行结果如下：

 6、计算水源涵养

结合前期算出的产水量，计算水源涵养："V_249_2020_compare1.tif" * "TWI_03_compare1.tif" * "K_300_compare1.tif" * "wyield_2020.tif"

运行结果如下： 

 OK，今天的分享就到这里了，整个水源涵养的计算就结束了，我们下期开始新的生态系统服务功能的计算！让我们下期再见！！！