Final answer :
I hope this will be helpful.
[label,n] = bwlabel(maskofVron4Inve);
stats = regionprops(label,'Area','BoundingBox','Centroid','Eccentricity','EquivDiameter','Perimeter');%'WeightedCentroid')
stats2Struct=struct2table(stats);
CenterXY=stats2Struct.Centroid(:,:);
for i=1:length(CenterXY)
clearvars nN2
CenterXYasRef=stats2Struct.Centroid(i,:);
[nN,dis] = knnsearch(CenterXY,CenterXYasRef,'k',5);
%sort the distances;
%// Sort the distances
Make_Dis_Tab(:,1)=nN;
Make_Dis_Tab(:,2)=dis;
indx_dis=find(Make_Dis_Tab(:,2)>50 & Make_Dis_Tab(:,2)<200);
Make_Dis_Tab2=Make_Dis_Tab(indx_dis,:);
nN2=Make_Dis_Tab2(:,1);
dis2=Make_Dis_Tab2(:,2);
%// Get the indices of the closest distances
hold on
%end
imshow(Im1);hold on
scatter(CenterXY(:,1),CenterXY(:,2))
plot(stats2Struct.Centroid(i,1), stats2Struct.Centroid(i,2),'y--*','MarkerSize',15);
hold on
line(CenterXY(nN2,1),CenterXY(nN2,2),'color','y','marker','o','linestyle','none','markersize',10)
%waitfor()
pause(1.2);
close
end