isosurfaces1.py 493 B

123456789101112131415161718
  1. """Interactively cut a set of isosurfaces from a volumetric dataset"""
  2. from vedo import dataurl, show, BoxCutter, Volume
  3. # generate an isosurface the volume for each value
  4. values = [0.1, 0.25, 0.4, 0.6, 0.75, 0.9]
  5. isos = Volume(dataurl+'quadric.vti').isosurface(values) # Mesh
  6. plt = show(isos, __doc__, axes=1, interactive=False)
  7. cutter = BoxCutter(isos)
  8. plt.add(cutter)
  9. plt.interactive()
  10. plt.close()
  11. cmesh = cutter.get_cut_mesh(invert=False).clean()
  12. print("Get the cut mesh:\n", cmesh)