Plot to avoid certain regions in 2D plot

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 件のコメント

KSSV
KSSV 2018 年 11 月 2 日
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?

サインインしてコメントする。

 採用された回答

KSSV
KSSV 2018 年 11 月 2 日
編集済み: KSSV 2018 年 11 月 2 日

1 投票

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')

1 件のコメント

Ananya Malik
Ananya Malik 2018 年 11 月 2 日
Thanks a lot.. Exact solution as always (y).

サインインしてコメントする。

その他の回答 (1 件)

madhan ravi
madhan ravi 2018 年 11 月 2 日

0 投票

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

1 件のコメント

Ananya Malik
Ananya Malik 2018 年 11 月 2 日
編集済み: Ananya Malik 2018 年 11 月 2 日
Thank you for the ans however, this is giving me 1D array of results. Cannot be plotted.

サインインしてコメントする。

カテゴリ

ヘルプ センター および File ExchangePolar Plots についてさらに検索

製品

質問済み:

2018 年 11 月 2 日

コメント済み:

2018 年 11 月 2 日

Community Treasure Hunt

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

Start Hunting!

Translated by