ct图像重建加速的几种方法-计算机系统应用.pdf

2009 年 第 4 期 计 算 机 系 统 应 用 ① CT 图像重建加速的几种方法 Research Progress on CT Image Reconstruction Acceleration 邹永宁 谭 辉 黄 亮 (重庆大学 ICT 研究中心 重庆 400044) 摘 要： CT 图像重建速度是衡量CT 系统的重要指标之一。论文针对滤波反投影算法综述了当前国内外的重建 加速技术，阐述了普通PC 上的加速方法，通过构建集群系统进行并行重建的方法，采用多个DSP 连 接成并行计算模块的加速方法，基于现场可编程门阵列(FPGA)的并行重建方案，采用 GPU 图像卡实 现图像重建，采用细胞宽带引擎(CBE)的图像重建方案。在应用性能方面分析了各种加速方法的优缺点。 关键词 ： 计算机断层成像 图像重建 加速技术 计算机断层成像(Computed Tomography, CT) 调用三角函数将非常耗时。因此可以在反投影前先行 作为一种现代成像技术在医疗和工业领域中都有着广 确定三角函数的值，将所有计算出的值存入一张参数 泛的应用。CT 图像重建是 CT 系统软件的重要组成之 表中，反投影时根据几何位置直接查参数表。采用几 一，重建速度是 CT 系统的重要指标之一。长期以来， 何参数表方法可明显地减少反投影时间，加速比至少 许多学者提出了各具特色的重建加速方法。本文列举 在 10 倍以上。这种方法不会消耗太多内存，也不会 了用于图像重建加速的各种技术，并对它们进行了比 与其他方法冲突，在实际应用中基本上都要采用。 较分析。 1.2 对称法 CT 图像重建算法多种多样，但通常情况下采用最 利用反投影算法的对称性，只计算出部分图像数 [3] 多的是滤波反投影算法，因此下面只讨论此类算法的 据并对称得到其它图像数据 。对于扇束 CT 扫描，具 加速技术。 有一定关系的四幅投影的反投影运算可同时完成，从 滤波反投影是把某一视角下测得的投影数据做一 而减少最外层循环次数到 N ／4(N 是投影的幅数) ，这 维傅里叶变换，对上述结果乘以滤波函数H (  )  ， 样计算量仅为原来的四分之一，达到优化的目的[4]。 再做傅里叶反变换，得到修正的投影数据，然后进行 最后综合所有重建结果生成一幅完整的图像。 反投影。将各个视角下的反投影值进行累加，得重建 1.3 Intel S