Calcul d'un intervalle de temps en Python

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:

---

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.

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:

---

1 2 3 4

self.pouls = self.delta/44100 self.pouls = str(self.textBrowser) print self.pouls

---

Autre chose, pour poster du code:
http://www.oqapy.eu/releases/balisescode.png

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 ?

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:

---

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)

---

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 ?

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.

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:

---

1 2	self.pouls = self.delta/44100

---

ça donne quoi comme résultat ?

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:

---

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

Citation:

---

Envoyé par aeiouy49

il ne reconnait pas attrs

---

Moi non plus, d'ailleurs.

et que donne:

Code:

---

1 2	print self.recordedSignal

---

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

Code:

---

1 2	print self.recordedSignal

---

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

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

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