Iterative Image Reconstruction for CT with Unmatched Projection Matrices Using the Generalized Minimal Residual Algorithm

Emil Y. Sidky,Per Christian Hansen,Jakob Sauer Jørgensen,Xiaochuan Pan
DOI: https://doi.org/10.48550/arxiv.2201.07408
2022-01-01
Abstract:The generalized minimal residual (GMRES) algorithm is applied to image reconstruction using linear computed tomography (CT) models. The GMRES algorithm iteratively solves square, non-symmetric linear systems and it has practical application to CT when using unmatched back-projector/projector pairs and when applying preconditioning. The GMRES algorithm is demonstrated on a 3D CT image reconstruction problem where it is seen that use of unmatched projection matrices does not prevent convergence, while using an unmatched pair in the related conjugate gradients for least-squares (CGLS) algorithm leads to divergent iteration. Implementation of preconditioning using GMRES is also demonstrated.
What problem does this paper attempt to address?