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.