Discussion :

[SpyGer]

Bonjour à tous,

Merci par avance de l'aide que vous pourrez m'apporter.
Voila mon problème.

J'ai ceci:



Il y a un flux neutronique que j'ai calculer pour chaque cellule. (Si vous regardez bien il y a un quadrillage de 17x17 cellules).

J'ai, sous forme de matrice 17x17, la valeur du flux sur chaque cellule. J'aimerais le représenter en 3D avec comme plan horizontal celui de la figure (l'assemblage) et comme ordonné la valeur du flux (correspondant à chaque cellule). Chaque cellule fait 1.265 cm de long et j'aimerais que la légende (x,y) soit en cm, centré en (0,0) qui correspond au centre de la cellule centrale.

J'ai essayé d'être le plus clair possible. J'ai déjà commencer par utiliser: surfc, bar3

En gros j'arrive à un truc du genre:





Mais ici je n'ai que du 16x16 car dans le fond la valeur du flux est mise au intersection sur ce graphique alors que moi je voudrais vraiment des "histogrammes" sur un plan 2D.

Merci d'avance.

[Jerome Briot]

Une solution :

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
X = round(rand(17)*15);
 
figure
 
h = bar3(X,1);
 
for n=1:numel(h)
 
    % Centrage des valeurs en (0,0)
    x = get(h(n),'xdata');
    x = x-ceil(size(X,2)/2);
    y = get(h(n),'ydata');
    y = y-ceil(size(X,1)/2);
 
    % Modification de la couleur des barres en fonction de la valeur Z
    cdata=get(h(n),'zdata');
    cdata=repmat(max(cdata,[],2),1,4);
 
    set(h(n),'cdata',cdata,'facecolor','flat','xdata',x,'ydata',y);
 
end
 
axis tight
 
xt = get(gca,'xtick');
yt = get(gca,'ytick');
 
set(gca,'xticklabel',num2str(xt(:)*1.2625),'yticklabel',num2str(yt(:)*1.2625))

[SpyGer]

Merci beaucoup ca fonctionne plutôt bien.

Juste un peu de misère pour l'échelle quand je met la vue juste au dessus:

Voici le résultats (différents de celui du haut car je n'ai pas les mêmes données).





En tout cas merci bien!

[Jerome Briot]

Juste un peu de misère pour l'échelle quand je met la vue juste au dessus:

Il suffit de relancer les dernière lignes du code après la rotation du graphique :

Code :

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

1
2
3
4
xt = get(gca,'xtick');
yt = get(gca,'ytick');
 
set(gca,'xticklabel',num2str(xt(:)*1.2625),'yticklabel',num2str(yt(:)*1.2625))

[SpyGer]

Merci beaucoup