基于可变阶变分模型的医用低剂量CT图像去噪

doi:10.13700/j.bh.1001-5965.2018.0775

北京航空航天大学学报 ›› 2019, Vol. 45 ›› Issue (9): 1757-1764.doi: 10.13700/j.bh.1001-5965.2018.0775

基于可变阶变分模型的医用低剂量CT图像去噪

王娜^1,2, 张权^1,2, 刘祎^1,2, 贾丽娜^1,2, 桂志国^1,2

1. 中北大学生物医学成像与影像大数据山西省重点实验室, 太原 030051;
2. 中北大学信息与通信工程学院, 太原 030051

收稿日期:2019-01-02 出版日期:2019-09-20 发布日期:2019-09-29
通讯作者: 桂志国 E-mail:gzgtg@163.com
作者简介:王娜,女,博士研究生。主要研究方向:图像处理;桂志国,男,博士,教授,博士生导师。主要研究方向:信号与信息处理、图像处理和识别、图像重建。
基金资助:
国家自然科学基金（61671413，61801438）；国家重大科学仪器设备开发专项（2014YQ24044508）；山西省归国学者基金（2016-085）

Medical low-dose CT image denoising based on variable order variational model

WANG Na^1,2, ZHANG Quan^1,2, LIU Yi^1,2, JIA Lina^1,2, GUI Zhiguo^1,2

1. Shanxi Provincial Key Laboratory for Biomedical Imaging and Big Data, North University of China, Taiyuan 030051, China;
2. School of Information and Communication Engineering, North University of China, Taiyuan 030051, China

Received:2019-01-02 Online:2019-09-20 Published:2019-09-29
Supported by:
National Natural Science Foundation of China (61671413,61801438); National Key Scientific Instrument and Equipment Development Project of China (2014YQ24044508); Shanxi Provincial Foundation for Returned Scholars (Main Program), China (2016-085)

摘要/Abstract

摘要： 为了降低患者的辐射风险，低剂量CT（LDCT）广泛用于临床诊断，但辐射剂量的减少在重建的LDCT图像中引入了斑点噪声和条纹伪影。为了提高LDCT图像的质量，提出了一种基于可变阶变分模型的后处理技术。所提出的变分模型使用边缘指示器控制变分阶数，根据图像的特征在一阶全变分（TV）正则项和二阶有界Hessian（BH）正则项之间交替变换。采用基于快速傅里叶变换（FFT）的分裂Bregman算法求解所提出的变分模型。该模型在保留高剂量CT（HDCT）图像相应结构的同时，有效抑制了斑点噪声和条纹伪影。重建的图像和实验数据表明，所提出的变分模型比现有的先进模型具有更好的质量。

关键词: 低剂量CT (LDCT), 图像降噪, 边缘指示器, 全变分(TV), 有界Hessian (BH), 快速傅里叶变换(FFT), 分裂Bregman算法

Abstract: Low-dose CT (LDCT) is widely used for clinical diagnosis to reduce radiation risk to patients. However, the radiation dose reduction introduces mottle noise and streak artifacts into the reconstructed LDCT images. In this paper, a post-processing technique is proposed based on variable order variational model to improve the LDCT image quality. The proposed variational model employs the edge indicator to control the order of variation, which can alternate between the first order total variation (TV) regularizer and second order bounded Hessian(BH) regularizer based on the image feature. Moreover, the proposed model is solved by split Bregman algorithm based on fast Fourier transform (FFT). The proposed model effectively suppresses mottle noise and streak artifacts, meanwhile preserving structure in reference to high-dose CT (HDCT) images. The reconstructed images and experimental data indicate that the proposed model has better quality than some existing state-of-the-art models.

Key words: low-dose CT (LDCT), image denoising, edge indicator, total variation (TV), bounded Hessian (BH), fast Fourier transform (FFT), split Bregman algorithm

中图分类号:

TP391

王娜, 张权, 刘祎, 贾丽娜, 桂志国. 基于可变阶变分模型的医用低剂量CT图像去噪[J]. 北京航空航天大学学报, 2019, 45(9): 1757-1764.

WANG Na, ZHANG Quan, LIU Yi, JIA Lina, GUI Zhiguo. Medical low-dose CT image denoising based on variable order variational model[J]. JOURNAL OF BEIJING UNIVERSITY OF AERONAUTICS AND A, 2019, 45(9): 1757-1764.

参考文献

[1] ZHU Y,ZHAO M,ZHAO Y,et al.Noise reduction with low dose CT data based on a modified ROF model[J].Optics Express,2012,20(16):17987-18004.
[2] CHEN Y,YIN X,SHI L,et al.Improving abdomen tumor low-dose CT images using a fast dictionary learning based processing[J].Physics in Medicine and Biology,2013,58(16):5803-5820.
[3] ZHANG C,ZHANG T,LI M,et al.Low-dose CT reconstruction via L₁ dictionary learning regularization using iteratively reweighted least-squares[J].BioMedical Engineering OnLine,2016,15(1):66.
[4] LEE D,LEE J,KIM H,et al.A feasibility study of low-dose single-scan dual-energy cone-beam CT in many-view under-sampling framework[J].IEEE Transactions on Medical Imaging,2017,36(12):2578-2587.
[5] CHEN Y,LIU J,HU Y,et al.Discriminative feature representation:An effective postprocessing solution to low dose CT imaging[J].Physics in Medicine and Biology,2017,62(6):2103-2131.
[6] CHEN Y,LIU J,XIE L,et al.Discriminative prior-prior image constrained compressed sensing reconstruction for low-dose CT imaging[J].Scientific Reports,2017,7(1):13868.
[7] FRUSH D P,DONNELLY L F,ROSEN N S.Computed tomography and radiation risks:What pediatric health care providers should know[J].Pediatrics,2003,112(4):951-957.
[8] BRENNER D J,HALL E J.Computed tomography-An increasing source of radiation exposure[J].New England Journal of Medicine,2007,357(22):2277-2284.
[9] CHEN Y,SHI L,YANG J,et al.Radiation dose reduction with dictionary learning based processing for head CT[J].Australasian Physical and Engineering Sciences in Medicine,2014,37(3):483-493.
[10] YANG Q,YAN P,ZHANG Y,et al.Low dose CT image denoising using a generative adversarial network with wasserstein distance and perceptual loss[J].IEEE Transactions on Medical Imaging,2018,37(6):1348-1357.
[11] LIU J,MA J,ZHANG Y,et al.Discriminative feature representation to improve projection data inconsistency for low dose CT imaging[J].IEEE Transactions on Medical Imaging,2017,36(12):2499-2509.
[12] HASAN A M,MELLI A,WAHID K A,et al.Denoising low-dose CT images using multi-frame blind source separation and block matching filter[J].IEEE Transactions on Radiation and Plasma Medical Sciences,2018,2(4):279-287.
[13] DIWAKAR M,KUMAR M.CT image denoising using NLM and correlation-based wavelet packet thresholding[J].IET Image Processing,2018,12(5):708-715.
[14] YOU C,YANG Q,SHAN H,et al.Structure-sensitive multi-scale deep neural network for low-dose CT denoising[J].IEEE Access,2018,6:41839-41855.
[15] LIU Y,ZHANG Y.Low-dose CT restoration via stacked sparse denoising autoencoders[J].Neurocomputing,2018,284:80-89.
[16] 罗立民,胡轶宁,陈阳.低剂量CT成像的研究现状与展望[J].数据采集与处理,2015,30(1):24-34.LUO L M,HU Y N,CHEN Y.Research status and prospect for low-dose CT imaging[J].Data Acquisition and Processing,2015,30(1):24-34(in Chinese).
[17] LU W,DUAN J,QIU Z,et al.Implementation of high-order variational models made easy for image processing[J].Mathematical Methods in the Applied Sciences,2016,39(14):4208-4233.
[18] DUAN J,QIU Z,LU W,et al.An edge-weighted second order variational model for image decomposition[J].Digital Signal Processing,2016,49:162-181.
[19] DUAN J,WARD W O C,SIBBETT L,et al.Introducing anisotropic tensor to high order variational model for image restoration[J].Digital Signal Processing,2017,69:323-336.
[20] CHEN Y,YANG Z,HU Y,et al.Thoracic low-dose CT image processing using an artifact suppressed large-scale nonlocal means[J].Physics in Medicine and Biology,2012,57(9):2667-2688.
[21] BUADES A,COLL B,MOREL J M.A non-local algorithm for image denoising[C]//IEEE Computer Society Conference on Computer Vision and Pattern Recognition.Piscataway,NJ:IEEE Press,2005,2:60-65.
[22] WANG J,LU H,WEN J,et al.Multiscale penalized weighted least-squares sinogram restoration for low-dose X-ray computed tomography[J].IEEE Transactions on Biomedical Engineering,2008,55(3):1022-1031.

基于可变阶变分模型的医用低剂量CT图像去噪

Medical low-dose CT image denoising based on variable order variational model

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 6

编辑推荐

Metrics

本文评价

[1]	孔志杰, 郝新红, 栗苹, 闫晓鹏. 调频引信谐波时序检测抗干扰方法及实现[J]. 北京航空航天大学学报, 2018, 44(3): 549-555.
[2]	冉达, 尹灿斌, 贾鑫. 频率步进SAR虚拟阵列模型成像方法[J]. 北京航空航天大学学报, 2018, 44(1): 82-88.
[3]	刘晓玉, 王之栎, 丁蕾, 王博文, 刘立坤. 反转轴间双端面气膜密封性能分析[J]. 北京航空航天大学学报, 2017, 43(3): 608-614.
[4]	谭浩, 许传龙, 张彪, 王式民. 基于近场散射的颗粒粒径分布测量[J]. 北京航空航天大学学报, 2017, 43(2): 381-387.
[5]	李龙军, 王布宏, 夏春和, 刘新波, 曹帅. 基于子阵激励能量匹配的多子阵交错阵列设计[J]. 北京航空航天大学学报, 2016, 42(11): 2395-2402.
[6]	王钢, 杨民. X射线动态数字图像降噪方法与快速实现[J]. 北京航空航天大学学报, 2010, 36(6): 741-744.