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