三角网格细分法在重建隧道三维结构中的应用

1、三角网格细分法在重建隧道三维结构中的应用(1)摘 要 曲面造型一直是计算机图形学的重要研究内容之一。曲面的显示效果涉及到数据的存储结构及对数据的操作方法等多方面的知识。网格细分是曲面的一种表示方法，但却是目前广泛使用的一种曲面造型方法。而三角网格细分的许多算法已经成功地应用于许多领域。本文介绍了曲面造型的相关理论，重点介绍了 Loop 细分算法及其在重建隧道三维结构中应用。关键词 曲面;网格;细分;逼近;隧道0 引言隧道的施工安全，常常与隧道的掌子面观察紧密联系。现行观察一般是人工纪录和掌子面状况描绘或者照像，这种描绘往往由于人为因素，使获得的数据缺乏一致性和可比性。而照片在冲洗过程中，由于相

2、片本身的分辨率也会导致冲照片在洗过程中出现信息畸变、丢失的现象，而且不能数字化和及时处理。因此在对隧道开挖过程中的掌子面图象进行处理后，就需要根据所获得的特征参数，重建、显示隧道的三维结构，并用以指导隧道的施工。但三维结构在重建过程中，重建的效率是十分重要的。本文利用曲面网格的相关理论，并结合三角网格细分算法来进行遂道三维结构的重建，以加快生成速度和提高隧道三维结构显示的质量。 1 曲面网格的构建空间曲面的构建一般采用参数方式。若用双参数构建曲面，其数学形式为：S=S(u,v)=x(u,v),y(u,v),z(u,v) (aub,c vd)当 u,v 在各自的定义域内变化时， S（u,v ）在

3、空间坐标系中变化，且O-UV坐标系中的任意一点 (u,v) 均与 O-XYZ坐标系中的点 (X,Y,Z) 呈一一映射关系。当 v 不变而 u 变化时，得到 u 线；反之，则得到 v 线。所有的 u 线和 v 线所形成的网叫参数曲线网。由曲线 p(u,c),p(u,d),p(a,v),p(b,v) 四条边界曲所围成的部分，被称为网格，它定义在矩形域（ (a ub,c vd) ）上。矩形定义域一般采用正方形（ 0u1,0 v1）。 Coons曲面， Bezier 曲面， B 样条曲面均可用这种方式来建构。曲面建构方法亦可用三角域方式来构建。选取三角形区域 （u,v,w ）|u 0,v 0,w0,u

4、 v w=1,并采用 B网来构建曲面。对于不共线的三点P0,P 1,P 2 可以构成一个三角形，该三角形所在的面上任一点P 可表示为 P（ u,v,w ）=uP0 vP 1 wP2 (u v w=1)若当 0u v w 1, 则点 P(u,v,w) 位于 P0，P1， P2 所构成的三角形之内。例如三角域上的 n 次 Bezier 曲面（三角域上的三次 Bezier 曲面如图 1）为（1）Bijk n (u,v,w) 是与点 P ijk 对应的基函数，且图 1 三角域上三次Bezier 曲面B 样条曲面亦可采用三角域来构建，但基函数形式与建构 Bezier 曲面的基函数不同。 2 网格的细分（

5、 Subdivision ） 2.1 曲面造型一直是计算机图形学的重要研究内容之一。网格细分(Subdivision)是曲面的一种表示方法, 是曲面造型的一个分支。该方法实际上是从一个被称为控制网格 ( 绝大多数网格数据可用数字化仪通过人工模型来输入 ) 的多面体开始 , 递归地计算新网格上的每个顶点 ( 这些顶点一般是其上一级细分网格上某几个顶点的加权平均 ), 从而产生新的小面来获得更好的图形显示。这种方法只要计算机的硬件条件允许，就可以不断进行下去。针对网格构建方法，产生了网格细分的许多算法（网格细分算法的思想可以追溯到 50 多年前 G. de Rham 用于生成二维光滑曲线的“砍角算

6、法”）。常见的有Catmull-Clark算法、 Loop 算法、蝶形算法、 Doo-Sabin 算法， 算法等。一般而言，网格的细分造型可分为两步 : 拓扑分裂和几何平均。拓扑分裂定义了细分增加的网格点的拓扑结构 , 几何平均确定网格点的几何位置。 2.2 三角网格细分法的特点 : 任意拓扑、可伸缩性、表示的一致性、数值稳定性、简单性。 正是由于三角网格具有上述特点 , 使得三角网格细分方法作为一种 CAGD技术具有很强的生命力。 3 经典的网格细分算法 -Loop 算法细分算法是美国犹他大学的 Charles Loop 于 1987 年在其 3.1 Loop硕士论文中提出的一种逼近型三角形

7、面分裂细分算法。它是基于B 样条逼近型算法的一种策略 , 应用于规则网格时可以产生C2 连续的曲面 , 在非正规点处则可1达到 C 连续。该算法使每个三角片在分裂后生成四个小的三角面片，如图2。图 2 一个三角面片上的 Loop 细分格式根据点生成方式的不同可将新生成的点分为现两类 :(1) 第一类控制点：由原三角形的边生成的对应控制点，如E1， E2 等；(2) 第二类控制点：由原三角形的顶点生成的对应控制点，如等。 3.2 对于规则的三角网格， Loop 细分格式如图 3 所示：(1) 第一类点（ 2）第二类点 图 3 规则网格上 Loop 细分格式下各顶点的权由原三角形的边生成的第一类控

8、制点由下式进行计算：(2)由原三角形的顶点生成的第二类控制点由下式进行计算：(3)V 3 是与第一类点相邻的两个三角形的顶点V其中， V 0、V 1、V 2,、是原始控制顶点 , Q 是与 V 有一条边相连的原始控制顶点的和的平均。3.3对任意三角网格 , 运用相同的规则生成第一类控制点和第二类控制点 ,但第二类控制点的权修改为如图 4 所示。图 4非规则网格上第二类控制点的权此时公式（ 3）应为：（ 4）其中 。当 n=3 时， , 当 n>3 时， ； Vi ( i=0,1,2, ,n -1) 是与 V 有一条边相连的原始控制顶点。 4 总结三角网格细分方案是目前广泛应用的曲面建构方

9、法 , 在隧道掌子面三维重建过程中如果能有效地利用三角网格细分方法 , 不但可以有效地减少数据的存储量 , 而且可以改善建模的生成速度和模型的显示效果。更重要的是不仅能够根据获得的特征参数，重建、显示隧道的三维结构，而且能够根据重建的隧道的三维结构，进行一定范围内的地质情况的预测，从而更加安全、高效地指导隧道的施工。 参考文献 1 Loop C. Smooth Subdivision Surfaces Based on Triangles. Master s thesis,Department of Mathematics,University of Utah,1987 2 Loop C. Triangle Mesh Subdivision with Bounded Curvature and the Convex Hull Property. Microsoft Corporation,2001 3张学辉 . 三角网格细分技术在三维地质模型中的应用. 科学技术与工程 J. 2006年 1 月第 1 期,84-85 4 王金生等 . 几种经典网格细分算法的比较 . 计算机应用研究 J. 第六期 , 139-14