Algebraic reconstruction technique explained

frame|right|Animated sequence of reconstruction steps, one iteration.

The algebraic reconstruction technique (ART) is an iterative reconstruction technique used in computed tomography. It reconstructs an image from a series of angular projections (a sinogram). Gordon, Bender and Herman first showed its use in image reconstruction;^[1] whereas the method is known as Kaczmarz method in numerical linear algebra.^{[2] [3]}

An advantage of ART over other reconstruction methods (such as filtered backprojection) is that it is relatively easy to incorporate prior knowledge into the reconstruction process.

ART can be considered as an iterative solver of a system of linear equations

A x = b , where:A is a sparse m x n matrix whose values represent the relative contribution of each output pixel to different points in the sinogram (m being the number of individual values in the sinogram, and n being the number of output pixels);x represents the pixels in the generated (output) image, arranged as a vector, and:b is a vector representing the sinogram. Each projection (row) in the sinogram is made up of a number of discrete values, arranged along the transverse axis. b is made up of all of these values, from each of the individual projections.^[4]

Given a real or complex matrix

A and a real or complex vector b , respectively, the method computes an approximation of the solution of the linear systems of equations as in the following formula,x ^{k+1} = x ^{k} + λ _{k} (b _{i} - \langle a _{i,} x ^{k} \rangle)/(\|a ^{2} _{i\|}) a ^{T} _{i}

where

i = k \bmod m + 1 , a _{i} is the i-th row of the matrix A , b _{i} is the i-th component of the vector b .λ _{k} is an optional relaxation parameter, of the range 0 < λ _{k} \leq 1 . The relaxation parameter is used to slow the convergence of the system. This increases computation time, but can improve the signal-to-noise ratio of the output. In some implementations, the value of λ _{k} is reduced with each successive iteration.^[4]

A further development of the ART algorithm is the simultaneous algebraic reconstruction technique (SART) algorithm.

References

1. Gordon. R. Bender, R. Herman, GT. Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and x-ray photography. Journal of Theoretical Biology. December 1970. 29. 3. 471–81. 5492997. 10.1016/0022-5193(70)90109-8. 1970JThBi..29..471G.
2. Book: Herman, Gabor T.. Fundamentals of computerized tomography : image reconstruction from projections. 2009. Springer. Dordrecht. 978-1-85233-617-2. 2nd.
3. Book: Natterer, F.. The mathematics of computerized tomography. 1986. B.G. Teubner. Stuttgart. 0-471-90959-9.
4. Book: Kak. Avinash. Slaney. Malcolm. Principles of Computerized Tomographic Imaging. limited. 1999. IEEE Press. New York. 978-0898714944. 276–277, 284.