Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new user3D Compton scattering imaging: study of the spectrum and contour reconstruction
120
0
0.0
(
0
)
Added by
Gael Rigaud
Publication date
2019
fields
Informatics Engineering
and research's language is
English
Authors
Gael Rigaud
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
The recent development of energy-resolving cameras opens the way to new types of applications and imaging systems. In this work, we consider computerized tomography (CT) with fan beam geometry and equipped with such cameras. The measured radiation is then a function of the positions of the source and detectors and of the energy of the incoming photons. Due to the Compton effect, the variations in energy (or spectrum) of the measurement are modeled in terms of scattering events leading to the so-called Compton scattering tomography (CST). We propose an analysis of the spectral data in terms of modelling and mapping properties which results in a general reconstruction strategy. Thanks to the supplementary information given by the energy, this joint CT-CST scanner makes accurate reconstructions of characteristics of the sought-for object possible for very few source positions and a small number of detectors. The general reconstruction strategy is finally validated on synthetic data via a total variation iterative scheme. We further show how the method can be extended to high energetic polychromatic radiation sources. Also illustrative, this work motivates the potential of combining conventional CT and Compton scattering imaging (CSI) with various architectures in 2D and 3D.
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.
Here we introduce a new reconstruction technique for two-dimensional Bragg Scattering Tomography (BST), based on the Radon transform models of [arXiv preprint, arXiv:2004.10961 (2020)]. Our method uses a combination of ideas from multibang control and microlocal analysis to construct an objective function which can regularize the BST artifacts; specifically the boundary artifacts due to sharp cutoff in sinogram space (as observed in [arXiv preprint, arXiv:2007.00208 (2020)]), and artifacts arising from approximations made in constructing the model used for inversion. We then test our algorithm in a variety of Monte Carlo (MC) simulated examples of practical interest in airport baggage screening and threat detection. The data used in our studies is generated with a novel Monte-Carlo code presented here. The model, which is available from the authors upon request, captures both the Bragg scatter effects described by BST as well as beam attenuation and Compton scatter.
We consider the two dimensional quantitative imaging problem of recovering a radiative source inside an absorbing and scattering medium from knowledge of the outgoing radiation measured at the boundary. The medium has an anisotropic scattering property that is neither negligible nor large enough for the diffusion approximation to hold. We present the numerical realization of the authors recently proposed reconstruction method. For scattering kernels of finite Fourier content in the angular variable, the solution is exact. The feasibility of the proposed algorithms is demonstrated in several numerical experiments, including simulated scenarios for parameters meaningful in optical molecular imaging.
We propose an iterative solution method for the 3D high-frequency Helmholtz equation that exploits a contour integral formulation of spectral projectors. In this framework, the solution in certain invariant subspaces is approximated by solving complex-shifted linear systems, resulting in faster GMRES iterations due to the restricted spectrum. The shifted systems are solved by exploiting a polynomial fixed-point iteration, which is a robust scheme even if the magnitude of the shift is small. Numerical tests in 3D indicate that $O(n^{1/3})$ matrix-vector products are needed to solve a high-frequency problem with a matrix size $n$ with high accuracy. The method has a small storage requirement, can be applied to both dense and sparse linear systems, and is highly parallelizable.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
