Characterization and suppression of edge and aliasing artefacts in iterative x-ray CT reconstruction.

For the purpose of obtaining x-ray tomographic images, statistical reconstruction (SR) provides a general framework with possible advantages over analytical algorithms such as filtered backprojection (FBP) in terms of flexibility, resolution, contrast and image noise. However, SR images may be seriously affected by some artefacts that are not present in FBP images. These artefacts appear as aliasing patterns and as severe overshoots in the areas of sharp intensity transitions ('edge artefacts'). We characterize this inherent property of iterative reconstructions and hypothesize how discretization errors during reconstruction contribute to the formation of the artefacts. An adequate solution to the problem is to perform the reconstructions on an image grid that is finer than that typically employed for FBP reconstruction, followed by a downsampling of the resulting image to a granularity normally used for display. Furthermore, it is shown that such a procedure is much more effective than post-filtering of the reconstructions. Resulting SR images have superior noise-resolution trade-off compared to FBP, which may facilitate dose reduction during CT examinations.

[1]  R. Lewitt Alternatives to voxels for image representation in iterative reconstruction algorithms , 1992, Physics in medicine and biology.

[2]  Gengsheng Lawrence Zeng,et al.  Unmatched projector/backprojector pairs in an iterative reconstruction algorithm , 2000, IEEE Transactions on Medical Imaging.

[3]  R. Leahy,et al.  Recent developments in iterative image reconstruction for PET and SPECT [Editorial] , 2000, IEEE Transactions on Medical Imaging.

[4]  Patrick Dupont,et al.  Reduction of metal streak artifacts in X-ray computed tomography using a transmission maximum a posteriori algorithm , 1999 .

[5]  Hakan Erdogan,et al.  Ordered subsets algorithms for transmission tomography. , 1999, Physics in medicine and biology.

[6]  Wojciech Zbijewski,et al.  Suppression of intensity transition artifacts in statistical x-ray computer tomography reconstruction through radon inversion initialization. , 2003, Medical physics.

[7]  Jeffrey A. Fessler,et al.  Statistical image reconstruction for polyenergetic X-ray computed tomography , 2002, IEEE Transactions on Medical Imaging.

[8]  K. Lange Convergence of EM image reconstruction algorithms with Gibbs smoothing. , 1990, IEEE transactions on medical imaging.

[9]  Klaus Mueller,et al.  Rapid 3-D cone-beam reconstruction with the simultaneous algebraic reconstruction technique (SART) using 2-D texture mapping hardware , 2000, IEEE Transactions on Medical Imaging.

[10]  Robert M. Lewitt,et al.  Practical considerations for 3-D image reconstruction using spherically symmetric volume elements , 1996, IEEE Trans. Medical Imaging.

[11]  G.T. Gullberg,et al.  Frequency domain implementation of the three-dimensional geometric point response correction in SPECT imaging , 1991, Conference Record of the 1991 IEEE Nuclear Science Symposium and Medical Imaging Conference.

[12]  L. Shepp,et al.  Maximum Likelihood Reconstruction for Emission Tomography , 1983, IEEE Transactions on Medical Imaging.

[13]  Ronald J. Jaszczak,et al.  Modeling the axial extension of a transmission line source within iterative reconstruction via multiple transmission sources , 2002, IEEE Transactions on Medical Imaging.

[14]  L. J. Thomas,et al.  Noise and Edge Artifacts in Maximum-Likelihood Reconstructions for Emission Tomography , 1987, IEEE Transactions on Medical Imaging.

[15]  R. Siddon Fast calculation of the exact radiological path for a three-dimensional CT array. , 1985, Medical physics.

[16]  M. Defrise,et al.  Iterative reconstruction for helical CT: a simulation study. , 1998, Physics in medicine and biology.

[17]  Freek J. Beekman,et al.  Accelerated iterative transmission CT reconstruction using an ordered subsets convex algorithm , 1998, IEEE Transactions on Medical Imaging.

[18]  B. De Man,et al.  Distance-driven projection and backprojection , 2002, 2002 IEEE Nuclear Science Symposium Conference Record.

[19]  F. Beekma,et al.  Ordered subset reconstruction for x-ray CT. , 2001, Physics in medicine and biology.

[20]  Simon R Cherry,et al.  Effect of phantom voxelization in CT simulations. , 2002, Medical physics.

[21]  K. Lange,et al.  EM reconstruction algorithms for emission and transmission tomography. , 1984, Journal of computer assisted tomography.

[22]  R M Lewitt,et al.  Multidimensional digital image representations using generalized Kaiser-Bessel window functions. , 1990, Journal of the Optical Society of America. A, Optics and image science.

[23]  Max A. Viergever,et al.  Dual matrix ordered subsets reconstruction for accelerated 3D scatter compensation in single-photon emission tomography , 1997, European Journal of Nuclear Medicine.

[24]  Klaus Mueller,et al.  Fast Implementation of Algebraic Methods for 3D Reconstruction from Cone-Beam Data , 1999, IEEE Trans. Medical Imaging.

[25]  Patrick Dupont,et al.  An iterative maximum-likelihood polychromatic algorithm for CT , 2001, IEEE Transactions on Medical Imaging.

[26]  H Guan,et al.  Computed tomography using algebraic reconstruction techniques (ARTs) with different projection access schemes: a comparison study under practical situations , 1996, Physics in medicine and biology.

[27]  Ge Wang,et al.  An iterative algorithm for X-ray CT fluoroscopy , 1998, IEEE Transactions on Medical Imaging.