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.

Recovering Fine Details from Under-Resolved Electron Tomography Data using HOTV Regularization

Low dose PET reconstruction with total variation regularization.

An Adaptive Total Variation Regularization Method for Electrical Resistance Tomography

Reweighted Anisotropic Total Variation Minimization for Limited-Angle CT Reconstruction

Mixed Total Variation and L1 Regularization Method for Optical Tomography Based on Radiative Transfer Equation

Low-dose CT reconstruction via edge-preserving total variation regularization

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

Sampling Limits for Electron Tomography with Sparsity-exploiting Reconstructions

Data-proximal complementary $\ell^1$-TV reconstruction for limited data CT

A Hybrid Regularization Method Combining Tikhonov with Total Variation for Electrical Resistance Tomography

Total Variation Regularization for Tomographic Reconstruction of Cylindrically Symmetric Objects

GPU-based Low Dose CT Reconstruction via Edge-preserving Total Variation Regularization

Better than the Total Variation Regularization

Data-Proximal Complementary l1-TV Reconstruction for Limited Data Computed Tomography

An improved total variation regularization method for electrical resistance tomography

Total variation based Fourier reconstruction and regularization for computer tomography

A nonlinear total variation based computed tomography (CT) image reconstruction method using gradient reinforcement

Total variation regularization for nonlinear fluorescence tomography with an augmented Lagrangian splitting approach

A Second-Order Generalized Total Variation with Improved Alternating Direction Method of Multipliers Algorithm for Electrical Impedance Tomography Reconstruction

Image reconstruction from few views by ℓ0-norm optimization

RTE-based Parameter Reconstruction with T V + L 1 Regularization