A Caltech Library Service

A Forward-Adjoint Operator Pair Based on the Elastic Wave Equation for Use in Transcranial Photoacoustic Computed Tomography

Mitsuhashi, Kenji and Poudel, Joemini and Matthews, Thomas P. and Garcia-Uribe, Alejandro and Wang, Lihong V. and Anastasio, Mark A. (2017) A Forward-Adjoint Operator Pair Based on the Elastic Wave Equation for Use in Transcranial Photoacoustic Computed Tomography. SIAM Journal on Imaging Sciences, 10 (4). pp. 2022-2048. ISSN 1936-4954. PMCID PMC5788322. doi:10.1137/16M1107619.

[img] PDF - Published Version
See Usage Policy.

[img] PDF - Accepted Version
See Usage Policy.


Use this Persistent URL to link to this item:


Photoacoustic computed tomography (PACT) is an emerging imaging modality that exploits optical contrast and ultrasonic detection principles to form images of the photoacoustically induced initial pressure distribution within tissue. The PACT reconstruction problem corresponds to an inverse source problem in which the initial pressure distribution is recovered from measurements of the radiated wavefield. A major challenge in transcranial PACT brain imaging is compensation for aberrations in the measured data due to the presence of the skull. Ultrasonic waves undergo absorption, scattering, and longitudinal-to-shear wave mode conversion as they propagate through the skull. To properly account for these effects, a wave-equation-based inversion method should be employed that can model the heterogeneous elastic properties of the skull. In this work, a forward model based on a finite-difference time-domain discretization of the three-dimensional elastic wave equation is established and a procedure for computing the corresponding adjoint of the forward operator is presented. Massively parallel implementations of these operators employing multiple graphics processing units are also developed. The developed numerical framework is validated and investigated in computer-simulation and experimental phantom studies whose designs are motivated by transcranial PACT applications.

Item Type:Article
Related URLs:
URLURL TypeDescription CentralArticle Paper
Wang, Lihong V.0000-0001-9783-4383
Anastasio, Mark A.0000-0002-3192-4172
Additional Information:© 2017, Society for Industrial and Applied Mathematics. Submitted: 12 December 2016. Accepted: 15 May 2017. Published online: 16 November 2017. The work of the authors was supported in part by NIH awards EB01696301 and 5T32EB01485505 and NSF award DMS1614305.
Funding AgencyGrant Number
NIH Predoctoral Fellowship5T32EB01485505
Subject Keywords:photoacoustic computed tomography, transcranial imaging, image reconstruction, elastic wave equation
Issue or Number:4
Classification Code:AMS Subject Headings 92C55, 94A08
PubMed Central ID:PMC5788322
Record Number:CaltechAUTHORS:20180118-145948329
Persistent URL:
Usage Policy:No commercial reproduction, distribution, display or performance rights in this work are provided.
ID Code:84399
Deposited By: George Porter
Deposited On:18 Jan 2018 23:40
Last Modified:15 Nov 2021 20:19

Repository Staff Only: item control page