基于多重网格的SIRT加速算法-CT理论与应用研究.PDF

第22卷 第 1 期 CT理论与应用研究 Vol.22, No.1 2013年1月（25-32） CT Theory and Applications Jan., 2013 乔灵博, 李亮, 陈志强. 基于多重网格的SIRT 加速算法[J]. CT 理论与应用研究, 2013, 22(1): 25-32. Qiao LB, Li L, Chen ZQ. Accelerated SIRT algorithm based on multigrid[J]. CT Theory and Applications, 2013, 22(1): 25-32. 基于多重网格的 SIRT 加速算法 a,b a,b a,b 乔灵博 ，李亮 ，陈志强 （清华大学 a.工程物理系；b.粒子技术与辐射成像教育部重点实验室，北京 100084） 摘要：在同时迭代重建算法，即 SIRT 的基础上，提出了一种基于多重网格的加速算法。SIRT 算法在对每个像素更新时，更新量与所有穿过该像素的射线所覆盖的其他像素有关。鉴于此， 通过运用多重网格的思想，对SIRT算法进行了分析，指出了其全局性的迭代特点，进而根据“细 网光滑、粗网校正”，实现较快的 SIRT 迭代重建。本文提出的新算法主要通过降低计算量来实 现加速效果,而与其他常见的加快迭代收敛速度的加速策略并不冲突。数值实验证明，该算法对 噪声抑制更好，同时在本文数值实验条件下，计算速度相比SIRT提高接近一倍。 关键词：SIRT；多重网格；加速 文章编号：1004-4140（2013）01-0025-08 中图分类号：TP301.6 文献标志码：A 常用的 CT（Computerized Tomography）重建算法一般可以分为解析重建算法和迭代重 建算法两类[1]。解析重建算法物理思想清晰、计算速度快，已经广泛应用于各类商业系统； 而迭代重建算法需要反复迭代，计算量庞大，执行速度慢一直成为限制其发展的重要因素。 随着近年来计算机硬件的快速发展以及各种加速技术的出现，迭代重建算法以其独特的优 势开始受到广泛关注。与解析重建算法相比，迭代重建算法对数据采集模式要求小，易于 加入各种约束条件和先验知识，适用面广泛。另一方面，得到相同质量的图像，迭代重建 算法需要更少的投影数据，这对减少受体的照射剂量有着重要意义。 本文关注迭代重建算法（Iterative Reconstruction Technique，IRT）中的同时迭代 重建算法（Simultaneous Iterative Reconstruction Technique，SIRT）。SIRT 是在代数 重建算法[2] （Algebraic Reconstruction Techniques，ART）的基础上提出的，旨在使重 建图像对测量误差及噪声不敏感[3]。SIRT 具有更好的稳定性，但是收敛速度慢、计算量 [4] [5] 较 ART 更大。1972 年，Gilbert 首次提出 SIRT 算法。1994 年，Hudson 等 通过将投影 数据分块提出了适合包括 EM（Expectation Maximization）、SIRT 等在内的有序子集加 速算法，可以使收敛速度提高一到两个数量级。Gregor 等[6]提出了一种基于本征值的近 似最佳松弛因子选择方法，采用该方法给出的近似最佳松弛因子，得到相同的图像质量 SIRT 的迭代步数减半。Guo 等[7]引入了 SIRT 前一次迭代结果对当前迭代数据的补偿，同 样达到了一定的加速效果。 注意到 SIRT全局性的迭代特点以及不同分辨率图像的相关性，本文提出了一种基于多 重网格 SIRT加速算法。通过运用“细网光滑、粗网校正”的策略，将迭代过程主要集中在 收稿日期：2012-08-30。 项目基金：国家自然科学基金；北京市自然科学基金资助项目。 26