A systematic 3-D magnetic particle imaging simulation model for quantitative analysis of reconstruction image quality

Abstract

Background and objectiveMagnetic particle imaging (MPI) is an emerging imaging technology in medical tomography that utilizes the nonlinear magnetization response of superparamagnetic iron oxide (SPIO) particles to determine the in vivo spatial distribution of nanoparticle contrast agents. The reconstruction image quality of MPI is determined by the characteristics of magnetic particles, the setting of the MPI scanner parameters, and the hardware interference of MPI systems. We explore a feasible method to systematically and quickly analyze the impact of these factors on MPI reconstruction image quality. MethodsWe propose a systematic 3-D MPI simulation model. The MPI simulation model has the capability of quickly producing the simulated reconstruction images of a scanned phantom, and quantitative analysis of MPI reconstruction image quality can be achieved by comparing the differences between the input image and output image. These factors are mainly classified as imaging parameters and interference parameters in our model. In order to reduce the computational time of the simulation model, we introduce GPU parallel programming to accelerate the processing of large complex matrix data. For ease of use, we also construct a reliable, high-performance, and open-source 3-D MPI simulation software tool based on our model. The efficiency of our model is evaluated by using OpenMPIData. To demonstrate the capabilities of our model, we conduct simulation experiments using parameters consistent with a real MPI scanner for improving MPI image quality. ResultsThe experimental results show that our simulation model can systematically and quickly evaluate the impact of imaging parameters and interference parameters on MPI reconstruction image quality. ConclusionsWe developed an easy-to-use and open-source 3-D MPI simulation software tool based on our simulation model incorporating all the stages of MPI formation, from signal acquisition to image reconstruction. In the future, our simulation model has potential guiding significance to practical MPI images.