Discussion :

[Polybead]

Bonjour à tous,


Je dipose d'un vecteur qui contient des coordonnées de points. J’aimerais calculer pour chacun d’entre eux la distance aux autres points et si celle-ci est inférieure à une distance 'd' alors placer cette distance dans un vecteur. Ce dernier serait donc constitué d’un nombre de lignes égal au nombre de points et le nombre de distances sur chaque ligne renseignera sur le nombre de points situés à une distance inférieure à 'd'.

je bloque un peu sur l'écriture du script...

Merci d'avance.

Piet

[Caro-Line]

As-tu commencer à coder quelque chose ?

Comment calcules-tu déjà la distance entre 2 points ?
Une fois que tu as cela il suffit de boucler sur ton vecteur, quelque chose du style :

Code :

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

1
2
3
4
5
6
7
distances=[];
for k=1:length(V)-1
  d=calcul_distance(V(k),V(k+1));
  if d < seuil
     distances(end+1)=d;
  end
end

[Polybead]

en fait, ce sont les distances euclidiennes qui m'intéressent
j'avais tapé les écritures suivantes mais je n'ai aucune idée si c'est correct

Code :

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

1
2
3
4
5
 for n=1:N
   if pdist(XY,'euclidean')<14
      Z(n)=pdist(XY,'euclidean')
   end
end

où XY est mon vecteur contenant les coordonées des points

merci pour ton aide, je vais tester cela ce soir

[Caro-Line]

Attention tu n'as pas mis d'indice à XY. Donc pdist(XY,'euclidean') aura toujours la même valeur à chaque itération de ta boucle et de plus donne un vecteur.
De plus la façon dont tu fais va donner un vecteur Z de taille 1*N avec des 0 quand la condition ne sera pas respectée et je n'avais pas l'impression que c'était ce que tu souhaitais.

J'ai regardé un peu la doc de pdist : elle a l'air de prendre en entrée une matrice plutôt qu'un vecteur...
De quelle taille est XY ?

Une idée :

Code :

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

1
2
Z=pdist(XY,'euclidean');%donne un vecteur tel que décrit dans la doc
Z(Z<14)=[];%enlève de Z toutes les valeurs inférieures à 14

[Polybead]

il y a environ 2800 points

[Polybead]

XY est de type (2x2800)

[Caro-Line]

XY est de type (2x2800)

Donc c'est bien une matrice
C'était juste un problème de vocabulaire, un vecteur n'est finalement qu'une matrice spécifique à 1 dimension : 1*N.
La méthode que je t'ai donnée devrait donc fonctionner.

[Polybead]

oui en effet, désolé j'avais vu que j'avais indiqué vecteur mais c'est en effet une matrice