Abstract:The paper presents a survey of mathematical problems, techniques, and challenges arising in the Thermoacoustic and Photoacoustic Tomography.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is about the mathematical modeling and challenges in Thermoacoustic Tomography (TAT). Specifically, the paper focuses on the following aspects:
1. **Uniqueness problem**: Research under which conditions the target function \(f(x)\) can be uniquely reconstructed from the measured data. For example, when using a specific measurement surface \(S\), whether the initial value \(f(x)\) can be uniquely determined from the measured data \(g(y, t)\).
2. **Inversion formula and algorithm**: Explore how to effectively reconstruct the target function \(f(x)\) from the measured data \(g(y, t)\). This includes specific mathematical formulas and algorithm implementation.
3. **Stability analysis**: Study the stability of the inversion process, that is, whether small errors in the data will significantly affect the accuracy of the reconstruction results.
4. **Influence of incomplete data**: Analyze how the quality of the reconstruction results will be affected when only part of the data is available.
5. **Spatial description of ideal data**: Describe the space of all possible "ideal" data \(g(y, t)\), that is, the conditions that these data should satisfy. This is very important for improving the inversion algorithm, supplementing incomplete data, and detecting and compensating for certain data errors.
### Specific mathematical problems
1. **Wave equation model**:
- Consider the wave equation:
\[
\begin{cases}
p_{tt} = c^2(x)\Delta p, & t\geq0, x\in\mathbb{R}^3 \\
p(x, 0) = f(x), \\
p_t(x, 0) = 0
\end{cases}
\]
- The goal is to recover the initial value \(f(x)\) from the measured data \(g(y, t)\).
2. **Spherical mean transform model**:
- When the sound speed \(c(x)\) is constant, the spherical mean transform can be used:
\[
(Rf)(x, r)=\int_{|y| = 1}f(x + ry)\,dA(y)
\]
- The measured data \(g(y, t)\) can be regarded as the restriction of \(p(x, t)\) on the measurement surface \(S\), that is:
\[
p(y, t) = g(y, t), \quad y\in S, t\geq0
\]
3. **Uniqueness set**:
- Research which sets \(S\) are uniqueness sets, that is, placing detectors on these sets can provide enough data to uniquely determine \(f(x)\).
- For example, a closed surface \(S\) is usually a uniqueness set.
4. **Non - uniqueness set**:
- Describe which sets \(S\) are non - uniqueness sets, that is, placing detectors on these sets cannot uniquely determine \(f(x)\).
- For example, a straight line or a plane is usually a non - uniqueness set.
### Conclusion
The paper elaborately discusses these problems in thermoacoustic tomography through mathematical methods, providing theoretical basis and guidance for practical applications. The solution of these problems is of great significance for improving the resolution and contrast of thermoacoustic tomography technology.