Discussion :

[magictol53]

Salut à tous,

Voici un script qui reprend une partie de la triangulation de delaunay. Ma question est dans le script.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
load croisement.txt;
m = [croisement(:,1),croisement(:,2),croisement(:,3),croisement(:,4)];
sample = m(randperm(200),:);
dt = DelaunayTri(sample(:,1),sample(:,2));
triplot(dt)
qrypts = [547148.4193	6238621.283 ]
triids = pointLocation(dt,qrypts);% j'identifie le triangle qui englobe qrypts
M = dt(triids,:)
M =
 
    46    24   109  % c'est ici où j'ai un doute:ces identifiants sont les trois sommets de mon triangle mais correspondent ils biens aux numeros de lignes de sample?

Merci d'avance

MG

[Jerome Briot]

Tu peux le vérifier graphiquement :

Code :

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

1
2
3
figure
plot(sample(dt([triids triids(1)],:),1),sample(dt([triids triids(1)],:),2),'b-',qrypts(1),qrypts(2),'g*-')
axis equal