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. The important property of the wavelet is that it minimizes the product of its standard deviations in the time and frequency domain (given by the variances defined below). 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 }xf(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 }kF(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 ${\frac {1}{2}}$ . 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^{-ik_{0}(x-x_{0})}

As opposed to other functions commonly used as bases in Fourier Transforms such as $\sin$ and $\cos$ , Gabor wavelets have the property that they are localized, 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 (unitary, angular-frequency convention) of a Gabor wavelet, which is also a Gabor wavelet:

F(k)=ae^{-(k-k_{0})^{2}a^{2}}e^{-ix_{0}(k-k_{0})}

An example wavelet is given here:

A Gabor wavelet with a = 2, x₀ = 0, and k₀ = 1

Time-causal analogue of the Gabor wavelet

When processing temporal signals, data from the future cannot be accessed, which leads to problems if attempting to use Gabor functions for processing real-time signals that depend upon the temporal dimension. A time-causal analogue of the Gabor filter has been developed in ^[4] based on replacing the Gaussian kernel in the Gabor function with a time-causal and time-recursive smoothing kernel referred to as the time-causal limit kernel. In this way, time-frequency analysis based on the resulting complex-valued extension of the time-causal limit kernel makes it possible to capture essentially similar transformations of a temporal signal as the Gabor wavelets can handle, and corresponding to the Heisenberg group, while carried out with strictly time-causal and time-recursive operations, see ^[4] for further details.

References

^ Lee, Tai S. (October 1996). "Image Representation Using 2D Gabor wavelets" (PDF). IEEE Transactions on Pattern Analysis and Machine Intelligence. 18 (10): 959–971. doi:10.1109/34.541406.
^ Daugman, John. Computer Vision Lecture Series (PDF). University of Cambridge.
^ ^a ^b Daugman, John. Information Theory Lecture Series (PDF). University of Cambridge.
^ ^a ^b Lindeberg, T. (23 January 2023). "A time-causal and time-recursive scale-covariant scale-space representation of temporal signals and past time". Biological Cybernetics: 1–39. arXiv:2202.09209. doi:10.1007/s00422-022-00953-6.

External links

MATLAB code for 2D Gabor wavelets and Gabor feature extraction

[imageRep-1] Lee, Tai S. (October 1996). "Image Representation Using 2D Gabor wavelets" (PDF). IEEE Transactions on Pattern Analysis and Machine Intelligence. 18 (10): 959–971. doi:10.1109/34.541406.

[CVRef-2] Daugman, John. Computer Vision Lecture Series (PDF). University of Cambridge.

[infRef-3] Daugman, John. Information Theory Lecture Series (PDF). University of Cambridge.

[Lin23-4] Lindeberg, T. (23 January 2023). "A time-causal and time-recursive scale-covariant scale-space representation of temporal signals and past time". Biological Cybernetics: 1–39. arXiv:2202.09209. doi:10.1007/s00422-022-00953-6.

[1]

[2]

[3]

[4]