Researcher profile

Leonid Kunyansky

Leonid Kunyansky contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2020arXiv

Microlocally accurate solution of the inverse source problem of thermoacoustic tomography

We consider the inverse source problem of thermo- and photoacoustic tomography, with data registered on an open surface partially surrounding the source of acoustic waves. Under the assumption of constant speed of sound we develop an explicit non-iterative reconstruction procedure that recovers the Radon transform of the sought source, up to an infinitely smooth additive error term. The source then can be found by inverting the Radon transform. Our analysis is microlocal in nature and does not provide a norm estimate on the error in the so obtained image. However, numerical simulations show that this error is quite small in practical terms. We also present an asymptotically fast implementation of this procedure for the case when the data are given on a circular arc in 2D.

preprint2015arXiv

A dissipative time reversal technique for photo-acoustic tomography in a cavity

We consider the inverse source problem arising in thermo- and photo-acoustic tomography. It consists in reconstructing the initial pressure from the boundary measurements of the acoustic wave. Our goal is to extend versatile time reversal techniques to the case of perfectly reflecting boundary of the domain. Standard time reversal works only if the solution of the direct problem decays in time, which does not happen in the setup we consider. We thus propose a novel time reversal technique with a non-standard boundary condition. The error induced by this time reversal technique satisfies the wave equation with a dissipative boundary condition and, therefore, decays in time. For larger measurement times, this method yields a close approximation; for smaller times, the first approximation can be iteratively refined, resulting in a convergent Neumann series for the approximation.

preprint2015arXiv

Inversion of the spherical means transform in corner-like domains by reduction to the classical Radon transform

We consider an inverse problem arising in thermo-/photo- acoustic tomography that amounts to reconstructing a function $f$ from its circular or spherical means with the centers lying on a given measurement surface. (Equivalently, these means can be expressed through the solution $u(t,x)$ of the wave equation with the initial pressure equal to $f$.) An explicit solution of this inverse problem is obtained in 3D for the surface that is the boundary of an open octet, and in 2D for the case when the centers of integration circles lie on two rays starting at the origin and intersecting at the angle equal to $π/N$, $N=2,3,4,...$. Our formulas reconstruct the Radon projections of a function closely related to $f$, from the values of $u(t,x)$ on the measurement surface. Then, function $f$ can be found by inverting the Radon transform.

preprint2013arXiv

Exterior/interior problem for the circular means transform with applications to intravascular imaging

Exterior inverse problem for the circular means transform (CMT) arises in the intravascular photoacoustic imaging (IVPA), in the intravascular ultrasound imaging (IVUS), as well as in radar and sonar. The reduction of the IPVA to the CMT is quite straightforward. As shown in the paper, in IVUS the circular means can be recovered from measurements by solving a certain Volterra integral equation. Thus, a tomography reconstruction in both modalities requires solving the exterior problem for the CMT. Numerical solution of this problem usually is not attempted due to the presence of "invisible" wavefronts, which results in severe instability of the reconstruction. The novel inversion algorithm proposed in this paper yields a stable partial reconstruction: it reproduces the "visible" part of the image and blurs the "invisible" part. If the image contains little or no invisible wavefronts (as frequently happens in the IVPA and IVUS) the reconstruction is quantitatively accurate. The presented numerical simulations demonstrate the feasibility of tomography-like reconstruction in these modalities.

preprint2011arXiv

2D and 3D reconstructions in acousto-electric tomography

We propose and test stable algorithms for the reconstruction of the internal conductivity of a biological object using acousto-electric measurements. Namely, the conventional impedance tomography scheme is supplemented by scanning the object with acoustic waves that slightly perturb the conductivity and cause the change in the electric potential measured on the boundary of the object. These perturbations of the potential are then used as the data for the reconstruction of the conductivity. The present method does not rely on "perfectly focused" acoustic beams. Instead, more realistic propagating spherical fronts are utilized, and then the measurements that would correspond to perfect focusing are synthesized. In other words, we use \emph{synthetic focusing}. Numerical experiments with simulated data show that our techniques produce high quality images, both in 2D and 3D, and that they remain accurate in the presence of high-level noise in the data. Local uniqueness and stability for the problem also hold.

preprint2011arXiv

Fast reconstruction algorithms for the thermoacoustic tomography in certain domains with cylindrical or spherical symmetries

We propose three fast algorithms for solving the inverse problem of the thermoacoustic tomography corresponding to certain acquisition geometries. Two of these methods are designed to process the measurements done with point-like detectors placed on a circle (in 2D) or a sphere (in 3D) surrounding the object of interest. The third inversion algorithm works with the data measured by the integrating line detectors arranged in a cylindrical assembly rotating around the object. The number of operations required by these techniques is equal to O(n^3 log n) and O(n^3 log^2 n) for the 3D techniques (assuming the reconstruction grid with n^3 nodes) and to O(n^2 log n) for the 2D problem with n-by-n discretizetion grid. Numerical simulations show that our methods are at least two orders of magnitude faster than the existing algorithms, without any sacrifice in accuracy or stability. The results of reconstructions from real measurements done by the integrating line detectors are also presented, to demonstrate the practicality of our methods.

preprint2010arXiv

Reconstruction of a function from its spherical (circular) means with the centers lying on the surface of certain polygons and polyhedra

We present explicit filtration/backprojection-type formulae for the inversion of the spherical (circular) mean transform with the centers lying on the boundary of some polyhedra (or polygons, in 2D). The formulae are derived using the double layer potentials for the wave equation, for the domains with certain symmetries. The formulae are valid for a rectangle and certain triangles in 2D, and for a cuboid, certain right prisms and a certain pyramid in 3D. All the present inversion formulae yield exact reconstruction within the domain surrounded by the acquisition surface even in the presence of exterior sources.

preprint2009arXiv

Mathematics of Photoacoustic and Thermoacoustic Tomography

This is the manuscript of the chapter for a planned Handbook of Mathematical Methods in Imaging that surveys the mathematical models, problems, and algorithms of the Thermoacoustic (TAT) and Photoacoustic (PAT) Tomography. TAT and PAT represent probably the most developed of the several novel ``hybrid'' methods of medical imaging. These new modalities combine different physical types of waves (electromagnetic and acoustic in case of TAT and PAT) in such a way that the resolution and contrast of the resulting method are much higher than those achievable using only acoustic or electromagnetic measurements.

preprint2009arXiv

Synthetic focusing in ultrasound modulated tomography

Several hybrid tomographic methods utilizing ultrasound modulation have been introduced lately. Success of these methods hinges on the feasibility of focusing ultrasound waves at an arbitrary point of interest. Such a focusing, however, is difficult to achieve in practice. We thus propose a way to avoid the use of focused waves through the so called synthetic focusing, i.e. by the reconstruction of the would-be response to the focused modulation from the measurements corresponding to realistic unfocused waves. Examples of reconstructions from simulated data are provided. This non-technical paper describes only the general concept, while technical details will appear elsewhere.