Morlet wavelet

Real-valued Morlet wavelet
Complex-valued Morlet wavelet

In mathematics, the Morlet wavelet (or Gabor wavelet)[1] is a wavelet composed of a complex exponential (carrier) multiplied by a Gaussian window (envelope). This wavelet is closely related to human perception, both hearing[2] and vision.[3]

History

In 1946, physicist Dennis Gabor, applying ideas from quantum physics, introduced the use of Gaussian-windowed sinusoids for time-frequency decomposition, which he referred to as atoms, and which provide the best trade-off between spatial and frequency resolution.[1] These are used in the Gabor transform, a type of short-time Fourier transform.[2] In 1984, Jean Morlet introduced Gabor's work to the seismology community and, with Goupillaud and Grossmann, modified it to keep the same wavelet shape over equal octave intervals, resulting in the first formalization of the continuous wavelet transform.[4] (See also Wavelet history)

Definition

The wavelet is defined as a constant \kappa_{\sigma} subtracted from a plane wave and then localised by a Gaussian window:[5]

\Psi_{\sigma}(t)=c_{\sigma}\pi^{-\frac{1}{4}}e^{-\frac{1}{2}t^{2}}(e^{i\sigma t}-\kappa_{\sigma})

where \kappa_{\sigma}=e^{-\frac{1}{2}\sigma^{2}} is defined by the admissibility criterion and the normalisation constant c_{\sigma} is:

c_{\sigma}=\left(1+e^{-\sigma^{2}}-2e^{-\frac{3}{4}\sigma^{2}}\right)^{-\frac{1}{2}}

The Fourier transform of the Morlet wavelet is:

\hat{\Psi}_{\sigma}(\omega) = c_\sigma \pi^{-\frac{1}{4}} \left( e^{-\frac{1}{2}(\sigma-\omega)^2} - \kappa_\sigma e^{-\frac{1}{2}\omega^{2}} \right)

The "central frequency" \omega_{\Psi} is the position of the global maximum of \hat{\Psi}_{\sigma}(\omega) which, in this case, is given by the solution of the equation:

(\omega_{\Psi}-\sigma)^{2}-1=(\omega_{\Psi}^{2}-1)e^{-\sigma\omega_{\Psi}}

The parameter \sigma in the Morlet wavelet allows trade between time and frequency resolutions. Conventionally, the restriction \sigma>5 is used to avoid problems with the Morlet wavelet at low \sigma (high temporal resolution).

For signals containing only slowly varying frequency and amplitude modulations (audio, for example) it is not necessary to use small values of \sigma. In this case, \kappa_{\sigma} becomes very small (e.g. \sigma>5 \quad \Rightarrow \quad \kappa_{\sigma}<10^{-5}\,) and is, therefore, often neglected. Under the restriction \sigma>5, the frequency of the Morlet wavelet is conventionally taken to be \omega_{\Psi}\simeq\sigma.

The wavelet exists as a complex version or a purely real-valued version. Some distinguish between the "real Morlet" vs the "complex Morlet".[6] Others consider the complex version to be the "Gabor wavelet", while the real-valued version is the "Morlet wavelet".[7][8]

Matlab function

[PSI,X] = morlet(LB,UB,N) returns values of the Morlet wavelet on an N point regular grid in the interval [LB,UB].

Output arguments are the wavelet function PSI computed on the grid X, and the grid X.

This wavelet has [-4 4] as effective support. Although [-4 4] is the correct theoretical effective support, a wider effective support, [-8 8], is used in the computation to provide more accurate results. [9]

Use in medicine

The Morlet wavelet transform method presented offers an intuitive bridge between frequency and time information which can clarify interpretation of complex head trauma spectra obtained with Fourier transform. The Morlet wavelet transform, however, is not intended as a replacement for the Fourier transform, but rather a supplement that allows qualitative access to time related changes and takes advantage of the multiple dimensions available in a free induction decay analysis.[10]

Application in the electrocardiogram

The application of the Morlet wavelet analysis in the electrocardiogram (ECG) is mainly to discriminate the abnormal heartbeat behavior. Since the variation of the abnormal heartbeat is a non-stationary signal, then this signal is suitable for wavelet-based analysis. [11]

See also

References

  1. ↑ 1.0 1.1 A Real-Time Gabor Primal Sketch for Visual Attention "The Gabor kernel satisfies the admissibility condition for wavelets, thus being suited for multi-resolution analysis. Apart from a scale factor, it is also known as the Morlet Wavelet."
  2. ↑ 2.0 2.1 Time-Frequency Dictionaries, Mallat
  3. ↑ J. G. Daugman. Uncertainty relation for resolution in space, spatial frequency, and orientation optimized by two-dimensional visual cortical filters. Journal of the Optical Society of America A, 2(7):1160–1169, July 1985.
  4. ↑ http://rocksolidimages.com/pdf/gabor.pdf
  5. ↑ John Ashmead (2012). "Morlet Wavelets in Quantum Mechanics". Quanta 1 (1): 58–70. doi:10.12743/quanta.v1i1.5.
  6. ↑ Matlab Wavelet Families - "Morlet Wavelet: morl" and "Complex Morlet Wavelets: cmor"
  7. ↑ Mathematica documentation: GaborWavelet
  8. ↑ Mathematica documentation: MorletWavelet
  9. ↑ "morlet". http://www.mathworks.com/''.
  10. ↑ http://cds.ismrm.org/ismrm-2001/PDF3/0822.pdf
  11. ↑ "A Tutorial of the Morlet Wavelet Transform". http://djj.ee.ntu.edu.tw/''.