python - 从 OpenCascade 表面获取 (t, c, k) 值

Question

我创建了一个用于在 Python 中创建和使用 b 样条曲面的库，利用并行 scipy.interpolate.RectBivariateSpline() 实例来保存结向量 (X, Y, Z)控制点网格，以及 u 和 v 中的度数(执行表面评估的 (t, c, k) 元组)。我还编写了一个 STEP 解析器来读取从 CAD 包导出的表面数据；我从文件中的 b_spline_surface_with_knots 实体中获取 (t, c, k) 值，并将它们填充到我自己的对象中。表面库对我来说工作得很好，但 STEP 解析器很痛苦，几乎每次我使用它时都会以各种方式失败。所以我尝试使用“真正的”STEP 解析器，如下所示:

from OCC.STEPControl import STEPControl_Reader
from OCC.IFSelect import IFSelect_RetDone, IFSelect_ItemsByEntity

step_reader = STEPControl_Reader()
status = step_reader.ReadFile('c:/LPT/nomdata/lpt3.stp')

if status == IFSelect_RetDone:  # check status
    failsonly = False
    step_reader.PrintCheckLoad(failsonly, IFSelect_ItemsByEntity)
    step_reader.PrintCheckTransfer(failsonly, IFSelect_ItemsByEntity)

    ok = step_reader.TransferRoot(1)
    _nbs = step_reader.NbShapes()
    aResShape = step_reader.Shape(1)
else:
    print("Error: can't read file.")
    sys.exit(0)

现在我有了这个 aResShape 对象，但是在 IPython 中(也不是谷歌搜索)没有多少戳和刺激它揭示了如何获得定义表面的 (t, c, k) 值。

有人能告诉我揭示这些值的方法吗？或者是否可能存在另一种不那么不透明的基于 Python 的 STEP 解析器？

Answer 1

这个问题有点老了，但以防万一其他人遇到类似的问题...

step_reader.Shape()的结果是一个TopoDS_Shape，它是一个拓扑实体，可以分为以下组件拓扑:

• 顶点——对应于几何中的一个点的零维形状；
• 边——对应于曲线的形状，并由每个末端的顶点限制；
• 连线——由顶点连接的边序列；
• 面 – 由闭合线界定的平面(2D 几何)或表面(3D 几何)的一部分；
• 壳——由线边界的某些边缘连接的面的集合；
• Solid – 3D 空间的一部分，由壳包围；
• 复合固体——固体的集合。

通常，您会使用 TopoDS_Shape::ShapeType() 方法查询它，以便知道那个形状是什么(顶点？边？，...)。

如果模型是由单个 b 样条曲面形成的，那么形状应该是 TopoDS_Face，您可以通过调用获得:

face = aResShape.Face();

一旦手边有了 TopoDS_Face，就可以像这样获取底层几何体 (Geom_Surface):

surface = BRepAdaptor_Surface(face).Surface().BSpline();

现在您已经可以访问底层几何图形，您可以调用该对象的方法，它们将为您提供所需的信息。

它们记录在这里:

• > https://www.opencascade.com/doc/occt-7.1.0/refman/html/class_geom___b_spline_surface.html

OpenCASCADE 文档可能看起来很困惑，但我认为您可能对这个主题感兴趣:

• > https://www.opencascade.com/doc/occt-7.0.0/overview/html/occt_user_guides__modeling_data.html#occt_modat_3

希望对您有所帮助。