基于参数估计的低剂量CT图像统计去噪

doi:10.3969/j.issn.1001-2400.2011.03.016

J4 ›› 2011, Vol. 38 ›› Issue (3): 99-106.doi: 10.3969/j.issn.1001-2400.2011.03.016

基于参数估计的低剂量CT图像统计去噪

张元科¹；张军英¹；卢虹冰²

(1. 西安电子科技大学计算机学院，陕西西安 710071；
2. 第四军医大学计算机应用教研室，陕西西安 710032)

收稿日期:2010-05-17 出版日期:2011-06-20 发布日期:2011-07-14
通讯作者: 张元科
作者简介:张元科(1979-)，男，讲师，西安电子科技大学博士研究生，E-mail: yuankezhang@163.com．
基金资助:
国家自然科学基金资助项目(61070137);国家自然科学基金重点资助项目(60933009);国家自然科学基金资助项目(30470490)

Parameter estimation and statistical noise reduction for low-dose CT sinogram

ZHANG Yuanke¹；ZHANG Junying¹；LU Hongbing²

(1. School of Computer Science and Technology, Xidian Univ., Xi'an 710071， China|
2. Dept. of Computer Application/BME, Fourth Military Medical Univ., Xi'an 710032， China)

Received:2010-05-17 Online:2011-06-20 Published:2011-07-14
Contact: ZHANG Yuanke

摘要/Abstract

摘要：

由于提高低剂量CT图像的信噪比是低剂量CT获得有效临床应用的关键，为此，提出了一种低剂量CT投影域的自适应统计降噪算法．针对低剂量CT投影图像的非平稳高斯噪声特性，采用EM算法自适应地估计图像模型中的参数，并在此基础上对图像进行最大后验概率估计，从而达到图像降噪的目的．在对参数的估计过程中，引入MCMC的吉布斯采样技术，并在算法中引入两项初始化技术，从而减少了参数估计过程中的计算量，加快了算法的收敛速度．对仿真投影数据以及真实投影数据的实验结果表明，与传统算法相比，该算法在抑制噪声及保持分辨率方面均具有明显优势．

关键词: 低剂量CT, 图像降噪, 参数估计, 最大后验概率估计

Abstract:

Improvement of the SNR of low-dose CT images is a crucial issue for the low-dose CT application. In this paper, we propose a novel adaptive statistical noise reduction algorithm for low-dose CT sinogram. The algorithm first adopts an EM algorithm to adaptively estimate the parameters of the image model based on the non-stationary Gaussian noise property in the low-dose CT projection data, and then uses the MAP estimation to restore the sinogram. In the parameters estimation procedure, a Gibbs sampler is used to handle the complicated computation problem. In addition, two initialization strategies are used in the algorithm to accelerate the convergence speed too. The effectiveness of the proposed algorithm is validated by both computer simulations and experimental studies. The advantage of the proposed approach over other methods is quantified by noise-resolution tradeoff curves.

Key words: low-dose CT, noise reduction, parameter estimation, maximum a posteriori estimation

中图分类号:

TP391.41

张元科；张军英；卢虹冰. 基于参数估计的低剂量CT图像统计去噪[J]. J4, 2011, 38(3): 99-106.

ZHANG Yuanke；ZHANG Junying；LU Hongbing. Parameter estimation and statistical noise reduction for low-dose CT sinogram[J]. J4, 2011, 38(3): 99-106.

参考文献

［1］ Yazdi M, Beaulieu L. Artifacts in Spiral X-ray CT Scanners: Problems and Solutions［J］. International Journal of Biological and Medical Sciences, 2009, 4(3): 135-139.
［2］ Hsieh J. Adaptive Streak Artifact Reduction in Computed Tomography Resulting from Excessive X-ray Photon Noise［J］. Med Phys, 1998(25): 2139-2147.
［3］ Lu Hongbing, Hsiao I, Li X, et al. Noise Properties of Low-Dose CT Projections and Noise Treatment by Scale Transformations ［C］//Nuclear Science Symposium Conference Record. IEEE Nuclear and Plasma Sciences Society. San Dieco: IEEE, 2001: 1662-1666.
［4］陈丽霞, 丁宣浩, 宋国乡. 基于总变分与小波变换的图像去噪算法［J］. 西安电子科技大学学报, 2008, 35(6): 1075-1079.
Chen Lixia, Ding Xuanhao, Song Guoxiang. Image Denoising Algorithm Based on Total Variation and Wavelet Transform［J］. Journal of Xidian University, 2008, 35(6): 1075-1079.
［5］卢孝强, 孙怡. 基于乘性正则化的有限角度CT重建算法［J］. 光学学报, 2010, 30(5): 1285-1290.
Lu Xiaoqiang, Sun Yi. Limited Angle Computed Tomography Teconstruction Algorithm Based on Multiplicative Regularization Method［J］. Acta Optica Sinica, 2010, 30(5): 1285-1290.
［6］ Duan Xinhui, Zhang Li, Xing Yuxiang, et al. Few-view Projection Reconstruction with an Iterative Reconstruction-reprojection Algorithm and TV Constraint［J］. IEEE Trans on Nucl Sci, 2009, 56(3): 1377-1382.
［7］ La Riviere P, Bian J, Vargas P A. Penalized-likelihood Sinogram Restoration for Computed Tomography［J］. IEEE Trans on Med Imag, 2006, 25(8): 1022-1036.
［8］ Li T, Li X, Wang J, et al. Nonlinear Sinogram Smoothing for Low-Dose X-ray CT［J］. IEEE Trans on Nucl Sci, 2004, 51(5): 2505-2513.
［9］ Wang J, Li T, Lu Hongbing, et al. Penalized Weighted Least-squares Approach to Sinogram Noise Reduction and Image Reconstruction for Low-Dose X-ray Computed Tomography ［J］. IEEE Trans on Med Imag, 2006, 25(10): 1272-1283.
［10］ Wang J, Lu Hongbing, Wen J, et al. Multiscale Penalized Weighted Least-squares Sinogram Restoration for Low-Dose X-Ray Computed Tomography［J］. IEEE Trans on Biomedical Engineering, 2008, 55(3): 1022-1031.
［11］ Lei T, Sewchand W. Statistical Approach to X-ray CT Imaging and Its Applications in Image Analysis—Part I: Statistical Analysis of X-ray CT Imaging ［J］. IEEE Trans on Med Imag,1992, 11(1): 53-61.
［12］ Geman S, Geman D. Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images ［J］. IEEE Trans on Pattern Analysis and Machine Intelligence, 1984(6): 721-741.
［13］ Saquib S, Bouman C A, Sauer K. ML Parameter Estimation for Markov Random Fields with Applications to Bayesian Tomography ［J］. IEEE Trans on Image Processing, 1998, 7(7): 1029-1044.
［14］ Dempster A P, Laird N M, Rubin D B. Maximum Likelihood from Incomplete Data via the E-M Algorithm［J］. J R Statist Soc B, 1977(39): 1-38.
［15］ Perez P. Markov Random Fields and Images［J］. CWI Quarterly, 1998, 11(4): 413-437.
［16］ Li S Z. Markov Random Field Modeling in Image Analysis［M］. 3rd ed. Berlin: Springer, 2009:194-198.
［17］ Besag J. On the Statistical Analysis of Dirty Pictures (with discussion)［J］. J Royal Statistical Soc: Series B, 1984, 48(3): 259-302.

[1]	荆福龙,张春杰,司伟建,王宇,焦淑红. 一种有效的CFCR域LFM信号分析方法[J]. 西安电子科技大学学报, 2019, 46(3): 102-108.
[2]	张俊林;刘明骞. 脉冲噪声衰落信道下MPSK信号的符号周期估计[J]. 西安电子科技大学学报, 2018, 45(5): 38-42+95.
[3]	陈世超;罗丰;雒梅逸香;胡冲. 采用梅林变换的帕累托分布参数估计方法[J]. 西安电子科技大学学报, 2018, 45(4): 18-22+117.
[4]	刘明骞;陈钱;吴芸;王江宏. Underlay认知下时频重叠信号的分量个数估计[J]. 西安电子科技大学学报, 2017, 44(6): 43-47+64.
[5]	胡冲;罗丰;张林让;陈帅霖;范一飞. 闭合表达式的Pareto分布参数估计范围扩展[J]. 西安电子科技大学学报, 2017, 44(2): 33-36+113.
[6]	韩勋;杜兰;刘宏伟. 窄带雷达观测下的锥体目标参数估计方法[J]. J4, 2015, 42(6): 43-48.
[7]	解虎;冯大政;虞泓波;袁明冬;聂卫科. 一种采用稀疏表示的快速空时自适应方法[J]. J4, 2015, 42(5): 55-62.
[8]	杨秋;张群;吴勇;罗迎;王敏. 海面运动舰船微动特征分析与参数估计方法[J]. J4, 2015, 42(5): 200-206.
[9]	王晓峰;张国毅;芮国胜. 低信噪比下多相码信号检测与参数估计方法[J]. J4, 2015, 42(1): 174-179.
[10]	段佳;邢孟道;张磊;王金伟;梁毅. 联合属性散射中心的极化目标重构新方法[J]. J4, 2014, 41(6): 18-24.
[11]	卢云龙;李明;闫琰. 利用调频率匹配的DRFM欺骗干扰检测方法[J]. J4, 2014, 41(5): 67-73+134.
[12]	段佳;张磊;邢孟道;梁毅. 合成孔径雷达目标特征提取新方法[J]. J4, 2014, 41(4): 13-19.
[13]	巩继玲;廖桂生;许京伟;朱圣棋. 一种参数估计的自适应辅助通道方法[J]. J4, 2014, 41(3): 8-13+78.
[14]	王凌;刘坚强;李国林;张真;杨继坤. 用旋转不变技术实现互耦和通道不一致的联合抑制[J]. J4, 2014, 41(1): 182-188.
[15]	张树银;郭英;齐子森;霍文俊. 共形阵列LFM信号多参数估计的传播算子算法[J]. J4, 2013, 40(4): 181-187.

基于参数估计的低剂量CT图像统计去噪

Parameter estimation and statistical noise reduction for low-dose CT sinogram

PDF (PC)

赞

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

Metrics

本文评价

推荐阅读 10