vedo/examples/volumetric/vol2points.py

"""Extract all image voxels as points"""
from vedo import *

v = Volume(dataurl+'vase.vti')

pts = v.topoints().print() # returns Points

scalars = pts.pointdata[0]
pts.cmap('afmhot_r', scalars).point_size(1)

show([(v,__doc__), pts], N=2, viewup='z', bg2='lightblue', axes=1).close()