A Local Approach to Resolution Analysis of Image Reconstruction in Tomography

We propose a novel approach to analyzing resolution of tomographic reconstruction. Instead of following a conventional approach to obtain a global accuracy estimate, we investigate how the reconstructed function f_ε approximates the singularities of the original object f. The data is a discretized 2D Radon transform of f. The object could be static or change with time (dynamic tomography). Suppose the step-sizes along the angular and affine variables are O(ε). We pick a point x₀, where f has a jump singularity, and obtain the leading singular behavior of f_ε in an O(ε)-neighborhood of x₀ as ε→0. It turns out that the limiting behavior of f_ε depends only on the data microlocally near the singularity being reconstructed. This significantly simplifies the analysis and allows us to investigate complicated settings, e.g., dynamic tomography. Also, our resolution analysis is algorithm-specific---the same approach can be used for analyzing and optimizing various linear reconstruction algorithms. We present the results of numerical experiments in the static and dynamic cases. These results demonstrate an excellent agreement between predicted and actual behaviors of f_ε near a jump discontinuity of f.