激光与光电子学进展, 2017, 54 (2): 021006, 网络出版: 2017-02-10   

平滑约束的OSEM代数重建算法 下载: 548次

Smooth Constrained OSEM Iteration Reconstruction Algorithm
作者单位
西北核技术研究所, 陕西 西安 710024
摘要
有序子集期望最大化(OSEM)迭代算法是近年来发展较快的一种迭代类算法。但该算法在迭代过程中容易产生条纹状伪影、金属伪影或者散射伪影。构造了平滑约束矩阵作为先验信息引入到重建迭代过程, 建立了一种平滑约束OSEM(SC-OSEM)迭代重建算法。分别将中值滤波、全变差最小(TVM)方法作为平滑约束条件, 通过数值模拟, 针对不完备理想投影数据、含金属不完备投影数据、含噪声不完备投影数据三种情况, 重建出了与原始模型一致性较好的计算机层析成像技术(CT)图像, 比单独OSEM迭代算法重建质量高, 并且发现中值滤波约束重建图像的整体噪声较小, TVM算法使金属边界更清晰, 表明SC-OSEM迭代重建算法是一种精度高、适应性较强的CT重建算法。
Abstract
Order subset expectation maximization (OSEM) iterative algorithm is rapidly developed in recent years. But this algorithm is easy to generate some striation artifact, metal artifact or scattering artifact during iterative process. So a new iterative method named smooth constrained OSEM (SC-OSEM) is built, which utilizes smoothing constrained matrix as priori information in OSEM reconstruction. Median filtering algorithm and total variation minimization algorithm (TVM) are introduced as smoothing constrained conditions. The images reconstructed by computer tomography are consistent well with the initial model for the situations of incomplete projection data with ideal, metal, and noise, whose reconstruction qualities are better than only OSEM iterative algorithm. Median filtering constraint makes reconstructed image denoised, while TVM makes metal boundary clearer. These results manifest that SC-OSEM iterative algorithm is an adaptable and high precision CT reconstruction method.
参考文献

[1] 谭 海, 王大东, 薛艳玲, 等. 显微CT血管系统三维结构的骨架细化算法并行化设计实现[J]. 光学学报, 2015, 35(11): 1117003.

    Tan Hai, Wang Dadong, Xue Yanling, et al. Parallelization of 3D thinning algorithm for extracting skeleton of micro-CT vasculature[J]. Acta Optica Sinica, 2015, 35(11): 1117003.

[2] 刘 浩, 高万荣, 陈朝良. 手持式牙齿在体谱域光学相干层析成像系统研究[J]. 中国激光, 2016, 43(2): 0204003.

    Liu Hao, Gao Wanrong, Chen Chaoliang. Handheld spectral domain optical coherence tomography for in vivo tooth imaging[J]. Chinese J Lasers, 2016, 43(2): 0204003.

[3] Shepp L A, Vardi Y. Maximum likelihood reconstruction for emission tomography[J]. IEEE Transactions on Medical Imaging, 1982, 1(2): 113-122.

[4] Lange K, Carson R. EM reconstruction algorithms for emission and transmission tomography[J]. Journal of Computer Assisted Tomography, 1984, 8(2): 306-316.

[5] Hudson H M, Larkin R S. Accelerated image reconstruction using ordered subsets of projection data[J]. IEEE Transactions on Medical Imaging, 1994, 13(4): 601-609.

[6] 李 建, 吴志芳, 安洪振, 等. 松弛系数-有序子集-最大似然期望算法在60Co集装箱CT检测中的应用[J]. 核电子学与探测技术, 2007, 27(1): 45-49.

    Li Jian, Wu Zhifang, An Hongzhen, et al. A relaxed OSEM algorithm applied in 60Co container CT image reconstruction[J]. Nuclear Electronic & Detection Technology, 2007, 27(1): 45-49.

[7] 余 维. 不完备投影数据的CT重建算法研究[D]. 重庆: 重庆大学, 2014.

    Yu Wei. Investigation of reconstruction algorithms for incomplete CT projections[D]. Chongqing: Chongqing University, 2014.

[8] 阙介民, 王燕芳, 孙翠丽, 等. 基于不完备投影数据重建的四种迭代算法比较研究[J]. CT理论与应用研究, 2012, 21(2): 169-178.

    Que Jiemin, Wang Yanfang, Sun Cuili, et al. Comparison of four iterative algorithm based on incomplete projection reconstruction[J]. CT theory and Application, 2012, 21(2): 169-178.

[9] Ko S J, LeeY H. Center weighted median filters and their applications to image enhancement[J]. IEEE Transactions on Circuits and System, 1991, 38(9): 984-993.

[10] Liu Shaoqiang, Chen Lei, Fan Xiaoping, et al. Combining pseudo-median filter and median filter to improve performance[C]. IEEE International Conference on Computer Science and Information Technology, 2010, 4: 513-517.

[11] Rudin L, Osher S, Fatemi E. Nonlinear total variation based noise removal algorithms[J]. Physica D: Nonlinear Phenomena, 1992, 60(1-4): 259-268.

[12] Persson M, Bone D, Elmqvist H. Total variation norm for three-dimensional iterative reconstruction in limited vies angle tomography[J]. Physics in Medicine and Biology, 2001, 46(3): 853-866.

[13] Candes E, Romberg J, Tao T. Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information[J]. IEEE Trans Inf Theory, 2006, 52(2): 489-509.

[14] Sidky E Y, Kao C M, Pan X C. Accurate image reconstruction from few-views and limited-angle data in divergent beam CT[J]. X-Ray Sci Technol, 2006, 14(2): 119-139.

[15] Sidky E Y, Pan X C. Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization[J]. Phys Med Biol, 2008, 53(17): 4777-4807.

[16] 马继明, 张建奇, 宋顾周, 等. 全变分约束迭代滤波反投影CT重建[J]. 光学学报, 2015, 35(2): 0234002.

    Ma Jiming, Zhang Jianqi, Song Guzhou, et al. Total variation constrained iterative filtered backprojection CT reconstruction method[J]. Acta Optica Sinica, 2015, 35(2): 0234002.

张艳红, 盛亮, 李奎念, 李阳, 彭博栋, 张美, 赵吉祯, 袁媛. 平滑约束的OSEM代数重建算法[J]. 激光与光电子学进展, 2017, 54(2): 021006. Zhang Yanhong, Sheng Liang, Li Kuinian, Li Yang, Peng Bodong, Zhang Mei, Zhao Jizhen, Yuan Yuan. Smooth Constrained OSEM Iteration Reconstruction Algorithm[J]. Laser & Optoelectronics Progress, 2017, 54(2): 021006.

本文已被 1 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!