Thermoacoustic tomography - implementation of exact backprojection formulas

The problem of image reconstruction in thermoacoustic tomography requires inversion of a generalized Radon transform, which integrates the unknown function over circles in 2D or spheres in 3D. The paper investigates implementation of the recently discovered backprojection type inversion formulas for the case of spherical acquisition in 3D. A numerical simulation of the data acquisition with subsequent reconstructions are made for the Defrise phantom as well as for some other phantoms. Both full and partial scan situations are considered. The results are compared with the implementation of a previously used approximate inversion formula.

[1]  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.

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

[3]  S. Jacques,et al.  Iterative reconstruction algorithm for optoacoustic imaging. , 2002, The Journal of the Acoustical Society of America.

[4]  Eric Todd Quinto,et al.  Singularities of the X-ray transform and limited data tomography , 1993 .

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

[6]  G. Beylkin.  The inversion problem and applications of the generalized radon transform , 1984 .

[7]  R. Kruger,et al.  Breast cancer in vivo: contrast enhancement with thermoacoustic CT at 434 MHz-feasibility study. , 2000, Radiology.

[8]  Eric Todd Quinto,et al.  Geometry of stationary sets for the wave equation in ℝn: The case of finitely supported initial data , 2001 .

[9]  V. Palamodov,et al.  Reconstruction from limited data of arc means , 2000 .

[10]  A. Aisen,et al.  Thermoacoustic CT with radio waves: a medical imaging paradigm. , 1999, Radiology.

[11]  Xiaochuan Pan,et al.  Data redundancy and reduced-scan reconstruction in reflectivity tomography , 2003, IEEE Trans. Image Process..

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

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

[14]  S. Patch,et al.  Thermoacoustic tomography--consistency conditions and the partial scan problem. , 2004, Physics in medicine and biology.

[15]  W. Joines,et al.  The measured electrical properties of normal and malignant human tissues from 50 to 900 MHz. , 1994, Medical physics.

[16]  P. Kuchment,et al.  On local tomography , 1995 .

[17]  M. E. Davison,et al.  The Ill-Conditioned Nature of the Limited Angle Tomography Problem , 1983 .

[18]  Lihong V. Wang,et al.  Reconstructions in limited-view thermoacoustic tomography. , 2004, Medical physics.

[19]  S. Rowland,et al.  Computer implementation of image reconstruction formulas , 1979 .

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

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

[22]  M. E. Davison,et al.  Tomographic reconstruction with arbitrary directions , 1981 .