## Leonid Kunyansky

**Primary Department**

**Department Affiliations**

Professor, Applied Mathematics - GIDP

Professor, BIO5 Institute

Professor, Mathematics

Work Summary

I develop mathematics of biomedical imaging. All modalities of tomography imaging rely heavily on mathematical algorithms for forming an image. I develop the theory and the algorithm enabling this technology.

Cox, B. T., Holman, B., & Kunyansky, L. (2013). Photoacoustic tomography in a reflecting cavity. Progress in Biomedical Optics and Imaging - Proceedings of SPIE, 8581.

Abstract:

Almost all known photoacoustic image reconstruction algorithms are based on the assumption that the acoustic waves leave the object (the imaged region) after a finite time. This assumption is fulfilled if the measurements are made in free space and reflections from the detectors are negligible. However, when the object is surrounded by acoustically hard detectors arrays (and/or by additional acoustic mirrors), the acoustic waves will bounce around in such a reverberant cavity many times (in the absence of absorption, forever). This paper proposes fast reconstruction algorithms for the measurements made from the walls of a rectangular reverberant cavity. The algorithms are tested using numerical simulations. © 2013 Copyright SPIE.

Bruno, O. P., & Kunyansky, L. A. (2001). Surface scattering in three dimensions: An accelerated high-order solver. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 457(2016), 2921-2934.

Abstract:

We present a new algorithm for the numerical solution of problems of acoustic scattering by surfaces in three-dimensional space. This algorithm evaluates scattered fields through fast, high-order, accurate solution of the corresponding boundary integral equation. The high-order accuracy of our solver is achieved through use of partitions of unity together with analytical resolution of kernel singularities. The acceleration, in turn, results from use of high-order equivalent source approximations, which allow for fast evaluation of non-adjacent interactions by means of the three-dimensional fast Fourier transform (FFT). Our acceleration scheme has dramatically lower memory requirements and yields much higher accuracy than existing FFT-accelerated techniques. The present algorithm computes one matrix-vector multiply in O(N^{6/5} log N) to O(N^{4/3} log N) operations (depending on the geometric characteristics of the scattering surface), it exhibits super-algebraic convergence, and it does not suffer from accuracy breakdowns of any kind. We demonstrate the efficiency of our method through a variety of examples. In particular, we show that the present algorithm can evaluate accurately, on a personal computer, scattering from bodies of acoustical sizes (ka) of several hundreds.

Kunyansky, L. A. (2004). Inversion of the 3D exponential parallel-beam transform and the Radon transform with angle-dependent attenuation. Inverse Problems, 20(5), 1455-1478.

Abstract:

The inversion problem for the 3D parallel-beam exponential ray transform is solved through inversion of a set of the 2D exponential Radon transforms with complex-valued angle-dependent attenuation. An inversion formula for the latter 2D transform is derived; it generalizes the known Kuchment-Shneiberg formula valid for real angle-dependent attenuation. We derive an explicit theoretically exact solution of the 3D problem which is valid for arbitrary closed trajectory that does not intersect itself. A simple reconstruction algorithm is described, applicable for certain sets of trajectories satisfying Orlov's condition. In the latter case, our inversion technique is as stable as the Tretiak-Metz inversion formula. Possibilities of further reduction of noise sensitivity are briefly discussed in the paper. The work of our algorithm is illustrated by an example of image reconstruction from two circular orbits.

Guillement, J. -., Jauberteau, F., Kunyansky, L., Novikov, R., & Trebossen, R. (2002). On single-photon emission computed tomography imaging based on an exact formula for the nonuniform attenuation correction. Inverse Problems, 18(6), L11-L19.

Abstract:

An exact formula approach to nonuniform attenuation correction in single-photon emission computed tomography (SPECT) was discussed. The formula admits a numerical implementation via a direct generalization of the filtered backprojection (FBP) algorithm. The formula can be used for fast computing of an efficient first approximation for more complicated SPECT reconstruction techniques.

Kunyansky, L. A. (2007). Explicit inversion formulae for the spherical mean Radon transform. Inverse Problems, 23(1), 373-383.

Abstract:

We derive explicit formulae for the reconstruction of a function from its integrals over a family of spheres, or for the inversion of the spherical mean Radon transform. Such formulae are important for problems of thermo- and photo-acoustic tomography. A closed-form inversion formula of a filtration-backprojection type is found for the case when the centres of the integration spheres lie on a sphere in surrounding the support of the unknown function. © 2007 IOP Publishing Ltd.