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()