Mayavi Points3D 绘图动画不会更新下一帧

Question

我正在写一个可视化的四个字形通过空间移动的代码。目前，mayavi窗口显示了初始位置，但不会更新到下一个位置。

    #Library decleration 
import numpy as np
from mayavi.mlab import *
....
#Inputting the intital positions into the storage vector
storage_position[0].append(v_1.theta)
storage_position[1].append(v_1.phi)
#Calculating the rest of the positions using the symmetry given
storage_position = Sym(storage_position)
#Plotting the intitial positions 
x_coord = x_trans( storage_position)
y_coord = y_trans(storage_position)
z_coord = z_trans( storage_position)

plt = points3d(x_coord, y_coord, z_coord)
msplt = plt.mlab_source
@mlab.animate(delay=100)
def anim(storage_position, storage_vort, no_vort ,x_coord, y_coord, z_coord):
    f = mlab.gcf()
    while True:
    #for i in range(10):      
        #Working out the hamiltonian
        #Hami(storage_position, storage_vort, 1 - 1, no_vort-1)
        transfer_vector = method(storage_position, storage_vort, 1 - 1, no_vort-1)
        storage_position[0].append(transfer_vector[0])
        storage_position[1].append(transfer_vector[1])
        storage_position = Sym(storage_position)

        x_coord = x_trans( storage_position)
        y_coord = y_trans(storage_position)
        z_coord = z_trans( storage_position)
        msplt.set(x_coord = x_coord, y_coord = y_coord, z_coord = z_coord)
        yield

anim(storage_position, storage_vort, no_vort - 1,x_coord, y_coord, z_coord)
mlab.show()

x_coord等是numpy向量，用于存储四个符号的x坐标。X_trans等是用动画每一步的新坐标更新每个向量的函数。

Answer 1

尝试替换这一行:

msplt.set(x_coord = x_coord, y_coord = y_coord, z_coord = z_coord)

与这个:

msplt.set(x=x_coord, y=y_coord, z=z_coord)

点集数据源(应该是类型为MGlyphSource的对象)知道x、y和z，但不支持x_coord、y_coord和z_coord作为关键字参数。

下面是一个完整的工作示例:

from mayavi import mlab
from numpy import array, cos, sin, cos
x_coord = array([0.0, 1.0, 0.0, -1.0])
y_coord = array([1.0, 0.0, -1.0, 0.0])
z_coord = array([0.2, -0.2, 0.2, -0.2])
plt = mlab.points3d(x_coord, y_coord, z_coord)
msplt = plt.mlab_source
@mlab.animate(delay=100)
def anim():
    angle = 0.0
    while True:
        x_coord = array([sin(angle), cos(angle), -sin(angle), -cos(angle)])
        y_coord = array([cos(angle), -sin(angle), -cos(angle), sin(angle)])
        msplt.set(x=x_coord, y=y_coord)
        yield
        angle += 0.1
anim()
mlab.show()

注意，当您运行此命令时，您可能会看到大量警告打印到控制台上。忽略它们:它们是无害的，而且是已知的问题。