Generate Random position that are not duplicate

Question

Riccardo Giaffreda 2020-2-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/507901-generate-random-position-that-are-not-duplicate

评论： Adam Danz 2021-12-4

Hello, I want to randomly chose N vector of two coordinates each (x and y, so N positions with thier coordinates) inside a circular area and then approximate the coordinates to its nearest integer (to have only positions such as 34;20 or 6;18 but not 6,27 ; 17,98). I did this with this code but I need the positions not to replicate (so I can't have two positions 10;15 for example). How can I do it?

N=30;
Radius= 50
x0 = 0;                   %Center of the cell in x direction
y0 = 0;                   %Center of the cell in y direction
%Set of Points
t = 2*pi*rand(N,1);
r = Radius*sqrt(rand(N,1));
x = x0 + r.*cos(t);
y = y0 + r.*sin(t);
x=round(x);
y=round(y);

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Temu 2020-2-27

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/507901-generate-random-position-that-are-not-duplicate#answer_417623

在 MATLAB Online 中打开

Hi Riccardo,

I would go for rejection sampling: just keep on generating until you have enough points.

N = 30;
Radius = 5;
%Set of Points
x1 = [];
while( size( x1, 1 ) < N )
    x1 = [x1; round( Radius .* ( rand( N, 1 ) * 2 - 1 )) + ...
        1j .* round( Radius .* ( rand( N, 1 ) * 2 - 1 ))];
    [~,ind1] = unique( x1 );
    x1 = x1(ind1);
end
x1 = x1(1:N);
x = real( x1 );
y = imag( x1 );

hth,

Temu

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Adam Danz 2020-2-27

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/507901-generate-random-position-that-are-not-duplicate#answer_417626

编辑：Adam Danz 2020-3-6

在 MATLAB Online 中打开

This efficient solution takes the following steps.

List all integer coordinates within the circle
Randomly choose N of those coordinates, without replacement (requires Stats & Machine Learning toolbox)

See inline comments for details.

% Set up inputs
N=30;
Radius= 50
x0 = 0;                   %Center of the cell in x direction
y0 = 0;                   %Center of the cell in y direction
% get all integer coordinates inside the square that frames the circle
xVec = floor(x0-Radius) : ceil(x0+Radius); 
yVec = floor(y0-Radius) : ceil(y0+Radius); 
[xGrid,yGrid] = meshgrid(xVec, yVec); 
% Eliminate coodinates outside of the circle
xy = [xGrid(:), yGrid(:)]; 
xyDist = pdist2(xy,[x0,y0]); 
xy(xyDist > Radius, :) = [];  % Use >= if you want to exclude points on circumf.
% xy now lists all integer coordinates inside the circle. 
% For visual inspection:  
clf(); 
plot(xy(:,1),xy(:,2), 'k.'); 
axis equal
% Now, select a random sample (N) without replication
% Requires Stats & Machine Learning Toolbox
xySelect = datasample(xy, N, 'Replace', false);
% Add chosen coordinates to the plot
hold on 
plot(xySelect(:,1),xySelect(:,2), 'ro')

The image below shows all possible integer coordinates within the cirlce (black dots) and the N randomly chosen coordinates (red circles).

2 个评论
显示无隐藏无

omar th 2021-12-4

how can I determine the position of each random points ?

Adam Danz 2021-12-4

If you plotted the points, you don't need to determine them since you've already got their coordinates.

If you didn't plot them and have the fig file, you can easily extract their coordinates from the file.

If this is a flat image (jpg, png, etc) then you can using image analysis tools to isolate the red points and the axis frame to estimate their position. Also see grabit from the file exchange.

请先登录，再进行评论。