Remove centroid locations of objects that are too close

6 visualizaciones (últimos 30 días)
Jason
Jason el 8 de Ag. de 2022
Comentada: Jason el 9 de Ag. de 2022
Hi, I have a list of centroid locations of spots in an image (xf, yf). The spots are randomly distributed. I want to remove the centroid locations that are too close to each other, or say within a distance of d. I've got to the following but not sure what to do next. Also is there a better way to do this ratrher than use loops?
%using a loop
M=[];
for j=1:length(xf)
for i=1:length(xf)
dx=xf(j)-xf(i);
dy=yf(j)-yf(i);
dr=sqrt(power(dx,2)+power(dy,2));
M(j,1)=xf(j);
M(j,2)=yf(j);
M(j,3)=xf(i);
M(j,4)=yf(i);
M(j,5)=dr;
end
end
  2 comentarios
Matt J
Matt J el 9 de Ag. de 2022
If two spots are less than d apart should both be removed? If not, by what criterion do you decide which spot remains?
Image Analyst
Image Analyst el 9 de Ag. de 2022
To me it sounds like both are to be removed. This is easily done with pdist2.

Iniciar sesión para comentar.

Respuesta aceptada

Matt J
Matt J el 9 de Ag. de 2022
Editada: Matt J el 9 de Ag. de 2022
A=[xf(:) yf(:)];
D=pdist2(A,A);
D(D==0)=inf; %retroactive EDIT
keep = all(D>d,1);
xf=xf(keep);
yf=yf(keep);
  7 comentarios
Jason
Jason el 9 de Ag. de 2022
Of course, its clear now, thankyou
Jason
Jason el 9 de Ag. de 2022
Just to add, when I had over 400 objects, akthough i do other operations, using the for loop took 110s.using the pdist2 approach took just under 2s. Remarkable!

Iniciar sesión para comentar.

Más respuestas (1)

yanqi liu
yanqi liu el 9 de Ag. de 2022
yes,sir,may be use pdist2, such as
xf = rand(1, 10);
yf = rand(1, 10);
%using a loop
M=[];
for j=1:length(xf)
for i=1:length(xf)
dx=xf(j)-xf(i);
dy=yf(j)-yf(i);
dr=sqrt(power(dx,2)+power(dy,2));
M(end+1,1)=xf(j);
M(end,2)=yf(j);
M(end,3)=xf(i);
M(end,4)=yf(i);
M(end,5)=dr;
end
end
% second method
A=[xf(:) yf(:)];
B=pdist2(A,A);
C=B';
C=C(:);
isequal(M(:,end),C)
ans = logical
1

Categorías

Más información sobre Timing and presenting 2D and 3D stimuli en Help Center y File Exchange.

Etiquetas

Productos


Versión

R2022a

Community Treasure Hunt

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

Start Hunting!

Translated by