FR 3038508 A1 20170113 FR 1556393 A 20150706 1 Procédé de caractérisation d'une activité électrophysiologique Description DOMAINE TECHNIQUE Le domaine technique de l'invention est la caractérisation d'une activité électrophysiologique par détection de champs magnétiques induits par cette activité. Une application est la caractérisation d'une activité cardiaque par magnétocardiographie. ART ANTERIEUR Jusqu'à présent, lorsqu'il s'agit d'évaluer l'activité du myocarde à des fins de diagnostic en cardiologie clinique, l'électrocardiographie (ECG) est la mesure la plus répandue. Cette technique consiste à détecter les potentiels électriques générés par l'activité électrique cardiaque. Mais cette activité électrique induit également des champs magnétiques. Bien qu'étant d'amplitude très faible - de l'ordre de 10 -10 T à 10-14T, à comparer avec l'amplitude du champ magnétique terrestre, d'environ 10-4 T - les fluctuations de ces champs magnétiques peuvent être enregistrées. C'est l'objet de la magnétocardiographie, fréquemment désignée par l'acronyme MCG qui, à l'instar de la magnétoencéphalographie (MEG) pour le cerveau, permet d'analyser l'activité électrique du myocarde. Il s'agit d'une technique d'analyse non invasive, sans contact et n'émettant aucune radiation, ni champ magnétique. Bien que relativement ancienne, cette technique a récemment bénéficié des progrès dans la détection de champs magnétiques, en particulier grâce aux évolutions dans le domaine de la supraconductivité ou de l'isolation magnétique. De ce fait, bien qu'elle demeure en majorité cantonnée aux travaux de recherche, ses applications cliniques se développent. Des études récentes ont montré que la magnétocardiographie pouvait présenter un intérêt clinique majeur pour le diagnostic précoce de certaines pathologies, par exemple les syndromes coronariens aigus, la myocardie ischémique, ou l'arythmie ventriculaire. Par ailleurs, la magnétocardiographie présente certains avantages par rapport à l'électrocardiographie , en particulier une sensibilité plus élevée, pour une spécificité similaire, dans le diagnostic des arythmies et des maladies coronariennes. Elle offre également la possibilité de construire des cartes spatio-temporelles de l'activité cardiaque, permettant de localiser rapidement une zone myocardique fonctionnelle anormale. Un autre avantage est que 3038508 2 la MCG constitue un enregistrement natif de l'activité cardiaque, ne nécessitant pas de mesure de référence. Par ailleurs, la MCG peut être appliquée au diagnostic de pathologies cardiaques du foetus, domaine dans lequel l'ECG n'est pas bien adaptée du fait du confinement des signaux électriques foetaux par les tissus utérins et abdominaux. 5 Toutefois, les signaux produits par MCG sont souvent bruités et des recherches permettant de remédier à cet inconvénient sont effectuées. Dans ce cadre, les inventeurs ont inventé un procédé simple et efficace, pouvant s'appliquer à des signaux issus de magnétomètres afin de permettre une caractérisation fiable de l'activité électrophysiologique étudiée. 10 EXPOSE DE L'INVENTION Un objet de l'invention est un procédé de caractérisation d'une activité électrophysiologique périodique, comprenant les étapes suivantes i) acquisition, par un premier magnétomètre, d'un premier signal, comportant une mesure d'un champ magnétique induit par ladite activité électrophysiologique 15 périodique ; ii) à partir dudit premier signal, extraction d'une composante d'intérêt représentative dudit champ magnétique induit ; iii) identification, dans ladite composante d'intérêt, d'une pluralité de signaux élémentaires successifs, chaque signal élémentaire correspondant à une période de ladite activité 20 électrophysiologique ; iv) prise en compte d'un signal élémentaire de référence ; v) pour chaque signal élémentaire identifié lors de l'étape iii), calcul d'un indicateur représentant une similitude entre ledit signal élémentaire et ledit signal élémentaire de référence ; 25 vi) en fonction de l'indicateur calculé pour chaque signal élémentaire, sélection d'au moins une partie desdits signaux élémentaires. De cette façon, lesdits signaux élémentaires sélectionnés sont considérés comme caractérisant ladite activité électrophysiologique. Selon un mode de réalisation, l'étape i) comporte également l'acquisition, par un deuxième 30 magnétomètre, d'un deuxième signal, mesurant, de préférence, un champ magnétique 3038508 3 représentatif d'un environnement dans lequel sont placés lesdits premier et deuxième magnétomètres. L'étape ii) peut alors comporter une décomposition dudit premier signal en une pluralité de composantes, dont au moins une composante d'intérêt, cette décomposition étant effectuée, 5 en utilisant ledit deuxième signal, par analyse en composantes indépendantes ou par mesure différentielle. Selon un mode de réalisation, l'étape i) comporte l'acquisition, par une pluralité de magnétomètres, de premiers signaux, chaque premier signal comportant une mesure d'un champ magnétique induit 10 par ladite activité électrophysiologique périodique ; - l'étape ii) est réalisée en utilisant lesdits premiers signaux. Selon un mode de réalisation, l'étape iv) est réalisée l'aide d'un ensemble de signaux élémentaires identifiés dans l'étape iii), le signal élémentaire de référence étant obtenu, par exemple, à partir d'une moyenne, ou d'une autre grandeur statistique telle une médiane, des 15 signaux élémentaires formant ledit ensemble. Ledit signal élémentaire de référence peut, par exemple, être une moyenne desdits signaux élémentaires dudit ensemble. Selon un mode de réalisation, le signal élémentaire de référence est préalablement stocké dans une mémoire. Selon un mode de réalisation, l'indicateur calculé lors de l'étape v) comporte, pour chaque signal 20 élémentaire, un coefficient de corrélation entre ledit signal élémentaire et le signal élémentaire de référence. Selon un mode de réalisation, l'indicateur calculé lors de l'étape v) résulte d'une classification de chaque signal élémentaire, ladite classification étant effectuée entre au moins une classe non représentative dudit signal de référence et une classe représentative du signal de référence. 25 Le procédé peut comporter une étape vii) de détermination d'un signal résultant, représentatif des signaux élémentaires sélectionnés lors de l'étape vi), ce signal caractérisant ladite activité électrophysiologique. Ce signal résultant peut par exemple être obtenu par une moyenne des signaux élémentaires sélectionnés, ou par une moyenne de ces signaux pondérée par l'indicateur de similitude associé à chaque signal élémentaire sélectionné, déterminé lors de 30 l'étape v). 3038508 4 L'activité électrophysiologique caractérisée peut être l'activité du myocarde. Un autre objet de l'invention est un support d'enregistrement d'informations, comportant des instructions pour l'exécution d'un procédé tel que décrit-ci-dessus, ces instructions étant aptes à être exécutées par un processeur, en particulier un micro-processeur. Ce micro-processeur 5 peut alors être relié à une mémoire dans laquelle lesdites instructions sont stockées. Un autre objet de l'invention est un dispositif de caractérisation d'une activité électrophysiologique périodique, comportant : - un premier magnétomètre, apte à détecter un champ magnétique induit par ladite activité élecro-physiologique ; 10 - un processeur, configuré pour mettre en oeuvre un procédé tel que précédemment décrit. Le dispositif peut également comporter un deuxième magnétomètre, apte à détecter un champ magnétique dans un environnement dans lequel est placé ledit dispositif. FIGURES 15 La figure 1 représente un dispositif apte à la mise en oeuvre de l'invention. La figure 2 représente les principales étapes d'un procédé selon un mode de réalisation de l'invention. La figure 3 représente les principales étapes d'un procédé selon une variante du mode de réalisation représenté sur la figure 2. 20 La figure 4A représente un résultat expérimental, montrant un signal représentatif d'une activité cardiaque obtenu suite à une acquisition dont la durée est de 10 minutes. La figure 4B représente un résultat expérimental, montrant un signal représentatif d'une activité cardiaque, similaire à celle analysée en lien avec la figure 4A, ce signal étant obtenu suite à une acquisition dont la durée est de 30 secondes. 25 EXPOSE DE MODES DE REALISATION PARTICULIERS La figure 1 représente un dispositif 1 de caractérisation d'une activité électrophysiologique périodique. Le terme activité électrophysiologique désigne une activité électrique d'un organe, par exemple le coeur ou le cerveau. Par périodique, on entend que ladite activité se reproduit 3038508 5 dans le temps, selon une période déterminée. Cette période peut en particulier être comprise entre 100 ms à 5 minutes. Dans l'exemple représenté, on cherche à caractériser l'activité électrique du myocarde, ce terme désignant un tissu musculaire composant le coeur d'un sujet examiné. Pour ce faire, le dispositif 5 comporte un premier magnétomètre 10, disposé à proximité de la peau P dudit sujet, de façon à minimiser la distance entre l'organe O (dans cet exemple le coeur) et le magnétomètre 10. Le premier magnétomètre 10 n'est pas forcément au contact de la peau. Il peut être disposé à quelques millimètres de cette dernière, par exemple 5 mm, et de préférence à une distance inférieure à 1 cm. 10 Le premier magnétomètre 10 est apte à acquérir un premier signal li, représentatif de l'intensité d'un premier champ magnétique Bi, selon un ou plusieurs axes. Dans l'exemple considéré, on détermine l'intensité du champ magnétique selon un axe 57> orienté selon une direction normale à la peau P du sujet examiné. Un tel magnétomètre est par exemple un magnétomètre à SQUID (acronyme de Superconducting Quantum Interference Device) ou un magnétomètre à pompage 15 optique. Par exemple, il s'agit d'un magnétomètre à pompage optique à Hélium 4, ce type de magnétomètre étant décrit dans la demande W02013092383. Le premier magnétomètre 10 est relié à un circuit électronique 14 de mise en forme, assurant l'amplification du premier signal et/ou sa numérisation, puis il est transmis à un processeur 20, soit par liaison filaire, soit par liaison sans-fil, ou encore par l'intermédiaire d'un support 20 d'enregistrement de type CDROM, clef USB... Le processeur 20 est relié à une mémoire 22 comportant des instructions, ces dernières pouvant être exécutées par le processeur 20, pour mettre en oeuvre le procédé représenté sur les figures 2 ou 3, et décrit ci-après. Ces instructions peuvent être sauvegardées sur un support d'enregistrement, lisible par le processeur, de type disque dur, CDROM ou autre type de mémoire. Le processeur peut être relié à une unité 25 d'affichage 24, par exemple un écran. Le dispositif 1 comporte de préférence un deuxième magnétomètre 12, tenu à distance du premier magnétomètre 10. Comme précédemment décrit, le premier magnétomètre a pour fonction de détecter une intensité d'un champ magnétique généré par l'activité physiologique à caractériser. Mais le premier signal mesuré par le premier magnétomètre comporte également 30 une composante parasite liée au champ magnétique ambiant de l'environnement dans lequel est placé le dispositif 1. Le deuxième magnétomètre a pour fonction de ne détecter que, ou du moins très majoritairement, ce champ magnétique ambiant. Il est alors tenu à une distance 3038508 6 suffisante de l'organe générant l'activité électrophysiologique, de telle sorte que son exposition au champ magnétique induit par cette activité est négligeable. Ce deuxième magnétomètre 12 est par exemple placé à une distance supérieure à 10 cm du premier magnétomètre 10. Il est relié au circuit électronique de mise en forme 12 ainsi qu'au 5 processeur 20. Le deuxième magnétomètre est apte à acquérir un deuxième signal 12, représentatif de l'intensité d'un deuxième champ magnétique B2, selon le même axe, ou les mêmes axes, que ceux considérés lors de la mesure du premier champ magnétique. Ainsi, le premier champ magnétique Bi, dont l'intensité est mesurée par le premier magnétomètre 10, est tel que Bi = Ba + Benv où : 10 Ba représente le champ magnétique induit par l'activité électrique de l'organe O. Ce champ est de faible amplitude, typiquement comprise entre 10 -10 T à 10-14 T lorsque l'organe est un coeur humain. Ben, représente le champ magnétique ambiant, dans l'environnement dans lequel sont placés le premier magnétomètre et le deuxième magnétomètre. 15 Dans cet exemple, le deuxième champ magnétique B2, dont l'intensité est mesurée par le deuxième magnétomètre 12, est tel que B2 5:z-% - env , La figure 2 représente les principales étapes d'un procédé selon l'invention. Etape 100 : acquisition des données. Durant cette étape, chaque magnétomètre acquiert un signal dont l'intensité est représentative du champ magnétique au niveau dudit magnétomètre, 20 et en particulier une de ses composantes selon un ou plusieurs axes. Dans cet exemple, où deux magnétomètres sont utilisés, on acquiert deux signaux |z(t) et I2(t) représentatifs de l'intensité respective, à l'instant t, du premier champ magnétique Bi et du deuxième champ magnétique B2 , selon l'axe 57> précédemment évoqué. Etape 110 : prétraitement des données. Le signal délivré par chaque magnétomètre fait de 25 préférence l'objet d'un prétraitement, préalable à l'étape 120 d'extraction d'une composante d'intérêt. Il s'agit tout d'abord d'une numérisation du signal, la fréquence d'échantillonnage temporel étant par exemple de 2 kHz. Ce prétraitement peut également comporter une sous-étape visant à centrer chaque signal par rapport à sa valeur moyenne. Autrement dit, on soustrait de chaque signal, échantillonné en temps, une estimation de la valeur moyenne de ce 30 signal. 3038508 7 Le prétraitement peut également comprendre l'application de filtres fréquentiels, de manière à sélectionner certaines bandes spectrales. Par exemple, on peut appliquer un filtre passe-bande, dont la bande passante s'étend entre 0.5 Hz et 45 Hz, ainsi qu'un filtre de type Notch, pour rejeter une plage spectrale étroite, de quelques Hz, voire moins, située autour de la fréquence 5 50 Hz. Les étapes 110 à 170 peuvent être mises en oeuvre en post-traitement, c'est-à-dire suite à l'acquisition de l'ensemble des données, ou durant l'acquisition des données. Dans ce dernier cas, on dispose d'une période d'apprentissage, au cours de laquelle les données acquises sont utilisées pour établir des grandeurs nécessaires au traitement des données, en l'occurrence la 10 moyenne, la variance, ou la matrice W décrite en lien avec l'étape 120. Le prétraitement peut également comporter une sous-étape de blanchiment des données, c'est-à-dire traiter chaque signal li et 12, de façon à minimiser la covariance entre les signaux provenant de magnétomètres différents, et de rendre la variance de chaque signal égale à 1. /V) Ala suitede cette étapede pétratement on disposed'unvecteurde mesure 1() = I,/(t 15 dont chaque composante est une variable temporelle centrée réduite. Précisons que cette sous-étape de blanchiment est facultative. Elle est réalisée si l'étape 120 comporte la mise en oeuvre d'un algorithme de type Analyse en Composantes Indépendantes, comme décrit ci-dessous. On comprend que ces opérations de pré-traitement effectuées lors de l'étape 110 sont adaptées 20 à la méthode d'extraction d'une composante d'intérêt appliquée au cours de l'étape 120. Etape 120 : extraction d'une composante d'intérêt. Cette étape vise à décomposer chaque terme du vecteur de mesure en une combinaison linéaire de composantes statistiquement indépendantes l'une de l'autre. Une première composante Ci est représentative du champ magnétique Ba induit par l'activité électrophysiologique, et une deuxième composante C2 est 25 représentative du champ magnétique ambiant Ben, préalablement décrit. Cette étape peut être réalisée par un algorithme de type Analyse en Composantes Indépendantes, fréquemment désigné par le terme ICA (Independent Component Analysis). Cette méthode, connue de l'homme du métier, permet une décomposition du vecteur de mesure 1(t) telle que : 1(t) = M x C(t), M étant une matrice de pondération et C(t) étant un vecteur de composantes. Ainsi, 30 IlAt)1 pita M1,21 [Ci (t)1 i P12,1 M2,2] x [C2 (t)] 3038508 8 M1,1, M1,2, M2,1 et M2,2 étant les termes de la matrice de pondération M, Ci(t) et C2(t) désignant lesdites composantes. L'algorithme d'analyse en composantes indépendantes vise à déterminer une matrice d'extraction W, telle que C = WTx I. Ce type d'algorithme est aisément mis en oeuvre à l'aide de 5 logiciels de calculs mathématiques, par exemple le logiciel bien connu Matlab®. Précisons que l'étape 120 ne met pas nécessairement en oeuvre une analyse en composantes indépendantes. Par exemple, il est possible d'effectuer des mesures différentielles à partir des signaux de chaque magnétomètre, de façon à déterminer la première composante C1. Cette dernière peut alors être obtenue sur la base d'une différence entre les signaux délivrés par 10 chaque magnétomètre, en l'occurrence une soustraction du deuxième signal 12 au premier signal Cependant, une analyse en composantes indépendantes est préférable, car elle permet une meilleure discrimination de la composante d'intérêt par rapport au bruit. Ainsi, l'étape 120 permet d'extraire, à partir des premiers et deuxième signaux 11 et 12, une 15 composante d'intérêt, en l'occurrence la première composante C1, représentative de l'activité électrophysiologique à caractériser. Étape 130 : segmentation temporelle. L'activité électrique physiologique étant périodique, la composante d'intérêt Ci extraite lors de l'étape 120 l'est également. L'étape 130 est une étape visant à identifier, sur ladite composante d'intérêt, des périodes temporelles, chaque période 20 temporelle correspondant à une réalisation de l'activité électrique physiologique. Dans cet exemple, appliqué à l'activité du myocarde, à chaque période comprend un maximum local, correspondant à la contraction du muscle ventriculaire. Cette contraction génère trois ondes, communément désignées par les termes respectifs Q, R, S, ou complexe QRS, l'onde R étant celle présentant un pic d'amplitude le plus élevé, et cela dans chaque période. Aussi, en 25 appliquant un seuillage en amplitude sur la composante d'intérêt Ci, il est possible d'identifier des instants tl?, 0, chaque instant tn correspondant au pic d'une onde R d'une n'n'e période, n étant un entier compris entre 1 et N, N désignant le nombre de périodes comprises dans la composante d'intérêt analysée. La valeur du seuil peut être déterminée à priori, ou manuellement, par exemple en fonction 30 d'une représentation graphique du décours temporel de la première composante sur l'écran 24, ou automatiquement, à partir d'une estimation de l'écart type a de la première composante Ci, le seuil correspondant alors à une pondération de cet écart type, par exemple 26 ou 36. 3038508 9 A partir de chaque instant th?, on peut définir un instant de début de période th? - ati et un instant de fin de période th? + 8t2 , ati et 8t2 correspondant à des durées prédéterminées. Par exemple, ati = 8t2 = 300 ms. Ces durées sont déterminées a priori, en se basant par exemple sur le rythme cardiaque du sujet examiné. 5 On extrait ensuite un signal élémentaire Sn, représentant une eme période de Ci, tel que Sn(t = t1 = O) = - ati) et .S",(t = T) = Ci(t.g + 8t2), avec : T = ati + 8t2, T désignant la durée d'un signal élémentaire, déterminée en fonction de la période de l'activité cardiaque, t1 étant l'instant initial de chaque signal élémentaire. 10 Ainsi, chaque signal élémentaire Sn, est défini selon une plage temporelle comprise entre t = =0 et t = T. T est dans cet exemple égal à 600 ms. Cette étape 130 permet de segmenter la composante d'intérêt Ci selon une séquence de N signaux élémentaires Sn, chaque signal élémentaire étant représentatif d'une période de rang n de l'activité électrophysiologique à caractériser. 15 Etape 140 : obtention d'un signal élémentaire de référence. Cette étape consiste à obtenir un signal élémentaire de référence Sref, considéré comme représentatif de l'activité électrophysiologique étudiée. Ce signal élémentaire peut être stocké dans la mémoire 22. Il peut s'agir d'un signal établit sur la base d'un modèle ou d'un signal préalablement mémorisé, correspondant au sujet examiné. 20 De préférence, en lien avec la variante représentée sur la figure 3, le signal élémentaire de référence Sref est obtenu au cours d'une étape 135, en considérant un ensemble de P signaux élémentaires Sn identifiés dans la composante d'intérêt C1. Cet ensemble peut être par exemple constitué par les P premiers signaux élémentaires identifiés, ou par l'ensemble des signaux élémentaires identifiés, P étant un entier inférieur ou égal à N. Le signal élémentaire de 25 référence est alors obtenu en effectuant une moyenne de ces P signaux élémentaires, de telle sorte que : S' f(t) = mean (SS(t)). où mean désigne l'opérateur moyenne. P est de préférence supérieur à 10, voire supérieur à 30. La composante d'intérêt peut comporter une séquence N signaux élémentaires, dont les P 30 premiers signaux élémentaires, dans l'ordre chronologique, sont utilisés pour établir le signal de 3038508 10 référence S'f, les signaux élémentaires suivants étant utilisés pour caractériser l'activité électrophysiologique. Cela permet une caractérisation de l'activité électrophysiologique simultanément à l'acquisition des mesures. P peut être égal à N, auquel cas le signal de référence Sref est progressivement mis à jour au fur et à mesure que le nombre de signaux élémentaires 5 augmente. Il est donc réactualisé en permanence. En variante, la caractérisation est réalisée en post-traitement, auquel cas l'ensemble des signaux élémentaires peuvent être utilisés pour établir le signal de référence. D'autres grandeurs statistiques que la moyenne peuvent être utilisées pour calculer le signal de référence, en utilisant par exemple l'opérateur médiane au lieu de l'opérateur moyenne. 10 Etape 150 : calcul d'un indicateur exprimant une similitude entre chaque signal élémentaire Sr, et le signal de référence Sref. La similitude peut être une corrélation. L'indicateur peut alors être un coefficient de corrélation rn, exprimé par exemple selon l'expression : ti=T rn = ti-t1 V Sn(ti) X Sn(T - ti)ere f (ti) X Sref (7' ti) Sn(ti) X Sref (T - ti) où t, désigne chaque instant résultant de l'échantillonnage temporel des signaux Sn et Sref, ces 15 signaux étant échantillonnés entre les instants ti = ti et ti =T. Le dénominateur de cette expression est un terme facultatif de normalisation. D'autres méthodes que le calcul d'un coefficient de corrélation sont envisageables. On peut par exemple classifier chaque signal élémentaire de façon à ne retenir que ceux les plus corrélés au signal élémentaire de référence. La classification permet de classer chaque signal élémentaire 20 Sn selon au moins une classe représentative du signal de référence Sref et une classe non représentative de ce signal de référence. Une telle classification peut par exemple résulter de l'application d'un algorithme de réduction de la dimensionnalité de chaque signal élémentaire, par exemple par analyse en composantes principales. Il est alors possible d'établir une distance entre chaque signal élémentaire Sn et le signal de référence Sref, dans un espace défini par des 25 vecteurs de base déterminés par l'analyse en composantes principales. Une telle distance constitue un paramètre de classification, selon lequel un signal élémentaire est, ou non, considéré comme appartenant à une classe représentative du signal de référence. A partir de cette classification, l'indicateur de similitude rn associé au signal élémentaire Sn est établi à partir de la distance, auquel cas plus il est faible, plus le signal élémentaire est considéré comme 30 proche du signal de référence. 3038508 11 Alternativement, l'indicateur de similitude peut prendre la valeur 1 lorsque le signal élémentaire est classifié dans la classe représentative du signal de référence ; il prend la valeur 0 dans le cas contraire. Aussi, d'une façon générale, l'indicateur de similitude représente une ressemblance entre le 5 signal élémentaire considéré et le signal de référence. Cet indicateur peut alors être : un coefficient de corrélation, exprimant une corrélation entre le signal élémentaire Su considéré et le signal de référence S'f. une appartenance à une classe représentative du signal de référence S'f, suite à une classification du signal élémentaire Su considéré. 10 Étape 160 : sélection des signaux élémentaires Su en fonction de la valeur de l'indicateur calculé lors de l'étape 150, afin de ne retenir que les signaux élémentaires les mieux corrélés avec le signal de référence, et considérés comme caractérisant l'activité électrophysiologique étudiée. En effet, bien qu'un seuillage en amplitude soit préalablement appliqué au cours de l'étape 130, les inventeurs ont considéré qu'il était préférable de réaliser une sélection, sur la base de 15 l'indicateur de similitude avec le signal de référence. Cela permet d'améliorer la caractérisation de l'activité électrophysiologique, en rejetant des signaux considérés comme non représentatifs de cette activité, par exemple des signaux présentant des artéfacts. La sélection peut se faire sur la base d'un seuillage par rapport à la valeur du coefficient de corrélation ru défini ci-dessus, les signaux élémentaires 4, pour lesquels le coefficient de 20 corrélation ru dépasse un certain seuil rthreshold , étant sélectionnés. La valeur du seuil rthreshold peut être prédéterminée. Elle peut être également ajustée automatiquement, de façon à retenir un certain pourcentage de signaux élémentaires. Étape 170 : caractérisation de l'activité électrophysiologique sur la base des signaux élémentaires sélectionnés. Durant ces étapes, les signaux élémentaires sélectionnés, désignés 25 par Su* peuvent être moyennés, de façon à fournir un signal résultant Suutput que l'on considère comme caractérisant l'activité électrophysiologique. Ainsi, Suutput (0 = me an (Sn* (0), mean désignant l'opérateur moyenne. Le calcul du signal résultant Suutput peut également être réalisé par une moyenne des signaux élémentaires sélectionnés 4, pondérée par l'indicateur de similitude qui leur est associé ru, de façon à moduler le poids de chaque signal élémentaire 30 sélectionné, dans le calcul du signal résultant Suutp ut, en fonction de cet indicateur. 3038508 12 Essais réalisés : Des essais comparatifs ont été effectués, afin de mettre en évidence les avantages engendrés par l'invention. Des acquisitions de signaux magnétiques générés par le myocarde ont été réalisés en mettant en oeuvre, sur un sujet sain, un dispositif tel que représenté sur la figure 1. 5 Le premier magnétomètre est placé à 5 mm de la peau du sujet et le deuxième magnétomètre est placé à une distance de 10 cm du premier magnétomètre. Les essais sont réalisés à l'intérieur d'une enceinte magnétiquement isolée. Lors d'un premier essai, le temps d'acquisition est de 10 minutes, la fréquence d'échantillonnage étant de 2.056 kHz. 10 La figure 4A représente un signal résultant S'tput , caractérisant l'activité électrophysiologique, ce dernier étant établi à partir d'une moyenne obtenue : courbe 1: à partir de 88 signaux élémentaires Sn, identifiés lors du seuillage en amplitude appliqué lors de l'étape 130, sans mise en oeuvre des étapes 140, 150 et 160. Autrement dit, les signaux élémentaires identifiés lors de l'étape 130 sont tous considérés comme 15 représentatifs de l'activité électrophysiologique. courbe 2 : à partir de 119 signaux élémentaires sélectionnés 4, en mettant en oeuvre le procédé représenté sur la figure 3, le signal de référence étant déterminé à partir de l'ensemble des signaux élémentaires sélectionnés. L'axe des abscisses désigne le temps t, l'unité étant la seconde. 20 II est précisé que la valeur du seuil appliqué lors de l'étape 130 a été ajustée entre les acquisitions correspondant respectivement à la courbe 1 et à la courbe 2, de façon à obtenir un signal résultant calculé sur la base d'un nombre comparable de signaux élémentaires. La courbe 2 est mieux résolue. Les amplitudes du complexe QRS sont plus marquées, ce qui permet une estimation plus précise de sa durée. Par ailleurs, les ondes P et T sont plus nettes 25 que sur la courbe 1. L'activité électrophysiologique est donc caractérisée avec plus de précision. La figure 4B représente les résultats obtenus en ne considérant que 30 secondes d'enregistrement. Les courbes sont similaires à celles représentées en lien avec la figure 4A, la seule différence étant le nombre de signaux élémentaires considérés : sur la courbe 1, on a moyenné 11 signaux élémentaires S', identifiés durant 30 secondes, sans mise en oeuvre des 30 étapes 140, 150, 160. Sur la courbe 2, on a moyenné 14 signaux sélectionnés 4, après mise en oeuvre du procédé représenté sur la figure 3. 3038508 13 On remarque que la courbe 2, mettant en oeuvre l'invention, permet une meilleure définition du complexe QRS ainsi que des ondes P et T, en comparaison avec la courbe 1. Ce résultat est remarquable par le faible nombre de signaux élémentaires considérés pour obtenir une représentation satisfaisante de l'activité électrophysiologique. Cela laisse envisager un temps 5 d'acquisition de quelques dizaines de secondes, comparable au temps d'acquisition nécessaire à l'obtention d'un électrocardiogramme exploitable. Outre une meilleure caractérisation de l'activité électrophysiologique, l'invention permet une réduction significative du temps d'acquisition, ce dernier pouvant se limiter à quelques dizaines de secondes. 10 Bien que décrite en lien avec la caractérisation de l'activité du myocarde, l'invention peut être appliquée à d'autres organes, en particulier l'électroencéphalographie ou la caractérisation de l'activité musculaire en réponse à une stimulation, par exemple une stimulation magnétique trans-crânienne du cortex moteur. Par ailleurs, bien que l'exemple décrit soit réalisé sur la base de deux magnétomètres, 15 l'utilisation d'un nombre plus élevés de magnétomètres est tout à fait envisageable, l'étape 120 permettant l'extraction d'au moins une composante d'intérêt à partir des différents signaux mesurés. On peut par exemple disposer d'une pluralité de magnétomètres, chaque magnétomètre mesurant un champ magnétique induit par l'activité électrophysiologique. L'étape 120 peut être réalisée à partir des signaux enregistrés par ces magnétomètres, afin 20 d'extraire au moins une composante d'intérêt représentative de l'activité électrophysiologique à caractériser. REVENDICATIONS1. Procédé de caractérisation d'une activité électrophysiologique périodique, comprenant les étapes suivantes : i) acquisition, par un premier magnétomètre (10), d'un premier signal (la), comportant une mesure d'un champ magnétique induit (Ba) par ladite activité électrophysiologique périodique ; ii) à partir dudit premier signal, extraction d'une composante d'intérêt (C1) représentative du champ magnétique induit (Ba) ) ; iii) identification, dans ladite composante d'intérêt, d'une pluralité de signaux élémentaires successifs (Sn), chaque signal élémentaire correspondant à une période de ladite activité électrophysiologique ; iv) prise en compte d'un signal élémentaire de référence (Sref); v) pour chaque signal élémentaire (Sn) identifié lors de l'étape iii), calcul d'un indicateur (rn) représentant une similitude entre ledit signal élémentaire et ledit signal élémentaire de référence (Sref); vi) en fonction de l'indicateur (rn) calculé pour chaque signal élémentaire (Sn), sélection d'au moins une partie desdits signaux élémentaires, lesdits signaux élémentaires sélectionnés (4) étant considérés comme caractérisant de ladite activité électrophysiologique. 2. Procédé selon la revendication 1, dans lequel l'étape i) comporte également l'acquisition, par un deuxième magnétomètre (12), d'un deuxième signal (12), mesurant un champ magnétique ( 1 - représentatif d'un environnement dans lequel sont placés lesdits premier et deuxième ,- env . magnétomètres. 3. Procédé selon la revendication 2, dans lequel l'étape ii) comporte, une décomposition dudit premier signal (la, l'a) en une pluralité de composantes (C1, C2), dont au moins ladite composante 3038508 15 d'intérêt (C1), cette décomposition étant effectuée, en utilisant ledit deuxième signal (12, l'2), par analyse en composantes indépendantes ou par mesure différentielle. 4. Procédé selon l'une quelconque des revendications 1 à 3, dans lequel l'étape iv) est réalisée l'aide d'un ensemble de signaux élémentaires (Su) identifiés au cours de l'étape iii), le signal 5 élémentaire de référence (S'f ) étant obtenu à partir d'une moyenne des signaux élémentaires formant ledit ensemble. 5. Procédé selon la revendication 4, dans lequel ledit signal élémentaire de référence (S'f ) est une moyenne desdits signaux élémentaires dudit ensemble. 6. Procédé selon l'une quelconque des revendications 1 à 3, dans lequel le signal élémentaire de 10 référence (S'f) est préalablement stocké dans une mémoire (22). 7. Procédé selon l'une quelconque des revendications précédentes, dans lequel l'indicateur calculé lors de l'étape v) comporte, pour chaque signal élémentaire (4,), un coefficient de corrélation (rn) entre ledit signal élémentaire et le signal élémentaire de référence (Sref). 8 Procédé selon l'une quelconque des revendications 1 à 6, dans lequel l'indicateur calculé lors 15 de l'étape v) résulte d'une classification de chaque signal élémentaire (Su), ladite classification étant effectuée entre au moins une classe non représentative dudit signal de référence (Sref) et une classe représentative dudit signal de référence. 9. Procédé selon l'une quelconque des revendications précédentes, comportant une étape vii) de détermination d'un signal résultant (S'tput), calculé en fonction des signaux élémentaires 20 (SU) sélectionnés lors de l'étape vi), ce signal résultant caractérisant alors ladite activité électrophysiologique. 10. Procédé selon l'une quelconque des revendications précédentes, dans lequel ladite activité électrophysiologique est l'activité du myocarde. 11. Support d'enregistrement d'informations, comportant des instructions pour l'exécution d'un 25 procédé selon l'une quelconque des revendications 1 à 10, ces instructions étant aptes à être exécutées par un processeur (20). 12. Dispositif (1) de caractérisation d'une activité électrophysiologique périodique, comportant : 3038508 16 - un premier magnétomètre (10), apte à détecter un champ magnétique (Ba) induit par ladite activité élecro-physiologique ; - un processeur (20), configuré pour mettre en oeuvre un procédé selon l'une des revendications 1 à 8. 5 13. Dispositif selon la revendication 12, comportant également un deuxième magnétomètre (12), apte à détecter un champ magnétique (Benn) d'un environnement dans lequel est placé ledit dispositif (1).