Using Low-Rank Tensors for the Recovery of MPI System Matrices

In Magnetic Particle Imaging (MPI), the system matrix plays an important role, as it encodes the relationship between particle concentration and the measured signal. Its acquisition requires a time-consuming calibration scan, whereas its size leads to a high memory-demand. Both of these aspects can be limiting factors in practice. In order to reduce measurement time, compressed sensing exploits the knowledge that the MPI system matrix has a sparse representation in a suitably chosen domain. In this work we demonstrate that the rows of the system matrix allow a representation as low-rank tensors. We show that such an approximation leads to a denoising of the system matrix while introducing only a negligible bias. As an application, we develop a new matrix recovery method exploiting aforementioned low rank property in addition to sparsity in the DCT domain. Experiments show that the proposed matrix recovery method yields system matrices with reduced error when compared to a standard compressed sensing recovery.

[1]  Huazhong Shu,et al.  Domain Progressive 3D Residual Convolution Network to Improve Low-Dose CT Imaging , 2019, IEEE Transactions on Medical Imaging.

[2]  Tobias Knopp,et al.  Discriminating nanoparticle core size using multi-contrast MPI , 2019, Physics in medicine and biology.

[3]  Thorsten M. Buzug,et al.  Artifact free reconstruction with the system matrix approach by overscanning the field-free-point trajectory in magnetic particle imaging , 2016, Physics in medicine and biology.

[4]  B Gleich,et al.  First experimental evidence of the feasibility of multi-color magnetic particle imaging , 2015, Physics in medicine and biology.

[5]  Yoram Bresler,et al.  Transform Learning for Magnetic Resonance Image Reconstruction: From Model-Based Learning to Building Neural Networks , 2019, IEEE Signal Processing Magazine.

[6]  Daniel K Sodickson,et al.  Low‐rank plus sparse matrix decomposition for accelerated dynamic MRI with separation of background and dynamic components , 2015, Magnetic resonance in medicine.

[7]  Jin Liu,et al.  3D Feature Constrained Reconstruction for Low-Dose CT Imaging , 2018, IEEE Transactions on Circuits and Systems for Video Technology.

[8]  Anru Zhang,et al.  Sparse and Low-Rank Tensor Estimation via Cubic Sketchings , 2018, IEEE Transactions on Information Theory.

[9]  Jin Liu,et al.  Discriminative Feature Representation to Improve Projection Data Inconsistency for Low Dose CT Imaging , 2017, IEEE Transactions on Medical Imaging.

[10]  David L. Donoho,et al.  De-noising by soft-thresholding , 1995, IEEE Trans. Inf. Theory.

[11]  Matthias Graeser,et al.  Magnetic Particle Imaging for Real-Time Perfusion Imaging in Acute Stroke. , 2017, ACS nano.

[12]  Tobias Knopp,et al.  Reconstruction of the Magnetic Particle Imaging System Matrix Using Symmetries and Compressed Sensing , 2015 .

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

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

[15]  Guangming Shi,et al.  Compressive Sensing via Nonlocal Low-Rank Regularization , 2014, IEEE Transactions on Image Processing.

[16]  B Gleich,et al.  Three-dimensional real-time in vivo magnetic particle imaging , 2009, Physics in medicine and biology.

[17]  Tobias Knopp,et al.  Sparse Reconstruction of the Magnetic Particle Imaging System Matrix , 2013, IEEE Transactions on Medical Imaging.

[18]  Tobias Knopp,et al.  Magnetic particle imaging for in vivo blood flow velocity measurements in mice , 2018, Physics in medicine and biology.

[19]  Tobias Knopp,et al.  MPIReco.jl: Julia Package for Image Reconstruction in MPI , 2019 .

[20]  Huazhong Shu,et al.  Artifact Suppressed Dictionary Learning for Low-Dose CT Image Processing , 2014, IEEE Transactions on Medical Imaging.

[21]  Jörn Borgert,et al.  Magnetic Particle Imaging (MPI): Experimental Quantification of Vascular Stenosis Using Stationary Stenosis Phantoms , 2017, PloS one.

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

[23]  Tobias Knopp,et al.  OpenMPIData: An initiative for freely accessible magnetic particle imaging data , 2019, Data in brief.

[24]  Jochen Franke,et al.  Magnetic Particle Imaging: A Resovist based Marking Technology for Guide Wires and Catheters for Vascular Interventions. , 2016, IEEE transactions on medical imaging.

[25]  Bernhard Gleich,et al.  Simultaneous magnetic particle imaging (MPI) and temperature mapping using multi-color MPI , 2016 .

[26]  Tobias Knopp,et al.  Magnetic Particle / Magnetic Resonance Imaging: In-Vitro MPI-Guided Real Time Catheter Tracking and 4D Angioplasty Using a Road Map and Blood Pool Tracer Approach , 2016, PloS one.

[27]  Robert H. Halstead,et al.  Matrix Computations , 2011, Encyclopedia of Parallel Computing.

[28]  Tobias Knopp,et al.  Towards accurate modeling of the multidimensional magnetic particle imaging physics , 2019, New Journal of Physics.

[29]  Jong Chul Ye,et al.  ${k}$ -Space Deep Learning for Accelerated MRI , 2020, IEEE Transactions on Medical Imaging.

[30]  Bruce R. Rosen,et al.  Image reconstruction by domain-transform manifold learning , 2017, Nature.

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

[32]  J. Suykens,et al.  Nuclear Norms for Tensors and Their Use for Convex Multilinear Estimation , 2011 .

[33]  Bernhard Gleich,et al.  Analysis of a 3-D System Function Measured for Magnetic Particle Imaging , 2012, IEEE Transactions on Medical Imaging.

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

[35]  Erik G. Larsson,et al.  The Higher-Order Singular Value Decomposition: Theory and an Application [Lecture Notes] , 2010, IEEE Signal Processing Magazine.

[36]  Omer Burak Demirel,et al.  Fully automated gridding reconstruction for non-Cartesian x-space magnetic particle imaging , 2019, Physics in medicine and biology.

[37]  Zhi-Quan Luo,et al.  Guaranteed Matrix Completion via Non-Convex Factorization , 2014, IEEE Transactions on Information Theory.

[38]  Jonathan I. Tamir,et al.  T2 shuffling: Sharp, multicontrast, volumetric fast spin‐echo imaging , 2017, Magnetic resonance in medicine.

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

[40]  Wei Chen,et al.  Nonconvex Robust Low-Rank Tensor Reconstruction via an Empirical Bayes Method , 2019, IEEE Transactions on Signal Processing.

[41]  Tobias Knopp,et al.  Direct Image Reconstruction of Lissajous-Type Magnetic Particle Imaging Data Using Chebyshev-Based Matrix Compression , 2017, IEEE Transactions on Computational Imaging.

[42]  Thorsten M. Buzug,et al.  Model-Based Reconstruction for Magnetic Particle Imaging , 2010, IEEE Transactions on Medical Imaging.

[43]  T Knopp,et al.  Human-sized magnetic particle imaging for brain applications , 2018, Nature Communications.

[44]  Emine Ulku Saritas,et al.  Fast System Calibration With Coded Calibration Scenes for Magnetic Particle Imaging , 2019, IEEE Transactions on Medical Imaging.

[45]  Yipeng Liu,et al.  Smooth robust tensor principal component analysis for compressed sensing of dynamic MRI , 2020, Pattern Recognit..

[46]  Wotao Yin,et al.  Bregman Iterative Algorithms for (cid:2) 1 -Minimization with Applications to Compressed Sensing ∗ , 2008 .

[47]  Mariya Doneva,et al.  Compressed sensing reconstruction for magnetic resonance parameter mapping , 2010, Magnetic resonance in medicine.

[48]  S. Frick,et al.  Compressed Sensing , 2014, Computer Vision, A Reference Guide.

[49]  Zemin Zhang,et al.  Exact Tensor Completion Using t-SVD , 2015, IEEE Transactions on Signal Processing.

[50]  Andrzej Cichocki,et al.  Smooth PARAFAC Decomposition for Tensor Completion , 2015, IEEE Transactions on Signal Processing.

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

[52]  Bernhard Gleich,et al.  Tomographic imaging using the nonlinear response of magnetic particles , 2005, Nature.

[53]  Patrick W. Goodwill,et al.  The X-Space Formulation of the Magnetic Particle Imaging Process: 1-D Signal, Resolution, Bandwidth, SNR, SAR, and Magnetostimulation , 2010, IEEE Transactions on Medical Imaging.

[54]  Yanjun Li,et al.  Image Recovery via Transform Learning and Low-Rank Modeling: The Power of Complementary Regularizers , 2020, IEEE Transactions on Image Processing.

[55]  Alfred Mertins,et al.  On the Representation of Magnetic Particle Imaging in Fourier Space , 2020 .

[56]  Eero P. Simoncelli,et al.  Image quality assessment: from error visibility to structural similarity , 2004, IEEE Transactions on Image Processing.

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

[58]  Martin Möddel,et al.  Viscosity quantification using multi-contrast magnetic particle imaging , 2018, New Journal of Physics.

[59]  B Gleich,et al.  Weighted iterative reconstruction for magnetic particle imaging , 2010, Physics in medicine and biology.

[60]  Thorsten M Buzug,et al.  Toward cardiovascular interventions guided by magnetic particle imaging: First instrument characterization , 2013, Magnetic resonance in medicine.

[61]  Anand Rangarajan,et al.  Image Denoising Using the Higher Order Singular Value Decomposition , 2013, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[62]  Yi Ma,et al.  Robust principal component analysis? , 2009, JACM.

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

[64]  René M. Botnar,et al.  High‐dimensionality undersampled patch‐based reconstruction (HD‐PROST) for accelerated multi‐contrast MRI , 2019, Magnetic resonance in medicine.

[65]  Ce Zhu,et al.  Image Completion Using Low Tensor Tree Rank and Total Variation Minimization , 2019, IEEE Transactions on Multimedia.

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

[67]  Sabine Van Huffel,et al.  Multi-structural Signal Recovery for Biomedical Compressive Sensing , 2013, IEEE Transactions on Biomedical Engineering.

[68]  Stephen P. Boyd,et al.  Distributed Optimization and Statistical Learning via the Alternating Direction Method of Multipliers , 2011, Found. Trends Mach. Learn..

[69]  Jieping Ye,et al.  Tensor Completion for Estimating Missing Values in Visual Data , 2009, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[70]  Bernhard Gleich,et al.  Signal encoding in magnetic particle imaging: properties of the system function , 2009, BMC Medical Imaging.