Plot matrices (X,Y,Z et la valeur)

**Jimx59** · 02/05/2013, 18h25

Bonjour,

Voilà mon problème :
J'ai 4 matrices (5x9x4) :
- Les 3 premières sont des coordonnées X,Y et Z
- La dernière contient des valeurs prises à ces mêmes coordonnées

Je souhaite donc visualiser en 3D ma surface de mesurage (dans mon cas j'ai un parallélépipède) avec un code de couleur qui varie en fonction de la valeur contenue dans la 4ème matrice.

J'ai testé plot3, surface, contour3, colormap mais ces commandes servent à "ploter" 3 matrices, les 2 premières étant les coordonnées X et Y et la 3ème contenant la valeur.

Merci d'avance pour votre aide

**Jimx59** · 03/05/2013, 13h17

Un détail qui a son importance :
Les données que je possède sont les valeurs sur la surface de mon volume, en quelque sorte l'intérieur est creux.
Je veux donc juste visualiser l'extérieur d'une boite où la couleur de celle-ci dépend de la valeur à ce point.

Merci

**Jerome Briot** · 03/05/2013, 13h55

Pourrais-tu attacher un fichier mat (mis dans une archive zip ou rar) contenant ces matrices pour faire quelques tests ?

**duf42** · 03/05/2013, 19h29

Bonjour,

La fonction surf permet de spécifier 4 matrices où la 4ème matrice définit la couleur de chaque point.
Ca peut être une piste pour ton problème.

Duf

**Jimx59** · 14/05/2013, 10h52

Bonjour,

Merci pour votre réponse et désolé du retard, je teste la fonction surf et si ça ne marche pas, j'enverrai les fichiers mat.

Bonne journée.

**Jimx59** · 15/07/2013, 14h24

Bonjour,

Encore désolé, j'ai dû travailler sur autre chose et j'ai mis ce problème de coté.

A Dut, j'ai mis en PJ mon fichier de données. La 1ère colonne contient les X, la 2ème les Y et la 3ème les Z, représentant donc les coordonnées de mes points de mesures. En faisant un plot3 c'est plus parlant.
Les colonnes suivantes sont mes données à ces points.

J'aimerai donc visualiser la colonne 4 en fonction de X,Y et Z.

A duf42, j'ai testé la fonction surf. Mais apparemment cette fonction dessine Z en fonction de X et Y, et C permet d'ajuster la couleur.

Dans mon cas j'ai un parallélipipède qui représente ma surface de mesurage, j'ai plusieurs points de mesures aux coordonnées (X,Y,Z) qui sont sur la surface et pas à l'intérieur.
Je voudrai donc qu'aux valeurs les plus élevées, la zone autour du point soit plus rouge qu'aux points les moins élevées qui seraient plus bleus par exemple.

Merci d'avance.

**Jimx59** · 15/07/2013, 14h34

duf42, je viens de penser que la fonction surf pourrait être une solution en effet, mais ca nécessiterait des "bidoullages" que je ne sais pas faire.

Mon parallélépipède est constitué de 5 plans. Et, si je comprends bien la fonction surf, je peux dessiner 1 plan comme je le souhaite. il suffit que je garde les coordonnées d'un seul plan, puis dans Z je mets mes données à ces points.

Je peux ensuite faire de même pour les 4 plans suivants.

Le dernier problème serait ensuite de représenter l'ensemble de ces plans sur un même graphique et positionner l'ensemble de manière correcte. Mais ça, je ne sais pas le faire.

Qu'en pensez-vous?

**Jimx59** · 16/07/2013, 12h08

Bonjour,

J'ai testé la fonction surf sur un seul plan.

Ce n'est pas ce que je recherche, même si ce n'est pas loin.
Le résultat est des pics de couleur rouge aux valeurs élevées et des creux bleues aux valeurs faibles. Les couleurs sont bonnes mais je ne veux pas de pics et creux.

Je ne sais pas quoi essayer d'autres là

**Jerome Briot** · 16/07/2013, 12h31

La grille de mesure n'est pas régulièrement distribuée sur le domaine.

Est-ce normal ou bien dû au fait que toutes les mesures ne sont pas contenues dans le fichier que tu nous as donné ?

**Jerome Briot** · 16/07/2013, 12h53

Si les données fournies sont complètes, voici un début de solution (largement améliorable) :

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
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
load dir
 
x = dir_huile(:,1);
y = dir_huile(:,2);
z = dir_huile(:,3);
C = dir_huile(:,4); 
 
figure(1)
clf
plot3(x, y, z,'r*')
hold on
 
% MIN X
idx = find(x == min(x));
 
subx = x(idx);
suby = y(idx);
subz = z(idx);
subC = C(idx);
 
npoints = 30;
subxi = linspace(min(subx),max(subx),npoints);
subyi = linspace(min(suby),max(suby),npoints);
subzi = linspace(min(subz),max(subz),npoints);
 
subxi = ones(npoints)*min(x);
[subyi,subzi] = meshgrid(subyi,subzi);
 
subCi = griddata(suby,subz,subC,subyi,subzi);
 
surf(subxi,subyi,subzi,subCi)
 
% MAX X
idx = find(x == max(x));
 
subx = x(idx);
suby = y(idx);
subz = z(idx);
subC = C(idx);
 
npoints = 30;
subxi = linspace(min(subx),max(subx),npoints);
subyi = linspace(min(suby),max(suby),npoints);
subzi = linspace(min(subz),max(subz),npoints);
 
subxi = ones(npoints)*max(x);
[subyi,subzi] = meshgrid(subyi,subzi);
 
subCi = griddata(suby,subz,subC,subyi,subzi);
 
surf(subxi,subyi,subzi,subCi)
 
% MIN Y
idx = find(y == min(y));
 
subx = x(idx);
suby = y(idx);
subz = z(idx);
subC = C(idx);
 
npoints = 30;
subxi = linspace(min(subx),max(subx),npoints);
subyi = linspace(min(suby),max(suby),npoints);
subzi = linspace(min(subz),max(subz),npoints);
 
subyi = ones(npoints)*min(y);
[subxi,subzi] = meshgrid(subxi,subzi);
 
subCi = griddata(subx,subz,subC,subxi,subzi);
 
surf(subxi,subyi,subzi,subCi)
 
% MAX Y
idx = find(y == max(y));
 
subx = x(idx);
suby = y(idx);
subz = z(idx);
subC = C(idx);
 
npoints = 30;
subxi = linspace(min(subx),max(subx),npoints);
subyi = linspace(min(suby),max(suby),npoints);
subzi = linspace(min(subz),max(subz),npoints);
 
subyi = ones(npoints)*max(y);
[subxi,subzi] = meshgrid(subxi,subzi);
 
subCi = griddata(subx,subz,subC,subxi,subzi);
 
surf(subxi,subyi,subzi,subCi)
 
% MAX Z
idx = find(z == max(z));
 
subx = x(idx);
suby = y(idx);
subz = z(idx);
subC = C(idx);
 
npoints = 30;
subxi = linspace(min(subx),max(subx),npoints);
subyi = linspace(min(suby),max(suby),npoints);
subzi = linspace(min(subz),max(subz),npoints);
 
subzi = ones(npoints)*max(z);
[subxi,subyi] = meshgrid(subxi,subyi);
 
subCi = griddata(subx,suby,subC,subxi,subyi);
 
surf(subxi,subyi,subzi,subCi)
 
axis equal vis3d

Le pas de la grille d'interpolation a été choisi arbitrairement.

**Jimx59** · 16/07/2013, 15h16

Au vu de la PJ, ca a l'air d'être très prometteur

L'ensemble des mesures est présent dans le fichier que je t'ai donné. En fait, j'ai des micros placés en quinconce, ce qui cause des trous en mettant les données sous forme de matrice.
J'essayais de palier ce problème en utilisant la fonction griddata (c'était plutot mal parti pour le moment ^^').

Je vais regarder en détails le code que tu as fait. Je te remercie beaucoup pour ton aide.

Je tiens au courant.

**Jimx59** · 17/07/2013, 10h23

Bonjour Dut,

J'ai testé ton code, et regardé de plus près ce qu'il faisait en détails.

Apparemment c'est exactement la réponse à mon problème.

Si j'ai bien compris, tu commences par rechercher le min de l'axe X afin de construire la grille [Y;Z] en Xmin. Tu définis un nombre de points (ici 30) pour créer une grille plus fine que ma grille de mesures.

Tu utilises "meshgrid" pour créer cette grille de nouvelles coordonnées puis "griddata" pour interpoler mes données et remplir une nouvelle grille de données.

La fonction "surf" permet ensuite de "ploter" le résultat.

Tu réalises ça pour Xmin Xmax Ymin Ymax et Zmax, la fonction "hold on" permettant de tout mettre sur une seule figure.

Je me permets d'expliquer pour qu'un débutant comme moi puisse comprendre

Merci pour votre aide à duf42 et à toi.

**Jerome Briot** · 17/07/2013, 10h27

Exactement.

A toi maintenant d'améliorer ce code en adaptant le pas de la grille et le type d'interpolation utilisée dans griddata à ton problème physique

Il est également possible de le compacter car de nombreuses lignes sont similaires.

Plot matrices (X,Y,Z et la valeur) [Débutant]

MATLAB

Discussions similaires

Partager

Partager