Discussion :

[Memboro]

Bonjour,

J'aurais également une question concernant les indices. En fait, j'ai une matrice nx3 dans laquelle est stockée environ 140000 points. Il s'agit d'une surface, donc par exemple, différents x peuvent avoir la même coordonnée y.

Je souhaite en fait trouver l'indice le plus proche d'une certaine valeur de y (jusque là, rien de trop compliqué), cependant je veux que cet indice corresponde également à une certaine valeur minimum pour x car je veux utiliser x et z par la suite.

J'ai donc essayé ce code, mais il me renvoie un index avec des valeurs un peu aléatoire, mais croissantes. J'ai essayé de varier les exigences (en mettant des conditions où je suis certain qu'il y a des points) mais ça ne marche toujours pas.

Avez-vous une suggestion?

Code :

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

1
2
3
4
5
6
7
8
9
 
% recherche de l'indice avec la valeur la plus proche de y_ref/2
for l = 1 : length(v(:,1))
    if v(l,1) > 5
        index =  find(abs(v(:,2)-y_ref/2) < 10);
    else 
        index = 0;
    end
end

Ici, v(:,1) correspond à x, v(:,2) à y et v(:,3) à z.

[le fab]

salut

ce qui répond le mieux à ta question :

Code :

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

[~, index] = min(abs(v(v(:,1)>5,2)-y_ref/2));

[ol9245]

je ne sais aps si cela s'approche de ta question que je n'ai pas comprise à 100%. le code ci-dessous recherche le point le plus proche, en (x,y) du point (0.5, 0.5), parmi les points donnés dans le tableau data.

Code :

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

1
2
3
4
data = rand(1000, 3) ;
point = [0.5, 0.5] ;
[unused, index] = min(sum((bsxfun(@minus, data(:,1:2), point)).^2, 2)) ;
disp(data(index, :))

[Memboro]

Bonjour,

Merci pour votre aide. Cependant, ce n'est pas tout à fait ce que je recherche.

En fait j'ai une grande matrice avec des coordonnées en x, y et z. Je veux trouver l'indice avec le point qui est le plus proche de y_ref/2, mais avec des conditions qu'en cet indice, les points x et z soient compris par exemple entre 9 et 15. Est-ce possible?

Merci pour votre aide!

[le fab]

Code :

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

[~, index] = min(abs(v(v(:,1)>=9 & v(:,1)<=15 & v(:,3)>=9 & v(:,3)<=15,2)-y_ref/2));

[Memboro]

Merci, par contre ça ne fonctionne pas...

Si j'ai bien compris, il recherche les indices possibles avec les conditions pour x et z, et ensuite, il prend l'indice parmi ceux sélectionnés avec la valeur la plus proche de y_ref/2?

J'ai contrôlé plusieurs fois manuellement et il y a plusieurs indices qui pourraient correspondre, donc le problème ne vient pas du fait qu'il n'en trouve pas...

Une idée d'où peut provenir ce problème?

[le fab]

en effet je suis allé trop vite

Code :

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

1
2
3
4
5
6
7
8
idxContraintes = v(:,1)>=9 & v(:,1)<=15 & v(:,3)>=9 & v(:,3)<=15;
[~, idxMinVlimite] = min(abs(v(idxContraintes ,2)-y_ref/2));
if ~isempty(idxMinVlimite)
  index = find(idxContraintes);
  index = index(idxMinVlimite);
else
 %pas de solutions
end

sinon tu as bien compris le raisonnement

[Memboro]

Génial!

Merci beaucoup pour ton aide!

[le fab]

en plus propre et plus juste :

Code :

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

1
2
3
4
5
6
7
idxContraintes = find(v(:,1)>=9 & v(:,1)<=15 & v(:,3)>=9 & v(:,3)<=15);
if ~isempty(idxContraintes)
  [~, idxMinVlimite] = min(abs(v(idxContraintes ,2)-y_ref/2));
  index = idxContraintes (idxMinVlimite);
else
 %pas de solutions
end

[Memboro]

Merci!!!