Topology and Geometry in OpenCascade-Face-卡核

Topology and Geometry in OpenCascade-Face

摘要Abstract：本文简要介绍了几何造型中的边界表示法（BRep），并结合程序说明OpenCascade中的边界表示的具体实现，即拓朴与几何的联系。对具有几何信息的拓朴结构顶点（vertex）、边（edge）、面（face）进行了详细说明。本文仅对面（Face）进行说明。

关键字Key Words：OpenCascade、BRep、Topology、Geometry、Face

一、引言 Introduction

面（Face）由一个外环和若干个内环来表示，内环完全在外环之内。根据环（Wire）的定义，在面上沿环的方向前进，左侧总在面内，右侧总在面外。面有方向性，一般用其外法矢方向作为该面的正向。若一个面的法矢向外，称为正向面；反之，称为反向面。面的形状（surface）由面的几何信息来表示，可以是平面或曲面，平面可用平面方程来描述，曲面可以用控制多边形或型值点来描述（NURBS曲面），也可用曲面方程（隐式、显式或参数形式）来描述。对于参数曲面，通常在其二维参数域上定义环，这样就可以由一些二维的有向边来表示环，集合运算中对面的分割也可在二维参数域上进行。

OpenCascade中的面结构如下图所示：

Figure 1.1 OpenCascade Face

根据其类图可知，面除了其参数曲面mySurface外，还包含显示曲面的剖分（由三角形组成）。当面在着色显示模式下时会计算出面的三角形。可视化的算法是由BRepMesh::Mesh()来为每个面三角剖分后用来显示。

二、面 Face

面（Face）是用来描述三维实体边界的拓朴实体。面是由底层的曲面及一个或多个环（Wire）来描述。例如，一个圆柱体包含三个面：底面、顶面和侧面。每个面都是无限的（无界的）（Geom_Plane和Geom_CylindricalSurface），通过边界来限定无限的面得到面。即用位于Geom_Circle上的边形成的环，限定出底面和顶面。侧面包含四条边：其中两条边与顶面和底面其享，剩下的两条边是缝合边（seam edge），参看之前的讨论。限定侧面的环包含缝合边两次，两个缝合边具有不同的朝向。

2.1 曲面 Surface

让我们简要回顾一下什么是曲面（surface）。如果你在高中时学过数学分析（mathematical analysis），那么你可能对这个概念已经熟记于心。如果没学过，那么可能需要读些文章来自学。在wikipedia上有关于参数曲面的简单例子。

曲面将二维参数空间{u, v}映射到三维空间。如下图所示：

Figure 2.1 Map parameter space{u, v} into 3D space

参数区间可以有界，也可以是无界的，也可以只在一个方向上有界。如平面（Plane）的参数空间是无界的；NURBS是有界的；圆柱面（Cylindrical Surface）在U方向上有界（U∈[0，2π]），在V方向上是无界的。

Geom_Surface::Value()返回一个对应于参数空间中的参数（U，V）的空间点（X，Y，Z）。例如：地球上的任意一点都由纬度（V）和经度（U）表示，但是在世界坐标系中可以看作三维点（假如地球的中心定义为原点）。这个函数是纯虚函数，所有派生自Geom_Surface的类都有对这个函数的实现，使计算各种曲面上对应参数的空间点的方式统一。

让我们回想一下，边必须具有三维曲线（3D Curve）和曲面空间中的参数曲线（pcurve）。而OpenCascade要求面的边界（环wire）在三维和二维空间中必须是闭合的。因此，圆柱的侧面是采用前面我们讨论的那样来描述的。

2.2 朝向 Orientation

面的朝向表示面的法向与曲面法向之间的关系（Face orientation shows how face normal is aligned with its surface normal）。若面的朝向是TopAbs_FORWARD（向前），则面的法向与曲面的法向一致；若面的朝向是TopAbs_REVERSED（反向），则面的法向与曲面的法向相反。面的法向表示材质的位置，材料位于面的背面。在正确描述的实体中，所有面的法向都是向外的，如下图所示：

Figure 2.2 In a correct solid body all face normals go outward

面上的材料是由边的朝向确定的。方向是由曲面不是面（）的法向和边的微分的叉积确定。若边的朝向向前，则边的导数等于它的三维曲线的导数；若边的朝向反向，则边的导数与它的三维曲线的导数相反。也许考虑边的参数曲线会更易于理解：假如边是向前的，材料在它的左侧，假如是反向的，则材料在它的右侧。这让我想起了格林公式中对平面区域的边界曲线正向的规定：对平面区域D的边界曲线L，我们规定L的正向如下：当观察者沿L的这个方向行走时，D内在他附近处的那一部分总在他的左边。

Figure 2.3 Orientation of the edge

当用积分区域的概念来理解边的朝向时好像要容易些。把参数空间看作积分区域D，其中蓝色的边表示材料在参数曲线的右侧；红色的边表示材料在边的左侧。

2.3 容差 Tolerance

面的容差的几何意义是包围面的一个具有厚度的板。如下图所示：

Figure 2.4 Face Tolerance

与边和顶点的容差相比，建模算法中用到面的容差的情况相对要少。通常都使用默认值Precision::Confusion()。通常情况下，OpenCascade要求注意以下条件：

当边位于面上，顶点位于边上时，面的容差<=边的容差<=顶点的容差。

2.4 三角剖分 Triangulation

除了面的参数表示，为了面的显示，需要对面进行三角剖分，剖分得到的三角形也保存在面的数据结构中。当在着色渲染模式下时，可视化算法内部调用了BRepMesh::Mesh()来对每个面进行三角剖分，并将三角剖分得到的三角形加到面的数据结构中。

2.5 附加位置 Additional location

与边和顶点不同，面有附加位置信息（TopLoc_Location），它是面（BRep_TFace）的成员变量。所以，在使用底层曲面或三角剖分得到的三角形时，不要忘了将其考虑进去。

三、示例程序 Example Code

3.1 底层创建面并访问其数据 Creating a face bottom-up and accessing the data

与创建边和顶点一样，需要类BRep_Builder和Brep_Tool来从底层创建面和访问面中的数据。代码如下所示：

BRep_Builder aBuilder;
TopoDS_Face aFace;
aBuilder.MakeFace (aFace, aSurface, Precision::Confusion());
...
TopLoc_Location aLocation;
Handle(Geom_Surface) aSurf = BRep_Tool::Surface (aFace, aLocation);
gp_Pnt aPnt = aSurf->Value (aU, aV).Transformed (aLocation.Transformation());
//or to transform a surface at once
//Handle(Geom_Surface) aSurf = BRep_Tool::Surface (aFace);
//gp_Pnt aPnt = aSurf->Value (aU, aV);
Handle(Poly_Triangulation) aTri = BRep_Tool::Triangulation (aFace, aLocation);
aPnt = aTri->Nodes.Value (i).Transformed (aLocation.Transformation());

一定要考虑面的附加位置信息。