python : fill cuboid with spheres

Question

我想知道如何在长方体中创建球体。这个长方体 (50,50,50) 包含有很多点，我想要的只是获取球体(半径 = 2)内的点位置及其值。

非常感谢任何想法/建议/代码/信息..

Answer 1

创建一个以其中心 (i,j,k) 和长方体 (x0,y0,z0) 为中心的球体:

points_in_sphere = (x0 - i)**2 + (y0 - j)**2 + (z0 - k)**2 <= (radius)**2

提取点的坐标:

x1,y1,z1 = np.argwhere(points_sphere)[:,0], np.argwhere(points_sphere)[:,1], np.argwhere(points_sphere)[:,2]

打印点值:

print cubiod[x1,y1,z1]