Brevet FR3107122A1 — publié le 20210813

L’invention porte sur un procédé de détermination d’un couplage entre magnétomètres d’un réseau de N magnétomètres, par exemple à pompage optique, où chaque magnétomètre comprend un système d’annulation de champ apte à être activé pour faire fonctionner le magnétomètre en champ nul. Ce procédé comprendune première phase (P1) au cours de laquelle les N magnétomètres sont séparés en N-1 magnétomètres dont le système d’annulation de champ est désactivé et un magnétomètre de mesure dont le système d’annulation de champ est activé. Cette première phase comprend : - la génération (GENj) par les magnétomètres d’une pluralité de champs magnétiques de référence d’amplitudes connues et de directions distinctes, - la mesure (MESi), par le magnétomètre de mesure, du champ magnétique ambiant sur une pluralité d’axesde mesure; - la détermination (CALCij) de coefficients de couplage entre le magnétomètre de mesure et chacun des N magnétomètres à partir de ladite mesure et desdites amplitudes connues. Figure pour l’abrégé: Figure 2 Réseau de magnétomètres fonctionnant en champ nul et procédé associé de calibration des couplages inter-magnétomètres Le domaine de l’invention est celui del’imagerie des champs biomagnétiques, l’invention portant plus particulièrement sur les réseaux de magnétomètres vectoriels utilisés notamment en magnétocardiographie ou en magnétoencéphalographie. L’imagerie des champs magnétiques générés par différents organes du corps humain, notamment le cœur et le cerveau, exploite des matrices de capteurs magnétiques (magnétomètres) pour identifier des informations pertinentes autant du point de vue de la recherche médicale que de celui du diagnostic de différentes pathologies. Ainsi, l’imagerie des courants cérébraux permet l’étude de la dynamique cérébrale rapide de manière non invasive. Elle sert également au diagnostic préopératoire de l’épilepsie, permettant de localiser les zones épileptogènes de manière nettement plus précise que l’électroencéphalographie. Pour réaliser une reconstruction des sources de champ magnétique (courants circulant dans le cœur ou dans le cerveau par exemple), il est nécessaire de disposer de mesures de champ magnétique selon un ou différents axes de mesure ainsi que de la position des capteurs qui les enregistrent. Actuellement, les techniques de magnétoencéphalographie mettent en œuvre souvent plusieurs centaines des capteurs magnétiques de type SQUID qui présentent des niveaux de bruits intrinsèques très réduits mais requièrent un refroidissement cryogénique pour fonctionner. Par ailleurs, ces techniques d’imagerie sont habituellement pratiquées à l’intérieur d’enceintes magnétiquement blindées pour s’affranchir des champs magnétiques extérieurs qui pourraient perturber les mesures. La taille de ces enceintes doit être suffisamment grande pour contenir le récipient cryogénique ce qui aboutit à des fortes contraintes de coût ainsi que sur l’architecture des bâtiments contenant ces dispositifs, contraintes qui nuisent à leur démocratisation. Alternativement, il est possible d’utiliser des magnétomètres à pompage optique qui présentent des niveaux de bruits intrinsèques semblables mais qui n’ont pas besoin d’être refroidis. Ainsi les capteurs peuvent être disposés plus proches de la peau du patient ce qui permet d’améliorer l’amplitude des signaux ainsi que la résolution spatiale des mesures. Le blindage magnétique peut également être réduit, ce qui permet d’envisager une plus forte diffusion de ces techniques dans les centres hospitaliers. Par ailleurs, les magnétomètres à pompage optique permettent d’adapter la position des capteurs à la surface du corps de chaque patient en particulier. Cela est un avantage, mais requiert une calibration de la position des capteurs qui doit être répétée à chaque ajustement du réseau à la morphologie d’un nouveau patient. Il est donc important de pouvoir calibrer efficacement la position des différents magnétomètres les uns par rapport aux autres, ainsi que la position relative de ce réseau par rapport à des éléments de référence positionnés sur le patient. Une bonne partie des paramètres qui sont recherchés pour localiser et calibrer les magnétomètres à pompage optique (notamment le gain et la linéarité) peuvent être stabilisés (au lieu de devoir être mesurés en continu) par une opération des magnétomètres en boucle fermée, c’est-à-dire en rétroagissant sur les bobines de chaque magnétomètre pour créer des champs dits de compensation de manière à ce que chaque magnétomètre opère dans un champ magnétique total (champ ambiant + champ de compensation) nul selon ses différents axes de mesure. Ce mode d’opération en boucle fermée a été utilisé pour des mesures de magnétocardiographie et magnétoencéphalographie comme cela est par exemple rapporté dans la publication suivante : E. Labyt et al., «Magnetoencephalography With Optically Pumped 4 He Magnetometers at Ambient Temperature», IEEE Transactions on Medical Imaging, vol. 38, n° 1, p. 90-98, janv. 2019. Cependant, cette publication met uniquement en œuvre deux capteurs situés suffisamment loin l’un de l’autre pour ne pas être affectés de façon significative par les champs de compensation de l’autre capteur. Or dans le but de réaliser des mesures médicales utiles, il est nécessaire de mettre en œuvre plusieurs dizaines voire centaines de ces capteurs, ce qui oblige à les rapprocher les uns des autres. Un effet indésirable qui apparaît alors est celui d’un couplage entre magnétomètres disposés proches les uns des autres qui sont les uns sujets aux champs de compensation des autres. Ainsi, il n’est plus possible d’obtenir la valeur du champ magnétique préexistant à l’opération des magnétomètres à partir d’une simple lecture des sorties des capteurs. Des procédures permettent toutefois de récupérer ces valeurs de champ à partir de la lecture des sorties des capteurs et d’une matrice qui traduit les couplages entre les différentes bobines de compensation et les magnétomètres du réseau. Le brevet EP 3343 240 B1 décrit ainsi une manière de corriger les couplages dans un réseau de magnétomètres opérant en boucle fermée, ainsi qu’une technique de détermination de la matrice de couplage qui peut être déclinée en deux variantes: l’une consiste à mesurer les couplages avec tous les magnétomètres opérant en boucle ouverte, l’autre de les mesurer avec les différents magnétomètres opérant successivement en boucle fermée. La première de ces méthodes (calibration en boucle ouverte) aboutit aux vrais coefficients de couplage uniquement dans le cas où les opérations sont menées dans un environnement où les champs résiduels sont quasiment nuls. Or les inventeurs ont pu observer expérimentalement que la présence de champs de quelques nanotesla aboutit à la mesure de couplages importants (jusqu’à 17%) là où les couplages réels sont inférieurs à 1%. Cet artefact semble provenir d’un effet de deuxième ordre par lequel les magnétomètres à résonance paramétrique et à effet Hanle sont sensibles non seulement aux champs selon leurs axes de mesure nominaux (par exemple B z pour l’axe z), mais également à de termes croisés tels que B x .B y . Par ailleurs, pour le cas des magnétomètres à résonances paramétriques et ainsi qu’il est décrit dans le brevet susmentionné, le couplage entre les bobines d’un capteur et l’élément sensible d’un autre capteur proche du premier aboutit également à un inconvénient supplémentaire. Dans cette configuration de magnétomètre, des champs radiofréquence sont appliqués selon un ou plusieurs axes du capteur qui définissent les axes géométriques selon lesquels sont mesurées les différentes composantes du champ magnétique. Les couplages entre deux magnétomètres dont les axes ne sont pas orientés de la même façon aboutit donc à une modification de leurs axes de mesure respectifs ce qui peut avoir un impact très négatif sur la justesse de la reconstruction des sources de champ magnétique. L’invention vise à proposer une méthode de détermination des couplages entre les différents magnétomètres d’un réseau de magnétomètres à pompage optique, par exemple à effet Hanle ou à résonances paramétriques, méthode qui s’affranchisse des inconvénients susmentionnés. A cet effet, l’invention concerne un procédé de détermination d’un couplage entre magnétomètres d’un réseau de N magnétomètres où chaque magnétomètre comprend un système d’annulation de champ apte à être activé pour faire fonctionner le magnétomètre en champ nul. Ce procédé comprendune première phase au cours de laquelle les N magnétomètres sont séparés en N-1 magnétomètres dont le système d’annulation de champ est désactivé et un magnétomètre de mesure dont le système d’annulation de champ est activé. Cette première phase comprend les étapes suivantes: - génération, par chacun des N magnétomètres, d’une pluralité de champs magnétiques de référence d’amplitudes connues et de directions distinctes, - mesure, par le magnétomètre de mesure, du champ magnétique ambiant sur une pluralité d’axesde mesure; - détermination de coefficients de couplage entre le magnétomètre de mesure et chacun des N magnétomètres, ladite détermination comprenant: - la détection d’une contribution d’un champ magnétique de référence dans la mesure du champ magnétique ambiant sur l’un des axes de mesure; - le calcul d’un ratio entre l’amplitude de ladite contribution du champ magnétique de référenceet l’amplitude connue dudit champ magnétique de référence. Certains aspects préférés mais non limitatifs de ce procédé sont les suivants : lesdites étapes de la première phase sont réitérées N en adoptant à chaque itération un nouveau magnétomètre du réseau en tant que magnétomètre de mesure; chaque champ magnétique de référence est porteur d’une information qui lui est propre et la détection d’une contribution d’un champ magnétique de référence comprend l’identification de l’information propre audit champ magnétique de référence; l’information propre à un champ magnétique de référence est une fréquence caractéristique, par exemple une fréquence de base élevée à une puissance d’un nombre donné; les magnétomètres sont des magnétomètres à résonance paramétriques qui comprennent chacun un système d’excitation apte à être activé pour induire des champs radiofréquence d’excitation des résonances paramétriques et au cours de la première phase ledit système d’excitation des N magnétomètres est activé; il comprend une deuxième phase identique à la première phase (P1) sauf en ce que le système d’excitation des N-1 magnétomètres autres que le magnétomètre de mesure est désactivé; la deuxième phase comporte N itérations, chaque itération utilisant un magnétomètre de mesure différent parmi les N magnétomètres. D'autres aspects, buts, avantages et caractéristiques de l’invention apparaîtront mieux à la lecture de la description détaillée suivante de formes de réalisation préférées de celle-ci, donnée à titre d'exemple non limitatif, et faite en référence aux dessins annexés sur lesquels : est un schéma d’un magnétomètre appartenant à un réseau de magnétomètresselon l’invention; est un schéma illustrant différentes étapes d’un procédé de calibration conforme à l’invention. L’invention porte sur un appareil de mesure de champ magnétique, plus particulièrement un appareil doté de N (entier naturel au moins égal à deux) magnétomètres vectoriels agencés en réseau. Cet appareil trouve en particulier application dans le domaine médical pour l’imagerie des champs biomagnétiques, notamment en magnétoencéphalographie ou en magnétocardiographie. L’invention porte plus particulièrement sur un appareil configuré pour permettre la mise en œuvre du procédé de calibration des magnétomètres décrit ci-après. Ce procédé permet d’estimer l’impact de la mise en réseau des magnétomètres, en venant mesurer la perturbation que le fonctionnement de l’un des magnétomètres vient introduire à la mesure réalisée par un autre magnétomètre. Pour appliquer ce procédé, ainsi que plus généralement pour pouvoir opérer avec succès un réseau de magnétomètres en boucle fermée, il est préférable que les magnétomètres de ce réseau soient capables de mesurer les trois composantes du champ magnétique. En effet, dans le cas contraire, une opération en boucle fermée ne permet pas de garantir l’absence de champ résiduel en tout magnétomètre du réseau, condition nécessaire à une mesure des couplages exempte des artefacts qui ont été décrits plus haut. Sans que cela soit limitatif de l’invention, les magnétomètres du réseau sont de préférence des magnétomètres à pompage optique, par exemple des magnétomètres à effet Hanle ou des magnétomètres à résonances paramétriques. La description qui suit prend ainsi l’exemple des magnétomètres à pompage optique. En référence à la , chacun des magnétomètres 20 du réseau comprend une cellule 1 remplie d’un gaz atomique, par exemple de l’hélium-4 ou un gaz alcalin, soumis à un champ magnétique ambiant dont la projection sur trois axes de coordonnées rectangulaires x, y, z définit trois composantes. Le champ magnétique ambiant se décompose ainsi en trois composantes Bx, By et Bz chacune selon l’un des axes de mesure du magnétomètre x, y et z. La cellule 1 est éclairée par une source de pompage optique 2,3 agencée pour émettre en direction de la cellule 1 un faisceau de lumière, par exemple un faisceau laser, accordé à une longueur d’onde de pompage (ce faisceau est ainsi également désigné par faisceau de pompe). La longueur d’onde de pompage est calée sur une raie de transition atomique, par exemple sur la raie D 0 à 1083 nm dans le cas de l’hélium-4. Le faisceau de lumière peut être émis par une source laser 2 et être polarisé linéairement au moyen d’un polariseur rectiligne 3 intercalé entre la source laser 2 et la cellule 1 ou directement intégré à la source laser 2. Le faisceau de lumière se propage selon une direction de propagation confondue avec l’axe x et est polarisé linéairement selon l’axe z. Dans le cas où l’élément sensible est l’hélium-4, le magnétomètre 20 comporte par ailleurs un système de décharge haute fréquence (HF), comprenant un générateur HF 4 et des bobines de surtension 5, pour amener les atomes du gaz atomique dans un état énergisé où ils sont aptes à subir la transition atomique lorsqu’ils sont éclairés par le faisceau de lumière, typiquement dans l’état métastable 2 3 S 1 . Lorsque le magnétomètre 20 est un magnétomètre à résonances paramétriques, il comprend également un système d’excitation apte à être activé pour induire des champs radiofréquenced’excitation des résonances paramétriques. Ce système d’excitation des résonances paramétriques comporte un générateur de radiofréquence 8 qui alimente des bobines d’Helmholtz 7 d’axes orthogonaux qui entourent la cellule afin de générer un champ magnétique d’excitation des résonances paramétriques, également désigné par champ radiofréquence d’excitation. Ce circuit d’excitation vient plus particulièrement générer un champ magnétique radiofréquence présentant deux composantes orthogonales à la direction de polarisation et oscillant chacune à sa propre fréquence d’oscillation, à savoir une composante B ω cosωt selon l’axe x oscillant à la pulsation ω (avec par exemple ω=2π.3000 kHz) et une composante B Ω cosΩt selon l’axe y oscillant à la pulsation Ω (avec par exemple Ω=2π.16 kHz). Ces composantes conduisent à des résonances à chacune des fréquences d’oscillation Ω/2π, ω/2π et à un inter-harmonique des fréquences d’oscillation (ω ± Ω)/2π, ces résonances étant associées aux valeurs du champ ambiant dans les directions x, y et z respectivement. Le magnétomètre 20 comprend par ailleurs un photodétecteur 6 agencé pour recevoir le faisceau de lumière ayant traversé la cellule et un circuit de détection de résonances paramétriques 9 configuré, lorsque le magnétomètre 20 est un magnétomètre à résonances paramétriques, pour réaliser une détection synchrone à un harmonique de chacune des fréquences d’oscillation d’un signal électrique délivré par le photodétecteur et une détection synchrone à un inter-harmonique des fréquences d’oscillation du signal électrique délivré par le photodétecteur. Le dispositif 9 comporte trois voies de détection: une première voie Vx pour la détection du signal à Ω /2π (axe x), une second voie Vy pour la détection du signal à ω/2π (axe y), et une troisième voie Vz pour la détection du signal à Ω±ω/2 π (axe z). Le signal sur chacune des première et deuxième voies Vx, Vy est d’abord amplifié puis filtré avec un filtre passe-bande correspondant à la fréquence centrale adéquate (i.e. correspondant à celle du champ RF appliqué). Le signal obtenu est alors multiplié par un signal de référence et traité par un détecteur synchrone DSx, DSy. La troisième voie Vz utilise deux détections synchrones en série, l’une à ω/2π au moyen du détecteur DSy de la deuxième voie Vy et l’autre à Ω /2π au moyen d’un détecteur synchrone DSz. Le magnétomètre 20 comprend également un système d’annulation de champ apte à être activé pour faire fonctionner le magnétomètre en champ nul. Ce système peut prendre la forme d’un circuit d’asservissement 10 du magnétomètre en boucle fermée. Un tel circuit d’asservissement comprend trois voies d’asservissement Wx, Wy, Wz chacune couplée à une sortie d’une voie de détection correspondante Vx, Vy, Vz. Chacune des voies d’asservissement Wx, Wy, Wz exploite la sortie de la voie de détection correspondante en tant que signal d’erreur pour venir constamment réajuster un champ de compensation. Chacune des voies d’asservissement Wx, Wy, Wz comprend un intégrateur Ix, Iy, Iz configuré pour délivrer un signal de compensation et un générateur de courant GCx, GCy, GCz piloté par le signal de compensation pour injecter du courant dans l’une des bobines de Helmholtz 7 afin de générer un champ magnétique de compensation BCx, BCy, BCz opposé à une composante du champ ambiant Bx, By, Bz. La mesure des courants circulant dans les bobines 7 permet de déduire les champs qu’il est nécessaire d’appliquer pour annuler les différentes composantes du champ ambiant, et donc de disposer de la valeur de ces différentes composantes. Le magnétomètre 20 est par ailleurs doté d’un circuit de génération de champs magnétiques de référence et l’appareil de mesure (i.e. le réseau de magnétomètres) comporte quant à lui un calculateur configuré pour piloter les magnétomètres du réseau et déterminer des coefficients de couplage entre ceux-ci conformément au procédé qui sera décrit par la suite. Le circuit de génération de champs magnétiques de référence peut comprendre un générateur de signaux de référence 11 connecté à chacune des bobines d’Helmholtz 7 en parallèle du circuit de connexion de la bobine au générateur 8 du système d’excitation des résonances paramètres et au circuit d’asservissement 10 du magnétomètre en boucle fermée. Cette connexion du générateur de signaux de référence 11 à chacune des bobines d’Helmholtz 7 peut être réalisée par l’intermédiaire d’une résistance de polarisation de valeur plus élevée que l’impédance de la bobine à basse fréquence ( Du fait de la mise en réseau des magnétomètres et en utilisation usuelle de ce réseau, le champ magnétique mesuré au niveau d’un magnétomètre intègre des contributions de l’ensemble des autres magnétomètres. Ainsi, le champ magnétique B Gi généré par un magnétomètre M i (en réalité par une bobine du magnétomètre M i ) induit un champ B Mj au niveau d’un magnétomètre M j (en réalité sur chacun des axes de mesure du magnétomètre M j ) selon un coefficient de couplage . On cherche à déterminer, lors d’une opération de calibration préalable à une utilisation usuelle du réseau de magnétomètres, une matrice de couplage constitué des différents coefficients de couplage. Cette matrice est de taille 3N x 3N, où N correspond au nombre de magnétomètres dans le réseau et le chiffre 3 illustre les trois axes des magnétomètres. On a ainsi , avec les indices i et j variant de 1 à 3N, chaque composante d’un champ généré par un magnétomètre étant vue sur chacun des axes de mesure des magnétomètres (on a ainsi trois mesures B Mj intégrant chacune la contribution de chacun des trois champs B Gi ). Lors d’une phase de mesure ultérieure à cette phase de calibration, il est possible, en ayant recours à une inversion de la matrice de couplage, de remonter pour chaque magnétomètre au champ magnétique réel, tel qu’il serait mesuré en l’absence du réseau. Les différentes étapes de cette opération de calibration sont décrites par la suite en référence à la . Cette opération exploite les bobines 7 des magnétomètres pour venir générer, au moyen du générateur de signaux de référence 11, des champs magnétiques d’amplitude connue. Le champ résultant est tour à tour mesuré par l’un des magnétomètres du réseau, afin d’établir la matrice de couplage et ses 3N x 3N coefficients. L’opération de calibration comprend une première phase P1 qui comporte N itérations, chaque itération utilisant un magnétomètre de mesure différent parmi les N magnétomètres. Au cours de chaque itération de la première phase P1, les N magnétomètres sont séparés en N-1 magnétomètres dont le système d’annulation de champ est désactivé et un magnétomètre de mesure dont le système d’annulation de champ est activé. Lorsque les magnétomètres sont du type à résonances paramétriques, leur système d’excitation des résonances paramétriques est activé au cours de cette première phase. Ainsi, le magnétomètre de mesure fonctionne en champ nul et réalise une mesure du champ (son système d’annulation de champ est activé), tandis que tous les autres ne mesurent pas et ne compensent donc pas le champ en leur position (leur système d’annulation de champ est désactivé), leur rôle se cantonnant à générer sur leurs bobines des champs de référence ainsi que le cas échéant les radiofréquences correspondantes à leur fonctionnement habituel (leur système d’excitation des résonances paramétriques est activée). Cette manière d’opérer avec le magnétomètre de mesure en boucle fermée garantit que la caractérisation des couplages est réalisée avec ce magnétomètre de mesure opérant en champ nul selon ses trois axes de mesure, ce qui évite l’apparition de nombreux artefacts procédant d’effets inter-axes qui sont nuisibles à la caractérisation des couplages. Par ailleurs, le fait que les circuits d’excitation des résonances paramétriques soient activés au cours de cette caractérisation des couplages permet de mesurer ces couplages en prenant en compte les éventuels désalignements des axes de mesure résultant des couplages entre les radio-fréquences des différents magnétomètres lorsque ceux-ci sont du type à résonances paramétriques. Ces désalignements sont également présents pendant l’opération du réseau de capteurs et il est donc nécessaire de les prendre en compte dans la détermination de la matrice de couplage. Chaque itération de la première phase P1 comprend les étapes suivantes: - génération GENj, par chacun des N magnétomètres, d’une pluralité de champs magnétiques de référence d’amplitudes connues et de directions distinctes, - mesure MESi, par le magnétomètre de mesure, du champ magnétique ambiant sur une pluralité d’axesde mesure; - détermination CALCij de coefficients de couplage entre le magnétomètre de mesure et les N magnétomètres, à partir de ladite mesure du champ ambiant et des amplitudes connues des champs de référence. Au cours de l’étape de génération GENj, le circuit de génération de champs magnétiques de référence de chacun des N magnétomètres peut ainsi être activé de telle manière que chacun des N magnétomètres génère sur chacun de ces trois axes (trois directions distinctes) un champ de référence d’amplitude connue. Et au cours de l’étape de mesure MESi, le magnétomètre de mesure peut venir mesurer le champ magnétique ambiant sur chacun de ses trois axesde mesure. Il peut arriver que l’amplitude d’un champ de référence soit trop forte et provoque des effets indésirables comme par exemple des courants de Foucault dans des éléments métalliques situés à proximité du réseau. Il est alors possible de décroître cette amplitude, par exemple d’un ordre de grandeur par rapport à l’amplitude des autres champs de référence, puis de prendre en compte cette réduction lors du calcul des coefficients de couplage. La détermination CALCij des coefficients de couplage comprend la détection d’une contribution de chacun des champs magnétiques de référence dans la mesure du champ magnétique ambiant sur l’un des axes de mesure. Cette détection est suivie, pour chacun des champs magnétiques de référence, du calcul du coefficient de couplage sous la forme d’un ratio entre l’amplitude de ladite contribution du champ magnétique de référenceet l’amplitude connue dudit champ magnétique de référence. Les coefficients de couplage mesurés peuvent, pour chacun des axes de mesure, être normalisés par l’auto-couplage de de cet axe, i.e. le ratio entre l’amplitude de la contribution sur cet axe du champ de référence généré sur cet axe et l’amplitude connue de ce champ de référence. Etant donné que les coefficients de couplages sont symétriques (i est couplé à j de la même manière que j à i), il est possible d’utiliser les coefficients de couplage calculés dans les deux sens pour identifier d’éventuels problèmes dans le procédé de calibration ou encore de ne calculer que la moitié de ces coefficients. Lors de l’étape de génération GENij, le circuit de génération de champs magnétiques de référence de chacun des magnétomètres peut venir générer des champs de référence chacun porteur d’une information qui lui est propre (i.e. une information propre à un axe de ce magnétomètre), cette information pouvant être identifiée dans la mesure du champ magnétique ambiant réalisée sur l’un des axes du magnétomètre de mesure afin d’isoler la contribution du champ magnétique de référence porteur de cette information. Cette information propre à un champ magnétique de référence peut être une fréquence caractéristique, les signaux de référence générés par le générateur de signaux 11 pouvant par exemple être composées de sinusoïdes situées à des fréquences propres à chaque axe de chaque magnétomètre. Ces fréquences peuvent par exemple être choisies à des harmoniques correspondant une fréquence de base élevée aux puissances d’un nombre donné, soit avec f la fréquence de base et p ledit nombre donné. Ainsi, si la fréquence de base est de 1 Hz et si ledit nombre donné est 2, les fréquences caractéristiques sont des fréquences de 1, 2, 4, 8, 16, 32… Hz. Ce choix des fréquences caractéristiques offre l’avantage que la présence de tout harmonique ne correspondant pas à une puissance dudit nombre donné constitue un bon moyen de diagnostic d’un effet inter-axe non souhaité. A titre d’exemple, un premier magnétomètre peut générer sur son axe x un champ de référence , sur son axe yun champ de référence et sur son axe z un champ de référence , tandis qu’un deuxième magnétomètre génère sur son axe x un signal , sur son axe yun champ de référence et sur son axe z un champ de référence , etc. La détection de la contribution d’un champ magnétique de référence dans la mesure du champ magnétique ambiant sur l’un des axes de mesure du magnétomètre de mesure peut alors comprendre une détection synchrone à la fréquence caractéristique de ce champ de référence. Alternativement, cette information propre à un axe d’un magnétomètre (à un champ de référence généré par ce magnétomètre) peut correspondre à la durée de phases haute et basse de signaux de référence prenant la forme de créneaux. Dans une autre variante, les signaux de référence forment une base orthogonale telle que celle utilisée dans les techniques d’émission radio par étalement de spectre, où chaque magnétomètre signe son émission sur l’un de ses axes par une séquence caractéristique d’émission dans différentes bandes de fréquence, ces bandes étant incluses dans la bande de signaux mesurables par les magnétomètres du réseau. L’étape de mesure MESi est par exemple réalisée séquentiellement sur chacun des axesdu magnétomètre de mesure, avec par exemple une mesure tout d’abord réalisée selon un premier axe du magnétomètre de mesure pendant un temps suffisamment long pour permettre de détecter de façon satisfaisante la contribution de chacun des champs magnétiques de référence dans la mesure du champ magnétique ambiant sur ce premier axe, par exemple pour distinguer des champs de référence présentant des fréquences caractéristiques proches. Un temps de 10 secondes s’avère ainsi très largement suffisant pour éviter toute superposition de raies lorsque la fréquence de base est comme dans l’exemple précédent de 1 Hz. Une procédure similaire est ensuite suivie sur le second axe, puis sur le troisième axe du magnétomètre de mesure. A l’issue d’une itération de la première phase, on procède, comme illustré par le bloc CHANi sur la , à la sélection d’un nouvel magnétomètre de mesure parmi les N magnétomètres pour une nouvelle itération de la première phase. Les opérations précédemment décrites sont ainsi réitérées de manière à être au total réalisées N fois en utilisant à chaque fois un magnétomètre de mesure différentparmi les N magnétomètres En référence à la , le procédé de calibration selon l’invention peut comprendre une deuxième phase P2 qui comporte également N itérations, chaque itération utilisant, comme illustré par le bloc CHANi, un magnétomètre de mesure différent parmi les N magnétomètres. La deuxième phase peut être mise en œuvre avant ou après la première phase. Alternativement, il est possible de mélanger itérations de la première phase et itérations de la deuxième phase, par exemple en faisant suivre une itération de la première phase par une itération de la deuxième phase utilisant le même magnétomètre de mesure. On a vu précédemment que le fait d’activer les radiofréquences d’excitation des résonances paramétriques pendant la première phase permet de prendre en compte les éventuels désalignements des axes de mesure résultant des couplages entre les radiofréquences des différents magnétomètres. En revanche, ces couplages sont nuisibles à l’estimation des orientations et positions des différents magnétomètres du réseau par une méthode de triangulation ou par d’autres méthodes plus élaborées telles qu’une estimation par la méthode des moindres carrées prenant en compte la configuration géométrique théorique des champs magnétiques créés par chaque bobine de chaque magnétomètre. C’est pour cette raison que pour le cas d’un magnétomètre à résonances paramétriques le procédé selon l’invention peut comprendre une deuxième phase P2 identique à la première phase P1 sauf en ce que, à chacune de ses itérations, le système d’excitation des résonances paramétriques des N-1 magnétomètres autres que le magnétomètre de mesure est désactivé. Ainsi, au cours d’une itération de la deuxième phase, les N magnétomètres sont séparés en N-1 magnétomètres dont le système d’excitation est désactivé et le système d’annulation de champ est désactivé et un magnétomètre de mesure dont le système d’excitation est activé et le système d’annulation de champ est activé. En référence à la , cette deuxième phase comprend les étapes suivantes: - génération GENSEj, par chacun des N magnétomètres, d’une pluralité de champs magnétiques de référence d’amplitudes connues et de directions distinctes, - mesure MESi, par le magnétomètre de mesure, du champ magnétique ambiant sur une pluralité d’axesde mesure; - détermination CALCij* de coefficients de couplage entre le magnétomètre de mesure et les N magnétomètres, à partir de ladite mesure du champ ambiant et des amplitudes connues des champs de référence. La matrice de couplage obtenue à l’issue de la deuxième phase P2 ne prend pas en compte les désalignements des axes résultants du couplage entre les radiofréquences des magnétomètres autres que celui de mesure et celui-ci. La matrice des couplages déterminée dans cette deuxième phase sert uniquement à la calibration des positions et orientations des différents magnétomètres. Ainsi les neuf coefficients de couplage entre deux magnétomètres peuvent être utilisés pour établir les trois angles et trois distances entre ces deux magnétomètres comme décrit par exemple dans le brevet EP 3343240 B1. Une comparaison entre les matrices déterminées par chacune des premières et deuxièmes phases permet également de mesurer les désalignements d’axes résultant des couplages radiofréquences. L’invention n’est pas limitée au procédé précédemment décrit mais s’étend comme indiqué précédemment à un appareil de mesure de champ magnétique configuré pour permettre la mise en œuvre de ce procédé. Cet appareil comprend un calculateur et un réseau de N magnétomètres où chaque magnétomètre comprend un système d’annulation de champ apte à être activé pour faire fonctionner le magnétomètre en champ nul. Les N magnétomètres sont notamment aptes à mettre en œuvre les étapes suivantes alors qu’ils sont séparés en N-1 magnétomètres dont le système d’annulation de champ est désactivé et un magnétomètre de mesure dont le système d’annulation de champ est activé: - génération (GENj), par chacun des N magnétomètres, d’une pluralité de champs magnétiques de référence d’amplitudes connues et de directions distinctes, - mesure (MESi), par le magnétomètre de mesure, du champ magnétique ambiant sur une pluralité d’axesde mesure. Le calculateur est quant à lui notamment configuré pour détecter une contribution d’un champ magnétique de référence dans la mesure du champ magnétique ambiant sur l’un des axes de mesureet pour déterminer un coefficient de couplage entre le magnétomètre de mesure et l’un des magnétomètres en calculant un ratio entre l’amplitude de ladite contribution du champ magnétique de référenceet l’amplitude connue dudit champ magnétique de référence. Procédé de détermination d’un couplage entre magnétomètres (20) d’un réseau de N magnétomètres où chaque magnétomètre comprend un système d’annulation de champ (10) apte à être activé pour faire fonctionner le magnétomètre en champ nul, le procédé étant caractérisé en ce qu’il comprendune première phase (P1) au cours de laquelle les N magnétomètres sont séparés en N-1 magnétomètres dont le système d’annulation de champ est désactivé et un magnétomètre de mesure dont le système d’annulation de champ est activé, ladite première phase comprenantles étapes suivantes: - génération (GENj), par chacun du magnétomètre de mesure et des N-1 magnétomètres dont le système d’annulation de champ est désactivé, d’une pluralité de champs magnétiques de référence d’amplitudes connues et de directions distinctes, - mesure (MESi), par le magnétomètre de mesure, du champ magnétique ambiant sur une pluralité d’axesde mesure; - détermination (CALCij) de coefficients de couplage entre le magnétomètre de mesure et chacun du magnétomètre de mesure et des N-1 magnétomètres dont le système d’annulation de champ est désactivé, ladite détermination comprenant: - la détection d’une contribution d’un champ magnétique de référence dans la mesure du champ magnétique ambiant sur l’un des axes de mesure; - le calcul d’un ratio entre l’amplitude de ladite contribution du champ magnétique de référenceet l’amplitude connue dudit champ magnétique de référence. Procédé selon la revendication 1, dans lequel lesdites étapes de la première phase sont réitérées N en adoptant à chaque itération un nouveau magnétomètre du réseau en tant que magnétomètre de mesure. Procédé selon l’une des revendications 1 et 2, dans lequel chaque champ magnétique de référence est porteur d’une information qui lui est propre et dans lequel la détection d’une contribution d’un champ magnétique de référence comprend l’identification de l’information propre audit champ magnétique de référence. Procédé selon la revendication 3, dans lequel l’information propre à un champ magnétique de référence est une fréquence caractéristique. Procédé selon la revendication 4, dans lequel la fréquence caractéristique d’un champ magnétique de référence est une fréquence de base élevée à une puissance d’un nombre donné. Procédé selon l’une des revendications 1 à 5, dans lequel les magnétomètres sont des magnétomètres à résonance paramétriques qui comprennent chacun un système d’excitation (8) apte à être activé pour induire des champs radiofréquence d’excitation des résonances paramétriques et dans lequel au cours de la première phase ledit le système d’excitation des N magnétomètres est activé. Procédé selon la revendication 6, comprenant en outre une deuxième phase (P2) identique à la première phase (P1) sauf en ce que le système d’excitation des N-1 magnétomètres autres que le magnétomètre de mesure est désactivé. Procédé selon la revendication 7, dans lequel dans lequel la deuxième phase comporte N itérations, chaque itération utilisant un magnétomètre de mesure différent parmi les N magnétomètres. Appareil de mesure de champ magnétique, comprenantun calculateur et un réseau de N magnétomètres vectoriels où chaque magnétomètre comprend un système d’annulation de champ (10) apte à être activé pour faire fonctionner le magnétomètre en champ nul, les N magnétomètres étant aptes à mettre en œuvre les étapes suivantes alors qu’ils sont séparés en N-1 magnétomètres dont le système d’annulation de champ est désactivé et un magnétomètre de mesure dont le système d’annulation de champ est activé: - génération (GENj), par chacun du magnétomètre de mesure et des N-1 magnétomètres dont le système d’annulation de champ est désactivé, d’une pluralité de champs magnétiques de référence d’amplitudes connues et de directions distinctes, - mesure (MESi), par le magnétomètre de mesure, du champ magnétique ambiant sur une pluralité d’axesde mesure; et le calculateur étant configuré pour détecter une contribution d’un champ magnétique de référence dans la mesure du champ magnétique ambiant sur l’un des axes de mesureet pour déterminer un coefficient de couplage entre le magnétomètre de mesure et l’un des magnétomètres en calculant un ratio entre l’amplitude de ladite contribution du champ magnétique de référenceet l’amplitude connue dudit champ magnétique de référence.