等值线等值面生成.pptVIP

Marching Cubes (MC)方法 3．MC方法存在的问题 (1) MC方法构造的三角面片是三维等值面的近似表示 首先，在MC方法中，等值面与体元边界的交点是基于函数值在体元边界上呈线性变化这一假设而求出的。当数据密度高、体元很小时，这一假设接近于实际情况。但是，在稀疏数据中，体元较大，如果仍然认为函数值在体元边界上呈线性变化，将会产生较大误差。这时，需要根据不同的具体情况对函数值沿体元边界的变化作其它适当的假设，才能较准确地求出等值面。 其次，即使函数值沿体元边界作线性变化这一假设符合实际，那么通过线性插值求得的交点位置是准确的，但是，将体元中同一个面上两条相邻边上的交点简单地用直线连接起来也是一种近似(如下图所示)。 第32页/共60页 第33页/共60页 为了说明这一问题，需要引入当体元各边界上函数值均为线性变化时的等值面模型。如图3.5所示，P(x, y, z)为小体元中的任意点，体元中的数据沿x, y, z三个方向均是线性变化的。如果点P1，P2为点P沿y轴在立方体两个面上的投影，P11、P12、P21、P22分别为P1，P2点沿z轴在立方体平面上的投影。设V为y轴上的坐标分量，f为函数值，那么，通过三次线性插值，可得： (1) 其中P1，P2两点的值可由P11，P12和P21，P22插值求得，而P11、P12、P21、P22四个点的值又可以由它们所在体元内的一条边上的两个顶点插值得到。这样，通过三次线性插值运算，就可以求得P(x, y, z)点的函数值，（1）式可具体展开为： 第34页/共60页 其中系数ai ( i= 0, …7 )取决于体元8个顶点处的函数值，如果给定的等值面的值为Ft，那么，等值面就被定义为满足如下方程的点的集合 (2) 改变Ft的值，就可以得到不同等值面的表达式。 由上述等值面方程可以方便地求出某等值面与体元边界面的交线方程。不失一般性，设某边界面所在平面的方程为z=z0，代入方程式(2)，可得 (3) 上式可进一步表示为： 第35页/共60页 （4） 显然，上述方程表示的是一条双曲线，即等值面与体元中某一个面的交线是一条双曲线或其中的一支。如果用一条直线来表示这条双曲线，则会引起误差（如图所示）。如果体元很小，这一误差是可以忽略不记的。对于稀疏的三维数据，这种近似引起的误差是难以接受的，可通过自适应剖分算法将三角形按给定的逼近精度递归地分成子三角形，使这些子三角形的顶点满足方程(3)，且子三角形与等值面的最大距离小于给定的容差。 第36页/共60页 Marching Cubes (MC)方法 (2) 连接方式上的二义性 Marching Cubes方法可以看成是二维等值线网格序列法在三维空间中的推广。在网格序列法中，如果矩形网格单元4个顶点中有两个顶点的值大于等值线的值，另两个顶点的值小于等值线的值，且这两个顶点交叉分布，那么等值线的连接就出现二义性。同理，在Marching Cubes方法中，如果正六面体单元的6个矩形表面中出现与此相同的情形，那么该正六面体单元中的等值面连接必然会出现二义性，这样的面称为二义性面，包含1个以上的二义性面的体元，即为具有二义性的体元，如下图所示。事实上，在Marching Cubes方法的15种模式中，第3、6、7、10、12、13等6种情况是具有二义性的。 第37页/共60页 第38页/共60页 Marching Cubes (MC)方法 4．用双曲线渐近线方法判别和消除二义性 MC方法存在连接方式的二义性问题在该方法提出后不久，就由M.J.Durst提