Fourier-based reconstruction via alternating direction total variation minimization in linear scan CT

Abstract In this study, we consider a novel form of computed tomography (CT), that is, linear scan CT (LCT), which applies a straight line trajectory. Furthermore, an iterative algorithm is proposed for pseudo-polar Fourier reconstruction through total variation minimization (PPF-TVM). Considering that the sampled Fourier data are distributed in pseudo-polar coordinates, the reconstruction model minimizes the TV of the image subject to the constraint that the estimated 2D Fourier data for the image are consistent with the 1D Fourier transform of the projection data. PPF-TVM employs the alternating direction method (ADM) to develop a robust and efficient iteration scheme, which ensures stable convergence provided that appropriate parameter values are given. In the ADM scheme, PPF-TVM applies the pseudo-polar fast Fourier transform and its adjoint to iterate back and forth between the image and frequency domains. Thus, there is no interpolation in the Fourier domain, which makes the algorithm both fast and accurate. PPF-TVM is particularly useful for limited angle reconstruction in LCT and it appears to be robust against artifacts. The PPF-TVM algorithm was tested with the FORBILD head phantom and real data in comparisons with state-of-the-art algorithms. Simulation studies and real data verification suggest that PPF-TVM can reconstruct higher accuracy images with lower time consumption.

[1]  Stanley Osher,et al.  Development and Optimization of Regularized Tomographic Reconstruction Algorithms Utilizing Equally-Sloped Tomography , 2010, IEEE Transactions on Image Processing.

[2]  Aleksandra Pizurica,et al.  Split-Bregman-based sparse-view CT reconstruction , 2011 .

[3]  J. Miao,et al.  Three-dimensional imaging of dislocations in a nanoparticle at atomic resolution , 2013, Nature.

[4]  A. H. Andersen Algebraic reconstruction in CT from limited views. , 1989, IEEE transactions on medical imaging.

[5]  Junfeng Yang,et al.  Alternating Direction Algorithms for 1-Problems in Compressive Sensing , 2009, SIAM J. Sci. Comput..

[6]  E. Sidky,et al.  Accurate image reconstruction from few-views and limited-angle data in divergent-beam CT , 2009, 0904.4495.

[7]  G. Herman,et al.  Image reconstruction from linograms: implementation and evaluation. , 1988, IEEE transactions on medical imaging.

[8]  Peter Cloetens,et al.  Low-dose x-ray phase-contrast and absorption CT using equally sloped tomography , 2010, Physics in medicine and biology.

[9]  E. Candès,et al.  Stable signal recovery from incomplete and inaccurate measurements , 2005, math/0503066.

[10]  Wotao Yin,et al.  Bregman Iterative Algorithms for (cid:2) 1 -Minimization with Applications to Compressed Sensing ∗ , 2008 .

[11]  Zhiqiang Chen,et al.  An Extrapolation Method for Image Reconstruction from a Straight-line Trajectory , 2006, 2006 IEEE Nuclear Science Symposium Conference Record.

[12]  H. Tuy AN INVERSION FORMULA FOR CONE-BEAM RECONSTRUCTION* , 1983 .

[13]  J. Miao,et al.  Electron tomography at 2.4-ångström resolution , 2012, Nature.

[14]  Yin Zhang,et al.  An efficient augmented Lagrangian method with applications to total variation minimization , 2013, Computational Optimization and Applications.

[15]  Li Zhang,et al.  An Improved Form of Linogram Algorithm for Image Reconstruction , 2008, IEEE Transactions on Nuclear Science.

[16]  Bin Yan,et al.  Edge guided image reconstruction in linear scan CT by weighted alternating direction TV minimization. , 2014, Journal of X-ray science and technology.

[17]  Jeffrey A. Fessler,et al.  A Splitting-Based Iterative Algorithm for Accelerated Statistical X-Ray CT Reconstruction , 2012, IEEE Transactions on Medical Imaging.

[18]  Xiaochuan Pan,et al.  Volume image reconstruction from a straight-line source trajectory , 2005, IEEE Nuclear Science Symposium Conference Record, 2005.

[19]  Junfeng Yang,et al.  A New Alternating Minimization Algorithm for Total Variation Image Reconstruction , 2008, SIAM J. Imaging Sci..

[20]  E. Sidky,et al.  Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization , 2008, Physics in medicine and biology.

[21]  J. Miao,et al.  Equally sloped tomography with oversampling reconstruction , 2005 .

[22]  Xiao Han,et al.  Optimization-based reconstruction of sparse images from few-view projections , 2012, Physics in medicine and biology.

[23]  Bin Yan,et al.  Image reconstruction based on total-variation minimization and alternating direction method in linear scan computed tomography , 2013 .

[24]  Li Zhang,et al.  Direct filtered-backprojection-type reconstruction from a straight-line trajectory , 2007 .

[25]  M. Austin,et al.  Quadrature Overlapped Raised-Cosine Modulation , 1981, IEEE Trans. Commun..

[26]  Junfeng Yang,et al.  A Fast Alternating Direction Method for TVL1-L2 Signal Reconstruction From Partial Fourier Data , 2010, IEEE Journal of Selected Topics in Signal Processing.

[27]  Xiao Han,et al.  Optimization-based image reconstruction from sparse-view data in offset-detector CBCT. , 2013, Physics in medicine and biology.

[28]  Liang Li,et al.  A few-view reweighted sparsity hunting (FRESH) method for CT image reconstruction. , 2013, Journal of X-ray science and technology.

[29]  Tarkeshwar Singh,et al.  Fan-beam reconstruction from a straight line of source points , 1993, IEEE Trans. Medical Imaging.

[30]  Emmanuel J. Candès,et al.  Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information , 2004, IEEE Transactions on Information Theory.

[31]  Bruce D. Smith Image Reconstruction from Cone-Beam Projections: Necessary and Sufficient Conditions and Reconstruction Methods , 1985, IEEE Transactions on Medical Imaging.

[32]  E.J. Candes,et al.  An Introduction To Compressive Sampling , 2008, IEEE Signal Processing Magazine.

[33]  Xiaochuan Pan,et al.  Accurate image reconstruction in circular cone-beam computed tomography by total variation minimization: a preliminary investigation , 2006, 2006 IEEE Nuclear Science Symposium Conference Record.

[34]  Jeffrey A. Fessler,et al.  Parallel MR Image Reconstruction Using Augmented Lagrangian Methods , 2011, IEEE Transactions on Medical Imaging.

[35]  G. Herman,et al.  Linograms in Image Reconstruction from Projections , 1987, IEEE Transactions on Medical Imaging.

[36]  Tom Goldstein,et al.  The Split Bregman Method for L1-Regularized Problems , 2009, SIAM J. Imaging Sci..

[37]  Ronald R. Coifman,et al.  A Framework for Discrete Integral Transformations I-The Pseudopolar Fourier Transform , 2008, SIAM J. Sci. Comput..

[38]  I. Daubechies,et al.  An iterative thresholding algorithm for linear inverse problems with a sparsity constraint , 2003, math/0307152.