def StaLta(self, slta):
#get values for sta and lta
tas = slta.split(' ')
sta = eval(tas[0])
lta = eval(tas[1])
#first copy data stream. copy.data will be cft
self_cft=self.copy()
# for each trace of self_cft stream apply Sta/Lta
# add sta lta values into the header
for i in range(len(self_cft)):
cft = recSTALTA(self_cft[i].data, \
int(sta * self_cft[i].stats.sampling_rate), \
int(lta * self_cft[i].stats.sampling_rate))
# substitution data-->cft
self_cft[i].data=cft
# add sta lta values into the header
self_cft[i].stats['Sta'] = sta
self_cft[i].stats['Lta'] = lta
return self_cft
def test_recSTALTAC(self):
"""
Test case for ctypes version of recSTALTA
"""
nsta, nlta = 5, 10
c1 = recSTALTA(self.data, nsta, nlta)
self.assertAlmostEqual(c1[99], 0.80810165)
self.assertAlmostEqual(c1[100], 0.75939449)
self.assertAlmostEqual(c1[101], 0.91763978)
self.assertAlmostEqual(c1[102], 0.97465004)
def test_recSTALTAC(self):
"""
Test case for ctypes version of recSTALTA
"""
nsta, nlta = 5, 10
c1 = recSTALTA(self.data, nsta, nlta)
self.assertAlmostEqual(c1[99], 0.80810165)
self.assertAlmostEqual(c1[100], 0.75939449)
self.assertAlmostEqual(c1[101], 0.91763978)
self.assertAlmostEqual(c1[102], 0.97465004)