Morlet wavelet explained

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 and vision.[2]

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.[3] 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]

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

-1
4
\pi
-1t2
2
e

(ei\sigma-\kappa\sigma)

where

\kappa\sigma

-1\sigma2
2
=e
is defined by the admissibility criterion,and the normalisation constant

c\sigma

is:

c\sigma

-\sigma2
=\left(1+e
-3\sigma2
4
-2e
-1
2
\right)

The Fourier transform of the Morlet wavelet is:

\hat{\Psi}\sigma(\omega)=c\sigma

-1
4
\pi

\left(

-1(\sigma-\omega)2
2
e

-\kappa\sigma

-1\omega2
2
e

\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 positive solution to:

\omega\Psi=\sigma

1
1-
-\sigma\omega\Psi
e

which can be solved by a fixed-point iteration starting at

\omega\Psi=\sigma

(the fixed-point iterations converge to the unique positive solution for any initial

\omega\Psi>0

).

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\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]

Uses

Use in medicine

In magnetic resonance spectroscopy imaging, the Morlet wavelet transform method offers an intuitive bridge between frequency and time information which can clarify the 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.[9]

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

Use in music

The Morlet wavelet transform is used in pitch estimation and can produce more accurate results than Fourier transform techniques.[10] The Morlet wavelet transform is capable of capturing short bursts of repeating and alternating music notes with a clear start and end time for each note.

A modified morlet wavelet was proposed to extract melody from polyphonic music.[11] This methodology is designed for the detection of closed frequency. The Morlet wavelet transform is able to capture music notes and the relationship of scale and frequency is represented as the follow:

fa={fc\overa x T}

where

fa

is the pseudo frequency to scale

a

,

fc

is the center frequency and

T

is the sampling time.

Morlet wavelet is modified as described as:

\Psi(t)=e-|{t|}cos(2\pit)

and its Fourier transformation:

F[\Psi(t)]={1\over{4\pi2f2+1}}[\delta(f-2\pi)+\delta(f+2\pi)]

Application

See also

References

  1. http://homepages.inf.ed.ac.uk/rbf/CAVIAR/PAPERS/05-ibpria-alex.pdf A Real-Time Gabor Primal Sketch for Visual Attention
  2. [John Daugman|J. G. Daugman]
  3. Book: Mallat, Stephane . A Wavelet Tour of Signal Processing, The Sparse Way . September 18, 2009 . Time-Frequency Dictionaries.
  4. Web site: Archived copy . 2012-05-12 . 2013-06-09 . https://web.archive.org/web/20130609045720/http://www.rocksolidimages.com/pdf/gabor.pdf . dead .
  5. John Ashmead . Morlet Wavelets in Quantum Mechanics . Quanta . 2012 . 1 . 1 . 58–70 . 10.12743/quanta.v1i1.5 . 1001.0250 . 73526961 .
  6. Web site: Matlab Wavelet Families . https://web.archive.org/web/20190810051900/https://www.mathworks.com/help/wavelet/ug/wavelet-families-additional-discussion.html . 2019-08-10 . live.
  7. Mathematica documentation: GaborWavelet
  8. Mathematica documentation: MorletWavelet
  9. http://cds.ismrm.org/ismrm-2001/PDF3/0822.pdf
  10. Kumar. Neeraj. Kumar. Raubin. 2020-01-29. Wavelet transform-based multipitch estimation in polyphonic music. Heliyon. 6. 1. e03243. 10.1016/j.heliyon.2020.e03243. free . 2405-8440. 7000807. 32042974. 2020Heliy...603243K .
  11. Kumar . Neeraj . Kumar . Raubin . Murmu . Govind . Sethy . Prabira Kumar . 2021-02-01 . Extraction of melody from polyphonic music using modified morlet wavelet . Microprocessors and Microsystems . 80 . 103612 . 10.1016/j.micpro.2020.103612 . 0141-9331.
  12. Modified Stacked Autoencoder Using Adaptive Morlet Wavelet for Intelligent Fault Diagnosis of Rotating Machinery. Haidong. Shao . Min . Xia . Jiafu . Wan . W. de Silva. Clarence. February 2022. IEEE/ASME Transactions on Mechatronics. 27 . 24–33 . 10.1109/TMECH.2021.3058061 .
  13. Designing of Morlet wavelet as a neural network for a novel prevention category in the HIV system. Sabir. Zulqurnain . Umar . Muhammad . Asif Zahoor Raja . Muhammad . Mehmet Baskonus. Haci. Wei. Gao. 2022. International Journal of Biomathematics. 15 . 4 . 10.1142/S1793524522500127 .
  14. Numerical computing to solve the nonlinear corneal system of eye surgery using the capability of Morlet wavelet artificial neural networks.. Wang, B. O. J. F. Gomez-Aguilar . Zulqurnain Sabir . Muhammad Asif Zahoor Raja. Wei-Feng Xia . H. A. D. I. Jahanshahi. Madini O. Alassafi. Fawaz E. Alsaadi. 2022. Fractals. 30 . 5 . 2240147–2240353 . 10.1142/S0218348X22401478 . free. 2022Fract..3040147W .
  15. LOS/NLOS Identification for Indoor UWB Positioning Based on Morlet Wavelet Transform and Convolutional Neural Networks. Z. Cui. Y. Gao . J. Hu . S. Tian. J. Cheng. March 2021. IEEE Communications Letters. 25 . 3 . 879–882 . 10.1109/LCOMM.2020.3039251 .
  16. Morlet Wavelet Filtering and Phase Analysis to Reduce the Limit of Detection for Thin Film Optical Biosensors. Simon J. Ward. Rabeb Layouni . Sofia Arshavsky-Graham . Ester Segal. Sharon M. Weiss. 2021. ACS Sensors. 6 . 8 . 2967–2978 . 10.1021/acssensors.1c00787 . 34387077 . 8403169.