extract_cylinder#

Import Packages#

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()
[03:26:29] ERROR    PARDISO solver not installed, run `pip install pypardiso`. Otherwise,          _workspace.py:56
                    simulations will be slow. Apple M chips not supported.                                         

Generate Image for Testing#

im = ps.generators.blobs([100, 100, 100])

Apply tool and visualize#

This tool is useful for extracting a cylindrical core from an otherwise cubic image. This might be helpful if a cylindrical domain is needed, or to test a workflow on a tomogram that has has not been cropped yet.

cyl = ps.tools.extract_cylinder(im=im, r=40, axis=2)

fig, ax = plt.subplots(1, 2, figsize=[8, 4]);
ax[0].imshow(im[..., 50]);
ax[0].axis(False)
ax[1].imshow(cyl[..., 50]);
ax[1].axis(False);
../../../_images/b75e77c2bb0600f712d6d4c556701ad8bc9a3971dd8da8140bb0388ecb154e67.png