Discussion :

[yAwa_]

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 :

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 :

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
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

[Jerome Briot]

Pour commencer, tu te trompes dans l'utilisation des indices des boucles FOR-END

Tu fais :

Code :

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

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 :

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

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 :

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

1
2
3
4
5
6
 
if  t2==1
 
else
 
end

[mr_samurai]

Salut,

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

Code :

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

help imerode

++

[yAwa_]

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 :

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
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!

[Jerome Briot]

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

Code :

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

[t1,t2,I,J]=size(IBinaire);

[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é

[Jerome Briot]

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 :

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

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

[Amarth]

pour une erosion/dilatation tu te complique beaucoup la vie

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

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
 
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

[yAwa_]

Oui merci Amarth, effectivement ça marche!
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