A Kernel Method for CT Reconstruction: a Fast Implementation using   Circulant Matrices
  Modern tomography involves gathering projection data from multiple directions and feeding them into a software algorithm for tomographic reconstruction. We focus our study on image reconstruction from Radon data in the setting of Computerized Tomography (CT), a non-invasive medical procedure that uses an X-rays equipment to produce cross-sectional images of the body. The detectors of the machine measure the X-ray projections through the sample, producing a so-called sinogram, which is actually the Radon transform of the attenuation coefficient. We provide a kernel-based reconstruction algorithm adapted to these kind of data, i.e. the Algebraic Reconstruction Technique (ART). The novel idea of this work is to reduce the complexity of the ART, through a faster implementation of the matrix-vector product, based on the storage of the involved matrix as a circulant matrix. We provide numerical results for both artificially generated (phantom data) and real data. The aim is to study the behaviour of the method with respect to its image reconstruction accuracy and its computational efficiency. In addition, we try to solve the trade-off between accuracy and efficiency and the one between accuracy and numerical stability. 
Maria Angela Narduzzo
math.NA