体绘制光线投射加速算法研究与实现.docVIP

体绘制光线投射加速算法研究与实现 　　摘要：本文在光线透射法的各个环节提出了加速算法，包括简化光照模型，法向量明暗度查找表，减少投射光线，加快求交及采样，简化三维几何变换及快速插值。 　　关键词：体绘制 光线投射 明暗度 三维变换 　　中图分类号：TP391文献标识码：A 文章编号：1007-9416（2011）05-0045-03 　　 　　1、引言 　　三维重建技术能充分利用CT、 MRI等医学图像体数据，采用面绘制或体绘制的成像算法，可根据需要得到任意视角透视的三维投影图像，便于医生对人体内部结构进行观察。 　　医学体数据三维重建方法主要分为面绘制和体绘制两大类[1][2]。面绘制基本思想：首先由三维空间数据场构造出曲面、平面等几何图元，通过几何单元拼接拟合物体表面来描述物体三维结构, 是基于二维图象边缘或轮廓线提取，然后再由传统的计算机图形学及硬件实现画面绘制。这类方法可以产生比较清晰的等值面图像，可以利用图形硬件实现绘制，绘制速度快，但不能完整的反映数据域中所包含的信息，并且这种方法需要对体数据进行两分类，即判别每一个体素是否在当前绘制的面上，因此在处理复杂的、边界模糊的人体组织时，经常出现分类上的错误。 　　体绘制基本思想：直接对所有的体数据进行明暗处理而在屏幕上产生具有三维效果的二维图像，即直接将体素投影到显示平面。这类方法不丢失数据所包含的信息，能产生三维数据场的整体图像，清楚的显示内部结构，无须进行分割，有利于保留三维医学图像的细节信息，并且图像质量高、便于并行处理。但需对所有体素进行处理，加大了计算开销，也难于利用图形硬件实现，限制了绘制速度。 　　鉴于体绘制方法的种种优点，此种方法以成为三维重建技术研究中的热点，其中最有代表性的是光线透射法，虽然它的成像质量较高，但速度较慢，所以如何改进此算法，提高利用光线投射法进行三维重建的成像速度具有重要意义。 　　2、光线投射体绘制算法基本原理 　　光线投射算法[1][2]的流程如图1所示。 　　该算法假定三维空间数据f(xi,yj,zk)分布在均匀网格或规则网格的网格点上（为简化起见，在插图中三维空间数据用f(i,j,k)表示。流程图中的数据预处理包括原始数据格式的转换、图像滤波及断层间插值等功能。接着，进行数据值分类（即体数据分割），其目的是根据数据值的不同，正确地将其分为若干类并给每类数据赋予不同的颜色值和不透明度值，以求正确地表示多种物质的不同分布或单一物质的不同属性。然后是重采样，即从屏幕上的每一个象素点根据设定的观察方向发出一条射线，这条射线穿过三维数据场，沿着这条射线选择K个等距的采样点，并由距离某一采样点最近的8个数据点的颜色值和不透明度值作三次线性插值，求出该采样点的不透明度值及颜色值。显而易见，在作重新采样前，需要将具有颜色值及不透明度值的三维数据场由物体空间坐标转换为相应的图像空间坐标。这一算法的最后一步是图像合成，即将每条射线上各采样点的颜色值和不透明度值由前向后或由后向前加以合成，即可得到发出该射线的象素点处的颜色值。重新采样和图像合成是按屏幕上每条扫描线的每个象素逐个进行的，因而这一算法又称为图像空间扫描的体绘制算法。光线投射算法中体素的不透明度在0和1之间变化，可以表示透明度不同的半透明体。1表示完全不透明。 　　3、光线透射算法的加速 　　3.1简化光照模型及法向量明暗度查找表 　　对医学图像进行体绘制过程中，要进行真实感处理，包括颜色，光照，明暗色调和纹理等，图1中给数据点赋颜色值其实就是一种真实感处理，如果绘制的是彩色图像，要给体数据中的每个体素赋RGB颜色值，如果仅绘制灰度图象，就是我们所说的明暗度处理。体绘制中的明暗度处理要模拟光照射到物体表面情况下的反射和投射现象。 　　综合了源-衰减模型、朗伯漫反射余弦定律和Phong镜面反射模型之后，我们得出计算每个体素v(x,y,z)明暗度的公式： 　　式中：Kf为常数;f(x,y,z)为体素v(x,y,z)的密度值;Ip为光源强度;Kd为漫反射系数;Ks为镜面反射系数;q为决定表面光洁度的因子，q越大，光洁度越高;d(x,y,z)为体素v(x,y,z)至观察者的距离;K1,K2为常数;L为光源方向的单位法向量;R为镜面反射方向的单位向量;V为观察方向的单位向量;N(x,y,z)为体素v（x,y,z）处的表面单位法向量。利用灰度梯度法估计多灰度值体数据中一点的法向量,即 　　（1）式中R矢量可通过点积N?L得到矢量L在法向量方向上的投影，R+L=(2N?L)N,于是，镜面反射矢量可由式：R=(2N?L)N-L一计算出来。直接运用式(1)进行明暗处理，计算量太大。我们假设光源在无限远，光源发出的光被近似看作平行光线，则L近似为常向量；又由于我们作的投