Experiments with the nonlinear and chaotic behaviour of the multiplicative algebraic reconstruction technique (MART) algorithm for computed tomography

Among the iterative reconstruction algorithms for tomography, the multiplicative algebraic reconstruction technique (MART) has two advantages that make it stand out from other algorithms: it confines the image (and therefore the projection data) to the convex hull of the patient, and it maximizes entropy. In this paper, we have undertaken a series of experiments to determine the importance of MART nonlinearity to image quality. Variants of MART were implemented aiming to exploit and exaggerate the nonlinear properties of the algorithm. We introduce the Power MART, Boxcar Averaging MART and Bouncing MART algorithms. Power MART is linked to the relaxation concept. Its behaviour is similar to that of the chaos of a logistic equation. There appears to be an antagonism between increasing nonlinearity and noise in the projection data. The experiments confirm our general observation that regularization as a means of solving simultaneous linear equations that are underdetermined is suboptimal: it does not necessarily select the correct image from the hyperplane of solutions, and so does not maximize the image quality:x-ray dose ratio. Our investigations prove that there is scope to optimize CT algorithms and thereby achieve greater dose reduction.