Article de reference

Wavelet transform

An example of the 2D discrete wavelet transform that is used in JPEG2000 In mathematics , a wavelet series is a representation of a square-integrable ( real - or complex -valued...

An example of the 2D discrete wavelet transform that is used in JPEG2000

In mathematics, a wavelet series is a representation of a square-integrable (real- or complex-valued) function by a certain orthonormalseries generated by a wavelet. This article provides a formal, mathematical definition of an orthonormal wavelet and of the integral wavelet transform.

Definition

A function orthonormal wavelet if it can be used to define a Hilbert basis, that is, a complete orthonormal system for the Hilbert space of square-integrable functions on the real line.

The Hilbert basis is constructed as the family of functions

If, under the standard inner product on

Completeness is satisfied if every function

with convergence of the series understood to be convergence in norm. Such a representation of wavelet series. This implies that an orthonormal wavelet is self-dual.

The integral wavelet transform is the integral transform defined as wavelet coefficients

Here, binary dilation or dyadic dilation, and binary or dyadic position.

Principle

The fundamental idea of wavelet transforms is that the transformation should allow only changes in time extension, but not shape, imposing a restriction on choosing suitable basis functions. Changes in the time extension are expected to conform to the corresponding analysis frequency of the basis function. Based on the uncertainty principle of signal processing,

where

The higher the required resolution in time, the lower the resolution in frequency has to be. The larger the extension of the analysis windows is chosen, the larger is the value of

When

  1. Bad time resolution
  2. Good frequency resolution
  3. Low frequency, large scaling factor

When

  1. Good time resolution
  2. Bad frequency resolution
  3. High frequency, small scaling factor

In other words, the basis function

This shows that wavelet transformation is good in time resolution of high frequencies, while for slowly varying functions, the frequency resolution is remarkable.

Another example: The analysis of three superposed sinusoidal signals

Wavelet compression

Wavelet compression is a form of data compression well suited for image compression (sometimes also video compression and audio compression). Notable implementations are JPEG 2000, DjVu and ECW for still images, JPEG XS, CineForm, and the BBC's Dirac. The goal is to store image data in as little space as possible in a file. Wavelet compression can be either lossless or lossy.

Method

First a wavelet transform is applied. This produces as many coefficients as there are pixels in the image (i.e., there is no compression yet since it is only a transform). These coefficients can then be compressed more easily because the information is statistically concentrated in just a few coefficients. This principle is called transform coding. After that, the coefficients are quantized and the quantized values are entropy encoded and/or run length encoded.

A few 1D and 2D applications of wavelet compression use a technique called "wavelet footprints".

Evaluation

Requirement for image compression

For most natural images, the spectrum density of lower frequency is higher. As a result, information of the low frequency signal (reference signal) is generally preserved, while the information in the detail signal is discarded. From the perspective of image compression and reconstruction, a wavelet should meet the following criteria while performing image compression:

  • Being able to transform more original image into the reference signal.
  • Highest fidelity reconstruction based on the reference signal.
  • Should not lead to artifacts in the image reconstructed from the reference signal alone.

Requirement for shift variance and ringing behavior

Wavelet image compression system involves filters and decimation, so it can be described as a linear shift-variant system. A typical wavelet transformation diagram is displayed below:

The transformation system contains two analysis filters (a low pass filter

Biorthogonal filters for wavelet image compression
Length Filter coefficients Regularity
Wavelet filter 1 H0 9 .852699, .377402, -.110624, -.023849, .037828 1.068
G0 7 .788486, .418092, -.040689, -.064539 1.701
Wavelet filter 2 H0 6 .788486, .047699, -.129078 0.701
G0 10 .615051, .133389, -.067237, .006989, .018914 2.068

By observing the impulse responses of the two filters, we can conclude that the second filter is less sensitive to the input location (i.e. it is less shift variant).

Another important issue for image compression and reconstruction is the system's oscillatory behavior, which might lead to severe undesired artifacts in the reconstructed image. To achieve this, the wavelet filters should have a large peak to sidelobe ratio.

So far we have discussed about one-dimension transformation of the image compression system. This issue can be extended to two dimension, while a more general term - shiftable multiscale transforms - is proposed.

Derivation of impulse response

As mentioned earlier, impulse response can be used to evaluate the image compression/reconstruction system.

For the input sequence

On the other hand, to reconstruct the signal x(n), we can consider a reference signal

To obtain the overall L level analysis/synthesis system, the analysis and synthesis responses are combined as below:

Finally, the peak to first sidelobe ratio and the average second sidelobe of the overall impulse response

Using a wavelet transform, the wavelet compression methods are adequate for representing transients, such as percussion sounds in audio, or high-frequency components in two-dimensional images, for example an image of stars on a night sky. This means that the transient elements of a data signal can be represented by a smaller amount of information than would be the case if some other transform, such as the more widespread discrete cosine transform, had been used.

Limitations

While wavelet transforms offer theoretical advantages, their practical limitations have effectively limited wavelet compression to analyzing localized changes and transient signals. Despite decades of research, wavelet-based compression systems for common multimedia like audio and video do not consistently match the efficiency and perceptual quality of current Discrete Cosine Transform-based systems.

For one-dimensional data like audio or ECGs, wavelets excel at representing and compressing transient signals—sudden, isolated events such as a drum hit in music or the sharp peaks in a heart rhythm. For example, the discrete wavelet transform has been successfully applied for the compression of electrocardiograph (ECG) signals. However, for smooth, periodic signals, which make up much of typical audio, harmonic analysis in the frequency domain with Fourier-related transforms achieve better compression and sound quality. Compressing data that has both transient and periodic characteristics may be done with hybrid techniques that use wavelets along with traditional harmonic analysis. For example, the Vorbisaudio codec primarily uses the modified discrete cosine transform to compress audio (which is generally smooth and periodic), however allows the addition of a hybrid wavelet filter bank for improved reproduction of transients.

For higher-dimensional data, wavelet compression faces significant challenges. In video, for instance, modern compression techniques such as intra coding and motion compensation (predicting parts of an image based on what's next to it spatially and temporally) and mixed and dynamic block sizes become incredibly complex with wavelets because of their overlapping nature. This complexity translates to more processing power and slower speed, making them less practical for widespread use. Furthermore, while wavelets might score well on traditional measures such as PSNR, DCT blocks create a perception of sharpness that wavelets often lack, requiring higher bitrates to achieve similar subjective quality.

Comparison with Fourier transform and time-frequency analysis

TransformRepresentationInput
Fourier transform
Time–frequency analysis
Wavelet transform

Wavelets have some slight benefits over Fourier transforms in reducing computations when examining specific frequencies. However, they are rarely more sensitive, and indeed, the common Morlet wavelet is mathematically identical to a short-time Fourier transform using a Gaussian window function. The exception is when searching for signals of a known, non-sinusoidal shape (e.g., heartbeats); in that case, using matched wavelets can outperform standard STFT/Morlet analyses.

Other practical applications

The wavelet transform can provide us with the frequency of the signals and the time associated to those frequencies, making it very convenient for its application in numerous fields. For instance, signal processing of accelerations for gait analysis, for fault detection, for the analysis of seasonal displacements of landslides, for design of low power pacemakers and also in ultra-wideband (UWB) wireless communications.

  1. Discretizing of the

    Applied the following discretization of frequency and time:

    Leading to wavelets of the form, the discrete formula for the basis wavelet:

    Such discrete wavelets can be used for the transformation:

  2. Implementation via the FFT (fast Fourier transform)

    As apparent from wavelet-transformation representation (shown below)

    where

    and as already mentioned in this context, the wavelet-transformation corresponds to a convolution of a function

    • Fourier-transformation of signal
    • Selection of a discrete scaling factor
    • Scaling of the wavelet-basis-function by this factor
    • Multiplication with the transformed signal YFFT of the first step
    • Inverse transformation of the product into the time domain results in
    • Back to the second step, until all discrete scaling values for
    There are many different types of wavelet transforms for specific purposes. See also a full list of wavelet-related transforms but the common ones are listed below: Mexican hat wavelet, Haar Wavelet, Daubechies wavelet, triangular wavelet.
  3. Fault detection in electrical power systems.
  4. Locally adaptive statistical estimation of functions whose smoothness varies substantially over the domain, or more specifically, estimation of functions that are sparse in the wavelet domain.

Time-causal wavelets

For processing temporal signals in real time, it is essential that the wavelet filters do not access signal values from the future as well as that minimal temporal latencies can be obtained. Time-causal wavelets representations have been developed by Szu et al and Lindeberg, with the latter method also involving a memory-efficient time-recursive implementation.

Synchro-squeezed transform

Synchro-squeezed transform can significantly enhance temporal and frequency resolution of time-frequency representation obtained using conventional wavelet transform.