Technical Note: Convergence analysis of a polyenergetic SART algorithm. Issue 7 (11th June 2015)

• Record Type: Journal Article
• Authors: Humphries, T.
• Is Part Of: Medical physics. Volume 42:Issue 7(2015)
• Page Start: 4007
• Page End: 4014
• Publication Date: 2015-06-11
• Subjects: computerised tomography -- convergence of numerical methods -- image reconstruction -- iterative methods -- Jacobian matrices -- medical image processing; Computed tomography -- Reconstruction -- Matrix theory -- Numerical approximation and analysis; Computerised tomographs -- Biological material, e.g. blood, urine; Haemocytometers -- Digital computing or data processing equipment or methods, specially adapted for specific applications -- Image data processing or generation, in general; computed tomography -- polyenergetic CT -- beam hardening -- algebraic reconstruction technique -- ART -- SART -- nonlinear fixed point iteration; Jacobians -- Computed tomography -- Medical image reconstruction -- Nonlinear dynamics -- Image reconstruction -- Medical X‐ray imaging -- Eigenvalues -- Linear equations -- X‐ray spectra; Medical physics -- Periodicals; Medical physics; Geneeskunde; Natuurkunde; Toepassingen; Biophysics; Periodicals; Periodicals; Electronic journals; 610.153
• DOI: 10.1118/1.4922204
• Languages: English
• ISSNs: 0094-2405

Abstract

Abstract : Purpose: The authors analyze a recently proposed polyenergetic version of the simultaneous algebraic reconstruction technique (SART). This algorithm, denoted polyenergetic SART (pSART), replaces the monoenergetic forward projection operation used by SART with a postlog, polyenergetic forward projection, while leaving the rest of the algorithm unchanged. While the proposed algorithm provides good results empirically, convergence of the algorithm was not established mathematically in the original paper. Methods: The authors analyze pSART as a nonlinear fixed point iteration by explicitly computing the Jacobian of the iteration. A necessary condition for convergence is that the spectral radius of the Jacobian, evaluated at the fixed point, is less than one. A short proof of convergence for SART is also provided as a basis for comparison. Results: The authors show that the pSART algorithm is not guaranteed to converge, in general. The Jacobian of the iteration depends on several factors, including the system matrix and how one models the energy dependence of the linear attenuation coefficient. The authors provide a simple numerical example that shows that the spectral radius of the Jacobian matrix is not guaranteed to be less than one. A second set of numerical experiments using realistic CT system matrices, however, indicates that conditions for convergence are likely to be satisfied in practice. Conclusions: Although pSART is not mathematically guaranteed to converge, their numerical experiments indicate that it will tend to converge at roughly the same rate as SART for system matrices of the type encountered in CT imaging. Thus, the authors conclude that the algorithm is still a useful method for reconstruction of polyenergetic CT data.