土地利用适宜性（LUS）中的生态适宜性（Suitability Assessment）即生态系统服务（ESs）评价教程

我们的研究利用生态系统服务和权衡综合估价(InVEST)模型来分析研究区的生态系统动态。

Annual Water Yield（年产水量模块）

界面如下：

工作空间：F:\LUS\Invest\Annual Water Yield\xianhai

最好是英文路径，避免报错

降水

首先下载数据，我是从GEE上下载的，代码如下：
https://code.earthengine.google.com/9946915f9793f7fc5eecc8573977387b

将下载的数据至arcgis pro里处理一下
Invest模型需要每一个数据的尺度范围一致，且坐标系为投影坐标
首先进行掩膜处理

然后进行投影变换，同时将分辨率统一，这里设置为250m，下面数据都是这样处理

蒸发量

下载：https://code.earthengine.google.com/7a7d38b6a2557b0cc12e090952bfc7a4
你自己下载其他数据的时候请注意其单位，可以参考Invest的在线指南

下载下来的数据同样arcgis里面进行掩膜、投影

根系限制层的深度

根系限制层深度可以选取土壤深度图代替
这里我用的HWSD_V2_SMU数据库，可以从官网下载，也可以从GEE下载
链接：https://code.earthengine.google.com/9c2f8dcb1a681286548ce5a016fc17d9
从GEE下载的波段如下

这里我们选择波段14进行处理

植物可用水分含量

选择HWSD的波段15进行处理
如果你是从官网下的联合国粮农组织(FAO)和维也纳国际应用系统研究所(IIASA)所构建的世界土壤数据库(Harmonized World Soil Database)
官方网址：https://www.fao.org/soils-portal/data-hub/soil-maps-and-databases/harmonized-world-soil-database-v12/en/
需要处理一下：

土地使用/土地覆盖

用的ESA WorldCover 10m v200
链接：https://code.earthengine.google.com/94fd7e089a8d36e620e36f7b7644c299

生物物理表

根据文献自己查，这里我根据查找的文献得出的表如下：

Z值

这是一个地区性参数，每一个地方都不一样，调整这个参数是作为复杂的，选取合适的参数对结果很重要，可以参看其他论文的设置，最关键的是要学会自己矫正，在视频https://www.bilibili.com/video/BV1PG4y1P73T/?spm_id_from=333.337.search-card.all.click&vd_source=d97cf16c69f59fa1b3a8e5790ca4083d（2：00：01）详细讲解如何调整的步骤。

流域

可以用你的研究区的矢量边界
注意，要加一个字段ws_id代表每个流域的唯一标识符。

运行

收集完成所有数据后，尽量对所有栅格数据进行空间参考，像元分辨率，边界范围进行统一，避免出行错误，这里也是一个十分繁琐的工作。完成这些后可以点击 RUN运行你会得到这样一个文件夹：wyield.tif 就是我们想要的结果！

在ArcGIS中，我们采用分位数分类法将生态系统服务功能划分为三个等级：不重要、中等重要和重要。四个主要的生态系统功能被赋予最高的重要性等级，确立了评估生态系统服务重要性的关键标准。在合并这些功能评估结果时，我们遵循‘遇高则高’的原则，使用‘取最大值’（MAXIMUM）的镶嵌运算符，确保任何一个区域在任何一个生态系统服务功能中被评定为高重要性时，都能在最终的生态系统服务功能重要性评价中得到相应的高等级保护。

分位数分类法分成3类，1为不重要，2为中等重要，3为重要。

Carbon Storage and Sequestration（碳储存和封存模块）

界面如下：

工作空间：F:\LUS\Invest\Carbon Storage and Sequestration\xianhai

现状LULC

使用上面处理后的lulc即可

碳库

同样，根据文献获取，注意要跟LULC的编号对应起来

运行

tot_c_cur.tif是我们要的结果

同样，去arcgis里重分类

Habitat Quality（生境质量）

界面如下：

###工作空间：F:\LUS\Invest\Habitat Quality\xianhai

现状土地覆盖

使用上面处理之后的即可

威胁表格

此表格如下所示，具体数据查阅相关文献

第一列分了两个威胁源，你可以分更多，比如道路，水体等等，最后一列是对应的栅格图像
此步骤我们利用上面的LULC数据重分类得到两个，耕地和建设用地，重分类成0和1的二值图
以耕地为例：
我们LULC下载的是ESA worldcover数据，其属性表如下：

所以我们重分类，将40设为1，其他设为0，命名为gd_c.tif，_c表示现在
如果考虑多LULC情景状态的生境质量，威胁源的命名要注意加后缀以区分：“_b”(基准)，“_c”((当前)，或 “_f”(未来)表明威胁地图的时间段

建筑用地将值为50的设为1，其他设为0

敏感性表格

这里要注意后两列名称要跟威胁源的名称对应起来，参数查文献

半包和常数

根据用户手册，我们选择默认值0.5

运行

其中的quality_c即为最终计算的生境质量值，其余数据为中间结果的输出，将最后结果导入Arcgis中打开并进行重分类

Sediment Delivery Ratio (SDR) （泥沙输移比/土壤保持模块）

界面如下：

工作空间：F:\LUS\Invest\Sediment Diversity Ratio\xianhai

数字高程模型

链接：https://code.earthengine.google.com/bc695301a98287dcfcf671bc8a1fc9ac

数据再arcgis中处理

侵蚀性

也就是所说的降水侵蚀性因子 ®
这个因子我们可以通过年降水量数据计算得到 此 降雨侵蚀力数据，对于此数据的计算有很多方法，这里使用下面公式

年降水量数据在第一个模块中已经下载过了，所以只需要用 栅格计算器计算即可
公式：R_n = 0.053 * Float(“Annual precipitation.tif”) ** 1.655

土壤侵蚀性

下载：https://esdac.jrc.ec.europa.eu/content/global-soil-erosion#tabs-0-description=0

我们用这个进行处理：

如果使用HWSD的话可以参考下面：

士地使用/土地覆盖

使用上面处理后的

生物物理表

参数查文献

流域

前面添加字段之后的

阈值流量累积(像元数量)

默认值：1000
汇流阈值（1000-5000）对结果的影响非常大，结果需要和实测结果校准，可以先用1000试一下

Borselli K参数

反映水文连通性的参数，默认为2

最大SDR值

反映土地上的泥沙最大有多大的概率会被带走，默认0.8

Borselli lC0参数

模型校正系数，反映泥沙输送的比例，表示泥沙有多少能进入河流。通常赋为0.5。

最大L值

LS因子中斜率长度参数（L）的最大允许值，默认参数值为122。

运行

得到的记过我们需要计算一下土壤保持量SD
公式如下：

所以：

• 1、将rkls和usle数据加载到ArcGIS Pro中。
• 2、加载好数据后，在【地理处理】处打开【栅格计算器】，在对话框中输入： “rkls.tif” - “usle.tif” 点击【运行】即可得到土壤保持量数据。

进行重分类

生态适宜性评价

至此我们有了四个模块的结果，四个模块分别是

根据前面所说，用arcgis的镶嵌至新栅格，选择最大值镶嵌，即可得到最终的生态适宜性

然后掩膜一下

处理后可能会有一个值为0的，这个值我们可以在属性表里删了，或者重分类将0和1分为一列1即可
这样就得到了三个不同适应性范围的结果了