CT图像重建加速的研究进展研究.pdfVIP

中国体视学学会cT理论与应用分会、中国工程物理研究院 2008年7月 应用电子学研究所、(CT理论与应用研究》编辑部论文集 CT图像重建加速的研究进展 邹永宁‘’6，谭辉～6，黄亮～。 (重庆大学a．光电技术及系统教育部重点实验室；b．ICT研究中心，重庆400044) 摘要：cT图像重建速度是衡量CT系统的重要指标之一。本文分析了cT系统常用的滤波反投影 算法的加速技术，重建的加速方法有很多，首先在普通Pc上主要有：几何参数表法，对称法， Intel SIMD指令优化法。其次通过构建集群系统可以进行并行重建，构建集群非常灵活，可以 根据应用情况使用不同的架构和节点数量，有刀片机，机架式并行机，也可以用普通PC机组建。 采用专门设计的硬件直接进行重建也可以极大地缩短重建时间，方法之一就是采用多个DSP按 一定的拓扑结构连接成并行计算模块，充分利用DSP的长流水线和并行能力，方法二是基于现 场可编程门阵列(FPGA)的方案，FPGA在顺序执行方面快，并且每一工作时钟周期内能够并行 完成多个逻辑运算，FPGA具有可重配置性，硬件方面表现出极大的灵活性。采用FP(；A实现cT 图像重建可利用其丰富的存储资源，使得投影数据的读取本地化，提高读取速度，解决存储器 带宽限制问题。目前GPU发展迅速，除了实现图形显示还广泛用于其它领域的通用计算，采用 GPU实现图像重建具有显著的优势，GPU是流处理器，具有一定的并行性，一般采用纹理映射的 方法实现反投影过程。IBM为分布计算开发了一种优化和新颖的通用目标架构，即细胞宽带引擎 5122图像，可以达到lO倍的加速比。使用双细胞系统(DcBS)和4个PS3节点集群进行锥柬重 建实验，两者每秒都可以处理多达30个5122的投影数据并以32位的浮点精度反投影到5123的 象素体上。 关键词；口；图像重建；加速技术 引言 计算机断层成像(Computed 域中都有着广泛的应用。CT图像重建是cT系统软件的重要组成之一，重建速度是cT系统 的重要指标之一。 长期以来，重建加速的研究日益受到工程界的重视，已有许多学者提出了各具特色的 方法。本文列举了用于图像重建加速的各种技术，并对它们进行了比较分析。CT图像重建 算法多种多样，但一般情况下采用的是滤波反投影算法(或称为卷积反投影算法，频域中 滤波等价于空域中卷积过程)，因此下面只讨论此类算法的加速技术。滤波反投影算法顾名 思义主要包括2个过程：滤波和反投影。除此之外的预处理过程计算耗时较少，可以不用 考虑加速。计算时间最长的是反投影，滤波次之，我们主要研究这两部分的加速。有些CT 扫描方式不同，比如扇束CT和锥束CT，相应地它们的滤波反投影算法也有些不一样。但是 当采用相同的加速技术实现重建算法时，其方法和效果是相似的，我们就不细致的区分了。 邹永宁等：cT图像重建加速的研究进展 277 1普通PC上的加速技术 1．1几何参数表法 在反投影过程中要频繁地使用三角函数，其输入变量由射线扫描时的几何位置确定， 故三角函数的值可以在反投影前预先确定。将所有计算出的值存入一张参数表中，反投影 时根据几何位置直接查参数表。几何参数表方法可明显地减少反投影时间，加速比至少在 10倍以上【I】。这种方法不会消耗太多内存，也不会与其他方法冲突，应尽量采用。 1．2对称法 利用反投影算法的对称性，只计算出部分图像数据并对称得到其它图像数据比1。对于扇 束CT扫描，具有一定关系的四幅投影的反投影运算可同时完成，从而减少最外层循环次数 到N／4(N是投影的幅数)，这样计算量仅为原来的1／4，达到优化的目的∞1。最后综合所 有重建结果生成一幅完整的图像。 1．3 Intel SIMD指令优化法 入了8个全新的128位)吲寄存器，能够一次执行4个紧缩单精度浮点数操作，同时引入 了m64与m128数据之间的转换指令，求倒数指令，求平方根指令等。这些指令非常适合滤 波和反投影这样浮点运算量大的应用场合H。。 Intel SIMD指令优化技术可以使图像的重建速度提高约5倍哺1。此技术适合细粒度的 并行运算，主要优化反投影和卷积循环算法的最内层代码