def slice_traverse(self, viewpoint=None):
if not hasattr(self.ds.index, "grid"):
raise NotImplementedError
for node in kd_traverse(self.tree.trunk, viewpoint=viewpoint):
grid = self.ds.index.grids[node.grid - self._id_offset]
dds = grid.dds
gle = grid.LeftEdge.in_units("code_length").ndarray_view()
nle = get_left_edge(node)
nre = get_right_edge(node)
li = np.rint((nle - gle) / dds).astype('int32')
ri = np.rint((nre - gle) / dds).astype('int32')
dims = (ri - li).astype('int32')
sl = (slice(li[0], ri[0]), slice(li[1],
ri[1]), slice(li[2], ri[2]))
gi = grid.get_global_startindex() + li
yield grid, node, (sl, dims, gi)
def slice_traverse(self, viewpoint = None):
if not hasattr(self.ds.index, "grid"):
raise NotImplementedError
for node in kd_traverse(self.tree.trunk, viewpoint=viewpoint):
grid = self.ds.index.grids[node.grid - self._id_offset]
dds = grid.dds
gle = grid.LeftEdge.in_units("code_length").ndarray_view()
nle = get_left_edge(node)
nre = get_right_edge(node)
li = np.rint((nle-gle)/dds).astype('int32')
ri = np.rint((nre-gle)/dds).astype('int32')
dims = (ri - li).astype('int32')
sl = (slice(li[0], ri[0]),
slice(li[1], ri[1]),
slice(li[2], ri[2]))
gi = grid.get_global_startindex() + li
yield grid, node, (sl, dims, gi)
def traverse(self, viewpoint=None):
for node in kd_traverse(self.tree.trunk, viewpoint=viewpoint):
yield self.get_brick_data(node)
def traverse(self, viewpoint=None):
for node in kd_traverse(self.tree.trunk, viewpoint=viewpoint):
yield self.get_brick_data(node)