Brevet FR2622990A1 — publié le 19890512

Dispositif de stéréolocalisation spatiale de sources d'activités cérébrales La présente invention a pour objet un dispositif de stéréolocalisation spatiale de sources d'activités cérébrales, c'est-àdire un dispositif permettant de localiser dans le volume cérébral, par rapport à un système de coordonnées spatiales, une ou plusieurs sources à l'origine d'activités électriques cérébrales spécifiques. Un domained'application particulier de l'invention est la localisation spatiale d'une activité cérébrale pathologique, par exemple un foyer épileptique. Il est connu de réaliser une cartographie cérébrale à partir de signaux électroencéphalographiques (signaux EEG) recueillis au moyen d'électrodes implantées en différents endroits du scalp d'un patient. La cartographie cérébrale consiste à établir des cartes illustrant la répartition spatiale de la puissance d'émission à différentes fréquences par des nuances de gris ou des couleurs différentes. Un dispositif d'analyse spectrale permettant d'obtenir de bons résultats en cartographie cérébrale a déjà été décrit dans la demande de brevet français No 87 07 600. Dans cette demande, il est notamment décrit un dispositif d'estimation d'une matrice interspectrale de signaux EEG recueillis par des électrodes, par mise en oeuvre, sur les signaux EEG, d'un algorithme récursif de blanchiment et décorrélation spatiale avec optimisation du compromis biais-variance pour chacun des éléments de la matrice. La cartographie cérébrale permet de fournir des informations précieuses quant à l'existence de sources d'activités cérébrales et à leur emplacement approximatif. Toutefois, en particulier lorsqu'une exérèse doit être pratiquée, il est nécessaire de localiser avec précision la source d'activité cérébrale pathologique. A l'heure actuelle, cette localisation est réalisée dans la majorité des cas par implantation d'électrodes au sein de la masse cérébrale, à travers la boîte cranienne, ce qui est particulièrement traumatisant. Aussi, la présente invention e-t-elle pour but de fournir un dispositif permettant, à partir d'une matrice interspectrale de signaux EEG recueillis, non seulement l'identification de sources d'activités cérébrales, mais aussi la localisation spatiale précise de ces sources au sein du volume cérébral. Ainsi, l'invention vise à fournir un dispositif grâce auquel, en particulier, la localisation de foyers épileptiques peut être réalisée de façon non traumatisante, ce dispositif n'étant pas limité à cette application mais devant permettre, plus généralement, la localisation de toute source d'activité cérébrale pathologique ou non. Conformément à l'invention, le dispositif de stéréolocalisation comporte: (a) un ensemble d'électrodes destinées à être disposées en des emplacements prédéterminés du scalp pour recueillir des signaux EEG, (b) des moyens convertisseurs pour convertir les signaux EEG recueillis sous forme numérique, et (c) une unité de traitement comprenant: - des moyens d'estimation d'une matrice interspectrale des signaux EEG recueillis, des moyens de mémorisation dans lesquels sont enregistrées des informations relatives à un modèle de propagation et permettant de caractériser au moins un vecteur source définissant dans le domaine fréquentiel des relations fonctions de paramètres de position de la ou chaque source, entre un signal- émis par une source dipolaire et les signaux disponibles sur le scalp, le modèle de propagation étant un modèle à trois discontinuités comprenant le cerveau, le crâne et le scalp, et - des moyens de traitement de la matrice interspectrale estimée pour, d'une part, identifier les sources d'activités cérébrales contribuant aux signaux EEG recueillis et déterminer leur nombre et, d'autre part, estimer les coordonnées de chaque source identifiée. De préférence, les moyens de traitement de la matrice interspectrale comprennent des moyens d'estimation des vecteurs propres correspondant aux valeurs propres de la matrice non liées aux sources identifiées, et des moyens de détermination des coordonnées de position de sources pour lesquelles le vecteur source correspondant est orthogonal auxdits vecteurs propres. Le traitement de la matrice interspectrale pour l'identification et la localisation des sources est effectué sur une matrice débarrassée du bruit. Une matrice interspectrale estimée non bruitée peut être obtenue par la technique du moyennage, dans les cas o l'activité cérébrale est provoquée par des stimulations, celles-ci étant répétées pour faire émerger le signal par rapport au bruit. Dans d'autres cas (notamment pour des sources d'activités cérébrales pathologiques), la technique du moyennage n'est pas toujours utilisable. Comme indiqué plus loin de-façon détaillée, on peut parvenir à estimer une matrice interspectrale débarrassée de son bruit à partir de la matrice de cohérence associée à la matrice interspectrale estimée. 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 annexés sur lesquels: - la figure 1 illustre très schématiquement un mode de réalisation d'un dispositif de stéréolocalisation conforme à l'invention, - la figure 2 est un organigramme montrant de façon générale les phases successives mises en oeuvre par le dispositif conforme à l'invention, - la figure 3 illustre un modèle de propagation d'ondes dans un milieu composite avec trois discontinuités de conductivité, et - les figures 4A à 4E illustrent la localisation d'une source par mise en oeuvre de la méthode du projecteur orthogonal sur les données fournies par la matrice interspectrale débarrassée du bruit, au moyen du dispositif conforme à. l'invention. Dans le dispositif illustré par la figure 1, les signaux EEG sont recueillis au moyen d'électrodes de scalp 10, puis amplifiés au moyen d'amplificateurs 12, multiplexés par un multiplexeur 13 et convertis sous forme numérique au moyen d'un os05 convertisseur analogique-numérique 14. Avantageusement, les amplificateurs 12 sont du type de celui décrit dans la demande de brevet français No 2 564 657. Les signaux numérisés sont traités au moyen d'une unité de traitement 20 qui, compte-tenu du volume de traitements à effectuer, est de préférence du type à architecture multiprocesseur organisée autour d'un bus commun 30. L'unité de traitement 20 comprend: - un processeur maitre 21 assurant le séquencement des opérations et le calcul final de la matrice interspectrale des signaux EEG recueillis, - un transformateur de Fourier rapide 22, - un processeur 23 effectuant des calculs de périodogrammes et interpériodogrammes, des prémoyennages et des transformations logarithmiques directes et inverses, intervenant dans le calcul de la matrice interspectrale, - un processeur 24 effectuant des opérations de lissage optimal, - un processeur 25 de calcul de la matrice interspectrale débarrassée du bruit, - un processeur 26 de calcul de projecteur orthogonal associé à des vecteurs sources, - un processeur de localisation 27 effectuant le calcul de coordonnées de sources d'activités cérébrales identifiées, - des mémoires mortes (ROM) 28a destinées à la mémorisation des programmes et d'informations représentatives de vecteurs sources pour un modèle de propagation prédéterminé entre source et électrodes, - des mémoires vives (RAM) 28b pour la mémorisation temporaire de données, et - un circuit d'entrées/sorties 29. Différents périphériques sont connectés à l'unité de traitement 20 dont, notamment, une imprimante 16 et une console de visualisation 18. Les phases principales du programme de localisation de sources d'activités cérébrales mis en oeuvre par le dispositif décrit ci-avant sont indiquées sur la figure 2. Estimation de la matrice interspectrale La phase 101 d'estimation de la matrice interspectrale des signaux EEG recueillis est mise en oeuvre de la façon décrite dans la demande de brevet français No 87 07 600 déjà citée. Essentiellement, l'estimation de la matrice interspectrale est réalisée par mise en oeuvre d'un algorithme itératif de blanchiment et décorrélation sur les transformées de Fourier des signaux EEG recueillis. Le blanchiment d'un signal est effectué après estimation de son autospectre par la méthode de lissage du Log-périodogramme après détermination d'une fenêtre de pondération optimale. Le passage au Logpériodogramme est de préférence précédé par un prémoyennage du périodogramme. La décorrélation de deux signaux est effectuée après estimation de leur cohérence par la méthode de lissage de l'interpériodogramme des deux signaux blanchis après détermination d'une fenêtre de pondération optimale. De par ses propriétés, et pour un modèle de propagation donné entre sources d'activités cérébrales et électrodes disposées sur le scalp, la matrice interspectrale estimée et débarrassée de son bruit permet de déterminer le nombre de sources décorrélées et leur localisation dans le volume cérébral par rapport à un système de coordonnées. Les signaux EEG recueillis ont un très faible rapport signal/bruit. Une utilisation de la matrice interspectrale estimée pour identifier et localiser des sources de façon sOre n'est possible que dans la mesure o la matrice estimée est non bruitée. L'élimination du bruit peut résulter de la mise en oeuvre de la technique du moyennage, c'est-à-dire en effectuant la moyenne de résultats obtenus par suite de stimulations répétées des sources. Lorsque la technique du moyennage n'est pas applicable, par exemple lorsque l'activitité des sources n'est pas liée à des stimulations extérieures contrôlées et n'a pas le caractère répétitif nécessaire pour faire émerger le signal par rapport au bruit par moyennage, il est nécessaire de débarrasser la matrice interspectrale estimée du bruit (phase 102). Estimation de la matrice interspectrale débarrassée du bruit Dans l'hypothèse de bruits décorrélés de puissance quelconque sur les électrodes, hypothèse vérifiéedans la majorité des applications neurophysiologiques, la matrice interspectrale Ix peut se mettre sous la forme de la somme d'une matrice signal y de e rang r inférieur à K et d'une matrice de bruit W diagonale mais non scalaire o O O rX = ir + ae x yE r étant la dimension du sous-espace sources (nombres de sources o décorrélées) et K le nombre d'électrodes, tandis que - la notation o signifie que l'on se trouve dans le domaine fréquentiel, - la notation x se rapporte aux signaux EEG recueillis, - la notation y se rapporte à la composante des signaux EEG provenant de la ou des sources à identifier et localiser, et - la notation e se rapporte à la composante de bruit des signaux EEG. L'estimation de la matrice interspectrale débarassée du o O bruit coneiste à rechercher une estimée Iy de ry à partir de y y l'estimée rx obtenue précédemment dans la phase 101 (la notation x signifie qu'il s'agit.d'une estimée). On peut montrer qu'il suffit, pour déterminer une matrice d'ordre K de façon unique à partir de ses éléments non diagonaux, que le rang r0 de cette matrice soit inférieur ou égal à (K/2 - 1) et que tout mineur (ro, ro) soit de rang ro. Mais, en pratique, en raison du bruit et du biais -35 introduit par l'estimation de la matrice interspectrale, on ne peut affirmer qu'il existe une matrice de rang ro dont les éléments non diagonaux sont identiques à ceux de x. De la sorte, et bien qu'en principe seuls les éléments diagonaux de x sont affectés par le X, bruit, on ne peut rechercher une estimée de ir'y directement à partir Y des éléments non diagonaux de Wx. x Une approche par minimisation d'un critère quadratique portant sur la différence entre les éléments non diagonaux de Wx et ZO x Wy montre que l'on peut se ramener à un traitement de la matrice de cohérence C associée à t, c'est-à-dire à la matrice dont les éléments C.z. sont donnés par: cD = j()12 J- s..,l= J & les A. étant les éléments de rx i étant l'indice de ligne et j l'indice de colonne. En effet, un critère quadratique de la forme Q = '"(_..:],,,C,j), dans leqbel w(i, j) est la pondération, c'est-à-dire: O i;) =z Yg O) peut se ramener, en remplaçant la pondération par la valeur approchée (.),: 2 2 rt-jô -ÈO Q=z 2 i/ |i Q/ i i i;,j r7 --r o7 L _, - V o QOP, QO à Q représentant la grandeur Q normalisée de manière à faire apparaître une erreur relative. Le traitement mis en oeuvre sur la matrice de cohérence estimée x revient donc à rechercher, dans un ensemble lC (r)t de matrices définies non négatives de rang r inférieur ou égal à ro, o la matrice Cm (r) qui minimise le critère. Elcij(r)- cx, QXP(r)= i !IccxI- Cette recherche est effectuée à l'aide d'un algorithme de type gradient matriciel à pas de convergence optimal. D'une façon générale, le processus consiste à rechercher, par valeur croissante de r, dans la suite des IQ (r)t, c'est- à-dire les minima de Qop (r), le rang ro correspondant à la plus faible valeur de r pour laquelle QoP (r) est voisin d'une borne QOF (r) définie par i. QOg(r) *L X i;j U* Le étant la largeur équivalente de la fenêtre de lissage fournie par la méthode optimale d'estimation spectrale des éléments C,, de o ' la matrice Cx, comme décrit dans la demande de brevet français No 87 07 600 déjà citée. Pour un faible rapport signal sur bruit, la borne Qop (r) tend vers une valeur Q P (r) définie par: ?O-(r) = [(l- r)/(K-, > (]). Si C(ro) est la matrice qui correspond au minimum Qop (ro), la matrice i' débarrassée du bruit sera estimée par la 0 y matrice r d'éléments: O Y Yxij =7xii YYjj óij(tO) L'algorithme de type gradient matriciel pour la recherche du minimum sera maintenant décrit. La différentielle dQ P (r) du critère QOP (r), qui est nulle au minimum de ce critère, s'exprime par dQ( c [(r) -.xj) dij(r)+(ci(r) Xi dcj(r)] [ i lxijl o ce qui, compte-tenu du caractère hermitien des matrices C (r) et X donne: óx'donne 2. 2(cji(r)- c xij)dcji(r) dQ P(r)=, A Wj Xii Si l'on désigne par 1(r) la matrice définie par les éléments hii(r)=O et (r) rl-c pour i différent de j, et en raisonnant sur le numérateur Q (r) de Qop (r) pour alléger les notations, on obtient o dQ(r) = 2 t ( H(r) dC(r)}. ce qui, la matrice L(r) étant factorisée sous la forme (r) =(r) (r), conduit à la relation cf(r) =11(r) + (r), conduit à le relation: dQ(r) = 2tr { H(r) (r) d(r)+} + 2tr H(r) dM(r) M la notation "tr" signifiant "trace de". La différentielle dQ(r), et, par suite la différentielle DQ P (r) sont nulles siHo(r)M(r) = 0. Pour trouver le minimum ainsi défini, on forme un algorithme de type gradient tel que o o o EN(r) = - H(r) M(r) avec 0 Pour j suffisamment petit, l'accroissement SQ(r) de Q(i peut être assimilé à sa différentielle, ce qui donne: o o +, 3Q(r) =-4 gtr H M(r) M(r) '} . La trace defl k(r) (r)*+Dl représente la somme des carrés des modules des éléments de la matrice9l H dr). Elle est donc positive ou nulle et ne peut être nulle que si tous les éléments de la matrice IHM(r)sont nuls, c'est-à-dire pour un extrémum de Q(r) ou, de façon équivalente, de Q0 P(r). Pour une suite de matrices tMk(r)3, l'algorithme peut donc être défini par: a O O O Mk+l(r) = Mk(r) - Hk kMk(r) Si l'on considère la suite correspondante des valeurs Qk(r) du critère, pour Mk suffisamment petit, l'on obtient: _ o o 1 o Qk+l(r) - Qk(r) = -4 k.tr(I Ekk(r) Mk(r) k}'k Cette variation est négative tant que la matrice nM(r) ne se réduit pas à la matrice nulle. -Les valeurs [Qk (r)} du critère forment donc une suite décroissante-qui converge vers un minimum local de Q(r)ou Q0 P(r). La vitesse de convergence vers le minimum peut être accélérée en recherchant le pas de convergence optimal, c'est-à-dire, dans la direction du gradient, le Mk qui donne le minimum du critère. Le Mk est déterminé de façon approchée en faisant l'hypothèse d'une variation localement quadratique. A partir de la relation ci-dessus donnant M!R+(r), on peut écrire: 0 0 + 0O+ 0 0 ', 0 0 O+ O+ Mk+l =.- O++ kM HM)M+ f. Mk 1M+lM is k k 9 k k + k k Ek) +, kMkMk k- En négligeant le terme enlk relatif au carré de A A l'erreur, on a: k +l k+l k] Oa O Mk+IlMk: - c =Mk Mk - C-I Dk O O O O+ O ' O+ O avec Dk= Hk Mk Mk + Mk Mk Hk'. Par IUie on o l o 0 0 0 Hk+l(gk)=]k ' k Dk pouri j; o 2 02 o o o o 02 Hk+l(90k)= Hk gk k k - gk Dk Hk+gk2 Dk; H 2 I00' tr{[ k+l(gk)]) 2tr {[ Hk- -k Dk] Dk) De la dernière relation, en annulant la dérivée, il vient:o o tr(Hk Dk) ilk 2 'tr( k) u'{Dk} ce qui permet de calculer leó&k optimal à chaque itération de k. On décrira maintenant un processus permettant d'initialiser l'algorithme afin d'éviter une éventuelle convergence vers une valeur erronnée par suite d'un mauvais choix de conditions o initiales (c'est-à-dire les cij(r). Ce processus consiste en une méthode simplifiée servant à estimer, par valeur inférieure, l'ordre de grandeur du rang recherché et à fournir une estimée de R4(r) pour initialiser l'algorithme précédent. La matrice C x étant définie non négative peut se décomposer de façon unique sous le forme d'un produit de matrices triangulaires avec éléments positifs sur la diagonale. Dans le cas d'une matrice de rang r inférieur à K, cette décomposition prend la forme suivante: =B B P2= BIB2 BF = =, Pu-- B2 B! È z2 Les éléments non diagonaux de Cx étant connus, il suffit de déterminer les r premiers éléments diagonaux pour que la matrice soit complètement définie. En effet, le bloc diagonal é1 étant alors supposé connu, on peut, par un algorithme classique, le décomposer en un produit tel que IP = a l. Ainsi, le blanc B2est déterminé par la ni relation B2 P [BD A ce qui permet de calculer les éléments diagonaux manquants de x à partir du produit matriciel: O O o0+ 74 = B2 B2. La détermination du bloc diagonal et du rang deC X débarrassée du bruit, donc la détermination des r premiers éléments diagonaux est réalisée comme suit. On désigne par w(j)l les vecteurs colonnes de C et par j10tw(j) les vecteurs formés par les (K -.r - 1) dernières composantes des [ (j)}, et par 3 la matrice formée par r vecteurs: ^ A {w(1), w(2)..., w(0) En général, le vecteur W(r + 1) ne s'exprimera pas comme combinaison linéaire de ces r vecteurs du fait que la matriceCest biaisée et bruitée. On définit alors r A w (r+1) = 50 G)@, j=I o les (j) sont déterminés en minimisant le critère A 2 sp w(r+l)- w(r+1)l QS(r)=^ I w (r+l)lI Si la matrice à est de rang r, le vecteur paramètre correspondant au minimum est donné par A 0+ o1 0+ È (r) :[B+ B] B w(r+1) En considérant la matrice de rang r, les paramètres { j sont les estimées des coefficients du (r + 1)ème vecteur colonne de la matrice de cohérence "débarrassée" du bruit en fonction des r premiers. Pour les indices de ligne {h} o h est inférieur ou égal à r, les éléments diagonaux a sont donc x[h estimés à partir de la relation: t r - Cxh(r+l) = ç Cxhj + (h) C Xhh o l'on supposera que Z(h) est différent de zéro. Pour une première estimation de l'ordre, on peut raisonner sur la valeur du critère QSP(r). Pour une matrice non perturbée QSP(r) s'annulerait pour r supérieur ou égal à ro. En présence d'erreur statistique, et comme, plus haut, pour le critère Q P(r), on peut déterminer, en l'absence de biais, un seuil supérieur QsP(r) pour l'espérance du majorant du minimum du critère QSP(r) Ce seuil supérieur est donné par: O (r) - (K-r-l) (l+l1 (r)l2) /Le > Q(r Iw(r+1)I A faible rapport signal sur bruit, cette borne tend vers K-2r-1 r QO(r) K-r-l Q(r) Dans la mesure o Q5P(r) est toujours un majorant du- seuil d'arrêt, le rang estimé à partir de ce critère (qui ne fait intervenir qu'une sous-matrice de c x) correspond bien à une estimation par défaut. Le bloc ' de la matrice Cx étant déterminé comme 1. o indiqué au sousparagraphe précédent et les matrices î3 et 1B2, dont l'ensemble définit la matrice M(r), étant calculées, on oz 0 0 o remarque que les matrices Cx et (r) ( C(r) = (r)M+(r)) ne diffèrent que par le bloc (P B) P4 (P4 =]B2]B2). La matrice C(r) étant déterminée entièrement à partir des (r + 1) premières colonnesdecles lignes et les colonnes sont ordonnées de façon à grouper dans ces colonnes les éléments pour lesquels le coefficient de cohérence est le plus-grand. Pour ce classement il suffit, par exemple, de s'intéresser à la sousmatrice inférieure et de: - rechercher dans cette sous-matrice le coefficient de cohérence maximum cx. et, par une permutation de l'ordre des indices de la matrice'de le positionner sur la dernière ligne de la première colonne, - combiner la première colonne aux suivantes de façon à faire apparaître des zéros sur la dernière ligne (à l'exception de l'élément de la première colonne), - rechercher dans la matrice ainsi obtenue, privée de sa dernière ligne et de sa première colonne, l'élément maximum et, par une nouvelle permutation de l'ordre des indices de la matrice, positionner cet élément sur l'avant dernière ligne de la deuxième colonne, et ainsi de suite. Ce classement conduit à une matrice triangulaire. Les éléments diagonaux de cette matrice sont rangés, de la gauche vers la droite, par valeur décroissante. Le passage par la recherche de combinaison entre les éléments de la matrice, dans ce classement, se justifie par le fait qu'il est nécessaire de s'assurer que la grandeur des coefficients de cohérence ne résulte pas de la proportionnalité de deux colonnes (ou de deux lignes) de la matrice Cx- Le processus d'estimation de la matrice interspectrale - débarrassée du bruit décrit ci-avant suppose seulement que les bruits sont décorrélés. Si l'on se place de plus dans l'hypothèse bruit blanc, l'estimation du sous-espace source, peut encore être améliorée. Dans le cas o tous les éléments de la matrice sont utilisés, le critère à minimiser s'écrit alors: K A 2 _k2(i) Qop() =tr C(r)-Cx== Kr ( tr(c.} È 2(i) i=l o C(r) est une matrice de cohérence définie non négative de rang inférieur ou égal à r et Cx la matrice de cohérence estimée -et x bruitée déduite de la matrice de covariance estimée et bruitée par la relation: Y ^ 7xi t _ _ _ _ _ _ cxiji t'; Sous l'hypothèse bruit blanc, la matrice C(r) peut s'écrire C(r) = M(r) M(r) + 2(r) D- o M(r) M +(r) représente la matrice de cohérence débarrassée du 2 - bruit au rang r, G (r) la puissance du bruit blanc au rang r et ID x la matrice diagonale définie par: I YIll A x = _0 Yxkk Dans ces conditions, il est possible de déterminer le f2(r) qui minimise Q' P(r), trA - (r) (} Dx tr {C.(r) - M (r) MV+() D x ^2 tr Dx) Ainsi, si l'on pose A,, A A1 Cx = &xt (r) Dx le rang r optimal est défini en minimisant le critère: lcij(r) - c^, 12 c xij Q P(r) =12 i;jxij au lieu du critère Q P (r). Le reste du processus décrit précédemment reste inchangé et s'applique directement. "2 Le calcul itératif sur r et sur C(r)permet d'identifier le sous-espace source en déterminant ro, d'estimer la puissance du bruit blanc et, de plus, fournit une estimation de la matrice de covariance débarrassée du bruit. L'estimation de la matrice interspectrale débarrassée du bruit permet d'évaluer le nombre ro de sources décorrélées (cette approche pouvant s'étendre au cas de sources corrélées) et de fournir, par le calcul des valeurs propres associées, une information relative à la puissance d'émission de ces sources. Mais l'estimation de la matrice interspectrale débarrassée du bruit permet aussi la localisation des sources, dans la mesure o un modèle de propagation entre sources-et électrodes a été défini. Définition d'un modèle de propagation Le modèle utilisé ici est un modèle à trois discontinuités de conductivité électrique schématisé par trois sphères concentriques de rayons respectifs a1, a2, a3 représentant le cerveau avec le, liquide céphalo-rachidien, le crâne et le scalp (figure 3). On se place dans l'hypothèse d'une charge ponctuelle placée en un point M à distance b du centre O des sphères. Les o discontinuités de conductivité sont très importantes, puisque l'on peut atteindre un rapport de 100 entre la conductivité Q2 du crâne et les conductivités Q' et q3 du cerveau et du scalp qui sont voisines l'une de l'autre, et puisque la conductivité r4 à l'extérieur du scalp est très différente de celle du scalp. Du fait de ces importantes discontinuités de conductivité, la contribution en sortie des ondes transmises intermédiaires peut être négligée. On ne considère donc, en sortie, que la contribution des ondes Dnm ayant subi n réflexions entre les deux couches 31, 32 du côté interne et m réflexions entre les deux couches 32, 33 du côté externe. La position Mnm sur le segment OM de l'image (source virtuelle) dont est issue l'onde Dnm est donnée par 1 1. ( a +m - t Quant à l'amplitude C nm de l'onde Dnm, elle est donnée par C: ( rn)ZC} i (C) (CO n avec -., _). o(0), ,q _ __ e cri__ __r_ _ c o -- (a..+ 5)Cî2 a.)Co la sphère pour un milieu homogène. '(S, OMnm) s'exprime par: O4=) l/ Ir2 avec r2=1+ 2q cosA La double somme infinie donnant le potentiel VR (d' b) peut être simplifiée par approximation par une intégrale continue de la forme tf R C cc' ex? oc o '- avec o i C jI - 1i'=a a;O (. Loc CO+ Log C) c. ai â'a?;e a -30 L'intégrale ci-dessus peut être calculée par la méthode des trapèzes en discrétisant la variable t, ce qui revient à placer des charges régulièrement espacées entre le centre 0 de la sphère et la position M de la charge réelle avec des poids donnés par la fonction de pondération: Le modèle décrit ci-avant pour une charge ponctuelle est facilement étendu à une source dipolaire qui offre une meilleure représentation d'une source d'activités cérébrales. Le passage à la source dipolaire est obtenu numériquement en sommant les potentiels 1 i créés par deux charges l et (o 1 caractérise le moment dipolaire et e leur distance) placées à distance (b - e) et b du centre des sphères avec e très petit devant b, et en remarquant que, dans le cas d'une sphère conductrice le potentiel à la surface de celle-ci peut s'interpréter comme la somme du potentiel de deux dipôles dans l'espace libre de même conductivité soit, finalement, à partir de la charge en Mnm: A titre d'exemple, en calculant par la méthode des trapèzes l'intégrale définie plus haut donnant une approximation du potentiel résultant'R (y, b), en utilisant cinq trapèzes (charges virtuelles placées à 1/8, 1/4, 3/8, 1/2 et 3/4 de la distance OMo à partir de 0, on obtientpour le potentiel créé par un dipôle réel: 35.9 c* xn et udp(.*2o _ 6Avé lû{n [3 p A)xXpex t)] to(,4p+ ep çs3ep 93 ] ) +[txp( X À N (Ai 2)P[) X 9Ai 08 (,_y avec sec= ci Le modèle de propagation qui vient d'être décrit permet donc de définir un modèle de vecteur source g i) dont les composantes représentent les valeurs des potentiels recueillis (à une constante près) sur les K électrodes en fonction des coordonnées de la source, coordonnées exprimées par un vecteur de positiona relatif à la source. En d'autres termes, le vecteur source A (%) caractérise la fonction de transfert dans le domaine fréquentiel entre la position de la source et les signaux EEG recueillis sur les électrodes. Le modèle de vecteur source (e)correspondant au modèle de propagation sélectionné est enregistré en mémoire ROM du dispositif de traitement pour être utilisé aux fins de localisation des sources identifiées dans la matrice interspectrale débarrassée du bruit estimée y. Localisation des sources De préférence, la localisation des sources est effectuée par mise en oeuvre de la méthode du projecteur orthogonal, ou encore la méthode du goniomètre généralisée (phase 103). Ayant constaté que pour le modèle dipolaire, le vecteur source est indépendant de la fréquence, on localisera la source, non pour chaque plan fréquentiel de la matrice 'y, mais pour une matrice I'o obtenue par moyennage sur les différents plans la fréquentiels. La matrice I[ypossède (K - ro) valeurs propres qui sont y voisines de zéro (si l'estimation était parfaite, ces (K - r) valeurs propres seraient nulles). A ces (K ro) valeurs propres correspondent des vecteurs propres{} qui sont orthogonaux aux vecteurs sources {Ci} qui forment le sous-espace source de dimension ro. On désigne pare(p) la fonctionnelle obtenue en projetant le modèle du vecteur source;(ksur les valeurs proprestI, avec normalisation par p(c): lY E (of) I vjt étant le conjugué transposé de v. A chaque fois que le vecteur de position cQ sera égal au vecteur de position s d'une source, la fonctionnelle ?gp) doit etre nulle à cause des fluctuations statistiques. C'est donc la recherche des minima de cette fonctionnelle (ou les maxima de son inverse) qui donnera les vecteurs de position recherchés (i des sources. On notera que, dans l'expression ci-dessus de (), la normalisation par I|(e)I est introduite pour tenir compte du fait que la norme des vecteurs sourcesn'est pas obligatoirement constante. La recherche de minima de e(c) est effectuée en balayant l'ensemble des valeurs du vecteur de position. En variante, la localisation des sources pourra être effectuée en utilisant le sous-espace source. Une méthode peut consister alors à minimiser la distance entre a [ r au,t 5:= Xi et {i étant les valeurs propres et vecteurs propres du sous-espace source, et t est la densité spectrale de puissance de la source i. Toutefois, par rapport à la méthode du projecteur orthogonal, cette méthode présente le désavantage de faire porter la recherche des paramètres de position sur tous les paramètres à la fois. On notera que la méthode de localisation par projecteur orthogonal décrite -ci-avant en liaison avec le modèle de propagation envisagé impose une disposition des électrodes dans un même plan et suppose une orientation radiale du dipôle représentant la source. Le choix du modèle dipolaire radial est ici justifié en raison de l'orientation radiale des cellules pyramidales au niveau du cortex. La disposition des électrodes dans un plan permet une -9 t2622990 localisation uniquement dans ce plan. Une localisation en trois dimensions est possible en effectuant des localisations dans plusieurs plans parallèles en liaison avec une cartographie donnant dans chaque plan les niveaux de puissance d'émission. La localisation en trois dimensions peut être directement obtenue en élaborant un modèle de propagation plus complexe faisant intervenir un paramètre de position supplémentaire et une orientation quelconque du (ou des) dipôles. Exemple Les figures 4A à 4E se rapportent à un exemple de mise en oeuvre du dispositif selon l'invention pour la localisation d'un foyer épileptique. La figure 4A montre les emplacements des 16 électrodes utilisées. La figure 4B montre les signaux analogiquas EEG recueillis sur les 16 électrodes. @ Une matrice interspectrale ix a été estimée à partir des x signaux EEG recueillis par les-électrodes coplanaires No 4, 10, 9, 13, 14 et 8. La diagonalisation de la matrice de cohérence x associée à la matrice interspectrale estimée Ix', dans la bande delta (0,1 à 4 Hz) fait apparaître six valeurs propres non nulles (autant que de capteurs) dont deux paraissent prépondérantes. Une matrice ry débarrassée du bruit a été estimée à partir de Ox comme décrit plus haut. De l'estimation de ry, il x ressort que la matrice a un rang égal à deux. La diagonalisation de la matrice de cohérence Cy associée à Oy fait apparaître deux y -y valeurs propres distinctes auxquelles peuvent être associées deux sources décorrélées. La localisation des deux sources a été effectuée par recherche des maxima de la fonction 1/e(e), ce qui est équivalent à la recherche des minima de e(e), en utilisant la méthode du projecteur orthogonal avec deux vecteurs source élaborés à partir du modèle de propagation à trois discontinuités de conductivité, comme décrit plus haut. La figure 4C montre la cartographie réalisée à partir de la matrice interspectrale estimée Ix'? dans les bandes alpha lent (7,5 à 10 Hz), bêta rapide (18 à 40 Hz), thêta (4 à 7,5 Hz), bêta lent (12,5 à 18 Hz), delta (0,1 à 4 Hz) et alpha rapide (10 à 12,5 Hz). Les cartographies dans la bande thêta et delta indiquent la présence d'un foyer épileptique dans la région temporale droite. La figure 4D montre la variation de la fonction l/e(e) dans le plan défini par les électrodes sélectionnées. Deux maxima de la fonction 1/e(p) apparaissent sur la figure 4D, les positions de ces maxima étant données sur la figure 4E par une croix et un astérisque. Sur la figure 4E, l'axe vertical passant par le centre O correspond à celui des électrodes No 14 et 10. La croix correspond à la localisation du foyer épileptique, tandis que l'astérisque reflète vraisemblablement une "source" relative à des mouvements oculaires. REVENDICATIONS 1. Dispositif de stéréolocalisation spatiale de sources d'activités cérébrales, caractérisé en ce qu'il comporte: (a) un ensemble d'électrodes (10) destinées à être disposées en des emplacements prédéterminés du scalp pour recueillir des signaux EEG, (b) des moyens convertisseurs (14) pour convertir les signaux EEG recueillis sous forme numérique, et (c) une unité de traitement comprenant: - des moyens d'estimation optimale d'une matrice interspectrale des signaux EEG recueillis, - des moyens de mémorisation dans lesquels sont enregistrées des informations relatives à un modèle de propagation et permettant de caractériser au moins un vecteur source définissant des relations, fonctions de paramètres de position de la ou chaque source, dans le domaine fréquentiel, entre un signal émis par une source dipolaire et les signaux disponibles sur le scalp, le modèle de propagation étant un modèle à trois discontinuités comprenant le cerveau, le crâne et le scalp, et - des moyens de traitement de la matrice interspectrale estimée pour, d'une part, identifier les sources d'activités cérébrales contribuant aux signaux.EEG recueillis et déterminer leur nombre et, d'autre part, estimer les coordonnées de chaque source identifiée. 2. Dispositif de stéréolocalisation selon la revendication 1, caractérisé en ce que les moyens de traitement de la matrice interspectrale comprennent des moyens d'estimation, sur la matrice interspectrale débarrassée du bruit, des vecteurs propres correspondant aux veleurs propres de la matrice non liées aux sources identifiées, et des moyens de détermination des coordonnées de position de source pour lesquelles le vecteur source correspondant est orthogonal auxdits vecteurs propres. 3. Dispositif de stéréolocalisation selon l'une quelconque des revendications 1 et 2, caractérisé en ce que les moyens de traitement de la matrice interspectrale comprennent des moyens d'estimation de la matrice interspectrale débarrassée du bruit, à partir de la matrice de cohérence associée à la matrice interspectrale. IO