如何均匀地填充由 3d 点列表给出的 convexHull 的体积？

Question

我有一个由 scipy.spatial.ConvexHull 生成的 3d 点列表给出的凸包。

我希望整个体积以某种有效的方式均匀地填充 3d 点。

如果我有一个快速的方法来知道一个点是在凸壳内部还是外部，我可以矢量化以某种分辨率遍历每个体素并返回它的中心，如果它是 "inside" 或什么都不返回对于那个体素，如果它是 "outside".

示例：

对于 2d 点，如 scipy 的 ConvexHull 示例，

我想计算红线内均匀分布的点列表。

如何才能更有效地做到这一点？

Answer 1

这是一些（未优化的）示例代码，演示了 .equations 属性的使用，该属性包含凸包的表面法线。

import numpy as np
from scipy import spatial

## set up example ##
# create sample points
np.random.seed(17)
samples = np.random.uniform(-5,5,(100,2))
# pick convex subset
outline = samples[(samples**2).sum(1)<4]
outline = spatial.ConvexHull(outline)
# choose tolerance for detection of boundary points
eps = 1e-9

## classify ##
outside = (outline.equations@np.c_[samples, np.ones(len(samples))].T > eps).any(0)
inside = (outline.equations@np.c_[samples, np.ones(len(samples))].T < -eps).all(0)
boundary = ~(outside|inside)

## plot ##
import pylab
closed = np.tile(outline.points[outline.vertices],(2,1))
closed = closed[:len(closed)//2+1]
pylab.plot(*closed.T,'b')
pylab.plot(*samples[outside].T,'og')
pylab.plot(*samples[inside].T,'or')
pylab.plot(*samples[boundary].T,'oy')
pylab.show()

如何均匀地填充由 3d 点列表给出的 convexHull 的体积？

How to uniformly fill a volume of a convexHull given by a list of 3d points?

python

numpy

scipy

convex-hull