Invest模型 ——生境质量计算

        上一期推文中，介绍了Invest模型的安装以及基本使用方法。本期向大家介绍使用Invest模型的生境质量计算。

        生境质量指生存资源的可获得性，高的生境质量可以看作其结构、功能的完整并在历史动态范围内。生境质量取决于一个生境对人类土地利用和这些土地利用强度的可接近性。

        生境质量模型如下图所示。其中红色边框圈出区域为必须输入数据。

工作空间：为自定义文件路径，这里就不赘述了，但是其中路径应该均为英文路径。

土地利用情况：土地利用数据要在Arcgis进行预处理。将nodata数据重分类为0，操作流程如下：

【system Toolboxes】-【Spatial Analyst Tools】-【Reclass】-【Reclassify】

重新分类后结果如下所示：

 胁迫因子：此处输入数据为CSV文件数据。文件格式如下图所示，其中【CUR_PATH】列的数据，在Arcgis中根据土地利用数据将水体、人造地表、裸地、水田、工矿用地等胁迫源进行提取；【MAX_DIST】为最大影响距离，此数值可以通过阅读与自己相同研究区的文献进行查找；【WEIGHT】影响权重，同样数值的确定查阅相关论文或者视自己的研究，采用主成分分析、专家评价等进行赋值；【THREAT】是胁迫源，这一列的名字要与采用Arcgis进行胁迫源提取的栅格图层名称去掉_c情况下一致；【DECAY】影响类型，此处可以为线性影响、指数影响等，更多的类型可以去查阅使用手册。

胁迫源提取：首先看好每种胁迫源重分类后的代码是什么，然后根据代码进行提取。【system Toolboxes】-【Spatial Analyst Tools】-【Extract】-【Extract by Attributes】

        如进行水体的提取，提取后需要将其重分类，有水体的地方分类为1，无水体的地方分类为0。将所有胁迫源均像水体一样进行提取。

注意：胁迫源命名应该为xxx_c.tif（命名务必遵守规则，主要是为了模型的数据识别）

水体最终分类结果如下图所示：

敏感性：此处输入文件同样为CSV格式文件。

    其中【LULC】列的编号一定要和重分类后的编号保持一致【NAME】列的名称可以随便命名；【HABITAT】列为生境适宜性，同样查找文献来确定，人造地表、裸地等可能会低一点，临界水体等要高一点。最后五列的数值代表每一类的土地受到威胁源影响的权重是多少。

半变异系数：根据自定义值0.5进行，详细可查找使用手册。

将上述5个必须输入参数输入至模型，运行后，便可以得到结果。最终的结果数据为：

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

---

获取上述数据，请关注公众号GIS攻略，后台回复0509供大家练习。