找到离山顶最近的Minima

Question

我使用scipy.singal.find_peaks查找峰值和最小值，使用：

import numpy as np
from scipy.signal import find_peaks

x=np.array([9.8,57.,53.,37.,24.,19.,16.,15.,13.,13.,12.,12.,11.,11.,11.,11.,11.,11.,10.,13.,13.,13.,15.,13.,12.,14.,15.,14.,51.,34.,24.,20.,24.,22.,18.,57.,63.,38.,27.,28.,31.,33.,94.,71.,48.,40.,43.,39.,31.,27.,22.,21.,20.,19.,18.,18.,19.,20.,20.,49.,62.,48.,43.,34.,33.,28.,26.,26.,24.,23.,23.,26.,27.,70.,97.,57.,46.,68.,82.,59.,49.,37.,40.,45.,36.,33.,28.,22.,23.,284.,524.,169.,111.,148.,98.,68.,50.,38.,30.,28.])
peaks, _ = find_peaks(x)
mins, _ =find_peaks(x*-1)

看上去像是：

​

​

现在，我感兴趣的是找到最接近每个峰值的最小值。这样我就能接受他们之间的区别了。

在查看了find_peaks文档之后，peak_prominece的论点似乎就是我要找的东西。

prominences = peak_prominences(x, peaks)[0]
contour_heights = x[peaks] - prominences

然后看起来是：

​

​

经过检查，peak_prominences发现了峰值之前的最小值。对于我的应用程序，我想要最近的峰值，不管它是在前面还是后面。

如何使用mins来定义peak_prominence计算的参数wlen。？

由于mins包含极小值的索引，我如何使用它来定义wlen？基本上，我必须找到约束每个峰值的指数(即peaks[i])。

是否有更好的方法可以通过使用mins和peaks来实现这一点？

Answer 1

现在，我感兴趣的是找到最接近每个峰值的最小值。这样我就能接受他们之间的区别了。

下面是你要找的东西吗？

closest_mins = [mins[np.argmin(np.abs(x-mins))] for x in peaks]
difference = x[peaks]-x[closest_mins]
print(difference)

47.3. 1.37。4. 45。54.3. 44。51.36.8. 413。37.

下面是用虚线表示的peaks、mins和peaks-closest mins对的图。请注意，有一些mins最接近多个peaks。

plt.plot(x)
plt.plot(peaks, x[peaks],'o', label = 'peaks')
plt.plot(mins, x[mins],'s', label = 'mins')
plt.plot(closest_mins, x[closest_mins],'*', label = 'closest mins')

for p, m in zip(peaks, closest_mins):
    plt.plot([p,m], [x[p], x[m]], 'k', dashes = (4,1))
plt.legend();

​

​