L1 data fitting for robust reconstruction in magnetic particle imaging: quantitative evaluation on Open MPI dataset
share
Magnetic particle imaging is an emerging quantitative imaging modality, exploiting the unique nonlinear magnetization phenomenon of superparamagnetic iron oxide nanoparticles for recovering the concentration. Traditionally the reconstruction is formulated into a penalized least-squares problem with nonnegativity constraint, and then solved using a variant of Kaczmarz method. In order to achieve good performance, a preprocessing step of frequency selection to remove the deleterious influences of highly noisy measurements is often adopted. In this work, we propose a complementary approach to frequency selection, by viewing highly noisy measurements as outliers, and employing the l1 data fitting, one popular approach from robust statistics. When compared with the standard approach, it is easy to implement with a comparable computational complexity. Experiments with a public domain dataset, i.e., Open MPI dataset, show that it can give accurate reconstructions, and is less prone to noisy measurements, which is clearly illustrated by the quantitative (PSNR / SSIM) and qualitative comparisons with the standard approach.
READ FULL TEXT