Topology and Geometry in OpenCascade

Topology and Geometry in OpenCascade-Edge eryar@163.com 摘要 Abstract: 本文简要介绍了几何造型中的边界表示法 ( BRep), 并结合程序说明 OpenCascade 中的边界表示的具体实现, 即拓朴与几何的联系 对具有几何信息的拓朴结构 顶点 (vertex) 边 (edge) 面 (face) 进行了详细说明 本文只对拓朴边数据进行说明, 着 重介绍了 OpenCascade 中两个种特别的边缝合边 (seam edge) 和退化边 (degenerated edge) 关键字 Key Words:OpenCascade BRep Topology Edge Geometry 一 引言 Introduction 边 (Edge) 是两个邻面 ( 对正则形体而言 ) 或多个邻面 ( 对非正则形体而言 ) 的交集 边有方向, 它由起始顶点和终止顶点来界定 ; 边的形状由边的几何信息来表示, 可以是直线, 也可以是曲线, 曲线边可用一系列控制点或型值点来描述, 也可以用显示 隐式或参数方程 来描述 边 (Edge) 是边界表示法中的重要结构, 因为边界表示法 (BRep) 是用形体的边界来 描述形体的一种方法 BRep 认为形体是由有限数量的边界表面 ( 平面或曲面 ) 构成, 而每 个表面又由若干边界边与顶点构成, 所有的单元面构成了形体的边界, 形体的边界将形体和周围的环境分隔开来 边界表示法不仅详细记录了构成形体的面 边方程的系数和顶点坐标值的几何信息, 而 且描述了这些几何元素之间的拓朴信息, 即体 面 边 顶点的组成关系等 在保证对形体 面的定义确定并且无二义性的前提下, 它允许根据形体的拓朴结构 面表示的方便性等因素 确定一个面是以一个整体表示, 还是以几个部分之和进行表示 在 OpenCascade 中边包含了一系列的曲线, 其结构如下图所示 : BRep_TEdge +Standard_Real mytolerance +Standard_Integer myflags +BRep_ListOfCurveRepresentation mycurves Figure 1.1 BRep_TEdge members 其中, 包含一系列的曲线由下面几种 : Figure 1.2 BRep_CurveRepresentation class diagram

二 边 Edge 边 (edge) 是对应于一维对象 - 曲线的拓朴实体 边可以是面的边界 ( 如长方体的 12 条边之一 ); 也可以只是一条不属于任何面的 悬空 边 (floating edge), 想像一下在构建锥形体或扫掠体之前的轮廓线 ; 面的边可以被两个或更多面共享, 或者只属于一个面 如下图所示 : Figure 2.1 Model used to illustrate Edge 在上图中用不同的颜色将不同类型的边区别开来 : 红色 : 不属于任何面的悬空边 (floating edge); 绿色 : 只属于一个面的自由边 (free edge); 黄色 : 属于两个或多个面的共享边 (shared edge); 边 Edge 包含几种几何表示 : 三维空间中的曲线 C(t), 由 Geom_Curve 实现 这是边的基本表示方式 ; 曲线 P(t) 为二维空间的参数曲线, 用来表示属于面的边, 通常被称为 pcurves, 由类 Geom2d_Curve 实现 ; 多段线 (Polygonal) 由一组三维点表示, 且由类 Poly_Polygon3D 实现 多段线 ( Polygonal ) 也可由一组三角剖分面上点的索引来表示, 且由类 Poly_PolygonOnTriangulation 实现 他们的表示都可以使用前面提到的类 BRep_Tool 来获取 例如 : Standard_Real afirst, alast, apfirst, aplast; Handle(Geom_Curve) acurve3d = BRep_Tool:: ::Curve (anedge, afirst, alast); Handle(Geom2d_Curve) apcurve aface, apfirst, aplast); = BRep_Tool:: ::CurveOnSurface (anedge, 边必须有曲面上的参数曲线 (pcurves), 除了平面以外 边中所有曲线必须一致, 即朝向 相同 这样边上的点可以使用任意表示方式计算得到, 如曲线 C(t), 可以用 [first, last]

区间上的 t 来计算 ; 也可根据 u 在区间 [first1, last1] 上取得曲面 S1(P1x(u), P1y(u)) 上的点 Pi, 这里 Pi 是曲面 Si 上的参数曲线 pcurve 上的一点 1. 边的标志位 Edge Flags 边中的标志位有三种 : static const Standard_Integer ParameterMask = 1; static const Standard_Integer RangeMask = 2; static const Standard_Integer DegeneratedMask = 4; 这里对前两种标志位进行说明 : RangeMask same range :(BRep_Tool::SameRange()) 取值区间相同, 即几何表示的曲线参数取值区间相同 ; ParameterMask same parameter :(BRep_Tool::SameParameter()) 参数相同, 即当 C (t)=s1(p1x(t), P1y(t)) 时, 对于同样的参数 t,c(t) 和曲面 S1 上的点 P1(t) 相同 即 任何边上的点都对应参数曲线上相同的参数值 许多算法假定设置了这两个标志位, 因此建议你注意这种情况, 一定要设置这些标志位 2. 边的容差 Tolerance 边的容差 (Tolerance) 是其三维曲线和其他任何表示方式之间的最大偏差 其几何意义 就是以容差为半径沿边的一个包含边的三维曲线及其他任何形式表示的管子 如下图所示 : Figure 2.2 Edge Tolerance 3. 特殊类型的边 Special edge types 在 OpenCascade 有两种特别类型的边, 他们是 :

缝合边 (seam edge): 即在同一个面上出现两次的边 ( 如 : 在同一个面上具有 2 个参数 曲线 ); 退化边 (degenerated edge): 这种边位于曲面的奇异点处, 在三维空间中退化为一个点 ; 球面中这两种类型的边都有 缝合边位于经线 (U iso-lines), 参数为 0 和 2*PI 退化边位于南北极点, 对应于纬线 (V iso-lines), 参数为 -PI/2 和 PI/2 因为球面的参数方程为 : 当参数 u 取 0 和 2*PI 时, 球面的参数方程计算如下 : S(0, v) = P + r cos( v) (cos(0) Dx + sin(0) Dy) + r sin( v) Dz S(2π, v) = P + r cos( v) (cos(2π ) Dx + sin(2π ) Dy) + r sin( v) Dz S(0, v) = P + r cos( v) Dx + r sin( v) Dz S(2π, v) = P + r cos( v) Dx + r sin( v) Dz 从计算结果可以看出, 缝合边是位于 Dx 和 Dz 所在平面上的半圆弧 当参数 v 取 -PI/2 和 PI/2 时, 球面的参数方程计算如下 : v [ π / 2, π / 2] v [ π / 2, π / 2] π π π S( u, ± ) = P + r cos( ± ) (cos( u) Dx + sin( u) Dy) + r sin( ± ) Dz 2 2 2 π S( u, ± ) = P ± r Dz 2 从计算结果可以看出, 曲面上的两个边分别退化为两个点 即 v 取 -PI/2 和 PI/2 时球面的两个退化边分别位于南北极点上 如下图所示 : Figure 2.3 seam edge and degenerated edge of sphere

另外例子 : 环形体 (torus) 圆柱体 (cylinder) 圆锥体 (cone) 环形体 torus 有两条 缝合边 (seam edge), 对应于它的参数空间的边界 ; 圆柱体 (cylinder) 有一条缝合边 (seam edge) 圆锥 (cone) 顶点处为退化边 (degenerated edge) 检查边是否是缝化边或退化边, 可以使用函数 BRep_Tool::IsClosed() 和 BRep_Tool::Degenerated() 4. 边的朝向 边的朝向为正向 (forward edge orientation) 意味着边的逻辑方向与曲线的方向相同 反 向 (reversed) 意味着逻辑方向与曲线方向相反 所以, 缝合边 (seam-edge) 在一个面中总 是有两个朝向 : 一个反向 (reversed), 一个正向 (forward)

三 示例程序 Example Code 以边界表示 BRep 表示的球面为例, 说明其边的类型 程序代码如下所示 : /* * Copyright (c) 2013 eryar All Rights Reserved. * * File : Main.cpp * Author : eryar@163.com * Date : 2013-08-24 16:11 * Version : 1.0v * * Description : Demonstrate seam edge and degenerated edge of sphere. * */ #include <iostream> // OpenCascade Library. #define WNT #include <TopoDS.hxx> #include <TopoDS_Edge.hxx> #include <TopExp_Explorer.hxx> #include <BRepPrimAPI_MakeSphere.hxx> #pragma comment(lib, "TKernel.lib") #pragma comment(lib, "TKMath.lib") #pragma comment(lib, "TKBRep.lib") #pragma comment(lib, "TKPrim.lib") #pragma comment(lib, "TKTopAlgo.lib") /** * @breif Dump orientation types. * Orientation definitaion: * enum TopAbs_Orientation * TopAbs_FORWARD, * TopAbs_REVERSED, * TopAbs_INTERNAL, * TopAbs_EXTERNAL * ; */ std::string dumporientation(const TopAbs_Orientation& orient) std::string strtype; switch (orient) case TopAbs_FORWARD: strtype = "TopAbs_FORWARD"; break; case TopAbs_REVERSED: strtype = "TopAbs_REVERSED"; break; case TopAbs_INTERNAL:

strtype = "TopAbs_INTERNAL"; break; case TopAbs_EXTERNAL: strtype = "TopAbs_EXTERNAL"; break; return strtype; /** * @breif Dump edge information. */ void processedge(const TopoDS_Edge& edge, const TopoDS_Face& face) Standard_Real dtolerance = BRep_Tool::Tolerance(edge); Standard_Boolean bisgeometric = BRep_Tool::IsGeometric(edge); Standard_Boolean bissameparameter = BRep_Tool::SameParameter(edge); Standard_Boolean bissamerange = BRep_Tool::SameRange(edge); Standard_Boolean bisdegenerated = BRep_Tool::Degenerated(edge); Standard_Boolean bisclosed = BRep_Tool::IsClosed(edge, face); TopAbs_Orientation norientation = edge.orientation(); // Dump edge info. std::cout<<"====== Edge Info ======="<<std::endl; std::cout<<"tolerance: "<<dtolerance<<std::endl; std::cout<<"orientation: "<<dumporientation(norientation)<<std::endl; std::cout<<"geometric: "<<(bisgeometric?"true":"false")<<std::endl; std::cout<<"same Parameter: "<<(bissameparameter?"true":"false")<<std::endl; std::cout<<"same Range: "<<(bissamerange? "True":"False")<<std::endl; std::cout<<"degenerated edge: "<<(bisdegenerated?"true":"false")<<std::endl; std::cout<<"seam edge: "<<(bisclosed? "True":"False")<<std::endl; // Dump vertex of the edge. for (TopExp_Explorer vertexitr(edge, TopAbs_VERTEX); vertexitr.more(); vertexitr.next()) const TopoDS_Vertex& avertex = TopoDS::Vertex(vertexItr.Current()); gp_pnt pnt = BRep_Tool::Pnt(aVertex); std::cout<<"vertex: ("<<pnt.x()<<", "<<pnt.y()<<", "<<pnt.z()<<")"<<std::endl; int main(void) Standard_Integer nspherefacecount = 0; Standard_Integer nsphereedgecount = 0; TopoDS_Shape sphere = BRepPrimAPI_MakeSphere(1.0); for (TopExp_Explorer faceitr(sphere, TopAbs_FACE); faceitr.more(); faceitr.next())

const TopoDS_Face& aface = TopoDS::Face(faceItr.Current()); ++nspherefacecount; for (TopExp_Explorer edgeitr(aface, TopAbs_EDGE); edgeitr.more(); edgeitr.next()) const TopoDS_Edge& aedge = TopoDS::Edge(edgeItr.Current()); processedge(aedge, aface); ++nsphereedgecount; std::cout<<"sphere face count: "<<nspherefacecount<<std::endl; std::cout<<"sphere edge count: "<<nsphereedgecount<<std::endl; return 0; 程序运行结果如下所示 : ====== Edge Info ======= Tolerance: 1e-007 Orientation: TopAbs_REVERSED Geometric: True Same Parameter: True Same Range: True Degenerated edge: True Seam edge: False Vertex: (6.12323e-017, -1.49976e-032, 1) Vertex: (6.12323e-017, -1.49976e-032, 1) ====== Edge Info ======= Tolerance: 1e-007 Orientation: TopAbs_FORWARD Geometric: True Same Parameter: True Same Range: True Degenerated edge: False Seam edge: True Vertex: (6.12323e-017, -1.49976e-032, 1) Vertex: (6.12323e-017, -1.49976e-032, -1) ====== Edge Info ======= Tolerance: 1e-007 Orientation: TopAbs_FORWARD Geometric: True Same Parameter: True Same Range: True Degenerated edge: True Seam edge: False

Vertex: (6.12323e-017, -1.49976e-032, -1) Vertex: (6.12323e-017, -1.49976e-032, -1) ====== Edge Info ======= Tolerance: 1e-007 Orientation: TopAbs_REVERSED Geometric: True Same Parameter: True Same Range: True Degenerated edge: False Seam edge: True Vertex: (6.12323e-017, -1.49976e-032, 1) Vertex: (6.12323e-017, -1.49976e-032, -1) Sphere face count: 1 Sphere edge count: 4 Press any key to continue... 从运行结果可以看, 当球的边为退化边时, 边的两个顶点的坐标值相同 退化边位于球 的南北极点上 缝合边为连接两个退化边的曲线 根据遍历顺序, 第一条边为退化边 (degenerated edge), 其朝向为反向 (reversed); 第二条边为缝合边 (seam-edge), 其朝向为正向 (forward); 第三条边为退化边, 其朝向为正向 (forward); 第四条边为缝合边, 其朝向为反向 (reversed) 由上可见, 缝合边有两个朝向, 一个正向一个反向

四 结论 Conclusion 对与几何相关的拓朴边 (edge) 的类的属性数据进行详细说明 并结合程序代码详细说明边的标志位 (myflags) 属性的意义, 从参数方程出发, 理解缝合边 (seam-edge) 和退化边 (degenerated edge), 即标志位中 DegeneratedMask 的意义 五 参考资料 1. Roman Lygin, OpenCascade notes, opencascade.blogspot.com 2. 孙家广等. 计算机图形学. 清华大学出版社 3. OpenCascade source code.