扇束重建中OSEM迭代算法与角度采样的模拟研究.pdfVIP

扇束重建中OSEM 迭代算法与角度采样的模拟研究 刘力 印胤 吴晓锋 （中国科学院高能物理研究所，北京 100039） 摘要：本文通过模拟生成扇束断层扫描成像中不同角度采样数目和不同角度采样范围的扇束 投影数据，对比研究了不同角度采集条件下图像重建OSEM （有序子集最大似然法）算法的 成像效果，验证了在扇束采集与重建中的规律，显示了OSEM 迭代重建算法相对于FBP （滤 波反投影法）的优越性，并探讨了如何控制角度采样来获取较好的重建效果，同时尽量减少 采集与处理的时间。 关键词：CT 图像重建 OSEM FBP 中图分类号： 在 CT (Computerized Tomography) 断层扫描与核医学 ECT （Emission Computerrized Tomography ）断层扫描中，扇束扫描是一种常用采集方式，其重 建方法就有解析算法和迭代算法。有序子集最大期望值方法（Ordered Subsets Expectation Maximization, OSEM ）是一种迭代重建算法，它具有较好的重建图像 质量和较短的计算时间，已应用在商业 ECT 的重建和CT 图像重建的研究中。 本文通过研究仿真 Phantom 模型数据分析比较了在不同角度采集条件下的扇束 OSEM 重建结果，并分别在无噪声和有噪声的情况下同FBP 重建图像进行对比， 显示了OSEM 迭代重建算法的特点和相对于FBP （滤波反投影法）的优越性。 1 CT 模型与OSEM 重建算法 本文模拟了一个扇束CT 原型，见图1： Y D FOV S γ max O X R dso 图1 模拟CT 原型示意图 其中被测物体范围 （FOV ）是一个以O 为圆心，半径为R 的圆（单层平面）， S 为点X 光源，D 为圆弧状等角间距排列的探测器 （bin ）阵列，以S 为其圆心； dso 为S-O 距离，γ 为扇束半张角，一个扇束覆盖整个FOV 。假设整个系统保 max 持物体不动，光源－探测器系统绕O 点进行360o 旋转采样。 我们选取图像大小为 128 ×128，以相邻两象素间的距离为长度单位， dso=247.3，γ 为15o ，D 中bin 的个数为128。 max 用i、j 标记象素，n 标记扇束的转角，m 标记一个扇束内不同张角的射线束 ray ）；断层密度分布记为f (i, j ) ，投影数据集合记为d (n, m) ，第(i，j )象素在扇 （ 束n 时对第m 个bin 的 （衰减）贡献因子记为G(i, j , n, m) ，若选用最大似然期望 法（Maximum Likelihood Expectation Maximization ，MLEM ）进行投影数据d 的 重建，其迭代公式可表示为[1] ： (k ) f (k −1) (i, j ) G(i, j , n, m)d (n, m) f (i, j ) ∑ (k −1) （1） ∑G(i, j , n, m) n ,m ∑G(i , j