Schnelle Wavelet-Transformation

Die schnelle Wavelet-Transformation, englisch fast wavelet transform, ist ein effizientes Verfahren zur Berechnung einer diskreten Wavelet-Transformation. Sie kann mit der Anwendung der schnellen Fourier-Transformation zur Berechnung der Koeffizienten einer Fourier-Reihe verglichen werden.

Konstruktion

Analyse-Filterbank, $g=a^{-},h=b^{-}$

Rekursive Anwendung einer Analyse-Filterbank

Ein gegebenes kontinuierliches Signal $f$ wird zunächst durch orthogonale Projektion auf einen Unterraum $V_{-J}$ einer orthogonalen Multiskalenanalyse in eine zeitdiskrete Koeffizientenfolge $s^{(-J)}$ umgewandelt. Je größer $J$ ist, desto genauer ist die dadurch erzielte Approximation. In vielen Fällen ist es ausreichend,

$s_{n}^{(-J)}:=2^{-J/2}\,f(n/2^{J})$

zu setzen. Nun wird rekursiv aus jedem Tiefpasssignal $s^{(k)}$ ein neues Tiefpasssignal

$s^{(k+1)}={\frac {1}{\sqrt {2}}}(\downarrow 2)(a_{-}*s^{(k)})$

und das Hochpasssignal

$d^{(k+1)}={\frac {1}{\sqrt {2}}}(\downarrow 2)(b_{-}*s^{(k)})$

erzeugt. Zusammen bilden diese eine Analyse-Filterbank, die Operationen darin werden weiter unten erklärt.

Nach $M$ Schritten der Rekursion ergeben sich die Folgen

$d^{(-J+1)},\dots ,d^{(-J+M)}$   und   $s^{(-J+M)}$.

Das Ziel dieser Transformation ist, dass die $d^{(k)}$ „dünn“ besetzt sind und sich daher gut komprimieren lassen.

Sind die Filter $a$ und $b$ ausreichend frequenzselektiv, war das Ausgangssignal bandbeschränkt und wurde dem WKS-Abtasttheorem entsprechend die erste Koeffizientenfolge $s_{n}^{(-J)}$ gewonnen, so enthält das erste Tiefpassergebnis alle Signalbestandteile bis zur halben Nyquist-Frequenz, das Hochpassergebnis die darüberliegenden, beide Male mit einer der Bandbreite entsprechenden Abtastrate.

Analyse und Synthese

Der Fischgrätenzerlegung in der Multiskalenanalyse entspricht eine aus dem Tiefpass $a$ und dem Hochpass $b$ zusammengesetzte zeitdiskrete Filterbank, es wird ein zeitdiskretes Signal $x$ aufgeteilt in ein hohes Band $b^{-}*x$ und ein tiefes Band $a^{-}*x$ (Faltung von Folgen). Danach werden beide Signale heruntergetaktet (englisch downsampling) zu

$s=(\downarrow 2)(a^{-}*x)$   und   $d=(\downarrow 2)(b^{-}*x)$.

Mit $a^{-}$ sei dabei die zeitinvertierte Folge

$a^{-}=\{\dots ,a_{2},a_{1},a_{0},a_{-1},a_{-2},\dots \}$

bezeichnet. Das Heruntertakten einer Folge bedeutet, dass eine neue Folge aus den Gliedern mit geradem Index gebildet wird,

$(\downarrow 2)(y)=\{\dots ,y_{-4},y_{-2},y_{0},y_{2},y_{4},\dots \}$.

Alle diese Operationen zusammengefasst ergibt sich eine gliedweise Berechnungsvorschrift der Analyse-Filterbank

$s_{n}=\sum _{k}a_{k}\cdot x_{2n+k}$   und   $d_{n}=\sum _{k}b_{k}\cdot x_{2n+k}$.

Aus der Orthogonalität ergibt sich, dass das Ausgangssignal $x$ zurückgewonnen werden kann, zuerst werden die Tiefpass- und Hochpassanteile $s$ und $d$ in der Abtastrate hochgerechnet, dies wird als Upsampling bezeichnet, mit den Skalierungs- und Waveletmasken gefaltet und dann addiert,

$2x=a*(\uparrow 2)s+b*(\uparrow 2)d$

oder koeffizientenweise

$2x_{n}=\sum _{k}a_{n-2k}\cdot s_{k}+\sum _{k}b_{n-2k}\cdot d_{k}$.

Der Übergang von $x$ zu $(s,d)$ heißt Analyse, der inverse Synthese. Es ist ersichtlich, dass die Transformierte $(s,d)$ eines endlichen Signals nun etwa genauso viele Samples wie das Signal $x$ selbst hat, also genauso viel Information enthält.

Erweiterungen

Es ist nicht erforderlich, dass die Folgen in der Analyse-Filterbank mit denen in der Synthese-Filterbank wie oben übereinstimmen, nur ist dann nicht garantiert, dass die Kombination beider Filterbänke das Ausgangssignal rekonstruiert. Ist dies doch der Fall, spricht man von vollständiger Rekonstruktion (englisch perfect reconstruction) oder von Biorthogonalität der Wavelet-Basen.