Non-Equispaced System Matrix Acquisition for Magnetic Particle Imaging Based on Lissajous Node Points

Magnetic Particle Imaging (MPI) is an emerging technology in the field of (pre)clinical imaging. The acquisition of a particle signal is realized along specific sampling trajectories covering a defined field of view (FOV). In a system matrix (SM) based reconstruction procedure, the commonly used acquisition path in MPI is a Lissajous trajectory. Such a trajectory features an inhomogeneous coverage of the FOV, i.e. the center region is sampled less dense than the regions towards the edges of the FOV. Conventionally, the respective SM acquisition and the subsequent reconstruction do not reflect this inhomogeneous coverage. Instead, they are performed on an equispaced grid. The objective of this work is to introduce a sampling grid that inherently features the aforementioned inhomogeneity by using node points of Lissajous trajectories. Paired with a tailored polynomial interpolation of the reconstructed MPI signal, the entire image can be recovered. It is the first time that such a trajectory related non-equispaced grid is used for image reconstruction on simulated and measured MPI data and it is shown that the number of sampling positions can be reduced, while the spatial resolution remains constant.

[1]  Achiya Dax,et al.  On Row Relaxation Methods for Large Constrained Least Squares Problems , 1993, SIAM J. Sci. Comput..

[2]  Jochen Franke,et al.  On the formulation of the image reconstruction problem in magnetic particle imaging , 2013, Biomedizinische Technik. Biomedical engineering.

[3]  Marco Vianello,et al.  Trivariate Polynomial Approximation on Lissajous Curves , 2015 .

[4]  David P. Dobkin,et al.  The quickhull algorithm for convex hulls , 1996, TOMS.

[5]  Justin J. Konkle,et al.  Projection Reconstruction Magnetic Particle Imaging , 2013, IEEE Transactions on Medical Imaging.

[6]  Thorsten M. Buzug,et al.  Reconstruction Enhancement by Denoising the Magnetic Particle Imaging System Matrix Using Frequency Domain Filter , 2015, IEEE Transactions on Magnetics.

[7]  Thorsten M. Buzug,et al.  Undersampling the system matrix of a single sided MPI-scanner , 2015, 2015 5th International Workshop on Magnetic Particle Imaging (IWMPI).

[8]  D. Sandwell.  BIHARMONIC SPLINE INTERPOLATION OF GEOS-3 AND SEASAT ALTIMETER DATA , 1987 .

[9]  Bernhard Gleich,et al.  Perspectives on clinical magnetic particle imaging , 2013, Biomedizinische Technik. Biomedical engineering.

[10]  Christian Kaethner,et al.  A survey on bivariate Lagrange interpolation on Lissajous nodes , 2015 .

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

[12]  B Gleich,et al.  Fast reconstruction in magnetic particle imaging , 2012, Physics in medicine and biology.

[13]  Volker Rasche,et al.  Resampling of data between arbitrary grids using convolution interpolation , 1999, IEEE Transactions on Medical Imaging.

[14]  Thorsten M. Buzug,et al.  Bivariate Lagrange interpolation at the node points of non-degenerate Lissajous curves , 2014, Numerische Mathematik.

[15]  Patrick Vogel,et al.  Traveling Wave Magnetic Particle Imaging , 2014, IEEE Transactions on Medical Imaging.

[16]  Mark A Griswold,et al.  Magnetic Particle Imaging Tracers: State-of-the-Art and Future Directions. , 2015, The journal of physical chemistry letters.

[17]  Georges Voronoi.  Nouvelles applications des paramètres continus à la théorie des formes quadratiques. Deuxième mémoire. Recherches sur les parallélloèdres primitifs. , 1908 .

[18]  Marco Vianello,et al.  Bivariate Lagrange interpolation at the Padua points: The generating curve approach , 2006, J. Approx. Theory.

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

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

[21]  Wolfgang Erb.  Bivariate Lagrange interpolation at the node points of Lissajous curves - the degenerate case , 2016, Appl. Math. Comput..

[22]  Thorsten M. Buzug,et al.  Compressed Sensing of the System Matrix and Sparse Reconstruction of the Particle Concentration in Magnetic Particle Imaging , 2015, IEEE Transactions on Magnetics.

[23]  T Knopp,et al.  Joint reconstruction of non-overlapping magnetic particle imaging focus-field data , 2015, Physics in medicine and biology.

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

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

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

[27]  P. Hansen.  Discrete Inverse Problems: Insight and Algorithms , 2010 .

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

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

[30]  Patrick W. Goodwill,et al.  Multidimensional X-Space Magnetic Particle Imaging , 2011, IEEE Transactions on Medical Imaging.

[31]  Bo Zheng,et al.  Projection X-Space Magnetic Particle Imaging , 2012, IEEE Transactions on Medical Imaging.

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

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

[34]  Thorsten M. Buzug,et al.  Efficient gradient field generation providing a multi-dimensional arbitrary shifted field-free point for magnetic particle imaging , 2014 .

[35]  Thorsten M. Buzug,et al.  Axially Elongated Field-Free Point Data Acquisition in Magnetic Particle Imaging , 2015, IEEE Transactions on Medical Imaging.

[36]  B Gleich,et al.  Experimental results on fast 2D-encoded magnetic particle imaging , 2008, Physics in medicine and biology.

[37]  B Gleich,et al.  A simulation study on the resolution and sensitivity of magnetic particle imaging , 2007, Physics in medicine and biology.

[38]  Marco Vianello,et al.  Bivariate Lagrange interpolation at the Padua points: Computational aspects , 2008 .