《第四讲三维标量场数据可视化》由会员分享，可在线阅读，更多相关《第四讲三维标量场数据可视化（104页珍藏版）》请在金锄头文库上搜索。
1、7(tnee Q皿。成况三维标量场可视化三维空间数据场方法三维空间数据场与二维数据场不同，它是对三维空间中的采样，表示了一 个三维空间内部的详细信息，这类数据场最典型的医学CT采样数据，每 个CT的照片实际上是一个二维数据场，照片的灰度表示了某一片物体的 密度。将这些照片按一定的顺序排列起来，就组成了一个三维数据场。此 外，用上规模计算机计算的航天飞机周围的密度分布也是一个三维数据场 岛例子。三维空间数据场方法主要分为:1抽取表面信息的可视化方法（面绘制）：分为断层间的构造等值面值面生成2.直接体绘制方法（体绘制）：光线投射.投影方法其它体绘制方法抽取素面信息的可视化方该（面條制丿断层间的构造
2、等值面（断层间表面重构）如CT采样数据场这样的三维数据，可以看成是由一些二 维数据场按一定顺序排列组成的，各断层数据之间有很大 的相关性。断层数据广泛存在于医学、生物、地质、无损 探伤等应用领域，其各断层间相互平行，每一断层与实体 的交线就是实体在该断层的轮廓线。如果先在各层之间找 出物体的边界线，再利用断层之间的连贯性，就可以从一 系列断面上的轮廓线中推导出相应物体的空间几何结构。斷点间的构造普值而（断层间表面重构丿主要应用领域是医疗成像 不同的扫描技术包括： 囚CT （计算机断层扫描）扫描仪给平均 一个地区价值-而 不是在一个点的值MRI MRI （核磁共振成像） 多个二维切片的三维图像切
3、片切片间的空隙:丁 切片三二瞬层间的构造等值面（断层间表面重构丿CTMRISPECT单光子发射计算机断层成像术断层数据的表面重构是从一系列断面上的轮廓线 推导出实栋的空间几何结构相邻层的轮廓线位于相互平行的两个平面图51实体与断层瞬层间的构造等值面（断层间表面重构丿如果在相邻两层，各自只有一条轮廓线，其三维重构问题称为单轮廓线重构问题如果在相邻两层，有多条轮廓线，其三维重构问 题称为多轮廓线重构问题图7- 1单轮廉线重构问题示意图图7.2多轮螂线重构问题示意图瞬层间的构造等值面（断层间表面重构丿在一个断层中找出物体的轮廓线可以利用上面介绍的等值线方法。 找到所有轮廓线后，第二步是在各个相邻的轮
4、廓线之间构造出物体的 表面，然后进行绘制。物体的表面可以用三角面片拼接出来，拼接的 方法如图所示，就是在相邻的两层上找出三个点，其中两个点在同一 层，另一个点在另一层。在拼接过程中，一次加入一条边，就可以组 成一个三角面片，但加入一条边有两种选择，例如图中P1Q2和P2Q1 ，如果选择不恰当，则拼接出的表面比较乱，也不光滑。最简单的选 择方法可以采用贪心方法，就是每次选择一条较短的边加入，这样可 以保证构造出的衰面比较光滑。面绘制算法将三维数据场中具有某种共同属性的采样点按其空间位 置连接起来，构成一张连续表面，然后对抽取出的表面 进行绘制等值面算法等值面：在一给定三维数据场中， 有空间点的集
5、合三维标量场可视化中最常用Marching Cubes 方法采样值均为某一给定值的所著值面笛咸（面條制丿署值面唆啟（面條制丿构造物体的表面也可以采用等值面的方法。等值面可以看成是等值 线的三维扩展。等值面的构造也就是等值线构造方法的三维扩展，最 典型的就是Marching Cube方法。二维数据场的基本单元是矩形,在 三维空间的基本单元是一个小立方体。如果我们找出每个小立方体中 的等值面，这些等值面也就构成了整个物体的表面。表面绘制方法绘制的动物骨骼表面数据集囚适用于三维规则标量场每一立方体单元称为一个体素（voxel）,数据场的数据值分布 在体素的8个顶点上典型代表：CT数据、MRI数据13
6、Marching Cube方法体数据体素每个体素转换 单位立方体顶点no数据增强-最近邻插值100100数据增强-三线性插值三线性插值:ffx V Z)=U(l-x)(l-y)(l-z) + fiooxO -y)(i-z) + foioO-x)y(l-z) + fooiO -x)(l-y)z + f110xy(l-z) + 加(1 -y)z + foiidy fxyz数据增强-三线性插值*- .一. _ .口 表示的数值由如下表示100From University of Bonn建立等值面简单起见，以数值为零的水平等值面为例，并记正顶点为立方体有八,均为正或均为负。共计28 = 256不同情
7、况以下两种最简单此时，立方体内部无等值面部分等值面建立- 一个顶点为正的情况逆线性插值的交点逆线性插值艰据线性插值公式t处的函数 值：f(x*) = f 1 + t ( f2 f 1 )根据f*反求t = (f*-fl)/(f2-fl)Xi t x2xlx2等值面建立- 一个顶点为正的情况连接交叉点，形成了一个三角形 等值面的一部分等值面建立正顶点在对应角处等值面建立通过分析这种方式所有256种情况，尽管它很相差 很大，很多情况下是相似的。3S例如：顶点都为正或为负的2种情况，都没有等值面16种情况，其中一个顶点从所有其他符号相反实际上，仅有15拓扑不同的情况。256种可能的情况可以归纳为这1
8、5个基于典型个案 互补的基础上（顶点为正和负的互换）和旋转对称 这样做的好处是 为便于实施我们 只需要代码15种情况代替256种思想：基于“分治(divide-and-conquer) 思想将整个 数据场的等值面抽取分解到每一个体素中去完成拥算法概述读入三维规则标量场0对于每一体素依据所需抽取的等值面的属性值，确定其8个顶点的状态对于体素的每一条边，依据顶点状态，判别它是否与等值面有交 点。若交点存在，则求出交点在求出了当前体素的所有边与等值面的交点后，依据一定的准则 将这些交点连接成三角形，作为等值面位于该体素内部分的近似 表示，并进行真实感绘制当处理完所有体素后，即完成了整个数据场的等值面
9、抽取与绘 制确定体素顶点状态设所需抽取的等值面的属性值为C。若某顶点U所存贮的数据值大于（或等于）C。,则认为U在等值 面外侧（或位于其上）,并记其状态值为1反乙若V所存贮的数据值小于卬贝0认为V在等值面内侧，并 记其状态值为0#确定体素顶点状态Example : 5个顶点均位于外侧，记为10111100Case = v8|v7|v6|v5|v4|v3|v2|v13判别体素的边与等值面是否有交对于某一条边E （其顶点为巾和v2）,若巾和卩2的状态值相同, 则边E位于等值面的外侧（或内侧），边E不与等值面相交； 反之，若巾和的状态值不同，边E必定与等值面相交 若边E与等值面有交点，可通过线性插值
10、计算出交点（C。2）*岭+（卩一C。）*/（62）# JSu健将体素各边与等值面的交点连接成三角形取决于体素每一顶点的状态值分布情况 存在着2*种不同情况每一体素有8个顶点每一顶点有两种状态值基于体素顶点状态翻转对称性和旋转对称性,将上述256种组合 情形减少到15种翻转对称性：如果体素各顶点的状态值0和1互换，所含等值面的 拓扑结构（即交点连接关系）不变旋转对称性：体素旋转后，所含等值面的拓扑结构不变#15种等值面连接模式0 1 210 11 128933第1步：8个顶点相对的等值面的值分类字节索引；1 +ve;O -ve1000001VI V2 V3 V4 V5 V6 V7 V8代码标识边
11、的相交：VIV4; VIV5; V2V3; V2V6; V5V8; V7V8; V4V8 第2步：查找表标识出等值面交点个数 例如：00000000 10000000 01000000Configuration 0Configuration 1Configuration 1 11000001Configuration 6 11111111Con figuration 0第3步：沿边缘的逆线性插值将找到的交汇点第4步：规范的配置将决定如何创建等值面片（0, 1, 2 ,3或4个三角形）第5步：通三角形的显示渲染必風4算法找等值面优点g等值面提取边效果较好定义为三角形在3D表面-著名的渲染技术提供
12、照 明，阴影和观看硬件支扌#缺点只显示数据切片 有二义性算法Marching cubes 算法缺点第3步：三角形是选择 不唯一?拥存在问题15种连接情形中，有些存在二义性，可能导致所生成的相邻体 素的等值面之间不连续，从而使最终生成的等值面存在空洞”39396)(a)6(b)6i(C)39算法二义性- .二二义性问题，可以判断交叉口一双值寻找鞍点连接V八歧汶原因，因为： 三线性插值是线性沿边缘 就成了一个双线性函数正确的拓扑e构，我 必须加入1个正确的久点、产生歧义的不同三角面存在问题及改进改进方法之一：增加连接模式，使其能与相邻体素的状态相匹 配以消除“空洞”10c12c13c発存在问题及改进改进方法之二：将六面体体素分解为四面体单元，并将等值面 抽取限制在四面体单元中进行424242等值面空洞由于二义性三角面的存在 ,会造成空洞出现。孔洞三线性插值三线性插值函数： f(x,y ,z) = f()oo( 1 -x)( 1 -y)( 1 -z)+fio()x(ly)(lz) +f010(l-x)y(l- z) + fooi(lx)(ly)z +f1J0xy(l- z)+fioiX(ly)z +fOn(l-x)yz +fmxyz相对复杂。淀例如,等值面f(x,y,z) = 0
