Discussion :

[conanran]

Bonjour,

J’ai une matrice (A) 3D qui prend des valeurs à chaque instant t, mon but consiste à afficher les résultats en forme d’une cartographie en 3D pour chaque instant.
En 2 D j’utilise la fonction « imagesc » comme suit.

Code :

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

1
2
3
4
for n=1:length(variable.simulation_pas)
    imagesc('CData',A(:,:,n),'XData','YData');
 pause(0.02);
end

Je n’arrive pas à exploiter la fonction pour faire un affichage en 3D (comme l’image suivante)

[Jerome Briot]

Commence par faire une recherche avec le terme "voxel" dans le File Exchange

[conanran]

Merci Jerome Briot
Merci pour votre réponse, par contre je ne voix pas comment je peux l'utiliser dans mon cas, en effet, comment je peux tracer plusieurs voxels sans utiliser le boucle for et aussi comment je change la couleur selon l'intensité ( valeur de A).

[Jerome Briot]

Télécharge Voxel Image

Ouvre le fichier voxel_image.m et modifie la line 82 comme ceci :

Code :

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

patch('Vertices',vert,'Faces',fac,'FaceColor','flat','FaceVertexCData',color,'FaceAlpha',alpha,'Edgecolor',edgec);

Ensuite tu fais 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
18
19
% Valeurs arbitraires à afficher
A = rand(5,4,8);
 
% Couleur = valeur nomrmalisée de chaque voxel dans l'intervalle [1 64]
% 64 = nombre de couleurs dans la palette de couleurs de la figure
color = floor(1 + 63 * (A(:)-min(A(:)))./(max(A(:))-min(A(:))))
% Chaque voxel possède 6 faces de la même couleur
color = kron(color(:), ones(6,1));
 
% Génération de la grille des voxels
[x,y,z] = ndgrid(1:size(A,1),1:size(A,2),1:size(A,3));
 
figure
 
% Création des voxels
voxel_image([x(:) y(:) z(:)], [1 1 1], color);
 
view(3)
axis equal tight vis3d

[conanran]

Merci beaucoup Jerome Briot
colorbar.

Juste je veux savoir s'il est possible de changer le type de color map, je veux le changer par "jet" et ajouter un color bar.

Code :

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

1
2
3
4
5
6
c=colorbar; colormap(jet);
clim_max=Vmax; 
clim_min=Vmin;
ax=gca;
ax.CLim=[clim_min clim_max];
c.Label.String = 'Name';

Merci

[Jerome Briot]

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
18
19
20
21
22
23
24
25
26
27
28
29
% Valeurs arbitraires à afficher
A = rand(5,8,8);
 
% Couleur = valeur nomrmalisée de chaque voxel dans l'intervalle [1 num_color]
% num_color = nombre de couleurs dans la palette de couleurs de la figure
num_color = 64;
color = floor(1 + (num_color-1) * (A(:)-min(A(:)))./(max(A(:))-min(A(:))));
% Chaque voxel possède 6 facesd de la même couleur :
color = kron(color(:), ones(6,1));
 
% Génération de la grille des voxels
[x,y,z] = ndgrid(1:size(A,1),1:size(A,2),1:size(A,3));
 
figure
% Palette de couleur jet
colormap(jet(num_color))
 
% Création des voxels
voxel_image([x(:) y(:) z(:)], [1 1 1], color);
 
% Ajout colorbar
c = colorbar;
% Modification des étiquettes de la colorbar en fonction des valeurs de A
yt = get(c, 'ytick');
v = (yt-1)/(num_color-1)*(max(A(:))-min(A(:)))+min(A(:));
set(c, 'yticklabel', sprintf('%.2f\n', v))
 
view(3)
axis equal tight vis3d

[conanran]

Merci beaucoup Jerome Briot

c'est excellent !

comment je peux fixer les limites de colarbar, entre 0 et 10 dans votre exemple?

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
c = colorbar;
% Modification des étiquettes de la colorbar en fonction des valeurs de A
yt = get(c, 'ytick');
v =linspace(0,10,6);
% v = (yt-1)/(num_color-1)*(max(A(:))-min(A(:)))+min(A(:));
set(c, 'yticklabel', sprintf('%.2f\n', v))
 
% c=colorbar('Ticks',v);
c.Label.String = 'Name';
clim_max=10; 
clim_min=0;
CLims=[clim_min clim_max]; caxis(CLims);

Merci

[Jerome Briot]

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
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
% Valeurs arbitraires à afficher
A = rand(5,8,8);
 
% Couleur = valeur nomrmalisée de chaque voxel dans l'intervalle [1 num_color]
% num_color = nombre de couleurs dans la palette de couleurs de la figure
num_color = 64;
color = floor(1 + (num_color-1) * (A(:)-min(A(:)))./(max(A(:))-min(A(:))));
% Chaque voxel possède 6 facesd de la même couleur :
color = kron(color(:), ones(6,1));
 
% Génération de la grille des voxels
[x,y,z] = ndgrid(1:size(A,1),1:size(A,2),1:size(A,3));
 
figure
 
% Palette de couleur jet 
colormap(jet(num_color))
 
% Création des voxels
voxel_image([x(:) y(:) z(:)], [1 1 1], color);
 
clim_max = 10; 
clim_min = 0;
caxis([clim_min clim_max]);
 
% Ajout colorbar
c = colorbar;
% Modification des étiquettes de la colorbar en fonction des valeurs de A
yt = get(c, 'ytick');
v = (yt-1)/(num_color-1)*(max(A(:))-min(A(:)))+min(A(:));
set(c, 'yticklabel', sprintf('%.2f\n', v))
 
view(3)
axis equal tight vis3d

[conanran]

Merci
par contre, le but est de fixer les valeurs de colarbar entre 0 et 10.