Conjugate gradient Mojette reconstruction

Iterative methods are now recognized as powerful tools to solve inverse problems such as tomographic reconstruction. In this paper, the main goal is to present a new reconstruction algorithm made from two components. An iterative algorithm, namely the Conjugate Gradient (CG) method, is used to solve the tomographic problem in the least square (LS) sense for our specific discrete Mojette geometry. The results are compared (with the same geometry) to the corresponding Mojette Filtered Back Projection (FBP) method. In the fist part of the paper, we recall the discrete geometry used to define the projection M and backprojection M* operators. In the second part, the CG algorithm is presented within the context of the Mojette geometry. Noise is then added onto these Mojette projections with respect to the sampling and reconstructions are performed. Finally the Toeplitz block Toeplitz (TBT) character of M*M is demonstrated.