Borrar filtros
Borrar filtros

Plot to avoid certain regions in 2D plot

1 visualización (últimos 30 días)
Ananya Malik
Ananya Malik el 2 de Nov. de 2018
Comentada: Ananya Malik el 2 de Nov. de 2018
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 comentario
KSSV
KSSV el 2 de Nov. de 2018
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?

Iniciar sesión para comentar.

Respuesta aceptada

KSSV
KSSV el 2 de Nov. de 2018
Editada: KSSV el 2 de Nov. de 2018
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')

Más respuestas (1)

madhan ravi
madhan ravi el 2 de Nov. de 2018
how about?
plot(loc(loc~=all(poly0.Vertices) & loc~=all(poly1.Vertices)))
  1 comentario
Ananya Malik
Ananya Malik el 2 de Nov. de 2018
Editada: Ananya Malik el 2 de Nov. de 2018
Thank you for the ans however, this is giving me 1D array of results. Cannot be plotted.

Iniciar sesión para comentar.

Categorías

Más información sobre Time and Frequency Domain Analysis en Help Center y File Exchange.

Productos

Community Treasure Hunt

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

Start Hunting!

Translated by