Jason Zalev,Michael C. Kolios
Abstract:Simulation involves predicting responses of a physical system. In this article, we simulate opto-acoustic signals generated in a three-dimensional volume due to absorption of an optical pulse. A separable computational model is developed that permits an order-of-magnitude improvement in computational efficiency over a non-separable model. The simulated signals represent acoustic waves, measured by a probe with a linear transducer array, in a rotated and translated coordinate frame. Light is delivered by an optical source that moves with the probe's frame. A spatio-temporal impulse response for rectangular element transducer geometry is derived using a Green's function solution to the acoustic wave equation. The approach permits fast and accurate simulation for a probe with arbitrary trajectory. For a 3D volume of $n^3$ voxels, computation is accelerated by a factor of $n$. This may potentially have application for opto-acoustic imaging, where clinicians visualize structural and functional features of biological tissue for assessment of cancer and other diseases.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: how to significantly improve the computational efficiency of three - dimensional photoacoustic signal simulation while ensuring accuracy. Specifically, the author has developed a separable computational model, which divides the processing into two steps, thus achieving a computational efficiency improvement of one order of magnitude higher than that of the non - separable model.
### Problem Background
Photoacoustic imaging is an emerging medical imaging technique. It forms images by detecting photoacoustic signals generated after biological tissues absorb near - infrared laser pulses. In order to accurately reconstruct these images, fast and accurate simulation of photoacoustic signals is required. However, the traditional three - dimensional photoacoustic signal simulation methods have high computational complexity, especially when using more accurate transducer models, which has become an important factor limiting its application.
### Main Contributions of the Paper
1. **Separable Mathematical Model**: The author proposes a new mathematical model, which can equate the 3D Green's function model to the cascade of two perpendicular 2D operations. This separability greatly improves the computational efficiency.
2. **Simulation of Arbitrary Trajectories**: This model can quickly and accurately simulate the data acquisition process when the linear array probe moves along an arbitrary trajectory, which is especially important for the situation where a handheld probe is freely scanned on the skin surface.
3. **Theoretical Proof and Implementation**: The author not only provides detailed mathematical derivations but also verifies the effectiveness of the proposed method through experiments and compares it with other existing 3D simulation methods.
### Formula Representation
The key formulas involved in the paper include:
- **Green's Function Solution**:
\[
g(\mathbf{x}, t)=\frac{1}{4\pi t}\delta(ct - \|\mathbf{x}\|)
\]
- **Pressure Pulse Response**:
\[
h_f(\mathbf{x}, t)=\left[\frac{\partial g(\mathbf{x}, t)}{\partial t}\right]\ast f(\mathbf{x})
\]
- **System Response Operator**:
\[
H_{A,b}^{f,\Phi}\{\eta\}(\mathbf{u}, t)=h_f(\mathbf{u}, t)\ast_u\left[\eta(A\mathbf{u}+b)\Phi(\mathbf{u})\right]
\]
where $\mathbf{x}$ represents the global coordinates, $\mathbf{u}$ represents the local coordinates of the probe, $A$ is the rotation matrix, $b$ is the translation vector, $f$ is the transducer aperture function, $\Phi$ is the optical energy distribution, and $\eta$ is the photoacoustic conversion efficiency.
### Conclusion
This research provides an efficient three - dimensional photoacoustic signal simulation method, which can significantly reduce the computation time while maintaining high precision. This provides strong support for real - time image reconstruction and analysis in photoacoustic imaging and helps to improve the accuracy and efficiency of clinical diagnosis.