Alexander Katsevich

Rima Alaifari, Michel Defrise, Alexander Katsevich

In limited data computerized tomography, the 2D or 3D problem can be reduced to a family of 1D problems using the differentiated backprojection (DBP) method. Each 1D problem consists of recovering a compactly supported function $f \in L^2(\mathcal F)$, where $\mathcal F$ is a finite interval, from its partial Hilbert transform data. When the Hilbert transform is measured on a finite interval $\mathcal G$ that only overlaps but does not cover $\mathcal F$ this inversion problem is known to be severely ill-posed [1]. In this paper, we study the reconstruction of $f$ restricted to the overlap region $\mathcal F \cap \mathcal G$. We show that with this restriction and by assuming prior knowledge on the $L^2$ norm or on the variation of $f$, better stability with Hölder continuity (typical for mildly ill-posed problems) can be obtained.

Alexander Katsevich

In this paper we study reconstruction of a function $f$ from its discrete Radon transform data in $\mathbb R^3$ when $f$ has jump discontinuities. Consider a conventional parametrization of the Radon data in terms of the affine and angular variables. The step-size along the affine variable is $ε$, and the density of measured directions on the unit sphere is $O(ε^2)$. Let $f_ε$ denote the result of reconstruction from the discrete data. Pick any generic point $x_0$ (i.e., satisfying some mild conditions), where $f$ has a jump. Our first result is an explicit leading term behavior of $f_ε$ in an $O(ε)$-neighborhood of $x_0$ as $ε\to0$. A closely related question is why can we accurately reconstruct functions with discontinuities at all? This is a fundamental question, which has not been studied in the literature in dimensions three and higher. We prove that the discrete inversion formula `works', i.e. if $x_0\not\in S:=\text{singsupp}(f)$ is generic, then $f_ε(x_0)\to f(x_0)$ as $ε\to0$. The proof of this result reveals a surprising connection with the theory of uniform distribution (u.d.). This is a new phenomenon that has not been known previously. We also present some numerical experiments, which confirm the validity of the developed theory.

Alexander Katsevich, Dimitri Rothermel, Thomas Schuster

In this article we present an improved exact inversion formula for the 3D cone beam transform of vector fields. It is well known that only the solenoidal part of a vector field can be determined by the longitudinal ray transform of a vector field in cone beam geometry. The exact inversion formula, as it was developed in A. Katsevich and T. Schuster, 'An exact inversion formula for cone beam vector tomography', Inverse Problems 29 (2013), consists of two parts. the first part is of filtered backprojection type, whereas the second part is a costly 4D integration and very inefficient. In this article we tackle this second term and achieve an improvement which is easily to implement and saves one order of integration. The theory says that the first part contains all information about the curl of the field, whereas the second part presumably has information about the boundary values. This suggestion is supported by the fact that the second part vanishes if the exact field is divergence free and tangential at the boundary. A number of numerical tests, that are also subject of this article, confirm the theoretical results and the exactness of the formula.