给定n个点或一个凸边形,求其最小外接矩形,可视化

这里写目录标题

  • 原理
  • 代码

原理

  1. 求n个点的最小外接矩形问题可以等价为先求这n个点的凸包,再求这个凸包的最小外接矩形。

    其中求凸包可以使用Graham-Scan算法
    需要注意的是,
    因为Graham-Scan算法要求我们从先找到凸包上的一个点,所以我们可以先对这些点进行排序,根据x/y坐标都可以。排序后的第一个点和最后一点一定在凸包上。
    同时,假设我们按照x排序,并且想要得到顺时针旋转的凸包边界,那么我们只能找到一半的边界,这是因为我们从左到右遍历这些点的时候只有上面的点才满足顺时针,下面的不可能满足顺时针(这里为了节省计算可以只对y坐标大于开始和结束点中最小y坐标的点进行计算)。
    为了找到全部的点,我们需要正向判断完后再反向相同方法判断一遍。
    这里为了判断是顺时针还是逆时针,我们使用的是两个方向叉乘的最后一项,从前一个方向转到后一个方向时:如果叉乘最后一项小于0,说明是顺时针旋转的(右手系)。
    在这里插入图片描述

  2. 给定一个凸多边形,求最小外接矩形的原理可以参考这篇博客https://blog.csdn.net/wang_heng199/article/details/74477738
    简单而言就是对凸包上的每条边构造矩形并比较面积。
    我们这里的做法是遍历每条边,对于每条边,我们根据它与水平方向的夹角计算其旋转矩阵,并将这些凸包上的所有点根据此旋转矩阵的逆矩阵/转置矩阵(旋转矩阵是正交矩阵)旋转至水平方向,此时我们可以很容易得到这些点水平方向和竖直方向的边界点,从而就有了一个外接矩形。

代码

导入包并且初始化随机点:

import numpy as np
import matplotlib.pyplot as plt
n = 100
points = np.random.randn(n, 2)

计算凸包并可视化

sorted_points = points[np.argsort(points[:,0])]def cal_rotate_orientation(cur,stack_end,stack_end_before):'''计算p1->p2的向量旋转到p2->p3的向量是顺时针还是逆时针'''p1,p2,p3 = sorted_points[stack_end_before], sorted_points[stack_end], sorted_points[cur]x2,y2 = p3 - p2x1,y1 = p2 - p1return x1*y2 - x2*y1    #叉乘最后一项Z轴方向,<0:顺时针,>0:逆时针def cal_convex_hull(index_list):result = list(index_list[:2])index_list = index_list[2:]for i in index_list:if (cal_rotate_orientation(i,result[-1],result[-2]) < 0):  #叉乘最后一项 ,符合顺时针# print(i,result[-1],result[-2])result.append(i)else:if len(result) == 2: result[-1] = i   continueresult.pop()while(cal_rotate_orientation(i,result[-1],result[-2]) > 0):   #不是顺时针,需要反复剔除result.pop()if len(result) < 2:breakresult.append(i)return result#n-1是最后一个点,因为它肯定在第一遍的结果里,所以两个结果合并时去除第一个
result = cal_convex_hull(range(0,n,1))+cal_convex_hull(range(n-1,-1,-1))[1:]
print(result)
#ploy是将结果按照顺序两两组合,形成边界
ploy = list(map(lambda i:result[i:i+2],range(0,len(result)-1)))plt.scatter(sorted_points[:, 0], sorted_points[:, 1])
plt.plot(sorted_points[result,0],sorted_points[result,1],color = 'red')
plt.show()

在这里插入图片描述

计算最小外接矩形
这里为了可视化,我们会将矩形转回来。

def get_rotate(cos,sin):#我们已知在原坐标系下的方向,要将其转化到目标系,用逆/转置return np.array([[cos,-sin],[sin,cos]])   def get_rotate2(theta):#我们已知在原坐标系下的方向,要将其转化到目标系,用逆/转置return np.array([[np.cos(theta),-np.sin(theta)],[np.sin(theta),np.cos(theta)]])row = 3
col = (len(ploy)+1)//3+1
fig, axs = plt.subplots(row, col, figsize=(15, 15))min_rect = [[],0,np.inf] #rect,angle,area
#假设是-90度 -> 90度
for idx, (i,j) in enumerate(ploy):x1,y1 = sorted_points[i]x2,y2 = sorted_points[j]tan_theta = (y2-y1)/(x2-x1)cos_theta = np.sqrt(1/(1+np.power(tan_theta,2)))sin_theta = cos_theta*tan_thetarotated_points = sorted_points@get_rotate(cos_theta,sin_theta)min_x,min_y = np.min(rotated_points,axis=0)max_x,max_y = np.max(rotated_points,axis=0)rect_xy = np.array([[min_x,min_y],[max_x,min_y],[max_x,max_y],[min_x,max_y]])rect_area = (max_x-min_x)*(max_y-min_y)if rect_area < min_rect[2]:min_rect = [rect_xy,np.arcsin(-sin_theta),rect_area]  #这个角度是取反,因为要转回去rect_xy = rect_xy@get_rotate(cos_theta,-sin_theta)  #转一下axs[idx//col,idx%col].scatter(points[:, 0], points[:, 1])plot_idx = [i for i in range(0,4)]+[0]axs[idx//col,idx%col].plot(rect_xy[plot_idx,0],rect_xy[plot_idx,1],color = 'orange')axs[idx//col,idx%col].set_xlim(-4,4)axs[idx//col,idx%col].set_ylim(-4,4)axs[idx//col,idx%col].set_aspect('equal')rect_xy = min_rect[0]@get_rotate2(min_rect[1])axs[(idx+1)//col,(idx+1)%col].scatter(points[:, 0], points[:, 1])
plot_idx = [i for i in range(0,4)]+[0]
axs[(idx+1)//col,(idx+1)%col].plot(rect_xy[plot_idx,0],rect_xy[plot_idx,1],color = 'red')
axs[(idx+1)//col,(idx+1)%col].set_xlim(-4,4)
axs[(idx+1)//col,(idx+1)%col].set_ylim(-4,4)
axs[(idx+1)//col,(idx+1)%col].set_aspect('equal')
plt.show()

在这里插入图片描述

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

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【加载自定义控制器 Objective-C语言】

一、接下来要做的操作,就比较关键了 1.前面,我们在从UI基础,然后到UI进阶,第一天,然后到目前为止, 所有的应用程序,新建的项目,启动以后,加载的控制器,是不是都是Main.storyboard里面带箭头儿的那个控制器, 然后呢,你也可以通过新建一个storyboard,然后呢,给它…

默默消失的6大APP,不少中年人估计都用过?

随着科技的发展&#xff0c;这些曾经风靡一时的APP逐渐消失在人们的视线中。然而&#xff0c;它们的消失并不代表它们的价值被完全抹去。相反&#xff0c;它们所代表的一段历史和人们的回忆将永远留存在人们心中。 第一个千千静听是我们青少年时代的伴侣&#xff0c;它陪伴我们…

Wnmp服务安装并结合内网穿透实现公网远程访问——“cpolar内网穿透”

文章目录 前言1.Wnmp下载安装2.Wnmp设置3.安装cpolar内网穿透3.1 注册账号3.2 下载cpolar客户端3.3 登录cpolar web ui管理界面3.4 创建公网地址 4.固定公网地址访问 前言 WNMP是Windows系统下的绿色NginxMysqlPHP环境集成套件包&#xff0c;安装完成后即可得到一个Nginx MyS…

MSF暴力破解SID和检测Oracle漏洞

暴力破解SID 当我们发现 Oracle 数据库的 1521 端口时,我们可能考虑使用爆破 SID(System Identifier)来进行进一步的探测和认证。在 Oracle 中,SID 是一个数据库的唯一标识符。当用户希望远程连接 Oracle 数据库时,需要了解以下几个要素:SID、用户名、密码以及服务器的 I…

【LeetCode:318. 最大单词长度乘积 | 模拟 位运算】

&#x1f680; 算法题 &#x1f680; &#x1f332; 算法刷题专栏 | 面试必备算法 | 面试高频算法 &#x1f340; &#x1f332; 越难的东西,越要努力坚持&#xff0c;因为它具有很高的价值&#xff0c;算法就是这样✨ &#x1f332; 作者简介&#xff1a;硕风和炜&#xff0c;…

【广州华锐互动】VR历史古城复原:沉浸式体验古代建筑,感受千年风华!

在科技日新月异的今天&#xff0c;虚拟现实&#xff08;VR&#xff09;技术已经成为了我们生活中不可或缺的一部分。从娱乐游戏到医疗健康&#xff0c;从教育培训到房地产销售&#xff0c;VR技术的应用领域日益广泛。而近年来&#xff0c;VR技术在文化遗产保护和古迹复原方面的…

2023年亚太杯APMCM数学建模大赛ABC题辅导及组队

2023年亚太杯APMCM数学建模大赛 ABC题 一元线性回归分析类 回归分析&#xff08;Regression Analysis)是确定两种或两种以上变量间相互依赖的定量关系的一种统计分析方法。   – 按涉及变量个数划分   • 一元回归分析   • 多元回归分析   – 按自变量和因变量之间关…

JS+CSS随机点名详细介绍复制可用(可自己添加人名)

想必大家也想拥有一个可以随机点名的网页&#xff0c;接下来我为大家介绍一下随机点名&#xff0c;可用于抽人&#xff0c;哈哈 <!DOCTYPE html> <html><head><meta charset"utf-8"><title></title><style>* {margin: 0;…

代码随想录 Day38 完全背包问题 LeetCode T70 爬楼梯 T322 零钱兑换 T279 完全平方数

前言 在今天的题目开始之前,让我们来回顾一下之前的知识,动规五部曲 1.确定dp数组含义 2.确定dp数组的递推公式 3.初始化dp数组 4.确定遍历顺序 5.打印dp数组来排错 tips: 1.当求取物品有限的时候用0-1背包,求取物品无限的时候用完全背包 结果是排列还是组合也有说法,当结果是组…

【Mybatis小白从0到90%精讲】12:Mybatis删除 delete, 推荐使用主键删除!

文章目录 前言XML映射文件方式推荐使用主键删除注解方式工具类前言 在实际开发中,我们经常需要删除数据库中的数据,MyBatis可以使用XML映射文件或注解来编写删除(delete)语句,下面是两种方法的示例。 XML映射文件方式 Mapper: int delete(int id);Mapper.xml:

抖音10月榜单有哪些看点?

10月20日&#xff0c;抖音双11好物节在抖音平台正式开启抢跑&#xff0c;据数据显示&#xff0c;截止10月31日平台多项双11销售增长记录再次被刷新。 *新抖双十一活动也已开启&#xff0c;最高可省30788元&#xff0c;活动详情&#x1f449; 抖音平台内大促氛围火爆&#xff0…

麒麟KYLINIOS软件仓库搭建02-软件仓库添加新的软件包

原文链接&#xff1a;麒麟KYLINIOS软件仓库搭建02-软件仓库添加新的软件包 hello&#xff0c;大家好啊&#xff0c;今天给大家带来麒麟桌面操作系统软件仓库搭建的文章02-软件仓库添加新的软件包&#xff0c;本篇文章主要给大家介绍了如何在麒麟桌面操作系统2203-x86版本上&…