scipy signal find_peaks_cwt not finding the peaks accurately?

后端 未结 3 1907
一整个雨季
一整个雨季 2020-12-14 02:20

I\'ve got a 1-D signal in which I\'m trying to find the peaks. I\'m looking to find them perfectly.

I\'m currently doing:

import scipy.signal as sign         


        
3条回答
  •  慢半拍i
    慢半拍i (楼主)
    2020-12-14 03:06

    Solved, solution:

    Filter data first:

      window = signal.general_gaussian(51, p=0.5, sig=20)
      filtered = signal.fftconvolve(window, data)
      filtered = (np.average(data) / np.average(filtered)) * filtered
      filtered = np.roll(filtered, -25)
    

    Then use angrelextrema as per rapelpy's answer.

    Result:

    enter image description here

提交回复
热议问题