A novel total variation based ring artifact suppression method for CBCT imaging with two-dimensional antiscatter grids.

PURPOSE Two-dimensional antiscatter grids (2DASG) for cone beam computed tomography (CBCT) is a new area of research to reduce scatter intensity, and consequently improve CBCT image quality. One of the challenges in implementation of 2DASGs is their septa shadows that are impinged on the projections. If these artifacts are not corrected, they may lead to ring artifacts in CBCT images. In this work, we present a novel method to suppress ring artifacts in FPD-based CBCT images. METHODS Briefly, our method first detects the locations of 2DASG's septa shadows in projections and then, reduces projection pixel values in septa shadows iteratively until a residual-based convergence criterion is met. To suppress the 2DASG's septa shadows, we developed a total variation minimization (TVM) formulation, referred to as adaptive-diffusive total variation minimization (adTVM), where the diffusivity regularization parameter was adaptively adjusted during each iteration based on the magnitude of the local pixel gradients. To test our method, we have acquired CBCT scans of phantoms using three 2DASG prototypes with different grid geometries. Projections were acquired with a linac mounted CBCT system, operated in offset detector geometry. These projections were then corrected in the following steps: first, projections were corrected using a gantry angle-specific gain correction map; next, projections were corrected by applying our adTVM method. CBCT images were reconstructed using FDK filtered backprojection algorithm. To evaluate adTVM's performance, pixel value statistics and contrast-to-noise ratio (CNR) were compared in CBCT images corrected with and without our adTVM method. RESULTS Without our adTVM method, all three 2DASG prototypes introduced ring artifacts with varying intensities in CBCT images. With our method, significant reduction in ring artifacts was observed in all test cases. Standard deviation of CT numbers was reduced by 7-74% in uniform density phantom CBCT images, CNR was increased by 8-67%, and CT number accuracy of contrast objects embedded in the phantom was preserved. CONCLUSION We propose a new method to suppress ring artifacts caused by the 2DASG's septa shadows in CBCT images. Our initial investigations indicated that adTVM method could substantially reduce such ring artifacts while preserving CT number accuracy and maintaining good spatial resolution. Therefore, our method may potentially play an important role in enabling the implementation of 2DASGs in flat panel detector based CBCT systems.

[1]  Chi-Wang Shu Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[2]  P. Moin,et al.  Fully Conservative Higher Order Finite Difference Schemes for Incompressible Flow , 1998 .

[3]  B. A. Carré,et al.  The Determination of the Optimum Accelerating Factor for Successive Over-relaxation , 1961, Comput. J..

[4]  Jeffrey H Siewerdsen,et al.  A simple approach to measure computed tomography (CT) modulation transfer function (MTF) and noise-power spectrum (NPS) using the American College of Radiology (ACR) accreditation phantom. , 2013, Medical physics.

[5]  Jan Sijbers,et al.  Reduction of ring artefacts in high resolution micro-CT reconstructions. , 2004, Physics in medicine and biology.

[6]  Stephen P. Boyd,et al.  An ADMM Algorithm for a Class of Total Variation Regularized Estimation Problems , 2012, 1203.1828.

[7]  D. Jaffray,et al.  Cone-beam computed tomography with a flat-panel imager: magnitude and effects of x-ray scatter. , 2001, Medical physics.

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

[9]  Carsten Raven,et al.  Numerical removal of ring artifacts in microtomography , 1998 .

[10]  Moyed Miften,et al.  Two‐dimensional antiscatter grid: A novel scatter rejection device for Cone‐beam computed tomography , 2018, Medical physics.

[11]  D. Conover,et al.  Cone beam volume CT image artifacts caused by defective cells in x-ray flat panel imagers and the artifact removal using a wavelet-analysis-based algorithm. , 2001, Medical physics.

[12]  Arvid Lundervold,et al.  Noise removal using fourth-order partial differential equation with applications to medical magnetic resonance images in space and time , 2003, IEEE Trans. Image Process..

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

[14]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[15]  S. Lele Compact finite difference schemes with spectral-like resolution , 1992 .

[16]  M. J. Turner,et al.  Digital removal of power frequency artifacts using a Fourier space median filter , 2005, IEEE Signal Processing Letters.

[17]  Soo Yeol Lee,et al.  Removal of ring artifacts in computed tomographic imaging using iterative center weighted median filter , 2010, Comput. Biol. Medicine.

[18]  P. Sweby High Resolution Schemes Using Flux Limiters for Hyperbolic Conservation Laws , 1984 .

[19]  Mirko Boin,et al.  Compensation of ring artefacts in synchrotron tomographic images. , 2006, Optics express.

[20]  B. Kavanagh,et al.  Transmission characteristics of a two dimensional antiscatter grid prototype for CBCT , 2017, Medical physics.

[21]  Truong Q. Nguyen,et al.  An Augmented Lagrangian Method for Total Variation Video Restoration , 2011, IEEE Transactions on Image Processing.

[22]  Dennis N. Assanis,et al.  Evaluation of various high-order-accuracy schemes with and without flux limiters , 1993 .