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)
plotter.plot_set(signal, ylabel = "potential [uV]")
print "Mean amplitude: %s [uV]" % emg_analyser.mean_amplitude(options.channel)
# plot full spectrum
freq, spec = emg_analyser.signal_spectrum(options.channel)
plotter.plot_channel(spec,
xaxis = freq,
xlabel = "frequency [Hz]",
ylabel = "magnitude [(uV s) ^ 2]",
title = "EMG spectrum")
if options.mode in ["trigger", "file"]:
trigger_moments = emg_analyser.get_trigger_times(seconds = True)
emg_moments = emg_analyser.get_emg_times(options.channel,
window = 10,
jump = 2.0,
min_sec_diff = 0.3,
seconds = True)
if options.mode == "trigger":
times = emg_analyser.get_time_differences(trigger_moments,
emg_moments,
max_allowed = 0.3,
full_info = True)
einthoven = data_preparator.einthoven_signal_set()
chann_list = [einthoven.get_channel(1) + einthoven.get_channel(3),
einthoven.get_channel(2)]
plotter.plot_channel_list(chann_list, title_list = ["Einthoven: I + III",
"Einthoven: II"])
## Analysing channel
# peaks
peaks = ekg_analyser.get_peaks(options.channel)
##Difference plot - beats
peaks_sec = ekg_analyser.get_peaks(options.channel, seconds = True)
y1 = ekg_analyser.inter_peak_diff(peaks, seconds = True)
plotter.plot_channel(y1,
xaxis = peaks_sec[1:],
ylabel = 'Roznica [s]',
point_plot = 'x')
##Difference plot - beats
beats = range(1, len(y1) + 1)
plotter.plot_channel(y1,
xaxis = beats,
xlabel = 'uderzenia serca [1]',
ylabel = 'Roznica [s]',
point_plot = 'x')
print "Peak number: ", len(peaks)
inter_peak = ekg_analyser.mean_inter_peak_time(peaks)
print "Mean inter-peak time: %s s" % inter_peak
print "Mean heart-rate: %s beat/min" % (60.0 / inter_peak)
print emo_grow
emo_grow.pop(3)
emo_grow.pop(4)
print emo_grow
print "Emo:"
print "Amplitude: %.2f | %.2f" % (np.mean(emo_amps), np.max(emo_amps))
print "Latency: %.2f | %.2f" % (np.mean(emo_lats), np.max(emo_lats))
print "Growth time: %.2f | %.2f" % (np.mean(emo_grow), np.max(emo_grow))
print "Half life: %.2f | %.2f" % (np.mean(emo_half), np.max(emo_half))
print "Neutral:"
print "Amplitude: %.2f | %.2f" % (np.mean(neu_amps), np.max(neu_amps))
print "Latency: %.2f | %.2f" % (np.mean(neu_lats), np.max(neu_lats))
print "Growth time: %.2f | %.2f" % (np.mean(neu_grow), np.max(neu_grow))
print "Half life: %.2f | %.2f" % (np.mean(neu_half), np.max(neu_half))
print "\n\n"
print emo_half, np.mean(emo_half)
print neu_half, np.mean(neu_half)
plotter.plot_set(signal, ylabel = "potential [uV]")
plotter.plot_channel(signal.get_channel(1),
triggers = triggers,
slopes = increase,
peaks = peaks,
ylabel = "potential [uV]",
title = "Hand electrode")
plotter.show()
