1基于Garland的边收缩算法的一种实现0引言随着科学技术的进步，在计算机图形学、虚拟现实、地理信息系统、医学图像系统等领域所构造和使用的模型越来越精细、越来越复杂。这些复杂的模型不但对计算机的存储容量、处理速度提出了很高的要求、而且成为实时绘制、网络传输的瓶颈。因此模型简化成为非常重要的研究课题。模型简化是指在保持原模型几何形状基本不变的前提下，采用适当的算法减少该模型的面片数、顶点数和边数。近年来，出现了很多有代表性的模型简化算法，其中Galand的基于二次误差度量的边收缩算法是目前最常采用且有效的算法。其基本思想是以顶点到相关三角形平面的距离的平方和为误差度量，通过重复的边收缩操作对模型进行简化。1算法分析与设计1.1基本概念：定义1三角网格是由三位空间中的三角形通过边和顶点连接而成的分段线性曲面。三角网格M可由顶点集V={v1,v2...vm}和三角集合F={f1.f2...fn}组成的二元组M=(V,F)来表示定义2对M种任一顶点vi，与顶点vi相关的三角形集合记作Planes(i)，与边(vi,vj)关联的三角形集合记作Planes(i,j)，所有与vi关联的边构成的集合Edge(i)。1.2基于二次误差度量的边收缩算法基于二次误差度量的边收缩算法是通过不断选择模型中的一条边进行收缩，达到对模型的简化。每收缩一条非边界边，模型减少2个三角形、1个顶点、三条边；收缩一条边界边，模型减少1个三角形、1个顶点、俩条边。1.2.1误差度量简化模型必须与原网格尽量相似，这取决于边收缩的顺序和边收缩后生成的新点的位置。如何选择合适的边进行收缩及如何生成新的顶点，有一个选择误差度量标准的问题。Garland算法以点到平面的距离为误差度量标准。设对边(vi,vj)进行收缩，则与(vi,vj)边相关联的三角形集合Planes(i,j)构成了原模型上的一个区域，设边收缩后生成的新位置v为[x,y,z,1]T，定义这次边收缩带来的新误差△(v)为v到三角形集合Planes(i,j)中每个三角形所在面的距离的平方和，表示三角形集合Planes(i,j)中的每个三角形所在面的平面方程ax+by+cz+d=0，且有1222cba。在根据点到平面ax+by+cz+d=0的距离公式为222dcbadczbyax以及等式运算,就有△(v)=2),()(vpjiPlanespT其中：p=(a,b,c,d)T该式可变换如下：2),()()]1,,,([)(jiPlanespTTTvppvzyxv),()(jiPlanespTvKpvvQvT其中：pK为4*4的对称矩阵，称为三角形的误差矩阵，它的定义如下Q称为该次收缩的误差矩阵，定义如下1.3本小组关于Garland算法的一种实现。为了简化算法实现的效率问题,使用本组关于Galand算法的一种实现:1、取边收缩的误差矩阵Q为边的俩顶点收缩的误差矩阵Qi和Qj之和，即(Q=Qi+Qj)。2、取新顶点位置为收缩边(i,j)中的一个，误差由公式QvvT计算得到，具体的选择由i收缩到j的所产生的误差和由j收缩到i所产生的误差值做比较得到的较小值决定。3、实际边收缩带来的结果是删除一个点，删除一条边，删除俩个与边关联的三角形，增加新边、三角形。具体算法步骤如下：步骤1根据给定的模型数据，计算出所有顶点的顶点误差矩阵Qi：为各平面对应的a、b、c关联矩pK阵的和:计算过程涉及求面的a、b、c(注意满足1222cba)，即为面的单位法向量,单位法向量的求法可根据平面俩点的叉乘在单位化即可。步骤2设计一关于结构体的优先队列来简化堆的操作，存放i，j，及其收缩带来的误差，代表从i折叠到j的误差,并设计一哈希表用于判定是否新边。3步骤3计算边折叠误差：根据步骤2，取边折叠从i到j和从j到i的教小者作为该边的误差入队并进入哈希散列。步骤4选择队列中的队头元素q出队(即代表从q.i折叠到q.j)。步骤5判断i，j同时相关的三角形的个数：如果个数大于1，那么删除与i，j同时关联的三角形，删除点i，进行边收缩(将剩余三角形中与i相关的三角形全部转移到与j相关),并往队列中增加由于折叠带来的新边(做与步骤3相同的操作)，否则转步骤4。步骤6如果边收缩达到给定的收缩哦边的个数，则结束收缩，否则转步骤4。2、实验分析2.1.1原始牛模型:5804个三角片2.1.2简化后，三角片2804简化时间2578ms4，三角片1404简化时间3235ms,，三角片704简化时间3500ms52.2..1原始兔子模型19859个三角片三角片3859简化时间40375ms62.2.2进一步简化...
