Ondelette Morlet

**meera** · 03/08/2006, 23h20

Salut
Je veux appliquer la transformation ondelettes sur une image.
J'ai choisi l'ondelette de Morlet et j'ai écrit ce code :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
36
37
38
39
40
41
42
43
44
45
figure;
y=imread ('lena.gif');      % image originale
imshow(y);
title('Image originale');
[m,n]=size(y);
 
 
clf
ii=sqrt(-1);
%n=length(y)-1;
 
 
pi=3.14;
p=fix(log(n)/log(2));
dx=1/n;
t=0:dx:1;
x=-1:dx:1;
%x=-4:dx:4;
%
% calcul de l'ondelette
  a0=0.5;
  %a0=1/2;
a0=a0^(1/2);
p=2*p;
W=zeros(p,n+1);
%for i=p-1:-1:0
for e=1:m 
for i=0:p-1   
 
%	disp('echelle')
	a=a0^i;
	xx=x/a; 
	%g=2*pi*((1-2*pi*xx.^2).*exp(-pi*xx.^2))./sqrt(a);
	%plot(x,g),title('ondelette chapeau mexicain')
                 g=exp(-xx.^2/2).*exp(ii*5*xx)./sqrt(a);
	    plot(x,real(g)),title('ondelette de Morlet-partie reelle ')
%   disp('appuyer sur une touche pour continuer')
%   pause
%
% calcul des coefficients d'ondelettes a l'echelle a
%
	wa=conv(y(e,:),g);
	W(i+1,1:n)=abs(wa(n+1:2*n));
end
end

mais, quand j'exécute, j'obtiens cette erreur :

??? Error using ==> filter
Function 'filter' not defined for variables of class 'uint8'.

Error in ==> C:\MATLAB6p1\toolbox\matlab\datafun\conv.m
On line 37 ==> c = filter(a, 1, b);

Error in ==> C:\MATLAB6p1\work\moretfin.m
On line 52 ==> wa=conv(y(e,

,g);

Est-ce que quelqu'un peut m'aider pour résoudre ce problème

**Kangourou** · 04/08/2006, 10h22

salut,

essaie de convertir tes donnees 'uint8' en donnees 'double', comme ca :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

res = conv(double(signal), f);

A+

**meera** · 04/08/2006, 14h11

j'ai essai de modifier :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

wa=conv(y(e,:),g);

par

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

wa=conv(double(y),g);

çàd j'ai supprimer la boucle de e

mais quand j'execute j'oubtient:

??? Error using ==> conv
A and B must be vectors.

Error in ==> C:\MATLAB6p1\work\moretfin.m
On line 53 ==> wa=conv(double(y),g);

d'autre part si je laisse la boucle et j'ajout Double :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

wa=conv(double(y(e,:)),g);

le matlab se bloque !!!!!!!!!!

**meera** · 04/08/2006, 21h17

j'ai modifier le ligne par celle çi:
code:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
 
 wa=conv(double(y(:)),g);
W(i+1,1:n)=abs(wa(n+1:2*n));

mais loes de l'execution il s'affiche:
citation :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
 
???  In an assignment  A(matrix,matrix) = B, the number of rows in B
and the number of elements in the A row index matrix must be the same.
 
Error in ==> C:\MATLAB6p1\work\moretfin.m
On line 55  ==> 	W(i+1,1:n)=abs(wa(n+1:2*n));

est quelqu'un sait commet je resolu ce pb??

**Matthieu Brucher** · 04/08/2006, 21h30

Normal, W(i+1;1:n) est du 1*n quand abs(wa(n+1:2*n)) est de taille n*1

**meera** · 04/08/2006, 21h59

oui oui
mais quand je convrtis les lignes par collonne
matlab se bloque

**Matthieu Brucher** · 04/08/2006, 23h46

Pas normal, tu l'as réinstallé ? Et c'est le cas sur plusieurs machines ?

**meera** · 05/08/2006, 13h32

non je ne sais pas , je l'essai sur ma machine seulement!!!
maintenant j'ai modifier une portion de mon code ,
précedament j'avais:

code:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
 
wa=conv(double(y(:)),g);
W(i+1,1:n)=abs(wa(n+1:2*n));

donc wa c'est une vecterur
moi j'aimerais qu'elle deviendra matrice pour pouvoir l'ajouter la marque apres
je l'ai modifié par:

code:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
 
wa(:)=conv(double(y(:)),g);
W(i+1,1:n)=abs(wa(n+1:2*n));

mais cet erreur s'affichera:

citation:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
 
???  In an assignment  A(:) = B, the number of elements in A and B
must be the same.
 
Error in ==> C:\MATLAB6p1\work\chapeaaaaaaaaaa.m
On line 60  ==> wa(:)=conv(double(y(:)),g);

je sais pas commet je resoudre!!!!!!!!!

**Matthieu Brucher** · 05/08/2006, 14h09

quelle est la taille de y, de g et de wa ?

**meera** · 05/08/2006, 14h18

la taille de y : 256*256
..............g : 1 * 513
pou la taille de wa qaund elle est vecteur est 66048 * 1
mais apres modification de code par :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
 
wa(e,:)=conv(y(e,:),g);
 W(i+1,1:n)=abs(wa(n+1:2*n));

la taille deviet wa : 256*256

et cet erreur s'affiche
citation:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
 
??? Error using ==> filter
Function 'filter' not defined for variables of class 'uint8'.
 
Error in ==> C:\MATLAB6p1\work\conv.m
On line 37  ==>     c = filter(a, 1, b);
 
Error in ==> C:\MATLAB6p1\work\chapeaaaaaaaaaa.m
On line 54  ==> wa(e,:)=conv(y(e,:),g);
 
quant j'ajout double devant y matlab se bloque:?

**Matthieu Brucher** · 05/08/2006, 14h23

L'erreur est là :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

wa(e,:)=conv(y(e,:),g);

Pas autrepart, donc vérifie les tailles. Je ne comprend pas ce que tu dis, une fois ça vaut 66048 et après, ça vaut 256*256, ce qui n'est pas possible.

**meera** · 05/08/2006, 15h02

code:
precedament j'ai tapé ce code :
Code :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
 
wa=conv(double(y(:)),g);
W(i+1,1:n)=abs(wa(n+1:2*n));

donc wa c'est une vecterur de taille 66048*1
moi j'aimerais que wa sooit une matrice
pour cela j'ai appliqué la convolution sur chaque ligne
comme vous voyez sur le code principale j'ai utilisé une boucle
(for for e=1:m )

Code :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
 
  wa(e,:)=conv(y(e,:),g);
   W(i+1,1:n)=abs(wa(n+1:2*n));

cette fois çi la taille de wa deviendra 256*256
et ça ce que je veux
car apres l'application de la transformation ondeletta j'ajoute wa + z avec z est un marque( matrice 256*256 rempli pseudo aleatoirement)

le pb c que matlab n'adet pas cette modification , il continu de me donner la meme erreur!!!!!!!

**Matthieu Brucher** · 05/08/2006, 15h05

wa est de taille 256*256, mais conv(y(e,

,g) est de taille 1*768, c'est pas trop possible de faire ce que tu fais.
Encore une fois, ça ne sert à rien de poster la deuxième ligne, elle nous embrouille et n'est pas responsable de ton erreur.

**meera** · 05/08/2006, 15h20

Ah oui vous avez raison
mais puisque wa: 256*256 et g :1 *768
Pourquoi je n'applique pas la convolution pour chaque ligne de wa
cette fois-ci j'aurais :
wa : 1*256
g:1*768
est-ce possible !

**Matthieu Brucher** · 05/08/2006, 15h23

En français, ça donne quoi, je n'ai pas compris ta question.

**meera** · 05/08/2006, 15h37

Ma question c'était :
si j'applique la convolution de y et g sur chaque ligne de wa
c'est à dire je considère que wa est une matrice de 256*256

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

wa=zeros(m,n);

pour pouvoir résoudre le problème de dimension
cette fois-ci pour chaque ligne (selon la boucle de ligne) j'aurais:
wa : 1*256
g : 1*768

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
for e=1:m
    wa(e,:)=conv(double(y(:)),g);
end

mais j'aurais cet erreur !

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
 
???  In an assignment  A(I) = B, the number of elements in B and
 I must be the same.
 
Error in ==> C:\MATLAB6p1\work\chapeaaaaaaaaaa.m
On line 55  ==>    wa(e)=conv(double(y(:)),g);

**Matthieu Brucher** · 05/08/2006, 15h39

Il te faut couper les éléments extérieurs, ne pas les utiliser. Tu stockes dans une variable temporaire le retour de la convolution et tu mets dans wa la partie qui t'intéresse.

**meera** · 05/08/2006, 15h52

j'ai pas bien compris
vous pouvez me montrer par un exemple

**Matthieu Brucher** · 05/08/2006, 16h00

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
p = conv(...);
wa(ceque tu veux) = p(la partie qui t'intéresse);

**meera** · 05/08/2006, 16h05

Ah ok.
Je vais essayer, et voir ce que cela me donne maintenant.
Merci

J'ai mis ce code:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
wa=zeros(m,n);
jj=conv(double(y(:)),g);
wa(:,e)=jj(1:e);
W(i+1,1:n)=abs(jj(n+1:2*n));     %  quand je change jj par wa MATLAB se bloque

j'obtiens cette erreur:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
 
??  In an assignment  A(matrix,matrix) = B, the number of rows in B
and the number of elements in the A row index matrix must be the same.
 
Error in ==> C:\MATLAB6p1\work\chapeaaaaaaaaaa.m
On line 58  ==>   W(i+1,1:n)=abs(jj(n+1:2*n));