As a pre-processing step, we resample waveforms using the Lanczos method. We do this here:
|
st_i.interpolate(sampling_rate=decimation_rate, method='lanczos', |
And here:
|
st_i.interpolate(sampling_rate=lowest_fs, method='lanczos', a=a_param) |
Currently we are not applying an anti-aliasing filter first — this can lead to bad results! We need to add code to apply an appropriate lowpass filter prior to the interpolate() method being called.