Recursive methods for computing the Abel transform and its inverse

Abstract

The Abel transform and its inverse appear in a wide variety of problems in which it is necessary to reconstruct axisymetric functions from line-integral projections. We present a new family of algorithms, principally for Abel inversion, that are recursive and hence computationally efficient. The methods are based on a linear, space-variant, state-variable model of the Abel transform. The model is the basis for deterministic algorithms, applicable when data are noise free, and least-squares-estimation (Kalman filter) algorithms, which accommodate the noisy data case. Both one-pass (filtering) and two-pass (smoothing) estimators are considered. In computer simulations, the new algorithms compare favorably with previous methods for Abel inversion.