嘿，这个问题在3D碰撞检测管线里太常见了！我用GLM做过不少类似的实现，正好给你捋清楚完整的计算方法，结合你提到的向量运算一步步拆解：

核心思路：向量投影 + 范围限制 #

本质上，我们要找的是点C在线段AB所在直线上的投影点，然后把这个投影点限制在线段AB的范围内，就得到了最近点D。

步骤分解（结合GLM向量操作） #

1. 定义基础向量
   先计算线段AB的方向向量，以及从线段起点A到点C的向量：

   glm::vec3 AB = B - A; // 线段AB的方向向量
   glm::vec3 AC = C - A; // 从A指向C的向量
   

2. 计算投影比例t
   t是向量AC在AB上的投影占AB长度的比例，用点积计算：

   float lengthSquaredAB = glm::dot(AB, AB); // 计算AB长度的平方（避免开根号，提升性能）
   

   这里要先判断lengthSquaredAB是否接近0（比如小于GLM的epsilon），如果是，说明A和B重合，最近点直接返回A即可，避免除以0的错误。

   正常情况下计算t：

   float t = glm::dot(AC, AB) / lengthSquaredAB;
   

3. 限制t的范围
   因为我们要的是线段上的点，不是无限长直线上的点，所以t必须限制在[0, 1]之间：

   • 当t < 0时，最近点是线段起点A
   • 当t > 1时，最近点是线段终点B
   • 当t在0-1之间时，最近点是投影点

   用GLM的clamp函数直接处理：

   t = glm::clamp(t, 0.0f, 1.0f);
   

4. 计算最近点D
   用t和向量AB计算出最终的最近点：

   glm::vec3 D = A + t * AB;
   

完整的GLM实现函数 #

把上面的步骤封装成可复用的函数：

#include <glm/glm.hpp>
#include <glm/gtc/constants.hpp> // 用于glm::epsilon

glm::vec3 getClosestPointOnSegment(const glm::vec3& segmentStart, const glm::vec3& segmentEnd, const glm::vec3& point) {
    glm::vec3 segmentDir = segmentEnd - segmentStart;
    float segmentLengthSq = glm::dot(segmentDir, segmentDir);

    // 处理线段长度为0的特殊情况（两点重合）
    if (segmentLengthSq < glm::epsilon<float>()) {
        return segmentStart;
    }

    glm::vec3 pointToStart = point - segmentStart;
    float t = glm::dot(pointToStart, segmentDir) / segmentLengthSq;

    // 限制t在线段范围内
    t = glm::clamp(t, 0.0f, 1.0f);

    // 计算并返回最近点
    return segmentStart + t * segmentDir;
}

关于你提到的叉积 #

你说的叉积确实有用，但更多是用来计算点到线段的距离，或者判断点相对于线段的位置：

• 点到直线AB的距离可以用glm::length(glm::cross(segmentDir, pointToStart)) / glm::length(segmentDir)
• 但只有当t在0-1之间时，这个距离才等于点到线段的距离；否则距离是点到A或点到B的直线距离

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