Discussion :

[agondel]

Bonjour,

J'ai créé un programme où à partir d'une fontion (eigenmode) que je trace sur une grille 100x100 grâce à des données en format .txt je trouve le watershed de cette fonction.

Il se passe d'autres choses dans mon programme mais ce n'est pas vraiment important à savoir pour ma question.

En effet, voici mon problème :
j'ai un tableau 100x100 de logical appelé 'tab'
j'ai un tableau 100x100 de double appelé 'eigenmode'
je voudrais simplement superposer 'tab' en rouge sur 'eigenmode' (en gardant l'échelle de eigenmode)

J'ai essayé avec imfuse mais je n'obtiens pas ce que je veux, car les couleurs affichées ne correspondent pas à l'échelle si on l'affiche

comment faire ?

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
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
clear all
 
delimiterIn = ' ';
 
%%Landscape
 
landscape = importdata('landscape_plaqueL_homogene.txt', delimiterIn);
 
 
[r, id] = unique(landscape, 'rows');
z = r(:,3);
x = r(:,1);
y = r(:,2);
 
n = 100; 
[X, Y] = meshgrid(linspace(min(x)+0.005,max(x)-0.005, n), linspace(min(y)+0.005, max(y)-0.005,n));
landscape = griddata(x,y,z,X,Y);
 
% figure
% imagesc(landscape)
% axis image
% set(gca,'ydir','normal')
 
 
water = watershed(-landscape);
valley = water==0;
m=landscape(valley);
[yy,xx]=find(valley);
 
% figure
% imagesc(m)
 
 
 
%%Valeur Propre
valpropre = importdata('EigenValP.txt', delimiterIn);
 
%%Eigenmode
%%eigenmode
eigenmode= importdata('test.txt', delimiterIn);
[re, id] = unique(eigenmode, 'rows');
 
 
 
for k=1:3
 
 
  u=m(:,1)*valpropre(k,1);
 
  if max(u)>1
      tab=valley;
      for i=1:269
 
        if m(i,1)*valpropre(k,1)>1
 
 
       tab(yy(i,1),xx(i,1))=0;  
 
 
        end
 
      end
 
       figure
       imagesc(tab)
       title([num2str(k)]) 
 
       ze = re(7340*(k-1)+1:7340*k,5);
       xe = re(7340*(k-1)+1:7340*k,3);
       ye = re(7340*(k-1)+1:7340*k,4);
 
       n = 100; 
       [Xe, Ye] = meshgrid(linspace(min(xe)+0.005,max(xe)-0.005, n), linspace(min(ye)+0.005, max(ye)-0.005,n));
       eigenmode = griddata(xe,ye,ze,Xe,Ye);  
 
 
 
       figure
       c=imfuse(eigenmode,tab);
       imagesc(c)
 
 
 
  end
 
 
end

malheureusement mes fichiers .txt sont trop grands pour pouvoir être joints...

[Jerome Briot]

La méthode la plus simple consiste à se servir de tab comme masque de transparence pour eigenmode :

Code :

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

1
2
3
figure
imagesc(eigenmode, 'alphadata', ~tab);
axis image

La couleur est ensuite donnée par la couleur de l'objet Axes situé en dessous de l'objet Image :

Code :

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

set(gca, 'color', 'r');

[agondel]

super merci beaucoup !

[agondel]

Maintenant j'enregistre les images superposées,

et le "tab" apparaît en blanc et non en rouge comme je voudrais, comment faire ?

[Jerome Briot]

Dans ce cas, c'est un peu plus compliqué. Une solution consiste à modifier la valeur des éléments de eigenmode correspondant aux valeurs non nulles de tab. La nouvelle valeur devra correspondre à la position de la couleur rouge dans la palette de couleurs

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
eigenmode = rand(5);
tab = rand(5)>0.8
 
m = min(eigenmode(:));
M = max(eigenmode(:));
 
numCol = 64;
cmap = jet(numCol);
 
step = (M-m)/numCol;
 
redIndex = ismember(cmap, [1 0 0], 'rows');
redIndex = find(redIndex);
 
redValue = m+redIndex*step;
 
eigenmode(tab) = redValue;
 
figure
colormap(cmap)
imagesc(eigenmode)
axis image