I have to plot some random points in 500 x 500 area by avoiding certain area (poly0 and poly1 in the example below) in the 2D-plot.
x = 350.5; y = 350.5; r = 50;
th = 0:pi/50:2*pi;
xunit = r * cos(th) + x;
yunit = r * sin(th) + y;
poly0 = polyshape(xunit(1:end-1),yunit(1:end-1));
plot(poly0);
hold on;
ploy1 = polyshape([100 100 200 200],[100 200 200 100]); %square
plot(poly1);
hold on
loc = randi([1,500],50, 2);
If I do this, random points overlap the poly0 and poly1 areas. Can anyone help? TIA.

1 Kommentar

x = 350.5; y = 350.5; r = 50;
th = 0:pi/50:2*pi;
xunit = r * cos(th) + x;
yunit = r * sin(th) + y;
poly0 = polyshape(xunit(1:end-1),yunit(1:end-1));
plot(poly0);
hold on;
poly1 = polyshape([100 100 200 200],[100 200 200 100]); %square
plot(poly1);
hold on
loc = randi([1,500],50, 2);
There was a typo error. What you want to avoid?

Melden Sie sich an, um zu kommentieren.

 Akzeptierte Antwort

KSSV
KSSV am 2 Nov. 2018
Bearbeitet: KSSV am 2 Nov. 2018

1 Stimme

Read about inpolygon
x = 350.5; y = 350.5; r = 50;
th = 0:pi/50:2*pi;
xunit = r * cos(th) + x;
yunit = r * sin(th) + y;
poly0 = polyshape(xunit(1:end-1),yunit(1:end-1));
plot(poly0);
hold on;
poly1 = polyshape([100 100 200 200],[100 200 200 100]); %square
plot(poly1);
hold on
loc = randi([1,500],50, 2);
% avoid poly0
idx = inpolygon(loc(:,1),loc(:,2),poly0.Vertices(:,1),poly0.Vertices(:,2)) ;
loc(idx,:) = NaN ;
% avoid poly1
idx = inpolygon(loc(:,1),loc(:,2),poly1.Vertices(:,1),poly1.Vertices(:,2)) ;
loc(idx,:) = NaN ;
plot(loc(:,1),loc(:,2),'.r')

Weitere Antworten (1)

madhan ravi
madhan ravi am 2 Nov. 2018

0 Stimmen

how about?
plot(loc(loc~=all(poly0.Vertices) & loc~=all(poly1.Vertices)))

1 Kommentar

Ananya Malik
Ananya Malik am 2 Nov. 2018
Bearbeitet: Ananya Malik am 2 Nov. 2018
Thank you for the ans however, this is giving me 1D array of results. Cannot be plotted.

Melden Sie sich an, um zu kommentieren.

Kategorien

Mehr zu 2-D and 3-D Plots finden Sie in Hilfe-Center und File Exchange

Produkte

Community Treasure Hunt

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

Start Hunting!

Translated by