Taille des coefficients de la décomposition en ondelette dwt2

Bonjour,
Petite question sur la taille des coefficients que la fonctions dwt2 de la Wavelet tlbx produit.

Voilà un bout de code où on charge une image de taille 256x256 et où on calcule l'approx et les détails de niveau 1 pour deux ondelettes.
Code:

1 2 3 4 5 6 7 8 clear all; % Chargement de l'image 'barbara', taille 256x256 load wbarb; clear map; % Single-level discrete 2-D wavelet transform [CA1,CH1,CV1,CD1] = dwt2(X,'haar'); whos [CA1,CH1,CV1,CD1] = dwt2(X,'db2'); whos
dont voici le résultat

Citation:

Name Size Bytes Class Attributes

CA1 128x128 131072 double
CD1 128x128 131072 double
CH1 128x128 131072 double
CV1 128x128 131072 double
X 256x256 524288 double

Name Size Bytes Class Attributes

CA1 129x129 133128 double
CD1 129x129 133128 double
CH1 129x129 133128 double
CV1 129x129 133128 double
X 256x256 524288 double

Pourquoi une taille de 129x129 des coefficients pour l'ondelette daubechie 2 ?

Théoriquement cela devrait être comme celle de haar, de taille 128x128.

De plus, en réalisant exactement le même type de décomposition mais cette fois-ci en utilisant l'outil graphique de matlab (wavemenu), on retrouve la bonne dimension :

http://img534.imageshack.us/img534/2126/wave1.png

Après avoir cliquer sur "Full size" :
http://img203.imageshack.us/img203/3092/wave2.png

Rien à ce sujet dans la doc de matlab, ni même sur google...

Bien sûr même problème avec wavedec2 (décomposition sur plusieurs niveaux) qui n'est qu'une boucle sur dwt2...

Bonjour,

Dans la doc de dtw2 :

Citation:

Let sx = size(X) and lf = the length of filters; then size(cA) = size(cH) = size(cV) = size(cD) = sa where sa = ceil(sx/2), if the DWT extension mode is set to symmetric padding. For the other extension modes, sa = floor((sx+lf-1)/2).

Donc vraisemblablement, avec les résultats de ton code, dans un cas c'est la formule ceil(sx/2), et dans l'autre floor((sx+lf-1)/2) pour définir la taille de ton résultat. Ce que je ne comprends pas, c'est que tu n'as pas modifié le DWT extension mode extension mode dans ton cas, donc pour moi dans les 2 cas, il aurait dû utiliser celui par défaut, mais ça dépend peut-être du filtre utilisé.

Dans cette discussion, dans la réponse on trouve :

Code:

dwtmode('per','nodisp');

Cela permet de préciser la façon dont doit être traité les distorsions sur les bords et de forcer à avoir une taille d'image égale à la moitié de la taille de l'image initiale. Donc tu peux essayer ce réglage avant d'appliquer tes filttres. Après je te conseille de regarder les réglages de dwtmode pour plus de renseignements :

Code:

help dwtmode

Bonjour,
Merci pour la réponse, ça marche en effet avec dwtmode.

Le mode par défaut chez moi est :

*******************************************************
** DWT Extension Mode: Symmetrization (half-point) **
*******************************************************

ce qui signifie que c'est floor((sx+lf-1)/2) qui s'applique et donc bien 128 pour haar et 129 pour db2 car les tailles (lf) des ondelettes varient.
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 sx = size(X) W = dbwavf('db1'); lf = length(W) floor((sx+lf-1)/2) sx = size(X) W = dbwavf('db2'); lf = length(W) floor((sx+lf-1)/2) sx = 256 256 lf = 2 ans = 128 128 sx = 256 256 lf = 4 ans = 129 129
Donc a priori dans wavemenu ils utilisent directement le mode 'periodization'...