English

Analysis of reconstruction from discrete Radon transform data in $\mathbb R^3$ when the function has jump discontinuities

Numerical Analysis 2019-03-21 v1

Abstract

In this paper we study reconstruction of a function ff from its discrete Radon transform data in R3\mathbb R^3 when ff 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 ϵ\epsilon, and the density of measured directions on the unit sphere is O(ϵ2)O(\epsilon^2). Let fϵf_\epsilon denote the result of reconstruction from the discrete data. Pick any generic point x0x_0 (i.e., satisfying some mild conditions), where ff has a jump. Our first result is an explicit leading term behavior of fϵf_{\epsilon} in an O(ϵ)O(\epsilon)-neighborhood of x0x_0 as ϵ0\epsilon\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 x0∉S:=singsupp(f)x_0\not\in S:=\text{singsupp}(f) is generic, then fϵ(x0)f(x0)f_{\epsilon}(x_0)\to f(x_0) as ϵ0\epsilon\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.

Keywords

Cite

@article{arxiv.1903.08216,
  title  = {Analysis of reconstruction from discrete Radon transform data in $\mathbb R^3$ when the function has jump discontinuities},
  author = {Alexander Katsevich},
  journal= {arXiv preprint arXiv:1903.08216},
  year   = {2019}
}

Comments

3 figures