Region-of-Interest reconstruction from truncated cone-beam projections
Region-of-Interest (ROI) tomography aims at reconstructing a region of interest inside a body using only x-ray projections intersecting with the goal to reduce overall radiation exposure when only a small specific region of the body needs to be examined. We consider x-ray acquisition from sources located on a smooth curve in verifying classical Tuy's condition. In this situation, the {\it non-trucated} cone-beam transform of smooth densities admits an explicit inverse ; however cannot directly reconstruct from ROI-truncated projections. To deal with the ROI tomography problem, we introduce a novel reconstruction approach. For densities in where is a bounded ball in , our method iterates an operator combining ROI-truncated projections, inversion by the operator and appropriate regularization operators. Assuming only knowledge of projections corresponding to a spherical ROI , given , we prove that if is sufficiently large our iterative reconstruction algorithm converges uniformly to an -accurate approximation of , where the accuracy depends on the regularity of quantified in the Sobolev norm . This result shows the existence of a critical ROI radius ensuring the convergence of the ROI reconstruction algorithm to -accurate approximations of . We numerically verified these theoretical results using simulated acquisition of ROI-truncated cone-beam projection data for multiple acquisition geometries. Numerical experiments indicate that the critical ROI radius is fairly small with respect to the support region~.