How to say the third selected point lies within the created geometry or not?

1 次查看(过去 30 天)
close all
clc
base_x = [0 50 50 20 0 0]
base_y = [0 0 10 10 10 0 ]
roof_x = [ 0 0 20 20 30 30]
roof_y =[ 10 50 10 50 50 10]
figure
plot(base_x, base_y,'-')
hold on
plot(roof_x, roof_y,'-')
axis([0 50 0 50])
grid on
[X,Y]=meshgrid(1:11);
[x] = ginput(2);
[x] = round(x);
fprintf ('the a1 coordinate value is %d,\n' , x(1,1));
fprintf ('the b1 coordinate value is %d ,\n' , x(1,2));
fprintf ('the a2 coordinate value is %d,\n' , x(2,1));
fprintf ('the b2 coordinate value is %d ,\n' , x(2,2));
How to say the third selected points a1,b1 a2,b2 lies within the created geometry ( base,roof) or not?
  4 个评论
Adam Danz
Adam Danz 2018-7-25
Kaleesh, I'm continuing your other thread (link above) here because this topic seems to be more relevant to the problems you're having.

请先登录,再进行评论。

采纳的回答

Adam Danz
Adam Danz 2018-7-24
编辑:Adam Danz 2018-7-25
If I understand your question correctly, you're trying to determine if a coordinate (x,y) lies within a polygon.
First try matlab's inpolygon().
Another way to achieve that is by "ray casting". The idea is to draw a line from the (x,y) coordinate and extend it rightward until it passes all of the possible polygon edges. Then count the number of times that ray intersects with an edge of the polygon. If that number is non-zero and odd, then the coordinate (x,y) is within the polygon.
To achieve that, you'll need to loop through each edge of your polygon and determine if each line segment intersects with the virtual ray. Here's a nice tutorial on how to detect line segment intersection. And here's more info on the ray casting technique.
  5 个评论
Ramesh Bala
Ramesh Bala 2018-7-25
编辑:Ramesh Bala 2018-7-25
yeah the bottom edge wasn't mentioned in the previous code,but there in the updated one
xv = [0 0 0 20 20 0 10 20 20 20 20 30 30 50 50 20 30 40 50]
yv = [0 0 50 10 0 0 0 0 0 10 50 50 10 10 0 0 0 0 0 ]
let me try to check with these values now Danke /Thanks Adam
Ramesh Bala
Ramesh Bala 2018-7-26
yeah the command worked ,even the previous one also works. for the further case just added frpintf by calling on && in with values to check

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Line Plots 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by