Higher Order Regularization using Harmonic Eigenfunctions for Model-Based Reconstruction in Magnetic Particle Imaging
Fuente:
arXiv
Gespeichert in:
| Hauptverfasser: | , , |
|---|---|
| Format: | Preprint |
| Veröffentlicht: |
2025
|
| Schlagworte: | |
| Online-Zugang: | |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| _version_ | 1866913980499361792 |
|---|---|
| author | März, Thomas Gapyak, Vladyslav Weinmann, Andreas |
| author_facet | März, Thomas Gapyak, Vladyslav Weinmann, Andreas |
| contents | Magnetic Particle Imaging (MPI) is a recent imaging modality where superparamagnetic nanoparticles are employed as tracers. The reconstruction task is to obtain the spatial particle distribution from a voltage signal induced by the particles. Generally, in computational imaging variational reconstruction techniques are common and rely on a mathematical model to describe the underlying physics. For the MPI reconstruction task we propose a model-based variational reconstruction technique which incorporates a higher order regularizer, where the regularizer is diagonalized by harmonic eigenfunctions. The proposed image reconstruction algorithm features two major stages: in the first stage, the core stage, the components of the MPI core response are reconstructed. This is the MPI-specific data approximation task which we formulate as a variational problem incorporating the higher order regularizer. The relationship between the particle distribution, the MPI core response and the measured data is given by a mathematical model which was introduced in our earlier research. According to this model the MPI core response is tied to the particle distribution by convolution. Therefore the outcome of the core stage yields the data for the second stage, the deconvolution stage, in which the final reconstructed image is produced by solving an ill-posed deconvolution problem in a robust way relying on earlier research. Interestingly, the quality of the final image depends significantly on the quality of the result of the core stage. A contribution is thus the enhancement of the core stage via higher order regularization. We provide a theoretical foundation for our approach and demonstrate its benefit with numerical examples. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2508_06306 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Higher Order Regularization using Harmonic Eigenfunctions for Model-Based Reconstruction in Magnetic Particle Imaging März, Thomas Gapyak, Vladyslav Weinmann, Andreas Numerical Analysis 65K10, 65R32, 65T40, 92C55 Magnetic Particle Imaging (MPI) is a recent imaging modality where superparamagnetic nanoparticles are employed as tracers. The reconstruction task is to obtain the spatial particle distribution from a voltage signal induced by the particles. Generally, in computational imaging variational reconstruction techniques are common and rely on a mathematical model to describe the underlying physics. For the MPI reconstruction task we propose a model-based variational reconstruction technique which incorporates a higher order regularizer, where the regularizer is diagonalized by harmonic eigenfunctions. The proposed image reconstruction algorithm features two major stages: in the first stage, the core stage, the components of the MPI core response are reconstructed. This is the MPI-specific data approximation task which we formulate as a variational problem incorporating the higher order regularizer. The relationship between the particle distribution, the MPI core response and the measured data is given by a mathematical model which was introduced in our earlier research. According to this model the MPI core response is tied to the particle distribution by convolution. Therefore the outcome of the core stage yields the data for the second stage, the deconvolution stage, in which the final reconstructed image is produced by solving an ill-posed deconvolution problem in a robust way relying on earlier research. Interestingly, the quality of the final image depends significantly on the quality of the result of the core stage. A contribution is thus the enhancement of the core stage via higher order regularization. We provide a theoretical foundation for our approach and demonstrate its benefit with numerical examples. |
| title | Higher Order Regularization using Harmonic Eigenfunctions for Model-Based Reconstruction in Magnetic Particle Imaging |
| topic | Numerical Analysis 65K10, 65R32, 65T40, 92C55 |
| url | https://arxiv.org/abs/2508.06306 |