Model-Based Reconstruction for MPI

A prerequisite for system matrix-based image reconstruction in Magnetic Particle Imaging (MPI) is the acquisition of a system matrix that describes the mapping between the MPI tracer and the measured signal. A common method for its acquisition is a time-consuming calibration procedure, during which the scanner is blocked for other uses. For this reason, a model-based approach that allows the system matrix to be obtained by simulation has great appeal. However, an accurate model that describes the magnetization behavior of the tracer, allows the identification of its parameters, and is computationally feasible has not yet been found.

In an ongoing collaboration with Hannes Albers and Tobias Kluth from the University of Bremen we investigate and refine a magnetization model based on the Néel rotation for the magnetic moments of the particles (see Kluth et al., 2019 and Albers et al., 2022). On the one hand, the identification of the parameters of the model is in focus, since these are unknown a priori, on the other hand, measured 2D MPI system matrices are describe with much higher accuracy than the current MPI models. Moreover, we are also interested in the limitations current MPI models in the context of fluid dynamics (see Möddel et al., 2023).

A comparison of the frequency components between a measured and model-based system matrix shows significant differences, as the simple and widely used equilibrium model is not able to fully capture the complex magnetization dynamics.

Project Publications

[134633]
Title: Towards accurate modeling of the multidimensional magnetic particle imaging physics.
Written by: T. Kluth, P. Szwargulski, and T. Knopp
in: <em>New Journal of Physics</em>. 10 (2019).
Volume: <strong>21</strong>. Number: (10),
on pages: 103032
Chapter:
Editor:
Publisher: {IOP} Publishing:
Series:
Address:
Edition:
ISBN:
how published:
Organization:
School:
Institution:
Type:
DOI: 10.1088/1367-2630/ab4938
URL: https://doi.org/10.1088%2F1367-2630%2Fab4938
ARXIVID:
PMID:

[www]

Note: article, openaccess, model-based

Abstract: The image reconstruction problem of the tomographic imaging technique magnetic particle imaging (MPI) requires the solution of a linear inverse problem. One prerequisite for this task is that the imaging operator that describes the mapping between the tomographic image and the measured signal is accurately known. For 2D and 3D excitation patterns, it is common to measure the system matrix in a calibration procedure, that is both, very time consuming and adds noise to the operator. The need for measuring the system matrix is due to the lack of an accurate model that is capable of describing the nanoparticles’ magnetization behavior in the MPI setup. Within this work we exploit a physical model that is based on Néel rotation for large particle ensembles and we find model parameters that describe measured 2D MPI data with much higher precision than state of the art MPI models. With phantom experiments we show that the simulated system matrix can be used for image reconstruction and reduces artifacts due to model-mismatch considerably.