计算机视觉——凸包计算

现在有一大堆点,然后你要找出一个可以围住这些点且面积最小的凸多边形,这个凸多边形称为凸包。

显而易见,如果要面积最小,那凸包的顶点势必得是这一大堆点的几个点,你也可以想成是用一条橡皮筋把这些点圈起来。

先把各个点按坐标从小到大排序,坐标相同则再按坐标从小到大排序,排序之后的点顺序会是由左至右、由下至上,这样一来我们就可以按这个顺序遍历这些点,这种往固定方向扫描的方式,称为扫描线。

先讨论一件事情:有一个凸多边形,它的顶点已经按逆时针顺序排好了,依次是 p 1 , p 2 , . . . , p n p_1, p_2, ..., p_n p1,p2,...,pn ,那么 p i p_i pi 和 p j p_j pj 与 p k p_k pk 的关系会是什么(令 i < j < k i < j < k i<j<k既然是凸多边形,那么它的边应该是要一直往同个方向转弯的,而如果将边逆时针排序,这个边的斜率也应该是一直往逆时针方向转弯,显然点也会是这样,因此:

再来,我们把凸包分成上下两部分:上凸包和下凸包,以极左和极右点分割,如果极左点或极右点有两个(最多只会有两个),上面那个属于上凸包,下面那个属于下凸包,否则极左点必属于下凸包,极右点必属于上凸包。

显然,上或下凸包中不会有 x 坐标相同的顶点,因此在上或下凸包中,每个顶点都能分别出左右的关系的,并且你会发现,如果把点按逆时针排序,在下凸包中的点也是由左而右排序的、在上凸包的点也是由右而左排序的。

这样一来左右关系就有用了,先用由左而右的扫描线把下凸包做出来,再用由右而左的扫描线把上凸包做出来,就可以得到整个凸包。

整个流程可以用一个栈来实现,在处理一个点的时候,我们尝试把它加进凸包里,此时这个点是 ( p ) ,栈顶的点是 ( top ) ,栈顶再往下一个点是 ( second ) ,把这些点代入刚刚的式子,符合条件或者栈大小小于 2 时就停止,否则就把栈顶 pop,然后继续重复,结束后就把目前处理中的点放入栈顶。

在做下凸包的时候,先从最左边且最下面的点开始做上述动作,做到最后,栈顶的点应该是最右边且最上面的点,把它 pop 掉,因为它应该属于上凸包;做上凸包的时候,从最右边且最上面的点开始做,最后栈顶会是最左且最下的点,把它 pop 掉后,这两个接起来就是完整的凸包。

因为要用到栈顶往下一个点,所以栈用 vector 来实现。

cpp 复制代码
#include <bits/stdc++.h>

#define mp(a, b) make_pair(a, b)
#define pb(a) push_back(a)
#define F first
#define S second

using namespace std;

template<typename T>
pair<T, T> operator-(pair<T, T> a, pair<T, T> b){
    return mp(a.F - b.F, a.S - b.S);
}

template<typename T>
T cross(pair<T, T> a, pair<T, T> b){
    return a.F * b.S - a.S * b.F;
}

template<typename T>
vector<pair<T, T>> getConvexHull(vector<pair<T, T>>& pnts){
    int n = pnts.size();
    sort(pnts.begin(), pnts.end());

    vector<pair<T, T>> hull;

    for(int i = 0; i < 2; i++){
        int t = hull.size();
        for(auto& pnt : pnts){
            while(hull.size() - t >= 2 && cross(hull.back() - hull[hull.size() - 2], pnt - hull[hull.size() - 2]) <= 0)
                hull.pop_back();
            hull.pb(pnt);
        }
        hull.pop_back();
        reverse(pnts.begin(), pnts.end());
    }

    return hull;
}

旋转卡尺

用旋转两条平行线、夹住一堆点,看在线上的点是哪些,就叫旋转卡尺。

旋转线、夹点感觉很麻烦,是不是要用什么角度的东西啊?其实不用,先来分析一下问题,用两条平行线夹一堆点,那么平行线只会碰到凸包上的点而已,所以不在凸包上的点都可以先忽略:

过一个点的直线有很多条,但是过一个线段的直线只有一条,所以先枚举线段,再去找和它平行的直线应该会夹到哪个点,这样问题就简单多了。要找平行线会碰到哪个点,显然离线段最远的点就是了。

不过算距离是另一个问题,听起来也很麻烦,但其实很简单。一个点距离一条直线的距离,等同于过该点在直线上作垂线段的长,而一开始选定的线段作为底、垂线段长作为高,那么就可以得一个平行四边形面积了,且底的长是固定的,只要枚举最远点,就等同于枚举高,而得出面积最大的,就是我们要求的最远点。

上图中,选定两个红色点所连成的线段为底,然后枚举各个顶点取高,得出蓝色垂线是最长的,因此蓝色点就是距离红色线段最远的点。

这就是旋转卡尺的基础应用------最远点对,找到距离每一线段最远的点,再取该点与线段两端点的距离取最大值,这样就可以得出所有点中最远的点对为何。

硬要这么做的方式,时间复杂度是 O ( n 3 ) O(n^3) O(n3), n n n 是凸包上点的数量(不计盖凸包的复杂度),枚举线段是 O ( n 2 ) O(n^2) O(n2) ,再枚举一个点要再乘上 n n n 。

这不够快,我们需要更有效率的方式。

仔细观察一下,点和线段的距离有一个规律------先渐大,到一个最大值,再渐小:

会发现它会呈现一个单峰函数,也就是一个先递增、再递减的函数,这样我们就可以三分搜找到最高点了,这样三分搜一次的复杂度是 O ( n ) O(n) O(n),再乘上点的数量,就是 O ( n 2 ) O(n^2) O(n2)。

这样子还是不够快,前面提到旋转卡尺是「旋转两条平行线」,刚才的动作都是旋转其中一条,再去搜寻另一条,那我们可不可以在旋转其中一条的同时,把另一条一起旋转?答案是:可以。

(以下的转都是指往同一个方向转)

先找到距离第一条边最远的点,过前者的线称为第一条平行线,过后者的称为第二条平行线,接下来我们转动第一条平行线,也就是把它转到第二条线段上,而第二条平行线不要动,会发现,第一条平行线离第二条平行线那个点近了一些,接着再转第二条平行线,也就是把它转到下一个点上,那么距离会变远。

也就是,可以在不重新来过的情況下,找到单峰函数的最高点,会发现这样就是把两条平行线绕一圈,因此这样的复杂度是 O ( n ) O(n) O(n) 。

原文地址:https://cp.wiwiho.me/convex-hull/

相关推荐
martian6651 小时前
【人工智能数学基础篇】——深入详解多变量微积分:在机器学习模型中优化损失函数时应用
人工智能·机器学习·微积分·数学基础
人机与认知实验室2 小时前
人、机、环境中各有其神经网络系统
人工智能·深度学习·神经网络·机器学习
黑色叉腰丶大魔王2 小时前
基于 MATLAB 的图像增强技术分享
图像处理·人工智能·计算机视觉
迅易科技5 小时前
借助腾讯云质检平台的新范式,做工业制造企业质检的“AI慧眼”
人工智能·视觉检测·制造
古希腊掌管学习的神6 小时前
[机器学习]XGBoost(3)——确定树的结构
人工智能·机器学习
ZHOU_WUYI6 小时前
4.metagpt中的软件公司智能体 (ProjectManager 角色)
人工智能·metagpt
靴子学长7 小时前
基于字节大模型的论文翻译(含免费源码)
人工智能·深度学习·nlp
AI_NEW_COME8 小时前
知识库管理系统可扩展性深度测评
人工智能
海棠AI实验室8 小时前
AI的进阶之路:从机器学习到深度学习的演变(一)
人工智能·深度学习·机器学习
hunteritself8 小时前
AI Weekly『12月16-22日』:OpenAI公布o3,谷歌发布首个推理模型,GitHub Copilot免费版上线!
人工智能·gpt·chatgpt·github·openai·copilot