上一期关于流速系数的计算,我采用的是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,今天的分享就到这里了,整个水源涵养的计算就结束了,我们下期开始新的生态系统服务功能的计算!让我们下期再见!!!