Advertisement
B0b0A

Untitled

Jan 28th, 2020
94
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.04 KB | None | 0 0
  1. #%% Opfati QRS
  2. kolku = int(fs/14)
  3. for k in range(len(opsezi)):
  4. j=0;
  5. for i in range(opsezi[k][0],opsezi[k][0]-kolku,-1):
  6. vrednost = np.round(ekg_filtered2[i],decimals=2)
  7. j = j + 1
  8. pom = np.round(ekg_filtered2[(opsezi[k][0]-kolku):(opsezi[k][0]-j)],decimals=2)
  9. najdeno=np.where(pom==vrednost)
  10. if (len(najdeno[0]) != 0):
  11. break;
  12. opsezi[k][0] = opsezi[k][0]-(pom.shape[0]-najdeno[0][-1])-j
  13. j=0;
  14. for i in range(opsezi[k][1],opsezi[k][1]+kolku,1):
  15. vrednost = np.round(ekg_filtered2[i],decimals=3)
  16. j = j + 1
  17. pom = np.round(ekg_filtered2[(opsezi[k][1]+j):(opsezi[k][1]+kolku)],decimals=3)
  18. najdeno=np.where(pom==vrednost)
  19. if (len(najdeno[0]) != 0):
  20. break;
  21. opsezi[k][1] = opsezi[k][1]+j+najdeno[0][0]
  22.  
  23. #%% Odvoj i nacrtaj QRS
  24. qrs = np.copy(ekg_filtered2)
  25. qrs[0:opsezi[0][0]] = None
  26. for i in range(len(opsezi)-1):
  27. qrs[opsezi[i][1]:opsezi[i+1][0]] = None
  28. qrs[opsezi[-1][1]:] = None
  29. plt.figure()
  30. plt.plot(qrs)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement