RESOLUTION OF 2 DIMENSIONAL RECONSTRUCTION OF FUNCTIONS WITH NONSMOOTH EDGES FROM DISCRETE RADON TRANSFORM DATA

Let f\epsilon, 0 < euro < euro0, be a family of functions in R2, and f rec \epsilon be a reconstruction of f\epsilon from its discrete Radon transform data. Here euro is both the data sampling rate and the parameter of the family. We study the resolution of reconstruction when f\epsilon has a jump discontinuity along a nonsmooth curve S\epsilon. The assumptions are that (a) S\epsilon is a family of O(euro)-size perturbations of a smooth curve S, and (b) S\epsilon is Ho"\lder continuous with some exponent \gamma E (0, 1]. Thus the size of the perturbation S-(sic).S\epsilon is of the same order of magnitude as the data sampling rate. We compute the discrete transition behavior (DTB) defined as the limit DTB(xv \) := lim\epsilon \rightarrow0 frec \epsilon (x0 + euroxv \), where x0 is generic. We illustrate the DTB by two sets of numerical experiments. In the first set, the perturbation is a smooth, rapidly oscillating sinusoid, and in the second, a fractal curve. The experiments reveal that the match between the DTB and reconstruction is worse as S\epsilon gets rougher. This is in agreement with the proof of the DTB, which suggests that the rate of convergence to the limit is O(euro\gamma /2). We then propose a new DTB, which exhibits an excellent agreement with reconstructions. Investigation of this phenomenon requires computing the rate of convergence for the new DTB. This, in turn, requires completely new approaches. We obtain a partial result along these lines and formulate a conjecture that the rate of convergence of the new DTB is O(euro1/2ln(1/euro)).