计算机视觉——凸包计算

news2024/11/16 9:03:54

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

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

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

先讨论一件事情:有一个凸多边形,它的顶点已经按逆时针顺序排好了,依次是 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 来实现。

#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;
}

旋转卡尺

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

旋转卡尺1

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

旋转卡尺2

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

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

旋转卡尺3

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

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

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

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

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

旋转卡尺4

会发现它会呈现一个单峰函数,也就是一个先递增、再递减的函数,这样我们就可以三分搜找到最高点了,这样三分搜一次的复杂度是 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/

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1992028.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

什么是云原生?(一)

1. 前言 停下手头的工作&#xff0c;让你的同事定义“云原生”一词。你很可能会得到几个不同的答案。 1.1 让我们从一个简单的定义开始&#xff1a; 云原生架构和技术是一种设计、构建和操作在云中构建并充分利用云计算模型的工作负载的方法。 1.2 云原生计算基金会给出了官方…

easyexcel使用教程--导入导出简单案例

java中操作excel的库是poi,但是poi有很多问题&#xff0c;使用复杂&#xff0c;内存占用大等&#xff0c;阿里开源了一个操作excel的库叫easyexcel,它基于poi做了很多优化&#xff0c;平时读写excel建议使用这个库 文档地址: 关于Easyexcel | Easy Excel 官网 写入excel 在…

操作无法完成错误0x000006ba?教你几种全面解析与解决方案指南

错误代码0x000006ba通常与 Windows 操作系统中的远程过程调用&#xff08;RPC&#xff09;服务有关。当你在尝试执行某些操作&#xff0c;如连接到网络共享、运行某些网络服务或使用依赖于 RPC 的应用程序时&#xff0c;可能会遇到这个错误。接下来就和大家聊聊操作无法完成错误…

“揭秘CentosChina爬虫项目:掌握Scrapy框架的必备技巧与数据库设计“

Centoschina 项目要求 爬取centoschina_cn的所有问题&#xff0c;包括文章标题和内容 数据库表设计 库表设计&#xff1a; 数据展示&#xff1a; 项目亮点 低耦合&#xff0c;高内聚。 爬虫专有settings custom_settings custom_settings_for_centoschina_cncustom_settin…

鸿蒙(API 12 Beta2版)媒体开发【使用OHAudio开发音频播放功能(C/C++)】

OHAudio是系统在API version 10中引入的一套C API&#xff0c;此API在设计上实现归一&#xff0c;同时支持普通音频通路和低时延通路。仅支持PCM格式&#xff0c;适用于依赖Native层实现音频输出功能的场景。 使用入门 开发者要使用OHAudio提供的播放或者录制能力&#xff0c…

根据 IP 地址配置子网示例(下挂 hub 接不同 vlan 终端)

我们一般根据端口配置子网比较简单&#xff0c;但是如果换接口&#xff0c;就又要到交换机上重新配置端口所属 vlan 了&#xff0c;紧急情况下&#xff0c;还是比较耽误时间的。但如果根据IP地址配置 vlan&#xff0c;则可以插在交换机上任意端口&#xff0c;排障时比较节省时间…

部分PC制造商不会帮助英特尔第13/14代酷睿延保 对用户来说可能是个问题

早前英特尔已经宣布大部分第 13/14 代酷睿桌面处理器将可以获得延长保修服务&#xff0c;即后续使用出现非人为损坏的问题后&#xff0c;用户都可以联系英特尔换一颗新处理器。然而英特尔承担责任不代表对用户来说就高枕无忧&#xff0c;例如只有盒装版处理器才能直接联系英特尔…

MySQL数据库分区

文章目录 1. 分区的基本概念2. 分区的类型3. 在 MySQL 中实现分区3.1 创建分区表3.2 插入数据3.3 查询数据3.4 修改分区3.5 查看分区信息 4. 具体示例4.1 范围分区 (Range Partitioning)4.2 列表分区 (List Partitioning)4.3 哈希分区 (Hash Partitioning)4.4 键分区 (Key Part…

Win32注册表操作

注册表的概念 注册表是一个存储计算机配置信息的数据库&#xff0c;用于存储计算机上的硬件、安装的软件、系统设置以及用户账户配置等重要信息。对注册表的编辑不当可能会影响计算机的正常运行。应用程序可以调用API函数来对注册表进行增、删等操作。 注册表结构 运行Regedi…

【Redis 进阶】Redis 典型应用 —— 分布式锁

一、什么是分布式锁 在一个分布式的系统中&#xff0c;也会涉及到多个节点访问同一个公共资源的情况&#xff0c;此时就需要通过锁来做互斥控制, 避免出现类似于 “线程安全” 的问题。 而 Java 的 synchronized 或者 C 的 std::mutex&#xff0c;这样的锁都是只能在当前进程…

Embedding技术之Word Embedding

Word Embedding是什么&#xff1f; Word Embedding——词嵌入是将单词映射为数值向量&#xff0c;以捕捉单词间的语义和句法关系&#xff0c;为自然语言处理任务提供有效的特征表示。——自然语言处理——Word2Vec、GloVe、FastText。 Word Embedding属于静态向量&#xff0c…

深度学习入门(四):激活函数与LSTM

激活函数 激活函数在神经网络中扮演着至关重要的角色。它们的主要功能是引入非线性因素&#xff0c;这使得神经网络能够学习和近似任何非线性函数&#xff0c;从而处理复杂的数据模式和决策边界。如果没有激活函数&#xff0c;即使网络拥有很多层&#xff0c;其表达能力仍然只…

青岛国真携手图扑软件共建青岛西海岸区一网统管平台

为深入贯彻关于垃圾分类的重要指示精神&#xff0c;积极响应住房和城乡建设部开展的全国城市生活垃圾分类宣传周活动&#xff0c;5 月 25 - 27 日&#xff0c;以“发挥行业协会新动能 助力垃圾分类新时尚”为主题的“ 2023 垃圾分类高峰论坛”在青岛西海岸成功举办。 青岛国真智…

使用影子凭证进行域权限维持

本文来源无问社区&#xff0c;更对实战内容&#xff0c;渗透思路可前往查看http://wwlib.cn/index.php/artread/artid/15293.html Microsoft 推出了 Windows Hello 企业版 &#xff08;WHfB&#xff09;&#xff0c;以使用基于密钥的信任模型替换传统的基于密码的身份验证。此…

【Material-UI】按钮组:尺寸与颜色详解

文章目录 一、按钮组概述1. 组件介绍2. 基本用法 二、按钮组的尺寸&#xff08;Sizes&#xff09;1. 小尺寸&#xff08;Small&#xff09;2. 中等尺寸&#xff08;Medium&#xff09;3. 大尺寸&#xff08;Large&#xff09; 三、按钮组的颜色&#xff08;Colors&#xff09;1…

gitea docker 快捷安装部署

前言 在前一篇博文&#xff08;什么是 Gitea&#xff1f;&#xff09;中&#xff0c;我们详细介绍了gitea的功能特性&#xff0c;以及其与其它git服务器之间的特性多维度对比。 在本文中&#xff0c;我们将详细介绍gitea的快捷安装部署&#xff0c;docker方式&#xff01; 1…

qt客户端与服务端通信

服务器要一处于监听状态 客户端主动连接服务器 服务器的ui界面 客户端的ui界面 很简陋 服务端listen 按钮的槽函数 QHostAddress::Any,port 监听任何端口 void MainWindow::on_listen_clicked() {if(ui->lineEdit->text().isEmpty()){return ;}int port ui->li…

【第九节】python中xml解析和json编解码

目录 一、Python XML 解析 1.1 什么是XML 1.2 Python 对 XML 的解析方法 1.3 SAX解析xml 1.4 xml.dom解析xml 1.6 ElementTree解析XML 二、Python编解码json 2.1 什么是json 2.2 使用json 库 2.3 使用第三方库Demjson 一、Python XML 解析 1.1 什么是XML XML&#x…

Python新手错误集锦(PyCharm)

# 自学Python&#xff0c;用Pycharm作环境。我这个手新到这时我学习的第一个编程软件&#xff0c;且本人专业是化学&#xff0c;以前对电脑最高级的使用是玩扫雷游戏。所以这里集合的错误都是小透明错误&#xff0c;大部分人请绕道。不断更新中...... 缩进错误 记住“indent”…

Jmeter之BeanShell使用(全网最详细的介绍)-第九天

一.什么是BeanShell Java写成的⼩型、免费的Java源代码解释器可以执⾏标准Java语句和表达式,完全符合java语法的java脚本语⾔&#xff08;需要会javase语⾔&#xff09;包括⼀些脚本命令&#xff0c;有⾃⼰的⼀些语法和⽅法&#xff0c;是⼀种松散类型的脚本语⾔&#xff08;这…