Time reversal for photoacoustic tomography based on the wave equation of Nachman, Smith and Waag

Abstract

The goal of \emph{photoacoustic tomography} (PAT) is to estimate an \emph{initial pressure function} $\phi$ from pressure data measured at a boundary surrounding the object of interest. This paper is concerned with a time reversal method for PAT that is based on the dissipative wave equation of Nachman, Smith and Waag\cite{NaSmWa90}. This equation has the advantage that it is more accurate than the \emph{thermo-viscous} wave equation. For simplicity, we focus on the case of one \emph{relaxation process}. We derive an exact formula for the \emph{time reversal image} $\I$, which depends on the \emph{relaxation time} ${\tau }_1$ and the \emph{compressibility} ${\kappa }_1$ of the dissipative medium, and show $\I(τ_1,κ_1)\toφ$ for ${\kappa }_1\to 0$. This implies that $\I=φ$ holds in the dissipation-free case and that $\I$ is similar to $\phi$ for sufficiently small compressibility ${\kappa }_1$. Moreover, we show for tissue similar to water that the \emph{small wave number approximation} $\I_0$ of the time reversal image satisfies $\I_0 = η_0 *_\x φ$ with $\hat η_0(|\k|)\approx const.$ for $|\k|<< \frac{1}{c_0\,τ_1}$. For such tissue, our theoretical analysis and numerical simulations show that the time reversal image $\I$ is very similar to the initial pressure function $\phi$ and that a resolution of $\sigma \approx 0.036\cdot mm$ is feasible (in the noise-free case).

References in corpus (1)

• Photoacoustic imaging in attenuating acoustic media based on strongly causal models

Cited by in corpus (3)

• Iterative Methods for Photoacoustic Tomography in Attenuating Acoustic Media
• Photoacoustic imaging taking into account thermodynamic attenuation
• Extension of the principle of least action with focus on dissipative equations