基于InVEST-MSPA-MCR模型的生态网络构建（二）

 ~~~继上文InVEST模型的识别结果，我们接着进行生态网络的构建~~~

        由于是第一次写博客，本意也只是记录并分享一下自身所学习到的知识，因此，如果在阅读时，发现有错误的地方请见谅，欢迎评论区指出！

ps：网络构建的构成内容比较多，使用到的软件比较多，因为是已经完成的项目，所以有些步骤可能省略了，主要是分享一个大致的流程~~因此，若需要深入学习相关软件，也可自行查询相关文章！

（2）MSPA识别潜在生态源地

        MSPA利用二值图像的图斑特点，将栅格图像分为了7个区（图），一般直接选取核心区作为生态源地。对于MSPA的计算，使用guidostoolbox完成，在进行MSPA分析前的，需要对数据进行一定的处理，这里我参考的是b站up主分享的MSPA使用经验：【ArcGIS应用系列-MSPA模型简要演示】https://www.bilibili.com/video/BV1Mp4y1v7Ha?vd_source=bcde76257b9d7c6771fad63c6a41f032。

 MSPA分区示意图：    

MSPA分区的生态学意义（主要是提取其中的核心区）：

        

随后，将INVEST模型获得的生境质量结果进行重分类，将生境质量0.6（这个阈值设置很重要，我这里参考了相关文章）以上的区域作为生态源地待选择区。背景取1，前景取2，nodata取0。

下图为重分类结果：

        包括一些参数的设置，如下，具体含义参考使用手册：

        

        计算结果就是这样的啦：

       

        随后使用一下arcgis的提取功能，将生态核心区提取。因为生态核心区太过琐碎，需要进行一定的筛选，筛选方法一般通过计算源地斑块之间的dpc指数，选取其中大于1（经验值，具体为什么也不知道）的斑块作为源地，dpc计算使用conefor 26完成（一个联合arcgis的插件），dpc指数是基于PC指数结果的，而PC指数的计算是需要基于一个距离阈值的，该距离阈值一般通过不同距离的下的NC指数、IIC指数、以及PC指数的变化来确定，这里选取2000-6500m之间各距离阈值，进行计算NC、IIC、PC指数，所使用软件为conefor 26，网址：Download Conefor | Conefor。

在这个界面找到适合你arcgis的版本

        首先在arcgis中使用插件计算斑块之间的距离连接线，生成相应的数据，随后在conefor26中选取要计算的指数，将各距离阈值下的NC、IIC、PC指数绘制成表格与图画，基于指数的变化趋势，确定5000m为合适的距离阈值。

插件长这样：

       

 工具界面如下：

  计算结果：     

主要利用的是两个txt文件，获得该两个txt文件后，打开conefor26独立软件，界面如下：

计算结果：

进行统计一下，就可以获得NC、IIC、PC指数随距离的变化规律啦~~：

        

       经过一顿胡乱分析，最终确定了5000m为阈值，计算核心区斑块斑块的dpc指数（用Fragstats 4.3软件计算），并利用arcgis通过表连接将斑块与dpc表相连，再筛选dpc>1的斑块作为生态源地（36个），经过连接COUNT属性表，发现在5000m为阈值下，所筛选的生态源地，其连接性路径数量总数为（634），占总的连接路径（2386）的26.5%。至此，重要生态源地识别结束。（这一段涉及的操作比较多，我后面有时间再补一下。。）

接下来将基于生态源地的分布，进行生态廊道的识别，主要还是MCR（最小累积模型的原理）

----》未完待续

9月12日-----------------------------------------------------------------------

博主近段时间有要事要忙，恐不能更新，读者可先另寻高文，以滋求识之思！