Computationally efficient fourier transform of nonequidistant sampled data

We present an approximate method of performing the Fourier transform of the data sampled in nonequidistant readouts. It is shown that the data can be recalculated as equidistant readouts by using a nonuniform convolution, i.e., convolution of a certain function whose form depends on the calculated element and the character of nonequidistance. Thus, this recalculation does not require calculation of the values of the initial data in intermediate readouts (unlike the linear approximation, spline, or other recalculations). Since the size of the kernel of this nonuniform convolution is about 9, the proposed method can be the basis for an efficient computational algorithm. Applicability of the proposed approach to spectral optical coherence tomography is demonstrated.

[1]  N. Lomb Least-squares frequency analysis of unequally spaced data , 1976 .

[2]  S. K. Mitra,et al.  The nonuniform discrete Fourier transform and its applications in filter design. II. 2-D , 1996 .

[3]  Teresa C. Chen,et al.  Ultrahigh-resolution high-speed retinal imaging using spectral-domain optical coherence tomography. , 2004, Optics express.

[4]  J. Fujimoto,et al.  Optical Coherence Tomography , 1991 .

[5]  Maciej Wojtkowski,et al.  Complex and Coherence Noise Free Fourier Domain Optical Coherence Tomography , 2008 .

[6]  W Drexler,et al.  Ultrahigh resolution Fourier domain optical coherence tomography. , 2004, Optics express.

[7]  Brett E. Bouma,et al.  Optical Coherence Tomography and Coherence Techniques III , 2007 .

[8]  J. Tukey,et al.  An algorithm for the machine calculation of complex Fourier series , 1965 .

[9]  H. Bartelt,et al.  Image formation by inversion of scattered field data: experiments and computational simulation. , 1979, Applied optics.

[10]  Peter Koch,et al.  Using nonequispaced fast Fourier transformation to process optical coherence tomography signals , 2009, European Conference on Biomedical Optics.

[11]  A. Macovski,et al.  Selection of a convolution function for Fourier inversion using gridding [computerised tomography application]. , 1991, IEEE transactions on medical imaging.

[12]  A. Fercher,et al.  Measurement of intraocular distances by backscattering spectral interferometry , 1995 .

[13]  Sébastien Vergnole,et al.  Experimental validation of an optimized signal processing method to handle non-linearity in swept-source optical coherence tomography. , 2010, Optics express.