Fast Trajectory-Independent Model-Based Reconstruction Algorithm for Multi-Dimensional Magnetic Particle Imaging
Journal:
arXiv
Published Date:
May 28, 2025
Abstract
Magnetic Particle Imaging (MPI) is a promising tomographic technique for
visualizing the spatio-temporal distribution of superparamagnetic
nanoparticles, with applications ranging from cancer detection to real-time
cardiovascular monitoring. Traditional MPI reconstruction relies on either
time-consuming calibration (measured system matrix) or model-based simulation
of the forward operator. Recent developments have shown the applicability of
Chebyshev polynomials to multi-dimensional Lissajous Field-Free Point (FFP)
scans. This method is bound to the particular choice of sinusoidal scanning
trajectories. In this paper, we present the first reconstruction on real 2D MPI
data with a trajectory-independent model-based MPI reconstruction algorithm. We
further develop the zero-shot Plug-and-Play (PnP) algorithm of the authors --
with automatic noise level estimation -- to address the present deconvolution
problem, leveraging a state-of-the-art denoiser trained on natural images
without retraining on MPI-specific data. We evaluate our method on the publicly
available 2D FFP MPI dataset ``MPIdata: Equilibrium Model with Anisotropy",
featuring scans of six phantoms acquired using a Bruker preclinical scanner.
Moreover, we show reconstruction performed on custom data on a 2D scanner with
additional high-frequency excitation field and partial data. Our results
demonstrate strong reconstruction capabilities across different scanning
scenarios -- setting a precedent for general-purpose, flexible model-based MPI
reconstruction.