Iterative image reconstruction for CT with unmatched projection matrices using the generalized minimal residual algorithm
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.
- Publication:
-
7th International Conference on Image Formation in X-Ray Computed Tomography
- Pub Date:
- October 2022
- DOI:
- 10.1117/12.2646511
- arXiv:
- arXiv:2201.07408
- Bibcode:
- 2022SPIE12304E..06S
- Keywords:
-
- Physics - Medical Physics
- E-Print:
- Accepted to the 2022 CT Meeting