3-D weighted CBCT iterative reconstruction with TV minimization from angularly under-sampled projection data

Abstract With increasing cone angle in cone beam computed tomography (CBCT), the CB artifacts due to cone angle in reconstructed images deteriorate, mainly because of increased severity of violating the data sufficiency condition (DSC). The CB artifacts in analytic reconstruction can be partially attributed to the data inconsistency between the conjugate rays that are 180° apart in their view angles; while those in iterative reconstruction result from the hypothesis that the error propagates preferentially along the rays with larger cone angles. The generation of CB artifacts in both sorts of reconstructions is closely related to the cone angle. In this paper, we analyze how the errors in the CB reconstructions are induced and propagated; meanwhile, inspired by the success of and in a way similar to the three-dimensional (3-D) weighted CB-FBP algorithm, we incorporate a similar 3-D weighting into the total-variation (TV) minimization based CB iterative algorithm for reconstructing images from angularly under-sampled projection data and suppressing the CB artifacts in the reconstructed images. Numerical simulation studies with the FORBILD thoracic phantom have been conducted. In the 3-D weighted CB iterative reconstruction, the CB artifacts can be effectively reduced, although there is an obvious distinction in the appearances of the CB artifacts in analytic and iterative reconstructions. In this work, we have developed a 3-D weighted CB iterative reconstruction with TV minimization, which does suppress the CB artifacts and improve the quality of the reconstructed images.

[1]  G. Wang,et al.  A general cone-beam reconstruction algorithm , 1993, IEEE Trans. Medical Imaging.

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

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

[4]  Bruce R. Whiting,et al.  Noise simulation in x-ray CT , 2005, SPIE Medical Imaging.

[5]  L. Feldkamp,et al.  Practical cone-beam algorithm , 1984 .

[6]  Shaojie Tang,et al.  Axial Cone-Beam Reconstruction by Weighted BPF/DBPF and Orthogonal Butterfly Filtering , 2016, IEEE Transactions on Biomedical Engineering.

[7]  Qiong Xu,et al.  Application of projection simulation based on physical imaging model to the evaluation of beam hardening corrections in X-ray transmission tomography , 2008 .

[8]  Tianye Niu,et al.  Three-Dimensional Weighting in Cone Beam FBP Reconstruction and Its Transformation Over Geometries , 2018, IEEE Transactions on Biomedical Engineering.

[9]  J. Hsieh,et al.  A three-dimensional weighted cone beam filtered backprojection (CB-FBP) algorithm for image reconstruction in volumetric CT under a circular source trajectory , 2005, Physics in medicine and biology.

[10]  Hengyong Yu,et al.  Compressed sensing based interior tomography , 2009, Physics in medicine and biology.

[11]  W. Eckelman,et al.  NCRP report no. 93: Ionizing radiation exposure of the population of the United States: National Council on Radiation Protection and Measurements, Bethesda, Maryland (1987). US$15.00 , 1988 .

[12]  Qiu Wang,et al.  A low dose simulation tool for CT systems with energy integrating detectors. , 2013, Medical physics.

[13]  J. Hsieh,et al.  A three-dimensional-weighted cone beam filtered backprojection (CB-FBP) algorithm for image reconstruction in volumetric CT—helical scanning , 2006, Physics in medicine and biology.

[14]  Katsuyuki Taguchi,et al.  A new weighting scheme for cone-beam helical CT to reduce the image noise. , 2004, Physics in medicine and biology.