Abstract
The correct formulation of the plane-wave decomposition of point-source seismograms involves a sequence of Fourier and Hankel transforms which can be evaluated in several ways. The procedure which is proposed here exploits the fact that the plane-wave response is bandlimited along the horizontal slowness axis. This property permits to expand the Hankel transform into a Fourier-Bessel series. In practice, this algorithm requires an interpolation in distance of the recorded dataset, but allows a considerable reduction of Bessel functions calculations. Numerical applications performed with synthetic and real data show that the Fourier-Bessel summation technique yields results which are equivalent to a discrete integration of the Hankel transform.
