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

5 次查看（过去 30 天）
Tamoor Shafique
Tamoor Shafique 2022-4-27
编辑： Matt J 2022-4-28
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
Bruno Luong 2022-4-28
编辑：Bruno Luong 2022-4-28
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更早的评论
Bruno Luong
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.

请先登录，再进行评论。

更多回答（1 个）

Matt J
Matt J 2022-4-27
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更早的评论
Tamoor Shafique
Tamoor Shafique 2022-4-28
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)
Any help would be welcomed

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Lighting, Transparency, and Shading 的更多信息

标签

产品


版本

R2021a

Community Treasure Hunt

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

Start Hunting!

Translated by