碰撞检测算法之闵可夫斯基差集法(Minkowski Difference)

news2024/12/23 8:56:40

在游戏开发和机器人路径规划乃至于现在比较火的自动驾驶中,我们常常需要确定两个物体是否发生碰撞,有一种通过闵可夫斯基差集法求是否相交的算法,下面将介绍一下

闵可夫斯基差集法的优势

闵可夫斯基差集法优势:

  • 可以处理复杂的形状和多边形/多面体。
  • 提供精确的接触点和分离向量信息。

AABB的优势:

  • 不像AABB一样简单快速,只需要几次比较就可以完成相交测试
  • 可生成AABBTree,快速排除不可能相交的物体

可以看出AABB更快,闵可夫斯基差法更准,所以可以拿AABBTree缩小检测目标范围,最后精准检测的时候可以用闵可夫斯基差集法

闵可夫斯基差集法图形求交

有两个相交的图形A和B

在A中取一个点A1,和B中的所有点求差

将图形连接起来

再在A中取一个点A2,和B中的所有点求差

取A3,和B中的所有点求差

(其实在这里,原点出现在了生成的图形中,就可以终止循环了,我后面优化部分会提到)

取A4,和B中的所有点求差

连接所有结果

可以看到原点在闵可夫斯基差集的图形中,则成两个图形存相交情况

如果没有交点,也就不存在焦点了

测试代码

import numpy as np

def minkowski_difference(A, B):
    """计算两个多边形A和B的闵可夫斯基差"""
    diff = []
    for a in A:
        for b in B:
            diff.append(a - b)
    return np.array(diff)

def is_origin_inside(polygon):
    """检查原点是否在多边形内部"""
    return true/false

def polygons_collide(A, B):
    """判断两个多边形是否碰撞"""
    # 计算闵可夫斯基差
    difference = minkowski_difference(A, B)
    
    # 检查原点是否在差集内
    return is_origin_inside(difference)

# 示例:定义两个多边形的顶点
polygon_A = np.array([[0, 0], [2, 0], [2, 2], [0, 2]])  # 正方形
polygon_B = np.array([[1, 1], [3, 1], [3, 3], [1, 3]])  # 另一个正方形,与A相交

# 检测碰撞
collision = polygons_collide(polygon_A, polygon_B)
print("Polygons collide:", collision)

优化:

①闵可夫斯基差的几何意义是所有可能从一个多边形内的任一点到另一个多边形内任一点的向量集合。而上面的方法生成的是一个点集,而不是这些点构成的实际几何形状(通常是凸包),因此,许多计算出来的点对于确定原点是否在闵可夫斯基差内部是没有必要的

我上边到A3遍历的时候,就已经知道相交了,就直接可以是否相交了,没必要再计算A4了

要避免不必要的计算,我们可以优化minkowski_difference函数,使其不生成完整的闵可夫斯基差集。我们可以在检测过程中动态地生成差异点,并且一旦发现原点位于闵可夫斯基差内,就可以立即返回结果,而不需要继续计算,简单优化下

import numpy as np

def minkowski_difference_check_origin(A, B):
    """检查原点是否在两个多边形A和B的闵可夫斯基差内"""
    for a in A:
        for b in B:
            diff = a - b
            if is_origin_inside_single_point(diff):
                return True
    return False

def is_origin_inside_single_point(point):
    """简单检查单个点是否非常接近原点"""
    return true/false

def polygons_collide(A, B):
    """判断两个多边形是否碰撞"""
    # 动态检查原点是否在闵可夫斯基差内
    return minkowski_difference_check_origin(A, B)

# 示例:定义两个多边形的顶点
polygon_A = np.array([[0, 0], [2, 0], [2, 2], [0, 2]])  # 正方形
polygon_B = np.array([[1, 1], [3, 1], [3, 3], [1, 3]])  # 另一个正方形,与A相交

# 检测碰撞
collision = polygons_collide(polygon_A, polygon_B)
print("Polygons collide:", collision)

当然了,这个碰撞算法不是最佳的,有更加好的GJK算法,是一种优化后的技术,它巧妙地利用了闵可夫斯基差的概念,但在实际应用中避免了直接计算其所有点所带来的高昂代价

后面我将会更新GJK算法的文章

补充

注意是Minkowski Difference(闵可夫斯基差集)不是Minkowski Distance(闵可夫斯基距离),虽然都来源于数学家赫尔曼·闵可夫斯基的工作,但它们是两个不同的概念,应用领域也有所不同,这里注意区分

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

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

相关文章

Python OCR 文字识别

一.引言 文字识别,也称为光学字符识别(Optical Character Recognition, OCR),是一种将不同形式的文档(如扫描的纸质文档、PDF文件或数字相机拍摄的图片)中的文字转换成可编辑和可搜索的数据的技术。随着技…

【系统】Windows11更新解决办法,一键暂停

最近的windows更新整的我是措不及防,干啥都要关注一下更新的问题,有的时候还关不掉,我的强迫症就来了,非得关了你不可! 经过了九九八十一难的研究之后,终于找到了一个算是比较靠谱的暂停更新的方法&#x…

复合翼与倾转旋翼飞行器:设计与控制算法对比

一、引言 复合翼(Compound Wing)和倾转旋翼(Tilt - Rotor)飞行器在现代航空领域均占据独特地位,二者在设计和控制算法方面展现出显著差异。这些差异在飞行模式切换、推进系统设计、控制算法复杂度以及飞行器稳定性等多…

空闲中断配合DMA

1.传统串口接收数据:来一个字节接受一个。 2.一次中断将一包数据存到缓冲区 3.DMA原理

三格电子——新品IE103转ModbusTCP网关

型号:SG-TCP-IEC103 产品概述 IE103转ModbusTCP网关型号SG-TCP-IEC103,是三格电子推出的工业级网关(以下简称网关),主要用于IEC103数据采集、DLT645-1997/2007数据采集,IEC103支持遥测和遥信,可…

HDLBits训练3

时间:2024.12.22 Hadd 代码 法一: module top_module( input a, b,output cout, sum );assign {cout,sum}ab; endmodule法二: 运行结果 Fadd 代码 法一: module top_module( input a, b, cin,output cout, sum );assign…

Qt之串口设计-线程实现(十二)

Qt开发 系列文章 - Serial-port(十二) 目录 前言 一、SerialPort 二、实现方式 1.创建类 2.相关功能函数 3.用户使用 4.效果演示 5.拓展应用-实时刷新 总结 前言 Qt作为一个跨平台的应用程序开发框架,在串口编程方面提供了方便易用…

STM32F407ZGT6-UCOSIII笔记12: 事件标志组

有时一个任务需要与多个事件同步,这就要用到事件标志组 本文学习与程序编写基于 正点原子的 STM32F1 UCOS开发手册 文章提供测试代码讲解、完整工程下载、测试效果图 目录 事件标志组: 定义与初始化事件标志组: #include "Public.h&quo…

聊一聊 C#前台线程 如何阻塞程序退出

一:背景 1. 讲故事 这篇文章起源于我的 C#内功修炼训练营里的一位朋友提的问题:后台线程的内部是如何运转的 ? ,犹记得C# Via CLR这本书中 Jeffery 就聊到了他曾经给别人解决一个程序无法退出的bug,最后发现是有一个 Backgrond…

JVM性能优化一:初识内存泄露-内存溢出-垃圾回收

本文主要是让你充分的认识到什么叫做内存泄露,什么叫做内存溢出,别再傻傻分不清了,别再动不动的升级服务器的内存了。 文章目录 1.基本概念1.1.内存泄露1.2.内存溢出1.3.垃圾回收1.4.内存泄露-垃圾回收-内存溢出三者的关系关系 2.代码示例2.…

为什么使用环形队列

1.看以下两种情况。第一种不会出现问题,当主流程读取次数比较慢时,数据会被覆盖。 2.扩大空间。不可取。 3.什么是队列

【WRF教程第3.6期】预处理系统 WPS 详解:以4.5版本为例

预处理系统 WPS 详解:以4.5版本为例 Geogrid/Metgrid 插值选项详解1. 插值方法的工作机制2. 插值方法的详细说明2.1 四点双线性插值(four_pt)2.2 十六点重叠抛物线插值(sixteen_pt)2.3 简单四点平均插值(av…

批量提取zotero的论文构建知识库做问答的大模型(可选)——含转存PDF-分割统计PDF等

文章目录 提取zotero的PDF上传到AI平台保留文件名代码分成20个PDF视频讲解 提取zotero的PDF 右键查看目录 发现目录为 C:\Users\89735\Zotero\storage 写代码: 扫描路径‘C:\Users\89735\Zotero\storage’下面的所有PDF文件,全部复制一份汇总到"C:\Users\89735\Downl…

Java模拟Mqtt客户端连接Mqtt Broker

Java模拟Mqtt客户端基本流程 引入Paho MQTT客户端库 <dependency><groupId>org.eclipse.paho</groupId><artifactId>org.eclipse.paho.mqttv5.client</artifactId><version>1.2.5</version> </dependency>设置mqtt配置数据 …

boost asio 异步服务器

boost网络框架使用方法 boost绑定 首先介绍io_context&#xff0c;可以理解为这是操作系统和应用层数据交互的桥梁。有了它不必关注内核态的缓冲区&#xff0c;只需要关注自己定义在用户态的缓冲区&#xff0c;因为它会通过桥梁运输到用户态的缓冲区。 boost::asio::io_contex…

图解HTTP-HTTP协议

HTTP HTTP是一种不保存状态&#xff0c;即无状态的协议。HTTP协议自身不对请求和响应之间的通信进行保存。为了保存状态因此后面也有一些技术产生比如Cookies技术。 HTTP是通过URI定位网上的资源&#xff0c;理论上将URI可以访问互联网上的任意资源。 如果不是访问特定的资源…

【Go】-限流器的四种实现方法

目录 关于限流和限流器 固定窗口限流器 滑动窗口限流器 漏桶限流器 令牌桶限流器 总结 关于限流和限流器 限流&#xff08;Rate Limiting&#xff09;是一种控制资源使用率的机制&#xff0c;通常用于防止系统过载和滥用。 限流器&#xff08;Rate Limiter&#xff09;是…

CTF_1

CTF_Show 萌新赛 1.签到题 <?php if(isset($_GET[url])){system("curl https://".$_GET[url].".ctf.show"); }else{show_source(__FILE__); }?> 和 AI 一起分析 1.if(isset($_GET[url]))检查GET请求中是否存在名为url的参数。 curl 2.curl…

[文献阅读] Unsupervised Deep Embedding for Clustering Analysis (无监督的深度嵌入式聚类)

文章目录 Abstract:摘要聚类深度聚类 KL散度深度嵌入式聚类(DEC)KL散度聚类软分配&#xff08;soft assignment&#xff09;KL散度损失训练编码器的初始化聚类中心的初始化 实验评估总结 Abstract: This week I read Unsupervised Deep Embedding for Clustering Analysis .It…

记录:virt-manager配置Ubuntu arm虚拟机

virt-manager&#xff08;Virtual Machine Manager&#xff09;是一个图形用户界面应用程序&#xff0c;通过libvirt管理虚拟机&#xff08;即作为libvirt的图形前端&#xff09; 因为要在Linux arm环境做测试&#xff0c;记录下virt-manager配置arm虚拟机的过程 先在VMWare中…