"""Interactively cut a set of isosurfaces from a volumetric dataset"""
from vedo import dataurl, show, BoxCutter, Volume

# generate an isosurface the volume for each value
values = [0.1, 0.25, 0.4, 0.6, 0.75, 0.9]
isos = Volume(dataurl+'quadric.vti').isosurface(values) # Mesh

plt = show(isos, __doc__, axes=1, interactive=False)

cutter = BoxCutter(isos)
plt.add(cutter)

plt.interactive()
plt.close()

cmesh = cutter.get_cut_mesh(invert=False).clean()
print("Get the cut mesh:\n", cmesh)