わかさぎのブログ

プログラミング、Atcoderの勉強とか

サビツキゴレイフィルタ(savgol)

def savgol(y,n):
#サビツキ-ゴレイフィルタ
    signed_y = signal.savgol_filter(y, n, 7, mode="mirror")

    return signed_y

結構お世話になっているフィルタ。

窓の長さは奇数にしないといけない。

参考にしたブログではデータ数/4 +1としていた。

n=len(list)

n=n-n%8

とかすればいい(8の倍数は4で割っても偶数で1足すと奇数になる)

あとはploymoninaiの次数とモードがパラメータ

 

 

scipy.signal.savgol_filter — SciPy v1.7.1 Manual

SciPy で Savitzky-Golay フィルタ | org-技術 (org-technology.com)