flood#

Floods each region with a specific value based on a specified statistical operation performed on values in that region.

Import packages#

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

im#

The distance transform can have statistical calculations performed

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

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

labels#

snow_partitioning can be used to create regions

regions = ps.filters.snow_partitioning(im, r_max=4, sigma=0.4)
labels = regions.regions

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

mode#

Various functions in scipy.ndimage.measurements are called to perform statistical calculation. The mode indicates which function to call.

x1 = ps.filters.flood(im=dt, labels=labels, mode='max')
x2 = ps.filters.flood(im=dt, labels=labels, mode='mean')
x3 = ps.filters.flood(im=dt, labels=labels, mode='sum')

fig, ax = plt.subplots(1, 3, figsize=[18, 18])
ax[0].imshow(x1)
ax[0].axis(False)
ax[0].set_title('mode = max')
ax[1].imshow(x2)
ax[1].axis(False)
ax[1].set_title('mode = mean')
ax[2].imshow(x3)
ax[2].axis(False)
ax[2].set_title('mode = sum');
../../../_images/bda9012fe88890e69a613fed175d4215df17761324a1575c760056b4ac956d90.svg