Deblurring in iterative reconstruction of half CBCT for image guided brain radiosurgery

A high spatial resolution iterative reconstruction algorithm is proposed for a half cone beam CT (HCBCT) geometry. The proposed algorithm improves spatial resolution by explicitly accounting for image blurriness caused by different factors, such as extended X-ray source and detector response. The blurring kernel is estimated using the MTF slice of the Catphan phantom. The proposed algorithm is specifically optimized for the new Leksell Gamma Knife Icon (Elekta AB, Stockholm, Sweden) which incorporates the HCBCT geometry to accommodate the existing treatment couch while covering down to the base-of-skull in the reconstructed field-of-view. Image reconstruction involves a Fourier-based scaling simultaneous algebraic reconstruction technique (SART) coupled with total variation (TV) minimization and non-local mean denoising, solved using a split Bregman separation technique that splits the reconstruction problem into a gradient based updating step and a TV-based deconvolution algorithm. This formulation preserves edges and reduces the staircase effect caused by regular TV-penalized iterative algorithms. Our experiments indicate that our proposed method outperforms the conventional filtered back projection and TV penalized SART methods in terms of line pair resolution and retains the favorable properties of the standard TV-penalized reconstruction.

[1]  Lei Xing,et al.  A Fourier-based compressed sensing technique for accelerated CT image reconstruction using first-order methods , 2014, Physics in medicine and biology.

[2]  Dong Yang,et al.  FDK Half-Scan with a Heuristic Weighting Scheme on a Flat Panel Detector-Based Cone Beam CT (FDKHSCW) , 2006, Int. J. Biomed. Imaging.

[3]  J Webster Stayman,et al.  Generalized least-squares CT reconstruction with detector blur and correlated noise models , 2014, Medical Imaging.

[4]  Pierrick Coupé,et al.  Author manuscript, published in "Journal of Magnetic Resonance Imaging 2010;31(1):192-203" DOI: 10.1002/jmri.22003 Adaptive Non-Local Means Denoising of MR Images with Spatially Varying Noise Levels , 2010 .

[5]  José M. Bioucas-Dias,et al.  A New TwIST: Two-Step Iterative Shrinkage/Thresholding Algorithms for Image Restoration , 2007, IEEE Transactions on Image Processing.

[6]  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.

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

[8]  L. Xing,et al.  Iterative image reconstruction for CBCT using edge-preserving prior. , 2008, Medical physics.

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

[10]  Lei Dong,et al.  Reducing metal artifacts in cone-beam CT images by preprocessing projection data. , 2007, International journal of radiation oncology, biology, physics.

[11]  David L Donoho,et al.  Compressed sensing , 2006, IEEE Transactions on Information Theory.

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

[13]  Luo Ouyang,et al.  Effects of the penalty on the penalized weighted least-squares image reconstruction for low-dose CBCT. , 2011, Physics in medicine and biology.

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

[15]  Jin Sung Kim,et al.  Fast compressed sensing-based CBCT reconstruction using Barzilai-Borwein formulation for application to on-line IGRT. , 2012, Medical physics.