思享工具箱

Topology and Geometry in OpenCascade-Vertex

摘要：
关键词：OpenCascade，BRep，Topology，Geometry 1.简介边界表示法也称为BRep表示法，是几何建模中最成熟、最明确的表示法。STEPISO-10303-42的相关资源：http://www.steptools.com/support/stdev_docs/express/step_irs/index.htmlFigure2.1TopologydatastructureinOpenCascadeTopoDS_Shape它由值控制，包含三个成员变量：myLocation、myOrient和myTShape。如下图所示：图2.3 OpenCascade中的拓扑数据结构从上面的类图中可以看出，只有三个拓扑对象具有几何表示数据：顶点、边和面，分别是BRep_TVertex、BRep_TEdge、BRep_face图3.1 VertexTolerance与其他几何库的全局精度不同。OpenCascade将公差作为局部属性。

Topology and Geometry in OpenCascade-Vertex

eryar@163.com

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

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

一、引言 Introduction

边界表示（Boundary Representation）也称为BRep表示，它是几何造型中最成熟、无二义的表示法。实体的边界通常是由面的并集来表示，而每个面又由它所在的曲面的定义加上其边界来表示，面的边界是边的并集，而边又是由点来表示的。 

边界表示的一个重要特征是描述形体的信息包括几何信息（Geometry）和拓朴信息（Topology）两个方面。拓朴信息描述形体上的顶点、边、面的连接关系，它形成物体边界表示的“骨架”。形体的几何信息犹如附着在“骨架”上的肌肉。例如，形体的某个面位于某一个曲面上，定义这一曲面方程的数据就是几何信息。此外，边的形状、顶点在三维空间中的位置（点的坐标）等都是几何信息，一般来说，几何信息描述形体的大小、尺寸、位置和形状等。 

在边界表示法中，边界表示就按照体－面－环－边－点的层次，详细记录构成形体的所有几何元素的几何信息及其相互连接的拓朴关系。这样，在进行各种运算和操作中，就可以直接取得这些信息。 

下图所示为由一条边连接的两个面组成的壳（shell）： 

wps_clip_image-13951

Figure 1.1 Structure of a shell formed from two faces 

上图所示的形状表示为TS， 面TF1和TF2，有七条边TE1~TE7和六个顶点TV1~TV6。 

环TW1引用边TE1~TE4；环TW2引用TE4~TE7 。边引用的顶点如下：TE1（TV1，TV4），TE2（TV1，TV2），TE3（TV2，TV3），TE4（TV3，TV4），TE5（TV4，TV5），TE6（TV5，TV6），TE7（TV3，TV6）。 

wps_clip_image-26435

Figure 1.2 Data structure of the shell formed from two faces connected at an edge 

注：OpenCascade中的这个数据结构中不包含“回溯引用（back references）”，即所有的引用只从复杂形状到简单形状。（Note that this data structure does not contain any “back references”. All references go from more comples underlying shapes to less complex ones.）有点有向图的意思。 

二、OpenCascade中的边界表示 BRep in OpenCascade

2.1 拓朴结构 TopoDS_Shape data structure

OpenCascade中的拓朴（topology）是根据STEP标准ISO-10303-42设计的。也许读一下这个标准中的有关概念还是很有帮助的。STEP ISO-10303-42的相关资源： 

http://www.steptools.com/support/stdev_docs/express/step_irs/index.html

wps_clip_image-23948

Figure 2.1 Topology data structure in OpenCascade 

wps_clip_image-14563

TopoDS_Shape由值控制，包含三个成员变量：myLocation、myOrient、myTShape。 

wps_clip_image-9340

Figure 2.2 TopoDS_Shape member fields 

2.2 拓朴与几何的联系 Connection with Geometry

现在我们来考虑一下拓朴结构与几何的关系。通过继承TopoDS包中的抽象的拓朴类实现了边界表示模型。如下图所示： 

wps_clip_image-67

Figure 2.3 Topology data structure in OpenCascade 

从上面的类图可以看出只有三种拓朴对象有几何表示数据：顶点（vertex）、边（edge）、面（face），分别为BRep_TVertex、BRep_TEdge、BRep_TFace。 

wps_clip_image-2791

Figure 2.4 TopoDS_TShape class diagram 

三、顶点 Vertex

顶点（vertex）的位置用几何点（gp_Pnt）来表示。点是几何造型中的最基本元素，自由曲线、曲面或其他形体均可用有序的点集表示。用计算机存储、管理、输出形体的实质就是对点集及其连接关系的处理。在正则形体定义中，不允许孤立点存在。 

顶点的另一个重要属性是容差（Tolerance），用来表示位置精度。顶点容差T的几何意义为以顶点为圆心半径为T的球。这个球必须包含所有与这个顶点相连的边的曲线的端点。 

wps_clip_image-25604

Figure 3.1 Vertex Tolerance 

与其他几何库有全局精度（global precision）不同，OpenCascade把容差作为局部属性（local properties）。由图2.4可知，容差是顶点、边、面的属性。这种方法有助于用更一般的方式来描述高精度的模型。如下图所示： 

wps_clip_image-13508

Figure 3.2 Vertex with different tolerance 

如果从底层来创建形状，最好的方法就是指定最小的允许误差。默认值Precision::Confusion()为1e-07。 

下面讨论顶点的朝向（orientation）属性。它没有直接的几何意义，但是根据约定，若顶点的朝向属性值为TopAbs_FORWARD，它就必须与表示边的曲线的参数值小的端部匹配。相应地，TopAbs_REVERSED的顶点与参数值大的端部匹配。例如，有条边位于圆弧上，圆弧半径为1且在Z＝0的平面上，起点为（1，0，0），向－Z轴向，曲线为逆时针方向。所以顶点（1，0，0）的朝向为TopAbs_FORWARD，顶点（0，1，0）的朝向为TopAbs_REVERSED。如下图所示： 

wps_clip_image-6315

Figure 3.3 Vertex Orientation attribute 

实现上图所示的程序代码如下所示： 

/*
*    Copyright (c) 2013 eryar All Rights Reserved.
*
*        File    : Main.cpp
*        Author  : eryar@163.com
*        Date    : 2013-08-17 21:46
*        Version : 1.0v
*
*    Description : Demonstrate how to build a edge bottom-up by BRepBuilderAPI_MakeEdge,
*                  and how to access vertex infomation by BRep_Tool.
*                  
*/

// OpenCascade library.
#define WNT
#include <gp_Circ.hxx>
#include <GC_MakeArcOfCircle.hxx>

#include <TopoDS_Edge.hxx>
#include <BRep_Tool.hxx>
#include <BRepBuilderAPI_MakeEdge.hxx>
#include <BRepBuilderAPI_MakeVertex.hxx>

#pragma comment(lib, "TKernel.lib")
#pragma comment(lib, "TKMath.lib")
#pragma comment(lib, "TKBRep.lib")
#pragma comment(lib, "TKGeomBase.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 attributes of the vertex.
*/
void DumpVertex(const TopoDS_Vertex& v)
{
    gp_Pnt p = BRep_Tool::Pnt(v);
    Standard_Real dTolerance = BRep_Tool::Tolerance(v);

    std::cout<<"Vertex position: ("<<p.X()<<", "<<p.Y()<<", "<<p.Z()<<")"<<std::endl;
    std::cout<<"Vertex Tolerance: "<<dTolerance<<std::endl;
    std::cout<<"Vertex orientation: "<<DumpOrientation(v.Orientation())<<std::endl;
    std::cout<<std::endl;
}

int main(int argc, char* argv[])
{
    gp_Circ circle;
    TopoDS_Edge edge;
    TopoDS_Vertex vertex1;
    TopoDS_Vertex vertex2;
    BRepBuilderAPI_MakeEdge edgeBuilder;

    circle.SetRadius(1.0);
    circle.SetAxis(gp::OZ());

    edgeBuilder.Init(GC_MakeArcOfCircle(circle, 0.0, M_PI/2.0, Standard_True));

    // Test single vertex.
    /*vertex1 = BRepBuilderAPI_MakeVertex(gp_Pnt(100.0, 200.0, 300.0));
    vertex2 = BRepBuilderAPI_MakeVertex(gp_Pnt(500.0, 600.0, 700.0));

    std::cout<<"Single vetex test: "<<std::endl;

    std::cout<<"Vertex 1 attributes: "<<std::endl;
    DumpVertex(vertex1);

    std::cout<<"Vertex 2 attributes: "<<std::endl;
    DumpVertex(vertex2);*/

    edge = edgeBuilder.Edge();
    vertex1 = edgeBuilder.Vertex1();
    vertex2 = edgeBuilder.Vertex2();

    std::cout<<"Test vertex belong to edge:"<<std::endl;

    std::cout<<"Vertex 1 attributes: "<<std::endl;
    DumpVertex(vertex1);

    std::cout<<"Vertex 2 attributes: "<<std::endl;
    DumpVertex(vertex2);

    return 0;
}

程序运行结果如下图所示： 

wps_clip_image-27439

Figure 3.4 Code example result 

BRep_Builder是从底层创建拓朴结构的类。如下代码所示为从底层创建顶点的示例：

gp_Pnt aPoint(100.0, 200.0, 300.0)
BRep_Builder aBuilder;
TopoDS_Vertex aVertex;

aBuilder.MakeVertex(aVertex, aPoint, Precision::Confusion());
aVertex.Orientation(TopAbs_REVERSED);

有一个方便的类也可用来创建顶点BRepBuilderAPI_MakeVertex，其内部也是使用了类BRep_Builder。所以，若想从底层创建拓朴结构，必须要熟悉BRep_Builder。 

BRep_Tool是用来访问拓朴结构中几何信息的工具，他的大部分的函数是静态的。如下代码所示为获取顶点的容差和几何点的方法：

Standard_Real aTolerance = BRep_Tool::Tolerance(aVertex);
gp_Pnt aPoint = BRep_Tool::Pnt(aVertex);

四、BRep文件中Vertex的数据

结合《BRep Format Description White Paper》中对<vertex data>的描述，及程序代码中对顶点数据的读取，分析OpenCascade的BRep表示中的顶点。 

wps_clip_image-14063

Figure 4.1 NBF-like definition of Vertex 

详细说明： 

<vertex data representation u parameter>u的使用方法说明如下： 

<vertex data representation data 1> 和参数u定义了三维曲线C上的点V的位置。参数u是曲线C上点V对应的参数：C（u）＝V。 

<vertex data representation data 2>和参数u定义了曲面上的二维曲线C上点V的位置。参数u是曲线C上点V对应的参数：C（u）＝V。 

<vertex data representation data 3>和参数u及<vertex data representation v parameter>v定义了曲面S上的点V：S（u，v）＝V。 

<vertex data tolerance>t定义如下所示： 

wps_clip_image-4091

读取Vertex部分的程序代码摘抄如下：

case TopAbs_VERTEX :
    {
      TopoDS_Vertex& V = TopoDS::Vertex(S);
      
      // Read the point geometry
      IS >> tol;
      IS >> X >> Y >> Z;
      myBuilder.MakeVertex(V,gp_Pnt(X,Y,Z),tol);
      Handle(BRep_TVertex) TV = Handle(BRep_TVertex)::DownCast(V.TShape());

      BRep_ListOfPointRepresentation& lpr = TV->ChangePoints();
      TopLoc_Location L;

      do {
        IS >> p1 >> val;
        
        Handle(BRep_PointRepresentation) PR;
        switch (val) {

        case 1 :
          {
            IS >> c;

//  Modified by Sergey KHROMOV - Wed Apr 24 13:59:09 2002 Begin
        if (myCurves.Curve(c).IsNull())
          break;
//  Modified by Sergey KHROMOV - Wed Apr 24 13:59:13 2002 End

            Handle(BRep_PointOnCurve) POC =
              new BRep_PointOnCurve(p1,
                                    myCurves.Curve(c),
                                    L);
            PR = POC;
          }
          break;

        case 2 :
          {
            IS >> pc >> s;

//  Modified by Sergey KHROMOV - Wed Apr 24 13:59:09 2002 Begin
        if (myCurves2d.Curve2d(pc).IsNull() ||
        mySurfaces.Surface(s).IsNull())
          break;
//  Modified by Sergey KHROMOV - Wed Apr 24 13:59:13 2002 End

            Handle(BRep_PointOnCurveOnSurface) POC =
              new BRep_PointOnCurveOnSurface(p1,
                                             myCurves2d.Curve2d(pc),
                                             mySurfaces.Surface(s),
                                             L);
            PR = POC;
          }
          break;

        case 3 :
          {
            IS >> p2 >> s;

//  Modified by Sergey KHROMOV - Wed Apr 24 13:59:09 2002 Begin
        if (mySurfaces.Surface(s).IsNull())
          break;
//  Modified by Sergey KHROMOV - Wed Apr 24 13:59:13 2002 End

            Handle(BRep_PointOnSurface) POC =
              new BRep_PointOnSurface(p1,p2,
                                      mySurfaces.Surface(s),
                                      L);
            PR = POC;
          }
          break;
        }

        if (val > 0) {
          IS >> l;
      if (!PR.IsNull()) {
        PR->Location(Locations().Location(l));
        lpr.Append(PR);
      }
        }
      } while (val > 0);
    }
    break;

从BRep文件中可知，大部分的Vertex只有如下数据：

0101101
*
Ve
1e-007
2 3 0
0 0

即只使用了BRep_Builder.MakeVertex()创建创建顶点（vertex），还可记录顶点的类型： 

1. 若顶点在三维空间中的曲线上Geom_Curve，则记录三维曲线的索引号及参数u； 

2. 若顶点在二维空间中的曲线上Geom2d_Curve，则记录二维曲线的索引号及参数u； 

3. 若顶点在曲面上Geom_Surface，则记录曲面的索引号及参数（u，v）； 

五、结论 Conclusion

结合博客“OpenCascade notes”及《BRep format Description white paper》对OpenCascade的拓朴结构中的顶点（vertex）的属性数据进行说明。结合程序说明了顶点的容差及朝向的意义及从底层创建顶点的方法。通过BRep_Tool的静态函数可以获取顶点的几何数据及其他属性。 

发现在BRep文件中还对顶点进行了分类：三维曲线上的点、二维曲线上的点和曲面上的点。 

六、参考资料 

1. OpenCascade notes: opencascade.blogspot.com 

2. 孙家广等. 计算机图形学. 清华大学出版社 

PDF Version: Topology and Geometry in OpenCascade-Vertex

标签：#opencascade#顶点坐标#拓朴 | 浏览：125 | 发布日期：

免责声明：文章转载自《Topology and Geometry in OpenCascade-Vertex》仅用于学习参考。如对内容有疑问，请及时联系本站处理。

上篇JS阻止默认行为C#堆和栈下篇

宿迁高防，2C2G15M，22元/月；香港BGP，2C5G5M，25元/月 雨云优惠码：MjYwNzM=

相关文章

OpenCASCADE动画功能2

OpenCASCADE动画功能2eryar@163.com OpenCASCADE是一个开发平台，主要提供三维曲面和实体建模、CAD数据交换及可视化等功能。OCCT最适用于开发三维建模CAD软件、加工制造或测量（CAM）软件及数值仿真软件。对于一些加工制造软件，需要简单的动画仿真功能。本文主要就来说说OCCT的动画功能。 OCCT从7.1.0版本引入了类A...

[原] OpenGL ES 学习笔记 (一)

1.OpenGL ES 的坐标系在屏幕上的分布 OpenGL ES 的坐标系{x, y, z} 通过图片的三维坐标系可以知道： - 它是一个三维坐标系 {x, y, z} - 三维坐标中心在正方体的几何中心 {0， 0， 0} - 整个坐标系是 [0, 1] 的点，也就是说 OpenGL 中只支持 0 ~ 1 的点 （这里所讲的 0 和 1 ，最好理...

使用 mesh 实现多边形裁剪图片！Cocos Creator！

和 mask 裁剪图片说拜拜，用上高性能的 mesh + shader 。文章底部获取完整代码！ 效果预览： 使用方法： 创建一个空节点 添加用户脚本组件 mesh-texture-mask 添加图片 添加修改多边形顶点坐标 实现原理 创建 mesh mesh 是什么？ mesh 是决定一个物体形状的东西。 例如在二维中可以是正方形、圆形、三角...

OpenCASCADE 参数曲面面积

OpenCASCADE 参数曲面面积 eryar@163.com Abstract. 本文介绍了参数曲面的第一基本公式，并应用曲面的第一基本公式，结合OpenCASCADE中计算多重积分的类，对任意参数曲面的面积进行计算。 Key Words. Parametric Curve, Parametric Surface, Gauss Integration,...

OpenGL实现多层绘制（Layered Rendering） [转]

http://blog.csdn.net/u010462297/article/details/50589991 引言 在某些情况下会需要用到多层绘制。FBO下有多个颜色挂接点（Color Attachment），可以用不同的挂接点挂接不同的纹理对象，实现绘制多张纹理（MRT），这在之前的文章里已经有所描述。但是有时候这种方法是不够好用的： - 当纹理非...

Direct3D11学习：（九）绘制基本几何体

转载请注明出处：http://www.cnblogs.com/Ray1024 一、概述 Direct3D中很多复杂的几何效果都是由基本的几何体组合而成的，这篇文章中，我们来学习集中常见的基本几何体的绘制方法。 二、准备工作 我们使用一个类来组织这些绘制基本几何体的代码，以方便我们以后的使用。GeometryGenerator是一个工具类，用于生成诸如网格、...

最新文章

随机推荐

思享工具箱导航