region_size#

Each voxel is replaced with the size of the region to which it belongs

Import packages#

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

im#

Generate a test image. Can be either 2D or 3D.

im = ps.generators.blobs(shape=[200, 200])

plt.figure(figsize=[6, 6])
plt.axis(False)
plt.imshow(im);
../../../_images/1de079eadab8afabad3cc0984209e6b9c38aa40f15e5cd4501f62583fbaa66ca.png

Apply region_size filter#

x = ps.filters.region_size(im=im)

plt.figure(figsize=[6, 6])
plt.axis(False)
plt.imshow(x);
../../../_images/7017677437696e29ec8c1cecae152a0e5e503bb26bd097896fbdab2f683426d3.png