ANALYSIS OF RECONSTRUCTION FROM DISCRETE RADON TRANSFORM DATA IN R3 WHEN THE FUNCTION HAS JUMP DISCONTINUITIES

In this paper we study reconstruction of a function f from its discrete Radon transform data in R-3 when f 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(epsilon(2)). Let f, denote the result of reconstruction from the discrete data. Pick any generic point x(0) (i.e., satisfying some mild conditions), where f has a jump. Our first result is an explicit leading term behavior of f(epsilon) in an O(epsilon)-neighborhood of x(0) as epsilon -> 0. 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 x(0) is not an element of S := singsupp(f) is generic, then f epsilon (x(0)) -> f(x(0) ) as epsilon -> 0. The proof of this result reveals a surprising connection with the theory of uniform distribution. 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.