在生境图中区分走廊与碎片


12

我想知道以前是否有人对这种问题感到沮丧:

我想基于距离为3像素(对我感兴趣的物种具有生物学意义)的同一类别(森林)的相邻像素来描绘片段。

我担心的是,例如,请参见所附图像,有时这些片段实际上是走廊,并且由于走廊和实际片段的接近性,它们经常被分组为同一片段。

在此处输入图片说明

我想知道是否有一种方法可以根据形状,周围像素的数量等将走廊与片段区分开来?

例如,在下面的框中,可能的走廊用红色框表示,而片段用绿色框表示。

我可以访问QGIS和R。

在此处输入图片说明


您可以添加一个图形,其中包含您要从示例中提取的走廊类型吗?
radouxju

Answers:


4

在开始任何分析之前,我强烈建议对您的数据应用过滤器,以清除“盐和胡椒”效应。任何算法都将与数据的当前结构模式作斗争。简单的多数表决可能会产生不希望的结果。一种更可靠的方法是应用筛分方法,其中可以指定最小映射单位。这可以通过GDAL中的gdal_sieve.py函数,raster > analysis > sieveQGIS中的函数或ArcGIS渐变度量工具箱中的sieve函数来完成。

这看起来可以用数学形态学运算符解决(例如,从图像中提取道路)。我想像一下,用一个膨胀运算符再加上一个闭合运算符可以弄清楚走廊。然后,您可以应用“打开”运算符,以删除道路,并对栅格进行差异化以将识别出的道路拉为单独的对象。这些类型的图像分解功能在MSPA和GUIDOS软件中有些自动化,但是再次受到数据不连续性的明显影响。

有一个用于MSPA的QGIS插件以及GRASS中的可用功能(可通过QGIS GUI使用)。MSPA和GUIDOS的问题之一是图像大小受到限制。不幸的是,在ESRI软件中,形态运算符仅在ArcScan扩展中可用。通过一些挖掘,您将找到其他软件选项以及通过具有自定义内核矩阵的栅格代数定义形态算子的方法。

另一种方法是边缘检测过滤方法,例如Sobal核算子。在ArcGIS Gradient Metrics Toolbox中以及spaceEco R程序包中都有一个sobal函数。R实现的优点是您可以返回运算符的梯度函数,而ArcGIS实现仅返回一阶函数(可能就是您所需要的)。我相信Orfeo工具箱(作为QGIS附件提供)在EdgeExtraction功能中具有sobal选项。


也许使用“盐和胡椒”过滤器可能会去除目标物种连接片段的重要小栖息地(尽管有时质量欠佳)。我会根据物种的扩散特征和数据的空间分辨率来决定使用这种过滤器。鸟类(通常具有较高的分散能力)可以使用那些小的斑块,而两栖动物(具有相对较低的分散能力)则可以不使用。只是一个想法..
卡莫

这始终是一个折衷,但是最小映射单元应该考虑这个给定的假设。但是,您还必须牢记,平衡是该模型表现出一定程度的不确定性,并且您将结果视为绝对值。老实说,为了支持您的主张,应将景观视为概率梯度而不是二项式过程。这种过滤是代表更多功能区域的一种长期接受的做法。除非模型在估算中包括考虑空间结构的术语,否则空间不确定性是现实。
Jeffrey Evans

5

这不是一个完整的解决方案,但请查看以下用于连通性分析的工具(第一个工具很适合您搜索的内容):

还可以考虑训练机器学习算法来对您的实例进行分类(走廊与片段)。您可以在补丁程序级别(例如补丁程序大小,周长面积比,圆比)和基于距离的要素(例如到碎片的距离)上提供空间属性进行分类。要计算分类所需的补丁程序级功能,您可以尝试FragStats(https://www.umass.edu/landeco/research/fragstats/fragstats.html)。

您还可以考虑使用更简单的“基于专家的规则系统”对每个实例进行分类。例如,走廊将比栖息地碎片等具有更高的周长比。

有关连接性分析的更多有趣内容:http : //conservationcorridor.org/corridor-toolbox/programs-and-tools/

但是,某些走廊没有“全像素连接”这一事实将是您首先需要解决的问题。我认为您必须定义某种基于距离的阈值标准,才能确定某个像素是否属于走廊的一部分。


0

看来这取决于片段的形状。如果宽度大于高度(或相反高度)的2(或3)倍,您可以称其为走廊吗?

您是否已经到达碎片的划界了?

By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.