Discussion :

[Wel Kol]

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

Voici ce que j'ai codé :

Code :

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

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

[Invité]

Bonjour,

Les valeurs du tableau edge sont utilisées en tant qu'indices de ligne de node qui est de taille 3x2.
Ces valeurs ne doivent pas donc dépasser cette taille, ce que font np1, np2 et np3.
D'après ce que je vois dans le code de cette fonction (prise sur le FEX si je ne me trompe pas), tu peux fournir un edge vide, pour que ce soit calculé automatiquement.
Le code en peu simplifié :

Code :

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

1
2
3
4
5
6
7
% ...
tol = 1e-3;
for i = 1:numtri
    nT = dt(i,:); % ième triangle
    node = [x11(nT) y11(nT)];
    [in, on] = inpoly(A, node, [], tol);
end

[Wel Kol]

Merci pour votre proposition. ça marche nickel. Avez-vous une idée de comment faire pour créer une matrice dans laquelle pour chacun des triangles, on associe les points qui sont à l'intérieur et au bord?

[Invité]

À partir de la sortie in, tu peux appliquer l'indexage logique dessus points = A(in,:) pour obtenir les points situés à l'intérieur du triangle, ou find pour obtenir leur indices : indices = find(in).
Leur quantité étant variable selon le triangle considéré, il te faudra stocker ces données dans un tableau de cellules.

[Wel Kol]

Merci!