首页 > 解决方案 > 在 pandas 和 scipy.signal 中查找局部最大值/最小值的索引

问题描述

我的目标是在 pandas 或 matplotlib 中找到函数的局部最大值和最小值的索引。

假设我们有一个噪声信号,其局部最大值和最小值已经绘制在以下链接中:

https://stackoverflow.com/a/50836425/15934571

这是代码(我只是从上面的链接复制并粘贴它):

import numpy as np
import matplotlib.pyplot as plt
import pandas as pd
from scipy.signal import argrelextrema

# Generate a noisy AR(1) sample

np.random.seed(0)
rs = np.random.randn(200)
xs = [0]
for r in rs:
    xs.append(xs[-1] * 0.9 + r)
df = pd.DataFrame(xs, columns=['data'])

n = 5  # number of points to be checked before and after

# Find local peaks

df['min'] = df.iloc[argrelextrema(df.data.values, np.less_equal,
                    order=n)[0]]['data']
df['max'] = df.iloc[argrelextrema(df.data.values, np.greater_equal,
                    order=n)[0]]['data']

# Plot results

plt.scatter(df.index, df['min'], c='r')
plt.scatter(df.index, df['max'], c='g')
plt.plot(df.index, df['data'])
plt.show()

所以,我不知道如何从这一点继续并在图上找到与获得的局部最大值和最小值相对应的索引。将不胜感激任何帮助!

标签: pandasnumpymatplotlibscipy

解决方案


我你使用df['min'].notna()它返回一个数据框,其中的行min不是nan。要查找局部最小值的索引,您可以使用该.loc方法。

df.loc[df['min'].notna()].index

您的示例的输出结果是:

Int64Index([0, 11, 21, 35, 43, 54, 67, 81, 105, 127, 141, 161, 168, 187], dtype='int64')

您可以使用类似的过程来查找局部最大值。


推荐阅读