How to plot points outside certain polygon area in MATLAB

Question

I have been looking around for a method to plot points outside a polygon's area(hexagon in my case). Here's the scenario that I want to achieve, I have a small hexagon located inside a big hexagon. The picture is as follows:

In the picture, I created a small hexagon (whom area indicated in pale red) and generate a random points (three in my case) inside it using inpolygon. Problem arise when I want to plot points (red triangles) in big hexagon (indicated in pale purple) without touching the small hexagon area. I have look around the net for this simple solutions 3 days to no avail.

I would really appreciate any helps or guidance I could get. Thank you so much!

My code is as follows:

clear
clc

bighexagon = 20;
smallhexagon = 4;

axis_min = 0;
axis_max = 40; 
axis([axis_min axis_max axis_min axis_max],'square');
hold on

L = linspace(30,390,7); 
bhex_x = bighexagon * (1+cosd(L))'; 
bhex_y = bighexagon*(1+sind(L))';

L2 = linspace(30,390,7); 
shex_x = smallhexagon * (1+cosd(L2))'; 
shex_y = smallhexagon * (1+sind(L2))';

plot(bhex_x,bhex_y,'LineWidth',3);

%---Move small hexagon into big hexagon
shex_vertices_x2(:,1) = shex_x + 16;
shex_vertices_y2(:,1) = shex_y + 16;
plot(shex_vertices_x2(:,1),shex_vertices_y2(:,1),'--k','LineWidth',3);


%---Plot points in small hexagon
no = 3;
point_x2 = (smallhexagon+20) - rand(1,9*no)*2*smallhexagon;
point_y2 = (smallhexagon+20) - rand(1,9*no)*2*smallhexagon;       

inside = inpolygon(point_x2,point_y2,shex_vertices_x2,shex_vertices_y2);

point_x2 = point_x2(inside);
point_y2 = point_y2(inside);

idx2 = randperm(length(point_x2));

point_x2 = point_x2(idx2(1:no));
point_y2 = point_y2(idx2(1:no));

plot(point_x2,point_y2,'ro','MarkerSize',1.5,'LineWidth',1, ...
'MarkerFaceColor','r');

%---Plot points in big hexagon
no2 = 4;
point_x = (bighexagon+20) - rand(1,9*no2)*2*bighexagon;
point_y = (bighexagon+20) - rand(1,9*no2)*2*bighexagon;

inside2 = inpolygon(point_x,point_y,bhex_x,bhex_y);

point_x = point_x(inside2);
point_y = point_y(inside2);

idx = randperm(length(point_x));

point_x = point_x(idx(1:no2));
point_y = point_y(idx(1:no2));

plot(point_x,point_y,'g^','MarkerSize',3,'LineWidth',3, ...
'MarkerFaceColor','g');

Answer 1

Pls check by adding following two lines after inside2 = inpolygon(point_x,point_y,bhex_x,bhex_y);

in1 = inpolygon(point_x,point_y,shex_vertices_x2,shex_vertices_y2);
inside2= logical(inside2-in1);

Answer 2

UPDATE

Thank to Hoki's suggestion, I finally able to work it out.

Note that I change this part inside the code:

validpoint = inpolygon(point_x,point_y,bhex_x,bhex_y) & ~inpolygon(point_x,point_y,shex_vertices_x2,shex_vertices_y2);

Hopefully, this clear the confusions and would help other users as well. I would like to thank Hoki and xenoclast for their help.

The code is as follows:

clear
clc

bighexagon = 20;
smallhexagon = 4;

axis_min = 0;
axis_max = 40; 
axis([axis_min axis_max axis_min axis_max],'square');
hold on

L = linspace(30,390,7); 
bhex_x = bighexagon * (1+cosd(L))'; 
bhex_y = bighexagon*(1+sind(L))';

L2 = linspace(30,390,7); 
shex_x = smallhexagon * (1+cosd(L2))'; 
shex_y = smallhexagon * (1+sind(L2))';

plot(bhex_x,bhex_y,'LineWidth',3);

%---Move small hexagon into big hexagon
shex_vertices_x2(:,1) = shex_x + 16;
shex_vertices_y2(:,1) = shex_y + 16;
plot(shex_vertices_x2(:,1),shex_vertices_y2(:,1),'--k','LineWidth',3);


%---Plot points in small hexagon
no = 3;
point_x2 = (smallhexagon+20) - rand(1,9*no)*2*smallhexagon;
point_y2 = (smallhexagon+20) - rand(1,9*no)*2*smallhexagon;       

inside = inpolygon(point_x2,point_y2,shex_vertices_x2,shex_vertices_y2);

point_x2 = point_x2(inside);
point_y2 = point_y2(inside);

idx2 = randperm(length(point_x2));

point_x2 = point_x2(idx2(1:no));
point_y2 = point_y2(idx2(1:no));

plot(point_x2,point_y2,'ro','MarkerSize',1.5,'LineWidth',1, ...
'MarkerFaceColor','r');

%---Plot points in big hexagon
no2 = 30;
point_x = (bighexagon+20) - rand(1,9*no2)*2*bighexagon;
point_y = (bighexagon+20) - rand(1,9*no2)*2*bighexagon;

%---As per Hoki's suggestion, it ensure the points are outside the small hexagon

validpoint = inpolygon(point_x,point_y,bhex_x,bhex_y) & ...
    ~inpolygon(point_x,point_y,shex_vertices_x2,shex_vertices_y2);

point_x = point_x(validpoint);
point_y = point_y(validpoint);

idx = randperm(length(point_x));

point_x = point_x(idx(1:no2));
point_y = point_y(idx(1:no2));

plot(point_x,point_y,'g^','MarkerSize',3,'LineWidth',3, ...
'MarkerFaceColor','g');

Answer 3

If your inner hexagon is defined by the vertices then you could use inpolygon (link) to test whether a given point is inside it or not.