使用ndimage.maximum_filter和skimage.peak_local_max查找图像峰值

时间:2019-05-15 11:33:27

标签: python scikit-image ndimage

我正在尝试查找给定图像的一些相对最大值。我了解有两种可能的方法,第一种使用scipy.ndimage.maximum_filter(),第二种使用skimage.feature.peak_local_max()

为了比较这两种方法,我从skimage shown here修改了一个示例,以便比较发现的峰。

from scipy import ndimage as ndi
import matplotlib.pyplot as plt
from skimage.feature import peak_local_max
from skimage import data, img_as_float

im = img_as_float(data.coins())

# use ndimage to find the coordinates of maximum peaks
image_max = ndi.maximum_filter(im, size=20) == im

j, i = np.where(image_max)
coordinates_2 = np.array(zip(i,j))

# use skimage to find the coordinates of local maxima
coordinates = peak_local_max(im, min_distance=20)

# display results
fig, axes = plt.subplots(1, 2, figsize=(8, 3), sharex=True, sharey=True)
ax = axes.ravel()

ax[0].imshow(im, cmap=plt.cm.gray)
ax[0].plot(coordinates_2[:, 0], coordinates_2[:, 1], 'r.')
ax[0].axis('off')
ax[0].set_title('Maximum filter')

ax[1].imshow(im, cmap=plt.cm.gray)
ax[1].autoscale(False)
ax[1].plot(coordinates[:, 1], coordinates[:, 0], 'r.')
ax[1].axis('off')
ax[1].set_title('Peak local max')

fig.tight_layout()

plt.show()

这给出了每种方法的下一个峰: image

我了解到size的参数maximum_filter不等于min_distance中的peak_local_max,但是我想知道是否有一种方法可以两者给出相同的结果。有可能吗?

关于stackoverflow的一些相关问题是:

Get coordinates of local maxima in 2D array above certain value

Peak detection in a 2D array

0 个答案:

没有答案