Discussion :

[ines ben alaya]

bonsoir,
j'ai 65 points uniformément répartis sur la sphère, je veux extraire 34 points (uniforme)à partir de 65 points.
pourriez vous m'aider à résoudre ce problème.
bien cordialement.

[tbc92]

Par curiosité, j'aimerais avoir la liste des 65 points d'origine.

Tu dis qu'ils sont uniformément répartis : C'est un maillage parfait, les distances 2-à-2 sont strictement égales, ou pas ? Et dans la 2nde option, ça veut dire que les 65 points sont presque uniformément répartis.

Et a priori, on est dans la 2nde option, car à ma connaissance, il n'existe pas de polyèdre régulier à 65 sommets.
Peux-tu déjà confirmer ce point.

[ines ben alaya]

Bonjour,
je vous remercie pour votre réponse.
vous trouverez ci-joint la liste des coordonnées des points.
je travaille dans le domaine de l'IRM de diffusion et l'uniformité des points est un critère très important pour pouvoir estimer le signal de diffusion.
j'ai utilisé cette application pour générer les 65 points.
Merci encore une fois pour votre aide.
Bien cordialement.

[ines ben alaya]

voici le lien de l'application
http://www.emmanuelcaruyer.com/q-space-sampling.php

[tbc92]

Bof...
J'ai regardé le fichier gradient65.
J'ai calculé les distances 2 à 2 entre tous les points. Et pour chaque point, j'ai regardé ses plus proches voisins. Je m'attendais à ce que, entre un point donné et ses 6 plus proches voisins, on obtienne 6 fois la même valeur, ou en tout cas 6 valeurs très proches les unes des autres. Pas de pot, ce n'est pas du tout le cas.
Par exemple, pour le 1er point (0.049, 0.919, 0.391) , ses pus proches voisins sont à des distances de 0.267, 0.327, 0.411, 0.428 ... ... Et on a des valeurs similaires pour les autres points.
Autrement dit, la définition de 'répartition uniforme' n'est pas du tout celle que j'avais imaginée.

En regardant le lien que tu as donné, on constate un truc : Quand on demande N valeurs, on a N points, puis s quand on sélectionne un N+x points, on a les mêmes N points, plus x nouveaux points.
Ca confirme que la définition de répartition uniforme n'est pas celle que j'attendais.

Mais, ça répond peut-être à la question initiale ?
Pour extraire 34 points parmi 65, il faudrait donc prendre les 34 premiers points du fichier ? ça paraît pas si mal, car si on prend le centre de gravité des N premiers points, on trouve toujours un point assez proche de (0,0,0)

[phryte]

Bonjour,

Je tirerais au hasard selon une distribution uniforme, 34 coordonnées parmi les 65.

En Matlab cela donne :

Code matlab :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
% Tira&ge de 34 points à partir de 65
M=load('gradient_65.txt');
plot3(M(:,1),M(:,2),M(:,3),'.r','MarkerSize',15)
hold on
[x,y,z]=sphere(34);
h=surfl(x,y,z)
set(h, 'FaceAlpha', 0.3)
shading interp
colormap(0.8*ones(2,3))
grid
axis equal
n=randperm(65,34);% permutation au hasard de 34 nombres parmi 65
M_34=M(n,:);% Nouvelle matrice de coordonnées
plot3(M_34(:,1),M_34(:,2),M_34(:,3),'.k','MarkerSize',15)