关于Python的FFT 滤波器

学生党少华

最近有一个需求，要对全球海温数据进行FFT滤波，尽管我自己写出来也搜到了好几个不同版本的FFT滤波器，但是效果都达不到origin软件中FFT滤波器的效果，哪位大神能写出来和origin效果一样的FFT滤波器，在下不胜感激，求代码！！

岫鹿森涌

同求

小徐同学

请问楼主解决了吗

学生党少华

小徐同学 发表于 2021-10-25 21:10
请问楼主解决了吗

没看悬赏金额还挂在那上面呢么~~

一大碗年糕

补一个可以试试https://geocat-comp.readthedocs. ... tml#fourier-filters

学生党少华

久远的问题 我自问自答吧。。。 butterworth效果也不差。 用它写了文章
##lowpass

from scipy.signal import butter, lfilter, filtfilt
def butter_lowpass(lowcut, fs, order=4):
    nyq = 0.5 * fs
    low = lowcut / nyq
    b, a = butter(order, low, btype='lowpass')
    return b, a

def butter_lowpass_filter(data, lowcut, fs, order=4):
    b, a = butter_lowpass(lowcut, fs, order=order)
    y = filtfilt(b, a, data)
    return y
##bandpass
def butter_bandpass(lowcut, highcut, fs, order=4):
    nyq = 0.5 * fs
    low = lowcut / nyq
    high = highcut / nyq
    b, a = butter(order, [low, high], btype='bandpass')
    return b, a

def butter_bandpass_filter(data, lowcut, highcut, fs, order=4):
    b, a = butter_bandpass(lowcut, highcut, fs, order=order)
    y = filtfilt(b, a, data)
    return y

雪域小丹

学生党少华 发表于 2023-9-19 19:12
久远的问题 我自问自答吧。。。 butterworth效果也不差。 用它写了文章
##lowpass

请问fs怎么选取？

学生党少华

雪域小丹 发表于 2023-10-27 08:14
请问fs怎么选取？

采样频率嘛 根据数据分辨率来看

紫藤萝

请问0.25度数据，fs是等于0.25吗？
lowcut和highcut是怎么定的呢？
滤波阶数为什么选择4呢？

学生党少华

紫藤萝 发表于 2023-10-31 00:23
请问0.25度数据，fs是等于0.25吗？
lowcut和highcut是怎么定的呢？
滤波阶数为什么选择4呢？

采样频率指的是时间分辨率；order具体原理 我也不是很熟悉，我自己看着效果来的。