✨你在我的四面体心房吗？——判断点是否在四面体内的图形学魔法

"他站在空间的边界线上，不知是否已进入了那个被四个顶点围起的微型宇宙。"

🌐 前言：图形学，不只是线和面

计算机图形学不仅仅是炫酷的模型和光影，它更像是空间中的逻辑诗人，衡量一切点与线、面与体之间的微妙关系。而今天我们要解开一个经典谜题：

如何判断一个点，是否真的"在"一个四面体内部？

别急，我们不靠玄学，而是靠 行列式 和 体积的符号，就像用四杯水测试一个陌生人是否喝醉------每杯代表一个角度，每次的反应就是一个提示。

🧠 数学直觉：四个体积，揭示一个真相

给定一个四面体，由点 A、B、C、D 构成，我们想判断点 P 是否处在这个四面体内部。我们将 P 与四面体的四个面分别组成 四个新的四面体：

PBCD
PACD
PABD
PABC

这些新四面体的体积，与原始四面体 ABCD 的体积，是判断 P 是否在其中的关键。

🧊 魔法原理：体积符号不变即在其中

我们把每个体积都看成一个有符号的量，如果 P 在四面体 ABCD 内部，那么这五个体积（原始体积 + 四个子体积）将拥有相同的符号。

换句话说：

全部为正，P 在里面。
全部为负，P 仍在里面（换个坐标系看而已）。
如果有异号，说明 P 不在四面体内部，它踩出界了。

🧮 编程实现：JS版的行列式魔术

下面我们不写公式，用向量和行列式函数来实现这一判断：

css 复制代码

// 计算 3D 向量叉积
function cross(a, b) {
  return [    a[1] * b[2] - a[2] * b[1],
    a[2] * b[0] - a[0] * b[2],
    a[0] * b[1] - a[1] * b[0]
  ];
}

// 计算 3D 向量点积
function dot(a, b) {
  return a[0] * b[0] + a[1] * b[1] + a[2] * b[2];
}

// 从四个点计算有符号体积（其实是六倍体积）
function signedVolume(a, b, c, d) {
  const ab = [b[0]-a[0], b[1]-a[1], b[2]-a[2]];
  const ac = [c[0]-a[0], c[1]-a[1], c[2]-a[2]];
  const ad = [d[0]-a[0], d[1]-a[1], d[2]-a[2]];
  return dot(cross(ab, ac), ad);
}

// 判断点 P 是否在四面体 ABCD 中
function isPointInTetrahedron(p, a, b, c, d) {
  const v0 = signedVolume(a, b, c, d);
  const v1 = signedVolume(p, b, c, d);
  const v2 = signedVolume(a, p, c, d);
  const v3 = signedVolume(a, b, p, d);
  const v4 = signedVolume(a, b, c, p);

  const sameSign = (x, y) => (x >= 0 && y >= 0) || (x <= 0 && y <= 0);

  return (
    sameSign(v0, v1) &&
    sameSign(v0, v2) &&
    sameSign(v0, v3) &&
    sameSign(v0, v4)
  );
}

🎲 举个栗子：点落入四面体的测试

ini 复制代码

const A = [0, 0, 0];
const B = [1, 0, 0];
const C = [0, 1, 0];
const D = [0, 0, 1];

const P_inside = [0.1, 0.1, 0.1];
const P_outside = [1, 1, 1];

console.log(isPointInTetrahedron(P_inside, A, B, C, D)); // true
console.log(isPointInTetrahedron(P_outside, A, B, C, D)); // false

🔍 拓展视角：为什么用行列式？

本质上，我们在用**混合积（Mixed Product）**来测体积，而混合积就等于行列式！你不必真的写出那个花里胡哨的 3x3 矩阵，但你心里要知道，每一个体积的背后，其实是：

行列式在说话，向量在歌唱。

🧬 总结：四面体内部，是空间中的忠诚试炼

判断一个点是否在四面体内部，看似只是几行代码，其实是对三维空间几何关系的深刻洞察。我们不只是处理点和面，而是在用向量和体积讲述"空间归属感"的故事。

所以下次当你写 isPointInTetrahedron() 的时候，不妨温柔地想一句：

"你，真的属于这个四面体吗？"

🚀 附加挑战

将函数扩展为支持浮点误差容忍（例如加一个 epsilon 参数）
可视化每个体积的大小，用颜色表示符号
把这个判断用于物理仿真、碰撞检测或网格剖分中