Rapid Three-Dimensional Photoacoustic Imaging Reconstruction for Irregularly Layered Heterogeneous Media

Photoacoustic imaging (PAI) is susceptible to speed of sound (SOS) differences in heterogeneous media which greatly reduce the resolutions and qualities of the imaging results. Several reconstruction methods have been reported to adapt for heterogenous media, but they are limited by specific deficiencies such as efficiency, accuracy, and model limitation problems. Among them, the plane wave model based on wavefield reconstruction is the most efficient and promising one for high-efficiency three-dimensional PAI. However, the classic plane wave model only suits for planar layered media, severely limiting its applications in practice. To this end, we modify the plane wave model to apply for irregularly layered heterogeneous media and propose a corresponding wavefield extrapolation to reconstruct photoacoustic image. This method employs split-step Fourier to compensate the SOS differences, extrapolates wavefields and reconstructs the image depth by depth. Furthermore, a floating discretization strategy is introduced to control and balance the efficiency and accuracy with a hyperparameter. The simulation and experiment results demonstrate that the proposed method can reconstruct the image with an equivalent resolution to time reversal’s and even have higher efficiency and robustness. To reconstruct a three-dimensional image with $50\times {50}\times 600$ pixels, the proposed method takes only 5.5 seconds using a laptop loaded with Intel(R) Core (TM) i7-8550U CPU @1.8GHz.