"""Lookup table mapping cube index to its deleteability"""
keep = np.logical_not(delete)
"""Lookup table mapping cube index to its non-deleteability"""
filename_non_removable = "PK12nr.npy"
"""Filename for the lookup table mapping a cube configuration to the non-removeability of the center pixel"""
non_removable = initialize_lookup_table(filename=filename_non_removable,
function=match_non_removable)
"""Lookup table mapping cube index to its non-removeability"""
consider = np.logical_not(non_removable)
"""Lookup table mapping cube index to whether it needs to be considered further"""
rotations = t3d.rotations12(base)
"""Rotations of the base cube for the sub-iterations"""
###############################################################################
### Skeletonization
###############################################################################
def skeletonize(binary,
points=None,
steps=None,
removals=False,
radii=False,
check_border=True,
delete_border=False,
return_points=False,
def _test():
import numpy as np
import ClearMap.ImageProcessing.Topology.Topology3d as top
from importlib import reload
reload(top)
label = top.cube_labeled()
top.print_cube(label)
# Test rotations
c = np.zeros((3, 3, 3), dtype=bool)
c[1, 0, 0] = True
top.print_cube(c)
cs = [top.rotate(c, axis=2, steps=r) for r in range(4)]
[top.print_cube(cc) for cc in cs]
reload(top)
l = top.cube_labeled()
rts = top.rotations6(l)
[top.print_cube(r) for r in rts]
reload(top)
b = top.cube_from_index(6)
i = top.cube_to_index(b)
print(i, 6)
us = np.zeros((3, 3, 3), dtype=int)
us[1, 1, 2] = 1
us[1, 0, 1] = 1
us[1, 2, 0] = 2
r12 = top.rotations12(us)
[top.print_cube(cc) for cc in r12]
#check configuration utlity
reload(top)
index = 11607
source = top.cube_from_index(index)
c = top.index_from_binary(source)
c[1, 1, 1] == index
x = np.random.rand(1500, 500, 500) > 0.6
c = top.index_from_binary(x)
import numpy as np
import ClearMap.ImageProcessing.Topology.Topology3d as top
#check fortran vs c order
x = np.random.rand(5, 5, 5) > 0.35
y = np.asanyarray(x, order='F')
ix = top.index_from_binary(x)
iy = top.index_from_binary(y)
ax = ix.array
ay = iy.array
#%% profile
import io
io.DEFAULT_BUFFER_SIZE = 2**32
import pstats, cProfile
import numpy as np
import ClearMap.ImageProcessing.Topology.Topology3d as top
x = np.ones((3000, 500, 1000), dtype=bool, order='F')
import ClearMap.IO.IO as io
import ClearMap.ParallelProcessing.DataProcessing.ArrayProcessing as ap
ap.write('test.npy', x)
y = io.as_source('test.npy')
z = io.create('resuly.npy', shape=y.shape, order='C', dtype='uint32')
cProfile.runctx(
"c =top.index_from_binary(y, method='!shared', sink=z, verbose=True, processes=None)",
globals(), locals(), "Profile.prof")
s = pstats.Stats("Profile.prof")
s.strip_dirs().sort_stats("time").print_stats()
import mmap
mmap.ACCESS_COPY