| 123456789101112131415161718192021 |
- """Density plot from a distribution of points in 3D"""
- import numpy as np
- from vedo import *
- settings.use_depth_peeling = True
- n = 3000
- p = np.random.normal(7, 0.3, (n,3))
- p[:int(n*1/3) ] += [1,0,0] # shift 1/3 of the points along x by 1
- p[ int(n*2/3):] += [1.7,0.4,0.2]
- pts = Points(p, alpha=0.5)
- vol = pts.density() # density() returns a Volume
- vol.cmap('Dark2').alpha([0.1,1])
- r = precision(vol.metadata['radius'][0], 2) # retrieve automatic radius value
- vol.add_scalarbar3d(title=f'Density (counts in r_s ={r})', italic=1)
- show(pts, vol, __doc__, axes=True).close()
|
