Discussion :

[ED2912]

Bonjour à tous et à toutes,

Je suis face à un obstacle: je trace une matrice 2D comportant un nombre important de points. En fait il s'agit de la forme d'un collimateur.

Ce que je souhaite c'est récupérer les coordonnées (x,y) des points qui sont à l'extérieur (de TOUS les points) de mon ensemble sachant que je connais la matrice (x,y) de tous les points formant mon ensemble. Donc le contour extérieur.

J'ai essayé convhull mais il ne m'offre qu'une approximation or je veux quelque chose d'ultra précis.

Ci-joint une image avec le contour de convhull.

Existe-t-il une solution me permettant de récupérer tous mes points?



Merci d'avance!

Ed

[Jean Dumoncel]

Bonjour,

si je comprends bien, tu ne possède pas le maillage ou une table de connexité de ces points? Ils ont l'air d'être placé sur une grille régulière, non?

[ED2912]

Exact. Je possède juste la matrice (x,y) de chaque points qui sont espacés du même intervalle dans les deux directions.

[Jerome Briot]

Quelques idées => Obtenir le contour d'une forme

[ED2912]

hmm oui j'avais pensé à implémenter un algo de ce style: je teste chaque ligne chaque colonne, dès que je rencontre un point je vérifie ces 8 voisins s'il en existe 7 ou moins je conserve le point s'il en a 8 je continue.
Mais à la base, j’eusse espéré qu'il existât une fonction toute faite ou une option de convhull qui me permette d'obtenir le résultat souhaité.

[Jean Dumoncel]

convhull permet d'obtenir l'enveloppe convexe, sans contrainte supplémentaire, tu ne pourras pas obtenir les zones concaves.

Comme ta grille est régulière, tu peux voir ton ensemble de points comme une image, et si tu as accès à la toolbox traitement d'image, tu peux utliser la fonction :

Code :

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

doc bwboundaries

qui d'ailleurs est basé sur l'algorithme de Moore qui est évoqué dans le lien pointé par Dut.

[ED2912]

Je ne suis pas sûr de comprendre ce que me retourne bwboundaries.
En tout cas il ne me retourne pas ce que je cherche.

[Jean Dumoncel]

Tu ne cherches pas les coordonnées des points de ton contour?

As-tu regardé l'exemple de la doc? Quel code as-tu utilisé?

[ED2912]

Voilà le code que j'ai utilisé mais je ne suis pas sûr de bien l'appliquer....
En fait si j'ai bien compris. L'argument de bwboundaries doit être une image binaire? Là en l’occurrence ce n'est pas ce que j'ai. En fait j'ai juste des points tracés en 2D issus de ma matrice colli_Matrix.

Code :

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

1
2
3
h=figure(1)
plot(colli_Matrix(:,1),colli_Matrix(:,2),'.');
B = bwboundaries(h);

[Jerome Briot]

Il "suffit" de passer du repère cartésien (x1...xn,y1...ym) au repère discret (1...mLigne,1...nColonne) comme ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
x = colli_Matrix(:,1);
y = colli_Matrix(:,2);
 
rows = 1+(y-min(y))/intervalle;
cols = 1+(x-min(x))/intervalle;
 
nrows = max(rows);
ncols = max(cols);
 
bw = zeros(nrows,ncols);
 
bw(rows,cols) = 1;

et ensuite :

Code :

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

B = bwboundaries(bw);

Le changement de repère inverse te donnera les valeurs (x,y) du "contour"

C'est l'idée...

[ED2912]

Et sinon si j'effectue le changement de repère comme tu le proposes, est-ce que contourc peut marcher aussi?

Bon en fait je récupère la matrice de 0 et de 1, sauf qu'en utilisant bwboundaries, je récupère quoi? Les indices de ma matrice pour lesquels bw(rows,cols)==1?

[ED2912]

J'avance un peu. Je récupère bien la matrice de mon contour mais c'est lors du retour au cartésien que je dois pas faire quelque chose de bien. Cf image et code:

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
xc = colli_Matrix(:,1);
yc = colli_Matrix(:,2);
 
xmin = min(xc);
ymin = min(yc);
 
B=bwboundaries(colliMat);
B{1,1}(1,:) = [];
 
rows = B{1,1}(:,1);
cols = B{1,1}(:,2);
 
x = (rows-1)*0.5+xmin;
y = (cols-1)*0.5+ymin;
 
figure(1)
plot(x,y,'r.');
hold on
 
plot(colli_Matrix(:,1),colli_Matrix(:,2),'.');

colliMat est la matrice de 0 et 1

[Jerome Briot]

Dans un premier temps :

Code :

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

1
2
rows = B{1,1}(:,2);
cols = B{1,1}(:,1);

La correspondance exacte est (x,y) <=> (colonne,ligne)