Discussion :

[lastrecrue]

je veux savoir comment crée une image couleur a partir de trois plan de couleur logiquement j'ai pensé a ça.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
function GA=Id(M)
I=imread(M);
R=I(:,:,1);
G=I(:,:,2);
B=I(:,:,3);
IN=rgb2gray(I);
[a,b]=size(IN);
R=floor(double(R));
G=floor(double(G));
B=floor(double(B));
GA=[R,G,B];
subplot(121);imshow(I);title('image _ M');
subplot(122);imshow(GA);title('image _ GR');

mais ça ne marche pas ça devais me donnée limage d'origine M sauf que non

[Jerome Briot]

Il n'est pas nécessaire de passer par les variables intermédiaires R, G et B ici :

Code :

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

1
2
3
4
5
6
7
8
9
function GA=Id(M)
I=imread(M);
IN=rgb2gray(I);
[a,b]=size(IN);
subplot(121);imshow(I);title('image _ M');
I(:,:,1)=floor(double(I(:,:,1)));
I(:,:,2)=floor(double(I(:,:,2)));
I(:,:,3)=floor(double(I(:,:,3)));
subplot(122);imshow(GA);title('image _ GR');

Dans le cas où ce serait vraiment nécessaire :

Code :

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

1
2
3
4
5
I=imread(M);
R=I(:,:,1);
G=I(:,:,2);
B=I(:,:,3);
GA=cat(3,R,G,B);

[lastrecrue]

j'ai lu le help du cat en fait cat(3,R,G,B) c'est la même chose que [R,G,B] et quand j'affiche avec cat ça me donne vraiment la même chose, donc il manque encore quelque chose...

[rostomus]

Au help je trouve ça:

Code :

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

cat(2, A, B) is the same as [A, B], and cat(1, A, B) is the same as [A; B].

et tu dis :

cat(3,R,G,B) c'est la même chose que [R,G,B]

relis le help s'il te plait

[ol9245]

Bonjour,
c'est intéressant les forums on apprend plein de choses. Je ne connaissais pas cat 3.

Jusqu'à présent je faisais :

Code :

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

rgb = reshape ([R G B], size(R,1), size(R,2), 3) ;

[lastrecrue]

vous voulez dire que ça marche cat(3,R,G,B) car pour moi ça ne veu pas marché?si non la fonction rgb = reshape ([R G B], size(R,1), size(R,2), 3) ; marche bien même si je suis pas très sure des résultats que ma gaussienne donne...

[rostomus]

Bon,

Je ne vois pas pourquoi ça ne marche pas pour toi. dis nous ce qui se passe, erreur?....?

sinon fais comme il faisait ol9245:

Code :

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

rgb = reshape ([R G B], size(R,1), size(R,2), 3) ;

[lastrecrue]

finalement si ça marche c'est moi qui a pas su m'en servir désole
voici mon code je crois que c'est pas vraiment ça le filtre gaussien????

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
27
function rgb=filtre_gausse(M)
I=imread(M);
R=I(:,:,1);
G=I(:,:,2);
B=I(:,:,3);
IN=rgb2gray(I);
[a,b]=size(IN);
R=floor(double(R));
G=floor(double(G));
B=floor(double(B));
GR=zeros(a-2,b-2);
GG=zeros(a-2,b-2);
GB=zeros(a-2,b-2);
for x=2:a-1;
    for y=2:b-1;
       GR(x-1,y-1)=uint8(double((R(x-1,y-1)+R(x,y-1)+R(x+1,y-1)+R(x-1,y)+R(x,y)+R(x+1,y)+R(x-1,y+1)+R(x,y+1)+R(x+1,y+1))/9)./255);
       GG(x-1,y-1)=uint8(double((G(x-1,y-1)+G(x,y-1)+G(x+1,y-1)+G(x-1,y)+G(x,y)+G(x+1,y)+G(x-1,y+1)+G(x,y+1)+G(x+1,y+1))/9)./255);
       GB(x-1,y-1)=uint8(double((B(x-1,y-1)+B(x,y-1)+B(x+1,y-1)+B(x-1,y)+B(x,y)+B(x+1,y)+B(x-1,y+1)+B(x,y+1)+B(x+1,y+1))/9)./255);
    end
end
R=GR;
G=GG;
B=GB;
rgb = reshape ([R G B], size(R,1), size(R,2), 3) ;
%rgb=cat(3,R,G,B);
subplot(121);imshow(I);title('image _ M');
subplot(122);imshow(rgb);title('image _filtré');

[Jerome Briot]

Pour bien montrer que les fonctions CAT et RESHAPE donnent le même résultat :

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
27
IM_0=rand(5,5,3);
 
R=IM_0(:,:,1);
G=IM_0(:,:,2);
B=IM_0(:,:,3);
 
IM_1=cat(3,R,G,B);
 
IM_2=reshape([R G B],size(IM_0));
 
 
figure
 
subplot(2,2,[1 2])
imagesc(IM_0)
title('Image originale')
axis image off
 
subplot(2,2,3)
imagesc(IM_1)
title('Image reconstituée avec CAT')
axis image off
 
subplot(2,2,4)
imagesc(IM_2)
title('Image reconstituée avec RESHAPE')
axis image off

[lastrecrue]

jai pas compri ca veut dire que c'est pas la meme chose ?

[Jerome Briot]

ca veut dire que c'est pas la meme chose ?

Le résultat est le même mais le temps d'exécution est différent :

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
clear all
 
IM=rand(5000,2500,3);
 
R=IM(:,:,1);
G=IM(:,:,2);
B=IM(:,:,3);
 
tic
X=cat(3,R,G,B);
toc
 
clear X
 
tic
XX=reshape([R G B],size(IM));
toc

Code :

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

1
2
3
4
5
6
7
8
elapsed_time =
 
    0.5000
 
 
elapsed_time =
 
    0.5620

Même si, sur une seule opération, la différence n'est pas très importante, la même opération à l'intérieur d'une structure itérative peut donner des écarts significatifs.