Local reconstruction analysis of inverting the Radon transform in the plane from noisy discrete data

Abstract

In this paper, we investigate the reconstruction error, Nrec(x), when a linear, filtered back-projection (FBP) algorithm is applied to noisy, discrete Radon transform data with sampling step size ε in two-dimensions. Specifically, we analyze Nrec(x) for x in small, O()-sized neighborhoods around a generic fixed point, x0, in the plane, where the measurement noise values, ηk,j (i.e., the errors in the sinogram space), are random variables. The latter are independent, but not necessarily identically distributed. We show, under suitable assumptions on the first three moments of the ηk,j, that the following limit exists: Nrec(;x0) = 0Nrec(x0+), for x in a bounded domain. Here, Nrec and Nrec are viewed as continuous random variables, and the limit is understood in the sense of distributions. Once the limit is established, we prove that Nrec is a zero mean Gaussian random field and compute explicitly its covariance. In addition, we validate our theory using numerical simulations and pseudo random noise.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…