Spectral filtering for the resolution of the Gibbs phenomenon in MPI applications

open3Polynomial interpolation on the node points of Lissajous curves using Chebyshev series is an e effective way for a fast image reconstruction in Magnetic Particle Imaging. Due to the nature of spectral methods, a Gibbs phenomenon occurs in the reconstructed image if the underlying function has discontinuities. A possible solution for this problem are spectral filtering methods acting on the coefficients of the interpolating polynomial. In this work, after a description of the Gibbs phenomenon in two dimensions, we present an adaptive spectral filtering process for the resolution of this phenomenon and for an improved approximation of the underlying function or image. In this adaptive filtering technique, the spectral filter depends on the distance of a spatial point to the nearest discontinuity. We show the effectiveness of this filtering approach in theory, in numerical simulations as well as in the application in Magnetic Particle Imaging.openDe Marchi, Stefano; Erb, Wolfgang; Marchetti, Francesco.DE MARCHI, Stefano; Erb, Wolfgang; Marchetti, Francesc

Spectral filtering for the reduction of the Gibbs phenomenon of polynomial approximation methods on Lissajous curves with applications in MPI

Polynomial interpolation and approximation methods on sampling points along Lissajous curves using Chebyshev series is an effective way for a fast image reconstruction in Magnetic Particle Imaging. Due to the nature of spectral methods, a Gibbs phenomenon occurs in the reconstructed image if the underlying function has discontinuities. A possible solution for this problem are spectral filtering methods acting on the coefficients of the approximating polynomial. In this work, after a description of the Gibbs phenomenon and classical filtering techniques in one and several dimensions, we present an adaptive spectral filtering process for the resolution of this phenomenon and for an improved approximation of the underlying function or image. In this adaptive filtering technique, the spectral filter depends on the distance of a spatial point to the nearest discontinuity. We show the effectiveness of this filtering approach in theory, in numerical simulations as well as in the application in Magnetic Particle Imaging

Shape-Driven Interpolation With Discontinuous Kernels: Error Analysis, Edge Extraction, and Applications in Magnetic Particle Imaging

Accurate interpolation and approximation techniques for functions with discontinuities are key tools in many applications, such as medical imaging. In this paper, we study a radial basis function type of method for scattered data interpolation that incorporates discontinuities via a variable scaling function. For the construction of the discontinuous basis of kernel functions, information on the edges of the interpolated function is necessary. We characterize the native space spanned by these kernel functions and study error bounds in terms of the fill distance of the node set. To extract the location of the discontinuities, we use a segmentation method based on a classification algorithm from machine learning. The results of the conducted numerical experiments are in line with the theoretically derived convergence rates in case that the discontinuities are a priori known. Further, an application to interpolation in magnetic particle imaging shows that the presented method is very promising in order to obtain edge-preserving image reconstructions in which ringing artifacts are reduced

A new 3D model for magnetic particle imaging using realistic magnetic field topologies for algebraic reconstruction

We derive a new 3D model for magnetic particle imaging (MPI) that is able to incorporate realistic magnetic fields in the reconstruction process. In real MPI scanners, the generated magnetic fields have distortions that lead to deformed magnetic low-field volumes (LFV) with the shapes of ellipsoids or bananas instead of ideal field-free points (FFP) or lines (FFL), respectively. Most of the common model-based reconstruction schemes in MPI use however the idealized assumption of an ideal FFP or FFL topology and, thus, generate artifacts in the reconstruction. Our model-based approach is able to deal with these distortions and can generally be applied to dynamic magnetic fields that are approximately parallel to their velocity field. We show how this new 3D model can be discretized and inverted algebraically in order to recover the magnetic particle concentration. To model and describe the magnetic fields, we use decompositions of the fields in spherical harmonics. We complement the description of the new model with several simulations and experiments.Comment: 27 pages, 11 figure, 3 table

System Characterization of a Human-Sized 3D Real-Time Magnetic Particle Imaging Scanner for Cerebral Applications

Since the initial patent in 2001, the Magnetic Particle Imaging (MPI) community has been striving to develop an MPI scanner suitable for human applications. Numerous contributions from different research fields, regarding tracer development, reconstruction methods, hardware engineering, and sequence design have been employed in pursuit of this objective. In this work, we introduce and thoroughly characterize an improved head-sized MPI scanner with an emphasis on human safety. The scanner is operated by open-source software that enables scanning, monitoring, analysis, and reconstruction, designed to be handled by end users. Our primary focus is to present all technical components of the scanner, with the ultimate objective to investigate brain perfusion imaging in phantom experiments. We have successfully achieved full 3D single- and multi-contrast imaging capabilities at a frame rate of 4 Hz with sufficient sensitivity and resolution for brain applications. To assess system characterization, we devised sensitivity, resolution, perfusion, and multi-contrast experiments, as well as field measurements and sequence analysis. The acquired images were captured using a clinically approved tracer and suitable magnetic field strengths, while adhering to the established human peripheral nerve stimulation thresholds. This advanced scanner holds potential as a tomographic imager for diagnosing conditions such as ischemic stroke or intracranial hemorrhage in environments lacking electromagnetic shielding. Furthermore, due to its low power consumption it may have the potential to facilitate long-term monitoring within intensive care units for various applications.Comment: 22 pages, 9 figure

Spectral filtering for the resolution of the Gibbs phenomenon in MPI applications by Lissajous sampling

Polynomial interpolation and approximation methods on sampling points along Lissajous curves using Chebyshev series is an e-ffective way for a fast image reconstruc-tion in Magnetic Particle Imaging. Due to the nature of spectral methods, a Gibbs phenomenon occurs in the reconstructed image if the underlying function has discon-tinuities. A possible solution for this problem are spectral filtering methods acting on the coefficients of the approximating polynomial. In this work, after a description of the Gibbs phenomenon and classical filtering techniques in one and several dimensions, we present an adaptive spectral filtering process for the resolution of this phenomenon and for an improved approximation of the underlying function or image. In this adaptive filtering technique, the spectral filter depends on the distance of a spatial point to the nearest discontinuity. We show the e-ffectiveness of this filtering approach in theory, in numerical simulations as well as in the application in Magnetic Particle Imaging.ope

TranSMS: Transformers for Super-Resolution Calibration in Magnetic Particle Imaging

Magnetic particle imaging (MPI) offers exceptional contrast for magnetic nanoparticles (MNP) at high spatio-temporal resolution. A common procedure in MPI starts with a calibration scan to measure the system matrix (SM), which is then used to set up an inverse problem to reconstruct images of the MNP distribution during subsequent scans. This calibration enables the reconstruction to sensitively account for various system imperfections. Yet time-consuming SM measurements have to be repeated under notable changes in system properties. Here, we introduce a novel deep learning approach for accelerated MPI calibration based on Transformers for SM super-resolution (TranSMS). Low-resolution SM measurements are performed using large MNP samples for improved signal-to-noise ratio efficiency, and the high-resolution SM is super-resolved via model-based deep learning. TranSMS leverages a vision transformer module to capture contextual relationships in low-resolution input images, a dense convolutional module for localizing high-resolution image features, and a data-consistency module to ensure measurement fidelity. Demonstrations on simulated and experimental data indicate that TranSMS significantly improves SM recovery and MPI reconstruction for up to 64-fold acceleration in two-dimensional imaging

Auto-diagnosis of time-of-flight for ultrasonic signal based on defect peaks tracking model

With the popularization of humans working in tandem with robots and artificial intelligence (AI) by Industry 5.0, ultrasonic non-destructive testing (NDT)) technology has been increasingly used in quality inspections in the industry. As a crucial part of handling ultrasonic testing results–signal processing, the current approach focuses on professional training to perform signal discrimination but automatic and intelligent signal optimization and estimation lack systematic research. Though the automated and intelligent framework for ultrasonic echo signal processing has already exhibited essential research significance for diagnosing defect locations, the real-time applicability of the algorithm for the time-of-flight (ToF) estimation is rarely considered, which is a very important indicator for intelligent detection. This paper conducts a systematic comparison among different ToF algorithms for the first time and presents the auto-diagnosis of the ToF approach based on the Defect Peaks Tracking Model (DPTM). The proposed DPTM is used for ultrasonic echo signal processing and recognition for the first time. The DPTM using the Hilbert transform was verified to locate the defect with the size of 2–10 mm, in which the wavelet denoising method was adopted. With the designed mechanical fixture through 3D printing technology on the pipeline to inspect defects, the difficulty of collecting sufficient data could be conquered. The maximum auto-diagnosis error could be reduced to 0.25% and 1.25% for steel plate and pipeline under constant pressure, respectively, which were much smaller than those with the DPTM adopting the cross-correlation. The real-time auto-diagnosis identification feature of DPTM has the potential to be combined with AI in future work, such as machine learning and deep learning, to achieve more intelligent approaches for industrial health inspection

Non-Equispaced System Matrix Acquisition for Magnetic Particle Imaging Based on Lissajous Node Points

Magnetic Particle Imaging (MPI) is an emerging technology in the field of (pre)clinical imaging. The acquisition of a particle signal is realized along specific sampling trajectories covering a defined field of view (FOV). In a system matrix (SM) based reconstruction procedure, the commonly used acquisition path in MPI is a Lissajous trajectory. Such a trajectory features an inhomogeneous coverage of the FOV, i.e. the center region is sampled less dense than the regions towards the edges of the FOV. Conventionally, the respective SM acquisition and the subsequent reconstruction do not reflect this inhomogeneous coverage. Instead, they are performed on an equispaced grid. The objective of this work is to introduce a sampling grid that inherently features the aforementioned inhomogeneity by using node points of Lissajous trajectories. Paired with a tailored polynomial interpolation of the reconstructed MPI signal, the entire image can be recovered. It is the first time that such a trajectory related non-equispaced grid is used for image reconstruction on simulated and measured MPI data and it is shown that the number of sampling positions can be reduced, while the spatial resolution remains constant