Présentation 3D de la bathymétrie d'une zone

[invite73486f18]

Bonjour tout le monde,
je suis nouveau sur ce forum et je voudrai de l'aide svp.
Je veux créer une présentation 3D de la bathymétrie d'une zone avec Matlab. Voilà ce que j'ai fais:

a=load('qc.txt');
x=a(:,1);
y=a(:,2);
z=a(:,3);
X=(min(x):0.0083:max(x));
Y=(min(y):0.0083:max(y));
[XI,YI]=meshgrid(X,Y);
Z=interp2(x,y,z,XI,YI);
mesh(x,y,w,XI,YI);

Un message d'erreur mentionnant "X and Y must be monotonic vectors or matrices produced by MESHGRID".
J'ai vérifié X et Y et ils sont tous les deux monotoniques.
Est ce que quelqu'un possède une solution svp
-----

[yoda1234]

Bonjour et bienvenue!

Merci d'utiliser les différentes balises à ta disposition: Dans ton cas il fallait employer la balise "Code" symbolisée par , je l'ai fait à ta place...
D'autre part j'ai déplacé ton message pour en faire un fil indépendant.

[invitea29b3af3]

Salut, bienvenue,

Tu pourrais mettre le fichier qc.txt en pièce jointe?

[invite73486f18]

Merci yoda1234 et excusez moi pour mon ignorance.
J'ai mis le fichier qc en pièce jointe sous format .zip
merci d'avance

[A voir en vidéo sur Futura]

[invitea29b3af3]

Re,

attention si tu regardes dans l'aide de interp2 (ZI = INTERP2(X,Y,Z,XI,YI)) les vecteurs X et Y auxquels se réfère l'erreur sont tes vecteurs x et y, mais x n'est pas strictement monotone, il est monotone par morceau. Fais plot(x(1:1000)) pour voir ça. Ou tape find(x==max(x)) et tu verras qu'il a 33 maximums.

[invite73486f18]

ah, merci fiatlux. Je vais essayer d'ordonner ça ! mille merci

[invitea29b3af3]

C'est rien.

[invite73486f18]

Grâce à vos précieuses commentaires j'ai réussi (plus ou moins) à faire une présentation 3D de la bathymétrie. J'ai eu recours à un tool qui s'appelle gridfit téléchargeable gratuitement. Voici le code

a=load('qc.txt');
x=a(:,1);
y=a(:,2);
z=a(:,3);
X=(min(x):0.0083:max(x));
Y=(min(y):0.0083:max(y));
g=gridfit(x,y,z,X,Y);
figure
colormap(hsv(256))
surf(X,Y,g);
camlight right;
lighting phong;
shading interp
%To limit axis
axis ([min(x) max(x) min(y) max(y) min(z) max(z)]);;

Merci pour tous ceux qui m'ont aidé