Abstract:This paper studies inverse problems in quantitative photoacoustic tomography with additional optical current data supplemented from diffuse optical tomography. We propose a three-stage image reconstruction method for the simultaneous recovery of the absorption, diffusion, and Grüneisen coefficients. We demonstrate, through numerical simulations, that: (i) when the Grüneisen coefficient is known, the addition of the optical measurements allows a more accurate reconstruction of the scattering and absorption coefficients; and (ii) when the Grüneisen coefficient is not known, the addition of optical current measurements allows us to reconstruct uniquely the Grüneisen, the scattering and absorption coefficients. Numerical simulations based on synthetic data are presented to demonstrate the effectiveness of the proposed idea.
Numerical Analysis,Optimization and Control,Computational Physics,Medical Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to simultaneously reconstruct multiple coefficients in the coupled photoacoustic tomography (PAT) and diffuse optical tomography (DOT). Specifically, the author proposes a three - stage image reconstruction method, aiming to simultaneously recover the absorption coefficient \(\sigma\), the scattering coefficient \(\gamma\) and the Grüneisen coefficient \(\Gamma\).
### Problem Background
1. **Diffuse Optical Tomography (DOT)**
- DOT uses near - infrared (NIR) light to detect the scattering and absorption characteristics of optically inhomogeneous media such as biological tissues.
- Its main advantage is that it can image the optical contrast of biological tissues, but its resolution is low and it is difficult to recover high - frequency information from the measured data.
2. **Photoacoustic Tomography (PAT)**
- PAT is a hybrid imaging method that combines ultrasonic imaging with diffuse optical imaging through the photoacoustic effect to achieve high - resolution imaging.
- The photoacoustic effect is caused by short - pulse NIR light irradiating an optically inhomogeneous medium, where some photons are absorbed and converted into thermal energy, and then an ultrasonic signal is generated.
### Research Objectives
- **Main Challenge**: In quantitative photoacoustic tomography (QPAT), reconstructing the optical properties inside the medium from the measured data is an inverse problem, which is usually divided into two stages: first, reconstruct the initial pressure field, and then reconstruct the optical coefficients.
- **Innovation**: This paper studies adding additional boundary optical current data from DOT in QPAT and proposes a three - stage reconstruction method to improve the reconstruction accuracy and be able to uniquely reconstruct the three optical coefficients when the Grüneisen coefficient is unknown.
### Solutions
1. **First Stage**: Reconstruct the scattering coefficient \(\gamma\) on the boundary.
2. **Second Stage**: Reconstruct the scattering coefficient \(\gamma\) and the absorption coefficient \(\sigma\) in the entire region.
3. **Third Stage**: Reconstruct the Grüneisen coefficient \(\Gamma\).
### Numerical Simulation Results
Through numerical simulation, the author shows the effectiveness of the proposed method:
- When the Grüneisen coefficient is known, adding optical measurement data can more accurately reconstruct the scattering and absorption coefficients.
- When the Grüneisen coefficient is unknown, adding optical current measurement data can uniquely reconstruct all three coefficients.
### Key Formulas
- Diffusion Equation (DOT):
\[
\frac{1}{c} \frac{\partial u}{\partial t}-\nabla\cdot(\gamma\nabla u)+\sigma u = 0\quad\text{in }\mathbb{R}^+\times\Omega
\]
- Boundary Conditions:
\[
u+\kappa\gamma\frac{\partial u}{\partial n}=g(t, x)\quad\text{on }\mathbb{R}^+\times\partial\Omega
\]
- Optical Current Data:
\[
J_g(t, x):=\left.\gamma\frac{\partial u}{\partial n}\right|_{(0, T]\times\partial\Omega}
\]
- Pressure Distribution (PAT):
\[
H_g(x)=\Gamma(x)\sigma(x)U(x)\quad\text{in }\bar{\Omega}
\]
### Summary
By introducing additional optical measurement data, this paper proposes a three - stage reconstruction method and successfully solves the inverse problem of simultaneously reconstructing multiple optical coefficients in the coupled PAT and DOT.