On the Adjoint Operator in Photoacoustic Tomography

Photoacoustic Tomography (PAT) is an emerging biomedical "imaging from coupled physics" technique, in which the image contrast is due to optical absorption, but the information is carried to the surface of the tissue as ultrasound pulses. Many algorithms and formulae for PAT image reconstruction hav...

Descripción completa

Guardado en:
Detalles Bibliográficos
Autores principales: Arridge, Simon R., Betcke, Marta M., Cox, Ben T., Lucka, Felix, Treeby, Brad E.
Formato: Artículo
Publicado: 2016
Materias:
Acceso en línea:https://arxiv.org/abs/1602.02027
http://arxiv.org/abs/1602.02027
Etiquetas: Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
_version_ 1879179258241220608
author Arridge, Simon R.
Betcke, Marta M.
Cox, Ben T.
Lucka, Felix
Treeby, Brad E.
author_browse Arridge, Simon R.
Betcke, Marta M.
Cox, Ben T.
Lucka, Felix
Treeby, Brad E.
author_facet Arridge, Simon R.
Betcke, Marta M.
Cox, Ben T.
Lucka, Felix
Treeby, Brad E.
author_sort Arridge, Simon R.
building arXiv
collection Acceso abierto
description Photoacoustic Tomography (PAT) is an emerging biomedical "imaging from coupled physics" technique, in which the image contrast is due to optical absorption, but the information is carried to the surface of the tissue as ultrasound pulses. Many algorithms and formulae for PAT image reconstruction have been proposed for the case when a complete data set is available. In many practical imaging scenarios, however, it is not possible to obtain the full data, or the data may be sub-sampled for faster data acquisition. In such cases, image reconstruction algorithms that can incorporate prior knowledge to ameliorate the loss of data are required. Hence, recently there has been an increased interest in using variational image reconstruction. A crucial ingredient for the application of these techniques is the adjoint of the PAT forward operator, which is described in this article from physical, theoretical and numerical perspectives. First, a simple mathematical derivation of the adjoint of the PAT forward operator in the continuous framework is presented. Then, an efficient numerical implementation of the adjoint using a k-space time domain wave propagation model is described and illustrated in the context of variational PAT image reconstruction, on both 2D and 3D examples including inhomogeneous sound speed. The principal advantage of this analytical adjoint over an algebraic adjoint (obtained by taking the direct adjoint of the particular numerical forward scheme used) is that it can be implemented using currently available fast wave propagation solvers.
doi_str_mv 10.1088/0266-5611/32/11/115012
first_indexed 2026-10-16T01:29:47Z
format Article
fullrecord <oai_dc:dc xmlns:oai_dc="http://www.openarchives.org/OAI/2.0/oai_dc/" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.openarchives.org/OAI/2.0/oai_dc/ https://www.openarchives.org/OAI/2.0/oai_dc.xsd"><identifier>arxiv-1602.02027</identifier><datestamp>2026-06-04</datestamp> <dc:title>On the Adjoint Operator in Photoacoustic Tomography</dc:title> <dc:creator>Arridge, Simon R.</dc:creator> <dc:creator>Betcke, Marta M.</dc:creator> <dc:creator>Cox, Ben T.</dc:creator> <dc:creator>Lucka, Felix</dc:creator> <dc:creator>Treeby, Brad E.</dc:creator> <dc:subject>Numerical Analysis</dc:subject> <dc:description>Photoacoustic Tomography (PAT) is an emerging biomedical "imaging from coupled physics" technique, in which the image contrast is due to optical absorption, but the information is carried to the surface of the tissue as ultrasound pulses. Many algorithms and formulae for PAT image reconstruction have been proposed for the case when a complete data set is available. In many practical imaging scenarios, however, it is not possible to obtain the full data, or the data may be sub-sampled for faster data acquisition. In such cases, image reconstruction algorithms that can incorporate prior knowledge to ameliorate the loss of data are required. Hence, recently there has been an increased interest in using variational image reconstruction. A crucial ingredient for the application of these techniques is the adjoint of the PAT forward operator, which is described in this article from physical, theoretical and numerical perspectives. First, a simple mathematical derivation of the adjoint of the PAT forward operator in the continuous framework is presented. Then, an efficient numerical implementation of the adjoint using a k-space time domain wave propagation model is described and illustrated in the context of variational PAT image reconstruction, on both 2D and 3D examples including inhomogeneous sound speed. The principal advantage of this analytical adjoint over an algebraic adjoint (obtained by taking the direct adjoint of the particular numerical forward scheme used) is that it can be implemented using currently available fast wave propagation solvers. </dc:description> <dc:description>submitted to "Inverse Problems"</dc:description> <dc:date>2016-02-05</dc:date> <dc:date>2016-08-01</dc:date> <dc:type>text</dc:type> <dc:identifier>https://arxiv.org/abs/1602.02027</dc:identifier> <dc:identifier>http://arxiv.org/abs/1602.02027</dc:identifier> <dc:identifier>doi:10.1088/0266-5611/32/11/115012</dc:identifier> </oai_dc:dc>
id arxiv-1602.02027
institution arXiv
last_indexed 2026-10-16T01:29:47Z
publishDate 2016
publishDateRange 2016
publishDateSort 2016
record_format oai_dc
spelling arxiv-1602.020272026-06-04 On the Adjoint Operator in Photoacoustic Tomography Arridge, Simon R. Betcke, Marta M. Cox, Ben T. Lucka, Felix Treeby, Brad E. Numerical Analysis Photoacoustic Tomography (PAT) is an emerging biomedical "imaging from coupled physics" technique, in which the image contrast is due to optical absorption, but the information is carried to the surface of the tissue as ultrasound pulses. Many algorithms and formulae for PAT image reconstruction have been proposed for the case when a complete data set is available. In many practical imaging scenarios, however, it is not possible to obtain the full data, or the data may be sub-sampled for faster data acquisition. In such cases, image reconstruction algorithms that can incorporate prior knowledge to ameliorate the loss of data are required. Hence, recently there has been an increased interest in using variational image reconstruction. A crucial ingredient for the application of these techniques is the adjoint of the PAT forward operator, which is described in this article from physical, theoretical and numerical perspectives. First, a simple mathematical derivation of the adjoint of the PAT forward operator in the continuous framework is presented. Then, an efficient numerical implementation of the adjoint using a k-space time domain wave propagation model is described and illustrated in the context of variational PAT image reconstruction, on both 2D and 3D examples including inhomogeneous sound speed. The principal advantage of this analytical adjoint over an algebraic adjoint (obtained by taking the direct adjoint of the particular numerical forward scheme used) is that it can be implemented using currently available fast wave propagation solvers. submitted to "Inverse Problems" 2016-02-05 2016-08-01 text https://arxiv.org/abs/1602.02027 http://arxiv.org/abs/1602.02027 doi:10.1088/0266-5611/32/11/115012
spellingShingle Numerical Analysis
Arridge, Simon R.
Betcke, Marta M.
Cox, Ben T.
Lucka, Felix
Treeby, Brad E.
On the Adjoint Operator in Photoacoustic Tomography
title On the Adjoint Operator in Photoacoustic Tomography
title_full On the Adjoint Operator in Photoacoustic Tomography
title_fullStr On the Adjoint Operator in Photoacoustic Tomography
title_full_unstemmed On the Adjoint Operator in Photoacoustic Tomography
title_short On the Adjoint Operator in Photoacoustic Tomography
title_sort on the adjoint operator in photoacoustic tomography
topic Numerical Analysis
url https://arxiv.org/abs/1602.02027
http://arxiv.org/abs/1602.02027
work_keys_str_mv AT arridgesimonr ontheadjointoperatorinphotoacoustictomography
AT betckemartam ontheadjointoperatorinphotoacoustictomography
AT coxbent ontheadjointoperatorinphotoacoustictomography
AT luckafelix ontheadjointoperatorinphotoacoustictomography
AT treebybrade ontheadjointoperatorinphotoacoustictomography