3D线段与四边形平面相交检测结果异常的原因排查

在3D场景中，我用两点P0、P1表示线段起止点，四点ABCD构成四边形平面，需要检测线段是否与平面相交。参考MATLAB的plane_line_intersect函数，示例运行正常，但代入自定义坐标后出现异常：线段端点P1移出平面范围时，检测仍返回相交状态。

MATLAB代码示例 #

>> AB = B-A
AB = 0.8660    0.5000         0
>> AD = D-A
AD = 0     0     1
>> n = cross(AB,AD)/sqrt(dot(cross(AB,AD),cross(AB,AD)))
n = 0.5000   -0.8660         0
>> [I,check]=plane_line_intersect(n,A,P0,P1)
I = 1.0961   44.5116    6.6948
check = 3

自定义坐标数据 #

A = [1312.901701, 2360.12542566832, 74.24415425756794];
B = [1312.901701, 2423.58274087539, 88.31230234462595];
C = [1312.901701, 2371.241313396465, 24.103624956470995];
D = [1312.901701, 2434.698628603535, 38.171773043529];

P0 = [1495.468628, 2261.038086, 161.329498];
P1 = [1153.250854, 2479.341797, -14.787056];

对应的JavaScript实现代码 #

const ab = subtract(B, A);
const ad = subtract(D, A);
const n = divide(cross(ab, ad), sqrt(dot(cross(ab, ad), cross(ab, ad))));

const intersects = planeLineIntersect(n, A, P0, P1);

function planeLineIntersect(n, v0, p0, p1) {
    const u = subtract(p1, p0);
    const w = subtract(p0, v0);
    
    const D = dot(n, u);
    const N = -dot(n, w);

    let check = 0;

    if (abs(D) < 0.0000001) {
        if (N === 0) {
            check = 2;
            return check;
        }
        check = 0;
        return check;
    }

    const sI = divide(N, D);
    const I = add(p0, multiply(sI, u)); // not used

    if (sI < 0 || sI > 1) {
        check = 3;
    } else check = 1;

    return check;
}

异常原因分析 #

当前实现的核心问题是：只检测了线段所在直线与无限平面的相交关系，没有验证交点是否落在四边形的有限区域内。

• 函数中的sI判断仅确认交点在线段P0-P1的范围内（sI∈[0,1]时返回check=1），但完全没判断这个交点是否位于ABCD构成的四边形内部。
• 当线段所在直线与无限平面的交点在线段范围内，但该交点处于四边形的外部（无限平面的其他区域）时，函数仍会返回“相交”（check=1），这就是你看到的异常现象——实际上线段和无限平面相交，但和四边形平面并没有相交。

解决方向 #

要实现准确的线段与四边形相交检测，需要在计算出交点I后，添加额外的判断逻辑：

1. 将交点I投影到四边形所在的平面坐标系中；
2. 判断该投影点是否落在四边形的边界内。

常用的判断方法包括：

• 把四边形拆分为两个三角形（如ABC和ACD），分别判断交点是否在任意一个三角形内部；
• 针对凸四边形，使用向量叉积判断点是否在四边形所有边的同一侧；
• 使用重心坐标法验证点是否在四边形范围内。

内容的提问来源于stack exchange，提问作者yoranus