Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userMathematical Analysis of the 1D Model and Reconstruction Schemes for Magnetic Particle Imaging
110
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Magnetic particle imaging (MPI) is a promising new in-vivo medical imaging modality in which distributions of super-paramagnetic nanoparticles are tracked based on their response in an applied magnetic field. In this paper we provide a mathematical analysis of the modeled MPI operator in the univariate situation. We provide a Hilbert space setup, in which the MPI operator is decomposed into simple building blocks and in which these building blocks are analyzed with respect to their mathematical properties. In turn, we obtain an analysis of the MPI forward operator and, in particular, of its ill-posedness properties. We further get that the singular values of the MPI core operator decrease exponentially. We complement our analytic results by some numerical studies which, in particular, suggest a rapid decay of the singular values of the MPI operator.
rate research
Read More
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.
We propose a new variational model for joint image reconstruction and motion estimation in spatiotemporal imaging, which is investigated along a general framework that we present with shape theory. This model consists of two components, one for conducting modified static image reconstruction, and the other performs sequentially indirect image registration. For the latter, we generalize the large deformation diffeomorphic metric mapping framework into the sequentially indirect registration setting. The proposed model is compared theoretically against alternative approaches (optical flow based model and diffeomorphic motion models), and we demonstrate that the proposed model has desirable properties in terms of the optimal solution. The theoretical derivations and efficient algorithms are also presented for a time-discretized scenario of the proposed model, which show that the optimal solution of the time-discretized version is consistent with that of the time-continuous one, and most of the computational components is the easy-implemented linearized deformation. The complexity of the algorithm is analyzed as well. This work is concluded by some numerical examples in 2D space + time tomography with very sparse and/or highly noisy data.
In this paper, we consider the imaging problem of terahertz (THz) tomography, in particular as it appears in non-destructive testing. We derive a nonlinear mathematical model describing a full THz tomography experiment, and consider linear approximations connecting THz tomography with standard computerized tomography and the Radon transform. Based on the derived models we propose different reconstruction approaches for solving the THz tomography problem, which we then compare on experimental data obtained from THz measurements of a plastic sample.
We consider in this work the numerical computation of transport coefficients for Brownian dynamics. We investigate the discretization error arising when simulating the dynamics with the Smart MC algorithm (also known as Metropolis-adjusted Langevin algorithm). We prove that the error is of order one in the time step, when using either the Green-Kubo or the Einstein formula to estimate the transport coefficients. We illustrate our results with numerical simulations.
3D Compton scattering imaging is an upcoming concept exploiting the scattering of photons induced by the electronic structure of the object under study. The so-called Compton scattering rules the collision of particles with electrons and describes their energy loss after scattering. Although physically relevant, multiple-order scattering was so far not considered and therefore, only first-order scattering is generally assumed in the literature. The purpose of this work is to argument why and how a contour reconstruction of the electron density map from scattered measurement composed of first- and second-order scattering is possible (scattering of higher orders is here neglected). After the development of integral representations for the first- and second-order scattering, this is achieved by the study of the smoothness properties of associated Fourier integral operators (FIO). The second-order scattered radiation reveals itself to be structurally smoother than the radiation of first-order indicating that the contours of the electron density are essentially encoded within the first-order part. This opens the way to contour-based reconstruction techniques when using multiple scattered data. Our main results, modeling and reconstruction scheme, are successfully implemented on synthetic and Monte-Carlo data.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
