Convolution d'une image et fonction

Version imprimable

11/10/2007, 18h20
samia_6

Convolution d'une image et fonction

bonjour,
je suis étudiante, je veux réaliser une application dans laquelle je calcule la convolution d'une image avec la fonction g et d'autre:
g(x)=1/(sqrt(pi*2)* exp(-x^2/2*sig^2), pour des valeurs différentes de sig
le problème c'est que je n'arrive pas à la faire j'ai essayé avec conv et conv2 mais hélas.
est ce que vous pouvez m'aider?
merci d'avance. sache que je suis préssé du temps.
11/10/2007, 18h29
Jerome Briot
Quel est ton problème ? La constructions du filtre gaussien ? Son application avec CONV2 ?

Si tu as accès à l'Image Processing Toolbox, regarde la fonction IMFILTER.

Si tu cherches des informations plus générales à propos des filtres :
- Introduction au traitement numérique d'image
- Les filtres usuels en traitement d'images
11/10/2007, 18h47
samia_6

désolée, mais coment faire la convolution d'une fonction avec une matrice.
11/10/2007, 18h57
samia_6

je vais m'expliquer:
je lis une image (A=imred) puis je veux convoluer A avec une fonction g(x), le premier problème c'est que g est en fonction d'un seul veriable, aussi g est menu d'un paramétre sigma , je doit calculer la convolution de A et g pour des différents valeur de sig et redissiner l'image.
merci d'avance

Exemple du filtre de flou uniforme (tiré du lien que j'ai donné ci-dessus) :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
figure
load mandrill
 
colormap gray
 
subplot(1,2,1)
imagesc(X)
axis image
title('Original')
 
% Filtre Flou uniforme
F=ones(3);
 
% Application du filtre
Xfiltr=conv2(X,F,'same');
 
subplot(1,2,2)
imagesc(Xfiltr)
axis image
title('Flou uniforme (noyau 3x3)')

11/10/2007, 19h09
samia_6

pardant mais je cherche pas le filtre, je cherche la convolution?
est ce que c'est la méme chose? et comment faire pour varier le paramétre sig?:oops:
11/10/2007, 19h19
Jerome Briot

Tu me sembles débuter en traitement d'images.
Il n'y a aucun problème à ça.

Par contre, tu devrais prendre le temps de lire les deux liens que je t'ai donné (ce n'est pas très long ;)).
11/10/2007, 19h21
samia_6

merci, je vais essayer de les implimenter. désolée pour le dérangement:D
12/10/2007, 02h54
paradize3

Hello,

ton filtre, m'a tout l'air d'etre un filtre gaussien. Pour l'appliquer a ton image tu as plusieurs possibilites plus ou moins efficaces (en terme de complexite). soit tu l'applique tel quel (x represente la distance = sqrt(dx^2+dy^2)), c'est l'approche generale que tu peux appliquer a n'importe quel filtre, mais aussi la moins efficace.

Soit tu le decomposes en deux (propriete du filtre Gaussien), t'appliques une fois suivant l'axe des x puis suivant l'axe des y. Simple est rapide.

Ou pour aller encore plus vite, tu peux encore utiliser une fft.

Salutations,

Gregoire
15/10/2007, 22h57
samia_6

et comment je le ferai pour x et y ?
15/10/2007, 23h23
samia_6

:cry:, car les resultat obtenu par un masque gaussien eyt la fonction imfilter ne me donne pas de bon résultat, la plupart des cas j'obtient une image blanche. sachant que j'ai essayé de l'appliquer sur une image binaire, est ce que c'est normal?
pour résoudre ce problème, que puis-je faire?
23/10/2007, 08h07
paradize3

Reprenons,

Le plus simple pour faire ta convolution c;est probablement d'utiliser la fonction conv2. En l'occurence tu veux faire une convolution de ton image (I) avec un filtre gaussien (appelons la matrice carree correspondante M).

Une premiere chose a decided c'est la taille de M (p.ex. 3x3, 5x5, ...), disons qu'elle est de taille NxN.

Tu peux remplir M simplement grace a ta fonction g(x, sigma). En l'occurence pour l'element M(i,j) tu peux donner une valeur du genre: g(sqrt((i-(N+1)/2)^2 + (j-(N+1)/2)^2), sigma).

Ensuite tu n'as plus qu'a effectuer la convolution.

Salutations,

Gregoire