Full-Wave Iterative Image Reconstruction in Photoacoustic Tomography With Acoustically Inhomogeneous Media

Existing approaches to image reconstruction in photoacoustic computed tomography (PACT) with acoustically heterogeneous media are limited to weakly varying media, are computationally burdensome, and/or cannot effectively mitigate the effects of measurement data incompleteness and noise. In this work...

Full description

Saved in:
Bibliographic Details
Published in:IEEE transactions on medical imaging 2013-06, Vol.32 (6), p.1097-1110
Main Authors: Chao Huang, Kun Wang, Liming Nie, Wang, L. V., Anastasio, M. A.
Format: Article
Language:English
Subjects:
Acoustic heterogeneity
Acoustics
Algorithms
Animals
Computational modeling
Computer Simulation
Haplorhini
Image Processing, Computer-Assisted - methods
Image reconstruction
iterative image reconstruction
Mathematical model
Models, Biological
optoacoustic tomography
Phantoms, Imaging
Photoacoustic Techniques - methods
photoacoustic tomography
Reconstruction algorithms
Skull - anatomy & histology
thermoacoustic tomography
Tomography - methods
Vectors
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Existing approaches to image reconstruction in photoacoustic computed tomography (PACT) with acoustically heterogeneous media are limited to weakly varying media, are computationally burdensome, and/or cannot effectively mitigate the effects of measurement data incompleteness and noise. In this work, we develop and investigate a discrete imaging model for PACT that is based on the exact photoacoustic (PA) wave equation and facilitates the circumvention of these limitations. A key contribution of the work is the establishment of a procedure to implement a matched forward and backprojection operator pair associated with the discrete imaging model, which permits application of a wide-range of modern image reconstruction algorithms that can mitigate the effects of data incompleteness and noise. The forward and backprojection operators are based on the k-space pseudospectral method for computing numerical solutions to the PA wave equation in the time domain. The developed reconstruction methodology is investigated by use of both computer-simulated and experimental PACT measurement data.
ISSN:0278-0062
1558-254X
DOI:10.1109/TMI.2013.2254496