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