Produces a set of images with each panel containing one saturation. This method can be applied for visualizing image-based invasion percolation algorithm ibip filter.

import matplotlib.pyplot as plt
import numpy as np
import porespy as ps
import inspect
<Signature (satn, im, bins=None, axis=0, slice=None, **kwargs)>


The input image is a Boolean image True values indicating the void voxels and False for solid. Let’s create a test image:

im = ps.generators.blobs(shape=[50,50,50])
fig, ax = plt.subplots(1, 1, figsize=[6, 6])
ax.imshow(im[:,:,20], origin='lower', interpolation='none')


The saturation image can be generated from ibip data using seq_to_satn method. The satn is the image of porous material where each voxel indicates the global saturation at which it was invaded. Voxels with 0 values indicate solid and and -1 indicate uninvaded. As the image is large, we will be visualizing a section of the image: Note that the regions that were uninvaded at sat=0.7 are now invaded at sat=0.75 and remaining regions are invaded at sat=1:

out = ps.filters.ibip(im=im)
inv_seq, inv_size = out.inv_sequence, out.inv_sizes
satn = ps.filters.seq_to_satn(seq=inv_seq)
fig, ax =ps.visualization.satn_to_panels(satn[:15,:15,1], im[:15,:15,1]);


Indicates for which saturations images should be made. By default all saturation values in the image are used. To visualize the image for a list of equally space values between 0 and 1, an int value can be passed as the input (the number of saturation points between [0,1]).

fig, ax =ps.visualization.satn_to_panels(satn, im, bins= 5);


If the image is 2D this variable is ignored. If the image is 3D, a 2D image is extracted at the specified slice taken along this axis. By default axis=0 indicating the slice is extracted at x axis. Let’s extract the slices at y axis (because slice is not passed as the input, the default slice is at the mid-point of the axis):

fig, ax =ps.visualization.satn_to_panels(satn, im, bins= 5, axis=1);


If the image is 2D this variable is ignored. If the image is 3D, a 2D image is extracted from this slice along the given axis. By default a slice at the mid-point of the axis is returned. Let’s extract a slice at y=4:

fig, ax =ps.visualization.satn_to_panels(satn, im, bins= 5, axis=1, slice=4);

Note that extracting a slice at y=25 gives the same output as we have seen in the axis example above (default slice at mid-point axis):

fig, ax =ps.visualization.satn_to_panels(satn, im, bins= 5, axis=1, slice=25);


Additional keyword arguments can be sent to the imshow function, such as interpolation.

fig, ax =ps.visualization.satn_to_panels(satn, im, bins= 5, axis=1, slice=25, interpolation=None);