Skip to content
Snippets Groups Projects
Commit 5dab9fc1 authored by Cain, Payton David's avatar Cain, Payton David Committed by Jeremy M Fee
Browse files

Account for even taps, update test

parent ac793232
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
etc/filter/coeffs.json
+ 1
− 1
View file @ 5dab9fc1
{"window": [7.2383512291811e-20, -1.2056376791655608e-06, -4.277196717501353e-06, -7.5805261904864375e-06, -8.59990999152702e-06, -4.352200290023471e-06, 8.043963361551641e-06, 3.0701315868917085e-05, 6.427827433229615e-05, 0.0001073700606602676, 0.0001560816599615354, 0.00020388616910009576, 0.00024185700619301848, 0.00025933209524235326, 0.0002450226563246427, 0.00018852107781718654, 8.209652210988845e-05, -7.739918781232675e-05, -0.0002867500325918265, -0.0005349914806186782, -0.0008027060600497537, -0.0010622553823292288, -0.0012791392276245408, -0.0014145673052782018, -0.0014291915650927083, -0.0012877883456455183, -0.0009645169835343987, -0.00044823496583773606, 0.0002527592710099408, 0.0011072361616178833, 0.0020594556416456693, 0.003030052274417012, 0.003919751576353979, 0.004616024837980743, 0.005002478119825618, 0.00497043638509675, 0.004431861432275871, 0.0033324707796725693, 0.0016637408771624192, -0.0005275870204826471, -0.0031338054694460423, -0.005985429793316386, -0.008856003607624372, -0.011472773771948486, -0.013532939805804448, -0.014724611301843292, -0.014751069175579001, -0.01335648103342966, -0.010350920336899746, -0.00563242456755335, 0.0007960763923794444, 0.008816820028521606, 0.01819561067428389, 0.028589246251964705, 0.03956189575747365, 0.050609556738525054, 0.06119097341299261, 0.07076276709804766, 0.07881608572184706, 0.08491186632556437, 0.08871184974287996, 0.09000279117740173, 0.08871184974287996, 0.08491186632556437, 0.07881608572184706, 0.07076276709804766, 0.06119097341299261, 0.05060955673852507, 0.03956189575747365, 0.028589246251964708, 0.01819561067428389, 0.008816820028521606, 0.0007960763923794446, -0.005632424567553352, -0.010350920336899747, -0.013356481033429665, -0.014751069175579008, -0.014724611301843292, -0.013532939805804448, -0.01147277377194849, -0.008856003607624379, -0.005985429793316389, -0.003133805469446045, -0.0005275870204826471, 0.0016637408771624203, 0.0033324707796725693, 0.004431861432275873, 0.004970436385096754, 0.005002478119825619, 0.0046160248379807445, 0.003919751576353979, 0.003030052274417014, 0.0020594556416456715, 0.0011072361616178837, 0.0002527592710099409, -0.00044823496583773606, -0.0009645169835343994, -0.00128778834564552, -0.001429191565092709, -0.0014145673052782036, -0.0012791392276245408, -0.0010622553823292297, -0.0008027060600497556, -0.0005349914806186782, -0.00028675003259182703, -7.739918781232675e-05, 8.209652210988856e-05, 0.00018852107781718695, 0.0002450226563246427, 0.000259332095242354, 0.00024185700619301848, 0.0002038861691000964, 0.0001560816599615354, 0.00010737006066026785, 6.427827433229639e-05, 3.0701315868917085e-05, 8.043963361551653e-06, -4.352200290023471e-06, -8.59990999152696e-06, -7.5805261904864375e-06, -4.277196717501353e-06, -1.2056376791654907e-06, 7.2383512291811e-20]}
\ No newline at end of file
{"window": [7.2383512291811e-20, -1.2056376791655608e-06, -4.277196717501353e-06, -7.5805261904864375e-06, -8.59990999152702e-06, -4.352200290023471e-06, 8.043963361551641e-06, 3.0701315868917085e-05, 6.427827433229615e-05, 0.0001073700606602676, 0.0001560816599615354, 0.00020388616910009576, 0.00024185700619301848, 0.00025933209524235326, 0.0002450226563246427, 0.00018852107781718654, 8.209652210988845e-05, -7.739918781232675e-05, -0.0002867500325918265, -0.0005349914806186782, -0.0008027060600497537, -0.0010622553823292288, -0.0012791392276245408, -0.0014145673052782018, -0.0014291915650927083, -0.0012877883456455183, -0.0009645169835343987, -0.00044823496583773606, 0.0002527592710099408, 0.0011072361616178833, 0.0020594556416456693, 0.003030052274417012, 0.003919751576353979, 0.004616024837980743, 0.005002478119825618, 0.00497043638509675, 0.004431861432275871, 0.0033324707796725693, 0.0016637408771624192, -0.0005275870204826471, -0.0031338054694460423, -0.005985429793316386, -0.008856003607624372, -0.011472773771948486, -0.013532939805804448, -0.014724611301843292, -0.014751069175579001, -0.01335648103342966, -0.010350920336899746, -0.00563242456755335, 0.0007960763923794444, 0.008816820028521606, 0.01819561067428389, 0.028589246251964705, 0.03956189575747365, 0.050609556738525054, 0.06119097341299261, 0.07076276709804766, 0.07881608572184706, 0.08491186632556437, 0.08871184974287996, 0.09000279117740173, 0.08491186632556437, 0.07881608572184706, 0.07076276709804766, 0.06119097341299261, 0.05060955673852507, 0.03956189575747365, 0.028589246251964708, 0.01819561067428389, 0.008816820028521606, 0.0007960763923794446, -0.005632424567553352, -0.010350920336899747, -0.013356481033429665, -0.014751069175579008, -0.014724611301843292, -0.013532939805804448, -0.01147277377194849, -0.008856003607624379, -0.005985429793316389, -0.003133805469446045, -0.0005275870204826471, 0.0016637408771624203, 0.0033324707796725693, 0.004431861432275873, 0.004970436385096754, 0.005002478119825619, 0.0046160248379807445, 0.003919751576353979, 0.003030052274417014, 0.0020594556416456715, 0.0011072361616178837, 0.0002527592710099409, -0.00044823496583773606, -0.0009645169835343994, -0.00128778834564552, -0.001429191565092709, -0.0014145673052782036, -0.0012791392276245408, -0.0010622553823292297, -0.0008027060600497556, -0.0005349914806186782, -0.00028675003259182703, -7.739918781232675e-05, 8.209652210988856e-05, 0.00018852107781718695, 0.0002450226563246427, 0.000259332095242354, 0.00024185700619301848, 0.0002038861691000964, 0.0001560816599615354, 0.00010737006066026785, 6.427827433229639e-05, 3.0701315868917085e-05, 8.043963361551653e-06, -4.352200290023471e-06, -8.59990999152696e-06, -7.5805261904864375e-06, -4.277196717501353e-06, -1.2056376791654907e-06, 7.2383512291811e-20]}
\ No newline at end of file
geomagio/algorithm/FilterAlgorithm.py
+ 12
− 1
View file @ 5dab9fc1
......@@ -2,6 +2,7 @@ from __future__ import absolute_import
from .Algorithm import Algorithm
import numpy as np
import scipy.signal as sps
from scipy.interpolate import interp1d
import sys
from numpy.lib import stride_tricks as npls
from obspy.core import Stream, Stats
......@@ -68,12 +69,22 @@ class FilterAlgorithm(Algorithm):
if data is None or data == "":
return
window = data["window"]
if len(window) % 2 == 0:
half = len(window) // 2
fill_val = np.average([window[half], window[half + 1]])
window1 = window[:half]
window1.append(fill_val)
window2 = window[half:]
window = window1 + window2
self.steps = [
{
"name": "name" in data and data["name"] or "custom",
"input_sample_period": self.input_sample_period,
"output_sample_period": self.output_sample_period,
"window": data["window"],
"window": window,
}
]
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment