基于曲面方程的三角形网格求交算法.docVIP

基于曲面方程的三角形网格模型求交方法 陈振1 汤军1,2 廖环宇1 明廷宝1 梁宏1 陈俊飞1 刘圆圆1 长江大学地球科学学院，湖北武汉 430100） 长江大学油气资源勘探教育部重点实验室，湖北武汉 430100） 摘 要 在应用GIS不规则三角网的地层建模、表达和分析的研究过程中，会涉及到大量复杂的三角网格模型的求交运算。本文提出了一种基于曲面方程的三角形网格模型求交算法。该算法首先通过三维点要素，插值拟合出两个曲线方程，然后求解方程组获得交线方程，再将该交线分别投影到三角形网格上，获得顺序的投影三角形组，最后对两个三角形组做求交运算。该算法能快速有效排除不相交三角形，同时对候选三角形和三角形交线依次记录，大大减少了运算时间。 关键字 三角网 表面模型 曲线方程 交线 投影 中国法分类号：TP ； 文献标志码：A1 引言 随着地质信息计算机模拟技术和可视化技术的发展，三维地层建模及可视化技术是空间地理信息化建设的关键技术之一，并成为数学地质、石油勘探、岩土工程，GIS与计算机技术等多个领域的研究前沿和热点[1]。 在涉及三维地质体模拟与分析的领域中, 表面模型是重要的三维数据模型之一。表面模型是用3D 空间实体的几何表面模拟、表达实体, 侧重于3D 空间实体的外部形态表示, 如地质层面、断层剖面、地形表面等。一般来说, 几何表面模型的表示主要是Rainald Lohner（1996）提出的曲面和多边形曲面两种方法[2]。参数曲面主要由特定的数学方程来确定, 便于编辑和控制, 多用于设计领域；多边形曲面则由分布在实体表面的离散点集生成的多边形面元的集合构成, 离散点一般来源于实际采集的三维空间点信息, 面元大多采用三角形面元, 表面用不规则三角网（TIN）表示。 多边形曲面广泛应用于工程领域的三维可视化和分析计算。 本文主要讨论两个曲面相交的问题，，传统的方法是将曲面模型所包含的所有三角形进行两两相交测试，如毕林、王李管等在2008年提出的网格模型布尔运算算法，但随着三维技术的不断发展，三维网格模型的复杂度不断增加，所包含的三角形数量也不断飙升[3]。使得剔除不相交三角形减少求交的计算量，提高不规则三角网求交算法效率的有效途径之一有潘振宽、崔树娟等提出的层次包围盒法[4]刘丽萍、王琰提出的包围盒编码法[5]，张少丽、王毅刚等提出的空间分解法[6]等。 层次包围盒法需要将模型对象建立最小包围盒，某情况下求交效率不太理想。包围盒编码法则是通过空间编码将三角形相交情况逻辑化，但在剔除不相交三角形上效果不理想。空间分解法虽然在整体上剔除不相交三角形效果较好，但在剔除同一空间网格的三角形时任需要借助其他的方法，有时需要用到包围盒法简化结果，使得算法复杂度增加，运算效率不太高。 针对以上方法的不足，如果能找到一种方法，不但能在三角形求交之前删减三角形数量，而且在获得两个三角形交线时顺序存放，那将能更快连接成交线，得到结果。 基于此思想，本文将提出一种新的曲面求交算法：基于曲面方程的三角形网格求交算法。首先通过组成三角形网格的三维点，插值拟合出两个求交三角形网格的曲线方程，求解方程组获得交线方程，将该交线分别投影到两个三角形网格上，获得顺序的投影三角形组，最后对两个三角形组做求交运算。 交线作为两个曲面相交部分，在表面模型的构建以及后续的空间分析中都起着至关重要的作用。因此，可靠快速的不规则三角网求交法，是三角网表面模型建模分析的基本算法之一。表面模型构建过程中，交线是不同表面相交的边界部分，一定会构成一个有序的线段链或者线段环。于是对于不规则三角网求交算法则主要包括三部分：（1）确定相交不规则三角形对；（2）求取对应三角形对的交线线段；（3）把各交线线段连接成链或者环[7]。 首先获取生成三角形的具有三维坐标的点要素k=1,2,3，…，N，再应用插值算法[8]将点要素拟合成面要素，并得出曲面方程z=fa(x,y)和z=fb(x,y)，联立方程，求解得到交线方程式y=g(x)和z=h(x)。 2.2.1 求取交线与三角形网格模型的一个交点，作为开始获取交线的起点 将求交的两个三角形网格模型记分别为V1、V2，获取中的V1，求得其三角形总数k，并通过顶点坐标求得每个三角形面的参数方程 z=fm(x,y)（x1xx2,y1yy2,m=1,2,…,k），算法流程图如图1所示： 图1 获取起点和初始三角形流程图 Figure 1 Flow chart of obtaining the starting point and initial triangle 2.2.2 从△P1P2P3开始求得交线投影的三角形组N[n] 从交点p(x,y,z)开始，做交线在该△P1P2P3面上的投影，投影线必经过△P1P2P3的一条边上一点A