Premium
A separable quadratic surrogate total variation minimization algorithm for accelerating accurate CT reconstruction from few‐views and limited‐angle data
Author(s) -
Hou Xiaowen,
Teng Yueyang,
Kang Yan,
Qi Shouliang
Publication year - 2018
Publication title -
medical physics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.473
H-Index - 180
eISSN - 2473-4209
pISSN - 0094-2405
DOI - 10.1002/mp.12692
Subject(s) - algorithm , regularization (linguistics) , iterative reconstruction , mathematics , image quality , mathematical optimization , quadratic equation , norm (philosophy) , separable space , minification , optimization problem , computer science , artificial intelligence , image (mathematics) , mathematical analysis , geometry , political science , law
Purpose Total variation minimization (TVM) is a popular and useful method for accurate CT reconstruction from few‐views and limited‐angle data. However, the optimization procedure of previous TVM‐based algorithms is very time‐consuming. The purpose of this paper was to accelerate the high image quality CT reconstruction from few‐views and limited‐angle data. Method A new optimization algorithm based on the optimization transfer principle is proposed. The proposed algorithm uses TVM as the regularization term of the cost function that ensures the quality of the CT image reconstructed from few‐views and limited‐angle data. Additionally, half of the square of the difference between the original projection and the forward projection is used as the fidelity term. We then proved that the regularization term, i.e., 2‐norm TV, is a convex function. Based on the convexity of the regularization and fidelity terms, a separable quadratic surrogate (SQS) function was proposed to substitute the regularization and fidelity terms. The solution of the cost function can be obtained by minimizing the SQS function and building the next SQS function at the minimum point. Results Both numerical simulations and simulations using real experimental data showed that the proposed algorithm reconstruct high‐quality CT image from few‐views and limited‐angle data. The differences between the image reconstructed by the proposed algorithm and the images reconstructed by the previous algorithms are very small. However, the proposed algorithm required less than 1/10 time of the computational time of the previous algorithms. The image quality is not assessed in a rigorous way. Conclusion The proposed algorithm can greatly accelerate the accurate CT reconstruction from few‐views and limited‐angle data relative to previous algorithms.