Identifier les points les plus proches

Bonjour,

J'ai une matrice qui contient les cordonnées des points. J'aimerais pour chaque point de la matrice identifier les points les plus proches. J'ai utilisé la fonction dsearchn. Ci-après le code que j'ai écrit :

Code:

---

1 2 3 4 5 6 7

num = size(CG,1);   for h=1:num checkpoint(1,1)=CG(h,1); checkpoint(1,2)=CG(h,2); k=dsearchn(checkpoint,CG); end

---

Ça marche mais je n'obtiens pas de bons résultats.
Pourriez-vous me faire des suggestions?
Merci

Bonjour,

Code:

---

checkpoint(1,1)=CG(h,1); checkpoint(1,2)=CG(h,2);

---

est inutile au sein de la boucle. Il suffit de faire checkpoint= CG; en dehors, mais ce n'est même pas la peine ici. Il ne faut pas appeler dsearchn au sein de la boucle, car dsearchn est vectorisé. Il suffit de faire

Code:

---

k=dsearchn(CG,CG);

---

Peut être que je ne comprends pas bien votre solution. Je pense que j'ai besoin d'une boucle pour identifier pour chacun des points de la matrice les points voisins.
Si je mets directement ce code

Code:

---

k=dsearchn(CG,CG);

---

J'obtiens comme résultat tous les points de la matrice.

doc dsearchn:

Citation:

---

k = dsearchn(X,T,XI) returns the indices k of the closest points in X for each point in XI

---

Autant pour moi, ce que j'ai dit est vrai, mais du coup le plus proche voisin est toujours lui-même ! Du coup, autant abandonner cette fonction et faire le travail soi même:
1/ calculer une matrice de distance (via boucle ou bsxfun)
2/ mettre l'infini dans la diagonale
3/ prendre le min sur chaque colonne via min()

Sinon, il faut effectivement exclure le point courant à chaque fois

Code:

---

1 2 3 4 5 6

k=nan(1,num); for h=1:num CG_sansx=CG; CG_sansx(h,:)=[]; % suppression du point courant k(h)=dsearchn(CG_sansx,CG(h)); % distance entre point courant et les autres end

---

Bonjour,

Merci pour vos propositions. J'ai utilisé votre proposition 2. Quand je lance le calcul, ce message d'erreur suivant apparaît :

??? Error using ==> dsearchn at 61
X and XI should have the same column dimension. Error in ==> test at 72
k(h)=dsearchn(CG_sansx,CG(h));

Avec le le code ci-après :

Code:

---

1 2 3 4 5 6

num = size(CG,1); for i=1:num point(1,1:2)=CG(i,1:2); CG(i,:)=[]; k(i)=dsearchn(CG,point); end

---

Le message d'erreur suivant apparaît :" ??? Index exceeds matrix dimensions. Error in ==> test at 71 point(1,1:2)= CG(i,1:2);"
Pourriez-vous m'expliquer comment résoudre ce problème?

Code:

---

1 2 3 4 5 6 7 8

CG= rand(5,3); % random data num = size(CG,1); for i=1:num point=CG(i,:); others=CG; others(i,:)=[]; k(i)=dsearchn(others,point); end

---

Attention, cela renvoie l'indice dans la sous matrice others.

Le problème était que tu détruisais les lignes de CG dans la boucle sans les restaurer.