next up previous [pdf]

Next: Continuous case and seismic Up: Interpolation with Fourier basis Previous: Continuous Fourier basis

Discrete Fourier basis

Assuming that the range of the variable $x$ is limited in the interval from $-N$ to $N$, the discrete Fourier basis (Fast Fourier Transform) employs a set of orthonormal periodic functions

\begin{displaymath}
\psi_k (x) = \frac{1}{\sqrt{2N}} e^{i \pi \frac{k}{N} x} \;,
\end{displaymath} (23)

where the discrete frequency index $k$ also ranges, according to the Nyquist sampling criterion, from $-N$ to $N$. The interpolation function is computed from equation (15) to be
$\displaystyle W (x, n)$ $\textstyle =$ $\displaystyle \frac{1}{2 N} \sum_{k=-N}^{N-1} e^{i \pi \frac{k}{N} (x-n)} =
\fr...
...1 + e^{i \pi \frac{x-n}{N}} + \cdots + e^{i \pi \frac{2N-1}{N} (x-n)}
\right] =$  
    $\displaystyle \frac{1}{2 N} e^{- i \pi (x-n)}
\frac{e^{2i \pi (x-n)} - 1}{e^{i ...
... - e^{-i \pi (x-n)}}{e^{i \pi \frac{x-n}{2 N}} -
e^{- i \pi \frac{x-n}{2 N}}} =$  
    $\displaystyle e^{-i \pi \frac{x-n}{2 N}}
\frac{\sin \left[\pi (x - n)\right]}{2N \sin\left[\pi (x - n)/2N\right]} \;.$ (24)

An interpolation function equivalent to (24) has been found by Muir (Popovici et al., 1993,1996; Lin et al., 1993). It can be considered a tapered version of the sinc interpolant (21) with smooth tapering function

\begin{displaymath}
\frac{\pi (x - n)/2N}{\tan\left[\pi (x - n)/2N\right]}\;.
\end{displaymath}

Unlike most other tapered-sinc interpolants, Muir's interpolant (24) satisfies not only the obvious property (2), but also properties (3) and (19), where the interpolation function $W (x, n)$ should be set to zero for $x$ outside the range from $n - N$ to $n+N$. The form of this function is shown in Figure 11.

ma-sinc
ma-sinc
Figure 11.
The left plots show the sinc interpolation function. Note the slow decay in $x$. The middle shows the effective tapering function of Muir's interpolation; the right is Muir's interpolant. The top is for $N=2$ (5-point interpolation); the bottom, $N=6$ (13-point interpolation).
[pdf] [png] [sage]

The development of the mathematical wavelet theory (Daubechies, 1992) has opened the door to a whole universe of orthonormal function bases, different from the Fourier basis. The wavelet theory should find many useful applications in geophysical data interpolation, but exploring this interesting opportunity would go beyond the scope of the present work.

The next section carries the analysis to the continuum and compares the mathematical interpolation theory with the theory of seismic imaging.


next up previous [pdf]

Next: Continuous case and seismic Up: Interpolation with Fourier basis Previous: Continuous Fourier basis

2014-02-21