Discussion :

[Wel Kol]

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 :

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

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

[shaiHulud]

Bonjour,

Code :

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

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 :

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

k=dsearchn(CG,CG);

[Wel Kol]

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 :

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

k=dsearchn(CG,CG);

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

[shaiHulud]

doc dsearchn:

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

[shaiHulud]

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 :

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

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

[Wel Kol]

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));

[Wel Kol]

Avec le le code ci-après :

Code :

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

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?

[shaiHulud]

Code :

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

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.

[shaiHulud]

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