Facebook
From Harmless Gibbon, 6 Years ago, written in Plain Text.
Embed
Download Paste or View Raw
Hits: 244
  1. import pyedflib
  2. import numpy as np
  3. import matplotlib.pyplot as plt
  4. from scipy import signal
  5. import pylab
  6.  
  7. f = pyedflib.EdfReader("n7.edf")
  8. n = f.signals_in_file
  9. signal_labels = f.getSignalLabels()
  10. sigbufs = np.zeros((n, f.getNSamples()[0]))
  11. for i in np.arange(n):
  12.      sigbufs[i, :] = f.readSignal(i)
  13. print (sigbufs)
  14. print(signal_labels)
  15.  
  16. sygnal = []
  17. sygnal = sigbufs[0]
  18. fs = 10000
  19.  
  20.  
  21. for i in range(0,len(sigbufs)):
  22.         print(sigbufs[i])
  23.         sygnal[0:90000]=sigbufs[i][10000:100000]
  24.         fs = f.getSampleFrequency(i)
  25.         f, t, Sxx = signal.spectrogram(sygnal, fs)
  26.         plt.pcolormesh(t, f, Sxx)
  27.         plt.ylabel('Frequency [Hz]')
  28.         plt.xlabel('Time [sec]')
  29.         plt.show()
  30.