Dynamic MR Image Reconstruction From Highly Undersampled (k, t)-Space Data Exploiting Low Tensor Train Rank and Sparse Prior

Dynamic magnetic resonance imaging (dynamic MRI) is used to visualize living tissues and their changes over time. In this paper, we propose a new tensor-based dynamic MRI approach for reconstruction from highly undersampled (k, <inline-formula> <tex-math notation="LaTeX">$t$ </tex-math></inline-formula>)-space data, which combines low tensor train rankness and temporal sparsity constraints. Considering tensor train (TT) decomposition has superior performance in dealing with high-dimensional tensors, we introduce TT decomposition and utilize the low rankness of TT matrices to exploit the inner structural prior information of dynamic MRI data. First, ket augmentation (KA) scheme is used to permute the 3-order (k, <inline-formula> <tex-math notation="LaTeX">$t$ </tex-math></inline-formula>)-space data to a high order tensor and low rankness of each TT matrix is enforced with different weights. To reduce the computational complexity, we replace the nuclear norm of TT matrices with the minimum Frobenius norm of two factorization matrices to avoid singular value decomposition. Secondly, the <inline-formula> <tex-math notation="LaTeX">$l_{1}$ </tex-math></inline-formula> norm of the Fourier coefficients along the temporal dimension is added as a sparsity constraint to further improve the reconstruction. Lastly, an effective algorithm based on the alternating direction method of multipliers (ADMM) is developed to solve the proposed optimization problem. Numerous experiments have been conducted on three dynamic MRI data sets to estimate the performance of our proposed method. The experimental results and comparisons with several state-of-the-art imaging methods demonstrate the superior performance of the proposed method.

[1]  Zhi-Pei Liang,et al.  SPATIOTEMPORAL IMAGINGWITH PARTIALLY SEPARABLE FUNCTIONS , 2007, 2007 4th IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[2]  Tamara G. Kolda,et al.  Tensor Decompositions and Applications , 2009, SIAM Rev..

[3]  Craig H Meyer,et al.  Block LOw-rank Sparsity with Motion guidance ( BLOSM ) for accelerated dynamic MRI , 2012 .

[4]  V. Aggarwal,et al.  Efficient Low Rank Tensor Ring Completion , 2017, 2017 IEEE International Conference on Computer Vision (ICCV).

[5]  Volkan Cevher,et al.  An SVD-free Approach to a Class of Structured Low Rank Matrix Optimization Problems with Application to System Identification , 2013, CDC 2013.

[6]  Di Guo,et al.  A Fast Low Rank Hankel Matrix Factorization Reconstruction Method for Non-Uniformly Sampled Magnetic Resonance Spectroscopy , 2017, IEEE Access.

[7]  Michael Lustig,et al.  k-t SPARSE: High frame rate dynamic MRI exploiting spatio-temporal sparsity , 2006 .

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

[9]  Felix J. Herrmann,et al.  Optimization on the Hierarchical Tucker manifold – Applications to tensor completion , 2014, Linear Algebra and its Applications.

[10]  Misha Elena Kilmer,et al.  Third-Order Tensors as Operators on Matrices: A Theoretical and Computational Framework with Applications in Imaging , 2013, SIAM J. Matrix Anal. Appl..

[11]  E. Tyrtyshnikov,et al.  TT-cross approximation for multidimensional arrays , 2010 .

[12]  J. Ballani,et al.  Black box approximation of tensors in hierarchical Tucker format , 2013 .

[13]  Yang Wang,et al.  Real-time dynamic MRI using parallel dictionary learning and dynamic total variation , 2017, Neurocomputing.

[14]  Justin P. Haldar,et al.  Low-rank approximations for dynamic imaging , 2011, 2011 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[15]  Junzhou Huang,et al.  An efficient algorithm for dynamic MRI using low‐rank and total variation regularizations , 2018, Medical Image Anal..

[16]  Justin P. Haldar,et al.  Low rank matrix recovery for real-time cardiac MRI , 2010, 2010 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[17]  Yin Zhang,et al.  Solving a low-rank factorization model for matrix completion by a nonlinear successive over-relaxation algorithm , 2012, Mathematical Programming Computation.

[18]  Mathews Jacob,et al.  Recovery of Damped Exponentials Using Structured Low Rank Matrix Completion , 2017, IEEE Transactions on Medical Imaging.

[19]  Zhi-Pei Liang,et al.  Accelerated High-Dimensional MR Imaging With Sparse Sampling Using Low-Rank Tensors , 2016, IEEE Transactions on Medical Imaging.

[20]  Ivan Oseledets,et al.  Tensor-Train Decomposition , 2011, SIAM J. Sci. Comput..

[21]  Pablo A. Parrilo,et al.  Guaranteed Minimum-Rank Solutions of Linear Matrix Equations via Nuclear Norm Minimization , 2007, SIAM Rev..

[22]  A. Manduca,et al.  Local versus Global Low-Rank Promotion in Dynamic MRI Series Reconstruction , 2010 .

[23]  Huiqian Du,et al.  Dynamic MRI Reconstruction Using Tensor-SVD , 2018, 2018 14th IEEE International Conference on Signal Processing (ICSP).

[24]  Eric L. Miller,et al.  Tensor-Based Formulation and Nuclear Norm Regularization for Multienergy Computed Tomography , 2013, IEEE Transactions on Image Processing.

[25]  Ashraf A. Kassim,et al.  Dynamic MRI reconstruction using low rank plus sparse tensor decomposition , 2016, 2016 IEEE International Conference on Image Processing (ICIP).

[26]  Emmanuel J. Candès,et al.  A Singular Value Thresholding Algorithm for Matrix Completion , 2008, SIAM J. Optim..

[27]  Joos Vandewalle,et al.  A Multilinear Singular Value Decomposition , 2000, SIAM J. Matrix Anal. Appl..

[28]  M. Kilmer,et al.  Factorization strategies for third-order tensors , 2011 .

[29]  Jong Chul Ye,et al.  Improved k–t BLAST and k–t SENSE using FOCUSS , 2007, Physics in medicine and biology.

[30]  Zhi-Pei Liang,et al.  Spatiotemporal Imaging with Partially Separable Functions , 2007, 2007 Joint Meeting of the 6th International Symposium on Noninvasive Functional Source Imaging of the Brain and Heart and the International Conference on Functional Biomedical Imaging.

[31]  Jong Chul Ye,et al.  Motion Adaptive Patch-Based Low-Rank Approach for Compressed Sensing Cardiac Cine MRI , 2014, IEEE Transactions on Medical Imaging.

[32]  Nathan Srebro,et al.  Learning with matrix factorizations , 2004 .

[33]  Marko Filipovic,et al.  Tucker factorization with missing data with application to low-$$n$$n-rank tensor completion , 2015, Multidimens. Syst. Signal Process..

[34]  Ivan V. Oseledets,et al.  Approximation of 2d˟2d Matrices Using Tensor Decomposition , 2010, SIAM J. Matrix Anal. Appl..

[35]  Liqing Zhang,et al.  Tensor Ring Decomposition , 2016, ArXiv.

[36]  Eugene E. Tyrtyshnikov,et al.  Breaking the Curse of Dimensionality, Or How to Use SVD in Many Dimensions , 2009, SIAM J. Sci. Comput..

[37]  Ce Zhu,et al.  Low Rank Tensor Completion for Multiway Visual Data , 2018, 1805.03967.

[38]  Justin P. Haldar,et al.  Further development of image reconstruction from highly undersampled (k, t)-space data with joint partial separability and sparsity constraints , 2011, 2011 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[39]  Bing Chen,et al.  Hybrid CS-DMRI: Periodic Time-Variant Subsampling and Omnidirectional Total Variation Based Reconstruction , 2017, IEEE Transactions on Medical Imaging.

[40]  Vin de Silva,et al.  Tensor rank and the ill-posedness of the best low-rank approximation problem , 2006, math/0607647.

[41]  José Ignacio Latorre,et al.  Image compression and entanglement , 2005, ArXiv.

[42]  M. Lustig,et al.  Fast dynamic 3D MR spectroscopic imaging with compressed sensing and multiband excitation pulses for hyperpolarized 13C studies , 2011, Magnetic resonance in medicine.

[43]  Minh N. Do,et al.  Efficient Tensor Completion for Color Image and Video Recovery: Low-Rank Tensor Train , 2016, IEEE Transactions on Image Processing.

[44]  Jong Chul Ye,et al.  A General Framework for Compressed Sensing and Parallel MRI Using Annihilating Filter Based Low-Rank Hankel Matrix , 2015, IEEE Transactions on Computational Imaging.

[45]  David Atkinson,et al.  Dynamic MR Image Reconstruction–Separation From Undersampled (${\bf k},t$)-Space via Low-Rank Plus Sparse Prior , 2014, IEEE Transactions on Medical Imaging.

[46]  Zhi-Pei Liang,et al.  High-Resolution Cardiovascular MRI by Integrating Parallel Imaging With Low-Rank and Sparse Modeling , 2013, IEEE Transactions on Biomedical Engineering.