A series solution and a fast algorithm for the inversion of the spherical mean Radon transform

An explicit series solution is proposed for the inversion of the spherical mean Radon transform. Such an inversion is required in problems of thermo- and photo-acoustic tomography. Closed-form inversion formulae are currently known only for the case when the centres of the integration spheres lie on a sphere surrounding the support of the unknown function, or on certain unbounded surfaces. Our approach results in an explicit series solution for any closed measuring surface surrounding a region for which the eigenfunctions of the Dirichlet Laplacian are explicitly known—such as, for example, cube, finite cylinder, half-sphere etc. In addition, we present a fast reconstruction algorithm applicable in the case when the detectors (the centres of the integration spheres) lie on a surface of a cube. This algorithm reconstructs 3D images thousands times faster than backprojection-type methods.

[1]  Rakesh,et al.  The range of the spherical mean value operator for functions supported in a ball , 2006 .

[2]  Stephen J. Norton,et al.  Reconstruction of a two‐dimensional reflecting medium over a circular domain: Exact solution , 1980 .

[3]  Eric Todd Quinto,et al.  The Radon Transform, Inverse Problems, and Tomography , 2006 .

[4]  Peter Kuchment,et al.  A Range Description for the Planar Circular Radon Transform , 2006, SIAM J. Math. Anal..

[5]  Frank Natterer,et al.  Mathematical methods in image reconstruction , 2001, SIAM monographs on mathematical modeling and computation.

[6]  R. Kruger,et al.  Photoacoustic ultrasound (PAUS)--reconstruction tomography. , 1995, Medical physics.

[7]  E. T. Quinto,et al.  Range descriptions for the spherical mean Radon transform. I. Functions supported in a ball , 2006, math/0606314.

[8]  Minghua Xu,et al.  Erratum: Universal back-projection algorithm for photoacoustic computed tomography [Phys. Rev. E 71, 016706 (2005)] , 2007 .

[9]  Stephen J. Norton,et al.  Ultrasonic Reflectivity Imaging in Three Dimensions: Exact Inverse Scattering Solutions for Plane, Cylindrical, and Spherical Apertures , 1981, IEEE Transactions on Biomedical Engineering.

[10]  John A. Fawcett,et al.  Inversion of N-dimensional spherical averages , 1985 .

[11]  Markus Haltmeier,et al.  Inversion of Spherical Means and the Wave Equation in Even Dimensions , 2007, SIAM J. Appl. Math..

[12]  Minghua Xu,et al.  Time-domain reconstruction for thermoacoustic tomography in a spherical geometry , 2002, IEEE Transactions on Medical Imaging.

[13]  E. T. Quinto,et al.  Range descriptions for the spherical mean , 2006 .

[14]  P. Kuchment.  Generalized Transforms of Radon Type and Their Applications , 2005 .

[15]  Lihong V Wang,et al.  Universal back-projection algorithm for photoacoustic computed tomography , 2005, SPIE BiOS.

[16]  Peter Kuchment,et al.  On the injectivity of the circular Radon transform , 2005 .

[17]  P. Bérard,et al.  Spectres et groupes cristallographiques. II : domaines sphériques , 1980 .

[18]  L. Kunyansky,et al.  Explicit inversion formulae for the spherical mean Radon transform , 2006, math/0609341.

[19]  Eric Todd Quinto,et al.  Injectivity Sets for the Radon Transform over Circles and Complete Systems of Radial Functions , 1996 .

[20]  R A Kruger,et al.  Thermoacoustic computed tomography--technical considerations. , 1999, Medical physics.

[21]  Rakesh,et al.  Determining a Function from Its Mean Values Over a Family of Spheres , 2004, SIAM J. Math. Anal..