初学计算几何（四）——初识凸包

最新推荐文章于 2022-07-01 19:49:22 发布

weixin_30500289

最新推荐文章于 2022-07-01 19:49:22 发布

阅读量94

点赞数

原文链接：http://www.cnblogs.com/chenxiaoran666/p/ConvexHull.html

版权

前言

学习过了多边形的一些相关内容，总算可以开始学习凸包了。

这篇博客主要介绍如何在给出的点集中求出凸包。

关于凸包面积可以参考初学计算几何（三）——多边形的简单操作中的多边形面积的求法。

排序

求凸包的第一步便是将点集按照这个方法进行排序：

inline bool operator < (Vector A,Vector B) {return fabs(A.x-B.x)>eps?A.x<B.x:A.y<B.y;}

这应该还是比较好理解的吧，就是以\(x\)坐标为第一关键字，以\(y\)坐标为第二关键字进行排序。

接下来的步骤

显然，排完序后得到的第一个点和最后一个点肯定在凸包内。

每次要加入一个新的点，如果已加入点数大于\(1\)，我们将当前点\(p\)与最后加入的点\(S_n\)比较，如果\(\vec{S_{n-1},S_n}\)不在\(\vec{S_{n-1},p}\)左边，我们就可以将\(S_n\)弹出。

不断重复该过程，直到无法继续弹出了，我们再将点\(p\)加入凸包中。

然后倒着执行一遍类似的操作即可。

具体实现见代码：

inline ConvexHull GetConvexHull(Polygon S)
{
    register int i,t;register ConvexHull res;
    for(sort(S.p+1,S.p+S.n+1),i=1;i<=S.n;++i)//先排一遍序
    {
        while(res.n>1&&dcmp(Cro(res.p[res.n]-res.p[res.n-1],S.p[i]-res.p[res.n-1]))<=0) --res.n;//将不满足条件的点弹出
        res.p[++res.n]=S.p[i];//将当前点加入凸包中
    }
    for(t=res.n,i=S.n-1;i;--i)//倒着执行一遍类似的操作
    {
        while(res.n>t&&dcmp(Cro(res.p[res.n]-res.p[res.n-1],S.p[i]-res.p[res.n-1]))<=0) --res.n;
        res.p[++res.n]=S.p[i];
    }
    return res;
}