The Reconstruction of Magnetic Particle Imaging: Current Approaches Based on the System Matrix

A greedy regularized block Kaczmarz method for accelerating reconstruction in magnetic particle imaging

Objective.Magnetic particle imaging (MPI) is an emerging medical tomographic imaging modality that enables real-time imaging with high sensitivity and high spatial and temporal resolution. For the system matrix reconstruction method, the MPI reconstruction problem is an ill-posed inverse problem that is commonly solved using the Kaczmarz algorithm. However, the high computation time of the Kaczmarz algorithm, which restricts MPI reconstruction speed, has limited the development of potential clinical applications for real-time MPI. In order to achieve fast reconstruction in real-time MPI, we propose a greedy regularized block Kaczmarz method (GRBK) which accelerates MPI reconstruction.Approach.GRBK is composed of a greedy partition strategy for the system matrix, which enables preprocessing of the system matrix into well-conditioned blocks to facilitate the convergence of the block Kaczmarz algorithm, and a regularized block Kaczmarz algorithm, which enables fast and accurate MPI image reconstruction at the same time.Main results.We quantitatively evaluated our GRBK using simulation data from three phantoms at 20 dB, 30 dB, and 40 dB noise levels. The results showed that GRBK can improve reconstruction speed by single orders of magnitude compared to the prevalent regularized Kaczmarz algorithm including Tikhonov regularization, the non-negative Fused Lasso, and wavelet-based sparse model. We also evaluated our method on OpenMPIData, which is real MPI data. The results showed that our GRBK is better suited for real-time MPI reconstruction than current state-of-the-art reconstruction algorithms in terms of reconstruction speed as well as image quality.Significance.Our proposed method is expected to be the preferred choice for potential applications of real-time MPI.

GPU-accelerated parallel image reconstruction strategies for magnetic particle imaging

Objective. Image reconstruction is a fundamental step in magnetic particle imaging (MPI). One of the main challenges is the fact that the reconstructions are computationally intensive and time-consuming, so choosing an algorithm presents a compromise between accuracy and execution time, which depends on the application. This work proposes a method that provides both fast and accurate image reconstructions.Approach. Image reconstruction algorithms were implemented to be executed in parallel ingraphics processing units(GPUs) using the CUDA framework. The calculation of the model-based MPI calibration matrix was also implemented in GPU to allow both fast and flexible reconstructions.Main results. The parallel algorithms were able to accelerate the reconstructions by up to about6,100times in comparison to the serial Kaczmarz algorithm executed in the CPU, allowing for real-time applications. Reconstructions using the OpenMPIData dataset validated the proposed algorithms and demonstrated that they are able to provide both fast and accurate reconstructions. The calculation of the calibration matrix was accelerated by up to about 37 times.Significance. The parallel algorithms proposed in this work can provide single-frame MPI reconstructions in real time, with frame rates greater than 100 frames per second. The parallel calculation of the calibration matrix can be combined with the parallel reconstruction to deliver images in less time than the serial Kaczmarz reconstruction, potentially eliminating the need of storing the calibration matrix in the main memory, and providing the flexibility of redefining scanning and reconstruction parameters during execution.

A fused LASSO operator for fast 3D magnetic particle imaging reconstruction

Objective.Magnetic particle imaging (MPI) is a promising imaging modality that leverages the nonlinear magnetization behavior of superparamagnetic iron oxide nanoparticles to determine their concentration distribution. Previous optimization models with multiple regularization terms have been proposed to achieve high-quality MPI reconstruction, but these models often result in increased computational burden, particularly for dense gridding 3D fields of view. In order to achieve faster reconstruction speeds without compromising reconstruction quality, we have developed a novel fused LASSO operator, total sum-difference (TSD), which effectively captures the sparse and smooth priors of MPI images.Methods.Through an analysis-synthesis equivalence strategy and a constraint smoothing strategy, the TSD regularized model was solved using the fast iterative soft-thresholding algorithm (FISTA). The resulting reconstruction method, TSD-FISTA, boasts low computational complexity and quadratic convergence rate over iterations.Results.Experimental results demonstrated that TSD-FISTA required only 10% and 37% of the time to achieve comparable or superior reconstruction quality compared to commonly used fused LASSO-based alternating direction method of multipliers and Tikhonov-based algebraic reconstruction techniques, respectively.Significance.TSD-FISTA shows promise for enabling real-time 3D MPI reconstruction at high frame rates for large fields of view.

Machine Learning and Deep Learning Applications in Magnetic Particle Imaging

In recent years, magnetic particle imaging (MPI) has emerged as a promising imaging technique depicting high sensitivity and spatial resolution. It originated in the early 2000s where it proposed a new approach to challenge the low spatial resolution achieved by using relaxometry in order to measure the magnetic fields. MPI presents 2D and 3D images with high temporal resolution, non-ionizing radiation, and optimal visual contrast due to its lack of background tissue signal. Traditionally, the images were reconstructed by the conversion of signal from the induced voltage by generating system matrix and X-space based methods. Because image reconstruction and analyses play an integral role in obtaining precise information from MPI signals, newer artificial intelligence-based methods are continuously being researched and developed upon. In this work, we summarize and review the significance and employment of machine learning and deep learning models for applications with MPI and the potential they hold for the future. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY: Stage 1.

Recent Metal Nanotheranostics for Cancer Diagnosis and Therapy: A Review

In recent years, there has been an increasing interest in using nanoparticles in the medical sciences. Today, metal nanoparticles have many applications in medicine for tumor visualization, drug delivery, and early diagnosis, with different modalities such as X-ray imaging, computed tomography (CT), magnetic resonance imaging (MRI), positron emission tomography (PET), etc., and treatment with radiation. This paper reviews recent findings of recent metal nanotheranostics in medical imaging and therapy. The study offers some critical insights into using different types of metal nanoparticles in medicine for cancer detection and treatment purposes. The data of this review study were gathered from multiple scientific citation websites such as Google Scholar, PubMed, Scopus, and Web of Science up through the end of January 2023. In the literature, many metal nanoparticles are used for medical applications. However, due to their high abundance, low price, and high performance for visualization and treatment, nanoparticles such as gold, bismuth, tungsten, tantalum, ytterbium, gadolinium, silver, iron, platinum, and lead have been investigated in this review study. This paper has highlighted the importance of gold, gadolinium, and iron-based metal nanoparticles in different forms for tumor visualization and treatment in medical applications due to their ease of functionalization, low toxicity, and superior biocompatibility.

Gradient-Based Pulsed Excitation and Relaxation Encoding in Magnetic Particle Imaging

Magnetic particle imaging (MPI) is a radiation-free vessel- and target-imaging modality that can sensitively detect nanoparticles. A static magnetic gradient field, referred to as a selection field, is required in MPI to provide a field-free region (FFR) for spatial encoding. The image resolution of MPI is closely related to the size of the FFR, which is determined by the selection field gradient amplitude. Because of the limitations of existing gradient coil hardware, the image resolution of MPI cannot satisfy the clinical requirements of human in vivo imaging. Pulsed excitation has been confirmed to improve the image resolution of MPI by breaking down the 'relaxation wall.' This work proposes the use of a pulsed waveform magnetic gradient from magnetic resonance imaging to further improve the image resolution of MPI. Through alignment of the gradient direction along the field-free line (FFL), each location on the FFL is able to have a unique excitation field strength that generates a specific relaxation-induced decay signal. Through excitation of nanoparticles on the FFL with many gradient profiles, a high-resolution, one-dimensional (1D) image can be reconstructed on the FFL. For larger magnetic nanoparticles, simulation results revealed that a pulsed excitation field with a greater flat portion generates a 1D bar pattern phantom image with a higher correlation and spatial resolution. With parallel FFL and gradient coil movements, high-resolution, two-dimensional (2D) Shepp-Logan phantom and brain vessel maps were reconstructed through repetition of the spatially resolved measurement of magnetic nanoparticles on the FFL.

Recent developments of the reconstruction in magnetic particle imaging

Magnetic particle imaging (MPI) is an emerging molecular imaging technique with high sensitivity and temporal-spatial resolution. Image reconstruction is an important research topic in MPI, which converts an induced voltage signal into the image of superparamagnetic iron oxide particles concentration distribution. MPI reconstruction primarily involves system matrix- and x-space-based methods. In this review, we provide a detailed overview of the research status and future research trends of these two methods. In addition, we review the application of deep learning methods in MPI reconstruction and the current open sources of MPI. Finally, research opinions on MPI reconstruction are presented. We hope this review promotes the use of MPI in clinical applications.