这更像是一个几何问题，而不是一个编程问题。对于x，y，z空间中的任何给定点，x^2 + z^2 < (0.15/0.23 * y)^2将位于您定义的圆锥内的点，所以您可以这样做：

data(bunny)

pts <- cbind(c(0, 0), c(0, 0.23), c(0, 0))
radii <- c(0.0, 0.15)
cone <- cylinder3d(pts, radii, sides = 100)

cone_bunny <- bunny[bunny[,1]^2 + bunny[,3]^2 < (0.15/0.23 * bunny[,2])^2,]

points3d(cone_bunny, col="pink", alpha = 0.5)
shade3d(cone, col = "lightblue", alpha = 0.3)

字符串
x1c 0d1x的数据
我不确定你是否在寻找一个更通用的3D形状相交的解决方案-这要困难得多，但许多简单的几何相交可以从简单的代数计算。

增编

在评论中指出，你可能不想要一个实际的圆锥体，而是一个20边的金字塔，如你的例子所示。这也可以通过数学解决方案来解决：

cone <- cylinder3d(pts, radii, sides = 20) 

lims <- (0.15/0.23 * bunny[,2])^2/cos(atan2(bunny[,3], bunny[,1]) %% 
         (2*pi)/20 - pi/20)

cone_bunny <- bunny[bunny[,1]^2 + bunny[,3]^2 < lims,]

points3d(cone_bunny, col="pink", alpha = 0.5)
shade3d(cone, col = "lightblue", alpha = 0.3)

型

的