def process_data(self):
run_param = self.run_param
data = self.data
process_param = self.process_param
# Applying band pass filter.
data['freq_sample'] = 1.0/data['dt']*1000 # kHz
nyq = 0.5 * data['freq_sample']
low = process_param['freq_low'] / nyq * 1000
high = process_param['freq_high'] / nyq * 1000
b, a = butter(process_param['order'], [low, high], btype='band')
if self.process_param['filter'] == 'filtfilt':
data['LFP'] = filtfilt(b, a, data['LFP'], axis=1)
elif self.process_param['filter'] == 'lfilter':
data['LFP'] = lfilter(b, a, data['LFP'], axis=1)
else:
raise ValueError("process_param['filter'] is not a valid string.")
# Run the rest of the processing.
SphereRand.process_data(self)
