Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNumerical inversion of a broken ray transform arising in single scattering optical tomography
275
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The article presents an efficient image reconstruction algorithm for single scattering optical tomography (SSOT) in circular geometry of data acquisition. This novel medical imaging modality uses photons of light that scatter once in the body to recover its interior features. The mathematical model of SSOT is based on the broken ray (or V-line Radon) transform (BRT), which puts into correspondence to an image function its integrals along V-shaped piecewise linear trajectories. The process of image reconstruction in SSOT requires inversion of that transform. We implement numerical inversion of a broken ray transform in a disc with partial radial data. Our method is based on a relation between the Fourier coefficients of the image function and those of its BRT recently discovered by Ambartsoumian and Moon. The numerical algorithm requires solution of ill-conditioned matrix problems, which is accomplished using a half-rank truncated singular value decomposition method. Several numerical computations validating the inversion formula are presented, which demonstrate the accuracy, speed and robustness of our method in the case of both noise-free and noisy data.
rate research
Read More
The significance of the broken ray transform (BRT) is due to its occurrence in a number of modalities spanning optical, x-ray, and nuclear imaging. When data are indexed by the scatter location, the BRT is both linear and shift invariant. Analyzing the BRT as a linear system provides a new perspective on the inverse problem. In this framework we contrast prior inversion formulas and identify numerical issues. This has practical benefits as well. We clarify the extent of data required for global reconstruction by decomposing the BRT as a linear combination of cone beam transforms. Additionally we leverage the two dimensional Fourier transform to derive new inversion formulas that are computationally efficient for arbitrary scatter angles. Results of numerical simulations are presented.
We lay the foundations for a new fast method to reconstruct the electron density in x-ray scanning applications using measurements in the dark field. This approach is applied to a type of machine configuration with fixed energy sensitive (or resolving) detectors, and where the X-ray source is polychromatic. We consider the case where the measurements in the dark field are dominated by the Compton scattering process. This leads us to a 2D inverse problem where we aim to reconstruct an electron density slice from its integrals over discs whose boundaries intersect the given source point. We show that a unique solution exists for smooth densities compactly supported on an annulus centred at the source point. Using Sobolev space estimates we determine a measure for the ill posedness of our problem based on the criterion given by Natterer (The mathematics of computerized tomography SIAM 2001). In addition, with a combination of our method and the more common attenuation coefficient reconstruction, we show under certain assumptions that the atomic number of the target is uniquely determined. We test our method on simulated data sets with varying levels of added pseudo random noise.
We study a new class of Radon transforms defined on circular cones called the conical Radon transform. In $mathbb{R}^3$ it maps a function to its surface integrals over circular cones, and in $mathbb{R}^2$ it maps a function to its integrals along two rays with a common vertex. Such transforms appear in various mathematical models arising in medical imaging, nuclear industry and homeland security. This paper contains new results about inversion of conical Radon transform with fixed opening angle and vertical central axis in $mathbb{R}^2$ and $mathbb{R}^3$. New simple explicit inversion formulae are presented in these cases. Numerical simulations were performed to demonstrate the efficiency of the suggested algorithm in 2D.
We present the numerical solution of two-point boundary value problems for a third order linear PDE, representing a linear evolution in one space dimension. The difficulty of this problem is in the numerical imposition of the boundary conditions, and to our knowledge, no such computations exist. Instead of computing the evolution numerically, we evaluate the solution representation formula obtained by the unified transform of Fokas. This representation involves complex line integrals, but in order to evaluate these integrals numerically, it is necessary to deform the integration contours using appropriate deformation mappings. We formulate a strategy to implement effectively this deformation, which allows us to obtain accurate numerical results.
We propose a numerical method to approximate the scattering amplitudes for the elasticity system with a non-constant matrix potential in dimensions $d=2$ and $3$. This requires to approximate first the scattering field, for some incident waves, which can be written as the solution of a suitable Lippmann-Schwinger equation. In this work we adapt the method introduced by G. Vainikko in cite{V} to solve such equations when considering the Lame operator. Convergence is proved for sufficiently smooth potentials. Implementation details and numerical examples are also given.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
