Utilisation de la syntaxe longue inpoly
Bonjour,
J'utilise la fonction inpoly pour identifier les points appartenant à chacun des triangles de mon maillage obtenu à partir de la triangulation de Delaunay. Quand j'utilise la syntaxe courte inpoly(p, node), ça marche mais lorsque j'utilise la syntaxe longue inpoly(p, node, edge, tol), le message d'erreur suivant apparaît :
??? Error using ==> inpoly at 96, Invalid EDGE
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
| x1 = A(:,1);
y1 = A(:,2);
x11 = A(1:47:end,1);
y11 = A(1:47:end,2);
dt = DelaunayTri(x11,y11);
numtri = size(dt,1);
for i = 1:numtri
% Indexation linéaire des points des triangles
nT =dt(i,:); % ième triangle
np1=dt(i,1);
np2=dt(i,2);
np3=dt(i,3);
%Coordonnées des points au sommet des triangles
x_p1=x11(np1);
y_p1=y11(np1);
x_p2=x11(np2);
y_p2=y11(np2);
x_p3=x11(np3);
y_p3=y11(np3);
p1=[x_p1,y_p1];
p2=[x_p2,y_p2];
p3=[x_p3,y_p3];
% Utilisation de inpoly
p=[x1,y1];
node=[p1(1) p1(2); p2(1) p2(2); p3(1) p3(2)];
edge=[np1 np2; np2 np3; np3 np1];
tol=1e-3;
[in, on]=inpoly(p, node, edge, tol);
end |
Pourriez-vous m'aider à résoudre ce problème provenant de "edge".
Merci
