非轴对齐系统中Xbim构件距离计算及碰撞检测技术咨询

一、替换AABB为定向包围盒（OBB）优化粗筛选 #

针对长距离非轴对齐风管，轴对齐包围盒（AABB）会包含大量无关空间，导致粗筛精度极低。改用定向包围盒（OBB），完全贴合构件自身坐标系，能大幅减少候选构件数量：

• 从风管构件（如IfcDuctSegment）获取LocalPlacement对应的全局变换Matrix3D，这个矩阵包含了构件的旋转、平移、缩放信息
• 获取构件局部坐标系下的边界盒（通过GetBoundingBox()得到局部Min/Max值）
• 将局部边界盒的8个顶点通过Matrix3D变换到全局坐标系，构建出与风管自身对齐的OBB
• 可给OBB设置小范围的扩展阈值（比如5cm），避免因模型精度问题漏检，再判断目标点是否在扩展后的OBB范围内，快速排除无关风管

二、精确计算点到风管的最小距离 #

基于已获取的Matrix3D、面、垂线等信息，分两种风管类型处理：

1. 规则线性风管（矩形/圆形截面，直线路径） #

这类风管在局部坐标系下路径通常沿某一轴（如Z轴），截面在XY平面内，计算步骤如下：

• 用构件Matrix3D的逆矩阵将全局坐标点变换到局部坐标系：var localPoint = Matrix3D.Invert(matrix).Transform(globalPoint);
• 计算局部点到风管路径直线的垂直距离：利用向量叉乘公式Vector3D.CrossProduct(localPoint - lineStart, lineDir).Length / lineDir.Length，其中lineStart是风管局部起点，lineDir是路径方向向量
• 判断垂直距离与截面尺寸的关系：
  • 若垂直距离 ≤ 截面半宽（矩形）/半径（圆形），再检查点在路径上的投影是否在风管长度范围内：同时满足则点在风管内部，距离为0
  • 若垂直距离 > 截面尺寸，计算“垂直距离 - 截面尺寸”得到点到风管侧壁的距离；若投影超出风管长度范围，额外计算点到风管两端面的距离，取所有结果中的最小值

2. 复杂弯曲风管 #

对于带弯头、变径的复杂风管，需拆解为基础几何单元计算：

• 通过Xbim的IModelGeometry获取风管的三角化网格或NURBS曲面数据
• 遍历每个几何单元（线段、三角面），计算点到单元的最小距离：
  • 点到线段的距离：用向量投影判断投影点是否在线段上，再计算垂直距离或端点距离
  • 点到三角面的距离：先判断点在面内的投影是否在三角面范围内，再计算垂直距离，否则取到三个顶点的最小距离
• 所有几何单元的最小距离即为点到风管的距离

三、XbimSpatialAnalyser的后续搜索逻辑说明 #

XbimSpatialAnalyser中AABB粗筛后的核心步骤是精确几何校验：

1. 对OctTree粗筛出的候选构件，调用几何相交/距离检测方法（如Intersects、DistanceTo）
2. 默认用AABB会导致候选集过大，替换为OBB后，能直接减少进入精确校验的构件数量
3. 可参考FindElementsWithinDistance方法的源码逻辑，将其中的AABB检测替换为OBB检测，再保留后续的精确几何计算流程

四、关键代码片段示例 #

// 获取风管的全局变换矩阵
var matrix = duct.ObjectPlacement?.GetMatrix();
if (matrix == null) return float.MaxValue;

// 将全局点转换到局部坐标系
var inverseMatrix = Matrix3D.Invert(matrix);
var localPoint = inverseMatrix.Transform(globalPoint);

// 假设风管局部路径沿Z轴，起点(0,0,0)，终点(0,0,length)，截面为矩形半宽w/2，半高h/2
var lineDir = new Vector3D(0, 0, 1);
var lineStart = new Point3D(0, 0, 0);
var lineEnd = new Point3D(0, 0, ductLength);

// 计算点到路径直线的垂直距离
var vecToStart = localPoint - lineStart;
var projLength = Vector3D.DotProduct(vecToStart, lineDir);
var projPoint = lineStart + lineDir * projLength;
var perpendicularDist = (localPoint - projPoint).Length;

// 判断投影点是否在风管长度范围内
bool isProjInRange = projLength >= 0 && projLength <= ductLength;

float minDistance = float.MaxValue;
if (isProjInRange)
{
    // 点在风管路径范围内，计算到侧壁的距离
    if (perpendicularDist <= Math.Max(w/2, h/2))
        minDistance = 0; // 点在内部
    else
        minDistance = perpendicularDist - Math.Max(w/2, h/2);
}
else
{
    // 投影超出范围，计算到两端点的距离
    var distToStart = (localPoint - lineStart).Length;
    var distToEnd = (localPoint - lineEnd).Length;
    minDistance = Math.Min(distToStart, distToEnd) - Math.Max(w/2, h/2);
    minDistance = Math.Max(minDistance, 0); // 避免负数（点在端面内部）
}

// 转换回全局坐标系下的距离（因局部和全局缩放一致，无需额外转换）
return minDistance;

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