Bonjour ! 

J'ai un champ vectoriel en 3D (xx,yy,zz) avec les coordonnées (X,Y,Z) crées par un meshgrid. J'aimerais pouvoir slicer mes données au milieu de chaque plan (xy,xz,zy). Voici ce que j'ai essayé pour xy :

points = 3
ext = 1
xrange = np.linspace(-ext,ext,points)
yrange = np.linspace(-ext*2,ext*2,points)
zrange = np.linspace(-ext*4,ext*4,points)
X,Y,Z = np.meshgrid(xrange, yrange, zrange)

xx = np.ones_like(X)
yy = np.ones_like(X)
zz = np.ones_like(X)

fig, ax = plt.subplots(1,figsize=(8,8))


# plane xy
ax.quiver(X[:,:,1],Y[:,:,1], 
              xx[:,:,1],yy[:,:,1])
ax.set_xlabel("x")
ax.set_ylabel("y")

De la même manière pour zy :

# plane zy
ax.quiver(Z[1,:,:],Y[1,:,:], 
              zz[1,:,:],yy[1,:,:])
ax.set_xlabel("z")
ax.set_ylabel("y")

et pour xz : 

# plane xz
ax.quiver(X[:,1,:],Z[:,1,:], 
              xx[:,1,:],zz[:,1,:])
ax.set_xlabel("x")
ax.set_ylabel("z")

Il semblerait que mon slice n'est pas correcte parce que le plot ne montre pas correctement les plans recherchées. Auriez- vous s'il-vous-plaît une idée de comment améliorer cela ? Merci beaucoup !

résolu? pas résolu ? si résolu, comment ?