计算机应用 ›› 2014, Vol. 34 ›› Issue (12): 3605-3608.

• 虚拟现实与数字媒体 • 上一篇    下一篇

基于变指数各向异性扩散和非局部的最大似然期望最大低剂量CT重建算法

张芳1,崔学英1,张权1,董婵婵2,刘祎1,孙未雅2,白云蛟1,桂志国1,3   

  1. 1. 电子测试技术国家重点实验室(中北大学),太原 030051
    2. 电子测试技术国家重点实验室(中北大学),太原 030051;
    3. 仪器科学与动态测试教育部重点实验室(中北大学),太原 030051
  • 收稿日期:2014-07-07 修回日期:2014-08-28 出版日期:2014-12-01 发布日期:2014-12-31
  • 通讯作者: 桂志国
  • 作者简介:张芳(1989-),女,山西朔州人,硕士研究生,主要研究方向:基于低剂量CT的图像重建;崔学英(1978-),女,山西临汾人,讲师,博士研究生,主要研究方向:图像处理、微分方程应用;张权(1974-),男,山西大同人,副教授,博士研究生,主要研究方向:图像处理与科学可视化;董婵婵(1988-),女,山东金乡人,硕士研究生,主要研究方向:基于偏微分的图像降噪;刘祎(1987-),女,河南睢县人,博士研究生,主要研究方向:图像处理、医学图像重建;孙未雅(1991-),女,河北遵化人,硕士研究生,主要研究方向:基于低剂量CT的图像重建;白云蛟(1990-),男,山西五台人,硕士研究生,主要研究方向:图像处理与重建;桂志国(1972-),男,河北蓟县人,教授,博士,主要研究方向:信号与信息处理、图像处理和识别、图像重建。
  • 基金资助:

    国家自然科学基金资助项目;山西省国际合作项目;国家自然科学基金资助项目;山西省研究生优秀创新项目;中北大学第十届研究生科技基金资助项目;山西省高等学校优秀青年学术带头人支持计划资助项目;中北大学2013年校科学基金资助项目

MLEM low-dose CT reconstruction algorithm based on variable exponent anisotropic diffusion and non-locality

ZHANG Fang1,CUI Xueying1,ZHANG Quan1,DONG Chanchan1, 1,SUN Weiya1,BAI Yunjiao1,GUI Zhiguo1,2   

  1. 1. National Key Laboratory for Electronic Measurement Technology (North University of China), Taiyuan Shanxi 030051, China;
    2. Key Laboratory of Instrumentation Science and Dynamic Measurement of Ministry of Education (North University of China), Taiyuan Shanxi 030051, China
  • Received:2014-07-07 Revised:2014-08-28 Online:2014-12-01 Published:2014-12-31
  • Contact: GUI Zhiguo

摘要:

针对低剂量计算机断层扫描(CT)重建图像发生严重衰退的问题,提出一种基于变指数和非局部的最大似然期望最大(MLEM)低剂量CT重建算法。该算法考虑了传统各向异性扩散中降噪不充分的缺点,把可以有效折中热传导和各向异性扩散(P-M)这两种模型的变指数,以及代替梯度检测边缘和细节的相似度函数运用到传统各向异性扩散中,从而达到所期望的效果。该算法在每次迭代中首先采用基本的MLEM算法对低剂量CT投影数据进行重建; 然后利用基于非局部的相似性测度以及变指数和模糊数学的理论对各向异性扩散的扩散函数进行改进,用改进后的各向异性扩散对重建图像进行降噪;最后使用中值滤波对图像进行处理从而消除脉冲噪声点。实验结果表明,所提出算法的均方绝对误差、归一化均方距离均比有序子集惩罚最小二乘(OS-PLS)、有序子集惩罚最大似然一步迟疑(OS-PML-OSL)、基于传统P-M、基于方差的算法小,获得了高达10.52的信噪比。该算法重建出的图像可以在有效消除噪声的同时较好地保持图像的边缘和细节信息。

Abstract:

Concerning the serious recession problems of the low-dose Computed Tomography (CT) reconstruction images, a low-dose CT reconstruction method of MLEM based on non-locality and variable exponent was presented. Considering the traditional anisotropic diffusion noise reduction is insufficient, variable exponent which could effectively compromise between heat conduction and anisotropic diffusion P-M models, and the similarity function which could detect the edge and details instead of gradient were applied to the traditional anisotropic diffusion, so as to achieve the desired effect. In each iteration, firstly, the basic MLEM algorithm was used to reconstruct the low-dose projection data. And then the diffusion function was improved by the non-local similarity measure, variable index and fuzzy mathematics theory, and the improved anisotropic diffusion was used to denoise the reconstructed image. Finally median filter was used to eliminate impulse noise points in the image. The experimental results show the proposed algorithm has a smaller numerical value than OS-PLS (Ordered Subsets-Penalized Least Squares), OS-PML-OSL (Ordered Subsets-Penalized Maximum Likelihood-One Step Late), and the algorithm based on the traditional PM, in the variance of Mean Absolute Error (MAE), and Normalized Mean Square Distance (NMSD), especially its Signal-to-Noise Ratio (SNR) is up to 10.52. This algorithm can effectively eliminate the bar of artifacts, and can keep image edges and details information better.

中图分类号: