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

Research output: Contribution to journalArticle

24 Citations (Scopus)

Abstract

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
Volume6
Issue number1
DOIs
StatePublished - Feb 2012

Fingerprint

Thermoacoustics
Spherical Symmetry
Reconstruction Algorithm
Tomography
Fast Algorithm
Detector
Detectors
Grid
Computational Grid
Line
Time Reversal
Inverse problems
Finite Difference
Inversion
Rotating
Inverse Problem
Circle
Numerical Simulation
Geometry
Computer simulation

Keywords

  • 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

@article{a87863075b944a48952c16c0a9c31da5,
title = "Fast reconstruction algorithms for the thermoacoustic tomography in certain domains with cylindrical or spherical symmetries",
abstract = "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.",
keywords = "Fast algorithms, Integrating detectors, Radon transform, Spherical means, Thermoacoustic tomography",
author = "Leonid Kunyansky",
year = "2012",
month = "2",
doi = "10.3934/ipi.2012.6.111",
language = "English (US)",
volume = "6",
pages = "111--131",
journal = "Inverse Problems and Imaging",
issn = "1930-8337",
publisher = "American Institute of Mathematical Sciences",
number = "1",

}

TY - JOUR

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

AU - Kunyansky, Leonid

PY - 2012/2

Y1 - 2012/2

N2 - 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.

AB - 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.

KW - Fast algorithms

KW - Integrating detectors

KW - Radon transform

KW - Spherical means

KW - Thermoacoustic tomography

UR - http://www.scopus.com/inward/record.url?scp=84857171095&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84857171095&partnerID=8YFLogxK

U2 - 10.3934/ipi.2012.6.111

DO - 10.3934/ipi.2012.6.111

M3 - Article

AN - SCOPUS:84857171095

VL - 6

SP - 111

EP - 131

JO - Inverse Problems and Imaging

JF - Inverse Problems and Imaging

SN - 1930-8337

IS - 1

ER -