快速球体绘图与python

Question

好的，我在stackoverflow中的第一篇文章。我必须在同一图中绘制具有不同半径的30000个球体。我尝试了以下代码，它使用了matplotlib.pyplot

def plot_sphere(ax, x, y, z , r , resolution=30, **kwargs):

    u = np.linspace(0, 2 * np.pi, resolution)
    v = np.linspace(0, np.pi, resolution)

    xx = r * np.outer(np.cos(u), np.sin(v)) + x
    yy = r * np.outer(np.sin(u), np.sin(v)) + y
    zz = r * np.outer(np.ones(np.size(u)), np.cos(v)) + z

    ax.plot_surface(xx, yy, zz, rstride=1, cstride=1 , linewidth=0, **kwargs)

然而，它需要大量内存并且需要数小时才能绘制。有没有办法优化这个？我也对其他可能使用mayavi2或opengl等其他库的解决方案持开放态度。谢谢！