Discrete tomography by convex–concave regularization and D.C. programming Original Research Article

We present a novel approach to the tomographic reconstruction of binary objects from few projection directions within a limited range of angles. A quadratic objective functional over binary variables comprising the squared projection error and a prior penalizing non-homogeneous regions, is supplemented with a concave functional enforcing binary solutions. Application of a primal-dual subgradient algorithm to a suitable decomposition of the objective functional into the difference of two convex functions leads to an algorithm which provably converges with parallel updates to binary solutions. Numerical results demonstrate robustness against local minima and excellent reconstruction performance using five projections within a range of image.