计算机应用 ›› 2020, Vol. 40 ›› Issue (6): 1793-1798.DOI: 10.11772/j.issn.1001-9081.2019111955

• 虚拟现实与多媒体计算 • 上一篇    下一篇

基于Chambolle-Pock算法框架的高阶TV图像重建算法

席雅睿, 乔志伟, 温静, 张艳娇, 杨雯晶, 闫慧文   

  1. 山西大学 计算机与信息技术学院,太原030006
  • 收稿日期:2019-11-18 修回日期:2020-01-07 出版日期:2020-06-10 发布日期:2020-06-18
  • 通讯作者: 乔志伟(1977—)
  • 作者简介:席雅睿(1993—),女,山西临汾人,硕士研究生,主要研究方向:医学图像重建、图像处理.乔志伟(1977-), 男,山西临汾人,教授,博士,主要研究方向:医学图像重建、信号处理、高性能计算.温静(1995—),女,山西介休人,硕士研究生,主要研究方向:医学图像处理、深度学习.张艳娇(1995—),女,山西长治人,硕士研究生,主要研究方向:图像处理、深度学习.杨雯晶(1995—),女,山西晋城人,硕士研究生,主要研究方向:医学图像处理、深度学习.闫慧文(1996—),女,山西晋城人,硕士研究生,主要研究方向:医学图像重建、图像处理.
  • 基金资助:
    山西省重点研发计划项目(201803D421012);山西省留学人员科技活动项目(RSC1622)。

High order TV image reconstruction algorithm based on Chambolle-Pock algorithm framework

XI Yarui, QIAO Zhiwei, WEN Jing, ZHANG Yanjiao, YANG Wenjing, YAN Huiwen   

  1. School of Computer and Information Technology, Shanxi University, Taiyuan Shanxi 030006
  • Received:2019-11-18 Revised:2020-01-07 Online:2020-06-10 Published:2020-06-18
  • Contact: QIAO Zhiwei, born in 1977, Ph. D., professor. His research interests include medical image reconstruction, signal processing, high performance computing.
  • About author:XI Yarui, born in 1993, M. S. candidate. Her research interests include medical image reconstruction, image processing.QIAO Zhiwei, born in 1977, Ph. D., professor. His research interests include medical image reconstruction, signal processing, high performance computing.WEN Jing, born in 1995, M. S. candidate. Her research interests include medical image processing, deep learning.ZHANG Yanjiao, born in 1995,M. S. candidate. Her research interests include image processing, deep learning.YANG Wenjing, born in 1995, M. S. candidate. Her research interests include medical image processing, deep learning.YAN Huiwen, born in 1996, M. S. candidate. Her research interests include medical image reconstruction, image processing.
  • Supported by:
    Shanxi Provincial Key Research and Development Plan (201803D421012), the Returnees Science and Technology Activities Project of Shanxi Province (RSC1622).

摘要: 传统的总变差(TV)最小算法是一种基于压缩感知(CS)的经典迭代重建算法,可以从稀疏数据或含噪数据中高精度地重建图像。然而,TV算法在重建分段常数特征不明显的图像时可能会引入块状伪影,通过研究得出,在图像去噪中使用高阶总变差(HOTV)能有效压制TV模型引入的块状伪影。鉴于此,提出了一种HOTV图像重建模型及其Chambolle-Pock(CP)求解算法。具体来说,以二阶梯度构建二阶TV范数,进而设计了一种数据保真约束的二阶TV最小重建模型,并推导出了相应的CP算法。在理想数据投影和含噪数据投影条件下,分别采用基于波浪背景的Shepp-Logan模体、灰度渐变模体以及真实CT图像模体进行重建实验,并进行定性和定量分析。理想数据投影的重建结果表明,和传统TV算法相比,HOTV算法能有效压制块状伪影并提高重建精度。含噪数据投影的重建结果表明,HOTV算法和TV算法均有良好的抗噪能力,但HOTV算法的保边性能更好且抗噪性更强。在重建分段常数特征不明显而灰度波动特征明显的图像时,HOTV算法是一种比TV算法更优的重建算法。所提HOTV算法可以被推广到各种扫描模式下的CT重建及其他成像模态中。

关键词: 高阶总变差, 约束优化, 压缩感知, 图像重建, Chambolle-Pock算法

Abstract: The traditional Total Variation (TV) minimization algorithm is a classical iterative reconstruction algorithm based on Compressed Sensing (CS), and can accurately reconstruct images from sparse and noisy data. However, the block artifacts may be brought by the algorithm during the reconstruction of image having not obvious piecewise constant feature. Researches show that the use of High Order Total Variation (HOTV) in the image denoising can effectively suppress the block artifacts brought by the TV model. Therefore, a HOTV image reconstruction model and its Chambolle-Pock (CP) solving algorithm were proposed. Specifically, the second order TV norm was constructed by using the second order gradient, then a data fidelity constrained second order TV minimization model was designed, and the corresponding CP algorithm was derived. The Shepp-Logan phantom in wave background, grayscale gradual changing phantom and real CT phantom were used to perform the image reconstruction experiments and qualitative and quantitative analysis under ideal data projection and noisy data projection conditions. The reconstruction results of ideal data projection show that compared to the traditional TV algorithm, the HOTV algorithm can effectively suppress the block artifacts and improve the reconstruction accuracy. The reconstruction results of noisy data projection show that both the traditional TV algorithm and the HOTV algorithm have good denoising effect but the HOTV algorithm is able to protect the image edge information better and has higher anti-noise performance. The HOTV algorithm is a better reconstruction algorithm than the TV algorithm in the reconstruction of image having not obvious piecewise constant feature and obvious grayscale fluctuation feature. The proposed HOTV algorithm can be extended to CT reconstruction under different scanning modes and other imaging modalities.

Key words: High Order Total Variation (HOTV), constrained optimization, Compressed Sensing (CS), image reconstruction, Chambolle-Pock (CP) algorithm

中图分类号: