Abstract:Photoacoustic imaging (PAI) suffers from inherent limitations that can degrade the quality of reconstructed results, such as noise, artifacts, and incomplete data acquisition caused by sparse sampling or partial array detection. In this study, we proposed a new optimization method for both two-dimensional (2D) and three-dimensional (3D) PAI reconstruction results, called regularized iteration method with shape prior. The shape prior is a probability matrix derived from the reconstruction results of multiple sets of random partial array signals in a computational imaging system using any reconstruction algorithm, such as Delay-and-Sum (DAS) and Back-Projection (BP). In the probability matrix, high-probability locations indicate high consistency among multiple reconstruction results at those positions, suggesting a high likelihood of representing the true imaging results. In contrast, low-probability locations indicate higher randomness, leaning more towards noise or artifacts. As a shape prior, this probability matrix guides the iteration and regularization of the entire array signal reconstruction results using the original reconstruction algorithm (the same algorithm for processing random partial array signals). The method takes advantage of the property that the similarity of the object to be imitated is higher than that of noise or artifact in the results reconstructed by multiple sets of random partial array signals of the entire imaging system. The probability matrix is taken as a prerequisite for improving the original reconstruction results, and the optimizer is used to further iterate the imaging results to remove noise and artifacts and improve the imaging fidelity. Especially in the case involving sparse view which brings more artifacts, the effect is remarkable. Simulation and real experiments have both demonstrated the superiority of this method.

Photoacoustic Imaging Reconstruction Algorithm Based on a Combined First and Second Order Total Variation Regularization

A photoacoustic image reconstruction method using total variation and nonconvex optimization

Photoacoustic Imaging Reconstruction Using Combined Nonlocal Patch and Total-Variation Regularization for Straight-Line Scanning

A Photoacoustic Imaging Reconstruction Method Based on Directional Total Variation with Adaptive Directivity

Total Variation Based Gradient Descent Algorithm for Sparse-View Photoacoustic Image Reconstruction

Total Variation-Based Limited-View Photoacoustic Reconstruction Method with Segmentation-Based Regularization

Simultaneous Algebraic Reconstruction Technique Based on Total Variation for Photoacoustic Image

High Total Variation-Based Method for Sparse-View Photoacoustic Reconstruction

Efficient L1/2 Regularization-Based Reconstruction for Photoacoustic Imaging Using Adaptively Iterative Thresholding Algorithm

An Efficient Compensation Method for Limited-View Photoacoustic Imaging Reconstruction Based on Gerchberg–Papoulis Extrapolation

Multispectral interlaced sparse sampling photoacoustic tomography based on directional total variation

Photoacoustic Iterative Optimization Algorithm with Shape Prior Regularization

An Improved Filtered Back-Projection Algorithm for Photoacoustic Tomography

An Iterative Reconstruction Method for Photoacoustic Tomography in Heterogeneous Media Based on Point Source Response

Model‐based quantitative photoacoustic tomography with directional total variation

Effect of Iterative Reconstruction Integrating Sart and Fbp on Photoacoustic Imaging

Efficient Block-Sparse Model-Based Algorithm for Photoacoustic Image Reconstruction

Reconstruction Of High Quality Photoacoustic Tomography With A Limited-View Scanning

Accelerated model‐based iterative reconstruction strategy for sparse‐view photoacoustic tomography aided by multi‐channel autoencoder priors

Photo-acoustic tomographic image reconstruction from reduced data using physically inspired regularization

Low Dose PET Image Reconstruction with Total Variation Using Alternating Direction Method.