MATLAB Delaunay算法提取离散点边界的方法

2020-01-06 19:59:48刘景俊

最近在项目进行中遇到要提取离散点边界的问题,像我这样的对于matlab不是特别熟练的朋友一开始肯定摸不着头脑,到底选用哪种算法可以有效地提取到所有已知点的轮廓线呢。本人经过大量的文献搜索及代码实验找到了几个效果比较好的轮廓提取代码,在这里做个总结,并且希望能够对遇到同样问题的朋友有所启发。

关于离散点边界提取的三种方法:

1.Convhull 离散点集获得边界

2.Alpha Shape算法检测边缘点

3.Delaunay 三角剖分算法

前两种方法在之前的博客中已经做了总结这里就不展开了,现在主要介绍第三种算法。

该算法的总体思路如下:

1、利用 delaunay 函数,对所有数据点进行 Delaunay 三角剖分处理,delaunay 函数的返回值是一个 N * 3 的矩阵,其中 N 为剖分出的三角形个数,3 为每个三角形的三个端点的序号。

2、根据 triangles 矩阵,提取出所有 delaunay 三角剖分时所连接的边,依次扫描 triangles 矩阵的每一行,将 delaunay 三角剖分时所连接的边添加到一个新的矩阵中,最后构成一个 M * 2 的矩阵,其中 M 是一共所连接的边的条数。

3、显然,最小凸多边形上的边应该仅在以上矩阵中出现一次,因此,将以上矩阵中那些出现次数超过一次的边全部去掉,最后保留的便是最小凸多边形的边。

4、根据最小凸多边形的边,很容易得到构成最小凸多边形的结点的顺序,从而解决问题。