reduce_peaks

Broad or elongated peaks are replaced with single peaks located at the centre of original voxels

Import packages

import matplotlib.pyplot as plt
import numpy as np
from edt import edt

import porespy as ps

ps.visualization.set_mpl_style()
np.random.seed(0)

peaks

Find peaks on the test image

im = ps.generators.blobs(shape=[200, 200])
dt = edt(im)
peaks = ps.filters.find_peaks(dt, r_max=4)

plt.figure(figsize=[6, 6])
plt.axis(False)
plt.imshow(peaks / im);

Apply reduce_peaks filter

x = ps.filters.reduce_peaks(peaks=peaks)
plt.figure(figsize=[6, 6])
plt.axis(False)
plt.imshow(x / im);