Discussion :

[Jerome Briot]

Objectifs :

• ce bref article montre comment obtenir des contours colorés à partir d'un tableau logique (0 ou 1) comme celui retourné par la fonction EDGE.
• une méthode de superposition des contours en couleur sur une image en niveaux de gris est également présentée.

Public :

• la première partie s'adresse aux utilisateurs de l'Image Processing Toolbox (utilisation des fonction IMSHOW et IND2RGB)
• la seconde partie donne une alternative aux utilisateurs ne possédant pas cette toolbox.

N'hésitez pas à donner vos avis et à apporter des commentaires sur cet article à la suite de ce message

Vous pouvez également noter cet article en utilisant l'outil de notation de cette discussion dans la barre de menu en haut à droite

[Jerome Briot]

Prenons l'exemple issu de la documentation de la fonction EDGE

L'exemple montre comment obtenir les contours à l'aide des filtres de Prewitt et de Canny.

Code :

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

1
2
3
I = imread('circuit.tif');
BW1 = edge(I,'prewitt');
BW2 = edge(I,'canny');

Les trois images correspondante sont obtenues comme ceci :

Code :

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

1
2
3
4
5
6
7
8
figure
imshow(I)
 
figure
imshow(BW1)
 
figure
imshow(BW2)

[Jerome Briot]

Il peut être utile d'obtenir les images des contours issus des filtres en couleur

Pour ce faire, il suffit d'imposer une palette de 2 couleurs à la fonction IMSHOW.
La palette (Colormap en anglais) de M couleurs est un tableau Mx3 dont chaque ligne contient les valeurs R,G et B de chaque couleur.

Par exemple, pour obtenir les contours issus du filtre de Prewitt (tableau BW1) en rouge sur un fond noir :

Code :

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

1
2
3
4
5
figure
cmap = [0 0 0   % Couleur RGB pour les valeurs 0 de BW1 (fond)
        1 0 0]; % Couleur RGB pour les valeurs 1 de BW1 (contour)
 
imshow(BW1,'colormap',cmap)

Ou pour obtenir les contours issus du filtre de Canny (tableau BW2) en vert sur un fond magenta :

Code :

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

1
2
3
4
5
figure
cmap = [1 0 1   % Couleur RGB pour les valeurs 0 de BW2 (fond)
        0 1 0]; % Couleur RGB pour les valeurs 1 de BW2 (contour)
 
imshow(BW2,'colormap',cmap)

[Jerome Briot]

On peut aussi vouloir superposer les contours obtenus en couleur sur l'image initiale.
Dans ce cas, la technique précédente ne fonctionne plus. Il faut passer du tableau logique 2D à un tableau 3D type RGB (couleurs vraies)

Par exemple, pour superposer les contours issus du filtre de Prewitt (tableau BW1) en rouge sur l'image initiale (tableau I) :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
figure
 
imshow(I)
 
hold on
 
cmap = [0 0 0   % Couleur RGB pour les valeurs 0 de BW1 (fond)
        1 0 0]; % Couleur RGB pour les valeurs 1 de BW1 (contour)
 
BW1RGB = ind2rgb(BW1,cmap);
h = imshow(BW1RGB);
set(h,'alphadata',BW1)

ou pour superposer les contours issus du filtre de Canny (tableau BW2) en vert sur l'image initiale (tableau I) :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
figure
 
imshow(I)
 
hold on
 
cmap = [1 1 1   % Couleur RGB pour les valeurs 0 de BW2 (fond)
        0 1 0]; % Couleur RGB pour les valeurs 1 de BW2 (contour)
 
BW2RGB = ind2rgb(BW2,cmap);
h = imshow(BW2RGB);
set(h,'alphadata',BW2)

[Jerome Briot]

Pour les utilisateurs qui ne possèdent pas l'Image Processing Toolbox

Il est possible de colorer les contours, en supposant que vous avez un tableau de contour BW (0 et 1) correspondant à celui retourné par EDGE, comme ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
 
figure
 
cmap = [0 0 0   % Couleur RGB pour les valeurs 0 de BW (fond)
        1 0 0]; % Couleur RGB pour les valeurs 1 de BW (contour)
 
colormap(cmap);
 
imagesc(BW)
axis image

ou comme ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
figure
 
cmap = [0 0 0   % Couleur RGB pour les valeurs 0 de BW1 (fond)
        1 0 0]; % Couleur RGB pour les valeurs 1 de BW1 (contour)
 
RGB(:,:,1) = BW*cmap(2,1)+~BW*cmap(1,1);
RGB(:,:,2) = BW*cmap(2,2)+~BW*cmap(1,2);
RGB(:,:,3) = BW*cmap(2,3)+~BW*cmap(1,3);
 
image(RGB)
axis image

[Jerome Briot]

Pour les utilisateurs qui ne possèdent pas l'Image Processing Toolbox

Il est possible de superposer les contours en couleur à l'image initiale, en supposant que vous avez un tableau de contour BW (0 et 1) correspondant à celui retourné par EDGE et que I est le tableau de l'image initiale, comme ceci :

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
figure
colormap(gray)
imagesc(I);
axis image
 
hold on
 
cmap = [0 0 0   % Couleur RGB pour les valeurs 0 de BW1 (fond)
        1 0 0]; % Couleur RGB pour les valeurs 1 de BW1 (contour)
 
RGB(:,:,1) = BW*cmap(2,1)+~BW*cmap(1,1);
RGB(:,:,2) = BW*cmap(2,2)+~BW*cmap(1,2);
RGB(:,:,3) = BW*cmap(2,3)+~BW*cmap(1,3);
 
image(RGB,'alphadata',BW)
 
axis image