在球体表面创建一组随机点

Question

我编写了一个代码来在同心球表面上生成一组随机点，在每个球上我会根据我想要的同心球表面的数量生成一定数量的点，如下面编写的代码所示，您可以尝试运行它，也许可以更好地理解手头的问题：

class point:
    def __init__(self, x, y, z):
        self.x = x
        self.y = y
        self.z = z

    R_list=(0.001 , 0.002, 0.003 , 0.004 ,0.005)
    num_points = 1000/len(R_list)
    LD1 = 0.05
    LD23 = 0.025
    obj_points = []
    condition = True
    for R in R_list:
        while len(rpoints) < num_points:
           theta = np.random.uniform(-math.pi, math.pi)
           phi = np.random.uniform(-math.pi, math.pi)
           pi = R * math.cos(theta) * math.sin(phi) + (LD1 / 2), R * math.sin(theta) * math.sin(phi) + (
                    LD23 / 2), R * math.cos(phi) + (LD23 / 2)
           obj_points.append(point(pi[0], pi[1], pi[2]))
    
      # plott array
    lt = len(obj_points)
    Vis_arr = np.zeros((lt, 3))
    n = 0
    for n in range(0, lt-1 ):
        Vis_arr[n][0] = obj_points[n].x
        Vis_arr[n][1] = obj_points[n].y
        Vis_arr[n][2] = obj_points[n].z
    
    # scattering of points
    fig = plt.figure()
    ax = fig.add_subplot(111, projection='3d')
    # scattering of points
    ax.scatter3D(Vis_arr[:, 0], Vis_arr[:, 1], Vis_arr[:, 2], marker="o")
    # graph output fully no evolution
    plt.show()

但是当我运行代码时，一个随机点突然出现：

我怀疑这是创建点数组时的错误，但我无法检测到问题的确切原因。

Answer 1

您的代码中有两个错误：

while len(rpoints) < num_points:

应该是：

while len(obj_points) < num_points:

和

for n in range(0, lt-1 ):

应该是：

for n in range(0, lt ):