Topology and Geometry in OpenCascade

Size: px

Start display at page:

Download "Topology and Geometry in OpenCascade"

复侍郎
5 years ago
Views:

1 Topology and Geometry in OpenCascade-Edge 摘要 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.

2 二边 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]

3 区间上的 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), 参数为

4 缝合边 (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 π S( u, ± ) = P ± r Dz 2 从计算结果可以看出, 曲面上的两个边分别退化为两个点即 v 取 -PI/2 和 PI/2 时球面的两个退化边分别位于南北极点上如下图所示 : Figure 2.3 seam edge and degenerated edge of sphere

5 另外例子 : 环形体 (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)

6 三示例程序 Example Code 以边界表示 BRep 表示的球面为例, 说明其边的类型程序代码如下所示 : /* * Copyright (c) 2013 eryar All Rights Reserved. * * File : Main.cpp * Author : eryar@163.com * Date : :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") /** 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:

7 strtype = "TopAbs_INTERNAL"; break; case TopAbs_EXTERNAL: strtype = "TopAbs_EXTERNAL"; break; return strtype; /** 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())

8 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: ( e-017, e-032, 1) Vertex: ( e-017, e-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: ( e-017, e-032, 1) Vertex: ( e-017, e-032, -1) ====== Edge Info ======= Tolerance: 1e-007 Orientation: TopAbs_FORWARD Geometric: True Same Parameter: True Same Range: True Degenerated edge: True Seam edge: False

9 Vertex: ( e-017, e-032, -1) Vertex: ( e-017, e-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: ( e-017, e-032, 1) Vertex: ( e-017, e-032, -1) Sphere face count: 1 Sphere edge count: 4 Press any key to continue... 从运行结果可以看, 当球的边为退化边时, 边的两个顶点的坐标值相同退化边位于球的南北极点上缝合边为连接两个退化边的曲线根据遍历顺序, 第一条边为退化边 (degenerated edge), 其朝向为反向 (reversed); 第二条边为缝合边 (seam-edge), 其朝向为正向 (forward); 第三条边为退化边, 其朝向为正向 (forward); 第四条边为缝合边, 其朝向为反向 (reversed) 由上可见, 缝合边有两个朝向, 一个正向一个反向

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

Topology and Geometry in OpenCascade

Topology and Geometry in OpenCascade Location and Orientaion eryar@163.com 摘要 Abstract: 本文简要介绍了几何造型中的边界表示法 (BRep), 并结合程序说明 OpenCascade 中的边界表示的具体实现, 即拓朴与几何的联系拓朴结构中的位置 (Location) 和朝向 (Orientation) 进行了详细说明