|
Sommaire Généralités
|
|
Ondelettes continues | |
|
Résolution en temps et en fréquences | |
|
Inversion de la transformée en ondelettes | |
|
Transformations en ondelettes discrètes |
|
Principe de l'analyse espace - échelle ou analyse multirésolution | |
|
Définition |
|
Pourquoi des algorithmes pyramidaux ? | |
|
Décomposition | |
|
Application aux images | |
|
Restitution |
|
Théorie sur les ondelettes à support compact | |
|
Les ondelettes à support compact d’Ingrid Daubechies |
|
Les ondelettes isotropes | |
|
Les ondelettes anisotropes |
On veut trouver une base orthonormale d'ondelettes dans laquelle il sera
possible de décomposer le signal
appartenant à
(
est l'espace des réels). Pour cela, on peut utiliser l'analyse multirésolution
dans
.
Une telle analyse consiste à décomposer le signal sur une gamme très étendue d'échelles,opération que l'on peut comparer à
une cartographie. A chaque échelle, le signal sera remplacé par
l'approximation la plus adéquate que l'on puisse y tracer. En allant des échelles
les plus grossières vers les échelles les plus fines, on accède à des représentations
de plus en plus précises du signal donné. L'analyse s'effectue en calculant ce qui diffère
d'une échelle à l'autre, c'est-à-dire les détails à une résolution donnée.
Ceux-ci permettent, en corrigeant une approximation encore assez grossière,
d'accéder à une représentation d'une qualité
meilleure.
Les ondelettes correspondent à des degrés de résolutions: elles sont définies par la différence entre 2 fonctions échelles consécutives. Elles permettent de représenter les détails gagnés lors du passage d'une échelle à l'échelle plus petite suivante.
La méthode AMR donne toutes les bases d'ondelettes connues à ce jour. L'AMR de l'espace des fonctions de carrés sommables
consiste à le découper en une suite croissante de sous-espaces vectoriels fermés
. Ces sous-espaces ont un certain nombre de propriétés : ils forment une suite
emboîtée, leur intersection est réduite à
et leur réunion est dense dans
.
Chaque sous-espace est l'ensemble de toutes les approximations possibles d'un même
signal à l'échelle associée au sous-espace. Le signal à analyser sera
approximé par une succession de projections orthogonales sur les sous-espaces
.
|
Propriété |
Interprétation |
|
Vj+1 est l'image de Vj par une dilatation d'un facteur 2 |
Il existe une grille fréquentielle sous-jacente en progression géométrique. |
|
Pour tout j, Vj+1 est un sous-espace de Vj |
Un signal basse résolution est aussi un signal à haute résolution. |
|
Vj est invariant par translation de 2j |
Il existe une grille temporelle sous-jacente par pas de 2j. Il aurait d'ailleurs suffi de donner la propriété pour j=0. |
|
L'intersection des Vj est réduite à 0 dans L2. |
A résolution minimale, on perd toute l'image. |
|
La réunion des Vj est dense dans L2. |
A résolution infinie, on reproduit parfaitement tous les signaux. |
|
Il existe une fonction q telle que les translations entières de q forment une base de Riesz de V0. |
Chaque résolution est engendrée par une base d'atomes translatés de 2j. Une base de Riesz est une frame de vecteurs indépendants. |
Il faut à présent trouver une base engendrant
Soit la fonction
Soit
pour toute fonction
De plus, l'ensemble des fonctions
La projection de Le passage d'une résolution
.
appartenant à
telle que
,
(
)
soit une base orthonormée dans
.
La fonction
est appelée fonction d'échelle ou fonction d'interpolation.
le sous-espace vectoriel engendré par
,
(
).
est défini à partir de
par simple changement d'échelle :
de
.
est une base de
.
L'échelle associée
à la résolution
est
.
Le choix d'un facteur 2 correspond à une analyse dyadique. Si on pose :
,
alors la famille
engendrée par translation (paramètre
)
et par dilation (paramètre
)
de
est encore une base orthonormale de
.
sur
,
appelée
,
s'écrit sous la forme :
à une résolution plus fine
se fait en projetant la fonction
appartenant à
sur le sous-espace
.
Or
peut s'écrire sous la forme :
D'où, après un changement de variable adéquat :
Si ,
on a :
Pour son algorithme, Mallat considère une suite croissante de grilles emboîtées
qui va de la grille la plus fine à la grille la plus grossière. Le signal à
analyser est échantillonné sur la grille fine. Étant donné un signal
original discret
(avec
)
son approximation
(avec
)
est obtenue par convolution de
et du filtre de réponse impulsionnelle
précédemment décrit en gardant un échantillon sur deux à la sortie du
filtre. Cette procédure est itérée un certain nombre de fois.
Il est possible d'associer à h, précédemment décrit, un
filtre discret H dont il est la réponse impulsionnelle. Si la
transformation de Fourier discrète de ce filtre est notée ,
alors
doit vérifier les relations suivantes :
La transformée de Fourier de la fonction échelle
associée au filtre
est alors donnée par :
Le signal original est donc filtré par une suite de filtres passe-bas en cascade dont la fréquence de coupure décroît de moitié lorsque l'on passe de la résolution j+1à la résolution j.
où Wj est le sous-espace complémentaire de Vj dans Vj+1. est donc orthogonal à W, pour tout k différent de j.
Pour calculer la projection d'une fonction f(x) sur Wj, il faut définir une fonction y(x), appelée ondelette, telle qu'une base de West construite à partir de la translation et dilatation de y(x).
Exemple:

