Möller-Trumbore算法-射线三角形相交算法

最新推荐文章于 2024-11-18 18:11:20 发布

zhanxi1992

最新推荐文章于 2024-11-18 18:11:20 发布

阅读量9.9k

点赞数 29

分类专栏：图形学文章标签：图形学几何学线性代数算法

本文链接：https://blog.csdn.net/zhanxi1992/article/details/109903792

版权

图形学专栏收录该内容

19 篇文章

订阅专栏

Möller-Trumbore算法

一、概述
二、准备知识
三、Möller-Trumbore 算法推导
- 推导过程
四、代码实现

一、概述

Möller-Trumbore 射线三角相交算法是一种快速计算射线与三角形在三个维度上的交点的方法，通过向量与矩阵计算可以快速得出交点与重心坐标，而无需对包含三角形的平面方程进行预计算。另外，它还应用于计算机图形学中以实现涉及三角形网格的光线跟踪计算。算法名字是以发明者 TomasMöller 和 Ben Trumbore 的名字来命名的。

二、准备知识

1、三阶方阵的行列式等于三个列向量的混合积
在这里插入图片描述

2、克莱姆法则（部分译克拉默法则）
如果一个线性方程组 $\boldsymbol{Ax} = \boldsymbol{c}$ , 其中 $\boldsymbol{A}$ 是可逆方阵， $\boldsymbol{x}$ ， $\boldsymbol{c}$ 都是列向量，那么方程有解，且 $\boldsymbol{x}$ 的每一个解
在这里插入图片描述

其中 $A_i$ 是被列向量取代了第 $i$ 列的矩阵。

三、Möller-Trumbore 算法推导

已知光线 $\boldsymbol{Ray} = \boldsymbol{O} + t\boldsymbol{D}$ （ $\boldsymbol{O}$ 为起点， $\boldsymbol{D}$ 为射线方向， $t$ 为时间）, 三角形三个顶点 $P_0$ ， $P_1$ ， $P_2$ 。光线与三角形相交时，可得如下等式：
在这里插入图片描述

则可以解
在这里插入图片描述

参数定义：
在这里插入图片描述

推导过程

从这里开始
在这里插入图片描述

括号展开，移项可得
在这里插入图片描述

观察一下上面的括号以及等式左边的内容，都是已知的点，因此点的加减可以用向量来表示，令
在这里插入图片描述

得到
在这里插入图片描述

也即
在这里插入图片描述

这是一个形如 $\boldsymbol{Ax} = \boldsymbol{c}$ 的等式，所以可以用克拉莫法则
在这里插入图片描述

由向量混合积可以得出
分母部分：
在这里插入图片描述

令
在这里插入图片描述

原式等于
在这里插入图片描述

分子部分：
在这里插入图片描述
令

原式等于
在这里插入图片描述

因此
在这里插入图片描述

另外的两个参数 $b_1$ ， $b_2$ 可以同样推出来
在这里插入图片描述

四、代码实现

bool rayTriangleIntersect(const Vector3f& v0, const Vector3f& v1, const Vector3f& v2, const Vector3f& orig,
                          const Vector3f& dir, float& tnear, float& u, float& v)
{
    bool isIn = false;
    Vector3f E1 = v1 - v0;
    Vector3f E2 = v2 - v0;
    Vector3f S = orig - v0;
    Vector3f S1 = crossProduct(dir, E2);
    Vector3f S2 = crossProduct(S, E1);
    float coeff = 1.0 / dotProduct(S1, E1); // 共同系数
    float t = coeff * dotProduct(S2, E2);
    float b1 = coeff * dotProduct(S1, S);
    float b2 = coeff * dotProduct(S2, dir);
    if (t >= 0 && b1 >= 0 && b2 >= 0 && (1 - b1 - b2) >= 0)
    {
        isIn = true;
        tnear = t;
        u = b1;
        v = b2;
    }

    return isIn;
}