Gabor wavelet

Gabor wavelets are wavelets invented by Dennis Gabor using complex functions constructed to serve as a basis for Fourier transforms in information theory applications. They are very similar to Morlet wavelets. They are also closely related to Gabor filters (see Gabor filter#Wavelet space). The important property of the wavelet is that it minimizes the product of its standard deviations in the time and frequency domain. Put another way, the uncertainty in information carried by this wavelet is minimized. However they have the downside of being non-orthogonal, so efficient decomposition into the basis is difficult. Since their inception, various applications have appeared, from image processing to analyzing neurons in the human visual system. [1] [2]

Minimal uncertainty property

The motivation for Gabor wavelets comes from finding some function f(x) which minimizes its standard deviation in the time and frequency domains. More formally, the variance in the position domain is:

(\Delta x)^2 = \frac {\int_{-\infty}^{\infty} (x-\mu)^2 f(x)f^{*}(x) \,dx} {\int_{-\infty}^{\infty} f(x)f^{*}(x) \, dx}

where f^{*}(x) is the complex conjugate of f(x) and \mu is the arithmetic mean, defined as:

\mu = \frac {\int_{-\infty}^{\infty} x f(x)f^{*}(x) \,dx} {\int_{-\infty}^{\infty} f(x)f^{*}(x)\,dx}

The variance in the wave number domain is:

(\Delta k)^2 = \frac {\int_{-\infty}^{\infty} (k-k_0)^2 F(k)F^{*}(k) \, dk} {\int_{-\infty}^{\infty} F(k)F^{*}(k) \, dk}

Where k_0 is the arithmetic mean of the Fourier Transform of f(x) , F(x) :

k_0 = \frac {\int_{-\infty}^{\infty} k F(k)F^{*}(k) \,dk} {\int_{-\infty}^{\infty} F(k)F^{*}(k) \,dk}

With these defined, the uncertainty is written as:

(\Delta x)(\Delta k)

This quantity has been shown to have a lower bound of \frac12 . The quantum mechanics view is to interpret (\Delta x) as the uncertainty in position and \hbar (\Delta k) as uncertainty in momentum. A function f(x) that has the lowest theoretically possible uncertainty bound is the Gabor Wavelet.[3]

Equation

The equation of a 1-D Gabor wavelet is a Gaussian modulated by a complex exponential, described as follows:[3]

f(x) = e^{-(x - x_0)^2/a^2}e^{-i k_0(x-x_0)}

As opposed to other functions commonly used as bases in Fourier Transforms such as \sin and \cos , Gabor wavelets have locality properties, meaning that as the distance from the center x_0 increases, the value of the function becomes exponentially suppressed. a controls the rate of this exponential drop-off and k_0 controls the rate of modulation.

It is also worth noting the Fourier transform of a Gabor wavelet, which is also a Gabor wavelet:

F(k) = e^{-(k - k_0)^2 a^2}e^{-i x_0(k-k_0)}

An example wavelet is given here:

See also

References

[1]

  1. ↑ 1.0 1.1 Lee, Tai S. (October 1996). "Image Representation Using 2D Gabor wavelets". IEEE Transactions on Pattern Analysis and Machine Intelligence (IEEE) 18 (10).
  2. ↑ Daugman, John. Computer Vision Lecture Series. University of Cambridge.
  3. ↑ 3.0 3.1 Daugman, John. Information Theory Lecture Series. University of Cambridge.