凸多边形间最小距离

给定两个非连接(比如不相交)的凸多边形 P 和 Q, 目标是找到拥有最小距离的点对 (p,q) (p 属于 P 且 q 属于Q)。 


事实上, 多边形非连接十分重要, 因为我们所说的多边形包含其内部。 如果多边形相交, 那么最小距离就变得没有意义了。 然而, 这个问题的另一个版本, 凸多边形顶点对间最小距离对于相交和非相交的情况都有解存在。 

回到我们的主问题: 直观的, 确定最小距离的点不可能包含在多边形的内部。 与最大距离问题相似, 我们有如下结论: 

两个凸多边形  P  和  Q  之间的最小距离由 多边形间的对踵点对 确立。 存在凸多边形间的三种 多边形间的对踵点对 , 因此就有三种可能存在的最小距离模式:
  1. “顶点-顶点”的情况
  2. “顶点-边”的情况
  3. “边-边”的情况

换句话说, 确定最小距离的点对不一定必须是顶点。 下面的三个图例表明了以上结论: 
     
给定结果, 一个基于旋转卡壳的算法自然而然的产生了: 
考虑如下的算法, 算法的输入是两个分别有 m 和 n 个顺时针给定顶点的凸多边形 P 和 Q。
  1. 计算 P 上 y 坐标值最小的顶点(称为 yminP ) 和 Q 上 y 坐标值最大的顶点(称为 ymaxQ)。
  2. 为多边形在 yminP 和 ymaxQ 处构造两条切线 LP 和 LQ 使得他们对应的多边形位于他们的右侧。 此时 LP和 LQ 拥有不同的方向, 并且 yminP 和 ymaxQ 成为了多边形间的一个对踵点对。
  3. 计算距离(yminP,ymaxQ) 并且将其维护为当前最小值。
  4. 顺时针同时旋转平行线直到其中一个与其所在的多边形的边重合。
  5. 如果只有一条线与边重合, 那么只需要计算“顶点-边”对踵点对和“顶点-顶点”对踵点对距离。 都将他们与当前最小值比较, 如果小于当前最小值则进行替换更新。 如果两条切线都与边重合, 那么情况就更加复杂了。 如果边“交叠”, 也就是可以构造一条与两条边都相交的公垂线(但不是在顶点处相交), 那么就计算“边-边”距离。 否则计算三个新的“顶点-顶点”对踵点对距离。 所有的这些距离都与当前最小值进行比较, 若小于当前最小值则更新替换。
  6. 重复执行步骤4和步骤5, 直到新的点对为(yminP,ymaxQ)。
  7. 输出最大距离。
旋转卡壳模式保证了所有的对踵点对(和所有可能的子情况)都被考虑到。 此外, 整个算法拥有现行的时间复杂度, 因为(除了初始化), 只有与顶点数同数量级的操作步数需要执行。 

最小距离和最大距离的问题表明了旋转卡壳模型可以用在不同的条件下(与先前的直径和宽度问题比较)。 这个模型可以应用于两个多边形的问题中。 
“最小盒子”问题( 最小面积外接矩形 )通过同一多边形上两个正交切线集合展示了另一种条件下旋转卡壳的应用。  

 

 

原文地址:http://cgm.cs.mcgill.ca/~orm/mind2p.html


  • 2
    点赞
  • 10
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
旋转卡壳法(Rotating Calipers Algorithm),也称为圆规算法,是一种用于计算二维凸多边形内最近点对之间最短距离的有效算法。在MATLAB中,你可以通过编写一个函数来实现这个过程。以下是基本步骤: 1. **初始化**:选择两个顶点作为初始“圆规”(通常取顶点集合的第一和最后一个元素),记为`p1`和`p2`。 2. **迭代**:对于每个顶点`p`,计算它与当前圆规对角线所围成的三角形的两边长度`d1`和`d2`。 - 如果`d1`+`d2`小于当前最短距离,更新最短距离,并调整圆规为新形成的对角线。 - 否则,将圆规顺时针移动到下一个顶点,除非即将形成闭合路径,这时需要逆时针回溯。 3. **结束条件**:当圆规再次回到起始位置时,循环结束,因为这意味着所有可能的对角线都已被考虑过。 4. **返回结果**:最短距离即为找到的最小距离。 在MATLAB中,可以创建一个函数,例如`minDistancePolyhedron`,接受凸多边形的顶点列表作为输入,然后按照上述逻辑遍历并计算。下面是一个简化的示例代码片段: ```matlab function d = rotatingCalipers(p) n = length(p); if n < 3 % 凸多边形最少要有三个顶点 error('Not a valid polygon'); end p1 = p(1); % 初始圆规 p2 = p(end); minDist = Inf; % 初始化最小距离 for i = 2:n-1 d1 = norm(p(i) - p1); d2 = norm(p(i+1) - p1); if d1 + d2 < minDist [d, ~] = min([minDist, d1, d2]); p2 = p(i+1); else if i == n-1 % 如果形成闭合路径,回溯到开始 p2 = p1; else p2 = p(i); end end end d = minDist; % 返回最小距离 end % 示例使用 vertices = [0 0; 1 0; 1 1; 0 1]; % 四边形 minDistance = rotatingCalipers(vertices); ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值