Motion-Compensated Cone Beam CT Using a Conjugate Gradient Least Squares Algorithm and EIT Imaging Motion Data

Cone Beam Computed Tomography (CBCT) is an imaging modality that has been used in image-guided radiation therapy (IGRT). For applications such as lung radiation therapy, CBCT images are greatly affected by the motion artefacts. This is mainly due to low temporal resolution of CBCT. Recently, a dual modality of Electrical Impedance Tomography (EIT) and CBCT has been proposed, in which the high temporal resolution EIT imaging system provides motion data to a motion compensated algebraic reconstruction technique (ART) based CBCT reconstruction software. High computational time associated with ART and indeed other variations of ART make it less practical for real applications. This paper develops a motion-compensated conjugate gradient least squares (CGLS) algorithm for CBCT. A motion-compensated CGLS offers several advantages over ART based methods; including possibilities for: explicit regularisation, rapid convergence, and parallel computations. This paper for the first time demonstrates motion compensated CBCT reconstruction using CGLS and reconstruction results are shown in limited data CBCT considering only a quarter of full data set. The proposed algorithm is tested using simulated motion data in generic motion compensated CBCT as well as measured EIT data in dual EIT-CBCT imaging.

[1]  David A. Jaffray,et al.  Respiration correlated cone-beam computed tomography and 4DCT for evaluating target motion in Stereotactic Lung Radiation Therapy , 2006, Acta oncologica.

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

[3]  Uwe Oelfke,et al.  Linac-integrated 4D cone beam CT: first experimental results , 2006, Physics in medicine and biology.

[4]  K. Lam,et al.  Uncertainties in CT-based radiation therapy treatment planning associated with patient breathing. , 1996, International journal of radiation oncology, biology, physics.

[5]  Per Christian Hansen,et al.  Rank-Deficient and Discrete Ill-Posed Problems , 1996 .

[6]  M. Soleimani,et al.  New iterative cone beam CT reconstruction software: Parameter optimisation and convergence study , 2010, Comput. Methods Programs Biomed..

[7]  Steve B. Jiang,et al.  Effects of intra-fraction motion on IMRT dose delivery: statistical analysis and simulation. , 2002, Physics in medicine and biology.

[8]  J. G. Alessi,et al.  The Algebraic Reconstruction Technique (ART) , 1997, Proceedings of the 1997 Particle Accelerator Conference (Cat. No.97CH36167).

[9]  Ming Jiang,et al.  Convergence of the simultaneous algebraic reconstruction technique (SART) , 2003, Conference Record of Thirty-Fifth Asilomar Conference on Signals, Systems and Computers (Cat.No.01CH37256).

[10]  Simon Rit,et al.  Comparison of Analytic and Algebraic Methods for Motion-Compensated Cone-Beam CT Reconstruction of the Thorax , 2009, IEEE Transactions on Medical Imaging.

[11]  Cristian T. Badea,et al.  CONE BASED 3D RECONSTRUCTION: A FDK - SART COMPARISON FOR LIMITED NUMBER OF PROJECTIONS , 2001 .

[12]  M. Jiang,et al.  Ordered-subset simultaneous algebraic reconstruction techniques (OS-SART) , 2004 .

[13]  M Soleimani,et al.  A motion-compensated cone-beam CT using electrical impedance tomography imaging. , 2011, Physiological measurement.

[14]  M Grass,et al.  Motion-compensated iterative cone-beam CT image reconstruction with adapted blobs as basis functions , 2008, Physics in medicine and biology.

[15]  Åke Björck,et al.  Numerical methods for least square problems , 1996 .

[16]  H. Shirato,et al.  Four-dimensional treatment planning and fluoroscopic real-time tumor tracking radiotherapy for moving tumor. , 2000, International journal of radiation oncology, biology, physics.

[17]  S. Rit,et al.  On-the-fly motion-compensated cone-beam CT using an a priori model of the respiratory motion. , 2009, Medical physics.

[18]  M. Goitein.  Organ and tumor motion: an overview. , 2004, Seminars in radiation oncology.

[19]  Steve B. Jiang,et al.  GPU-based iterative cone-beam CT reconstruction using tight frame regularization , 2010, Physics in medicine and biology.

[20]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[21]  A. H. Andersen.  Algebraic reconstruction in CT from limited views. , 1989, IEEE transactions on medical imaging.

[22]  G. Herman,et al.  Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and x-ray photography. , 1970, Journal of theoretical biology.