trim_small_clusters#

trim_small_clusters function is a filter that removes isolated voxels or clusters of a given size or smaller.

import numpy as np
import porespy as ps
import scipy.ndimage as spim
import matplotlib.pyplot as plt
import skimage
ps.visualization.set_mpl_style()

im#

This function works on 2D and 3D images.

im = ps.generators.blobs(shape=[500, 500], blobiness=5, porosity = 0.5)

plt.figure(figsize=[6,6]);
plt.imshow(im);
plt.axis(False);
../../../_images/bb27c72fc95a59dd6fe495c5344f1b65e10de1dd71af54cc713be81545b05737.svg

size#

The maximum size of a cluster to trim is varied trying 10, 100, and 500 voxels.

x1 = ps.filters.trim_small_clusters(im, size=10)
x2 = ps.filters.trim_small_clusters(im, size=100)
x3 = ps.filters.trim_small_clusters(im, size=500)

fig, ax = plt.subplots(1, 3, figsize=[18, 18]);
ax[0].imshow(x1);
ax[0].axis(False);
ax[0].set_title('size = 10', fontdict={'fontsize': 18});
ax[1].imshow(x2);
ax[1].axis(False);
ax[1].set_title('size = 100', fontdict={'fontsize': 18});
ax[2].imshow(x3);
ax[2].axis(False);
ax[2].set_title('size = 500', fontdict={'fontsize': 18});
../../../_images/7b06a91b899da9ee6a6550fc1ec47e89c7f6dcdee818147735e31a6a6e3d5337.svg