1234567891011121314151617 |
- """Cut a mesh with another mesh"""
- from vedo import dataurl, settings, Plotter, Volume, Ellipsoid
- settings.tiff_orientation_type = 4 # data origin is bottom-left
- vol = Volume(dataurl + "embryo.tif")
- iso = vol.isosurface(30, flying_edges=False).normalize().pos(0,0,0)
- emsh = Ellipsoid().scale(0.4).pos(2.8, 1.5, 1.5).wireframe()
- # make a working copy and cut it with the ellipsoid
- cut_iso = iso.clone().cut_with_mesh(emsh).c("gold").bc("t")
- plt = Plotter(N=2, axes=1)
- plt.at(0).show(iso, emsh, __doc__)
- plt.at(1).show(cut_iso, viewup="z")
- plt.interactive().close()
|