Backprojection Wiener deconvolution for computed tomographic reconstruction

Analytical CT reconstruction is popular in practice because of its computational efficiency, but it suffers from low reconstruction quality when an insufficient number of projections are used. To address this issue, this paper presents a new analytical method of backprojection Wiener deconvolution (BPWD). BPWD executes backprojection first, and then applies a Wiener deconvolution to the whole backprojected image. The Wiener filter is derived from a ramp filter, enabling the proposed approach to perform reconstruction and denoising simultaneously. The use of a filter after backprojection does not differentiate between real sampled projections and interpolated ones, introducing reconstruction errors. Therefore a weighted ramp filter was applied to increase the contribution of real sampled projections in the reconstruction, thus improving reconstruction quality. Experiments on synthetic data and real phase-contrast x-ray images showed that the proposed approach yields better reconstruction quality compared to the classical filtered backprojection (FBP) method, with comparable reconstruction speed.

[1]  Ivan Lee,et al.  Frequency-based image deblurring with periodic point spread function , 2014 .

[2]  Hengyong Yu,et al.  A General-Thresholding Solution for $l_{p} (0, 2015, IEEE Transactions on Image Processing.

[3]  M. Vannier,et al.  Why do commercial CT scanners still employ traditional, filtered back-projection for image reconstruction? , 2009, Inverse problems.

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

[5]  Pamela C. Cosman,et al.  Evaluating quality of compressed medical images: SNR, subjective rating, and diagnostic accuracy , 1994, Proc. IEEE.

[6]  Franz Franchetti,et al.  Performance analysis of the filtered backprojection image reconstruction algorithms , 2005, Proceedings. (ICASSP '05). IEEE International Conference on Acoustics, Speech, and Signal Processing, 2005..

[7]  G. Herman,et al.  Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and x-ray photography. , 1970, Journal of theoretical biology.

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

[9]  Terry M. Peters,et al.  Image reconstruction from finite numbers of projections , 1973 .

[10]  Steve B. Jiang,et al.  GPU-based fast cone beam CT reconstruction from undersampled and noisy projection data via total variation , 2010 .

[11]  Vladimir Y. Panin,et al.  Fully 3-D PET reconstruction with system matrix derived from point source measurements , 2006, IEEE Transactions on Medical Imaging.

[12]  A. Kak,et al.  Simultaneous Algebraic Reconstruction Technique (SART): A Superior Implementation of the Art Algorithm , 1984, Ultrasonic imaging.

[13]  Zhaoying Bian,et al.  Optimizing a Parameterized Plug-and-Play ADMM for Iterative Low-Dose CT Reconstruction , 2019, IEEE Transactions on Medical Imaging.

[14]  Predrag Sukovic,et al.  Penalized weighted least-squares image reconstruction for dual energy X-ray transmission tomography , 2000, IEEE Transactions on Medical Imaging.

[15]  Andreas Fouras,et al.  Live small-animal X-ray lung velocimetry and lung micro-tomography at the Australian Synchrotron Imaging and Medical Beamline. , 2015, Journal of synchrotron radiation.

[16]  Jing Wang,et al.  Nonlinear sinogram smoothing for low-dose X-ray CT , 2004, IEEE Transactions on Nuclear Science.

[17]  David W. Townsend,et al.  Image Reconstruction for Medical Applications , 1990, Advances in Computer Graphics.

[18]  S Suzuki,et al.  Comparison between an image reconstruction method of filtering backprojection and the filtered backprojection method. , 1988, Applied optics.

[19]  Jiang Hsieh,et al.  Computed Tomography: Principles, Design, Artifacts, and Recent Advances, Fourth Edition , 2022 .

[20]  Ivan Lee,et al.  Backprojection regularization with weighted ramp filter for tomographic reconstruction , 2015, 2015 37th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC).

[21]  Grant T. Gullberg,et al.  Can the backprojection filtering algorithm be as accurate as the filtered backprojection algorithm? , 1994, Proceedings of 1994 IEEE Nuclear Science Symposium - NSS'94.

[22]  Avinash C. Kak,et al.  Principles of computerized tomographic imaging , 2001, Classics in applied mathematics.

[23]  Jeffrey A. Fessler,et al.  Fast Variance Prediction for Iteratively Reconstructed CT Images With Locally Quadratic Regularization , 2017, IEEE Transactions on Medical Imaging.

[24]  C R Crawford,et al.  CT filtration aliasing artifacts. , 1991, IEEE transactions on medical imaging.

[25]  Alvin C. Silva,et al.  Innovations in CT dose reduction strategy: application of the adaptive statistical iterative reconstruction algorithm. , 2010, AJR. American journal of roentgenology.

[26]  W P Segars,et al.  Realistic CT simulation using the 4D XCAT phantom. , 2008, Medical physics.

[27]  M. Glas,et al.  Principles of Computerized Tomographic Imaging , 2000 .

[28]  Kentaro Uesugi,et al.  High‐resolution visualization of airspace structures in intact mice via synchrotron phase‐contrast X‐ray imaging (PCXI) , 2008, Journal of anatomy.

[29]  Gaohang Yu,et al.  Sparse-view x-ray CT reconstruction via total generalized variation regularization , 2014, Physics in medicine and biology.

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

[31]  Bruno De Man,et al.  An outlook on x-ray CT research and development. , 2008, Medical physics.

[32]  Kees Joost Batenburg,et al.  Improving Filtered Backprojection Reconstruction by Data-Dependent Filtering , 2014, IEEE Transactions on Image Processing.

[33]  Jingyan Xu,et al.  Is iterative reconstruction ready for MDCT? , 2009, Journal of the American College of Radiology : JACR.

[34]  J. Leipsic,et al.  State of the Art: Iterative CT Reconstruction Techniques. , 2015, Radiology.

[35]  P. Gilbert.  Iterative methods for the three-dimensional reconstruction of an object from projections. , 1972, Journal of theoretical biology.

[36]  Jeffrey A. Fessler,et al.  Alternating Dual Updates Algorithm for X-ray CT Reconstruction on the GPU , 2015, IEEE Transactions on Computational Imaging.

[37]  Jing Wang,et al.  Penalized weighted least-squares approach to sinogram noise reduction and image reconstruction for low-dose X-ray computed tomography , 2006, IEEE Transactions on Medical Imaging.

[38]  Jie Tang,et al.  Prior image constrained compressed sensing (PICCS): a method to accurately reconstruct dynamic CT images from highly undersampled projection data sets. , 2008, Medical physics.