P13A05 Développement d’un algorithme de traitement du signal pour des signaux otologiques électro-physiologique appelés potentiels évoqués auditifs » History » Version 9

axel BARRIEUX, 04/08/2021 11:15 AM

1 1 axel BARRIEUX
h1=. Projet GE2-GE3 2013 :
2 1 axel BARRIEUX
3 1 axel BARRIEUX
4 1 axel BARRIEUX
h2<. P13A05 Développement d’un algorithme de traitement du signal pour des signaux otologiques électro-physiologique appelés potentiels évoqués auditifs (PEA)
5 1 axel BARRIEUX
6 1 axel BARRIEUX
p((((. Recherche automatisée du seuil d'audition par le recueil des PEA
7 1 axel BARRIEUX
8 1 axel BARRIEUX
*Entreprise / Client* : ECHODIA / Thierry HASSOUN
9 1 axel BARRIEUX
*Auteurs* : Ibtissame BOUSSIF / Charles ONDET
10 1 axel BARRIEUX
*Responsable Projet* : Jacques LAFFONTTuteur industriel : Isabelle GOI
11 1 axel BARRIEUX
12 2 axel BARRIEUX
[[1. Résumé]]
13 2 axel BARRIEUX
[[2. Abstract]]
14 2 axel BARRIEUX
[[3. Introduction]]
15 2 axel BARRIEUX
[[4. Présentation du Sujet]]
16 2 axel BARRIEUX
[[5. Cahier des Charges]]
17 2 axel BARRIEUX
[[6. Developpement]]
18 1 axel BARRIEUX
19 2 axel BARRIEUX
[[7. Problématiques]]
20 1 axel BARRIEUX
21 7 axel BARRIEUX
[[1. Faisabilité]]
22 2 axel BARRIEUX
[[2. Etude Théorique]]
23 2 axel BARRIEUX
[[3. Solutions]]
24 1 axel BARRIEUX
25 2 axel BARRIEUX
[[8. Gestion de Projet]]
26 1 axel BARRIEUX
27 2 axel BARRIEUX
[[1. W.B.S.]]
28 2 axel BARRIEUX
[[2.Gantt]]
29 1 axel BARRIEUX
30 2 axel BARRIEUX
[[9. Notes d'application]]
31 2 axel BARRIEUX
32 2 axel BARRIEUX
[[1. sujet 1]]
33 2 axel BARRIEUX
[[2. sujet 2]]
34 2 axel BARRIEUX
35 2 axel BARRIEUX
[[10. Bilan]]
36 2 axel BARRIEUX
37 2 axel BARRIEUX
[[1. Etat d'avancement]
38 2 axel BARRIEUX
[[2. Analyse Critique]]
39 2 axel BARRIEUX
[[3. Perspectives]]
40 2 axel BARRIEUX
41 2 axel BARRIEUX
[[11. Bibliographie]]
42 2 axel BARRIEUX
43 2 axel BARRIEUX
p((. *Résumé*
44 2 axel BARRIEUX
45 2 axel BARRIEUX
ECHODIA est une société située à Clermont-Ferrand à l'école de médecine. Cette société fabrique des appareils médicaux ORL. Ces équipements destinés aux médecins ORL et du travail sont conçus et assemblés par ECHODIA. Cette société a créé son propre logiciel d’acquisition et traitement des signaux appelé ECHOSOFT.Le but de ce projet est d'automatiser la détection du seuil auditif. Cette méthode est basée sur le recueil des potentiels évoqués auditifs (PEA), réponse électrique du système auditif à une stimulation acoustique. Les PEA témoignent de la santé du nerf auditif. Un sujet normal présente trois ondes caractéristiques notées ondes I, III et V. La présence de l'onde V témoigne de l'audition du patient pour une puissance donnée.
46 2 axel BARRIEUX
47 2 axel BARRIEUX
Notre démarche s'est déroulée suivant deux axes. Le premier est une étude théorique basée sur une méthode itérative de recherche de seuil automatisée. Le deuxième axe consiste à implémenter cette méthode et réaliser des statistiques sur les paramètres des signaux tels que la corrélation des signaux, l’amplitude de l’onde V.
48 2 axel BARRIEUX
49 2 axel BARRIEUX
L’environnement du développement du logiciel ECHOSOFT est Neatbeans. Cet outil utilise le langage JAVA.
50 2 axel BARRIEUX
51 2 axel BARRIEUX
L’étude théorique a mené vers une machine à états basée sur la méthode itérative de la dichotomie. Cette dernière implémentée conduit à la détection du seuil en concordance avec le seuil réel.
52 2 axel BARRIEUX
53 2 axel BARRIEUX
Les études statistiques ont conduit à des résultats en fonction de la corrélation et de l’amplitude de l’onde V.
54 2 axel BARRIEUX
55 2 axel BARRIEUX
Mots Clés : PEA, Onde V, Dichotomie, JAVA, Seuil.
56 2 axel BARRIEUX
57 2 axel BARRIEUX
---
58 2 axel BARRIEUX
59 2 axel BARRIEUX
p((. *Abstract*
60 2 axel BARRIEUX
61 2 axel BARRIEUX
ECHODIA is a company which is located in Clermont-Ferrand at School of Medicine. This company produces ENT medical devices. Equipment for ENT and occupational doctors are designed and assembled by ECHODIA. This company has created its own acquisition and signal processing software called Echosoft.
62 2 axel BARRIEUX
63 2 axel BARRIEUX
The aim of this project is to automate the detection of auditory threshold. This method is based on the collection of auditory brainstem response (ABR) which is an electrical response of the auditory system at acoustical stimulation. The ABR reflect the health of the auditory nerve. A normal subject has three characteristic waves recorded waves I, III and V. The presence of wave V reflects the patient's hearing for a given power.
64 2 axel BARRIEUX
65 2 axel BARRIEUX
Our approach followed two axis. The first is a theoretical study based on an iterative search method automated threshold. The second axis is to implement this method and perform statistics on signal parameters such as the correlation of signals, the amplitude of the wave V.
66 2 axel BARRIEUX
67 2 axel BARRIEUX
The development environment of the Echosoft software is called Neatbeans. This tool uses the Java language.
68 2 axel BARRIEUX
69 2 axel BARRIEUX
The theoretical study has led to a finite state machine based on the iterative method of dichotomy. This method leads to the implemented detection threshold in accordance with the actual threshold.
70 2 axel BARRIEUX
71 2 axel BARRIEUX
Statistical studies have led to results of the correlation function and the wave V amplitude.
72 2 axel BARRIEUX
73 2 axel BARRIEUX
74 2 axel BARRIEUX
KEY WORDS:
75 2 axel BARRIEUX
76 2 axel BARRIEUX
ABR, wave V, Dichotomy, JAVA, Threshold.
77 2 axel BARRIEUX
78 2 axel BARRIEUX
---
79 2 axel BARRIEUX
80 2 axel BARRIEUX
p((. *Introduction*
81 2 axel BARRIEUX
82 2 axel BARRIEUX
ECHODIA (S.A.S ) est un équipementier dans le domaine du diagnostic médical ORL et Neurologie. Société créée en 2009 par Thierry Hassoun, elle compte 6 employés.
83 2 axel BARRIEUX
84 2 axel BARRIEUX
Cette entreprise produit et commercialise des appareils destinés aux médecins ORL. Ces appareils permettent plusieurs types de mesures tels que l'audiométrie (détection de sourdité), la tympanométrie (détection d'otites), les potentiels évoqués auditifs et les otoémissions (détection de défauts sur le nerf auditif).
85 2 axel BARRIEUX
86 2 axel BARRIEUX
La société ECHODIA a conçu un appareil portatif, l'ELIOS pour effectuer ces mesures. Voir image suivante.
87 2 axel BARRIEUX
88 3 axel BARRIEUX
p<. !https://forge.clermont-universite.fr/attachments/download/16322/Materiel_20140113093754_20140114160841.png!
89 3 axel BARRIEUX
90 3 axel BARRIEUX
*+Figure 1 : L'ELIOS et ECHOSOFT : acquisition et traitement des données.+*
91 3 axel BARRIEUX
92 3 axel BARRIEUX
---
93 3 axel BARRIEUX
94 3 axel BARRIEUX
p((. *Présentation du Sujet*
95 3 axel BARRIEUX
96 3 axel BARRIEUX
L'ELIOS permet de réaliser le dignostic suivant cette procédure:
97 3 axel BARRIEUX
98 4 axel BARRIEUX
p<. !https://forge.clermont-universite.fr/attachments/download/16323/Recueil_20140114161307_20140114162204.png!
99 4 axel BARRIEUX
100 4 axel BARRIEUX
p<. +*Figure 2*: Procédure de test.+
101 4 axel BARRIEUX
102 4 axel BARRIEUX
L'appareil émet un signal acoustique (niveau sonore compris entre 0 et 120 dB) à l'oreille du patient via un écouteur. L'activité du nerf auditif se répercute à des niveaux très précis du crâne et derrière l'oreille du patient. Des électrodes placées à ces endroits récupèrent le signal électro physiologique (PEA : Potentiels Evoqués Auditifs). Le boitier ECHO-DIF permet de faire l'amplification puis l'acquisition des signaux.
103 4 axel BARRIEUX
104 4 axel BARRIEUX
Ce type de test donne les résultats suivants :
105 4 axel BARRIEUX
106 5 axel BARRIEUX
p<. !https://forge.clermont-universite.fr/attachments/download/16324/ImagePEA_20140114163159_20140114163409.png!
107 5 axel BARRIEUX
108 5 axel BARRIEUX
+*Figure3* : PEA précoces obtenus après filtrage et moyennage.+
109 5 axel BARRIEUX
110 5 axel BARRIEUX
Le diagnostic actuel consiste à rechercher les trois pics caractéristiques de ce signal, ces trois pics apparaissent aux alentours des 1ms, 3ms et 5ms, et sont notés I, III et V.
111 5 axel BARRIEUX
Contexte :
112 5 axel BARRIEUX
113 5 axel BARRIEUX
p(((. • Mesure des PEA alternée, réponse de l’activité du nerf auditif.
114 5 axel BARRIEUX
• Moyennage et filtrage des données alternées : Réduction du bruit.
115 5 axel BARRIEUX
116 5 axel BARRIEUX
Diagnostic actuel :
117 5 axel BARRIEUX
118 5 axel BARRIEUX
p(((. • Démarrage à 80 dB.
119 5 axel BARRIEUX
• 17 clics/s.
120 5 axel BARRIEUX
• Décrémentation de 10 dB, pour trouver un seuil normal autour de 20 dB.
121 5 axel BARRIEUX
• 8 mesures : 10 minutes par oreille avec 30% de réjection.
122 5 axel BARRIEUX
123 5 axel BARRIEUX
L'objectif général consiste à détecter automatiquement le seuil auditif à partir de la détection de l’onde caractéristique notée V du signal PEA. Les PEA, signal électro-physiologique permettent de déterminer la santé du nerf auditif du patient.
124 5 axel BARRIEUX
125 5 axel BARRIEUX
---
126 5 axel BARRIEUX
127 5 axel BARRIEUX
p((. *Cahier des Charges*
128 5 axel BARRIEUX
129 5 axel BARRIEUX
130 5 axel BARRIEUX
Cahier des Charges
131 5 axel BARRIEUX
132 5 axel BARRIEUX
• Recherche du seuil (onde V)
133 5 axel BARRIEUX
• Machine à états : méthode itérative.
134 5 axel BARRIEUX
• Gamme de puissance : 80 à 0 dB, avec pas réglable
135 5 axel BARRIEUX
136 5 axel BARRIEUX
p((. - Précision du seuil : 10 dB.
137 5 axel BARRIEUX
138 5 axel BARRIEUX
p(((((. - Temps d’exécution total : diviser par 3.
139 5 axel BARRIEUX
• Corrélation à puissance égales sur une oreille.
140 5 axel BARRIEUX
• Stabilité des latences :
141 5 axel BARRIEUX
142 5 axel BARRIEUX
p(((. -retard par rapport à la normale
143 5 axel BARRIEUX
- retard à puissances égales
144 5 axel BARRIEUX
- droite et gauche
145 5 axel BARRIEUX
146 5 axel BARRIEUX
• Implantation de l’algorithme :
147 5 axel BARRIEUX
148 5 axel BARRIEUX
p(((. - En JAVA sur logiciel ECHOSOFT
149 5 axel BARRIEUX
150 5 axel BARRIEUX
---
151 5 axel BARRIEUX
152 5 axel BARRIEUX
p(((. *Developpement*
153 5 axel BARRIEUX
154 5 axel BARRIEUX
---
155 5 axel BARRIEUX
156 5 axel BARRIEUX
p(((. *Problématiques*
157 5 axel BARRIEUX
158 5 axel BARRIEUX
159 5 axel BARRIEUX
Le problème à résoudre consiste à préciser le seuil d’audition c'est-à-dire la puissance à partir de laquelle le patient entend.
160 5 axel BARRIEUX
La recherche du seuil sera traitée de manière itérative. Cette méthode sera implémentée sous forme de machines à états.
161 5 axel BARRIEUX
162 5 axel BARRIEUX
Plusieurs pistes seront analysées telles que la puissance de départ, l’amplitude de l’onde V de départ, la corrélation.
163 5 axel BARRIEUX
164 5 axel BARRIEUX
Un document technique sera rédigé.
165 5 axel BARRIEUX
166 5 axel BARRIEUX
L'interface graphique sera modifiée par l’implantation du nouvel algorithme et pourra être utilisée par les médecins ORL.
167 5 axel BARRIEUX
168 5 axel BARRIEUX
---
169 5 axel BARRIEUX
170 5 axel BARRIEUX
p(((. *Faisabilité*
171 5 axel BARRIEUX
172 5 axel BARRIEUX
---
173 5 axel BARRIEUX
174 5 axel BARRIEUX
p(((. *Etude Théorique*
175 5 axel BARRIEUX
176 5 axel BARRIEUX
h3<. La recherche automatique du seuil auditif : 
177 5 axel BARRIEUX
178 5 axel BARRIEUX
179 5 axel BARRIEUX
Nous avons utilisé lors du projet la méthode de la dichotomie pour détecter le seuil d’audition. En
180 5 axel BARRIEUX
effet, la dichotomie est, en algorithme, un processus itératif ou récursif de recherche. Le principe
181 5 axel BARRIEUX
consiste à couper en deux parties (pas forcement égales) un espace de recherche qui devient
182 5 axel BARRIEUX
restreint à l’une de ces deux parties.
183 5 axel BARRIEUX
184 5 axel BARRIEUX
Un test est dit positif si les ondes V sont placées sur les deux courbes issues de la répétabilité.
185 5 axel BARRIEUX
186 5 axel BARRIEUX
L'algorithme réalise un premier test à une puissance de stimulation de 80 dB. Si le test est positif, l'algorithme recherche la présence d'ondes V pour les tests réalisés avec une puissance de stimulation inférieure. Si un test est négatif, l'algorithme réalise un autre test à une puissance de stimulation supérieure, et ainsi de suite.
187 5 axel BARRIEUX
Le seuil d’audition est la puissance à partir de laquelle le patient entend. L'algorithme retourne le seuil qui correspond à la puissance du dernier test.
188 5 axel BARRIEUX
189 5 axel BARRIEUX
L’algorithme spécifié est représenté sur l'Illustration suivante :
190 5 axel BARRIEUX
191 6 axel BARRIEUX
p<. !https://forge.clermont-universite.fr/attachments/download/16325/FSM1_20140114163159_20140114164554.png!
192 6 axel BARRIEUX
193 6 axel BARRIEUX
+*Figure 4* : Dichotomie basique.+
194 6 axel BARRIEUX
195 6 axel BARRIEUX
Power = 80 dB
196 6 axel BARRIEUX
Si (onde V est détectée) {
197 6 axel BARRIEUX
Power est décrémentée d’un pas réglable ( Power = Power – pas) }
198 6 axel BARRIEUX
Sinon {
199 6 axel BARRIEUX
Power est incrémenté d’un pas réglable. ( Power = Power + pas) }
200 6 axel BARRIEUX
Seuil =Power
201 6 axel BARRIEUX
202 6 axel BARRIEUX
Il y a 8 valeurs possibles : le nombre maximum d’itérations est 4. Le temps maximal théorique, avec 30 %
203 6 axel BARRIEUX
réjection, à 17 c/s, pour 1000 stimulations est de :
204 6 axel BARRIEUX
205 6 axel BARRIEUX
Tmax = 4 * 1.3 * 1000 * (1/17) = 305 s = 5 minutes
206 6 axel BARRIEUX
207 6 axel BARRIEUX
Ensuite, nous avons généralisé notre dichotomie en fixant la puissance de départ. Cette méthode a
208 6 axel BARRIEUX
permis de diminuer le nombre d’itérations à 2, 3 ou 4 en fonction de la proximité entre le choix de la
209 6 axel BARRIEUX
puissance de départ et la valeur du seuil réel.
210 6 axel BARRIEUX
211 1 axel BARRIEUX
Dichotomie en fonction de l'amplitude de l'onde V pour le test de démarrage à 80 dB :
212 7 axel BARRIEUX
213 7 axel BARRIEUX
p<. !https://forge.clermont-universite.fr/attachments/download/16326/FSMondeV_20140106134554_20140106134613.png!
214 7 axel BARRIEUX
215 7 axel BARRIEUX
+*Figure5* : FSM en prenant en compte l'amplitude de l'onde V pour le test de démarrage à 80 dB.+
216 7 axel BARRIEUX
217 7 axel BARRIEUX
Il s'agit d'une autre piste de recherche, une dichotomie avec un étage de moins, afin
218 7 axel BARRIEUX
de diminuer le temps de diagnostic, et ainsi diminuer la fatigue du patient.
219 7 axel BARRIEUX
Une étude a été conduite sur les données des 25 étudiants de GE4 pour définir une relation entre les
220 7 axel BARRIEUX
amplitudes des ondes V de fortes puissances (80, 70 et 60 dB) afin d'estimer la validité de la
221 7 axel BARRIEUX
méthode. Cette étude a mené à des statistiques sur l’amplitude de l’onde V en fonction du seuil.
222 7 axel BARRIEUX
223 7 axel BARRIEUX
1er test : 80 dB
224 7 axel BARRIEUX
2ème test : 60 dB 40 dB 20 dB
225 7 axel BARRIEUX
3ème test : 70 dB 50 dB 30 dB 10 dB
226 7 axel BARRIEUX
227 7 axel BARRIEUX
Cette piste n’est pas probante par le manque de données et la grande dispersion des données. On ne peut déduire ni de normalité ni de relation
228 7 axel BARRIEUX
entre l'amplitude de l'onde V du test à 80 dB et le seuil d'audition du patient.
229 7 axel BARRIEUX
230 7 axel BARRIEUX
Nous avons par la suite généralisé notre dichotomie, en effet, suivant le principe exposé en illustration 14, le test de départ est réalisé pour une puissance de stimulation comprise dans la gamme 80 à 20 dB, et la recherche suivante est réalisée dans la gamme comprise entre la puissance de départ et 80 dB, dans le cas d'un test négatif pour le premier test, ou 10 dB, dans le cas d'un test positif pour le premier test.
231 6 axel BARRIEUX
232 8 axel BARRIEUX
p<. !https://forge.clermont-universite.fr/attachments/download/16327/DichoGenerale_20140114163159_20140114165007.png!
233 8 axel BARRIEUX
234 8 axel BARRIEUX
+*Figure+* : Dichotomie généralisée.+
235 8 axel BARRIEUX
236 8 axel BARRIEUX
Avec cette méthode, en fonction de la proximité entre la valeur de la puissance de départ et le seuil réel, la valeur du seuil d'audition est retournée entre 2 et 4 itérations.
237 8 axel BARRIEUX
238 8 axel BARRIEUX
h3<. La corrélation :
239 8 axel BARRIEUX
240 8 axel BARRIEUX
Pour valider un test pour une puissance donnée et pour une oreille donnée, deux mesures sont réalisées. Afin d’établir un paramètre de normalité, une étude sur la corrélation des courbes issues de la répétabilité a été réalisée. La corrélation est une méthode statistique qui permet de définir la similitude entre deux signaux.
241 8 axel BARRIEUX
La corrélation est calculée dans la zone des PEA précoces, c'est à dire pour les 10 premières millisecondes du signal.
242 8 axel BARRIEUX
243 8 axel BARRIEUX
Formule utilisée pour calculer la corrélation des deux signaux :
244 8 axel BARRIEUX
245 9 axel BARRIEUX
p<. !https://forge.clermont-universite.fr/attachments/download/16328/FormuleCorrelation_20140106142119_20140106142406.png!
246 9 axel BARRIEUX
247 9 axel BARRIEUX
+*Figure 7* : Formule de la corrélation.+
248 9 axel BARRIEUX
249 9 axel BARRIEUX
La figure représente les résultats de corrélation obtenus, pour l’oreille droite et gauche respectivement :
250 9 axel BARRIEUX
251 2 axel BARRIEUX
p<. !!