Locate elements of a vector inside a meshgrid

1 次查看（过去 30 天）
Alberto Belvedere
Alberto Belvedere 2020-10-6
评论： Alberto Belvedere 2020-10-6
I'm trying to build a matrix RES, with the same dimensions of X and Y, that has a '1' in the position pointed (with a tollerance tollX and tollY) by the i-th couple contained in V.
[X,Y]=meshgrid(0:1:3,-2:1:2);
V=[2.1 1.2;
0.2 0.7;
3.1 1.9;
1.6 -1];
tollX=0.5;
tollY=0.5;
RES=0;
for i=1:length(V)
RES=(V(i,1)<X+tollX).*(V(i,1)>X-tollX).*(V(i,2)<Y+tollY).*(V(i,2)>Y-tollY)+RES;
end
This "rough" solution works well with small meshgrids and V, but since i have to manage far bigger data sets i would like to vectorize and refine the code to get better performances.
P.S. I expected '&' operator to be faster in general than ' .* ', but this doesn't seem to be true, at least for my case.
To test this I simply changed the statement inside the for loop with this one:
RES=((V(i,1)<X+tollX)&(V(i,1)>X-tollX)&(V(i,2)<Y+tollY)&(V(i,2)>Y-tollY))|RES;

请先登录，再回答此问题。

回答（1 个）

Rik
Rik 2020-10-6
I would suggest using ismembertol, or consider functions like normxcorr2 from the image processing toolbox.
  3 个评论
显示 1更早的评论
Rik
Rik 2020-10-6
I'm not quite sure how implicit expansion (with bsxfun) would increase performance. I suspect ismembertol would be more efficient. I also don't necessarily see how find would be useful.
Alberto Belvedere
Alberto Belvedere 2020-10-6
Thanks, i'll do some tests to see which one performs better.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Resizing and Reshaping Matrices 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by