Computational Analysis and Improvement of SIRT

Iterative X-ray computed tomography (CT) algorithms have the potential for producing high-quality images but are computationally very demanding, especially when applied to high-resolution problems. Focusing on simultaneous iterative reconstruction technique (SIRT), we provide an eigenvalue based scheme for automatically determining a near-optimal value of the relaxation parameter. This accelerates the convergence rate of SIRT to the point where only half the number of iterations normally required is needed. We also modify the way SIRT uses preconditioning to solve a weighted least squares problem. The resulting algorithm, which we call PSIRT, is associated with a smaller memory footprint and calls for less data to be communicated in a distributed-memory implementation. Experimental residual norm and timing results are provided based on cone-beam micro-CT mouse data, including for an ordered subsets study.