| 1234567891011121314151617181920 |
- """Probe a Volume (voxel dataset) with lines"""
- from vedo import *
- vol = Volume(dataurl + "embryo.slc")
- lines = []
- for i in range(60): # probe scalars on 60 parallel lines
- step = (i - 30) * 2
- p1 = vol.center() + [-100, step, step]
- p2 = vol.center() + [ 100, step, step]
- ln = Line(p1, p2, res=100)
- lines.append(ln)
- lines = merge(lines)
- # Probe the Volume with the lines and add the scalars as pointdata
- lines.probe(vol)
- lines.lw(3).cmap('hot', vmin=0, vmax=110).add_scalarbar()
- print(lines)
- show(vol, lines, __doc__, axes=1).close()
|
