Discussion :

[ImaneN]

Bonsoir à tous,

Je voudrai implémenter sous matlab un filtre gaussien, cad une fonction qui admet comme argument une variance sigma et qui retourne une matrice carrée dont la taille est définie dans le corps de la fonction et dont les éléments sont des images de la fonction gaussienne de paramètre sigma en 2D. J'utiliserai alors cette fonction comme masque pour mettre du flou dans une image donnée.

Voici ce que j'ai écrit pour la fonction filtre:

Code :

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

1
2
3
4
5
6
7
function [f]=filtre_gaussien_2d(sigma)
N=1+2*ceil(3*sigma);
c=ceil(3*sigma)+1; %point central de la matrice (de taille impaire)
xint=(1-c):1: (N-c); 
x=exp(-(xint.^2)/(2*(sigma^2)));
f=(1/(2*pi*sigma^2))*(x*xt);
end

Malheureusement, lorsque je teste cette fonction en exécutant la commande :

Code :

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

filtre_gaussien_2d(0.5)

par exemple, on me retourne le message d'erreur suivant

??? Error: File: filtre_gaussien_2d.m Line: 19 Column: 5
The expression to the left of the equals sign is not a valid target for an assignment.

Je ne comprend pas pourquoi on me retounr 'Column 5'. A quoi ça correspond ? Et puis la ligne 19 de mon script filtre_gaussien_2d.m correspond à 'end'. Peut être on parle du signe égal qui se trouve juste à la ligne du dessus ?

Je n'arrive pas à voir où se trouve l'erreur. Si vous pouvez m'aider !

Bien à vous.

[Jerome Briot]

Le code que tu nous montre ne correspond pas au message d'erreur.

Le nom de la fonction n'est pas le même dans l'appel et dans le message d'erreur (filtre_gaussien_2d et gaussian_filter_2d).
Il n'y a pas 19 lignes dans le code que tu nous montre.

Pour le "Column: 5", je te renvoi vers la FAQ : Comment exploiter les messages d'erreurs ?

[ImaneN]

Très bonne remarque ! Je viens d'éditer mon post en conséquence.

Dans mon code original il y a effectivement 19 lignes car j'ai rajouté des commentaires, et la ligne 19 correspond au 'end' de la fonction.

[ImaneN]

Finalement, c'était la transposé du vecteur ligne x qui posait problème. Je pensais que la syntaxe était xt alors que c'est x'.
Novice en Matlab !

Bon, enfin je teste ma fonction pour différents sigma et je trouve qu'elle me retourne les matrices attendues
(sauf pour le cas sigma=0.5 où elle me retourne une matrice de taille 5 alors qu'en calculant N avec sigma=0.5, on obtient N=3, je n'ai pas encore compris pourquoi).

Je prends une image et je la convolue avec ma fonction filtre en faisant varier sigma. Je constate que plus sigma est grand plus les bords de l'image convoluée sont noirs.

Je me demande pourquoi et tout d'un coup, ça me fait penser à la manière dont on fait cette opération de convolution : on prend chaque pixel de l'image d'origine et on considère la sous-matrice centrée sur ce pixel et qui est de même taille que le filtre et on fait le produit de convolution bien connu des matrices carrées de même taille. Ce produit de convolution donne un réel. On met alors à jour la valeur du pixel considéré avec le réel obtenu.

Et ben cela soulève un petit problème : c'est quoi la sous matrice centrée sur un pixel situé au coin de l'image ? Par exemple pour le pixel au coin haut gauche, il n'y a pas de pixels de l'image d'origine en dessus ni gauche !

Comment on fait alors pour ces pixels ? Est ce que ça a une relation avec le fait que les bords de l'image convoluée deviennent de plus en plus sombre en agrandissant sigma ?