基于GPU的三维锥束CT图像的快速重建

基于GPU的三维锥束CT图像的快速重建

论文摘要

计算机断层成像技术(computed Tomography,简称为CT)已经广泛应用于人体组织成像、工业无损检测等领域。锥束CT也叫三维CT,由于在空间轴向分辨率、扫描速度以及射线利用效率等方面的优势,三维锥束CT逐渐得到了广泛的应用。锥束CT的重建算法大致分为迭代法和解析法两大类。解析法又分为精确重建和近似重建,近似算法重建图像的质量虽然没有精确重建的好,但是它具有执行效率高,机械运动简单等优点。Feldkamp、David和Kress提出了基于圆形扫描轨迹的实用近似重建算法(简称为FDK算法),后被推广到螺旋扫描轨迹,目前该算法是商用锥束CT机上通用的算法。相比解析重建法,迭代重建算法需要迭代运算,计算比较耗时,但是在某些特殊场合,例如投影数据不全或锥角很大,迭代重建算法具有明显的优势。迭代算法也是无损检测领域的研究热点,典型的迭代算法有SART、ML-EM、及OS-EM。三维图像重建运算量大,在很多情况下,重建所需时间难以满足实际应用的需求,CT图像重建的速度成了制约三维锥束CT普及的一个重要因素。为了满足这种实际要求,国内外发展出了很多图像重建的加速方法。这些加速方法可以分为两类:基于算法改进的加速方法、以及基于硬件的加速方法。硬件加速方法中又可分为:专用集成电路ASIC、FPGA门阵列、图形处理器GPU、以及PC集群并行计算等方法。基于GPU的3D快速CT重建算法,由于硬件发展迅速、实现成本低,实用性强,加速效果明显,已经成为一种重要的3D图像重建加速方法。为了充分利用GPU的并行处理能力以及提高数据传输效率,本文仔细分析了锥束CT扫描的FDK算法及GPU用于图像重建的特点,并在此基础上研究了一种结合使用GPU多重纹理(multitexture)加速三维圆周扫描锥束CT的FDK图像重建过程的方法。该方法采用多重纹理映射来提高反投影速度、减少中间数据存储量、减少浮点累加次数,使用顶点颜色通道来实现距离加权运算,采用扩展方法来增加并行反投影的纹理单元,从而提高重建速度。计算机实验结果表明,使用普通PC机重建尺寸为2563的图像,在保证数据精度为16bit浮点数的要求下,GPU反投影计算可以在10秒以内完成。与仅使用CPU的重建方法相比,GPU重建图像加速方法达到了较高的时间加速比。本文仔细分析了螺旋扫描轨迹下的FDK算法及GPU用于此算法时的特点,将用于圆周扫描轨迹下的GPU多重纹理加速FDK图像重建的方法,移植到螺旋扫描轨迹下的锥束CT图像的加速重建过程中,从而得到一种结合使用GPU多重纹理加速螺旋扫描轨迹下三维锥束CT的FDK图像重建过程的方法。计算机实验结果表明,使用普通PC机重建尺寸为2563的图像,旋转分度为512(两周)时,GPU反投影计算可以在16秒以内完成。与圆周扫描轨迹下的FDK图像重建时间(利用GPU进行加速)相比较,两种扫描轨迹下的GPU加速图像重建方法的加速效果是相当的。

论文目录

  • 中文摘要
  • 英文摘要
  • 1 绪论
  • 1.1 CT技术的发展与重建加速的意义
  • 1.1.1 CT技术的应用与发展
  • 1.1.2 CT重建加速的意义
  • 1.2 国内外研究现状
  • 1.2.1 图像重建的加速方法
  • 1.2.2 GPU加速图像重建的现状
  • 1.3 本文研究的目的和研究内容
  • 1.3.1 本文研究的目的
  • 1.3.2 本文研究的内容
  • 2 GPU加速CT重建的基础
  • 2.1 引言
  • 2.2 GPU用于CT重建的硬件基础
  • 2.2.1 GPU硬件的特点
  • 2.2.2 GPU实现投影计算的方法
  • 2.2.3 GPU实现反投影计算的方法
  • 2.3 二维扇形束投影重建算法
  • 2.3.1 投影与Radon变换
  • 2.3.2 傅立叶切片定理
  • 2.3.3 二维等间距扇形束投影重建算法
  • 2.4 锥束近似重建算法(FDK算法)
  • 2.4.1 锥束中心平面等间距扇束重建公式
  • 2.4.2 锥束重建公式
  • 2.5 锥束螺旋扫描轨迹的近似重建算法
  • 2.5.1 一般扫描轨迹下的锥束CT重建算法
  • 2.5.2 螺旋扫描轨迹下的锥束CT重建算法
  • 2.6 本章小结
  • 3 GPU加速FDK图像重建
  • 3.1 引言
  • 3.2 GPU实现FDK重建的步骤
  • 3.3 GPU反投影加速实现
  • 3.3.1 GPU多重纹理反投影方法
  • 3.3.2 距离加权系数实现方法
  • 3.3.3 GPU执行反投影的矩阵表示
  • 3.3.4 并行反投影纹理单元的扩展
  • 3.4 实验结果及分析
  • 3.4.1 GPU加速FDK图像重建的效果
  • 3.4.2 扩展并行纹理的加速效果
  • 3.5 GPU重建的误差分析
  • 3.5.1 GPU定点纹理图像重建的误差公式
  • 3.5.2 误差公式的一个应用
  • 3.6 本章小结
  • 4 GPU加速螺旋扫描轨迹FDK图像重建
  • 4.1 引言
  • 4.2 GPU实现螺旋轨迹FDK重建的步骤
  • 4.3 GPU加速螺旋扫描轨迹FDK重建的特点
  • 4.3.1 重建物体切片方式
  • 4.3.2 GPU反投影方式
  • 4.3.3 投影图像子集
  • 4.4 GPU反投影加速实现
  • 4.4.1 多重纹理映射
  • 4.4.2 距离加权计算
  • 4.4.3 并行纹理单元扩展
  • 4.5 实验结果及分析
  • 4.6 本章小结
  • 5 结论与展望
  • 5.1 研究工作总结
  • 5.2 存在问题及后续工作展望
  • 致谢
  • 参考文献
  • 附录
  • 相关论文文献

    标签:;  ;  ;  ;  ;  

    基于GPU的三维锥束CT图像的快速重建
    下载Doc文档

    猜你喜欢