Discussion :

[romainmill]

Bonsoir,

Voici mon problème : j'affiche une image en noir et blanc sur lequel j'aimerai afficher un champ de vitesse en couleur. Mon code initial est :

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
28
29
30
31
32
33
34
35
 
 
im=imread(strcat(num2str(lat1),'_',num2str(lat2),'.TIF'));
 
displ=load (strcat(num2str(lat1),'_',num2str(lat2),'.out'));%je charge mon fichier contenant mon champ de vitesse
 
%stockage des données dans des variables
x=displ(:,1);
y=displ(:,2);
dep=displ(:,3);
dx=displ(:,6);
dy=displ(:,7);
 
figure
imagesc(im,'XData',[664170 734790], 'YData', [3889890 3981270]);
 
axis on
 
set(gca,'xtick',664170:10000:734790);
set(gca,'ytick',3889890:10000:3981270);
axis([664170 734790 3889890 3981270])
 
hold on
 
%Calcul de la norme du deplacement
N=sqrt((dx.^2)+(dy.^2));
N=N.*60*(369/365);
 
X=[664170:15:734790]';
Y=[3889890:15:3981270]';
x_N=X(x);
y_N=Y(y);
 
 
b=scatter(x_N,y_N,60,N,'filled')

Mais mon champ de déplacement est afficher en gris, alors que j'aimerai avoir : colormap('jet') pour le champ de déplacement.
Du coup j'ai modifier mon code en ajoutant un freezecolors :

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
28
29
30
31
32
33
34
35
36
 
 
im=imread(strcat(num2str(lat1),'_',num2str(lat2),'.TIF'));
 
displ=load (strcat(num2str(lat1),'_',num2str(lat2),'.out'));
 
%stockage des données dans des variables
x=displ(:,1);
y=displ(:,2);
dep=displ(:,3);
 
figure
%image(im,'XData',[664170 734790], 'YData', [3889890 3981270]);
imagesc(im,'XData',[664170 734790], 'YData', [3889890 3981270]);
colormap('gray')
freezeColors
axis on
 
set(gca,'xtick',664170:10000:734790);
set(gca,'ytick',3889890:10000:3981270);
axis([664170 734790 3889890 3981270])
 
hold on
 
%Calcul de la norme du déplacement
N=sqrt((dx.^2)+(dy.^2));
N=N.*60*(369/365);
%N=N.*60*(401/365);
 
X=[664170:15:734790]';
Y=[3889890:15:3981270]';
x_N=X(x);
y_N=Y(y);
 
 
b=scatter(x_N,y_N,60,N,'filled')

Mais au niveau du scatter matlab plante et n'arrive pas à m"afficher le scatter sur l'image.

Au tout départ, cela marchait si j'utilisé :

Code :

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

1
2
 
im=ind2rgb(im,gray(256))

Mais cela prend beaucoup trop de mémoire à l'ordinateur.

Auriez vous des conseils à me donner pour pouvoir afficher mon champ de déplacement en couleur sur mon image en noir et blanc (peut etre utiliser autre chose que imagesc, conseils pour utilisation de freezecolors...)?

Merci d'avance

Romain

[Jerome Briot]

Peux-tu nous montrer ce que renvoie ceci ?

Code :

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

1
2
3
im = imread(strcat(num2str(lat1),'_',num2str(lat2),'.TIF'));
 
whos im

[romainmill]

Code :

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

1
2
3
4
 
  Name         Size                 Bytes  Class    Attributes
 
  im        6092x4708            28681136  uint8

[Jerome Briot]

Et quelle est la définition (résolution) du ou des écrans sur lequel tu comptes afficher l'image ?

[romainmill]

La résolution de mon écran est de : 1366x768.
En fait maintenant cette solution marche : ,mais bon ça rame à fond pour manipuler les images :

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
28
29
30
31
32
 
 
im=imread(strcat(num2str(lat1),'_',num2str(lat2),'.TIF'));
 
displ=load (strcat(num2str(lat1),'_',num2str(lat2),'.out'));
 
%stockage des données dans des variables
x=displ(:,1);
y=displ(:,2);
dep=displ(:,3);
 
figure
%image(im,'XData',[664170 734790], 'YData', [3889890 3981270]);
b=imagesc(im)
colormap('gray')
freezeColors(gca)
%Calcul de la norme du déplacement
N=sqrt((dx.^2)+(dy.^2));
N=N.*60*(369/365);
mask=imread(strcat(num2str(lat1),'_',num2str(lat2),'_15mask.TIF'));
M = [];
for n = 1:numel(x)
 
    if mask(y(n),x(n))>0
 
        M = [M ; x(n) y(n) N(n)];
    end
 
end
hold on
b=scatter(M(:,1),M(:,2),20,M(:,3),'filled')
colormap('jet')

[Jerome Briot]

La résolution de mon écran est de : 1366x768.

Tu conviendras donc qu'il est inutile de vouloir afficher 6092x4708 pixels sur un écran qui n'en comporte que 1366x768.

Si j'étais à ta place, je redimensionnerais la taille de l'image avant de l'afficher et j'utiliserais la technique que tu cites dans ton premier message :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
I = imread(...);
 
Ir = imresize(I,...);
 
Irgb = ind2rgb(Ir, ...);
 
figure
 
image(Irgb)
hold on
scatter(...)

[romainmill]

Ok ça à l'air de marcher en divisant la résolution par 4, merci !
Du coup il va falloir ré-échantilloner mon champ de déplacement avec reshape par exemple ?

[Jerome Briot]

Pas si tu utilises la ligne suivante :

Code :

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

image(im,'XData',[664170 734790], 'YData', [3889890 3981270]);

[romainmill]

Ouai en effet ca marche merci!