基于全变分最小化和快速一阶方法的低剂量ct有序子集-物理学报.pdfVIP

物理学报 Acta Phys. Sin. Vol. 63, No. 13 (2014) 138701 基于全变分最小化和快速一阶方法的低剂量CT 有序子集图像重建∗ 毛宝林 陈晓朝 孝大宇 范晟昱 滕月阳 康雁 1)(东北大学中荷生物医学与信息工程学院, 沈阳 110819) 2)(东北大学医学影像计算教育部重点实验室, 沈阳 110819) ( 2014 年3 月31 日收到; 2014 年5 月5 日收到修改稿) 低剂量计算机断层成像(computed tomography, CT) 具有减少X 射线对患者的伤害的优势. 本文主要针 对从不完备投影数据重建出高质量低剂量CT 图像的问题. 通常, 这个问题可以通过统计图像重建方法来实 现, 而统计重建算法需要非常多的迭代次数, 导致了巨大的计算时间压力, 以至于很难应用在实践中. 为解决 此问题, 本文提出一种有序子集重建算法, 该算法结合了全变分最小化和快速一阶方法以减少重建的迭代次 数, 采用Split Bregman 交替方向法求解上述优化问题, 利用投影到凸集合的方法加快迭代的收敛速率. 实验 结果表明, 在同样的迭代次数下, 本文提出的方法与基于有序子集的一阶方法相比较, 相对重建误差的下降速 度更快. 关键词: 全变分, 投影到凸集合, 交替方向法, 有序子集 PACS: 87.59.–e, 42.30.Wb, 81.70.Tx, 07.85.–m DOI: 10.7498/aps.63.138701 感知(compressed sensing, CS) 理论. CS 理论利用 1 引 言 信号在变换域稀疏的性质, 突破了奈奎斯特-香农 采样定理, 可以用较低的频率采集信号并重构出非 由于CT 的广泛使用, 导致患者吸收的X 射线 常好的原始信号. 受CS 理论的启发, 基于稀疏表 辐射剂量越来越高, 这对身体健康造成极大的隐 达或稀疏约束的CT 重建算法, 可以利用不完备投 1 5 患 . 在保证成像质量的前提下, 减少放射剂量具 影数据重建出高质量的CT 图像 . 2014 年, 杨富 有非常重要的现实意义和应用价值. 低剂量CT 通 6 强等 对CT 不完备投影数据重建算法做了综述 常采用脉冲调制的方式控制球管产生X 射线, 同时 性研究, 其中对压缩感知理论在CT 重建中的应用 减少探测器采集的投影数据, 以降低患者吸收的辐 进行了详细的讨论. 射剂量. 在减少探测器采集投影数的扫描方式下, 通常, 全变分(total variation, TV) 79 和 成像问题变为一个病态反问题. 滤波反投影方法重 范数是CT 重建中图像稀疏约束的两种常用方式. 建出的图像有伪影严重、视觉质量差等问题, 无法 对于带约束的稀疏最小化模型的求解, 出现很多 满足诊断的基本要求. 对于这种待求解的图像数据 新的方法, 最具典型的两类方法是交替方向法和一 远大于已获取的投影数据的病态反问题, 迭代重建 阶方法, 其中交替方向法可以分为乘子交替方向法 方法