フィルターのクリア

Plot to avoid certain regions in 2D plot

2 ビュー (過去 30 日間)
Ananya Malik
Ananya Malik 2018 年 11 月 2 日
コメント済み: Ananya Malik 2018 年 11 月 2 日
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 日
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 日
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.

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

カテゴリ

Help Center および File ExchangeTime and Frequency Domain Analysis についてさらに検索

製品

Community Treasure Hunt

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

Start Hunting!

Translated by