Vladyslav Gapyak

Vladyslav Gapyak, Thomas März, Andreas Weinmann

Model-based reconstruction approaches for the medical imaging modality Magnetic Particle Imaging (MPI) are typically based on the Langevin model, which assumes instantaneous alignment of the particles magnetic momenta with the applied field. Regarding the application to real data, Langevin model-based reconstruction methods require model transfer functions (MTF) obtained from calibrations to preprocess the data. There are also model-based reconstruction approaches that include relaxation effects and other particle-level dynamics. However, they are limited either to 1D or 1D-like scanning scenarios when considering real data, or are limited to simulated data in the case of multi-dimensional field-free point (FFP) MPI. Thus, fully model-based reconstructions from multi-dimensional FFP scanning data that incorporate relaxation effects without using an MTF have not yet been demonstrated. In this work, we incorporate relaxation effects by considering a multi-dimensional Debye model and provide reconstruction formulae. In particular, we show that the Debye model-based signal is the response of a linear time-invariant system with exponential memory applied to a Langevin model-based signal. We provide a reconstruction algorithm for the introduced multi-dimensional Debye model. To this end, we devise a relaxation adaption step. For the resulting relaxation-adapted Debye signal, we show that it can be expressed by the well-studied MPI core operator derived from the Langevin theory. This results in a three-stage algorithm with low additional cost over the Langevin model, as the relaxation adaption scales linearly in the input data. We provide numerical results for the proposed algorithmic approach. In particular, we obtain fully model-based reconstructions from real 2D MPI data without involving any specific MTF analogous to the Langevin model case.

Vladyslav Gapyak, Corinna Rentschler, Thomas März et al.

Objective: Magnetic particle imaging (MPI) is an emerging medical imaging modality which has gained increasing interest in recent years. Among the benefits of MPI are its high temporal resolution, and that the technique does not expose the specimen to any kind of ionizing radiation. It is based on the non-linear response of magnetic nanoparticles to an applied magnetic field. From the electric signal measured in receive coils, the particle concentration has to be reconstructed. Due to the ill-posedness of the reconstruction problem, various regularization methods have been proposed for reconstruction ranging from early stopping methods, via classical Tikhonov regularization and iterative methods to modern machine learning approaches. In this work, we contribute to the latter class: we propose a plug-and-play approach based on a generic zero-shot denoiser with an $\ell^1$-prior. Approach: We validate the reconstruction parameters of the method on a hybrid dataset and compare it with the baseline Tikhonov, DIP and the previous PP-MPI, which is a plug-and-play method with denoiser trained on MPI-friendly data. Main results: We offer a quantitative and qualitative evaluation of the zero-shot plug-and-play approach on the 3D Open MPI dataset. Moreover, we show the quality of the approach with different levels of preprocessing of the data. Significance: The proposed method employs a zero-shot denoiser which has not been trained for the MPI task and therefore saves the cost for training. Moreover, it offers a method that can be potentially applied in future MPI contexts.

Vladyslav Gapyak, Thomas März, Andreas Weinmann

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.