import numpy as np
from mayavi import mlab
x, y, z = np.ogrid[-10:10:20j,-10:10:20j,-10:10:20j]
s = np.sin(x*y*z)/(x*y*z)
mlab.contour3d(s)
mlab.show()
切平面
import numpy as np
from mayavi import mlab
x, y, z = np.ogrid[-10:10:20j,-10:10:20j,-10:10:20j]
s = np.sin(x*y*z)/(x*y*z)
mlab.pipeline.image_plane_widget(mlab.pipeline.scalar_field(s),# 数据的标量数据场
plane_orientation='x_axes',# 切平面的方向
slice_index=10)
mlab.pipeline.image_plane_widget(mlab.pipeline.scalar_field(s),
plane_orientation='y_axes',
slice_index=10)
mlab.outline()
mlab.show()
复合观测
import numpy as np
from mayavi import mlab
x, y, z = np.ogrid[-10:10:20j,-10:10:20j,-10:10:20j]
s = np.sin(x * y * z)/(x * y * z)
src = mlab.pipeline.scalar_field(s)
mlab.pipeline.iso_surface(src, contours=[s.min()+0.1* s