if (len(sys.argv)) == 2: filename = sys.argv[1] nanoparticles = [] botx, topx, boty, topy, botz, topz = \ nanoparticle.file_to_nanoparticle_list(filename, nanoparticles) if (botx >= topx) or (boty >= topy) or \ (boty >= topy): print "Error Invalid BOX" exit() scx, scy, scz, radius = \ nanoparticle.nanoparticle_list_to_arrays(nanoparticles) renderer = vtk.vtkRenderer() renWin = vtk.vtkRenderWindow() renWin.AddRenderer(renderer) iren = vtk.vtkRenderWindowInteractor() iren.SetRenderWindow(renWin) for selected_index in range(0, 1): nanop = nanoparticles[selected_index] renderer.AddActor(nanop.get_vtk_actor(True, 1.0)) for p in nanop.inside_point_grid(): renderer.AddActor(p.get_actor(0.1, rc=1.0, gc=0.0, bc=0.0)) renderer.AddActor(nanop.get_vtk_actor(opacity=0.5))
zmax = xmax = ymax = -10000.0 zmin = xmin = ymin = 10000.0 xmin, xmax, ymin, ymax, zmin, zmax = \ nanoparticle.file_to_nanoparticle_list(filename, nanoparticles) cube_actors = cube.cube_to_actors(xmin, ymin, zmin, \ xmax, ymax, zmax, 1.0, 1.0, 1.0) for a in cube_actors: ren.AddActor(a) Dz = (zmax - zmin) Dy = (ymax - ymin) Dx = (xmax - xmin) scx, scy, scz, radius = nanoparticle.nanoparticle_list_to_arrays(nanoparticles) meanr = radius.mean() meand = 2.0 * meanr # voglio fermarmi a circa 2 D dalla vetta visto che in cima avro' sempre una # densita' minore(ricorda la prima sfera che supera zmax ferma la procedura. dz = ((zmax - zmin) - (3.0 * meand)) / float(hw_many_planes + 1) poly_data_points = [] center = [] for iplane in range(hw_many_planes): zplane = zmin + meand + (iplane + 1) * dz