Fast 4D cone-beam reconstruction using the McKinnon-Bates algorithm with truncation correction and nonlinear filtering

A challenge in using on-board cone beam computed tomography (CBCT) to image lung tumor motion prior to radiation therapy treatment is acquiring and reconstructing high quality 4D images in a sufficiently short time for practical use. For the 1 minute rotation times typical of Linacs, severe view aliasing artifacts, including streaks, are created if a conventional phase-correlated FDK reconstruction is performed. The McKinnon-Bates (MKB) algorithm provides an efficient means of reducing streaks from static tissue but can suffer from low SNR and other artifacts due to data truncation and noise. We have added truncation correction and bilateral nonlinear filtering to the MKB algorithm to reduce streaking and improve image quality. The modified MKB algorithm was implemented on a graphical processing unit (GPU) to maximize efficiency. Results show that a nearly 4x improvement in SNR is obtained compared to the conventional FDK phase-correlated reconstruction and that high quality 4D images with 0.4 second temporal resolution and 1 mm3 isotropic spatial resolution can be reconstructed in less than 20 seconds after data acquisition completes.

[1]  M. V. van Herk,et al.  Respiratory correlated cone beam CT. , 2005, Medical physics.

[2]  J. Hornegger,et al.  Fast GPU-Based CT Reconstruction using the Common Unified Device Architecture (CUDA) , 2007, 2007 IEEE Nuclear Science Symposium Conference Record.

[3]  Guang-Hong Chen,et al.  Prior image constrained compressed sensing (PICCS) , 2008, SPIE BiOS.

[4]  Richard A. Robb,et al.  3-D Reconstruction of the Heart from Few Projections: A Practical Implementation of the McKinnon-Bates Algorithm , 1986, IEEE Transactions on Medical Imaging.

[5]  Roberto Manduchi,et al.  Bilateral filtering for gray and color images , 1998, Sixth International Conference on Computer Vision (IEEE Cat. No.98CH36271).

[6]  Klaus Mueller,et al.  A comparative study of popular interpolation and integration methods for use in computed tomography , 2006, 3rd IEEE International Symposium on Biomedical Imaging: Nano to Macro, 2006..

[7]  Markus Kowarschik,et al.  GPU-accelerated SART reconstruction using the CUDA programming environment , 2009, Medical Imaging.

[8]  Jean-Michel Morel,et al.  A non-local algorithm for image denoising , 2005, 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR'05).

[9]  Guang-Hong Chen,et al.  Streaking artifacts reduction in four-dimensional cone-beam computed tomography. , 2008, Medical physics.

[10]  Graeme C. Mc Kinnon,et al.  Towards Imaging the Beating Heart Usefully with a Conventional CT Scanner , 1981, IEEE Transactions on Biomedical Engineering.

[11]  Lei Xing,et al.  Enhanced 4D cone-beam CT with inter-phase motion model. , 2007, Medical physics.

[12]  Rüdiger Westermann,et al.  Acceleration techniques for GPU-based volume rendering , 2003, IEEE Visualization, 2003. VIS 2003..