Skip to content

Commit

Permalink
Update the call to scipy.signal.butter (#299)
Browse files Browse the repository at this point in the history
  • Loading branch information
Alex L. Urban committed Apr 9, 2019
1 parent a801697 commit 4096a70
Showing 1 changed file with 3 additions and 2 deletions.
5 changes: 3 additions & 2 deletions gwdetchar/omega/core.py
Expand Up @@ -70,8 +70,9 @@ def highpass(series, f_low, order=12, analog=False, ftype='sos'):
gwpy.timeseries.TimeSeries.filter
"""
corner = f_low / 1.5
wn = 2 * corner * series.dt.to('s').value
hpfilt = butter(order, wn, btype='highpass', analog=analog, output=ftype)
fs = series.sample_rate.to('Hz').value
hpfilt = butter(order, corner, btype='highpass', analog=analog,
output=ftype, fs=fs)
hpseries = series.filter(hpfilt, filtfilt=True)
return hpseries

Expand Down

0 comments on commit 4096a70

Please sign in to comment.