Abstract:BackgroundWith the development of low-dose computed tomography (CT), incomplete data reconstruction has been widely concerned. The total variation (TV) minimization algorithm can accurately reconstruct images from sparse or noisy data. PurposeHowever, the traditional TV algorithm ignores the direction of structures in images, leading to the loss of edge information and block artifacts when the object is not piecewise constant. Since the anisotropic information can facilitate preserving the edge and detail information in images, we aim to improve the TV algorithm in terms of reconstruction accuracy via this approach. MethodsIn this paper, we propose an adaptive-weighted high order total variation (awHOTV) algorithm. We construct the second order TV-norm using the second order gradient, adapt the anisotropic edge property between neighboring image pixels, adjust the local image-intensity gradient to keep edge information, and design the corresponding Chambolle-Pock (CP) solving algorithm. Implementing the proposed algorithm, comprehensive studies are conducted in the ideal projection data experiment where the Structural similarity (SSIM), Root Mean Square Error (RMSE), Contrast to noise ratio (CNR), and modulation transform function (MTF) curves are utilized to evaluate the quality of reconstructed images in statism, structure, spatial resolution, and contrast, respectively. In the noisy data experiment, we further use the noise power spectrum (NPS) curve to evaluate the reconstructed images and compare it with other three algorithms. ResultsWe use the 2D slice in the XCAT phantom, 2D slice in TCIA Challenge data and FORBILD phantom as simulation phantoms and use real bird data for real verification. The results show that, compared with the traditional TV and FBP algorithms, the awHOTV has better performance in terms of RMSE, SSIM, and Pearson correlation coefficient (PCC) under the projected data with different sparsity. In addition, the awHOTV algorithm is robust against the noise of different intensities. ConclusionsThe proposed awHOTV method can reconstruct the images with high accuracy under sparse or noisy data. The awHOTV solves the strip artifacts caused by sparse data in the FBP method. Compared with the TV method, the awHOTV can effectively suppress block artifacts and has good detail protection ability.

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

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

A photoacoustic image reconstruction method using total variation and nonconvex optimization

Simultaneous Algebraic Reconstruction Technique Based on Total Variation for Photoacoustic Image

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

Reweighted Anisotropic Total Variation Minimization for Limited-Angle CT Reconstruction

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

Enhancing finite element-based photoacoustic tomography using total variation minimization

Reconstruction of Sparse-View X-ray Computed Tomography Based on Adaptive Total Variation Minimization

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

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

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

Adaptive-weighted High Order TV Algorithm for Sparse-View CT Reconstruction

Low dose PET reconstruction with total variation regularization.

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

Sparse-view x-ray CT reconstruction via total generalized variation regularization

Sparse-view reconstruction for photoacoustic tomography combining diffusion model with model-based iteration

Enhanced Time-domain Photoacoustic Tomography through Total-variation Minimization

A fast iterative reconstruction method based on the selective total variation for sparse angular CT

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

Adaptive-weighted total variation minimization for sparse data toward low-dose x-ray computed tomography image reconstruction