OpenCascade 边界表示法(BRep)
转自 http://www.cppblog.com/eryar/archive/2013/08/20/202678.html
Topology and Geometry in OpenCascade-Vertex
摘要Abstract:本文简要介绍了几何造型中的边界表示法(BRep),并结合程序说明OpenCascade中的边界表示的具体实现,即拓朴与几何的联系。对具有几何信息的拓朴结构顶点(vertex)、边(edge)、面(face)进行了详细说明。本文只对顶点数据进行说明。
关键字Key Words:OpenCascade、BRep、Topology、Geometry
一、引言 Introduction
边界表示(Boundary Representation)也称为BRep表示,它是几何造型中最成熟、无二义的表示法。实体的边界通常是由面的并集来表示,而每个面又由它所在的曲面的定义加上其边界来表示,面的边界是边的并集,而边又是由点来表示的。
边界表示的一个重要特征是描述形体的信息包括几何信息(Geometry)和拓朴信息(Topology)两个方面。拓朴信息描述形体上的顶点、边、面的连接关系,它形成物体边界表示的“骨架”。形体的几何信息犹如附着在“骨架”上的肌肉。例如,形体的某个面位于某一个曲面上,定义这一曲面方程的数据就是几何信息。此外,边的形状、顶点在三维空间中的位置(点的坐标)等都是几何信息,一般来说,几何信息描述形体的大小、尺寸、位置和形状等。
在边界表示法中,边界表示就按照体-面-环-边-点的层次,详细记录构成形体的所有几何元素的几何信息及其相互连接的拓朴关系。这样,在进行各种运算和操作中,就可以直接取得这些信息。
下图所示为由一条边连接的两个面组成的壳(shell):
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)。
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
Figure 2.1 Topology data structure in OpenCascade
TopoDS_Shape由值控制,包含三个成员变量:myLocation、myOrient、myTShape。
Figure 2.2 TopoDS_Shape member fields
2.2 拓朴与几何的联系 Connection with Geometry
现在我们来考虑一下拓朴结构与几何的关系。通过继承TopoDS包中的抽象的拓朴类实现了边界表示模型。如下图所示:
Figure 2.3 Topology data structure in OpenCascade
从上面的类图可以看出只有三种拓朴对象有几何表示数据:顶点(vertex)、边(edge)、面(face),分别为BRep_TVertex、BRep_TEdge、BRep_TFace。
Figure 2.4 TopoDS_TShape class diagram
三、顶点 Vertex
顶点(vertex)的位置用几何点(gp_Pnt)来表示。点是几何造型中的最基本元素,自由曲线、曲面或其他形体均可用有序的点集表示。用计算机存储、管理、输出形体的实质就是对点集及其连接关系的处理。在正则形体定义中,不允许孤立点存在。
顶点的另一个重要属性是容差(Tolerance),用来表示位置精度。顶点容差T的几何意义为以顶点为圆心半径为T的球。这个球必须包含所有与这个顶点相连的边的曲线的端点。
Figure 3.1 Vertex Tolerance
与其他几何库有全局精度(global precision)不同,OpenCascade把容差作为局部属性(local properties)。由图2.4可知,容差是顶点、边、面的属性。这种方法有助于用更一般的方式来描述高精度的模型。如下图所示:
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。如下图所示:
Figure 3.3 Vertex Orientation attribute
实现上图所示的程序代码如下所示:
2 * Copyright (c) 2013 eryar All Rights Reserved.
3 *
4 * File : Main.cpp
5 * Author : eryar@163.com
6 * Date : 2013-08-17 21:46
7 * Version : 1.0v
8 *
9 * Description : Demonstrate how to build a edge bottom-up by BRepBuilderAPI_MakeEdge,
10 * and how to access vertex infomation by BRep_Tool.
11 *
12 */
13
14 // OpenCascade library.
15 #define WNT
16 #include <gp_Circ.hxx>
17 #include <GC_MakeArcOfCircle.hxx>
18
19 #include <TopoDS_Edge.hxx>
20 #include <BRep_Tool.hxx>
21 #include <BRepBuilderAPI_MakeEdge.hxx>
22 #include <BRepBuilderAPI_MakeVertex.hxx>
23
24 #pragma comment(lib, "TKernel.lib")
25 #pragma comment(lib, "TKMath.lib")
26 #pragma comment(lib, "TKBRep.lib")
27 #pragma comment(lib, "TKGeomBase.lib")
28 #pragma comment(lib, "TKTopAlgo.lib")
29
30 /**
31 * @breif Dump orientation types.
32 * Orientation definitaion:
33 * enum TopAbs_Orientation {
34
35 * TopAbs_FORWARD,
36
37 * TopAbs_REVERSED,
38
39 * TopAbs_INTERNAL,
40
41 * TopAbs_EXTERNAL
42
43 * };
44 */
45 std::string DumpOrientation(const TopAbs_Orientation& orient)
46 {
47 std::string strType;
48
49 switch (orient)
50 {
51 case TopAbs_FORWARD:
52 strType = "TopAbs_FORWARD";
53 break;
54
55 case TopAbs_REVERSED:
56 strType = "TopAbs_REVERSED";
57 break; 58 59 case TopAbs_INTERNAL: 60 strType = "TopAbs_INTERNAL"; 61 break; 62 63 case TopAbs_EXTERNAL: 64 strType = "TopAbs_EXTERNAL"; 65 break; 66 } 67 68 return strType; 69 } 70 71 /** 72 * @breif Dump attributes of the vertex. 73 */ 74 void DumpVertex(const TopoDS_Vertex& v) 75 { 76 gp_Pnt p = BRep_Tool::Pnt(v); 77 Standard_Real dTolerance = BRep_Tool::Tolerance(v); 78 79 std::cout<<"Vertex position: ("<<p.X()<<", "<<p.Y()<<", "<<p.Z()<<")"<<std::endl; 80 std::cout<<"Vertex Tolerance: "<<dTolerance<<std::endl; 81 std::cout<<"Vertex orientation: "<<DumpOrientation(v.Orientation())<<std::endl; 82 std::cout<<std::endl; 83 } 84 85 int main(int argc, char* argv[]) 86 { 87 gp_Circ circle; 88 TopoDS_Edge edge; 89 TopoDS_Vertex vertex1; 90 TopoDS_Vertex vertex2; 91 BRepBuilderAPI_MakeEdge edgeBuilder; 92 93 circle.SetRadius(1.0); 94 circle.SetAxis(gp::OZ()); 95 96 edgeBuilder.Init(GC_MakeArcOfCircle(circle, 0.0, M_PI/2.0, Standard_True)); 97 98 // Test single vertex. 99 /*vertex1 = BRepBuilderAPI_MakeVertex(gp_Pnt(100.0, 200.0, 300.0));100 vertex2 = BRepBuilderAPI_MakeVertex(gp_Pnt(500.0, 600.0, 700.0));101 102 std::cout<<"Single vetex test: "<<std::endl;103 104 std::cout<<"Vertex 1 attributes: "<<std::endl;105 DumpVertex(vertex1);106 107 std::cout<<"Vertex 2 attributes: "<<std::endl;108 DumpVertex(vertex2);*/109 110 edge = edgeBuilder.Edge();111 vertex1 = edgeBuilder.Vertex1();112 vertex2 = edgeBuilder.Vertex2();113 114 std::cout<<"Test vertex belong to edge:"<<std::endl;115 116 std::cout<<"Vertex 1 attributes: "<<std::endl;117 DumpVertex(vertex1);118 119 std::cout<<"Vertex 2 attributes: "<<std::endl;120 DumpVertex(vertex2);121 122 return 0;123 }
程序运行结果如下图所示:
Figure 3.4 Code example result
BRep_Builder是从底层创建拓朴结构的类。如下代码所示为从底层创建顶点的示例:
2
3 BRep_Builder aBuilder;
4
5 TopoDS_Vertex aVertex;
6
7 aBuilder.MakeVertex(aVertex, aPoint, Precision::Confusion());
8
9 aVertex.Orientation(TopAbs_REVERSED);
有一个方便的类也可用来创建顶点BRepBuilderAPI_MakeVertex,其内部也是使用了类BRep_Builder。所以,若想从底层创建拓朴结构,必须要熟悉BRep_Builder。
BRep_Tool是用来访问拓朴结构中几何信息的工具,他的大部分的函数是静态的。如下代码所示为获取顶点的容差和几何点的方法:
2
3 gp_Pnt aPoint = BRep_Tool::Pnt(aVertex);
4
四、BRep文件中Vertex的数据
结合《BRep Format Description White Paper》中对<vertex data>的描述,及程序代码中对顶点数据的读取,分析OpenCascade的BRep表示中的顶点。
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定义如下所示:
读取Vertex部分的程序代码摘抄如下:
2 // vertex
3 //---------
4
5 case TopAbs_VERTEX :
6 {
7 TopoDS_Vertex& V = TopoDS::Vertex(S);
8
9 // Read the point geometry
10 IS >> tol;
11 IS >> X >> Y >> Z;
12 myBuilder.MakeVertex(V,gp_Pnt(X,Y,Z),tol);
13 Handle(BRep_TVertex) TV = Handle(BRep_TVertex)::DownCast(V.TShape());
14
15 BRep_ListOfPointRepresentation& lpr = TV->ChangePoints();
16 TopLoc_Location L;
17
18 do {
19 IS >> p1 >> val;
20
21 Handle(BRep_PointRepresentation) PR;
22 switch (val) {
23
24 case 1 :
25 {
26 IS >> c;
27
28 // Modified by Sergey KHROMOV - Wed Apr 24 13:59:09 2002 Begin
29 if (myCurves.Curve(c).IsNull())
30 break;
31 // Modified by Sergey KHROMOV - Wed Apr 24 13:59:13 2002 End
32
33 Handle(BRep_PointOnCurve) POC =
34 new BRep_PointOnCurve(p1,
35 myCurves.Curve(c),
36 L);
37 PR = POC;
38 }
39 break;
40
41 case 2 :
42 {
43 IS >> pc >> s;
44
45 // Modified by Sergey KHROMOV - Wed Apr 24 13:59:09 2002 Begin
46 if (myCurves2d.Curve2d(pc).IsNull() ||
47 mySurfaces.Surface(s).IsNull())
48 break;
49 // Modified by Sergey KHROMOV - Wed Apr 24 13:59:13 2002 End
50
51 Handle(BRep_PointOnCurveOnSurface) POC =
52 new BRep_PointOnCurveOnSurface(p1,
53 myCurves2d.Curve2d(pc),
54 mySurfaces.Surface(s),55 L);56 PR = POC;57 }58 break;59 60 case 3 :61 {62 IS >> p2 >> s;63 64 // Modified by Sergey KHROMOV - Wed Apr 24 13:59:09 2002 Begin65 if (mySurfaces.Surface(s).IsNull())66 break;67 // Modified by Sergey KHROMOV - Wed Apr 24 13:59:13 2002 End68 69 Handle(BRep_PointOnSurface) POC =70 new BRep_PointOnSurface(p1,p2,71 mySurfaces.Surface(s),72 L);73 PR = POC;74 }75 break;76 }77 78 if (val > 0) {79 IS >> l;80 if (!PR.IsNull()) {81 PR->Location(Locations().Location(l));82 lpr.Append(PR);83 }84 }85 } while (val > 0);86 }87 break;
从BRep文件中可知,大部分的Vertex只有如下数据:
2 *
3 Ve
4 1e-007
5 2 3 0
6 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. 孙家广等. 计算机图形学. 清华大学出版社
OpenCascade 边界表示法(BRep)的更多相关文章
- OpenCASCADE PCurve of Topological Face
OpenCASCADE PCurve of Topological Face eryar@163.com Abstract. OpenCASCADE provides a class BRepBuil ...
- OpenCASCADE BRepTools
OpenCASCADE BRepTools eryar@163.com Abstract. OpenCASCADE BRepTools provides utilities for BRep data ...
- OpenCASCADE Make Primitives-Box
OpenCASCADE Make Primitives-Box eryar@163.com Abstract. By making a simple box to demonstrate the BR ...
- Topology and Geometry in OpenCascade-Face
Topology and Geometry in OpenCascade-Face eryar@163.com 摘要Abstract:本文简要介绍了几何造型中的边界表示法(BRep),并结合程序说明O ...
- OpenCASCADE BRep vs. OpenNURBS BRep
OpenCASCADE BRep vs. OpenNURBS BRep eryar@163.com Abstract. BRep short for Boundary Representation. ...
- OpenCascade Primitives BRep - Sphere
OpenCascade Primitives BRep - Sphere eryar@163.com Abstract. BRep is short for Boundary Representati ...
- OpenCascade Primitives BRep - Box
OpenCascade Primitives BRep - Box eryar@163.com Abstract. BRep is short for Boundary Representation. ...
- Topology Shapes of OpenCascade BRep
Topology Shapes of OpenCascade BRep eryar@163.com 摘要Abstract:通过对OpenCascade中的BRep数据的读写,理解边界表示法的概念及实现 ...
- Geometry Surface of OpenCascade BRep
Geometry Surface of OpenCascade BRep eryar@163.com 摘要Abstract:几何曲面是参数表示的曲面 ,在边界表示中其数据存在于BRep_TFace中, ...
随机推荐
- postgreSQL中跨库查询在windows下的实现方法
以下是在postgreSQL 8.1版本中的实践,其他版本类似: 1.将C:\Program Files\PostgreSQL\8.1\share\contrib下的dblink.sql复制到C:\P ...
- Deutsch lernen (10)
Dieser Weg Dieser Weg wird kein leichter sein. Dieser Weg wird steinig und schwer. Nicht mit vielen ...
- svn命令行批量删除和批量添加
svn命令行批量删除和批量添加 如果使用svn的命令行,例如在linux下的终端中使用,svn的添加命令是svn add,删除命令是svn del,但是缺乏批量的操作,如果我在资源管理器中,手动添加了 ...
- OpenCV:Adaboost训练时数据扩增
更准确的模型需要更多的数据,对于传统非神经网络机器学习方法,不同的特征需要有各自相符合的数据扩增方法. 1. 在使用opencv_traincascade.exe 过程中,图像读取在 classi ...
- centos6.5 安装Python3.6.0
首先安装python3.6可能使用的依赖 # yum install openssl-devel bzip2-devel expat-devel gdbm-devel readline-devel ...
- 15.5.2 【Task实现细节】骨架方法的结构
尽管骨架方法中的代码非常简单,但它暗示了状态机的职责.代码清单15-11生成的骨架方 法如下所示: [DebuggerStepThrough] [AsyncStateMachine(typeof(De ...
- [kuangbin带你飞]专题1-23题目清单总结
[kuangbin带你飞]专题1-23 专题一 简单搜索 POJ 1321 棋盘问题POJ 2251 Dungeon MasterPOJ 3278 Catch That CowPOJ 3279 Fli ...
- python-pycharm windows安装
pycharm_IDE安装 1. 首先先pycharm官网,或者直接输入网址:http://www.jetbrains.com/pycharm/download/#section=windows,下载 ...
- 2019-04-16 sql tran and try catch :
begin try begin tran tran_addresource -- 标记事务的开始 delete rp insert into Cube.ResourcePool(ResourceTyp ...
- JavaScript对接百度地图api实现地图标点功能
粗略的做了个地图标点功能 首先,去百度注册开发者账号,然后进入到百度地图开放平台 进入到控制台, 创建浏览器端应用,给个安全域名 然后去开发者文档JavaScript里面找地图展示文档,直接怼上去就行 ...