Efficient and accurate computation of spherical mean values at scattered center points

Spherical means are a widespread model in modern imaging modalities like photoacoustic tomography. Besides direct inversion methods for specific geometries, iterative methods are often used as reconstruction scheme such that each iteration asks for the efficient and accurate computation of spherical means. We consider a spectral discretization via trigonometric polynomials such that the computation can be done via nonequispaced fast Fourier transforms. Moreover, a recently developed sparse fast Fourier transform is used in the three dimensional case and gives optimal arithmetic complexity. All theoretical results are illustrated by numerical experiments.

[1]  Lexing Ying,et al.  Sparse Fourier Transform via Butterfly Algorithm , 2008, SIAM J. Sci. Comput..

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

[3]  Frank Natterer,et al.  Photo-acoustic inversion in convex domains , 2012 .

[4]  Lihong V. Wang,et al.  Biomedical Optics: Principles and Imaging , 2007 .

[5]  L. Milne‐Thomson.  A Treatise on the Theory of Bessel Functions , 1945, Nature.

[6]  Otmar Scherzer,et al.  Filtered backprojection for thermoacoustic computed tomography in spherical geometry , 2005, Mathematical Methods in the Applied Sciences.

[7]  Stefan Kunis,et al.  Using NFFT 3---A Software Library for Various Nonequispaced Fast Fourier Transforms , 2009, TOMS.

[8]  V. Ntziachristos,et al.  Model-based optoacoustic inversions with incomplete projection data. , 2011, Medical physics.

[9]  Markus Haltmeier,et al.  Spatial resolution in photoacoustic tomography: effects of detector size and detector bandwidth , 2010 .

[10]  Otmar Scherzer,et al.  Thermoacoustic computed tomography with large planar receivers , 2004 .

[11]  M. Haltmeier,et al.  Exact and approximative imaging methods for photoacoustic tomography using an arbitrary detection surface. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[12]  Eric Todd Quinto,et al.  Helgason's Support Theorem and Spherical Radon Transforms , 2008 .

[13]  Leonid Kunyansky,et al.  Reconstruction of a function from its spherical (circular) means with the centers lying on the surface of certain polygons and polyhedra , 2010, 1009.0288.

[14]  Peter Kuchment,et al.  Uniqueness of reconstruction and an inversion procedure for thermoacoustic and photoacoustic tomography with variable sound speed , 2007, 0706.0598.

[15]  MARKUS HALTMEIER,et al.  A Mollification Approach for Inverting the Spherical Mean Radon Transform , 2011, SIAM J. Appl. Math..

[16]  Otmar Scherzer,et al.  Exact reconstruction in photoacoustic tomography with circular integrating detectors II: Spherical geometry , 2010 .

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

[18]  P. Burgholzer,et al.  Thermoacoustic tomography with integrating area and line detectors , 2005, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[19]  Markus Haltmeier,et al.  Universal Inversion Formulas for Recovering a Function from Spherical Means , 2012, SIAM J. Math. Anal..

[20]  P. Kuchment,et al.  Mathematics of thermoacoustic tomography , 2007, European Journal of Applied Mathematics.

[21]  Leonid Kunyansky.  A series solution and a fast algorithm for the inversion of the spherical mean Radon transform , 2007 .

[22]  F. John.  Plane Waves and Spherical Means: Applied To Partial Differential Equations , 1981 .

[23]  Stefan Kunis,et al.  A Stable and Accurate Butterfly Sparse Fourier Transform , 2012, SIAM J. Numer. Anal..

[24]  Otmar Scherzer,et al.  A Reconstruction Algorithm for Photoacoustic Imaging Based on the Nonuniform FFT , 2009, IEEE Transactions on Medical Imaging.

[25]  P. Heywood.  Trigonometric Series , 1968, Nature.

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

[27]  M. Haltmeier,et al.  Exact Series Reconstruction in Photoacoustic Tomography with Circular Integrating Detectors , 2008 .

[28]  Eric Bonnetier,et al.  Small volume asymptotics for anisotropic elastic inclusions , 2012 .

[29]  Markus Haltmeier,et al.  Inversion of circular means and the wave equation on convex planar domains , 2012, Comput. Math. Appl..

[30]  Frank Filbir,et al.  Reconstruction from circular and spherical mean data , 2010 .