Paper detail

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

The goal of \emph{photoacoustic tomography} (PAT) is to estimate an \emph{initial pressure function} $φ$ 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} $τ_1$ and the \emph{compressibility} $κ_1$ of the dissipative medium, and show $\I(τ_1,κ_1)\toφ$ for $κ_1\to 0$. This implies that $\I=φ$ holds in the dissipation-free case and that $\I$ is similar to $φ$ for sufficiently small compressibility $κ_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 $φ$ and that a resolution of $σ\approx 0.036\cdot mm$ is feasible (in the noise-free case).