How can I draw many circles without overlapping?
5 ビュー (過去 30 日間)
古いコメントを表示
The following code allows overlapping. How can I draw circles without overlapping?
t = linspace(0,2*pi,100);
figure
cx = 100 + (rand(20, 1) .* 1000); %center
cy = 100 + (rand(20, 1) .* 1000); % center
r = 50; % radius
plot(r*sin(t)+cx,r*cos(t)+cy)
y = r*cos(t)+cy
x = r*sin(t)+cx
for i =0:20
fill(x',y','k')
%fill(r*sin(t)+cx,r*cos(t)+cy, 'k')
r = r+1;
end
i=i+1;
axis([0,1200,0,1200])
axis square
0 件のコメント
採用された回答
michio
2018 年 1 月 12 日
Interesting problem. One way is to keep generating a set of centers until every center is separated by more than the fixed radius from any other centers. Simple logic. Here's an example function that does just that. Inputs are the number of centers required and the radius.
function centers = getCentersWithoutOverlaps(numCenters, radius)
%
% these are parameters you can set
a = 0;
b = 10;
%
radius2 = radius^2;
centers = zeros(numCenters,2);
centers(1,:) = a + rand(1,2).*b; % first center
for ii=2:numCenters
tmp = a + rand(1,2).*b;
% keep generating a random number until
% the new center is separated by "radius" from any other center points.
while any(sum((centers(1:ii-1,:)-tmp).^2,2) < radius2)
tmp = a + rand(1,2).*b;
end
% if the condition is met, the new center is added to "centers"
centers(ii,:) = tmp;
end
%
end
Please save above as a function and the following is the example script that you can try.
t = linspace(0,2*pi,100);
figure
r = 1; % radius
centers = getCentersWithoutOverlaps(10, 2);
cx = centers(:,1);
cy = centers(:,2);
plot(r*sin(t)+cx,r*cos(t)+cy);
y = r*cos(t)+cy;
x = r*sin(t)+cx;
fill(x',y','k')
axis square
Note that you need to specify reasonable values for numCenters and radius. Otherwise the function will not find the non overlapping set of centers.
0 件のコメント
その他の回答 (2 件)
H Y
2018 年 1 月 13 日
2 件のコメント
Manisha Kumari
2021 年 10 月 26 日
Hi ,
in the below code when i am putting axis([0,100,0,100]) then the above code is not working fine. could you please resolve this issue.
clc
close all
clear all
t = linspace(0,2*pi,100);
r = 1; % radius
centers = getCentersWithoutOverlaps(10,2);
cx = centers(:,1);
cy = centers(:,2);
plot(r*sin(t)+cx,r*cos(t)+cy);
y = r*cos(t)+cy;
x = r*sin(t)+cx;
%axis([0,100,0,100])
fill(x',y','r')
axis([0,100,0,100])
Thank you
MFB
2019 年 4 月 20 日
@michio
you put different radius in function calling i.e. radius 2 in the function as input argument and in the script you put r = 1. Can you explain why???
1 件のコメント
michio
2019 年 4 月 22 日
the r=1 in the script is just for visualization purpose. It can be any value. If you set r to be smaller than 2, you can draw circles without overlapping.
参考
カテゴリ
Help Center および File Exchange で Modeling についてさらに検索
製品
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!