Discussion :

[jponsoda]

Bonjour,

J'essaie de créer une interpolation entre des points d'un scatter 3D en utilisant la triangulation de Delaunay.

Mon code est le suivant:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
x = xlsread('Database_PAZFLOR_DN450_rev06.xls','P01','A16:A872')
y = xlsread('Database_PAZFLOR_DN450_rev06.xls','P01','B16:B872')
z = xlsread('Database_PAZFLOR_DN450_rev06.xls','P01','C16:C872')
 
TRI=DelaunayTri(x,y,z)
A=[1:3958]';
F=TriScatteredInterp(TRI,A)
[qx,qy]=meshgrid(0:400,0:80);
qz= F(qx)
mesh(qx,qy,qz);
hold on;
plot3(x,y,z,'o');
grid on;

j'ai les erreurs suivantes:

??? Error using ==> TriScatteredInterp
The number of data point locations should equal the number of data point values.

Error in ==> Interpscatter at 55
F=TriScatteredInterp(TRI,A)

Je ne comprends pas!!!

en fait je sais qu'il faut un vecteur de type A dans mon code mais je ne sais pas à quoi il sert???

Je ne sais pas si j'utilise la bonne méthode: je voudrais en fait récupérer les coordonnées des points issues de la triangulation de delaunay en précisant moi même en quel abscisse et ordonnée je me place.

Merci d'avance pour votre aide.

[Dardie]

Bonjour,
le vecteur A représente les données à interpoler (A=f(x,y,z). Que cherches tu à faire exactement ?

Si j'ai bien compris, tu cherche à obtenir le z pour un (x,y) donné ?

Dans ce cas, il faut simplement appeler

Code :

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

1
2
3
4
5
6
7
8
 
 
F=Triscatteredinterp(x,y,z);
 
%vecteur de points à interpoler
[qx,qy]=meshgrid(mesabcisses,mesordonnées);
 
qz=F(qx,qy);

http://www.mathworks.com/access/help...terpclass.html

C'est très bien expliqué ici.

Si tu cherches à trouver des simplex dans la triangulation de delaunay, utilise tsearch ( http://www.mathworks.com/access/help.../tsearch.shtml )