Discussion :

[conanran]

Bonjour,
dans l'habitude, qauns j'ai une matrice 2D (Xmin*Ymin) de données Z et je veux afficher une cartographie j'utilise la fonction suivante

Code :

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

imagesc('CData', Z(:,:),'XData',(0:Xmax),'YData',(0:Ymax));

Maintenant, je dispose de données sous format X Y Z, c'est à dire 3 vecteurs colonne: abscisse - ordonnée - altitude.
Avez vous une idée comment je peux créer une carte de Z ? Sachant qu'on peut définir un pas x et pas y pour créer une grille.
merci

[Jerome Briot]

Tu trouveras quelques pistes ici : cartographie avec surf

[conanran]

Bonjour,
merci Jerome Briot
j'ai fait un essai selon les exemples données mais il m'affiche le message suivant :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
X=Zone_fine(:,1);
Y=Zone_fine(:,2);
Z=Zone_fine(:,3);
PasX =1;
Xi=min(Zone_fine(:,1)):PasX:max(Zone_fine(:,1));
PasY =1;
Yi=min(Zone_fine(:,2)):PasY:max(Zone_fine(:,2));
[Xi,Yi] = meshgrid(Xi,Yi); 
Zi = griddata(X,Y,Z,Xi,Yi);
figure
contourf(Xi,Yi,Zi)

[Jerome Briot]

Que renvoi ceci ?

Code :

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

whos Zone_fine

[conanran]

Voici le message :
......
ET pour faire une interpolation sur la grille, j'ai essayé :

Code :

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

1
2
F=scatteredInterpolant(X,Y,Z);
qz = F(Xi,Yi);

Mais il ne marche pas !
Merci!

[Jerome Briot]

Comme ceci :

Code :

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

1
2
3
X = double(Zone_fine(:,1));
Y = double(Zone_fine(:,2));
Z = double(Zone_fine(:,3));