Fact-checked by Grok 2 weeks ago

Wavelet transform

The wavelet transform is a mathematical technique used to decompose a function or signal into a set of basis functions called wavelets, which are localized in both time and frequency domains, enabling the analysis of non-stationary signals with varying frequency content over time.^[1] Unlike the Fourier transform, which provides global frequency information without temporal localization, the wavelet transform achieves multiresolution analysis by scaling and translating a mother wavelet function, allowing simultaneous examination of signal details at different resolutions.^[2] The development of wavelet theory traces its roots to early 20th-century work on harmonic analysis, including contributions from Alfred Haar in 1910 with the first wavelet basis, but gained modern prominence in the 1980s through the efforts of geophysicist Jean Morlet and mathematician Alex Grossmann, who introduced the continuous wavelet transform (CWT) for seismic signal analysis.^[3] Key advancements followed in the late 1980s and early 1990s, with Stéphane Mallat developing multiresolution analysis in 1989 and Ingrid Daubechies constructing compactly supported orthonormal wavelets in 1988, enabling efficient discrete implementations.^[2] There are two primary forms: the continuous wavelet transform (CWT), defined as W(a, b) = \frac{1}{\sqrt{a}} \int_{-\infty}^{\infty} f(t) \psi^*\left(\frac{t - b}{a}\right) dt, where a > 0 is the scale parameter, b is the translation parameter, and \psi is the mother wavelet satisfying admissibility conditions like zero mean; this form is ideal for theoretical analysis and visualization of time-scale representations.^[1] The discrete wavelet transform (DWT), often implemented via filter banks or the Mallat algorithm, discretizes scales and translations (typically dyadically, as powers of 2), facilitating fast computation through pyramid decomposition into approximation and detail coefficients, as pioneered by Daubechies' orthogonal wavelets.^[3]^[2] Wavelet transforms have broad applications across fields, including signal and image processing for compression (e.g., the JPEG 2000 standard), denoising, and feature extraction; biomedical engineering for ECG and EEG analysis; geophysics for seismic data interpretation; and control systems for motion tracking and fault detection.^[2] Their ability to handle transient phenomena and provide sparse representations has made them foundational in modern data analysis, with extensions like wavelet packets and lifting schemes addressing challenges in higher-dimensional or adaptive processing.^[3]

Fundamentals

Definition

The wavelet transform is a mathematical technique for decomposing a signal into components that reveal its behavior at different scales, achieved by representing the signal as a superposition of wavelets—localized, oscillatory functions generated by scaling and translating a single prototype function known as the mother wavelet. This approach provides simultaneous localization in both time and frequency domains, allowing for the analysis of non-stationary signals where features vary over time.^[1]^[4] In contrast to the Fourier transform, which decomposes signals into global sinusoidal basis functions and thus excels at frequency analysis but sacrifices temporal information for transient events, the wavelet transform uses compactly supported wavelets to capture localized, short-duration phenomena such as discontinuities or abrupt changes.^[5]^[4] The basic building block is the family of wavelets defined by

\psi_{a,b}(t) = \frac{1}{\sqrt{|a|}} \psi\left( \frac{t - b}{a} \right),

where \psi(t) denotes the mother wavelet, a \in \mathbb{R} \setminus \{0\} controls the dilation (scale), and b \in \mathbb{R} governs the shift (translation), ensuring the functions maintain constant energy across scales.^[1] The origins of wavelet theory trace back to Alfred Haar's 1910 introduction of the first explicit example of an orthogonal wavelet basis, though the modern continuous wavelet transform emerged in the 1980s through the foundational work of Alexandre Grossmann and Jean Morlet on decomposing functions into square-integrable wavelets, with further advancements by Stéphane Mallat in multiresolution frameworks.^[4]^[6]^[7] This entry presupposes basic familiarity with continuous-time signals and the Fourier transform as a global frequency decomposition tool.^[5]

Continuous Wavelet Transform

The continuous wavelet transform (CWT) provides a time-scale representation of a signal f \in L^2(\mathbb{R}) by correlating it with scaled and translated versions of a mother wavelet \psi. The transform is defined by the integral

W_f(a, b) = \int_{-\infty}^{\infty} f(t) \frac{1}{\sqrt{|a|}} \overline{\psi}\left( \frac{t - b}{a} \right) \, dt,

where a \neq 0 is the scale parameter controlling dilation, b \in \mathbb{R} is the translation parameter shifting the wavelet in time, and \overline{\psi} denotes the complex conjugate of the mother wavelet.^[8]^[1] This formulation arises from the inner product \langle f, \psi_{a,b} \rangle, where \psi_{a,b}(t) = \frac{1}{\sqrt{|a|}} \psi\left( \frac{t - b}{a} \right) is the dilated and translated wavelet, assuming \|\psi\|_2 = 1 and \int \psi(t) \, dt = 0 to ensure zero mean and oscillatory behavior.^[8]^[9] For the CWT to allow perfect reconstruction of the original signal, the mother wavelet must satisfy the admissibility condition

C_\psi = \int_{-\infty}^{\infty} \frac{|\hat{\psi}(\omega)|^2}{|\omega|} \, d\omega < \infty,

where \hat{\psi} is the Fourier transform of \psi, and \hat{\psi}(0) = 0.^[8]^[9] This condition, derived by Calderón and formalized by Grossmann and Morlet, ensures the wavelet has sufficient decay at low and high frequencies, preventing energy leakage and enabling invertibility; C_\psi serves as the admissibility constant in reconstruction.^[8] The original signal can be recovered from the CWT coefficients via the inverse formula

f(t) = \frac{1}{C_\psi} \int_{-\infty}^{\infty} \int_{0}^{\infty} W_f(a, b) \frac{1}{\sqrt{|a|}} \psi\left( \frac{t - b}{a} \right) \frac{da \, db}{a^2}.

^[8]^[9] This double integral over scale and translation reconstructs f pointwise, with the measure da \, db / a^2 reflecting the invariant structure of the affine group.^[9] The CWT exhibits several key properties that underpin its utility in analysis. It is linear, so W_{cf + g}(a, b) = c W_f(a, b) + W_g(a, b) for constants c and signals f, g.^[8] Shift invariance holds: translating the signal by \tau shifts the transform by \tau in the b variable, i.e., W_{f(\cdot - \tau)}(a, b) = W_f(a, b - \tau).^[8]^[1] Scale covariance adjusts accordingly: dilating the signal by \lambda > 0 transforms the CWT as W_{f(\lambda \cdot)}(a, b) = \lambda^{-1/2} W_f(\lambda a, \lambda b).^[8] Energy is preserved via a wavelet analog of Parseval's theorem:

\int_{-\infty}^{\infty} |f(t)|^2 \, dt = \frac{1}{C_\psi} \int_{-\infty}^{\infty} \int_{0}^{\infty} |W_f(a, b)|^2 \frac{da \, db}{a^2},

ensuring the transform maintains the signal's L^2 norm in the time-scale plane.^[8]^[9] Common examples include the Morlet wavelet, defined as \psi(t) = \pi^{-1/4} e^{i \omega_0 t} e^{-t^2/2} with central frequency \omega_0 \approx 6 to satisfy admissibility, which resembles a Gaussian-modulated plane wave and excels at detecting oscillatory features due to its balanced time-frequency localization.^[10]^[1] The Mexican hat wavelet, \psi(t) = (1 - t^2) e^{-t^2/2}, is the negative second derivative of a Gaussian, producing a shape with two lobes symmetric about zero; its Fourier transform \hat{\psi}(\omega) = -\omega^2 e^{-\omega^2/2} emphasizes mid-frequencies, making it suitable for detecting transients or edges in signals.^[1] Unlike the discrete wavelet transform, which samples scales and translations on a dyadic grid for computational efficiency and non-redundant representation, the CWT operates over continuous parameters, yielding a redundant but richly detailed transform ideal for theoretical analysis and visualization rather than fast numerical implementation.^[9]

Discrete Wavelet Transform

Mathematical Formulation

The discrete wavelet transform (DWT) arises as a dyadic discretization of the continuous wavelet transform, selecting scales a = 2^{-j} and translations b = k \cdot 2^{-j} for j, k \in \mathbb{Z}, which enables efficient computation on discrete signals while preserving the time-frequency localization properties.^[7] This leads to the wavelet coefficients defined as

c_{j,k} = 2^{j/2} \int_{-\infty}^{\infty} f(t) \psi\left(2^{j} t - k \right) \, dt,

where \psi is the mother wavelet and f is the input signal.^[7] This discretization forms the theoretical basis for the DWT, drawing from the continuous wavelet transform while adapting it to sampled data.^[7] In practice, the DWT is implemented via a two-channel filter bank, where the signal is decomposed using a low-pass scaling filter h (associated with the scaling function) and a high-pass wavelet filter g (associated with the wavelet), each followed by downsampling by a factor of 2.^[7] Starting from approximation coefficients a_j at scale j, the approximation coefficients at the next coarser scale j-1 are computed as

a_{j-1} = \sum_{l} h[l - 2m] \, a_j,

while the detail (wavelet) coefficients are

d_{j-1} = \sum_{l} g[l - 2m] \, a_j.

^[7] This recursive process builds a multiresolution decomposition, justified by the multiresolution analysis framework that underlies the dyadic structure.^[7] For perfect reconstruction of the original signal, the analysis and synthesis filters must satisfy conditions that eliminate aliasing and distortion, typically achieved through quadrature mirror filter designs where the low-pass filter H(z) and high-pass filter G(z) obey H(z) G(-z) - H(-z) G(z) = 2 z^{-l} for some integer delay l. The inverse DWT reconstructs the signal by upsampling the approximation and detail coefficients by 2, then applying the corresponding synthesis filters (often the time-reversed analysis filters for orthogonal bases) and summing the results.^[7] The filter bank implementation enables a fast pyramidal algorithm with computational complexity O(N) for an input signal of length N, as each level processes half the samples of the previous level, making the DWT suitable for real-time digital signal processing.^[7]

Multiresolution Analysis

Multiresolution analysis (MRA) provides the theoretical framework for decomposing signals in L^2(\mathbb{R}) into hierarchical approximations and details at dyadic scales, underpinning the discrete wavelet transform as its practical implementation. An MRA consists of a sequence of nested closed subspaces \{V_j\}_{j \in \mathbb{Z}} of L^2(\mathbb{R}) such that \cdots \subset V_{j} \subset V_{j+1} \subset \cdots, with \bigcap_j V_j = \{0\}, \bigcup_j V_j dense in L^2(\mathbb{R}), and a scaling invariance property: f \in V_j if and only if f(2 \cdot) \in V_{j+1}. Complementary detail subspaces W_j satisfy V_{j+1} = V_j \oplus W_j (orthogonal direct sum) and L^2(\mathbb{R}) = \bigoplus_j W_j.^[7] The scaling function \phi \in V_0 generates orthonormal bases for each V_j through integer translates and dyadic dilations: \phi_{j,k}(t) = 2^{j/2} \phi(2^j t - k), k \in \mathbb{Z}. This function obeys the two-scale dilation equation

\phi(t) = \sqrt{2} \sum_{k \in \mathbb{Z}} h_k \, \phi(2t - k),

where \{h_k\}_{k \in \mathbb{Z}} is a square-summable sequence of real coefficients satisfying \sum_k h_k = \sqrt{2} to preserve the integral of \phi.^[7] The wavelet function \psi, which spans the detail spaces W_j, is constructed from the scaling function via

\psi(t) = \sqrt{2} \sum_{k \in \mathbb{Z}} g_k \, \phi(2t - k),

with basis functions \psi_{j,k}(t) = 2^{j/2} \psi(2^j t - k), k \in \mathbb{Z}. In orthogonal MRAs, the high-pass coefficients relate to the low-pass ones by g_k = (-1)^k h_{1-k}.^[7] Essential properties include orthonormality of the scaling functions at level zero, \int_{-\infty}^{\infty} \phi(t - k) \phi(t - l) \, dt = \delta_{kl}, which propagates to all levels j due to the dilation structure. For the wavelet basis in orthogonal MRAs, completeness holds via Parseval's identity: for any f \in L^2(\mathbb{R}),

\|f\|_{L^2}^2 = \sum_{j \in \mathbb{Z}} \sum_{k \in \mathbb{Z}} |\langle f, \psi_{j,k} \rangle|^2,

ensuring exact reconstruction from wavelet coefficients.^[7] In more general MRAs, stability requires the collection \{\psi_{j,k}\}_{j,k \in \mathbb{Z}} to form a Riesz basis for L^2(\mathbb{R}), characterized by positive constants $0 < A \leq B < \infty such that

A \|f\|_{L^2}^2 \leq \sum_{j,k} |\langle f, \psi_{j,k} \rangle|^2 \leq B \|f\|_{L^2}^2

for all f \in L^2(\mathbb{R}); this bounded equivalence guarantees that coefficient sequences faithfully represent the signal without distortion or instability.^[11] The sequence \{h_k\} corresponds to the coefficients of the projection of \phi onto its own dilates and translates, serving as the impulse response of a low-pass filter in the equivalent filter bank representation of the MRA, while \{g_k\} defines the high-pass filter for extracting details.^[7] The Haar MRA illustrates these concepts simply, with scaling function \phi(t) = \chi_{[0,1)}(t) (the indicator function of [0,1)), satisfying h_0 = h_1 = 1/\sqrt{2} and all other h_k = 0. The corresponding wavelet is \psi(t) = \phi(2t) - \phi(2t - 1), a step function equal to 1 on [0, 1/2), -1 on [1/2, 1), and 0 elsewhere, with g_0 = 1/\sqrt{2}, g_1 = -1/\sqrt{2}.^[7]

Wavelet Construction and Properties

Mother Wavelet and Scaling

The mother wavelet, denoted as \psi(t), serves as the prototype function for generating wavelet bases in the wavelet transform. It is an oscillatory function in L^2(\mathbb{R}) with finite energy, satisfying \int_{-\infty}^{\infty} |\psi(t)|^2 \, dt = 1, and zero mean, \int_{-\infty}^{\infty} \psi(t) \, dt = 0.^[12] These properties ensure that \psi(t) is localized in both time and frequency domains, allowing it to capture transient features effectively without contributing to the overall average of the signal.^[12] The scaling function, or father wavelet, denoted as \phi(t), complements the mother wavelet by providing low-pass approximations in multiresolution analysis. It is typically a positive function that integrates to unity, \int_{-\infty}^{\infty} \phi(t) \, dt = 1, and forms an orthogonal basis for the approximation space V_0.^[12] This normalization ensures that \phi(t) represents the coarse-scale structure of signals, with \int_{-\infty}^{\infty} |\phi(t)|^2 \, dt = 1 for orthonormality.^[12] Wavelet bases are generated from the mother wavelet through dilation and translation operations, yielding \psi_{j,k}(t) = 2^{j/2} \psi(2^j t - k) for integers j, k \in \mathbb{Z}.^[12] Similarly, the scaling functions are \phi_{j,k}(t) = 2^{j/2} \phi(2^j t - k). These transformations scale the wavelet by factors of $2^j (with positive j corresponding to finer resolutions) and shift it by k, producing a family that spans L^2(\mathbb{R}). The resulting basis satisfies orthogonality conditions, \int_{-\infty}^{\infty} \psi_{j,k}(t) \psi_{m,n}(t) \, dt = \delta_{j m} \delta_{k n}, ensuring unique decompositions.^[12] Key properties of the mother wavelet include compact support, where \psi(t) is nonzero only on a finite interval, enabling efficient finite impulse response filters in implementations.^[12] Additionally, vanishing moments up to order M-1, defined by \int_{-\infty}^{\infty} t^m \psi(t) \, dt = 0 for m = 0, 1, \dots, M-1, allow wavelets to approximate polynomials well in scaling spaces, enhancing the transform's ability to represent smooth components.^[12] A simple example is the Haar mother wavelet, \psi(t) = 1 for $0 \leq t < 0.5, \psi(t) = -1 for $0.5 \leq t < 1, and \psi(t) = 0 otherwise, which has compact support on [0,1) and one vanishing moment.^[12] This discontinuous function exemplifies the basic construction, with its corresponding scaling function \phi(t) = 1 for $0 \leq t < 1 and 0 elsewhere.^[12]

Orthogonal Wavelet Families

Orthogonal wavelets must satisfy the two-scale dilation equation, ensuring the scaling function and wavelet form a multiresolution analysis, while orthogonality requires that the translates and dilates are orthonormal, and compact support implies finite impulse response (FIR) filters for computational efficiency. These conditions enable perfect reconstruction in discrete wavelet transforms without redundancy.^[13] Daubechies wavelets, introduced in 1988, provide the shortest filters achieving a specified number of vanishing moments for compactly supported orthogonal wavelets, maximizing regularity for given support length. Their construction involves solving the two-scale relation through spectral factorization of the polynomial P(z) = \sum_{k} |H(\omega + 2\pi k)|^2, where H(\omega) is the low-pass filter, chosen as a maximally flat half-band filter to ensure the desired vanishing moments. For the Daubechies D4 wavelet (with 2 vanishing moments and 4-tap filters), the low-pass coefficients are:

\begin{align*} h_0 &= \frac{1 + \sqrt{3}}{4\sqrt{2}}, \\ h_1 &= \frac{3 + \sqrt{3}}{4\sqrt{2}}, \\ h_2 &= \frac{3 - \sqrt{3}}{4\sqrt{2}}, \\ h_3 &= \frac{1 - \sqrt{3}}{4\sqrt{2}}. \end{align*}

The high-pass coefficients are derived by alternating signs and reversal: g_k = (-1)^k h_{3-k}.^[13] Biorthogonal wavelets, developed by Cohen, Daubechies, and Feauveau in 1992, relax strict orthogonality by using dual wavelet functions \psi and \tilde{\psi} that satisfy \langle \psi_{j,k}, \tilde{\psi}_{j',k'} \rangle = \delta_{j,j'} \delta_{k,k'}, enabling perfect reconstruction while allowing symmetric filters for reduced phase distortion.^[14] This duality permits separate analysis and synthesis filters, often with unequal lengths, balancing smoothness and support.^[14] Other notable families include Coiflets, constructed by Daubechies to balance vanishing moments between the scaling function and wavelet, achieving approximate shift-invariance and better polynomial approximation.^[13] Symlets, also by Daubechies, modify the Daubechies filters through least-asymmetric optimization to increase symmetry while preserving orthogonality and minimal support.^[13] A key trade-off in these families is that increasing the number of vanishing moments enhances smoothness and approximation order but extends filter length, potentially reducing time localization in the wavelet transform.

Applications

Signal and Image Compression

Wavelet-based compression leverages the discrete wavelet transform (DWT) as the core decomposition tool to represent signals and images in a multi-resolution framework, enabling efficient data reduction through transform coding. The process begins by applying the DWT to decompose the input into wavelet coefficients, which capture both spatial and frequency information. Small coefficients, often representing noise or fine details, are then thresholded using techniques such as hard thresholding (setting coefficients below a threshold to zero) or soft thresholding (shrinking coefficients toward zero), followed by quantization to further reduce precision and entropy coding to compress the resulting bitstream. This approach exploits the energy compaction property of orthogonal wavelet families, such as Daubechies wavelets, where most signal energy is concentrated in a few large coefficients.^[15] A seminal method in this domain is the Embedded Zerotree Wavelet (EZW) algorithm, introduced by Jerome M. Shapiro in 1993, which exploits inter-scale dependencies in the wavelet coefficients to achieve progressive coding. In EZW, the DWT coefficients are organized into a hierarchical tree structure called zerotrees, where a coefficient at a coarser scale and its descendants at finer scales are predicted to be insignificant (below the threshold) together if the parent is insignificant, allowing efficient scanning and encoding of significant coefficients in embedded bit planes from most to least significant bits. This zerotree concept significantly reduces redundancy and enables scalable compression, where partial bitstreams yield approximations at varying quality levels.^[15] Building on EZW, the Set Partitioning in Hierarchical Trees (SPIHT) algorithm, developed by Amir Said and William A. Pearlman in 1996, improves efficiency by using spatial orientation trees instead of simple zerotrees, incorporating sign coding and more refined partitioning of coefficient sets. SPIHT refines the significance map testing to avoid the need for explicit lists, reducing memory usage and computational overhead while maintaining embedded coding properties and outperforming EZW in compression ratios for the same distortion levels. Its tree-based structure better captures the statistical dependencies across scales and orientations, making it particularly effective for images with diverse textures. Performance of wavelet compression is typically evaluated using metrics such as Peak Signal-to-Noise Ratio (PSNR), which measures reconstruction fidelity, and compression ratio, indicating data reduction efficiency. For instance, at high compression ratios above 30:1, wavelet methods achieve higher PSNR values compared to Discrete Cosine Transform (DCT)-based JPEG, due to superior energy compaction that preserves edges and textures without introducing prominent block artifacts. For example, tests on the Lena image show wavelet coders like SPIHT attaining a PSNR of 37.12 dB at 0.5 bits per pixel, surpassing JPEG's approximately 33 dB under similar conditions.^[16] To address shift-invariance issues inherent in the decimated DWT, which can cause coefficient misalignment and degrade compression for shifted inputs, undecimated wavelet transforms are employed, retaining all coefficients without downsampling to maintain translation invariance. Additionally, ringing artifacts—oscillatory effects near edges akin to the Gibbs phenomenon—can arise from abrupt coefficient truncation but are minimized by selecting smooth mother wavelets with higher regularity, which reduce high-frequency oscillations in the basis functions.^[17]^[18] Despite these advances, wavelet compression has limitations, including higher computational costs for long filters in the DWT, which increase encoding time for large datasets, and potential block artifacts in tiled image processing, though less severe than in DCT methods. It is also less effective for very smooth signals, where energy spreads more evenly across coefficients, leading to poorer compaction compared to frequency-domain transforms optimized for stationary content.^[19] A prominent real-world application is the JPEG 2000 standard, which adopts the biorthogonal 9/7 wavelet (CDF 9/7) for lossy compression, combining it with tiered entropy coding to support resolutions from 0.25 to 16 bits per pixel while achieving superior visual quality over JPEG at equivalent rates. This choice of wavelet balances computational efficiency and reconstruction quality, enabling features like region-of-interest coding and error resilience.^[20]

Denoising and Other Signal Processing

Wavelet denoising techniques leverage the sparse representation properties of wavelet transforms to remove noise from signals while preserving important features. In the 1990s, David Donoho introduced wavelet shrinkage methods, which apply nonlinear thresholding to wavelet coefficients to suppress noise. These approaches decompose the noisy signal into wavelet coefficients, threshold those dominated by noise, and reconstruct the signal via inverse transform. To estimate the noise variance σ, one computes the median of the absolute values of the detail coefficients at the finest scale and divides by 0.6745, approximating the standard deviation for Gaussian noise. A common threshold is the universal value λ = σ √(2 log N), where N is the signal length, which ensures near-optimal performance for worst-case scenarios. Several thresholding strategies build on this foundation. VisuShrink employs the universal threshold across all coefficients, providing a simple, globally adaptive method that achieves minimax risk rates for signals in Besov spaces. SureShrink, developed by Donoho and Johnstone, uses a level-dependent approach based on Stein's unbiased risk estimate (SURE), applying subband-adaptive thresholds to minimize estimated mean squared error.^[21] BayesShrink adopts a model-based perspective, assuming a generalized Gaussian prior for signal coefficients and deriving a data-driven threshold per subband to minimize Bayesian risk under this model. Beyond denoising, wavelets enable edge detection by identifying singularities in signals. The modulus maxima method examines the absolute values of wavelet coefficients at fine scales; sharp edges correspond to localized maxima chains across scales, allowing detection of discontinuities without smoothing them excessively. This technique, rooted in multiscale analysis, locates edges by tracking these maxima, which persist along ridges in the wavelet transform plane. For feature extraction, wavelet packets extend the standard dyadic decomposition by allowing adaptive frequency partitioning. This creates a library of bases through full binary tree decompositions of both approximation and detail spaces, enabling finer time-frequency tilings. Best basis selection minimizes an information cost function, such as Shannon entropy, to identify the optimal subspace for representing signal features with maximal sparsity. In biomedical signal analysis, wavelet denoising facilitates QRS complex detection in electrocardiograms (ECGs) by isolating the characteristic peaks amid baseline wander and muscle artifacts.^[22] Similarly, in geophysical applications, these methods enhance seismic reflection data by attenuating random noise and ground roll, improving the visibility of subsurface reflectors for structural interpretation.^[23] Wavelets offer advantages in processing non-stationary signals through their sparse representations, which concentrate energy around transients and handle discontinuities more effectively than Fourier methods by localizing both time and frequency information. However, limitations include sensitivity to threshold selection, where suboptimal choices can increase mean squared error or lead to over-smoothing of genuine signal features.^[21]

Comparisons and Extensions

Comparison with Fourier Transform

The Fourier transform decomposes signals into global sinusoidal basis functions of the form e^{i \omega t}, providing excellent frequency resolution but no inherent time localization, making it ideal for stationary signals where frequency content is consistent over time.^[24] In contrast, the wavelet transform employs localized basis functions \psi_{a,b}(t) = \frac{1}{\sqrt{a}} \psi\left(\frac{t - b}{a}\right), where a is the scale and b the translation, allowing simultaneous time and frequency analysis with adaptive resolution.^[6] This localization arises because the Fourier transform of the mother wavelet \Psi(\omega) decays away from its center frequency, ensuring compact support in both domains.^[24] A key limitation of the Fourier transform for non-stationary signals is addressed by the short-time Fourier transform (STFT), which applies a fixed-width window to achieve time localization; however, this fixed window enforces a trade-off dictated by the Heisenberg uncertainty principle, where improved time resolution degrades frequency resolution uniformly across scales, lacking multi-resolution capability.^[25] Wavelets overcome this by using a variable window: wider for low frequencies to enhance frequency resolution and narrower for high frequencies to improve time resolution, enabling better capture of transients such as impulses or chirps in signals with varying frequency content.^[25] The development of wavelet transforms in the 1980s was largely motivated by these inadequacies of Fourier-based methods for analyzing non-stationary data, particularly in fields like seismology.^[6] In practical examples, the Fourier transform struggles with discontinuous signals, producing the Gibbs phenomenon—persistent oscillations near jump discontinuities due to the global nature of its basis, with overshoots reaching about 9% of the jump height regardless of truncation level. Wavelets, with their localized support, yield sparse coefficients that avoid such ringing, providing more accurate reconstructions around discontinuities. Time-frequency representations further illustrate the differences: the STFT spectrogram |\text{STFT}_f(t,\omega)|^2 offers constant resolution cells, limiting adaptability to frequency variations, whereas the wavelet scalogram |W_f(a,b)|^2 from the continuous wavelet transform provides scale-dependent resolution, yielding sharper depictions for signals with evolving frequencies like chirps.^[26] Synchrosqueezed variants can refine this further for even sharper reassignment.^[26] Hybrid approaches combining wavelets and Fourier transforms are effective for wideband signals, where wavelets handle local transients and Fourier captures global low-frequency components.^[27]

Time-Causal and Synchrosqueezed Variants

Time-causal wavelets emerged in the 1990s as adaptations of traditional wavelet transforms to handle streaming data in real-time applications, ensuring that the analysis relies solely on past and present signal values without accessing future data.^[28] These progressive wavelets feature infinite support extending to the left (past) but are designed with compact or decaying support to the right, preserving causality for temporal signals.^[29] Unlike standard discrete wavelet transforms (DWT), which may introduce non-causal filtering delays, time-causal variants avoid "peeking" into the future, making them suitable for online processing.^[28] The formulation of time-causal wavelets requires their Fourier transforms to be supported exclusively on positive frequencies, ensuring an analytic signal representation that aligns with causality principles.^[29] For instance, Cauchy wavelets, defined in the frequency domain as \hat{\psi}(\omega) = \frac{c^{1/2} \Gamma(c+1) e^{i \theta} }{(\omega + i c)^{c+1}} for \omega > 0 and zero otherwise (with c > 0 controlling decay), provide progressive support and are used for directional analysis in causal contexts.^[30] Similarly, the progressive Morlet wavelet modifies the standard Morlet by restricting its Gaussian-modulated exponential to positive frequencies, yielding \psi(t) = \pi^{-1/4} e^{i k_0 t} e^{-t^2/2} with the understanding that negative frequency components are suppressed.^[29] In applications, time-causal wavelets enable real-time signal processing in domains like control systems, where low-latency analysis of dynamic inputs is critical, such as in adaptive filtering for feedback loops.^[28] They outperform non-causal DWT in streaming scenarios by maintaining temporal order without buffering future samples.^[31] However, these variants often sacrifice orthogonality for causality, leading to higher redundancy compared to symmetric orthogonal families.^[29] The synchrosqueezed wavelet transform (SWT), introduced by Daubechies, Lu, and Wu in 2009, refines the continuous wavelet transform (CWT) by reassigning its energy to instantaneous frequencies, yielding sharper time-frequency representations.^[32] This post-processing step computes the instantaneous frequency as \omega(a,b) = -i \frac{\partial_b W_f(a,b)}{W_f(a,b)}, where W_f(a,b) is the CWT coefficient at scale a and time b, concentrating energy along ridges corresponding to signal components.^[32] SWT produces sharper ridges than the standard CWT, effectively reducing blurring in non-stationary signals like chirps, where frequency varies linearly with time.^[33] Key properties of SWT include its invertibility, allowing perfect reconstruction of the original signal from the reassigned representation, similar to CWT but with enhanced localization.^[33] This makes SWT valuable for mode decomposition in multicomponent signals, akin to empirical mode decomposition but with rigorous mathematical foundations.^[32] Within wavelet frameworks, SWT improves time-frequency resolution without increasing the redundancy inherent in CWT, unlike causal variants that trade symmetry for latency.^[33] Limitations arise from its computational intensity, as the partial derivative computation scales with the number of scales and times, often requiring optimized implementations for large datasets.^[34]

References

[1]
[PDF] THE CONTINUOUS WAVELET TRANSFORM: A TOOL FOR SIGNAL ...
In this article, the continuous wavelet transform is introduced as a signal processing tool for investigating time-varying frequency spectrum characteristics ...
[2]
[PDF] Wavelet theory and applications - Pure
Jan 1, 2005 · Since then a lot of research into wavelets and the wavelet transform is performed. This report gives an overview of the main wavelet theory.<|control11|><|separator|>
[3]
(PDF) A Review of Wavelet Analysis and Its Applications
This article reviews the development history of wavelet theory, from the construction method to the discussion of wavelet properties.
[4]
Time Frequency Analysis of Wavelet and Fourier Transform
Continuous wavelet transform. Morlet and Grossman formalized the continuous wavelet transform in 1984 in [21]. ... 1910 by Alfred Haar as an alternative to the ...Missing: seminal | Show results with:seminal
[5]
[PDF] The wavelet transform, time-frequency localization and signal analysis
The first procedure is the short-time or windowed Fourier transform, the second is the “wavelet transform,” in which high frequency components are studied with ...
[6]
Decomposition of Hardy Functions into Square Integrable Wavelets ...
Decomposition of Hardy Functions into Square Integrable Wavelets of Constant Shape. Authors: A. Grossmann and J. Morlet ...
[7]
[PDF] A Theory for Multiresolution Signal Decomposition: The Wavelet ...
The wavelet representation can also be interpreted as a texton decom- position where each texton is equal to a particular func- tion of the wavelet orthonormal ...
[8]
[PDF] 4.3 Wavelet Transforms
Theorem 4.4 (Calderón, Grossmann, Morlet). Let ψ ∈ L2(R) be a real function ... is called the wavelet admissibility condition. To guarantee that this ...
[9]
https://epubs.siam.org/doi/10.1137/1031126
[10]
[PDF] A Practical Guide to Wavelet Analysis - NOAA
Wavelet analysis is becoming a common tool for analyzing localized variations of power within a time series. By decomposing a time series into time–fre-.
[11]
[PDF] Discrete Wavelets and Multiresolution Analysis
Abstract. In this paper we present an elementary discussion of the discrete wavelet transform. A major problem is formed by the construction of an or-.
[12]
[PDF] Orthonormal bases of compactly supported wavelets
We conclude this paper with the plots of a few of the compactly supported wavelets constructed here. 2. Multiresolution Analysis and Image Decomposition and ...
[13]
Ten Lectures on Wavelets | SIAM Publications Library
This monograph contains 10 lectures presented by Dr. Daubechies as the principal speaker at the 1990 CBMS-NSF Conference on Wavelets and Applications.
[14]
Biorthogonal bases of compactly supported wavelets
Orthonormal bases of compactly supported wavelet bases correspond to subband coding schemes with exact reconstruction in which the analysis and synthesis ...<|separator|>
[15]
Embedded image coding using zerotrees of wavelet coefficients
Abstract: The embedded zerotree wavelet algorithm (EZW) is a simple, yet remarkably effective, image compression algorithm, having the property that the ...
[16]
[PDF] Comparison of DCT and Wavelet based Image Compression ...
So the Discrete Wavelet Transform (DWT) is applied to an image and the energy compaction performance of both Discrete Cosine Transform (DCT) and DWT is compared ...
[17]
Non-decimated Wavelet Transform for a Shift-invariant Analysis
Aug 7, 2025 · On the other hand, the Non-Decimated Discrete Wavelet Transform is shift-invariant because it eliminates the downsampling and, consequently, is ...
[18]
Blur and Ringing Artifact Measurement in Image Compression using ...
In the spatial domain this causes ripples or oscillations around sharp edges or contours in the image. This is also known as the Gibbs phenomenon. A method for ...
[19]
Wavelet Compression - an overview | ScienceDirect Topics
7 Standard wavelet transforms may require many coefficients to represent jumps or edges, which can lead to artifacts in compressed images.Missing: Limitations | Show results with:Limitations
[20]
[PDF] Image compression using wavelets and JPEG2000: a tutorial
This paper presents a tutorial on the discrete wavelet transform (DWT) and introduces its application to the new JPEG2000* image compression standard. We start ...
[21]
Adapting to Unknown Smoothness via Wavelet Shrinkage
We attempt to recover a function of unknown smoothness from noisy sampled data. We introduce a procedure, SureShrink, that suppresses noise by thresholding the ...Missing: seminal | Show results with:seminal
[22]
A Wavelet Denoising and Teager Energy Operator-Based Method ...
Apr 7, 2020 · The major aim of this work is to develop a novel method for QRS complex detection under various ECG signal morphologies as well as under different ECG ...
[23]
Physical Wavelet Frame Denoising | GEOPHYSICS - SEG Library
The new approach, physical wavelet frame denoising uses a wavelet frame that takes into account the characteristics of seismic data both in time and space.
[24]
[PDF] Wavelet transforms and orthonormal wavelet bases
Bases. INGRID DAUBECHIES. ABSTRACT. We introduce the wavelet transform and discuss its motivation as a time-frequency localization tool. We review the different ...
[25]
[PDF] Comparison of STFT and Wavelet Transform in Time - DiVA portal
Wavelet transform provides a tool for both time –frequency localization and ... As shows in Fig 3, the Fourier transform can represent the frequency localization, ...
[26]
[PDF] Comparison of different time-frequency representations. - Hal-Inria
Dec 8, 2014 · Continuous wavelet transform and scalograms. Similarly with STFT and spectrogram, scalogram is squared value of continuous wavelet trans-.
[27]
https://ieeexplore.ieee.org/document/6663389
[28]
Causal analytical wavelet transform - SPIE Digital Library
Sep 1, 1992 · The causal analytical wavelet transform employs exponentially decaying nonsinusoidal wideband transient bases of compact support.
[29]
WAVELET TRANSFORMS AND THEIR APPLICATIONS TO ...
They are well adapted to analyze causal signals, i.e. signals for which there is some action of causality. This is because progressive wavelets preserve the ...
[30]
[PDF] Phase retrieval for the Cauchy wavelet transform - arXiv
Apr 3, 2017 · Abstract. We consider the phase retrieval problem in which one tries to reconstruct a function from the modulus of its wavelet transform.
[31]
[2510.05834] Time-causal and time-recursive wavelets - arXiv
Oct 7, 2025 · This paper describes how a time-causal wavelet analysis can be performed based on concepts developed in the area of temporal scale-space theory, ...Missing: Weiss | Show results with:Weiss
[32]
Synchrosqueezed Wavelet Transforms: a Tool for Empirical Mode ...
Dec 12, 2009 · In this paper we describe a method that captures the flavor and philosophy of the EMD approach, albeit using a different approach in constructing the ...Missing: 2008 | Show results with:2008
[33]
An empirical mode decomposition-like tool - ScienceDirect.com
Synchrosqueezed wavelet transforms: An empirical mode decomposition-like tool. Author links open overlay panelIngrid Daubechies , Jianfeng Lu , Hau-Tieng Wu.
[34]
[PDF] SYNCHROSQUEEZING TRANSFORM AND ITS APPLICATIONS
I.Daubechies, J .Lu, and H.-T. Wu, “Synchrosqueezed wavelet transforms: an empirical mode decomposition-like.