Affiliations: [a] School of Mathematical Sciences, Chongqing Normal University, ChongQing, China
| [b] School of Biomedical Engineering and Imaging, Xianning Medical College, Hubei University of Science and Technology, Xianning, China
| [c] Key Laboratory of Optoeletronic and Intelligent Control, Hubei University of Science and Technology, Xianning, China
| [d] College of Computer Science, Chengdu University, Chengdu, China
Correspondence:
[*]
Corresponding authors. Chengxiang Wang, School of Mathematical Sciences, Chongqing Normal University, ChongQing, 401331, China. E-mails: 20190009@cqnu.edu.cn and Wei Yu, School of Biomedical Engineering and Imaging, Xianning Medical College, Hubei University of Science and Technology, Xianning, 437100, China. E-mail: cq_yuwei@163.com.
Abstract: BACKGROUND:Due to the incomplete projection data collected by limited-angle computed tomography (CT), severe artifacts are present in the reconstructed image. Classical regularization methods such as total variation (TV) minimization, ℓ0 minimization, are unable to suppress artifacts at the edges perfectly. Most existing regularization methods are single-objective optimization approaches, stemming from scalarization methods for multiobjective optimization problems (MOP). OBJECTIVE:To further suppress the artifacts and effectively preserve the edge structures of the reconstructed image. METHOD:This study presents a multiobjective optimization model incorporates both data fidelity term and ℓ0-norm of the image gradient as objective functions. It employs an iterative approach different from traditional scalarization methods, using the maximization of structural similarity (SSIM) values to guide optimization rather than minimizing the objective function.The iterative method involves two steps, firstly, simultaneous algebraic reconstruction technique (SART) optimizes the data fidelity term using SSIM and the Simulated Annealing (SA) algorithm for guidance. The degradation solution is accepted in the form of probability, and guided image filtering (GIF) is introduced to further preserve the image edge when the degradation solution is rejected. Secondly, the result from the first step is integrated into the second objective function as a constraint, we use ℓ0 minimization to optimize ℓ0-norm of the image gradient, and the SSIM, SA algorithm and GIF are introduced to guide optimization process by improving SSIM value like the first step. RESULTS:With visual inspection, the peak signal-to-noise ratio (PSNR), root mean square error (RMSE), and SSIM values indicate that our approach outperforms other traditional methods. CONCLUSIONS:The experiments demonstrate the effectiveness of our method and its superiority over other classical methods in artifact suppression and edge detail restoration.
