数字化放射图像重建技术在癌症放疗中的应用.docxVIP

数字化放射图像重建技术在癌症放疗中的应用 1 直接体绘制算法 通过模拟x-ligh-t检验原理，获得与数字ct数据相同的x-ligh-ct图像效果。主要原理是在特定的介质中模拟x-ligh-t信息时产生的能量损失，这些数据被光学吸收和色散成像所包围。在成像领域，ct数据被视为具有源衰减效应的颗粒集。通过适当的传输函数设计和x-ligh-t广告模拟，可以直接绘制drr。两种传输函数的结构不同：直接映射源的传输函数由体素梯度决定或由用户互动确定，但数字辐射图像源的传输函数必须符合物理规律，将体素值转换为与身体和组织的实际对应器官和组织的x-ligh-l。 按照处理顺序,直接体绘制算法大致可以分为三类:图像空间方法(如光线投射),物体空间方法(如Cell projection和Splatting)和频域空间方?法. 频域空间方法虽然能够利用频域快速积?分?将数据遍历的计算复杂度从O(N3)降低至O(N2logN),但其不支持透视投影,无法获得术内应用所必须的深度提示. Splatting算法利用预积分footprint表获得交互的绘制速度,但不适合于DRR. 光线投射算法虽然能产生高质量图像,在绘制速度及交互性能上却不尽如人意.为了使直接体绘制技术得到实用,人们在这些算法上做了各种改进.纯粹的基于软件的加速方法有Space leaping,Early ray termination,Ray template和Shear-warp facterization等.在为数不多的有关DRR的文献中,蔡文立等提出了一种结合Shear-warp和Splatting算法的混合加速技巧,在平行投影条件下能够将光线投射算法绘制速度提高到原先的5.6倍左右. LaRose提出一种称为“Transgraph”的基于图像绘制的方法,由于牺牲了精度,用该方法生成的图像有严重的走样现象;同时,他们建议采用纹理硬件加速绘制,但受图形硬件对帧缓冲区存取效率的限制,该方法在目前并不可行. Wang等提出一种基于圆柱调和映射快速DRR技术. 虽然他们宣称能够达到实时的绘制速度,但从实验数据来看,该算法在Dell Precision420工作站上以2.21 fps的速度绘制大小为330×330×83的体数据,远远没有达到实时或交互绘制的要求. 随着可编程图形硬件的发展,出现了一系列加速直接体绘制的算法,包括二维和三维纹理映射,大规模并行计算,CPU优化和专有硬件等.本文提出一种新的硬件实现光线投射算法,结合三维体纹理、可编程顶点和像素渲染器及最新render-to-texture功能,获得快速、高质量的数字化放射图像效果. 2 研究背景 2.1 drr中mi的概念 X-光线衰减的物理原理构成了DRR绘制的基础.可以用一个指数函数来模拟当X-光线穿过某一物体时,在距离X-光线入射点的s处,X-光线的强度衰减为其中,I0为入射点处的X-光线强度,μ为线性衰减系数,E表示在某点t处的光子能量.式(1)的离散形式为 考虑到X光片的曝光效果,DRR技术以负片的形式给出最终结果 I=I0(1?exp(?∑i=0NμiΔs))=I0(1?exp(?∑i=0NmiρΔs))(3)Ι=Ι0(1-exp(-∑i=0ΝμiΔs))=Ι0(1-exp(-∑i=0ΝmiρΔs))(3) 其中,mi表示质量衰减系数,ρ为物体密度.在DRR中,衰减系数传输函数可用一个分段线性函数表示,用于模拟各种人体器官组织(包括空气、脂肪、软组织和骨骼)的特性,对应的mi值为不同器官组织百分比的一个线性组合,定义为 miρ=∑j=0Nmijρj(4)miρ=∑j=0Νmijρj(4) 其中mij是对应灰度值ρj质量衰减系数的第j分量. 2.2 图形硬件干道的生成和模拟 DRR建立在下述两个任务基础上: Step1. 利用光线投射算法计算式(2); Step2. 根据式(3)模拟X-光线的负片曝光效果. 其中的关键在于Step1的计算,即光线投射.文献包括图1 a所示的三个基本步骤:经过建模、视点和透视变换,生成一系列视点光线;每根光线与体数据求交,按照从前至后或从后至前的顺序,沿着光线采样体数据;对每个采样点进行光照明计算并按顺序融合,获得对应像素的光亮度值.同时,如图1 b所示的图形硬件流水线也可简化为三步:在顶点渲染器中将三角形顶点首先从物体空间变换至视点空间,再经透视变换到屏幕空间;三角形光删化后的点在像素渲染器中进行纹理映射和渲染操作;最后的结果颜色值被输出到成像平面.光线投射和图形硬件绘制分别在图像空间和物体空间中进行,本质上而言,光线投射算法无法在图形硬件中实现.注意到,两者都可表示为三个阶段组成的流水线操作,类比图形硬件流水线中的三角形,当我们定义光线为一种绘制基本元素时,则可将光线投射算法映射到