TY - GEN

T1 - A geometric approach to quadratic optimization in computerized tomography

AU - Gordon, Dan

AU - Mansour, Rawia

PY - 2005

Y1 - 2005

N2 - 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 Reconstruction Technique), quadratic optimization (QUAD - a method that minimizes the L2-norm of the residual) and the relatively new component-averaging (CAV) and BICAV algorithms. 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 BICAV. 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 efficiency, and the achieved error measures.

AB - 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 Reconstruction Technique), quadratic optimization (QUAD - a method that minimizes the L2-norm of the residual) and the relatively new component-averaging (CAV) and BICAV algorithms. 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 BICAV. 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 efficiency, and the achieved error measures.

KW - Computerized tomography

KW - Image reconstruction

KW - Medical imaging

KW - Normal equations

KW - Quadratic optimization

UR - http://www.scopus.com/inward/record.url?scp=33644560292&partnerID=8YFLogxK

M3 - Conference contribution

AN - SCOPUS:33644560292

SN - 0889865140

T3 - Proceedings of the Eighth IASTED International Conference on Computer Graphics and Imaging, CGIM 2005

SP - 60

EP - 65

BT - Proceedings of the Eighth IASTED International Conference on Computer Graphics and Imaging, CGIM 2005

A2 - Hamza, M.H.

T2 - Eighth IASTED International Conference on Computer Graphics and Imaging, CGIM 2005

Y2 - 15 August 2005 through 17 August 2005

ER -