English

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

Numerical Analysis 2022-06-10 v1 Numerical Analysis

Abstract

Let ff be a function in R2\mathbb R^2, which has a jump across a smooth curve S\mathcal S with nonzero curvature. We consider a family of functions fϵf_\epsilon with jumps across a family of curves Sϵ\mathcal S_\epsilon. Each Sϵ\mathcal S_\epsilon is an O(ϵ)O(\epsilon)-size perturbation of S\mathcal S, which scales like O(ϵ1/2)O(\epsilon^{-1/2}) along S\mathcal S. Let fϵrecf_\epsilon^{\text{rec}} be the reconstruction of fϵf_\epsilon from its discrete Radon transform data, where ϵ\epsilon is the data sampling rate. A simple asymptotic (as ϵ0\epsilon\to0) formula to approximate fϵrecf_\epsilon^{\text{rec}} in any O(ϵ)O(\epsilon)-size neighborhood of S\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{\"o}lder continuous) Sϵ\mathcal S_\epsilon. In this paper we provide a full proof of this result, which says that the magnitude of the error between fϵrecf_\epsilon^{\text{rec}} and its approximation is O(ϵ1/2ln(1/ϵ))O(\epsilon^{1/2}\ln(1/\epsilon)). The main assumption is that the level sets of the function H0(,ϵ)H_0(\cdot,\epsilon), which parametrizes the perturbation SSϵ\mathcal S\to\mathcal S_\epsilon, are not too dense.

Keywords

Cite

@article{arxiv.2206.04545,
  title  = {Novel resolution analysis for the Radon transform in $\mathbb R^2$ for functions with rough edges},
  author = {Alexander Katsevich},
  journal= {arXiv preprint arXiv:2206.04545},
  year   = {2022}
}
R2 v1 2026-06-24T11:45:15.412Z