利用INVEST模型、ARCGIS计算产水量、水源涵养量教程、数据来源、数据处理方法（一）

PS: 水源涵养整体比较复杂，需要找的数据和计算的数据量也比较多，本人也是参考了很多大神的视频，然后整理了一些自己做的步骤，随手记一下。

主要参考了B站上面的视频1 InVEST模型产水量模块（Annual Water Yield）数据获取与操作_哔哩哔哩_bilibili

一、数据下载获取

1.1 invest软件下载

下载链接InVEST | Natural Capital Project (stanford.edu)​​​​​​   下载后安装打开 annual water yield模块

1.2 数据获取

1.2.1 降水量数据

来源：国家青藏高原科学数据中心 中国1km分辨率逐月降水量数据集 可直接下载国家青藏高原科学数据中心 (tpdc.ac.cn)

1.2.2 蒸散量数据

来源：国家青藏高原科学数据中心  中国1km逐月潜在蒸散发数据集（1901-2022）

国家青藏高原科学数据中心 (tpdc.ac.cn)

1.2.3  根系深度

1.2.4 植物可利用含水量

利用中国土壤数据库里的数据计算

世界土壤数据库（HWSD）的中国土壤数据集(v1.1)从国家青藏高原科学数据中心下载。国家青藏高原科学数据中心 (tpdc.ac.cn)

计算公式为：

1.2.5 土地利用数据

中科院地理所的30m分辨率的土地利用数据 可免费获取

1.2.6 流域数据

自己百度下载的

二、数据处理

写在前面：1，所有空间输入必须具有完全相同的投影坐标系(以米为线性单位)，而不是地理坐标系(以度为单位)。所以进行数据处理之前要保证各个因子均具有相同的投影坐标系（线性米单位）。
              2，栅格输入可能有不同的单元大小，它们将被重新采样以匹配土地利用数据栅格的单元大小。因此，所有模型结果将具有与土地利用栅格相同的单元大小。

2.1 降水量蒸散量数据处理

降水量与蒸散量数据处理方法相同，所以只介绍降水量处理方法

降水量要求为年平均降水量，单位:mm/年，因此需要下载该年份前后一年或者两年的数据取平均值。如：目标年份为2020年，则分别下载2019、2020、2021年的降水量，或者18/19/20/21/22年的降水量。 

1.2.1中介绍的网站上下载的数据为nc格式，需要用arcgis中工具“创建NetCDF栅格图层”转换，输入nc文件，参数默认，波段维度选time，点击确定。

创建好后，双击图层，选择符号系统，可以看到共有36个波段，我这里下载的是15-17年3年的所以是36个月，共有36个波段（如果是一年的话就是12个波段）。

然后进行裁剪（按掩膜提取工具）和投影（定义投影工具），像元大小为0.0083是没有投影过的，投影栅格时需选择相应的投影坐标，特别注意输出像元大小改为1000，这样就不用重采样了（后均统一）。

投影好的数据波段数36，像元大小1000*1000，证明处理好了。

至此预处理工作已完成，然后利用栅格计算器求年均降水量，得到年均降水量栅格。

这个的计算有三种方法：（1）打开.nc降水气象数据制作年降水量地图_ncl可以画降水分布图吗-CSDN博客根据这个链接里的方法将单波段的数据导入进来，用spatial analyst-局部分析-像元统计数据工具求sum，再用栅格计算器计算。这个比较适合只有一年的数据，不会太多。

（2）使用波段算术函数，窗口--影像分析--选择要计算的影像-----添加函数----右键点击影像栅格---插入函数---波段算术函数

方法选择用户定义，输入表达式，这里是36个月求和，计算出三年的总降水量。表达式为：B1+B2+B3+B4+B5+B6+B7+B8+B9+B10+B11+B12+B13+B14+B15+B16+B17+B18+B19+B20+B21+B22+B23+B24+B25+B26+B27+B28+B29+B30+B31+B32+B33+B34+B35+B36

得出总和后就只有一个波段了，然后用栅格计算器除以30，意思是求三年的平均值，除以10是单位为0.1mm，化成mm要除以10。

最终结果为单波段的年均降水量1000多mm，符合实际。完毕。

---

蒸散量需要的是年度蒸散量，不需要取平均值，因此只需要将12个月的波段加到一起再除以10就可以。利用像元统计工具求和，再用栅格计算器除以10得出蒸散量。

2.2 根系深度处理

根系深度用的基岩深度数据代替，数据是网上下载的

这个数据就是tif格式的，分辨率1km，投影裁剪一下就可以直接用 比较简单 

2.3 植物可利用含水量处理

这个主要是处理土壤数据库的数据，教程在这个链接 可参考基于GIS简单处理世界土壤数据库（HWSD）的中国土壤数据集-CSDN博客

处理好后 利用1.2.4里的公式在Excel里计算 

2.4  土地利用数据处理

这个数据由于使用的是数据产品 只需要裁剪一下 重采样一下分辨率 投影下就好

2.5 生物物理参数表

根据说明，下载这个表的样例，然后进行修改就可以了，根据提示进行相应的修改。

要注意的是这个表的格式需要为UTF-8   csv文件 ，如果不是的话用记事本打开另存一下。描述那一列的内容最好用英文 不然容易报错

2.6 Z系数

Z系数与产水量成反比，这个值是一个经验值，需要出了结果后根据统计量修正。初次做的话先填1.5默认值。

2.7 流域数据

这个是用研究区范围代替  自己准备就好

2.8 研究区范围数据

也是现成的三级流域矢量面数据  裁剪投影一下就好

---

---

这个数据处理好后最好是重新导出一下，让这些数据拥有相同的范围、分辨率。便于导入到invest模型。

数据处理好后，就是把处理好的导入到invest里面去跑就可以了。

跑完后，找到结果文件夹，需要的结果就是wyield文件，

这是这个字段的解释，就是我们需要的产水量结果，单位为mm

跑出来的结果要看产水量的均值，根据这个值去调整Z系数。

首先看降水量的均值，计算产水系数=产水量/年降水量

根据研究区的水资源公报调整。直接去百度搜索研究区的政府网站，下载相应的水资源公报。

例如：海口市的产水量/降水量=23.04/45.66=0.505就是产水系数

所以计算出来的结果也要保证产水系数接近于0.505，如果invest结果计算出的产水系数为0.3，那就调高产水量，进而降低Z值 （产水量与Z成反比）再去计算一遍 直到跑出来的结果合适。

就做完了