Determine if a point is inside a cube

Question

Lulu 2011-11-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/21429-determine-if-a-point-is-inside-a-cube

回答： sajjad nasiri 2018-7-18

Hi! There are given 8 vertices of a cube and a random point. How can I define that the point lies inside the cube? I wanted to create 2 vectors from the random point. These vectors should be parallel to X and Y axis (Z is not needed here). Then, if vectors intersect with cube sides 4 times, then the point lies inside a cube. But this seems to be too complex approach for this problem. Does someone know a simpler way? A code snipped would be very helpful. Thanks.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Honglei Chen 2011-11-16

Could you give a little more information regarding the issue? Is the cube really a right-angle equal-side-length cube? Is the orientation arbitrary or it has to be aligned with axes? Maybe an example can clarify these.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Image Analyst 2011-11-16

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/21429-determine-if-a-point-is-inside-a-cube#answer_28222

The way I'd do it is to call convhulln(). Basically it gives you the vertices of your set of 3D vertices that are on the convex hull. The convex hull is what you'd get if you basically wrapped a balloon around your points.

So all you have to do is to pass in your cube coordinates concatenated with your "test" point and see if it returns the cube coordinates. If it doesn't then the point is outside. If it does, then the cube is the convex hull and that must mean your test point is inside the cube. It's one line of code to call convhulln() and another to do the check to see if what it returns matches your cube coordinates. You should be able to do it. If you can't then it's your turn to supply us with code that generates cube coordinates and test points (one inside and one outside) so that we can do those two lines for you.

2 个评论
显示无隐藏无

Lulu 2011-11-16

Thanks. I wrote the following code according to your idea:

% Create a cube

xyz = [0 0 0; 1 1 1; 0 0 1; 0 1 0; 1 0 0; 1 0 1; 1 1 0; 0 1 1];

random_point = [1 2 1]; % lies outside the cube

xyz = [xyz; random_point];

K = convhulln(xyz);

len = length(xyz(:,1));

inside = any(any(K==len))

Please let me know if there is a better implementation of your idea.

Image Analyst 2011-11-17

Well K looks more complicated than the 2D case. I'd have to think about it. But you'd have to call convhulln twice (once for cube and once for cube + extra points) and use ismember to compare that each row of K1 is included in K2.

请先登录，再进行评论。