commande double et choix du seuil pour la détection de contour

**helaaa** · 15/07/2009, 18h30

bonjour,

j'ai quelque lignes de codes que je n'arrive pas à comprendre

Code :

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

1
2
3
4
5
6
7
8
9
10
I1=imread([FileName]);
[l c]=size(I1);
I3=double(I1);
I2=double(edge(I1,'canny',0.1,2));
I4=I3+I2*255;
for i=1:c;
    ind1=find(double(I4(:,i))>144,1); %  chercher le premier pixel qui dépasse la valeur seuil
    if ~isempty(ind1); indz1(i)=ind1; 
    end;  
end;

ce programme de contour doit normalement détecter la surface d'une image qui en haut est noir et en bas pleine d 'eau donc normalement ça aurait suffit de faire find(double(I2(:,i))==1,1); et donc de tracer la sortie du premier contour détecté par edge, or ce que je comprend pas c'est la sortie de double sur laquelle il définit un seuil ( d'ailleurs je sais pas sur quel critère il le choisit) et prends la première valeurs dépassant ce seuil

cette méthode qu'il utilise (ajouter la recherche d'un autre seuil a celle de la sortie de edge) marche très bien sur mes images, mieux même que celle ou je prends les premières valeurs de 1, mais je comprends pas pourquoi?

merci d'avance

**Jerome Briot** · 15/07/2009, 21h28

Il faut repartir de l'image initiale.

A priori les données sont stockées en uint8 donc la plage de valeur est [0 255]

Comme il est parfois difficile de faire des opérations entre les variables de cette classe sous MATLAB, il est nécessaire de les transformer en Double (leur valeur ne change pas lors de cette conversion)

EDGE renvoie un tableau logique composé de 0 (faux) et 1 (vrai). Comme pour l'image, il faut la transformer en Double pour pouvoir faire de opération avec.

Donc la ligne suivante sert à renforcer les contours détectés par EDGE en mettant la valeur des pixels correspondant à 255 (valeur maximale en uint8) :

Code :

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

I4=I3+I2*255;

Ensuite pour la valeur du seuil utilisée, elle dépend certainement de l'image de départ, non ?

A noter que pour la relation logique, il n'est pas nécessaire de passer une nouvelle fois en classe Double :

Code :

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

ind1=find(I4(:,i)>144,1)

**helaaa** · 15/07/2009, 21h52

merci pour la réponse

mais moi en fait je voulais savoir pourquoi il fait cette étape là pour la détection de contour

Code :

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

ind1=find(double(I4(:,i))>144,1); %  chercher le premier pixel qui dépasse la valeur seuil

normalement il fait ceci

Code :

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

ind1=find(double(I2(:,i))==1,1);

le test sur I4 va donner quoi? et le choix du seuil n'est pas bizarre vu que les contours dans I4 sont multipliés par 255 alors qu'il ne met que 144 comme valeur???

**Jerome Briot** · 15/07/2009, 21h58

Envoyé par helaaa

le choix du seuil n'est pas bizarre vu que les contours dans I4 sont multipliés par 255 alors qu'il ne met que 144 comme valeur???

Oui mais les valeurs de I3 qui ne correspondent pas à un contour issus de EDGE sont conservées dans I4

Donc dans I3 il y a des valeurs comprises entre 0 et 255, non ?

**Jerome Briot** · 15/07/2009, 22h02

Envoyé par Dut

Donc la ligne suivante sert à renforcer les contours détectés par EDGE en mettant la valeur des pixels correspondant à 255 (valeur maximale en uint8) :

Code :

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

I4=I3+I2*255;

Je crois que j'ai dit une bêtise la... théoriquement, certaines valeurs de I4 peuvent être supérieures à 255... car I2, I3 et I4 sont de classe Double

**helaaa** · 15/07/2009, 22h41

donc le I4 c'est l'image initiale à laquelle on ajoute le contour de edge multiplié par 255
I4 est déjà en double, pourquoi alors il met cette condition encore une 2ème fois en double? cela peut justifier sa valeur de seuil??