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

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);

Apply region_size filter

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

plt.figure(figsize=[6, 6])
plt.axis(False)
plt.imshow(x);