Enhancing Finite Element-Based Photoacoustic Tomography by Localized Reconstruction Method

Iterative reconstruction algorithm based on finite element (FE) modeling is a powerful approach in photoacoustic tomography (PAT). However, an iterative inverse algorithm using conventional FE meshing of the entire imaging zone is computationally demanding, which hinders this powerful tool in applications where quick image acquisition and/or a large image matrix is needed. To address this challenge, parallel computing techniques are proposed and implemented in the field. Here, we present an alternative approach for 2D PAT, which locoregionally reconstructs the region of interest (ROI) instead of the full imaging zone. Our simulated and phantom experimental results demonstrate that this ROI reconstruction algorithm can produce almost the same image quality as the conventional full zone-based reconstruction algorithm; however, the computation time can be significantly reduced without any additional hardware cost by more than two orders of magnitude (100-fold). This algorithm is further applied and validated in an in vivo study. The major vessel structures in a rat’s brain can be imaged clearly using our ROI-based approach, coupled with a mesh of 11,801 nodes. This novel algorithm can also be parallelized using MPI or GPU acceleration techniques to further enhance the reconstruction performance of FE-based PAT.