# create a renderwindowinteractor
iren = vtk.vtkRenderWindowInteractor()
iren.SetRenderWindow(renWin)
hw_many_planes = 100
hw_many_points = 1000
nanoparticles = []
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)
nanop.get_max_sphere()).get_volume()
#renderer.AddActor(nanop.get_vtk_actor(False, 0.8))
cx, cy, cz = nanop.get_center()
A, B, H = nanop.get_dimensions()
dm = max(B, A, H) / 2.0
maxbox_x = cx + dm
minbox_x = cx - dm
maxbox_y = cy + dm
minbox_y = cy - dm
maxbox_z = cz + dm
minbox_z = cz - dm
cube_actors = cube.cube_to_actors(minbox_x, minbox_y, minbox_z, \
maxbox_x, maxbox_y, maxbox_z, 1.0, 0.0, 0.0)
for a in cube_actors:
renderer.AddActor(a)
"""
p = point.point(cx + dm/2.0, cy-dm, cz+dm/4.0)
renderer.AddActor(p.get_actor(0.5, 1.0, 0.0, 0.0))
pp = nanop.project_point(p)
renderer.AddActor(pp.get_actor(0.5))
dist = nanop.get_distance(p)
renderer.AddActor(p.get_actor(dist))
"""
p = point.point(cx + dm, cy - dm / 2.0, cz + dm / 1.5)
renderer.AddActor(p.get_actor(0.5, 1.0, 0.0, 0.0))
MAX_NUM_OF_CUBE = NXN * NYN * NZN
#MAX_NUM_OF_CUBE = 10*NUM_OF_STARTING_CUBE
print "MAX_NUM_OF_CUBE: ", MAX_NUM_OF_CUBE
print >> sys.stderr, "Box limits: ", xmin, xmax, ymin, ymax, zmin, zmax
# 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.
visactors = False
actors = []
if (visactors):
actors.extend(cube.cube_to_actors (xmin, ymin, zmin, \
xmax, ymax, zmax))
cubcenterx = numpy.empty(0)
cubcentery = numpy.empty(0)
cubcenterz = numpy.empty(0)
cubradius = numpy.empty(0)
centers = []
cubes = []
globalindex = 0
print "First step... ", globalindex, " ", len(cubesets)
j = 0