Opération d'érosion

Version imprimable

Opération d'érosion

Bonsoir,
Voilà, alors je suis débutant en matlab et là j'ai un projet en traitement d'images où je dois effectuer divers effets comme l'érosion par exemple.

Je voudrais utiliser l'algorithme suivant :

Citation:

Pour tous les pixels (x, y) de image_1 Faire :
Début
Tab1(x, y) = Faux
Tab2(x, y) = Faux
Si image_1 >0 Alors Tab1(x, y) = VRAI
Fin

Pour tous les pixels (x, y) Faire
Début
Si Tab1(x-1, y-1) est VRAI
ET Tab1(x-1, y)
ET Tab1(x-1, y+1)
ET Tab1(x+1, y-1)
ET Tab1(x+1, y)
ET Tab1(x+1, y+1)
Alors Tab2(x, y) = Tab1(x, y)
Sinon Tab2(x, y) = FAUX
Fin

Pour tous les pixels (x, y) Faire
Début
Si Tab2(x, y) est VRAI alors image_1 = 255
Sinon image_1 = 0
Fin

Seulement voilà, je n'arrive pas à le traduire en matlab, je patauge complètement. Donc si vous pouviez m'y aider, ça serait sympa...

Après de nombreuses tentatives qui ne donnent strictement rien, je vais mettre le code le moins pourri que j'ai pondu... Je suis sûr de me planter dans l'utilisation des boucles for.
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 for I=1:i for J=1:j t1=0; t2=0; if IBinaire>0 t1=1; end if t1=1 & t1(i-1,j) & t1(i-1,j+1) & t1(i+1,j-1) & t1(i+1,j) & t1(i+1,j+1) t2=t1; else t2=0; end if t2=1 IBinaire=255; else IBinaire=0; end end end
Voilà, ça vaut ce que ça vaut, l'image a été binarisée au préalable.
Merci d'avance

Bonne soirée

13/05/2008, 07h43
Jerome Briot
Pour commencer, tu te trompes dans l'utilisation des indices des boucles FOR-END

Tu fais :
Code:

1 2 3 4 5 6 for I=1:i for J=1:j t(i,j) end end
Alors qu'il faut faire :
Code:

1 2 3 4 5 for i=1:I for j=1:J t(i,j) end end
En ayant pris soin de définir I et J au préalable avant la boucle

Ensuite, l'opérateur logique d'égalité est == sous MATLAB.

Donc :
Code:

1 2 3 4 5 6 if t2==1 else end
13/05/2008, 09h48
mr_samurai

Salut,

Une fonction qui fait l'érosion est déjà présente dans l'Image Processing Toolbox : IMERODE

Code:

help imerode

++

Oui merci, mais j'essaye d'écrire une fonction qui fait le boulot de imerode !

Alors, j'ai réussi à résoudre pas mal d'erreurs mais maintenant, je ne sais pas comment récupérer l'image traitée:

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
Im=imread('\\iutb-cpt\profetu$\l.lorinquer\Bureau\Projet_info\image.jpg');
seuil=128;
IBinaire = Im > seuil;
image(IBinaire);
[t1,t2]=size(IBinaire);
for i=1:I
   for j=1:J
        t1(i,j)=false;
        t2(i,j)=false;
        if (IBinaire > 0)
            t1(i,j)=true;
        else
            t1(i,j)=false;
        end
        if ((t1==true) & (t1(i-1,j-1))  & (t1(i-1,j)) & (t1(i-1,j+1)) & (t1(i+1,j-1)) & (t1(i+1,j)) & (t1(i+1,j+1)))
            t2(i,j)=t1(i,j);
        else
            t2(i,j)=false;
        end
        if (t2(i,j)==true);
            IBinaire=255;
        else
            IBinaire=0;
        end
   end
end

Je n'arrive pas à récupérer et à afficher l'image ayant subie la dilatation + l'erosion.

Voilà, merci d'avance!

20/05/2008, 15h13
Jerome Briot

Tu peux nous expliquer ce que représentent les variables "t1", "t2", "I" et "J" dans la ligne de code suivante ?

Code:

[t1,t2,I,J]=size(IBinaire);
20/05/2008, 15h46
yAwa_

Alors en fait, t1 et t2 sont deux tableaux de même taille que l'image binaire. L'idée de l'érosion c'est que l'on fait passer un élément structurant sur l'image à traiter.
Le centre de l’élément structurant se met sur le premier pixel à traiter, chaque partie de cet élément structurant devient alors momentanément égal à la valeur du pixel sur lequel est situé le pixel central. Après, le centre de l’élément structurant prend la plus petite valeur trouvée. Et ce principe se répète pour l’ensemble de l’image.

Pour I et J c'était une erreur de ma part, j'ai édité :aie:
20/05/2008, 16h07
Jerome Briot
Citation:

Envoyé par yAwa_

Alors en fait, t1 et t2 sont deux tableaux de même taille que l'image binaire.

Tiens c'est bizarre, ce n'est pas ce que dit la documentation :
Code:

1 2 [M,N] = SIZE(X) for matrix X, returns the number of rows and columns in X as separate output variables.

pour une erosion/dilatation tu te complique beaucoup la vie

tu a juste a faire un min/max sur ton element stucturant :

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
 
function [out]=Ero(in)
 
[Limg,Cimg,Bimg] = size(in);%recupere les infos de taille de l'image
out=zeros(Limg,Cimg);
 
 for i=2:Limg-1 %demarage a 2et fin a -1 pour eviter les erreurs due a l'effet de bord
    for j=2:Cimg-1
 
          temp(1)=in(i-1,j-1);
          temp(2)=in(i,j-1);
          temp(3)=in(i+1,j-1);
          temp(4)=in(i-1,j);
          temp(5)=in(i,j);
          temp(6)=in(i+1,j);
          temp(7)=in(i-1,j+1);
          temp(8)=in(i,j+1);
          temp(9)=in(i+1,j+1);
 
         %out(i,j)=min(min(temp)); %dilatation
 
         out(i,j)=max(max(temp));  %erosion
 
    end
end

voilou... c'est la version "je decompose tous" mais tu peu toujour modifier au besoin

31/05/2008, 18h09
yAwa_

Oui merci Amarth, effectivement ça marche! :D
Maintenant mon seul problème c'est que lorsque j'effectue ma dilatation ou mon erosion, l'effet est trop puissant. Disons en clair, que ça dilate ou erode trop!

Voilà, je voulais savoir si il y a un moyen de réduire cet effet, de fixer un "seuil" en quelque sorte.

Merci d'avance