The spherical Radon transform with centers on cylindrical surfaces

Abstract

Recovering a function from its spherical Radon transform with centers of spheres of integration restricted to a hypersurface is at the heart of several modern imaging technologies, including SAR, ultrasound imaging, and photo- and thermoacoustic tomography. In this paper we study an inversion of the spherical Radon transform with centers of integration restricted to cylindrical surfaces of the form Gamma x R-m, where r is a hypersurface in R-n. We show that this transform can be decomposed into two lower dimensional spherical Radon transforms, one with centers on r and one with a planar center-set in Rm+1. Together with explicit inversion formulas for the spherical Radon transform with a planar center-set and existing algorithms for inverting the spherical Radon transform with a center-set Gamma, this yields reconstruction procedures for general cylindrical domains. In the special case of spherical or elliptical cylinders we obtain novel explicit inversion formulas. For three spatial dimensions, these inversion formulas can be implemented efficiently by backprojection type algorithms only requiring Omicron(N-4/3) floating point operations, where N is the total number of unknowns to be recovered. We present numerical results demonstrating the efficiency of the derived algorithms.