Gaik Ambartsoumian

Gaik Ambartsoumian, Souvik Roy

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.

Gaik Ambartsoumian, Rim Gouia-Zarrad, Venkateswaran P. Krishnan et al.

We study inversion of the spherical Radon transform with centers on a sphere (the data acquisition set). Such inversions are essential in various image reconstruction problems arising in medical, radar and sonar imaging. In the case of radially incomplete data, we show that the spherical Radon transform can be uniquely inverted recovering the image function in spherical shells. Our result is valid when the support of the image function is inside the data acquisition sphere, outside that sphere, as well as on both sides of the sphere. Furthermore, in addition to the uniqueness result our method of proof provides reconstruction formulas for all those cases. We present a robust computational algorithm based on our inversion formula and demonstrate its accuracy and efficiency on several numerical examples.

Gaik Ambartsoumian, Sarah K. Patch

The problem of image reconstruction in thermoacoustic tomography requires inversion of a generalized Radon transform, which integrates the unknown function over circles in 2D or spheres in 3D. The paper investigates implementation of the recently discovered backprojection type inversion formulas for the case of spherical acquisition in 3D. A numerical simulation of the data acquisition with subsequent reconstructions are made for the Defrise phantom as well as for some other phantoms. Both full and partial scan situations are considered. The results are compared with the implementation of a previously used approximate inversion formula.