首页
学习
活动
专区
圈层
工具
发布
社区首页 >问答首页 >如何在MATLAB中绘制特定多边形区域外的点

如何在MATLAB中绘制特定多边形区域外的点
EN

Stack Overflow用户
提问于 2015-05-04 13:20:49
回答 3查看 472关注 0票数 2

我一直在寻找一种方法来绘制多边形区域外的点(在我的例子中是六边形)。这是我想要实现的场景,我在一个大六边形里面有一个小的六边形。情况如下:

在图片中,我创建了一个小的六边形(用浅红色表示的区域),并使用inpolygon在其中生成一个随机点(在我的例子中是三个)。当我想在大六边形(用浅紫色表示)中绘制点(红色三角形)时出现问题,而不接触小六边形区域。我已经环顾网络寻找这个简单的解决方案3天,但没有任何效果。

我真的很感激我能得到的任何帮助或指导。非常感谢!

我的代码如下:

代码语言:javascript
复制
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');
EN

回答 3

Stack Overflow用户

回答已采纳

发布于 2015-05-04 15:53:33

更新

感谢Hoki的建议,我终于想出了办法。

注意:我在代码中更改了这个部分:

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

希望这能消除困惑,也能帮助其他用户。我要感谢Hokixenoclast的帮助。

守则如下:

代码语言:javascript
复制
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');
票数 2
EN

Stack Overflow用户

发布于 2015-05-04 13:26:28

如果内六边形是由顶点定义的,那么您可以使用inpolygon (链接)来测试给定点是否在其中。

票数 1
EN

Stack Overflow用户

发布于 2015-05-04 16:09:10

请在inside2 = inpolygon(point_x,point_y,bhex_x,bhex_y);后添加以下两行检查

代码语言:javascript
复制
in1 = inpolygon(point_x,point_y,shex_vertices_x2,shex_vertices_y2);
inside2= logical(inside2-in1);
票数 1
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/30031241

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档