Source author record

Alexander Katsevich

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

math.CA math.NA math.FA math.SP Numerical Analysis math-ph math.MP physics.data-an Quantitative Methods

Catalog footprint

Research overview

12works

9topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2022arXiv

Novel resolution analysis for the Radon transform in $\mathbb R^2$ for functions with rough edges

Let $f$ be a function in $\mathbb R^2$ , which has a jump across a smooth curve $\mathcal S$ with nonzero curvature. We consider a family of functions $f_ε$ with jumps across a family of curves $\mathcal S_ε$ . Each $\mathcal S_ε$ is an $O(ε)$ -size perturbation of $\mathcal S$ , which scales like $O(ε^{-1/2})$ along $\mathcal S$ . Let $f_ε^{\text{rec}}$ be the reconstruction of $f_ε$ from its discrete Radon transform data, where $ε$ is the data sampling rate. A simple asymptotic (as $ε\to0$ ) formula to approximate $f_ε^{\text{rec}}$ in any $O(ε)$ -size neighborhood of $\mathcal S$ was derived heuristically in an earlier paper of the author. Numerical experiments revealed that the formula is highly accurate even for nonsmooth (i.e., only H{ö}lder continuous) $\mathcal S_ε$ . In this paper we provide a full proof of this result, which says that the magnitude of the error between $f_ε^{\text{rec}}$ and its approximation is $O(ε^{1/2}\ln(1/ε))$ . The main assumption is that the level sets of the function $H_0(\cdot,ε)$ , which parametrizes the perturbation $\mathcal S\to\mathcal S_ε$ , are not too dense.

preprint2021arXiv

Resolution analysis of inverting the generalized $N$ -dimensional Radon transform in $\mathbb R^n$ from discrete data

Let $\mathcal R$ denote the generalized Radon transform (GRT), which integrates over a family of $N$ -dimensional smooth submanifolds $\mathcal S_{\tilde y}\subset\mathcal U$ , $1\le N\le n-1$ , where an open set $\mathcal U\subset\mathbb R^n$ is the image domain. The submanifolds are parametrized by points $\tilde y\subset\tilde{\mathcal V}$ , where an open set $\tilde{\mathcal V}\subset\mathbb R^n$ is the data domain. The continuous data are $g={\mathcal R} f$ , and the reconstruction is $\check f=\mathcal R^*\mathcal B g$ . Here $\mathcal R^*$ is a weighted adjoint of $\mathcal R$ , and $\mathcal B$ is a pseudo-differential operator. We assume that $f$ is a conormal distribution, $\text{supp}(f)\subset\mathcal U$ , and its singular support is a smooth hypersurface $\mathcal S\subset\mathcal U$ . Discrete data consists of the values of $g$ on a lattice $\tilde y^j$ with the step size $O(ε)$ . Let $\check f_ε=\mathcal R^*\mathcal B g_ε$ denote the reconstruction obtained by applying the inversion formula to an interpolated discrete data $g_ε(\tilde y)$ . Pick a generic pair $(x_0,\tilde y_0)$ , where $x_0\in\mathcal S$ , and $\mathcal S_{\tilde y_0}$ is tangent to $\mathcal S$ at $x_0$ . The main result of the paper is the computation of the limit $f_0(\check x):=\lim_{ε\to0}ε^κ\check f_ε(x_0+ε\check x).$ Here $κ\ge 0$ is selected based on the strength of the reconstructed singularity, and $\check x$ is confined to a bounded set. The limiting function $f_0(\check x)$ , which we call the discrete transition behavior, allows computing the resolution of reconstruction.

preprint2020arXiv

On the spectral properties of the Hilbert transform operator on multi-intervals

Let $J,E\subset\mathbb R$ be two multi-intervals with non-intersecting interiors. Consider the following operator $A:\, L^2( J )\to L^2(E),\ (Af)(x) = \frac 1π\int_{ J } \frac {f(y)\text{d} y}{x-y},$ and let $A^\dagger$ be its adjoint. We introduce a self-adjoint operator $\mathscr K$ acting on $L^2(E)\oplus L^2(J)$ , whose off-diagonal blocks consist of $A$ and $A^\dagger$ . In this paper we study the spectral properties of $\mathscr K$ and the operators $A^\dagger A$ and $A A^\dagger$ . Our main tool is to obtain the resolvent of $\mathscr K$ , which is denoted by $\mathscr R$ , using an appropriate Riemann-Hilbert problem, and then compute the jump and poles of $\mathscr R$ in the spectral parameter $λ$ . We show that the spectrum of $\mathscr K$ has an absolutely continuous component $[0,1]$ if and only if $J$ and $E$ have common endpoints, and its multiplicity equals to their number. If there are no common endpoints, the spectrum of $\mathscr K$ consists only of eigenvalues and $0$ . If there are common endpoints, then $\mathscr K$ may have eigenvalues imbedded in the continuous spectrum, each of them has a finite multiplicity, and the eigenvalues may accumulate only at $0$ . In all cases, $\mathscr K$ does not have a singular continuous spectrum. The spectral properties of $A^\dagger A$ and $A A^\dagger$ , which are very similar to those of $\mathscr K$ , are obtained as well.

preprint2016arXiv

An improved exact inversion formula for cone beam vector tomography

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.

preprint2016arXiv

Reconstruction algorithms for a class of restricted ray transforms without added singularities

Let $X$ and $X^*$ denote a restricted ray transform along curves and a corresponding backprojection operator, respectively. Theoretical analysis of reconstruction from the data $Xf$ is usually based on a study of the composition $X^* D X$ , where $D$ is some local operator (usually a derivative). If $X^*$ is chosen appropriately, then $X^* D X$ is a Fourier Integral Operator (FIO) with singular symbol. The singularity of the symbol leads to the appearance of artifacts (added singularities) that can be as strong as the original (or, useful) singularities. By choosing $D$ in a special way one can reduce the strength of added singularities, but it is impossible to get rid of them completely. In the paper we follow a similar approach, but make two changes. First, we replace $D$ with a nonlocal operator $\tilde D$ that integrates $Xf$ along a curve in the data space. The result $\tilde D Xf$ resembles the generalized Radon transform $R$ of $f$ . The function $\tilde D Xf$ is defined on pairs $(x_0,Θ)\in U\times S^2$ , where $U\subset\mathbb R^3$ is an open set containing the support of $f$ , and $S^2$ is the unit sphere in $\mathbb R^3$ . Second, we replace $X^*$ with a backprojection operator $R^*$ that integrates with respect to $Θ$ over $S^2$ . It turns out that if $\tilde D$ and $R^*$ are appropriately selected, then the composition $R^* \tilde D X$ is an elliptic pseudodifferential operator of order zero with principal symbol 1. Thus, we obtain an approximate reconstruction formula that recovers all the singularities correctly and does not produce artifacts. The advantage of our approach is that by inserting $\tilde D$ we get access to the frequency variable $Θ$ . In particular, we can incorporate suitable cut-offs in $R^*$ to eliminate bad directions $Θ$ , which lead to added singularities.

preprint2016arXiv

Stability estimates for the regularized inversion of the truncated Hilbert transform

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.

preprint2015arXiv

Diagonalization of the Finite Hilbert Transform on two adjacent intervals

We study the interior problem of tomography. The starting point is the Gelfand-Graev formula, which converts the tomographic data into the finite Hilbert transform (FHT) of an unknown function $f$ along a collection of lines. Pick one such line, call it the $x$ -axis, and assume that the function to be reconstructed depends on a one-dimensional argument by restricting $f$ to the line. Let $Ω_1$ be the interval where $f$ is supported, and $Ω_2$ be the interval where the Hilbert transform of $f$ can be computed using the Gelfand-Graev formula. The equation we study is $H_1 f=g|_{Ω_2}$ , where $H_1$ is the FHT that integrates over $Ω_1$ and gives the result on $Ω_2$ , i.e. $H_1: L^2(Ω_1)\to L^2(Ω_2)$ . In the case of the interior problem the tomographic data are truncated, and $Ω_1$ is no longer a subset of $Ω_2$ . In this paper we consider the case when the intervals $Ω_1=(a_1,0)$ and $Ω_2=(0,a_2)$ are adjacent. Here $a_1 < 0 < a_2$ . First we find a differential operator $L$ that commutes with $H_1$ . Using the Titchmarsh-Weyl theory, we show that $L$ has only continuous spectrum and obtain two isometric transformations $U_1$ , $U_2$ , such that $U_2 H_1 U_1^*$ is the multiplication operator with the function $σ(λ)$ , $λ\geq(a_1^2+a_2^2)/8$ . Here $λ$ is the spectral parameter. Then we show that $σ(λ)\to0$ as $λ\to\infty$ exponentially fast. We also obtain the leading asymptotic behavior of the kernels involved in the integral operators $U_1$ , $U_2$ as $λ\to\infty$ . When the intervals are symmetric, i.e. $-a_1=a_2$ , the operators $U_1$ , $U_2$ are obtained explicitly in terms of hypergeometric functions.

preprint2014arXiv

Covariance Matrix Estimation for the Cryo-EM Heterogeneity Problem

In cryo-electron microscopy (cryo-EM), a microscope generates a top view of a sample of randomly-oriented copies of a molecule. The problem of single particle reconstruction (SPR) from cryo-EM is to use the resulting set of noisy 2D projection images taken at unknown directions to reconstruct the 3D structure of the molecule. In some situations, the molecule under examination exhibits structural variability, which poses a fundamental challenge in SPR. The heterogeneity problem is the task of mapping the space of conformational states of a molecule. It has been previously suggested that the leading eigenvectors of the covariance matrix of the 3D molecules can be used to solve the heterogeneity problem. Estimating the covariance matrix is challenging, since only projections of the molecules are observed, but not the molecules themselves. In this paper, we formulate a general problem of covariance estimation from noisy projections of samples. This problem has intimate connections with matrix completion problems and high-dimensional principal component analysis. We propose an estimator and prove its consistency. When there are finitely many heterogeneity classes, the spectrum of the estimated covariance matrix reveals the number of classes. The estimator can be found as the solution to a certain linear system. In the cryo-EM case, the linear operator to be inverted, which we term the projection covariance transform, is an important object in covariance estimation for tomographic problems involving structural variation. Inverting it involves applying a filter akin to the ramp filter in tomography. We design a basis in which this linear operator is sparse and thus can be tractably inverted despite its large size. We demonstrate via numerical experiments on synthetic datasets the robustness of our algorithm to high levels of noise.

preprint2014arXiv

On Sobolev instability of the interior problem of tomography

In this paper we continue investigation of the interior problem of tomography that was started in \cite{BKT2}. As is known, solving the interior problem {with prior data specified on a finite collection of intervals $I_i$ } is equivalent to analytic continuation of a function from $I_i$ to an open set ${\bf J}$ . In the paper we prove that this analytic continuation can be obtained with the help of a simple explicit formula, which involves summation of a series. Our second result is that the operator of analytic continuation is not stable for any pair of Sobolev spaces regardless of how close the set ${\bf J}$ is to $I_i$ . Our main tool is the singular value decomposition of the operator $\mathcal H^{-1}_e$ that arises when the interior problem is reduced to a problem of inverting the Hilbert transform from incomplete data. The asymptotics of the singular values and singular functions of $\mathcal H^{-1}_e$ , the latter being valid uniformly on compact subsets {of the interior of $I_i$ }, was obtained in \cite{BKT2}. {Using these asymptotics we can accurately measure the degree of ill-posedness of the analytic continuation as a function of the target interval ${\bf J}$ .} Our {last} result is the convergence of the asymptotic approximation of the singular functions {in the $L^2(I_i)$ sense}.

preprint2014arXiv

Singular value decomposition of a finite Hilbert transform defined on several intervals and the interior problem of tomography: the Riemann-Hilbert problem approach

We study the asymptotics of singular values and singular functions of a Finite Hilbert transform (FHT), which is defined on several intervals. Transforms of this kind arise in the study of the interior problem of tomography. We suggest a novel approach based on the technique of the matrix Riemann-Hilbert problem and the steepest descent method of Deift-Zhou. We obtain a family of matrix RHPs depending on the spectral parameter $λ$ and show that the singular values of the FHT coincide with the values of $λ$ for which the RHP is not solvable. Expressing the leading order solution as $λ\to 0$ of the RHP in terms of the Riemann Theta functions, we prove that the asymptotics of the singular values can be obtained by studying the intersections of the locus of zeroes of a certain Theta function with a straight line. This line can be calculated explicitly, and it depends on the geometry of the intervals that define the FHT. The leading order asymptotics of the singular functions and singular values are explicitly expressed in terms of the Riemann Theta functions and of the period matrix of the corresponding normalized differentials, respectively. We also obtain the error estimates for our asymptotic results.

preprint2013arXiv

Asymptotic analysis of the SVD for the truncated Hilbert transform with overlap

The truncated Hilbert transform with overlap $H_T$ is an operator that arises in tomographic reconstruction from limited data, more precisely in the method of Differentiated Back-Projection (DBP). Recent work [1] has shown that the singular values of this operator accumulate at both zero and one. To better understand the properties of the operator and, in particular, the ill-posedness of the inverse problem associated with it, it is of interest to know the rates at which the singular values approach zero and one. In this paper, we exploit the property that $H_T$ commutes with a second-order differential operator $L_S$ and the global asymptotic behavior of its eigenfunctions to find the asymptotics of the singular values and singular functions of $H_T$ .

preprint2013arXiv

Spectral analysis of the truncated Hilbert transform with overlap

We study a restriction of the Hilbert transform as an operator $H_T$ from $L^2(a_2,a_4)$ to $L^2(a_1,a_3)$ for real numbers $a_1 < a_2 < a_3 < a_4$ . The operator $H_T$ arises in tomographic reconstruction from limited data, more precisely in the method of differentiated back-projection (DBP). There, the reconstruction requires recovering a family of one-dimensional functions $f$ supported on compact intervals $[a_2,a_4]$ from its Hilbert transform measured on intervals $[a_1,a_3]$ that might only overlap, but not cover $[a_2,a_4]$ . We show that the inversion of $H_T$ is ill-posed, which is why we investigate the spectral properties of $H_T$ . We relate the operator $H_T$ to a self-adjoint two-interval Sturm-Liouville problem, for which we prove that the spectrum is discrete. The Sturm-Liouville operator is found to commute with $H_T$ , which then implies that the spectrum of $H_T^* H_T$ is discrete. Furthermore, we express the singular value decomposition of $H_T$ in terms of the solutions to the Sturm-Liouville problem. The singular values of $H_T$ accumulate at both $0$ and $1$ , implying that $H_T$ is not a compact operator. We conclude by illustrating the properties obtained for $H_T$ numerically.

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Alexander Katsevich

Research overview

Is this your work?

12 publications

Novel resolution analysis for the Radon transform in $\mathbb R^2$ for functions with rough edges

Resolution analysis of inverting the generalized $N$ -dimensional Radon transform in $\mathbb R^n$ from discrete data

On the spectral properties of the Hilbert transform operator on multi-intervals

An improved exact inversion formula for cone beam vector tomography

Reconstruction algorithms for a class of restricted ray transforms without added singularities

Stability estimates for the regularized inversion of the truncated Hilbert transform

Diagonalization of the Finite Hilbert Transform on two adjacent intervals

Covariance Matrix Estimation for the Cryo-EM Heterogeneity Problem

On Sobolev instability of the interior problem of tomography

Singular value decomposition of a finite Hilbert transform defined on several intervals and the interior problem of tomography: the Riemann-Hilbert problem approach

Asymptotic analysis of the SVD for the truncated Hilbert transform with overlap

Spectral analysis of the truncated Hilbert transform with overlap

See the researcher in context

Building this map preview

Alexander Katsevich

Research overview

Is this your work?

12 publications

Novel resolution analysis for the Radon transform in R2\mathbb R^2R2 for functions with rough edges

Resolution analysis of inverting the generalized NNN-dimensional Radon transform in Rn\mathbb R^nRn from discrete data

On the spectral properties of the Hilbert transform operator on multi-intervals

An improved exact inversion formula for cone beam vector tomography

Reconstruction algorithms for a class of restricted ray transforms without added singularities

Stability estimates for the regularized inversion of the truncated Hilbert transform

Diagonalization of the Finite Hilbert Transform on two adjacent intervals

Covariance Matrix Estimation for the Cryo-EM Heterogeneity Problem

On Sobolev instability of the interior problem of tomography

Singular value decomposition of a finite Hilbert transform defined on several intervals and the interior problem of tomography: the Riemann-Hilbert problem approach

Asymptotic analysis of the SVD for the truncated Hilbert transform with overlap

Spectral analysis of the truncated Hilbert transform with overlap

See the researcher in context

Building this map preview

Novel resolution analysis for the Radon transform in $\mathbb R^2$ for functions with rough edges

Resolution analysis of inverting the generalized $N$ -dimensional Radon transform in $\mathbb R^n$ from discrete data