可变约束OSEM图像重建算法仿真研究.pdfVIP

可变约束 OS-EM 图像重建算法仿真研究 吴朝霞 1 刘力 2 柴新禹 3 1 清华大学工程物理系，北京，100084 2 中国科学院高能物理研究所，北京，100856 3 上海交通大学生物医学工程系，上海，200030 摘要：有序子集最大期望值方法 （Ordered Subsets Expectation Maximization ，OS-EM ）具有 较高的重建图像质量和较短的计算时间，正逐步应用在正电子发射断层成像仪 （Position Emission Tomography，PET ）图像重建过程中。论文提出了一种改进的OS-EM 图像重建算 法——可变约束 OS-EM 迭代算法 （Variable Constraint OS-EM，VCOSEM ），对仿真 Phantom 模型在不同测量条件以及不同子集划分情况下的研究结果表明，该方法具有分辨率高，噪声 低的优点，能够提高重建图像质量。 关键词：可变约束 OS-EM PET 1．前言 正电子发射断层扫描仪（Positron Emission Tomography ，PET ）是当今最高层次的核医 学技术，也是当前医学界公认的最先进的大型医疗诊断成像设备，PET 是在分子水平上利用 影像技术研究人体心脑代谢和受体功能的唯一手段，是具有较高特异性的功能显像和分子显 像，除显示形态结构外，能够在分子水平上提供有关脏器及其病变的功能信息。PET 在肿瘤 学，心血管疾病学和神经系统疾病学研究中，以及新医药学开发研究等领域中已经显示出它 卓越的性能。图像重建方法是正电子发射断层成像扫描仪中的一个关键技术，随着先进的核 医学断层影像设备的广泛应用和计算机技术的迅速发展，图像重建方法的研究越来越受到人 们的重视。 传统的图像重建方法主要分为解析法和迭代法。解析法的优点是速度快，可用于临床实 时断层重建。但当测量噪声较大或采样不充分时，这类算法的成像效果不甚理想，且难以在 重建中引入各种校正和约束。迭代法最大优点之一是可以根据具体成像条件引入与空间几何 有关的或与测量值大小有关的约束和条件因子，但是迭代法收敛速度慢，运算时间长，运算 量大，而且重建图像会随着迭代次数的增加而趋于 “老化”甚至发散，出现 “硬边”现象和 “checkerboard ”效应等高频伪影，这些缺点极大地限制了它在临床中的应用[1,2] 。近年来一 种具有较高图像质量和较短计算时间的重建算法——有序子集最大期望值方法 （Ordered Subsets Expectation Maximization ，OS-EM ）正逐步应用在临床PET 重建技术中。本文提出 了一种改进的OS-EM 图像重建算法——可变约束OS-EM 迭代算法 （VCOSEM ），目的是在 加快收敛的同时，进一步抑制噪声，改善图像质量，并对仿真 Phantom 模型数据进行了深 入细致的研究。 2 ．可变约束OS-EM 迭代算法 有 序 子 集 最 大 期 望 值 法 是 在 最 大 似 然 期 望 法 （ Maximum Likelihood Expectation-maximization ，ML-EM ）的基础上发展起来的一种图像迭代重建算法。 传统的 ML-EM 方法计算式为[3] ： (k ) f (k −1) (i, j ) p (i, j , n, m)d (n, m) f (i, j ) ∑p (i, j , n, m) ∑ p (i , j , n, m)f (k −1) (i , j ) ( 1) n,m ∑ n ,m i,j ( ) ( ) ( ) 其中 i, j 为像素坐标，f i, j 为该射线所代表的放射性（强）浓度大小， n, m 为探测器 空间坐标， ( )