def draw_image(self):
w, h = (self.w, self.h)
clearws()
setwindow(0, self.w, 0, self.h)
setviewport(0, 1, 0, 1)
setbackgroundcolor(1, 1, 1, 0)
vertices, normals = triangulate(data, \
(1.0/64, 1.0/64, 1.0/128), (-0.5, -0.5, -0.5), self.isolevel)
mesh = createmesh(len(vertices)*3, vertices, normals, \
ones(vertices.shape))
drawmesh(mesh, 1, (0, 0, 0), (0, 0, 1), (0, 1, 0), (1, 1, 1),
(1, 1, 1))
center = spherical_to_cartesian(-2, pi * self.y / self.h + pi / 2,
pi * self.x / self.w)
up = spherical_to_cartesian(1, pi * self.y / self.h + pi,
pi * self.x / self.w)
cameralookat(center[0], center[1], -0.25 + center[2], 0, 0, -0.25,
up[0], up[1], up[2])
drawimage(0, self.w, 0, self.h, \
self.w, self.h, GR3_Drawable.GR3_DRAWABLE_GKS)
if self.export:
export("mri.html", 800, 800)
print("Saved current isosurface to mri.html")
self.export = False
clear()
deletemesh(c_int(mesh.value))
def draw_image(self):
if not self.needs_refresh:
return
self.needs_refresh = False
w, h = (self.w, self.h)
clearws()
setwindow(0, self.w, 0, self.h)
setviewport(0, 1, 0, 1)
setbackgroundcolor(1, 1, 1, 0)
vertices, normals = triangulate(data, \
(1.0/64, 1.0/64, 1.0/128), (-0.5, -0.5, -0.5), self.isolevel)
mesh = createmesh(len(vertices)*3, vertices, normals, \
ones(vertices.shape))
drawmesh(mesh, 1, (0,0,0), (0,0,1), (0,1,0), (1,1,1), (1,1,1))
center = spherical_to_cartesian(-2, pi*self.y/self.h+pi/2, pi*self.x/self.w)
up = spherical_to_cartesian(1, pi*self.y/self.h+pi, pi*self.x/self.w)
cameralookat(center[0], center[1], -0.25+center[2], 0, 0, -0.25, up[0], up[1], up[2])
drawimage(0, self.w, 0, self.h, \
self.w, self.h, GR3_Drawable.GR3_DRAWABLE_GKS)
if self.export:
export("mri.html", 800, 800)
print("Saved current isosurface to mri.html")
self.export = False
clear()
deletemesh(c_int(mesh.value))
def draw_image(self):
if not self.needs_refresh:
return
self.needs_refresh = False
gr.clearws()
gr.setwindow(0, self.w, 0, self.h)
gr.setviewport(0, 1, 0, 1)
gr3.setbackgroundcolor(1, 1, 1, 0)
vertices, normals = gr3.triangulate(data,
(1.0 / 64, 1.0 / 64, 1.0 / 128),
(-0.5, -0.5, -0.5), self.isolevel)
mesh = gr3.createmesh(
len(vertices) * 3, vertices, normals, np.ones(vertices.shape))
gr3.drawmesh(mesh, 1, (0, 0, 0), (0, 0, 1), (0, 1, 0), (1, 1, 1),
(1, 1, 1))
center = spherical_to_cartesian(-2,
np.pi * self.y / self.h + np.pi / 2,
np.pi * self.x / self.w)
up = spherical_to_cartesian(1, np.pi * self.y / self.h + np.pi,
np.pi * self.x / self.w)
gr3.cameralookat(center[0], center[1], -0.25 + center[2], 0, 0, -0.25,
up[0], up[1], up[2])
gr3.drawimage(0, self.w, 0, self.h, self.w, self.h,
gr3.GR3_Drawable.GR3_DRAWABLE_GKS)
if self.export:
gr3.export("mri.html", 800, 800)
print("Saved current isosurface to mri.html")
self.export = False
gr3.clear()
gr3.deletemesh(c_int(mesh.value))
def cavity_triangles(cavity_grid, cavity_indices, isolevel, step, offset,
discretization_grid):
cavity_triangles = []
cavity_surface_areas = []
grid = np.zeros(cavity_grid.shape, dtype=np.bool)
for cavity_index in cavity_indices:
grid = np.logical_or(grid, cavity_grid == -(cavity_index + 1))
views = []
for x, y, z in itertools.product(*map(xrange, (3, 3, 3))):
view = grid[x:grid.shape[0] - 2 + x, y:grid.shape[1] - 2 + y,
z:grid.shape[2] - 2 + z]
views.append(view)
grid = np.zeros(grid.shape, np.uint16)
grid[:, :, :] = 0
grid[1:-1, 1:-1, 1:-1] = sum(views) + 100
vertices, normals = triangulate(grid, (1, 1, 1), (0, 0, 0), 100 + isolevel)
discrete_vertices = np.array(np.floor(vertices + 0.5), dtype=np.int)
vertices *= np.tile(step, (vertices.shape[0], 3, 1))
vertices += np.tile(offset, (vertices.shape[0], 3, 1))
normals /= np.tile(step, (normals.shape[0], 3, 1))
cavity_surface_area = 0
for cavity_triangle, discrete_triangle in itertools.izip(
vertices, discrete_vertices):
any_outside = False
for vertex, discrete_vertex in itertools.izip(cavity_triangle,
discrete_triangle):
if discretization_grid[tuple(discrete_vertex)] != 0:
any_outside = True
break
if not any_outside:
v1, v2, v3 = cavity_triangle
a = v2 - v1
b = v3 - v1
triangle_surface_area = np.linalg.norm(np.cross(a, b)) * 0.5
cavity_surface_area += triangle_surface_area
return vertices, normals, cavity_surface_area
def cavity_triangles(cavity_grid,
cavity_indices,
isolevel, step, offset,
discretization_grid):
cavity_triangles = []
cavity_surface_areas = []
grid = np.zeros(cavity_grid.shape, dtype=np.bool)
for cavity_index in cavity_indices:
grid = np.logical_or(grid, cavity_grid == -(cavity_index + 1))
views = []
for x, y, z in itertools.product(*map(xrange, (3, 3, 3))):
view = grid[x:grid.shape[0] - 2 + x, y:grid.shape[1] - 2 + y, z:grid.shape[2] - 2 + z]
views.append(view)
grid = np.zeros(grid.shape, np.uint16)
grid[:, :, :] = 0
grid[1:-1, 1:-1, 1:-1] = sum(views) + 100
vertices, normals = triangulate(grid, (1, 1, 1), (0, 0, 0), 100 + isolevel)
discrete_vertices = np.array(np.floor(vertices + 0.5), dtype=np.int)
vertices *= np.tile(step, (vertices.shape[0], 3, 1))
vertices += np.tile(offset, (vertices.shape[0], 3, 1))
normals /= np.tile(step, (normals.shape[0], 3, 1))
cavity_surface_area = 0
for cavity_triangle, discrete_triangle in itertools.izip(vertices, discrete_vertices):
any_outside = False
for vertex, discrete_vertex in itertools.izip(cavity_triangle, discrete_triangle):
if discretization_grid[tuple(discrete_vertex)] != 0:
any_outside = True
break
if not any_outside:
v1, v2, v3 = cavity_triangle
a = v2 - v1
b = v3 - v1
triangle_surface_area = np.linalg.norm(np.cross(a, b)) * 0.5
cavity_surface_area += triangle_surface_area
return vertices, normals, cavity_surface_area
