使用mayavi库
import mayavi.mlab
import numpy as np
data = np.loadtxt('points9_label.txt', dtype=np.float64)
x = data[:, 0]
y = data[:, 1]
z = data[:, 2]
max = np.max(data[:, 3])
for i in data:
i[3] = float(i[3])/max
label = data[:, 3]
nodes = mayavi.mlab.points3d(x, y, z, mode="point")
nodes.glyph.scale_mode = 'scale_by_vector'
nodes.mlab_source.dataset.point_data.scalars = label
mayavi.mlab.show()