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

Version imprimable

02/05/2013, 17h25
Jimx59

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

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
03/05/2013, 12h17
Jimx59

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
03/05/2013, 12h55
Jerome Briot

Pourrais-tu attacher un fichier mat (mis dans une archive zip ou rar) contenant ces matrices pour faire quelques tests ?
03/05/2013, 18h29
duf42

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
14/05/2013, 09h52
Jimx59

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.
15/07/2013, 13h24
Jimx59

1 pièce(s) jointe(s)

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.
15/07/2013, 13h34
Jimx59

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?
16/07/2013, 11h08
Jimx59

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à :(
16/07/2013, 11h31
Jerome Briot

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é ?

1 pièce(s) jointe(s)

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

Code:

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.

16/07/2013, 14h16
Jimx59

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

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.
17/07/2013, 09h23
Jimx59

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 :mouarf:

Merci pour votre aide à duf42 et à toi. :ccool:
17/07/2013, 09h27
Jerome Briot

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.