In functional analysis, compactly supported wavelets derived from Legendre polynomials are termed Legendre wavelets or spherical harmonic wavelets.[1] Legendre functions have widespread applications in which spherical coordinate system is appropriate.[2] [3] [4] As with many wavelets there is no nice analytical formula for describing these harmonic spherical wavelets. The low-pass filter associated to Legendre multiresolution analysis is a finite impulse response (FIR) filter.
Wavelets associated to FIR filters are commonly preferred in most applications.[3] An extra appealing feature is that the Legendre filters are linear phase FIR (i.e. multiresolution analysis associated with linear phase filters). These wavelets have been implemented on MATLAB (wavelet toolbox). Although being compactly supported wavelet, legdN are not orthogonal (but for N = 1).[5]
Associated Legendre polynomials are the colatitudinal part of the spherical harmonics which are common to all separations of Laplace's equation in spherical polar coordinates.[2] The radial part of the solution varies from one potential to another, but the harmonics are always the same and are a consequence of spherical symmetry. Spherical harmonics
Pn(z)
2nd
\left(1-z2\right)
| d2y | |
| dz2 |
-2z
| dy | |
| dz |
+n(n+1)y=0.
Pn(\cos(\theta))
H(\omega)
|H(0)|=1
|H(\pi)|=0
|H(\omega)|
\nu=2n+1.
|H\nu(\omega)|=\left|
| ||||||||||
| P\nu\cos(0) |
\right|
Illustrative examples of filter transfer functions for a Legendre MRA are shown in figure 1, for
\nu=1,3,5.
-\pi<\omega<\pi
\nu
The low-pass filter transfer function is given by
H\nu
| ||||||
| (\omega)=-e |
P\nu\left(\cos\left(\tfrac{\omega}{2}\right)\right)
The transfer function of the high-pass analysing filter
G\nu(\omega)
H\nu(\omega)=-e-j
| \omega | |
| 2 |
Indeed,
|G\nu(0)|=0
|G\nu(\pi)|=1
A suitable phase assignment is done so as to properly adjust the transfer function
H\nu(\omega)
H\nu(\omega)=
| 1 | |
| \sqrt{2 |
The filter coefficients
\{hk\}k
| \nu | |
| h | |
| k |
=-
| \sqrt{2 | |
from which the symmetry:
| \nu | |
| {h | |
| k |
follows. There are just
\nu+1
Hn(\omega)
\nu
Table I - Smoothing Legendre FIR filter coefficients for
\nu=1,3,5
N
\nu=1(N=1) | \nu=3(N=2) | \nu=5(N=3) | |
h0 | -\tfrac{\sqrt{2}}{2} | -5\tfrac{\sqrt{2}}{16} | -63\tfrac{\sqrt{2}}{256} |
h1 | -\tfrac{\sqrt{2}}{2} | -3\tfrac{\sqrt{2}}{16} | -35\tfrac{\sqrt{2}}{256} |
h2 | -3\tfrac{\sqrt{2}}{16} | -30\tfrac{\sqrt{2}}{256} | |
h3 | -5\tfrac{\sqrt{2}}{16} | -30\tfrac{\sqrt{2}}{256} | |
h4 | -35\tfrac{\sqrt{2}}{256} | ||
h5 | -63\tfrac{\sqrt{2}}{256} | ||
N.B. The minus signal can be suppressed.
Legendre wavelets can be easily loaded into the MATLAB wavelet toolbox—The m-files to allow the computation of Legendre wavelet transform, details and filter are (freeware) available. The finite support width Legendre family is denoted by legd (short name). Wavelets: 'legdN'. The parameter N in the legdN family is found according to
2N=\nu+1
Legendre wavelets can be derived from the low-pass reconstruction filter by an iterative procedure (the cascade algorithm). The wavelet has compact support and finite impulse response AMR filters (FIR) are used (table 1). The first wavelet of the Legendre's family is exactly the well-known Haar wavelet. Figure 2 shows an emerging pattern that progressively looks like the wavelet's shape.
The Legendre wavelet shape can be visualised using the wavemenu command of MATLAB. Figure 3 shows legd8 wavelet displayed using MATLAB. Legendre Polynomials are also associated with windows families.[8]
Wavelet packets (WP) systems derived from Legendre wavelets can also be easily accomplished. Figure 5 illustrates the WP functions derived from legd2.