Thermoacoustic tomography with detectors on an open curve

An efficient reconstruction algorithm

Research output: Contribution to journalArticle

24 Citations (Scopus)

Abstract

Practical applications of thermoacoustic tomography require numerical inversion of the spherical mean Radon transform with the centers of integration spheres occupying an open surface. A solution of this problem is needed (both in 2D and in 3D) because frequently the region of interest cannot be completely surrounded by the detectors, as happens, for example, in breast imaging. We present an efficient numerical algorithm for solving this problem in 2D (similar methods are applicable in the 3D case). Our method is based on the numerical approximation of plane waves by certain single-layer potentials related to the acquisition geometry. After the densities of these potentials have been precomputed, each subsequent image reconstruction has the complexity of the regular filtration backprojection algorithm for the classical Radon transform. The performance of the method is demonstrated in several numerical examples: one can see that the algorithm produces very accurate reconstructions if the data are accurate and sufficiently well sampled; on the other hand, it is sufficiently stable with respect to noise in the data.

Original languageEnglish (US)
Article number055021
JournalInverse Problems
Volume24
Issue number5
DOIs
StatePublished - Oct 1 2008

Fingerprint

Thermoacoustics
Reconstruction Algorithm
Tomography
Efficient Algorithms
Radon Transform
Radon
Detector
Detectors
Curve
Spherical Means
Single Layer Potential
Numerical Inversion
Image Reconstruction
Region of Interest
Image reconstruction
Numerical Approximation
Numerical Algorithms
Plane Wave
Filtration
Imaging

ASJC Scopus subject areas

  • Signal Processing
  • Computer Science Applications
  • Applied Mathematics
  • Mathematical Physics
  • Theoretical Computer Science

Cite this

Thermoacoustic tomography with detectors on an open curve : An efficient reconstruction algorithm. / Kunyansky, Leonid.

In: Inverse Problems, Vol. 24, No. 5, 055021, 01.10.2008.

Research output: Contribution to journalArticle

@article{0b58d1748fcd4732a91b062f075be050,
title = "Thermoacoustic tomography with detectors on an open curve: An efficient reconstruction algorithm",
abstract = "Practical applications of thermoacoustic tomography require numerical inversion of the spherical mean Radon transform with the centers of integration spheres occupying an open surface. A solution of this problem is needed (both in 2D and in 3D) because frequently the region of interest cannot be completely surrounded by the detectors, as happens, for example, in breast imaging. We present an efficient numerical algorithm for solving this problem in 2D (similar methods are applicable in the 3D case). Our method is based on the numerical approximation of plane waves by certain single-layer potentials related to the acquisition geometry. After the densities of these potentials have been precomputed, each subsequent image reconstruction has the complexity of the regular filtration backprojection algorithm for the classical Radon transform. The performance of the method is demonstrated in several numerical examples: one can see that the algorithm produces very accurate reconstructions if the data are accurate and sufficiently well sampled; on the other hand, it is sufficiently stable with respect to noise in the data.",
author = "Leonid Kunyansky",
year = "2008",
month = "10",
day = "1",
doi = "10.1088/0266-5611/24/5/055021",
language = "English (US)",
volume = "24",
journal = "Inverse Problems",
issn = "0266-5611",
publisher = "IOP Publishing Ltd.",
number = "5",

}

TY - JOUR

T1 - Thermoacoustic tomography with detectors on an open curve

T2 - An efficient reconstruction algorithm

AU - Kunyansky, Leonid

PY - 2008/10/1

Y1 - 2008/10/1

N2 - Practical applications of thermoacoustic tomography require numerical inversion of the spherical mean Radon transform with the centers of integration spheres occupying an open surface. A solution of this problem is needed (both in 2D and in 3D) because frequently the region of interest cannot be completely surrounded by the detectors, as happens, for example, in breast imaging. We present an efficient numerical algorithm for solving this problem in 2D (similar methods are applicable in the 3D case). Our method is based on the numerical approximation of plane waves by certain single-layer potentials related to the acquisition geometry. After the densities of these potentials have been precomputed, each subsequent image reconstruction has the complexity of the regular filtration backprojection algorithm for the classical Radon transform. The performance of the method is demonstrated in several numerical examples: one can see that the algorithm produces very accurate reconstructions if the data are accurate and sufficiently well sampled; on the other hand, it is sufficiently stable with respect to noise in the data.

AB - Practical applications of thermoacoustic tomography require numerical inversion of the spherical mean Radon transform with the centers of integration spheres occupying an open surface. A solution of this problem is needed (both in 2D and in 3D) because frequently the region of interest cannot be completely surrounded by the detectors, as happens, for example, in breast imaging. We present an efficient numerical algorithm for solving this problem in 2D (similar methods are applicable in the 3D case). Our method is based on the numerical approximation of plane waves by certain single-layer potentials related to the acquisition geometry. After the densities of these potentials have been precomputed, each subsequent image reconstruction has the complexity of the regular filtration backprojection algorithm for the classical Radon transform. The performance of the method is demonstrated in several numerical examples: one can see that the algorithm produces very accurate reconstructions if the data are accurate and sufficiently well sampled; on the other hand, it is sufficiently stable with respect to noise in the data.

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

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

U2 - 10.1088/0266-5611/24/5/055021

DO - 10.1088/0266-5611/24/5/055021

M3 - Article

VL - 24

JO - Inverse Problems

JF - Inverse Problems

SN - 0266-5611

IS - 5

M1 - 055021

ER -