import numpy as np
import pyclesperanto_prototype as cle
import matplotlib.pyplot as plt
membranes = cle.imread("../../data/membranes.tif")
membranes
cle._ image
|
labels = cle.imread("../../data/membranes_labeled.tif").astype(np.uint32)
labels
cle._ image
|
for r in range(20):
print("radius =", r)
fig, ax = plt.subplots(1,3, figsize=(10,10))
cle.imshow(labels, labels=True, plot=ax[0])
cle.imshow(cle.opening_labels(labels, radius=r), labels=True, plot=ax[1])
cle.imshow(cle.smooth_labels(labels, radius=r), labels=True, plot=ax[2])
ax[0].set_title("Original")
ax[1].set_title("Opening r=" + str(r))
ax[2].set_title("????ing r=" + str(r))
plt.savefig("temp/" + str(r).zfill(2) + ".tif", dpi=300)
plt.show()
radius = 0
radius = 1
radius = 2
radius = 3
radius = 4
radius = 5
radius = 6
radius = 7
radius = 8
radius = 9
radius = 10
radius = 11
radius = 12
radius = 13
radius = 14
radius = 15
radius = 16
radius = 17
radius = 18
radius = 19