Resolution of 2D reconstruction of functions with nonsmooth edges from discrete Radon transform data
share
Let f be an unknown function in ℝ^2, and f_ϵ be its reconstruction from discrete Radon transform data, where ϵ is the data sampling rate. We study the resolution of reconstruction when f has a jump discontinuity along a nonsmooth curve 𝒮_ϵ. The assumptions are that (a) 𝒮_ϵ is an O(ϵ)-size perturbation of a smooth curve 𝒮, and (b) 𝒮_ϵ is Holder continuous with some exponent γ∈(0,1]. We compute the Discrete Transition Behavior (or, DTB) defined as the limit DTB(x̌):=lim_ϵ→0f_ϵ(x_0+ϵx̌), where x_0 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 𝒮_ϵ gets more rough. This is in agreement with the proof of the DTB, which suggests that the rate of convergence to the limit is O(ϵ^γ/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(ϵ^1/2ln(1/ϵ)).
READ FULL TEXT