Boundary pores during fibrous mat to network conversion

[kunal0kumar0sinha]

I am trying to convert a fibrous mat into a pore network model to simulate in openPNM. However, I am unable to label the boundary pores. I followed the steps suggested in "adding boundary pores". Code is shown below:

im = ps.generators.cylinders(shape=[200, 200, 200], r=5, ncylinders=200)
sigma = 0.4
dt = spim.distance_transform_edt(input=im)
dt1 = spim.gaussian_filter(input=dt, sigma=sigma)
peaks = find_peaks(dt=dt)
print('Initial number of peaks: ', spim.label(peaks)[1])
peaks = trim_saddle_points(peaks=peaks, dt=dt1)
print('Peaks after trimming saddle points: ', spim.label(peaks)[1])
peaks = trim_nearby_peaks(peaks=peaks, dt=dt)
peaks, N = spim.label(peaks)
print('Peaks after trimming nearby peaks: ', N)
regions = watershed(image=-dt, markers=peaks, mask=dt > 0)
regions = randomize_colors(regions)

After this I tried to add boundary pores

bds = ps.networks.add_boundary_regions(regions=regions*im, pad_width=10)
net = ps.networks.regions_to_network(bds)
pn = op.io.network_from_porespy(net)

However, I cannot locate boundary pores to do the simulation in PNM

Please help