Topology and Geometry in OpenCascade
Location and Orientaion
eryar@163.com
摘要Abstract：本文简要介绍了几何造型中的边界表示法（BRep），并结合程序说明OpenCascade中的边界表示的具体实现，即拓朴与几何的联系。拓朴结构中的位置（Location）和朝向（Orientation）进行了详细说明。
关键字Key Words：OpenCascade、BRep、Topology、Geometry、Location、Orientation
一、引言 Introduction
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
二、位置 Location
ToposDS_Shape有个TopLoc_Location的成员变量myLocation，该变量定义了子形状相对于该形状的偏移量。例如，环（wire）有一个位置变量，该变量沿着向量{0,0,10}移动，意味着环所有的边都沿着Z轴移动10个单位。下面以一个程序来具体说明：
1 /*
2 * Copyright (c) 2013 eryar All Rights Reserved.
3 *
4 * File : Main.cpp
5 * Author : eryar@163.com
6 * Date : 2013-09-26
7 * Version : V1.0
8 *
9 * Description : Shape location.
10 *
11 */
12
13 #define WNT
14 #include <Geom_Circle.hxx>
15
16 #include <TopoDS_Edge.hxx>
17 #include <TopoDS_Wire.hxx>
18 #include <TopoDS_Iterator.hxx>
19
20 #include <BRepBuilderAPI_MakeEdge.hxx>
21 #include <BRepBuilderAPI_MakeWire.hxx>
22
23 #pragma comment(lib, "TKernel.lib")
24 #pragma comment(lib, "TKMath.lib")
25 #pragma comment(lib, "TKG3d.lib")
26 #pragma comment(lib, "TKBRep.lib")
27 #pragma comment(lib, "TKTopAlgo.lib")
28
29 const std::string dumpShapeType(const TopAbs_ShapeEnum& type)
30 {
31 std::string strType("Shape");
32
33 switch (type)
34 {
35 case TopAbs_COMPOUND:
36 strType = "COMPOUND";
37 break;
38
39 case TopAbs_COMPSOLID:
40 strType = "COMPSOLID";
41 break;
42
43 case TopAbs_SOLID:
44 strType = "SOLID";
45 break;
46
47 case TopAbs_SHELL:
48 strType = "SHELL";
49 break;
50
51 case TopAbs_FACE:
52 strType = "FACE";
53 break;
54
55 case TopAbs_WIRE:
56 strType = "WIRE";
57 break;
58
59 case TopAbs_EDGE:
60 strType = "EDGE";
61 break;
62
63 case TopAbs_VERTEX:
64 strType = "VERTEX";
65 break;
66
67 default:
68 break;
69 }
70
71 return strType;
72 }
73
74 void dumpShapeLocation(const TopoDS_Shape& shape)
75 {
76 std::cout << "Shape Type: " << dumpShapeType(shape.ShapeType()) << std::endl;
77 shape.Location().ShallowDump(std::cout);
78
79 TopoDS_Iterator anItr(shape);
80
81 for (; anItr.More(); anItr.Next())
82 {
83 const TopoDS_Shape& aChild = anItr.Value();
84
85 dumpShapeLocation(aChild);
86 }
87 }
88
89 int main(void)
90 {
91 Handle_Geom_Curve aCircle = new Geom_Circle(gp::XOY(), 5.0);
92
93 TopoDS_Edge anEdge = BRepBuilderAPI_MakeEdge(aCircle);
94 TopoDS_Wire aWire = BRepBuilderAPI_MakeWire(anEdge);
95
96 std::cout << "Before transformation: " << std::endl;
97 dumpShapeLocation(aWire);
98
99 gp_Trsf trsf;
100 trsf.SetTranslation(gp_Vec(0, 0, 10));
101 TopLoc_Location location(trsf);
102
103 aWire.Location(location);
104
105 std::cout << "After transformation: " << std::endl;
106 dumpShapeLocation(aWire);
107
108 return 0;
109 }
程序结果如下所示：
1
Before transformation:
2Shape Type: WIRE
3TopLoc_Location : Identity
4
5Shape Type: EDGE
6TopLoc_Location : Identity
7
8Shape Type: VERTEX
9TopLoc_Location : Identity
10
11Shape Type: VERTEX
12TopLoc_Location : Identity
13
14After transformation:
15Shape Type: WIRE
16TopLoc_Location :
17 Exponent : 1
18 TopLoc_Datum3D 034100E8
19 ( 1, 0, 0, 0)
20 ( 0, 1, 0, 0)
21 ( 0, 0, 1, 10)
22
23
24Shape Type: EDGE
25TopLoc_Location :
26 Exponent : 1
27 TopLoc_Datum3D 034100E8
28 ( 1, 0, 0, 0)
29 ( 0, 1, 0, 0)
30 ( 0, 0, 1, 10)
31
32
33Shape Type: VERTEX
34TopLoc_Location :
35 Exponent : 1
36 TopLoc_Datum3D 034100E8
37 ( 1, 0, 0, 0)
38 ( 0, 1, 0, 0)
39 ( 0, 0, 1, 10)
40
41
42Shape Type: VERTEX
43TopLoc_Location :
44 Exponent : 1
45 TopLoc_Datum3D 034100E8
46 ( 1, 0, 0, 0)
47 ( 0, 1, 0, 0)
48 ( 0, 0, 1, 10)
49
50
51Press any key to continue . . .
三、朝向 Orientation
朝向（orientation）与位置（location）的工作原理相同。当将子对象从实体中分离出来时，父对象的朝向会影响到子对象的朝向。但是有个很重要的例外就是面（Face）上边（Edge）的朝向不遵守这个规则。在讨论面的朝向时，讨论过边的参数空间曲线（pcurve）的material问题。这个例外说的是计算面上边的朝向时不应该受到面的朝向的影响。即若要使用边的参数空间曲线（pcurve）就按下面的方式：
1 TopExp_Explorer aFaceExp (myFace.Oriented (TopAbs_FORWARD), TopAbs_EDGE);
2 for (; aFaceExp.More(); aFaceExp.Next())
3 {
4 const TopoDS_Edge& anEdge = TopoDS::Edge (aFaceExp.Current());
5 }
假如你研究得更为深入，这个例外也是可以理解的。让我们以一个具体例子来理解这个例外，从底层一步步的创建一个面：
1 Handle(Geom_Surface) aSurf = new Geom_Plane (gp::XOY());
2 //anti-clockwise circles if to look from surface normal
3 Handle(Geom_Curve) anExtC = new Geom_Circle (gp::XOY(), 10.);
4 Handle(Geom_Curve) anIntC = new Geom_Circle (gp::XOY(), 5.);
5 TopoDS_Edge anExtE = BRepBuilderAPI_MakeEdge (anExtC);
6 TopoDS_Edge anIntE = BRepBuilderAPI_MakeEdge (anExtC);
7 TopoDS_Wire anExtW = BRepBuilderAPI_MakeWire (anExtE);
8 TopoDS_Wire anIntW = BRepBuilderAPI_MakeWire (anIntE);
9 BRep_Builder aB;
10 TopoDS_Face aFace;
11 aB.MakeFace (aFace, aSurf, Precision::Confusion());
12 aB.Update (aFace, aSurf);
13 aB.Add (aFace, anExtW);
14 //material should lie on the right of the inner wire
15 aB.Add (aFace, anIntW.Reversed());
面默认的朝向是向前的（forward），让我们来遍历面的边（edge）和参数空间曲线（pcurve）。尽管我们没有显示地来添加它们，从原来的讨论中可知平面上的参数空间曲线可以默认计算（be computed on the fly）：
1 void TraversePCurves (const TopoDS_Face& theFace)
2 {
3 TopExp_Explorer anExp (theFace, TopAbs_EDGE);
4 for (; anExp.More(); anExp.Next())
5 {
6 Standard_Real aF = 0.0;
7 Standard_Real aL = 0.0;
8 const TopoDS_Edge& anEdge = TopoDS::Edge (anExp.Current());
9
10 Handle(Geom2d_Curve) aPCurve = BRep_Tool::CurveOnSurface (anEdge, theFace, aF, aL);
11 }
12 }
得到的参数空间曲线（pcurves）如下图所示，material在红线的左侧，在蓝线的右侧：
一切都很正确。现在设想一下，如果把面的朝向反向（reverse），然后再遍历边和参数空间曲线，会发生什么呢？
1 TopoDS_Face aRFace = TopoDS::Face (aFace.Reversed());
2
3 TraversePCurves (aRFace);
所有的边将会具有相反的朝向，对应边的参数空间曲线（pcurve）的material在外环的外侧，在内环的内侧，这明显是错误的！在前面讨论面的朝向时就说过面的朝向仅仅是面的逻辑朝向，而与其底层的曲面（surface）没有关系。在上面的例子中反转面aRFace只是一个法向为{0， 0， －1}的面。所以，要获得边的正确朝向，必须使用下面的方法来访问面中的边：
1 TopExp_Explorer anExp (theFace.Oriented (TopAbs_FORWARD), TopAbs_EDGE);
这样就确保面上的边具有正确的朝向，而与曲面（surface，注意在此不是face！）的法向没有关系。OpenCASCADE的算法对这种特殊的情况都做了处理，你也一定记得这样做。
四、结论 Conclusion
对拓朴形状（TopoDS_Shape）的位置（location）和朝向（orientation）进行深入理解，整个拓朴结构就变得清晰了，因为一个拓朴形状中除了子对象外，剩下就是位置和朝向成员变量了。
理解了拓朴结构后，对OpenCascade的模块ModelingData就有个较深刻地认识了。
五、参考资料 Reference
1. Roman Lygin, OpenCascade notes, opencascade.blogspot.com
2. 孙家广等. 计算机图形学. 清华大学出版社
3. OpenCascade source code.