Discussion :

[magictol53]

Salut à tous,

Je vous explique, j'ai calculé une triangulation à partir de coordonnées x et y structurées sous cette forme

  x       y
12.23   25.45
12.56   25.89
12.89   36.74
.......   ........

La fonction delaunay me renvoie alors les nœuds des triangles sous cette forme

56  89  87
78  56  89
..............

Mon objectif est de récupérer les coordonnées de ces nœuds pour ensuite appliquer une autre transformation à l'aide de ces nouvelles combinaisons de points.

Merci d'avance

[Jerome Briot]

Sous quelle forme souhaites-tu récupérer les données ?

Ta question n'est pas très claire...

[magictol53]

Salut,

Admettons que la triangulation me renvoie mon premier triangle sous cette forme avec les identifiants des points:

45 56 89

Je souhaite recuperer une matrice sous cette forme (ou en transposé, peu importe)


x y x y xy

soit mes coordonnées de départ mais en prenant en compte la combinaison renvoyée par la triangulation.

Merci

[Jerome Briot]

Essaie ceci :

Code :

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

1
2
3
4
5
6
7
8
x = rand(1,3)
y = rand(1,3)
 
T = [1 2 3 ; 3 2 1 ]
 
T = T.';
xy = [x(T(:)) ; y(T(:))];
xy = xy(:)

ou

Code :

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

1
2
3
4
5
6
7
8
9
x = rand(1,3)
y = rand(1,3)
 
T = [1 2 3 ; 3 2 1 ]
 
T = T.';
xy = zeros(1,2*numel(T));
xy(1:2:end) = x(T(:));
xy(2:2:end) = y(T(:))

[magictol53]

Re,

C'est à peu près ça merci. En fait ton code me renvoie une matrice d'une seule colonne telle que

x = le x du 1er point du 1er triangle
x = le x du 2e point du 1er triangle
x = ..............
x = le x du 1er point du 2e triangle
x
y = le y du 1e point du 1er triangle
y
y
y
y

Comment pourrais maintenant le scinder en 2 et reformer un tableau du type

x y
x y
x y

Merci beaucoup en tout cas

PierrOt

[magictol53]

Et pour le suivi:

Code :

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

1
2
xy(1:2:end) = x(T(:));
xy(2:2:end) = y(T(:))

Ces deux lignes renvoient l'erreur suivante:

??? Index exceeds matrix dimensions.

Merci

[Jerome Briot]

Les codes que j'ai fournis renvoie un vecteur car c'est ce que tu semblais demander dans ton précédent message :

x y x y xy

Sinon les deux codes mélangent bien les coordonnées x et y comme ceci [x1 y1 x2 y2 x3 y3 x4 y4]

Comment pourrais maintenant le scinder en 2 et reformer un tableau du type

x y
x y
x y

Ce n'est toujours pas clair pour moi.
Je vois un tableau 3x2 correspondant à un triangle.
Quelles sont les dimensions de ce tableau pour N triangles ?

[magictol53]

Admettons que j'ai calculé une triangulation. Je voudrais maintenant obtenir ceci


x y
x y
x y
x y
Les trois premières lignes doivent alors correspondre aux x et y des sommets de mon premier triangle et ainsi de suite...

En partant de ce que m'a renvoyée la triangulation c'ad:

56 32 97
87 36 98
..................

------> Ces id correspondent aux identifiants des points de mon jeu de données ( premiere ligne pour le premier triangle). Chaque point a des coordonnées x et y et mon jeu d'entrée se structure donc sous cette forme

1 x y
2 x y
3 x y
4
5
6

J'éspère être clair, merci pour votre aide

PierrOt

[Jerome Briot]

C'est une simple variation des solutions que j'ai données précédemment :

Code :

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

1
2
3
4
5
6
7
x = rand(1,3)
y = rand(1,3)
 
T = [1 2 3 ; 3 2 1 ]
 
T = T.';
xy = [x(T(:)) ; y(T(:))].'

ou

Code :

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

1
2
3
4
5
6
7
8
9
x = rand(1,3)
y = rand(1,3)
 
T = [1 2 3 ; 3 2 1 ]
 
T = T.';
xy = zeros(numel(T),2);
xy(:,1) = x(T(:));
xy(:,2) = y(T(:))

[magictol53]

Merci beaucoup, cela marche formidablement bien...

A bientôt

PierrOt