if not options.input_file:
print "Input file required; type -h for help."
exit(-1)
## construct main agents
bandstop_filter = BandstopFilter(options.frequency, [49, 51])
highpass_filter = HighpassFilter(options.frequency, 0.5)
data_preparator = DataPreparator(options, [bandstop_filter,
highpass_filter])
data_preparator.load_file()
ekg_analyser = EkgAnalyser(options.frequency, data_preparator.start_time)
plotter = SignalPlotter(data_preparator.prepare_timeline())
data_preparator.prepare_data()
data_preparator.apply_filters()
print "%s representation" % options.representation.capitalize()
print "Analysing channel %s\n" % options.channel
# plot signal
if options.representation == "raw":
signal = data_preparator.raw_signal_set()
elif options.representation == "einthoven":
signal = data_preparator.einthoven_signal_set()
elif options.representation == "goldberg":
signal = data_preparator.goldberg_signal_set()
ekg_analyser.load_signal_set(signal)
data = {}
if not options.input_file:
print "Input file required; type -h for help."
exit(-1)
## construct main agents
bandpass_filter = BandpassFilter(options.frequency, [20, 180])
data_preparator = DataPreparator(options,
channel_number,
max_multiplexed,
[bandpass_filter])
data_preparator.load_file()
plotter = SignalPlotter(data_preparator.prepare_timeline())
data_preparator.prepare_data()
if options.mode == "trigger":
data_preparator.apply_filters(exclude = [3])
elif options.mode == "simple":
data_preparator.apply_filters()
else:
data_preparator.apply_filters(exclude = [1])
print "Analysing channel %s\n" % options.channel
# plot signal
if options.mode == "trigger":
signal = data_preparator.trigger_signal_set()
elif options.mode == "simple":
signal = data_preparator.difference_signal_set()
else:
signal = data_preparator.trigger_signal_set(just_trigger = True)
emg_analyser.load_signal_set(signal)
