Discussion :

[aeiouy49]

Bonjour,

Je dois faire un projet (à rendre le 22 mai) dans lequel on me demande de calculer le pouls d'un patient à partir d'un signal d'électrocardiogramme.
Pour cela je me suis dit qu'il fallait calculer un deltaT qui correspondrait à l'intervalle entre deux pics ecg.
Ensuite la valeur du pouls sera la division entre deltaT et ma fréquence (44100).

Mon problème est le suivant: je n'arrive pas à transcrire ça sous langage python. j'ai essayé de multiples manières mais sans résultats.
Voici ce que j'ai fait :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
    self.duree = self.Durees_2.value()
    print self.duree
    self.recordedSignal=wavrecord(44100,self.duree)
    print self.recordedSignal
 
    # Normalisation
    self.recordedSignalMax = max(self.recordedSignal)
    self.recordedSignalMin = min(self.recordedSignal)
    if self.recordedSignalMax>self.recordedSignalMin:
      self.recordedSignal = self.recordedSignal/self.recordedSignalMax
    else:
      self.recordedSignal = self.recordedSignal/-(self.recordedSignalMin)
 
    plotWaveForm(self.widget,self.recordedSignal,44100,self.duree)                      
 
    #Calcul du pouls (bmp)
    # delta correspond à l'ecart entre deux pics soit t2-t1
    self.delta = self.recordedSignalMax(t2)-self.recordedSignal(t1)
    self.pouls = self.delta/44100
    self.pouls = str(self.textBrowser)
    print self.pouls

Si quelqu'un peut m'aider, ça serait gentil.

[VinsS]

Salut,

Difficile à déterminer comme ça, tu fais appel à des éléments et des widgets dont nous ignorons tout.

Je vois juste à la fin de ton code que tu redéfinis deux fois self.pouls

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
self.pouls = self.delta/44100
self.pouls = str(self.textBrowser)
print self.pouls

Autre chose, pour poster du code:

[aeiouy49]

Merci pour l'astuce pour insérer du code ^^

J'ai créer une interface à partir de QtCreator qui m'a généré un code et c'est à partir de celui que je peux récupérer certaine balise utiles.

Dans le cas d'un histogramme quelconque, comment fait on pour calculer un intervalle de temps entre deux pics les plus élevés ?

[aeiouy49]

J'ai trouvé ça pour calculer les pics max de mon signal ecg mais après je ne vois pas trop comment faire pour calculer l'intervalle de temps entre ces pics

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
 
# parameters
minheight = np.percentile(ecg, 0.9)
minwidth = sampling * (1 / 0.0015) # R-peak is about 15ms
 
# find R-peak using a translation of findpeak.m
findpeaks = vrp.ecg.findpeaks
peaks, locs = findpeaks(ecg, minheight, minwidth)
beat = zip(peaks, locs)

[VinsS]

44100 est la fréquence de l'électrocardiogramme ?

Si oui alors self.delta / 44100 te donne le temps en seconde entre deux battements et 44100 / self.delta le nombre de battements par seconde.

Je ne sais pas te dire comment tu peux trouver les pics, je ne vois pas comment tu reçois le signal de l'ecg. Tu es vraiment connectée à un appareil ou s'agit-il d'une fonction qui le simule ?

[aeiouy49]

oui 44100 c'est la fréquence de mon signal.

on a un système d'acquisition avec une carte électronique que l'on branche sur l'ordinateur. on relie la carte au patient avec des patchs spéciaux et à l'ordinateur avec un câble.
L'ordinateur reçoit les données et on doit afficher le signal ecg ainsi que le pouls sur l'interface.

[VinsS]

Dans ton code tu sembles avoir ce qui fait cela wavrecord(44100, self.duree)

Il n'y a pas de self devant, c'est un module importé ?

Et plus loin:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
self.pouls = self.delta/44100

ça donne quoi comme résultat ?

[aeiouy49]

waverecord est une fonction qui permet d'enregistrer un signal (celui qui vient de la carte ecg).

self.pouls est ma nouvelle variable pour afficher le pouls. Pour le moment elle n'affiche rien car mon code comprend des fautes.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
 
  def LancerAcquisition(self):
    self.duree = self.Durees_2.value()
    print self.duree
    self.recordedSignal=wavrecord(44100,self.duree)
    print self.recordedSignal
 
    # Normalisation
    self.recordedSignalMax = max(self.recordedSignal)
    self.recordedSignalMin = min(self.recordedSignal)
    if self.recordedSignalMax>self.recordedSignalMin:
      self.recordedSignal = self.recordedSignal/self.recordedSignalMax
    else:
      self.recordedSignal = self.recordedSignal/-(self.recordedSignalMin)
 
    # parameters
    sampling = (self.recordedSignal).attrs.get('sampling_frequency')
    minheight = np.percentile(self.recordedSignal, 0.9)
    minwidth = sampling * (1 / 0.0015) # R-peak is about 15ms
 
    # find R-peak using a translation of findpeak.m
    findpeaks = vrp.recordedSignal.findpeaks
    peaks, locs = findpeaks(self.recordedSignal, minheight, minwidth)
    beat = zip(peaks, locs)
    self.delta = minheight-minwidht
    self.pouls = 44100/self.delta
    plotWaveForm(self.widget,self.recordedSignal,44100,self.duree)

il ne reconnait pas attrs

[VinsS]

il ne reconnait pas attrs

Moi non plus, d'ailleurs.

et que donne:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
print self.recordedSignal

Puisque je vois que tu utilises min() et max()

[aeiouy49]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
print self.recordedSignal

Cette fonction permet d'afficher mon signal ecg dans la fenêtre ploteWaveForm

[VinsS]

Heu non, ça print dans le terminal, rien à voir avec un plotter.

[aeiouy49]

Ah, je ne savais pas. Merci.
J'arrive à visualiser mon interface mais le signal ecg ne s'affiche pas à cause de l'erreur "attrs".
Je vais continuer mes recherches pour en savoir plus sur ce type d'erreur et sa signification

Merci pour votre aide en tout cas