Fast reconstruction algorithms for the thermoacoustic tomography in certain domains with cylindrical or spherical symmetries

Research output: Contribution to journalArticle

25 Scopus citations


We propose three fast algorithms for solving the inverse problem of the thermoacoustic tomography corresponding to certain acquisition geometries. Two of these methods are designed to process the measurements done with point-like detectors placed on a circle (in 2D) or a sphere (in 3D) surrounding the object of interest. The third inversion algorithm works with the data measured by the integrating line detectors arranged in a cylindrical assembly rotating around the object. The number of operations required by these techniques is equal to O(n 3 log n) and O(n 3 log 2n) for the 3D techniques (assuming the reconstruction grid with n 3 nodes) and to O(n 2 log n) for the 2D problem with n × n discretizetion grid. Numerical simulations show that on large computational grids our methods are at least two orders of magnitude faster than the finite-difference time reversal techniques. The results of reconstructions from real measurements done by the integrating line detectors are also presented, to demonstrate the practicality of our algorithms.

Original languageEnglish (US)
Pages (from-to)111-131
Number of pages21
JournalInverse Problems and Imaging
Issue number1
Publication statusPublished - Feb 2012



  • Fast algorithms
  • Integrating detectors
  • Radon transform
  • Spherical means
  • Thermoacoustic tomography

ASJC Scopus subject areas

  • Analysis
  • Control and Optimization
  • Discrete Mathematics and Combinatorics
  • Modeling and Simulation

Cite this