A Wavelet Based Sparse Row-Action Method for Image Reconstruction in Magnetic Particle Imaging

PURPOSE Magnetic Particle Imaging (MPI) is a preclinical imaging technique capable of visualizing the spatio-temporal distribution of magnetic nanoparticles. The image reconstruction of this fast and dynamic process relies on efficiently solving an ill-posed inverse problem. Current approaches to reconstruct the tracer concentration from its measurements are either adapted to image characteristics of MPI but suffer from higher computational complexity and slower convergence or are fast but lack in the image quality of the reconstructed images. METHODS In this work we propose a novel MPI reconstruction method to combine the advantages of both approaches into a single algorithm. The underlying sparsity prior is based on an undecimated wavelet transform and is integrated into a fast row-action framework to solve the corresponding MPI minimization problem. RESULTS Its performance is numerically evaluated against a classical FISTA (Fast It-erative Shrinkage-Thresholding Algorithm) approach on simulated and real MPI data. The experimental results show that the proposed method increases image quality at significantly reduced computation times. CONCLUSIONS In comparison to state-of-the-art MPI reconstruction methods, our approach shows better reconstruction results and at the same time accelerates the convergence rate of the underlying row-action algorithm.

[1]  Bangti Jin,et al.  L1 data fitting for robust reconstruction in magnetic particle imaging: quantitative evaluation on Open MPI dataset , 2020, ArXiv.

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

[3]  Tao Zhang,et al.  A fundamental proof of convergence of alternating direction method of multipliers for weakly convex optimization , 2019, Journal of Inequalities and Applications.

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

[5]  T Knopp,et al.  Correction of linear system drifts in magnetic particle imaging , 2019, Physics in medicine and biology.

[6]  Tobias Knopp,et al.  Toward employing the full potential of magnetic particle imaging: exploring visualization techniques and clinical use cases for real-time 3D vascular imaging , 2019, Medical Imaging.

[7]  Bangti Jin,et al.  Enhanced reconstruction in magnetic particle imaging by whitening and randomized SVD approximation , 2019, Physics in medicine and biology.

[8]  Ion Necoara,et al.  Faster Randomized Block Kaczmarz Algorithms , 2019, SIAM J. Matrix Anal. Appl..

[9]  Zhi Wei Tay,et al.  Magnetic Particle Imaging-Guided Heating in Vivo Using Gradient Fields for Arbitrary Localization of Magnetic Hyperthermia Therapy. , 2018, ACS nano.

[10]  Jungwon Yoon,et al.  Real-Time Two-Dimensional Magnetic Particle Imaging for Electromagnetic Navigation in Targeted Drug Delivery , 2017, Sensors.

[11]  Xiaoming Yuan,et al.  Convergence Analysis of Douglas-Rachford Splitting Method for "Strongly + Weakly" Convex Programming , 2017, SIAM J. Numer. Anal..

[12]  T Knopp,et al.  Hybrid system calibration for multidimensional magnetic particle imaging , 2017, Physics in medicine and biology.

[13]  T M Buzug,et al.  Two dimensional magnetic particle spectrometry , 2017, Physics in medicine and biology.

[14]  Tolga Çukur,et al.  Comparison of System-Matrix-Based and Projection-Based Reconstructions for Field Free Line Magnetic Particle Imaging , 2017 .

[15]  Tobias Kluth,et al.  Improved image reconstruction in magnetic particle imaging using structural a priori information , 2017 .

[16]  Patrick W. Goodwill,et al.  Magnetic Particle Imaging: A Novel in Vivo Imaging Platform for Cancer Detection. , 2017, Nano letters.

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

[18]  Dirk A. Lorenz,et al.  Linear convergence of the randomized sparse Kaczmarz method , 2016, Mathematical Programming.

[19]  T Knopp,et al.  Online reconstruction of 3D magnetic particle imaging data , 2016, Physics in medicine and biology.

[20]  I. Bayram,et al.  Correction for "On the Convergence of the Iterative Shrinkage/Thresholding Algorithm with a Weakly Convex Penalty" , 2016, IEEE Trans. Signal Process..

[21]  Matthias Graeser,et al.  Magnetic particle imaging: current developments and future directions , 2015, International journal of nanomedicine.

[22]  Matthieu Kowalski,et al.  Thresholding RULES and iterative shrinkage/thresholding algorithm: A convergence study , 2014, 2014 IEEE International Conference on Image Processing (ICIP).

[23]  D. Needell,et al.  Randomized block Kaczmarz method with projection for solving least squares , 2014, 1403.4192.

[24]  Stephen J. Wright,et al.  An Asynchronous Parallel Randomized Kaczmarz Algorithm , 2014, ArXiv.

[25]  Dirk A. Lorenz,et al.  The Linearized Bregman Method via Split Feasibility Problems: Analysis and Generalizations , 2013, SIAM J. Imaging Sci..

[26]  Neal Parikh,et al.  Distributed Optimization and Statistical Learning via the Alternating Direction Method of Multipliers , 2011, Found. Trends Mach. Learn..

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

[28]  Patrick L. Combettes,et al.  Proximal Splitting Methods in Signal Processing , 2009, Fixed-Point Algorithms for Inverse Problems in Science and Engineering.

[29]  Jean-Luc Starck,et al.  Monotone operator splitting for optimization problems in sparse recovery , 2009, 2009 16th IEEE International Conference on Image Processing (ICIP).

[30]  B Gleich,et al.  Trajectory analysis for magnetic particle imaging , 2008, Physics in medicine and biology.

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

[32]  Otmar Scherzer,et al.  KACZMARZ METHODS FOR REGULARIZING NONLINEAR ILL-POSED EQUATIONS I: CONVERGENCE ANALYSIS , 2007 .

[33]  R. Vershynin,et al.  A Randomized Kaczmarz Algorithm with Exponential Convergence , 2007, math/0702226.

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

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

[36]  Per Christian Hansen,et al.  REGULARIZATION TOOLS: A Matlab package for analysis and solution of discrete ill-posed problems , 1994, Numerical Algorithms.

[37]  Stéphane Mallat,et al.  A Theory for Multiresolution Signal Decomposition: The Wavelet Representation , 1989, IEEE Trans. Pattern Anal. Mach. Intell..

[38]  H. Trussell,et al.  The Landweber iteration and projection onto convex sets , 1985, IEEE Trans. Acoust. Speech Signal Process..

[39]  C. Kaethner,et al.  Magnetic Particle Imaging , 2018 .

[40]  T. Knopp,et al.  Edge Preserving and Noise Reducing Reconstruction for Magnetic Particle Imaging , 2017, IEEE Transactions on Medical Imaging.

[41]  Xuemei Chen Kaczmarz Algorithm, row action methods, and statistical learning algorithms , 2017 .

[42]  Marc Teboulle,et al.  A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems , 2009, SIAM J. Imaging Sci..

[43]  Robert Tibshirani,et al.  The Elements of Statistical Learning: Data Mining, Inference, and Prediction , 2001, Springer Series in Statistics.

[44]  S. Mallat A wavelet tour of signal processing , 1998 .

[45]  B. Silverman,et al.  The Stationary Wavelet Transform and some Statistical Applications , 1995 .