Reconstruction of sparse-view X-ray computed tomography using adaptive iterative algorithms

Two-and-a-half order score-based model for solving 3D ill-posed inverse problems

Computed Tomography (CT) and Magnetic Resonance Imaging (MRI) are crucial technologies in the field of medical imaging. Score-based models demonstrated effectiveness in addressing different inverse problems encountered in the field of CT and MRI, such as sparse-view CT and fast MRI reconstruction. However, these models face challenges in achieving accurate three dimensional (3D) volumetric reconstruction. The existing score-based models predominantly concentrate on reconstructing two-dimensional (2D) data distributions, resulting in inconsistencies between adjacent slices in the reconstructed 3D volumetric images. To overcome this limitation, we propose a novel two-and-a-half order score-based model (TOSM). During the training phase, our TOSM learns data distributions in 2D space, simplifying the training process compared to working directly on 3D volumes. However, during the reconstruction phase, the TOSM utilizes complementary scores along three directions (sagittal, coronal, and transaxial) to achieve a more precise reconstruction. The development of TOSM is built on robust theoretical principles, ensuring its reliability and efficacy. Through extensive experimentation on large-scale sparse-view CT and fast MRI datasets, our method achieved state-of-the-art (SOTA) results in solving 3D ill-posed inverse problems, averaging a 1.56 dB peak signal-to-noise ratio (PSNR) improvement over existing sparse-view CT reconstruction methods across 29 views and 0.87 dB PSNR improvement over existing fast MRI reconstruction methods with × 2 acceleration. In summary, TOSM significantly addresses the issue of inconsistency in 3D ill-posed problems by modeling the distribution of 3D data rather than 2D distribution which has achieved remarkable results in both CT and MRI reconstruction tasks.

Reconstruction of Sparse-View X-ray Computed Tomography Based on Adaptive Total Variation Minimization

Sparse-view reconstruction has garnered significant interest in X-ray computed tomography (CT) imaging owing to its ability to lower radiation doses and enhance detection efficiency. Among current methods for sparse-view CT reconstruction, an algorithm utilizing iterative reconstruction based on full variational regularization demonstrates good performance. The optimized direction and number of computations for the gradient operator of the regularization term play a crucial role in determining not only the reconstructed image quality but also the convergence speed of the iteration process. The conventional TV approach solely accounts for the vertical and horizontal directions of the two-dimensional plane in the gradient direction. When projection data decrease, the edges of the reconstructed image become blurred. Exploring too many gradient directions for TV terms often comes at the expense of more computational costs. To enhance the balance of computational cost and reconstruction quality, this study suggests a novel TV computation model that is founded on a four-direction gradient operator. In addition, selecting appropriate iteration parameters significantly impacts the quality of the reconstructed image. We propose a nonparametric control method utilizing the improved TV approach as a solution to the tedious manual parameter optimization issue. The relaxation parameters of projection onto convex sets (POCS) are determined according to the scanning number and numerical proportion of the projection data; according to the image error before and after iterations, the gradient descent step of the TV item is adaptively adjusted. Compared with several representative iterative reconstruction algorithms, the experimental results show that the algorithm can effectively preserve edges and suppress noise in sparse-view CT reconstruction.

Fast alternating projection methods for constrained tomographic reconstruction

The alternating projection algorithms are easy to implement and effective for large-scale complex optimization problems, such as constrained reconstruction of X-ray computed tomography (CT). A typical method is to use projection onto convex sets (POCS) for data fidelity, nonnegative constraints combined with total variation (TV) minimization (so called TV-POCS) for sparse-view CT reconstruction. However, this type of method relies on empirically selected parameters for satisfactory reconstruction and is generally slow and lack of convergence analysis. In this work, we use a convex feasibility set approach to address the problems associated with TV-POCS and propose a framework using full sequential alternating projections or POCS (FS-POCS) to find the solution in the intersection of convex constraints of bounded TV function, bounded data fidelity error and non-negativity. The rationale behind FS-POCS is that the mathematically optimal solution of the constrained objective function may not be the physically optimal solution. The breakdown of constrained reconstruction into an intersection of several feasible sets can lead to faster convergence and better quantification of reconstruction parameters in a physical meaningful way than that in an empirical way of trial-and-error. In addition, for large-scale optimization problems, first order methods are usually used. Not only is the condition for convergence of gradient-based methods derived, but also a primal-dual hybrid gradient (PDHG) method is used for fast convergence of bounded TV. The newly proposed FS-POCS is evaluated and compared with TV-POCS and another convex feasibility projection method (CPTV) using both digital phantom and pseudo-real CT data to show its superior performance on reconstruction speed, image quality and quantification.