Analysis of resolution of tomographic-type reconstruction from discrete data for a class of distributions

Let $f(x)$, $x\in\mathbb R^2$, be a piecewise smooth function with a jump discontinuity across a smooth surface $\mathcal S$. Let $f{\Lambda\epsilon}$ denote the Lambda tomography (LT) reconstruction of $f$ from its discrete Radon data $\hat f(\alphak,pj)$. The sampling rate along each variable is $\sim\epsilon$. First, we compute the limit $f0(\check x)=\lim{\epsilon\to0}\epsilon f{\Lambda\epsilon}(x0+\epsilon\check x)$ for a generic $x0\in\mathcal S$. Once the limiting function $f0(\check x)$ is known (which we call the discrete transition behavior, or DTB for short), the resolution of reconstruction can be easily found. Next, we show that straight segments of $\mathcal S$ lead to non-local artifacts in $f{\Lambda\epsilon}$, and that these artifacts are of the same strength as the useful singularities of $f{\Lambda\epsilon}$. We also show that $f{\Lambda\epsilon}(x)$ does not converge to its continuous analogue $f\Lambda=(-\Delta)^{1/2}f$ as $\epsilon\to0$ even if $x\not\in\mathcal S$. Results of numerical experiments presented in the paper confirm these conclusions. We also consider a class of Fourier integral operators $\mathcal{B}$ with the same canonical relation as the classical Radon transform adjoint, and a class of distributions $g\in\mathcal{E}'(Zn)$, $Zn:=S^{{n-1}\times\mathbb} R$, and obtain easy to use formulas for the DTB when $\mathcal{B} g$ is computed from discrete data $g(\alpha{\vec k},p_j)$. Exact and LT reconstructions are particlular cases of this more general theory.