У меня есть матрица 50943x3, содержащая вершины поверхностной сетки. Я хочу найти индекс определенной вершины, используя ее координаты (x, y, z).
Я попробовал найти функцию Matlab, но она возвращает пустую матрицу 0 на 1.
Заранее спасибо,
Ваше здоровье
Ваша попытка, вероятно, не работает из-за ошибок округления с плавающей запятой. Подробнее об этом можно прочитать здесь. Вы можете изучить функцию eps или просто использовать этот пример :
% Your matrix
M = randn(50943 , 3);
% The coordinates you are looking for
P = [0,0,0];
% Distance between all coordinates and target point
D = sqrt(sum((M - bsxfun(@minus,M,P)).^2,2));
% Closest coordinates to target
[~ , pos] = min(D);
% Display result
disp(M(pos , :))
Попробуйте следующее:
mat = randi(30,50943,3);
vec = [1,2,3];
% R2106b+ code
ind = find(all(mat==vec,2));
% or: explicit expansion, works with all versions
ind = find(all(bsxfun(@eq,mat,vec),2));
Что он делает: == или eq проверяет, равны ли координаты (дает логическую матрицу [50943x3]) all возвращает true, только если все координаты равны find возвращает индекс всех ненулевых элементов
Это работает только для точного совпадения (отсюда целочисленные координаты, выбранные с помощью randi).
Поскольку ответ уже принят, я добавлю ответ @Zep, который предлагает решение для получения ближайшей точки, которая, похоже, и была той, которую изначально искали.
[min_dist,ind_nearest] = min(sum(bsxfun(@minus,mat,vec).^2,2)); % index to the nearest point
vec не существует в матрице mat, как получить ближайшую вершину? 19.11.2018 randi и целочисленными координатами) 19.11.2018 
pos = find(D<tol);, гдеtol— подходящий допуск, также подойдет вместо[~ , pos] = min(D);. Он успешно вернет 2 индекса, если есть две совпадающие точки, и пустой массив, если их нет, в то время какminвсегда будет возвращать один индекс. Это зависит от того, можно ли ожидать этих случаев или нет, и что следует делать в этом случае. 19.11.2018D = sqrt(sum((M - P).^2,2));19.11.2018(M-P)наbsxfun(@minus,M,P), а не использоватьrepmat. Производительность лучше, а код короче. И обратно совместим до версии R2007a. 19.11.2018