insert_cylinder

Import packages

import matplotlib.pyplot as plt

import porespy as ps

ps.visualization.set_mpl_style()

Generate image for testing

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

fig, ax = plt.subplots(figsize=[4, 4])
ax.imshow(im[..., 50])
plt.axis(False);

Demonstrate tool

im = ps.tools.insert_cylinder(im=im * 1.0, xyz0=[20, 20, 20], xyz1=[60, 60, 50], r=20)

fig, ax = plt.subplots(figsize=[4, 4])
ax.imshow(ps.visualization.xray(~(im > 1)))
plt.axis(False);