A Geometric Approach to Quadratic Optimization in Computerized Tomography

D. Gordon and R. Mansour (Israel)

Keywords

Computerized tomography, image reconstruction, medical imaging, normal equations, quadratic optimization.

Abstract

The problem of image reconstruction from projections in computerized tomography, when cast as a system of linear equations, leads to an inconsistent system. The problem is studied under very adverse conditions, consisting of low contrast images and a strongly underdetermined system, where the number of equations is only about 25% of the number of variables (fewer equations enable the use of less radiation). Various algorithms are examined with respect to their performance in such cases: ART (Algebraic Recon struction Technique), quadratic optimization (QUAD – a method that minimizes the L2-norm of the residual) and the relatively new component-averaging(CAV) and BICAV al gorithms. A variant of QUAD, called NQUAD, is obtained by normalizing the equations before applying QUAD. This has a geometric significance since the resulting system is independent of any particular algebraic representation of the equations – a property shared by ART, CAV and BI CAV. Experiments with phantom reconstructions show that NQUAD is always preferable to QUAD. Under the stated adverse conditions, NQUAD is much better than the other studied algorithms, in terms of image quality, runtime effi ciency, and the achieved error measures.

Important Links:



Go Back