Accelerating dynamic magnetic resonance imaging (MRI) for lung tumor tracking based on low-rank decomposition in the spatial-temporal domain: a feasibility study based on simulation and preliminary prospective undersampled MRI.

PURPOSE To evaluate a low-rank decomposition method to reconstruct down-sampled k-space data for the purpose of tumor tracking. METHODS AND MATERIALS Seven retrospective lung cancer patients were included in the simulation study. The fully-sampled k-space data were first generated from existing 2-dimensional dynamic MR images and then down-sampled by 5 × -20 × before reconstruction using a Cartesian undersampling mask. Two methods, a low-rank decomposition method using combined dynamic MR images (k-t SLR based on sparsity and low-rank penalties) and a total variation (TV) method using individual dynamic MR frames, were used to reconstruct images. The tumor trajectories were derived on the basis of autosegmentation of the resultant images. To further test its feasibility, k-t SLR was used to reconstruct prospective data of a healthy subject. An undersampled balanced steady-state free precession sequence with the same undersampling mask was used to acquire the imaging data. RESULTS In the simulation study, higher imaging fidelity and low noise levels were achieved with the k-t SLR compared with TV. At 10 × undersampling, the k-t SLR method resulted in an average normalized mean square error <0.05, as opposed to 0.23 by using the TV reconstruction on individual frames. Less than 6% showed tracking errors >1 mm with 10 × down-sampling using k-t SLR, as opposed to 17% using TV. In the prospective study, k-t SLR substantially reduced reconstruction artifacts and retained anatomic details. CONCLUSIONS Magnetic resonance reconstruction using k-t SLR on highly undersampled dynamic MR imaging data results in high image quality useful for tumor tracking. The k-t SLR was superior to TV by better exploiting the intrinsic anatomic coherence of the same patient. The feasibility of k-t SLR was demonstrated by prospective imaging acquisition and reconstruction.

[1]  Emmanuel J. Candès,et al.  Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information , 2004, IEEE Transactions on Information Theory.

[2]  A. Majumdar,et al.  An algorithm for sparse MRI reconstruction by Schatten p-norm minimization. , 2011, Magnetic resonance imaging.

[3]  L. Xing,et al.  4D cone beam CT via spatiotemporal tensor framelet. , 2012, Medical physics.

[4]  Cristian Lorenz,et al.  Impact of four-dimensional computed tomography pulmonary ventilation imaging-based functional avoidance for lung cancer radiotherapy. , 2011, International journal of radiation oncology, biology, physics.

[5]  Talissa A Altes,et al.  Evaluation of the reproducibility of lung motion probability distribution function (PDF) using dynamic MRI , 2007, Physics in medicine and biology.

[6]  E. Candès,et al.  Stable signal recovery from incomplete and inaccurate measurements , 2005, math/0503066.

[7]  B. Fallone,et al.  First MR images obtained during megavoltage photon irradiation from a prototype integrated linac-MR system. , 2009, Medical physics.

[8]  J G M Kok,et al.  Towards MRI-guided linear accelerator control: gating on an MRI accelerator , 2011, Physics in medicine and biology.

[9]  Hak Choy,et al.  International Journal of Radiation Oncology Biology Physics: Introduction , 2004 .

[10]  K. Sheng,et al.  The effect of respiratory motion variability and tumor size on the accuracy of average intensity projection from four-dimensional computed tomography: an investigation based on dynamic MRI. , 2008, Medical physics.

[11]  Hans-Ulrich Kauczor,et al.  Measurement of tumor diameter-dependent mobility of lung tumors by dynamic MRI. , 2004, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[12]  Emmanuel J. Candès,et al.  Exact Matrix Completion via Convex Optimization , 2009, Found. Comput. Math..

[13]  J Rottmann,et al.  Markerless EPID image guided dynamic multi-leaf collimator tracking for lung tumors , 2013, Physics in medicine and biology.

[14]  José M. Bioucas-Dias,et al.  Fast Image Recovery Using Variable Splitting and Constrained Optimization , 2009, IEEE Transactions on Image Processing.

[15]  Steve B. Jiang,et al.  Cine Cone Beam CT Reconstruction Using Low-Rank Matrix Factorization: Algorithm and a Proof-of-Principle Study , 2012, IEEE Transactions on Medical Imaging.

[16]  Hui Yan,et al.  A phantom study on the positioning accuracy of the Novalis Body system. , 2003, Medical physics.

[17]  Fang-Fang Yin,et al.  Four-dimensional magnetic resonance imaging (4D-MRI) using image-based respiratory surrogate: a feasibility study. , 2011, Medical physics.

[18]  D. Donoho,et al.  Sparse MRI: The application of compressed sensing for rapid MR imaging , 2007, Magnetic resonance in medicine.

[19]  Ke Sheng,et al.  Reproducibility of interfraction lung motion probability distribution function using dynamic MRI: statistical analysis. , 2008, International journal of radiation oncology, biology, physics.

[20]  Hans-Ulrich Kauczor,et al.  Imaging tumour motion for radiotherapy planning using MRI , 2006, Cancer imaging : the official publication of the International Cancer Imaging Society.

[21]  A. Majumdar.  Improved dynamic MRI reconstruction by exploiting sparsity and rank-deficiency. , 2013, Magnetic resonance imaging.

[22]  Steve B. Jiang,et al.  MRI-guided tumor tracking in lung cancer radiotherapy , 2011, Physics in medicine and biology.

[23]  R. Chartrand,et al.  Restricted isometry properties and nonconvex compressive sensing , 2007 .

[24]  Thomas E. Yankeelov,et al.  Real-Time Compressive Sensing MRI Reconstruction Using GPU Computing and Split Bregman Methods , 2012, Int. J. Biomed. Imaging.

[25]  K. Scheffler,et al.  Analysis and compensation of eddy currents in balanced SSFP , 2005, Magnetic resonance in medicine.

[26]  J G M Kok,et al.  Integrating a 1.5 T MRI scanner with a 6 MV accelerator: proof of concept , 2009, Physics in medicine and biology.

[27]  Angshul Majumdar,et al.  Exploiting rank deficiency and transform domain sparsity for MR image reconstruction. , 2012, Magnetic resonance imaging.

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

[29]  A Z Faranesh,et al.  Ultrafast pulse sequence techniques for cardiac magnetic resonance imaging. , 2000, Topics in magnetic resonance imaging : TMRI.

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

[31]  Mathews Jacob,et al.  Accelerated Dynamic MRI Exploiting Sparsity and Low-Rank Structure: k-t SLR , 2011, IEEE Transactions on Medical Imaging.