Discussion :

[magictol53]

Bonjour,

Je tente de chercher le plus proche voisin de chaque point parmi une liste qui contient ce même point. La fonction nearest neighbors semble bien appropriée
après une triangulation. En revanche j'ai un souci, lorsque je sors le pid qui doit me renvoyer l'indice de son voisin le plus proche, je me retrouve avec des correspondances 1-1, ce qui n'est pas bon puisque cela indique que son plus proche voisin est lui même...J'ai pourtant essayé avec une sous matrice mais cela ne marche pas également..

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
load built.txt;
x1 = built(:,1);
x2 = built(:,2);
h = [x1 x2];
for n= 1:193;
        V = [h(1:n-1,:);h(n+1:193,:)];
        dt = DelaunayTri(V(:,1),V(:,2));
        qrypts = h(n,1:2);
        pid(n) = nearestNeighbor(dt, qrypts);
end

Merci de votre réponse

MT

[FLB]

Il te manque une ligne de code :
lorsque tu retires le point de la liste, tu décales les indices de tous les points suivants :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
x1 = rand(10,1);
x2 = rand(10,1);
h = [x1 x2];
pid = zeros(10,1);
for n= 1:10;
        V = [h(1:n-1,:);h(n+1:10,:)];
        dt = DelaunayTri(V(:,1),V(:,2));
        qrypts = h(n,1:2);
        pid(n) = nearestNeighbor(dt, qrypts);
        if (pid(n)>=n) % Ne pas oublier qu'on a enlevé un point!
            pid(n) = pid(n)+1;
        end
end