Pour préciser le phénomène, nous chercherons maintenant à travailler sur une image synthétique contenant deux ondes concentriques de basse et haute fréquence (0.02 cycles/pixel et 0.2 cycle/pixel) pour cela nous créerons une image de taille 384*384 définie par :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
I(i,j)=0.5+0.25*cos(2*pi*f1*sqrt((i-193)^2+(j-193)^2)+0.25*cos(2*pi*f2*sqrt((i-193)^2)+(-193)^2))
alors, grossièrement j'ai crée le script foireux suivant :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
 
img =abs(I(384));
I(i,j)=(0.5+0.25*cos(2*pi*f1*sqrt((i-193)*(i-193)+(j-193)*(j-193))+0.25*cos(2*pi*f2*sqrt((i-193)*(i-193))+(-193)*(j-193))));
i=0.02;
j=0.2;
imshow (img)
mais j'ai l'impression de pas faire ce qu'il faut, forcément je me tape le code d'erreur
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
 
??? Index exceeds matrix dimensions.
 
Error in ==> C:\Documents and Settings\pat\Bureau\license\tdi\test.m
On line 1  ==> img =abs(I(384));
auriez vous une solution à m'apporter?