I want to find out points whose coordinates are within the volume formed by a cone shape.

22 次查看（过去 30 天）

显示更早的评论

Tamoor Shafique 2022-4-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1706600-i-want-to-find-out-points-whose-coordinates-are-within-the-volume-formed-by-a-cone-shape

编辑： Matt J 2022-4-28

采纳的回答： Bruno Luong

Hello,

I have plotted few points in a 3D volume as shown in the figure below:

I want to plot a 3D righ Cone from a few specefic points as shown:

These cones are 3D volumes and I want to find out points whose coordinates are within each cone. for example which point is within yello cone and which one in red and which one in purple cone.

Can some one help?

采纳的回答

Bruno Luong 2022-4-28

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1706600-i-want-to-find-out-points-whose-coordinates-are-within-the-volume-formed-by-a-cone-shape#answer_952710

编辑：Bruno Luong 2022-4-28

在 MATLAB Online 中打开

Assuming

the cone center is C: (3 x 1)
coneangle is the cone angle (0 to pi, pi corresponds to the extreme case where the cone takes the whole space)
the cone axis is N: (3 x 1), oriented toward the direction you consider to be in
Point P: (3 x 1) the point coordinates to be check

Check inside the half of the cone is

dot(P-C, N) >= cos(coneangle/2) * norm(N) * norm(P-C)

You can compute once

cos(coneangle/2) * norm(N)

if you have many points P to check.

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

Bruno Luong 2022-4-28

编辑：Bruno Luong 2022-4-28

The above code plots one part of the surface of the cone truncated at base hight = h, base radius r.

Why? Because the author wants it that way.

Tamoor Shafique 2022-4-28

I got it. sorry for being so dumb.

Thank you very much and much appreciated.

I have resolved this problem

请先登录，再进行评论。

更多回答（1 个）

Matt J 2022-4-27

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1706600-i-want-to-find-out-points-whose-coordinates-are-within-the-volume-formed-by-a-cone-shape#answer_952385

Draw a vector from the cone vertex to the point and compute the angle of this vector to the cone's axis. If the angle is less than the angular width of the cone, then the point is inside the cone. This assumes that the cone is not of finite length.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Tamoor Shafique 2022-4-28

在 MATLAB Online 中打开

Sorry one more problem:

Why is following not giving me a 3D Cone but a 2D Cone? How can I convert it into a 3D Cone?

function [X2,Y2,Z2] = cone_plot3(r,h,phi,theta,offset)
m = h/r;
[R,A] = meshgrid(linspace(0,r,50),linspace(0,2*pi,50));
X = R.*cos(A);
Y = R.*sin(A);
Z = m*R;
% Cone around the z-axis, point at the origin
X1 = X*cos(phi) - Z*sin(phi);
Y1 = Y;
Z1 = X*sin(phi) + Z*cos(phi);
% Previous cone, rotated by angle phi about the y-axis
X2 = X1*cos(theta) + offset(1);
Y2 = X1*sin(theta) + offset(2);
Z2 = Z1 + offset(3);
% Second cone rotated by angle theta about the z-axis
end

The Code is

clc
clear all
r=10;
h=30;
phi=pi/3;
theta=pi/4;
offset(1)=10;
offset(2)=10;
offset(3)=10;
[X3 Y3 Z3]=cone_plot3(r,h,phi,theta,offset)
surf(X3, Y3, Z3)