FR 2615978 A1 19881202 FR 8707599 A 19870529 DISPOSITIF D'ANALYSE DE SIGNAUX POUR LA DETECTION DE POTENTIELS EVOQUES CEREBRAUX La présente invention concerne un dispositif d'analyse de signaux pour la détection de potentiels évoqués cérébraux. Un potentiel évoqué est la réponse électrique transitoire du système nerveux central à une stimulation. Le potentiel évoqué étant le plus souvent de faible amplitude, il est noyé dans le bruit de fond de l'activité électrique cérébrale spontanée, mais se produit normalement à une latence fixe par rapport à la stimulation. L'enregistrement des potentiels évoqués impose de produire un stimulus sensoriel bref, par exemple auditif, d'enregistrer au moyen de capteurs les réponses électriques sur le scalp au niveau de régions du système nerveux central déterminées par l'anatomie des systèmes sensoriels impliqués et d'éliminer les signaux électriques parasites (bruit) se produisant de façon concomitante au potentiel évoqué afin de pouvoir déterminer un certain nombre de caractéristiques du potentiel évoqué, comme, notamment, la latence par rapport à la stimulation, l'amplitude... L'élimination du bruit pose un problème très difficile en raison du très faible rapport signal à bruit. Dans le cas par exemple de potentiels évoqués auditifs du tronc cérébral (PEATC), le rapport signal à bruit est de l'ordre de -15dB. Une technique actuellement mise en oeuvre pour l'élimi- nation du bruit est celle du moyennage qui met à profit la constance de la latence des potentiels évoqués par rapport à la stimulation. Le moyennage consiste à réaliser un certain nombre de stimulations puis à additionner les réponses obtenues afin de faire émerger le potentiel évoqué par rapport au bruit. Mais la faiblesse du rapport signal à bruit impose un nombre de stimulations élevé, couramment 2000 ou plus dans le cas de PEATC, entraînant de sérieux inconvénients. D'abord la durée des stimulations atteint plusieurs minutes, ce qui est pénible pour le patient. De plus, pour la validité de son interprétation, le résultat obtenu par la technique de moyennage impose la linéarité et la stationarité du milieu 2615 9 78 cérébral pendant toute La durée des stimulations, hypothèses qui ne sont en fait vérifiées que sur des durées de quelques secondes. Ainsi, la présente invention a-t-elle pour but de fournir un dispositif d'analyse destiné à la détection des potentiels évoqués et pouvant procurer des résultats fiables sans demander nécessairement la répétition d'un grand nombre de stimulations. Conformément à l'invention, ce but est atteint grâce à un dispositif comprenant des circuits de réception de signaux électriques recueillis au moyen d'un capteur en réponse à une stimulation et un dispositif de moyennage pour élaborer un signal représentatif, dans le domaine temporel, de la moyenne d'un certain nombre de signaux de réponse recueillis, dispositif dans lequel sont prévus en outre, conformément à l'invention: des moyens de calcul pour élaborer la transformée de Fourier du signal moyenné, et des moyens de traitement pour détecter des sinusoïdes exponentiellement amorties contenues dans la transformée de Fourier du signal moyenné et déterminer des paramètres des sinusoïdes détectées, afin de déduire desdits paramètres des informations caractérisant les ondes de potentiel évoqué contenues dans le signal moyenné. L'invention repose sur le constat que différentes ondes de potentiels évoqués cérébraux, notamment celles du potentiel évoqué auditif du tronc cérébral, peuvent être modélisées par une courbe représentative d'une fonction (lorentzienne) dont la transformée de Fourier directe est une sinusoïde exponentiellement amortie. Ainsi la détection et la localisation des potentiels évoqués peuvent alors se ramener à la détermination des paramètres de fréquences pures exponentiellement amorties. De ces paramètres sont ensuite déduites les caractéris- tiques recherchées des différentes ondes du potentiel évoqué. Ainsi, pour chaque sinusoïde exponentiellement amortie détectée, on détermine la fréquence, le maximum d'amplitude et le coefficient d'amortissement d'o l'on déduit respectivement la latence, l'amplitude et la largeur à mi-hauteur de l'onde correspondante du potentiel évoqué. 2615 9 78 Selon une particularité de L'invention, les moyens de traitement sont agencés pour détecter des sinusoîdes exponentielle- ment amorties et en déterminer les paramètres dans un signal comportant une composante importante de bruit. De la sorte, il n'est pas nécessaire, comme dans la technique antérieure, de réaliser un moyennage par accumulation d'un grand nombre de stimulations. En analyse spectrale paramétrique, on connaît des méthodes adaptées à la détermination des paramètres d'une somme de sinusoïdes exponentiellement amorties et de bruit. On pourra en particulier se référer à la méthode décrite par R. Kumaseran et D.W. Tufts dans un article intituLé "Estimating the parameters of exponentially damped sinusoids and poLe zero modeling in noise" I.E.E.E. A.S.S.P., Vol. 30, 1982, p. 833-840. Toutefois, cette méthode connue perd de son efficacité lorsque le rapport signal à bruit est faible, ce qui est généralement le cas dans le cadre de mise en oeuvre de l'invention. Aussi, de préférence, les moyens de traitement du dispositif conforme à l'invention sont agencés pour pouvoir effectuer de façon fiable une détection de sinusoîdes exponentiellement amorties dans un signal fortement bruité. A cet effet, et comme cela sera décrit de façon plus détaillée par la suite, les moyens de traitement comprennent des moyens d'estimation d'une matrice de covariance directe-inverse (progressive-rétrograde) et de séparation efficace des valeurs propres de la matrice liées au bruit. Il est alors possible de travailler à partir d'un signal moyenné en utilisant un nombre de réponses relativement peu élevé, par exemple moins de 200 au lieu de 2000 comme dans la pratique antérieure. L'invention sera mieux comprise à la lecture de la description faite ci-après, à titre indicatif mais non limitatif, en référence aux dessins annexes sur lesquels: - la figure 1 montre schématiquement les différentes ondes du potentiel évoqué auditif du tronc cérébral; 2615 9 78 - la figure 2 iLlustre la variation dans le temps de -La réponse électrique à une stimulation; - la figure 3 illustre la variation dans le temps d'un signal obtenu par moyennage des réponses électriques à 128 stimulations; - la figure 4 est un schéma par blocs d'un mode de réalisation d'un dispositif d'analyse conforme à l'invention; - les figures 5A et 5B sont des organigrammes montant les opérations effectuées par les moyens de traitement du dispositif conforme à l'invention; - la figure 6 illustre la séparation des paramètres liés au bruit dans la transformée de Fourier du signal moyenné, séparation telle que le dispositif conforme à l'invention permet de la réaliser; et - la figure 7 illustre un exemple de potentiels évoqués estimés respectivement par la technique antérieure et par le dispositif conforme à l'invention. La figure 1 représente le potentiel évoqué auditif du tronc cérébral. Il est caractérisé par un certain nombre d'ondes notées I à VII. Le praticien s'intéresse généralement aux ondes I à V en recherchant, pour chacune d'entre elles, sa latence t. (retard entre la stimulation et l'instant o l'on atteint le maximum d'amplitude) et, le cas échéant, son amplitude maximale a. et sa sa largeur à mi-hauteur bi, i étant alors un indice variant de 1 à 5. A titre d'exemple, la latence t3, l'amplitude maximale a3 et La largeur à mi-hauteur b3 de l'onde III sont indiquées sur la figure 1. L'invention est basée sur une modélisation de chaque onde par une courbe représentative d'une fonction particulière du temps, dite lorentzienne, qui s'exprime sous la forme: a. t étant le temps. + (t - t.)2 b. Dans sa partie utile, le potentiel évoqué peut donc être mis sous la forme d'une fonction du temps y (t) telle que: a. y(t) =. i=1 + (t-t)2 b. Comme déjà indiqué, la détection du potentiel évoqué est rendue particulièrement difficile en raison du bruit dans lequel il est noyé. La figure 2 représente un exemple de réponse éLectrique à une stimulation; il est impossible d'y reconnaître les ondes de la figure 1. La technique de moyennage permet une amélioration du rapport signal à bruit, mais son emploi n'est légitime que si les hypothèses de linéarité et de stationnarité du milieu cérébral sont vérifiées pendant toute la durée des stimulations. Cette durée doit alors être limitée à quelques secondes, ce qui limite en pratique le nombre de stimulations à moins de 200. La figure 3 illustre la variation dans le temps du signal obtenu après moyennage de 128 réponses. Ce signal ne peut être exploité directement pour rechercher les caractéristiques des ondes du potentiel évoqué, ondes que l'on peut simplement deviner. IL est donc nécessaire d'effectuer un traitement spécifique du signal moyenné. A titre indicatif, la courbe A de la figure 6 illustre un potentiel évoqué estimé après un moyennage effectué sur 2048 stimulations conformément à La pratique antérieure. Une exploitation directe de la courbe A est certes possible, mais sans garantie de fiabilité du résultat et au prix d'une épreuve pénible pour le patient. La figure 4 est un schéma par blocs d'un dispositif d'analyse conforme à l'invention. Une électrode de scalp 10 permet de recueillir les réponses électriques à des stimulations auxquelles un patient est soumis, par exemple des stimulations auditives constituées par des sons brefs produits avec une fréquence de répétition de plusieurs hertz au moyen d'un stimulateur 12. 2615 9 78 Un signal recueiLLi par l'électrode 10 (par exemple du type de celui illustré par la figure 2) est amplifié par un amplificateur 14 puis est converti sous forme numérique au moyen d'un convertisseur analogiquenumérique 16. L'information numérique obtenue est reçue par un circuit de traitement 20. De façon classique, le circuit de traitement 20 comprend une unité arithmétique et logique 22, des circuits de mémoire morte (ROM) 24 et vive (RAM) 26 et des circuits d'entrée/sortie 28 interconnectés par bus 30. La commande du stimulateur 12 est assurée par des signaux produits par l'unité de traitement afin de synchroniser les stimulations et les origines des temps des réponses électriques recueillies. Le circuit de traitement 20 fonctionne sous la commande de programmes stockés en mémoire morte. Les opérations réalisées par le circuit de traitement sont indiquées sur les organigrammes des figures 5A à 5B. Les phases successives de traitement sont indiquées de façon générale sur la figure 5A. Elles consistent dans des opérations de: - moyennage des réponses électriques à un nombre N de stimulations (phase 51), transformation de Fourier du signal moyenné (phase 52), détection et détermination des paramètres d'exponentielles complexes dans la transformée de Fourier du signal moyenné (phase 53), et - détermination des caractéristiques des ondes du potentiel évoqué (phase 54). Le moyennage est effectué en cumulant successivement les informations numériques correspondant aux signaux électriques recueillis en réponse aux différentes stimulations. Comme déjà indiqué, l'invention peut être mise en oeuvre en utilisant un nombre de stimulations relativement peu élevé, par exemple moins de 200. La transformation de Fourier est une opération classique. La modélisation des ondes du potentiel évoqué par des courbes 2615 9 78 lorentziennes représentatives de fait que ces ondes se retrouvent sous forme de sinusoïdes exponentiellement amorties dans la transformée de Fourier du signal moyenné. L'opération 53 de détection et de détermination des paramètres de sinusoïdes exponentiellement amorties dans La transformée de Fourier du signal moyenné comprend les phases indiquées dans l'organigramme de la figure 5B. Cet organigramme correspond à la mise en oeuvre d'une méthode originale permettant de détecter de façon fiable des sinusoïdes exponentiellement amorties noyées dans un bruit important. Dans une première phase,la matrice de covariance est estimée en cumulant les covariances élémentaires de P vecteurs progressifs et P vecteurs rétrogrades de dimension M. Chaque vecteur progressif est formé de M valeurs successives de la transformée de Fourier directe du signal moyenné, l'ensemble des P vecteurs progressifs étant obtenu par décalage successif d'un échantillon de la transformée de Fourier directe. Les vecteurs rétrogrades sont obtenus en retournant et en conjugant les - composantes des vecteurs progressifs. Pour s'assurer du fait que la matrice soit de rang M on doit utiliser un nombre total de vecteurs tel que 2 Pff M. Pour identifier R exponentielles complexes on doit choisir une dimension M supérieure à (2 R) + 1, 2 R étant la dimension de l'espace engendré par les R exponentielles complexes et les R exponentielles complexes rétrogrades associées (phase 531). La matrice de covariance est diagonalisée, la diagonalisation est une opération classique. Les 2 R plus grandes valeurs propres sont sélectionnées comme correspondant à l'espace signal. Celui-ci est identifié comme l'espace engendré par les 2 R vecteurs propres qui leur sont associés. Par exemple, dans le cas des PEATC, R est égal à 5 (phase 532). Un vecteur A est choisi orthogonal à l'espace signal de 2 M-i telle sorte que si K, Kzi, Kz2 +..., K est l'un quelconque des -, i estL'nqecnuds signaux à identifier, l'orthogonalité avec A entraîne: a+z +... + aM_1 zM-1 0. 0 M-i i 2615 78 c'est-à-dire que z. est zéro du polynôme A (z) associé à A: A (z) = a0 +a1z = a2z +. + aM_-Z 1. Si M > (2 R) + 1, il existe une infinité de vecteurs orthogonaux à L'espace signal. De façon originale, on introduit un critère de choix (norme minimaLe et contrainte de normalisation) qui permet de minimiser les erreurs dues au bruit résiduel, ce qui conduit à choisir la solution: A A e} +A0 e -i o A = u' U. i = 2 R+1 t = 1/2 argument (A00 A) et o A désigne le vecteur obtenu en inversant l'ordre des _UO ème composantes de A et en les conjugant, U. désigne le ime vecteur propre normé de la matrice de covariance et u li désigne la première composante du vecteur U. . La somme définissant A] porte bien sur les M -2R derniers vecteurs propres (correspondant aux plus petites valeurs propres), tous orthogonaux à l'espace signal (phase 533). Le choix d'un vecteur A orthogonal à l'espace signal assure, comme indiqué ci-dessus, que les paramètres zi caractéristiques des exponentielles figurent parmi les racines du polynômes A (z) associé à A. Grâce au mode original de choix du vecteur A décrit ci-dessus, on assure que l'ensemble des racines du polynôme est réparti de la façon suivante: - les R racines à déterminer zi, de module Izil relation z' = 1/z. i - les M - 2 R - 1 racines parasites supplémentaires de module égal à 1. La répartition des racines ainsi obtenue est indiquée figure 6. Elle permet donc d'identifier sans ambiguité parmi les M - 1 racines du polynôme A (z) les paramètres zi comme étant les R racines de module strictement inférieur à un, ceux liés au bruit étant sur le cercle de rayon unité (phase 534). 26 15 9 78 Les paramètres des exponentielles se déduisent de façon classique des paramètres complexes zi: amortissement Jzi |, fréquence argument (zi); et les amplitudes s'en déduisent par une méthode de moindres carrés (phase 535). REVENDICATIONS 1. Dispositif d'analyse de signaux pour La détection de potentiels évoqués cérébraux, comportant des circuits (10, 14) de réception de signaux électriques recueillis au moyen d'un capteur en réponse à une stimulation et un dispositif de moyennage pour élaborer un signal représentatif, dans le domaine temporel, de la moyenne d'un certain nombre de signaux de réponse recueillis, caractérisé en ce que sont prévus en outre: des moyens de calcul pour élaborer la transformée de Fourier du signal moyenné, et des moyens de traitement pour détecter des sinusoïdes exponentiellement amorties contenues dans la transformée de Fourier du signal moyenné et déterminer des paramètres des sinusoïdes détectées, afin de déduire desdits paramètres des informations caractérisant les ondes (I à V) de potentiel évoqué contenues dans le signal moyenné. 2. Dispositif selon la revendication 1, caractérisé en ce que les moyens de traitement comprennent des moyens pour successivement: estimer une matrice de covariance directe-inverse à partir des valeurs de la transformée de Fourier directe du signal moyenné; diagonaliser la matrice de covariance et sélectionner les plus grandes valeurs propres comme étant celles liées à l'espace signal; sélectionner un vecteur orthogonal à l'espace signal satisfaisant à une condition de norme minimale sous contrainte, de manière à imposer aux racines liées au bruit d'un polynôme associé au vecteur sélectionné de se trouver sur le cercle unité, et rechercher les racines dudit polynôme associé au vecteur sélectionné qui caractérisent des exponentielles amorties. 3. Dispositif selon la revendication 1, caractérisé en ce que les moyens de traitement sont agencés pour calculer l'un au moins des paramètres consistant dans la latence, l'amplitude et la largeur à mi-hauteur des ondes du potentiel évoqué à partir, respectivement, de la fréquence, du maximum d'amplitude et du coefficient d'amortissement des exponentielles amorties détectées. 4. Dispositif selon l'une quelconque des revendications 1 à 3, destiné à la caractérisation des ondes des potentiels évoqués il auditifs du tronc cérébraL.