English

Resolution analysis of inverting the generalized Radon transform from discrete data in $\mathbb R^3$

Numerical Analysis 2019-08-29 v2 Numerical Analysis

Abstract

A number of practically important imaging problems involve inverting the generalized Radon transform (GRT) R\mathcal R of a function ff in R3\mathbb R^3. On the other hand, not much is known about the spatial resolution of the reconstruction from discretized data. In this paper we study how accurately and with what resolution the singularities of ff are reconstructed. The GRT integrates over a fairly general family of surfaces Sy\mathcal S_y in R3\mathbb R^3. Here yy is the parameter in the data space, which runs over an open set VR3\mathcal V\subset\mathbb R^3. Assume that the data g(y)=(Rf)(y)g(y)=(\mathcal R f)(y) are known on a regular grid yjy_j with step-sizes O(ϵ)O(\epsilon) along each axis, and suppose S=singsupp(f)\mathcal S=\text{singsupp}(f) is a piecewise smooth surface. Let fϵf_\epsilon denote the result of reconstruction from the descrete data. We obtain explicitly the leading singular behavior of fϵf_\epsilon in an O(ϵ)O(\epsilon)-neighborhood of a generic point x0Sx_0\in\mathcal S, where ff has a jump discontinuity. We also prove that under some generic conditions on S\mathcal S (which include, e.g. a restriction on the order of tangency of Sy\mathcal S_y and S\mathcal S), the singularities of ff do not lead to non-local artifacts. For both computations, a connection with the uniform distribution theory turns out to be important. Finally, we present a numerical experiment, which demonstrates a good match between the theoretically predicted behavior and actual reconstruction.

Keywords

Cite

@article{arxiv.1908.04753,
  title  = {Resolution analysis of inverting the generalized Radon transform from discrete data in $\mathbb R^3$},
  author = {Alexander Katsevich},
  journal= {arXiv preprint arXiv:1908.04753},
  year   = {2019}
}
R2 v1 2026-06-23T10:46:36.652Z