Discussion :

[Valou_le_lapinou]

Bonjour, mon problème va sûrement paraître simple à vos yeux mais je début et suis bloqué.

J'ai un fichier Excel contenant des coordonnées géographiques : latitude, longitude, altitude.

Je voudrais faire un graphique 3D en utilisant surf ou mesh.

Quand je fais surf(X,Y,Z) X étant la latitude, Y les longitude et Z l'altitude.
J'ai testé tout un tas de truc pris à droite à gauche mais j'en reviens toujours à la même erreur : "Z must be a matrix, not a scalar or vector."

Je crois comprendre que Z doit être fonction de X et Y mais ça ne l'est pas du tout dans mon cas. Comment faire un surf/mesh avec 3 données indépendantes?

Merci d'avance à la communauté.

[le fab]

salut

en effet pour utiliser surf ou mesh il faut des données complètes (cad pour toutes les combinaisons des X et Y présents dans tes données tu dispose du Z correspondant)
ce qui s'écrit via une matrice Z de taille length(X)*length(Y)
et d'ailleurs il faut parfois (selon les fonctions appelées) construire des matrices X et Y avec la fonction meshgrid

là tu as des données "éparpillées", c'est à dire un ensemble de triplets X, Y, Z
alors il se peut que des données soient complète, auquel cas il te suffit de réorganiser X, Y, et Z en matrice (cf reshape)

mais il est plus probable que tes données soient réellement éparpillée (scattered)
auquel ça tu peux :
- appeler une fonction faite pour ça : scatter3 ou plot3 (mais ca ne trace pas des surface)
- reconstituer un ensemble de données complet par interpolation (cf scatteredInterpolant)
- reconstituer un ensemble de données complet par lissage (cf gridfit)
- ...

Fabien

[Valou_le_lapinou]

Ok, merci pour votre aide.

Je crois que mon problème n'est pas solvable je vais rester sur le scatter3.

Bonne continuation
Valou_le_lapinou

[le fab]

si il est solvable, mais il faut faire un peu de mise en forme des données
essayes ceci (construction des points manquants par interpolation)

Code :

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

1
2
3
4
5
f = scatteredInterpolant(x, y, Z);
[mx, my] = meshgrid(x,y);
mz = f(mx, my);
figure;
surf(mx, my, mz)

[Valou_le_lapinou]

Cela fonctionne, merci beaucoup !

Bonne continuation