A New Method of Solving Noisy Abel-Type Equations

Abstract

A new approach to solving noisy integral equations of the first kind is applied to the family of Abel equations. Such equations play a role in stereology (Wicksell's unfolding problem), medicine, engineering, and astronomy. The method is based on an expansion in an arbitrary orthonormal basis, coupled with exact inversion of the integral operator. The inverse appears in the Fourier coefficients of the expansion, where it can be carried over to the usually well-behaved basis elements in the form of the adjoint. This method is an alternative to Tikhonov regularization, regularization of the inverse of the operator itself, or a wavelet-vaguelette/singular-value decomposition. The method is particularly interesting in irregularity of the kernel, the input, or both. Because knowledge of the spectral properties of the operator is not required, the method is also of interest in regular cases where these spectral properties are not sufficiently known or are hard to deal with. For smooth input functions, the simple basis of trigonometric functions yields input estimators whose mean integrated squared error converges at the optimal rate for the entire family of Abel operators. This can be shown when smooth wavelets are used for Abel operators with index smaller than 1/2, and when the Haar wavelet is used for operators with index larger than 1/2.