导入STEP文件时边界框与三角化节点单位不匹配问题求助

使用PythonOCC 7.8.1.1从STEP文件导入轴对齐形状，几何外观显示正确，但计算出的边界圆柱高度符合预期，半径却远小于实际值。具体细节如下：

• 形状的边界框尺寸与STEP文件预期匹配（例如直径100对应坐标范围-50至+50）
• 通过BRep_Tool.Triangulation获取的三角化节点坐标范围仅为-8至+8，是边界框尺寸的1/6左右
• STEPControl_Reader.SystemLengthUnit()返回值为1.0，不确定该值的具体含义
• 最终计算出的圆柱半径约为8.2，远低于预期的50左右

请问该现象是否正常？如何解决？

原代码片段 #

step_reader = STEPControl_Reader()
status = step_reader.ReadFile(filename)
if status == IFSelect_RetDone:
    step_reader.TransferRoot(1)
    shape = step_reader.Shape(1)
    print(f"System Length Unit: {step_reader.SystemLengthUnit()}")  # 输出1.0
    # ...后续代码

# 为每个面获取三角化数据
tri = BRep_Tool.Triangulation(face, location)
# ...

# 计算圆柱半径的代码
result.radius = 0.0

for i in range(1, len(nodes) + 1):
    projNode = gp_Pnt(
        nodes.Value(i).X() * (1.0 - orientation.X()),
        nodes.Value(i).Y() * (1.0 - orientation.Y()),
        nodes.Value(i).Z() * (1.0 - orientation.Z()),
    )

    distance = projNode.Distance(projCenter)
    result.radius = max(distance, result.radius)
# ...

问题原因与解决方法 #

核心原因 #

BRep_Tool.Triangulation返回的节点坐标是面的局部坐标系下的坐标，而非全局坐标系。原代码未应用location参数对应的坐标转换，直接使用局部坐标计算半径，导致数值被缩小。

解决步骤 #

1. 将局部坐标转换为全局坐标：通过location.Transformation()获取转换矩阵，将三角化节点的局部坐标转换到全局坐标系
2. 修正投影节点计算逻辑：先完成坐标转换，再处理方向投影

修正后的代码片段 #

# 获取三角化数据时保留location
tri = BRep_Tool.Triangulation(face, location)
nodes = tri.Nodes()
# 获取坐标转换矩阵
transform = location.Transformation()

# 计算圆柱半径的代码
result.radius = 0.0

for i in range(1, len(nodes) + 1):
    # 先将局部点转换到全局坐标系
    global_node = nodes.Value(i).Transformed(transform)
    # 再处理方向投影
    projNode = gp_Pnt(
        global_node.X() * (1.0 - orientation.X()),
        global_node.Y() * (1.0 - orientation.Y()),
        global_node.Z() * (1.0 - orientation.Z()),
    )

    distance = projNode.Distance(projCenter)
    result.radius = max(distance, result.radius)

补充说明 #

• STEPControl_Reader.SystemLengthUnit()返回的1.0表示STEP文件使用的长度单位是米（OCC默认单位为米），如果你的STEP文件实际单位是毫米，需要额外做单位转换（乘以1000）。但从边界框尺寸匹配的情况看，单位问题不是本次半径偏小的原因。
• 若修正后仍有问题，需检查projCenter的计算是否基于全局坐标系，确保投影中心与节点坐标处于同一坐标系下。

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