Brevet FR3056761A1 — publié le 20180330

RÉSEAU DE MAGNÉTOMÈTRES VECTORIELS ET PROCÉDÉ ASSOCIÉ DE CALIBRATION DES COUPLAGES ENTRE MAGNÉTOMÈTRES DESCRIPTION DOMAINE TECHNIQUE Le domaine de l'invention est celui de l'imagerie des champs biomagnétiques, et concerne plus particulièrement les réseaux de magnétomètres vectoriels utilisés notamment en magnétocardiographie ou en magnétoencéphalographie. ÉTAT DE LA TECHNIQUE ANTÉRIEURE Le principe de mesure des magnétomètres atomiques vectoriels est basé sur l'utilisation d'au moins un harmonique parmi plusieurs générés au moyen d'une source d'excitation. Ces magnétomètres emploient une cellule remplie d'un gaz d'atomes, une source laser qui émet un faisceau polarisé ainsi qu'un photo-détecteur apte à délivrer un signal de sortie représentatif du faisceau ayant traversé la cellule. Une bobine entoure la cellule, alimentée par un générateur de fréquence pour générer un champ magnétique Û) d'excitation, sinusoïdal de fréquence 27r , perpendiculaire à la direction de polarisation des photons issus du laser s'ils sont polarisés rectilignement et perpendiculaire au faisceau laser s'ils sont polarisés circulairement. La projection du champ magnétique sur la direction du champ d'excitation sinusoïdale est ainsi mesurée. Le signal observée par le photo-détecteur contient des résonances (2q +1)co impaires aux fréquences 27r (avec q un entier positif ou nul) qui sont proportionnelles au champ magnétique à mesurer. Ces résonances impaires, dont seule la première est en réalité utilisée, s'annulent pour un champ magnétique total nul. L'observation de la première résonance permet ainsi de mesurer un champ magnétique Bo via l'asservissement d'un champ magnétique de compensation Bc dont l'amplitude est ajustée afin que la somme Bc + Bo soit maintenue à zéro en permanence. De cette manière, 3056761 2 la connaissance du champ de compensation permet de connaitre le champ magnétique ambiant Bo puisque Bc = -Bo. Cet asservissement en champ nul est représenté sur la figure 1. Le magnétomètre 1 comprend une source d'excitation 2 permettant la résonance, par 5 exemple une bobine entourant la cellule et alimentée par un générateur de fréquence, et un photo-détecteur 3 permettant de mesurer l'amplitude d'un signal de sortie à un harmonique de la fréquence d'oscillation de la source d'excitation. Cette amplitude, et la valeur de courant I, calculée à partir de cette amplitude au moyen d'un intégrateur 4, s'annulent pour un champ magnétique total BT nul. On applique pour ce faire le champ 10 magnétique de compensation Bc au moyen d'une bobine 5 alimentée par un courant Ic. La mesure du courant le injecté dans la bobine 5 permet de connaitre, via la fonction de transfert courant/champ magnétique de la bobine 5, la valeur du champ magnétique de compensation Bc et donc celle du champ ambiant Bo. On a représenté sur la figure 2 la courbe de résonance représentant 15 l'amplitude V du signal mesuré par le photo-détecteur 3 en fonction du champ magnétique total BT. Cette courbe présente une portion linéaire lorsque le champ magnétique total est faible et au voisinage de laquelle l'amplitude est directement proportionnelle au champ magnétique selon une pente de résonance. On comprend que l'asservissement en champ nul permet d'augmenter la dynamique de mesure en venant opérer le magnétomètre dans 20 cette portion linéaire, y compris lorsque le champ magnétique à mesurer n'est pas suffisamment faible. L'asservissement décrit en référence à la figure 1 est nécessaire pour chaque axe de mesure des magnétomètres. Il faut ainsi trois asservissements de ce type, afin de mesurer les trois composantes du champ ambiant Bo. Or la superposition au niveau d'un seul magnétomètre de plusieurs champs magnétiques d'excitation à des fréquences différentes pour couvrir les différents axes du magnétomètre dégrade la pente de résonance. Et dans le cadre d'un réseau de magnétomètres, les amplitudes et directions des champs d'excitation générés par un magnétomètre vont également modifier celles de leurs proches voisins. Ces modifications doivent être prises en compte pour générer des 3056761 3 champs magnétiques d'excitation ayant les bonnes amplitudes et directions pour chaque magnétomètre. Par ailleurs, toujours dans le cadre d'un réseau de magnétomètres, le champ magnétique de compensation Bc généré par les bobines d'un magnétomètre est 5 également vu par les autres magnétomètres du réseau, les plus proches en particulier. Il vient perturber les mesures effectuées par les autres magnétomètres du réseau et peut également conduire à des instabilités des magnétomètres. Ainsi, la mise en réseau de magnétomètres vient impacter les champs magnétiques d'excitation ainsi que la fonction de transfert des magnétomètres asservis. 10 EXPOSÉ DE L'INVENTION L'invention a pour objectif de quantifier cet impact lors d'une phase de calibration d'un réseau de magnétomètres vectoriels à pompage optique. Elle vise ainsi à proposer une solution pour estimer les couplages entre les magnétomètres du réseau. La connaissance de ces couplages permettrait alors, lors d'une phase de mesure, de corriger 15 l'amplitude et la direction des champs magnétiques d'excitation des résonances paramétriques, et donc d'améliorer la précision des magnétomètres. Cette connaissance permettrait également de compenser numériquement les biais induits par les champs magnétiques de compensation des magnétomètres employés en mode asservi. Afin de déterminer ces couplages qui vont permettre de limiter la 20 dégradation de la précision des magnétomètres liée à leur utilisation en réseau, l'invention propose un procédé de calibration d'un réseau de N magnétomètres vectoriels à pompage optique, qui comprend une étape de détermination d'une matrice de coefficients de couplage entre magnétomètres du réseau. Cette étape de détermination comporte les opérations suivantes : 25 génération séquentielle, par l'un des magnétomètres du réseau dit magnétomètre de référence, de m champs magnétiques de référence d'amplitudes connues et de directions distinctes ; mesure, par chacun des N magnétomètres du réseau, des m champs de référence suivant n axes ; 3056761 4 détermination de coefficients de couplage entre le magnétomètre de référence et chacun des N magnétomètres du réseau, ladite détermination comprenant le calcul de coefficients de couplage en effectuant les ratios entre les N mesures d'un des m champs de référence suivant l'un des n axes et l'amplitude dudit un 5 des m champs de référence. Certains aspects préférés mais non limitatifs de ce procédé sont les suivants : - la détermination des coefficients de couplage comprend en outre la mesure d'une incertitude portant sur chaque coefficient de couplage calculé et une validation 10 d'un coefficient de couplage calculé lorsque la mesure de l'incertitude portant sur le coefficient de couplage calculé est inférieure à un seuil ; - lorsque la mesure de l'incertitude portant sur un coefficient de couplage calculé est supérieure audit seuil, l'amplitude dudit un des m champs de référence est augmentée et on procède à la réitération du calcul du coefficient de couplage en 15 effectuant le ratio entre la mesure dudit un des m champs de référence suivant l'un des n axes et l'amplitude augmentée dudit un des m champs de référence ; - il comprend l'estimation d'une contribution au couplage des coefficients de couplage calculés dont la mesure d'incertitude est supérieure audit seuil, et l'amplitude dudit un des m champs de référence n'est augmentée (A++) que si 20 ladite estimation de la contribution au couplage est supérieure à un seuil ; - les opérations de génération séquentielle, de mesure et de détermination sont réitérées en adoptant à chaque itération un nouveau magnétomètre du réseau en tant que magnétomètre de référence ; - l'étape de détermination de la matrice de coefficients de couplage est réitérée, 25 après modification des amplitudes et directions des m champs magnétiques de référence générés par chacun des magnétomètres du réseau lorsqu'il est opéré comme magnétomètre de référence. 3056761 5 BRÈVE DESCRIPTION DES DESSINS 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 5 dessins annexés sur lesquels ; la figure 1, déjà décrite ci-dessus, illustre le principe de l'asservissement en champ magnétique nul mis en oeuvre en phase de mesure ; la figure 2, déjà décrite ci-dessus, illustre la portion linéaire de la courbe de résonance observable en présence d'un champ magnétique faible ; 10 la figure 3 est un schéma illustrant différentes étapes d'un procédé de calibration conforme à l'invention. EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS L'invention porte sur un appareil de mesure de champ magnétique, notamment un appareil doté de N (entier naturel au moins égal à deux) magnétomètres 15 vectoriels à pompage optique agencés en réseau et trouvant en particulier application dans le domaine médical pour l'imagerie des champs biomagnétiques, par exemple en magnétoencéphalographie ou en magnétocardiographie. Dans ces applications, les spécifications de précision sur la mesure magnétique faite par chaque magnétomètre du réseau, liées au besoin de localisation des signaux magnétiques qu'il mesure, peuvent être 20 de quelques millimètres. L'invention porte plus particulièrement sur un appareil configuré pour permettre la mise en oeuvre 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 qu'un champ généré par l'un des magnétomètres vient introduire 25 à la mesure réalisée par un autre magnétomètre. Le réseau de magnétomètres vectoriels à pompage optique est typiquement porté par une structure, par exemple conformable, destinée à être portée par un utilisateur. En magnétoencéphalographie, la structure prend la forme d'un casque à même d'être installée sur le crâne de l'utilisateur. 3056761 6 Du fait de la présence du réseau et en utilisation usuelle de celui-ci, 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 BG1 généré par un magnétomètre d'indice i (en réalité par une bobine d'un magnétomètre) induit un champ 5 Bm, au niveau d'un magnétomètre d'indice j (en réalité sur un axe de mesure d'un B. magnétomètre) selon un coefficient de couplage BGt On cherche à déterminer, lors d'une phase de calibration préalable à une utilisation usuelle du réseau de magnétomètres, la matrice de couplage constitué des différents coefficients de couplage. Cette matrice est de taille 3N*3N, où N correspond au 10 nombre de magnétomètres dans le réseau et le chiffre 3 illustre les trois axes des magnétomètres. On a ainsi BMI = Cu.BG1 , 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 Bm, intégrant chacune la contribution de chacun des trois champs BGi). 15 Lors d'une phase de mesure ultérieure à cette phase de calibration, on pourra, en ayant recours à une inversion de la matrice de couplage, 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 phase de calibration, ou phase d'estimation du couplage entre magnétomètres du réseau, sont décrites par la suite en 20 référence à la figure 3. Dans cette phase de calibration, les magnétomètres sont de préférence opérés en boucle ouverte, i.e. sans l'avertissement en champ nul décrit précédemment qui correspond à une opération en boucle fermée. En boucle ouverte, il n'y a pas de champ de compensation et les magnétomètres mesurent le champ magnétique ambiant Bo 25 directement via la mesure de l'amplitude V de la résonance et grâce à la connaissance de la pente de la portion linéaire de la courbe de résonance. Ainsi, l'amplitude du champ mesurée doit être suffisamment faible pour que la résonance s'opère dans sa portion linéaire, soit en pratique ±36 nT pour de l'hélium 4, muni d'un temps de relaxation de 1 ms 3056761 7 et d'un rapport gyromagnétique de 27c x 28.10-9 rad.s-1.T-1. Bien entendu, les magnétomètres peuvent supporter des champs plus grands que 36 nT, mais la mesure ne doit pas être faite dans ces conditions, sauf à prévoir un champ de compensation constant pour se ramener dans la portion linéaire de la résonance comme cela est décrit dans la 5 demande de brevet EP 2 998 758 A2. Dans la phase de calibration, les magnétomètres peuvent alternativement être opérés en boucle fermée, mais il est alors nécessaire pour éviter les interférences de les utiliser l'un après l'autre en phase de mesure d'un champ de référence, ce qui s'avère N fois plus long à réaliser. 10 On a vu précédemment que chaque magnétomètre du réseau comprend au moins une bobine, trois typiquement, qui permet de générer un champ magnétique d'amplitude et de direction connues. Dans la phase de calibration, chaque bobine de chaque magnétomètre génère tour à un champ magnétique d'amplitude connue, pendant que tous les magnétomètres (y compris celui qui dont la bobine génère le champ 15 magnétique) mesurent le champ magnétique résultant, afin d'établir la matrice de couplage et ses 3N x 3N coefficients. Le procédé de calibration du réseau de N magnétomètres vectoriels à pompage optique comprend ainsi une étape de détermination d'une matrice de coefficients de couplage entre magnétomètres du réseau, cette étape comportant les 20 opérations suivantes, mises en oeuvre pour l'un des magnétomètres du réseau, dit magnétomètre de référence, d'indice i. Le magnétomètre de référence est opéré pour réaliser la génération séquentielle GENi de m champs magnétiques de référence d'amplitudes connues et de directions distinctes, m étant un entier au moins égal à un, typiquement égal à trois. Ainsi 25 les bobines du magnétomètre de référence sont utilisées à tour de rôle pour générer un champ BGI. L'ensemble des N magnétomètres du réseau, et donc y compris le magnétomètre de référence, vient procéder à la mesure séquentielle MESj des m champs de référence suivant n axes, n étant un entier au moins égal à un, typiquement égal à trois. 3056761 8 A titre d'exemple, les magnétomètres mesurent simultanément le champ issu de la sinusoïde d'amplitude BG, par exemple à 10 Hz à l'aide d'une détection synchrone à 10 Hz, pour obtenir précisémentBmj . Le bruit des magnétomètres est typiquement de l'ordre 100 fT/VHz entre 0.1 et 300 Hz. Un filtre passe-bas d'ordre 4 à 1 Hz 5 peut être utilisé en sortie de détection synchrone pour obtenir un bruit efficace voisin de 100 fT. Puis, à partir des mesures réalisées simultanément par tous les axes de tous les magnétomètres du réseau lors du créneau temporel dédié à la génération d'un des m champs de référence par le magnétomètre de référence, on vient procéder à la 10 détermination CALCij de coefficients de couplage entre le magnétomètre de référence et chacun des N magnétomètres du réseau. Cette détermination comprend le calcul de B. coefficients de couplage C, = en effectuant les N ratios entre les N mesures Bm, dudit un des m champs de référence suivant l'un des n axes et l'amplitude BG1 dudit un des m champs de référence. Ce calcul est bien entendu effectué pour les mesures dudit un des m 15 champs de référence suivant chacun des n axes. Une fois ces coefficients calculés, on peut venir réitérer les opérations de génération séquentielle GENi, de mesure MESj et de détermination CALCij en adoptant à chaque itération un nouveau magnétomètre du réseau en tant que magnétomètre de référence, comme cela est illustré par le bloc incrémental « i++ » sur la figure 3. Le procédé 20 est arrêté lorsque tous les coefficients de couplage ont été calculés, i.e. après que chaque magnétomètre du réseau ait été utilisé en tant que magnétomètre de référence. Dans un mode de réalisation possible de l'invention, lors de l'opération de détermination CALCij des coefficients de couplage, si l'amplitude du champ mesuré Bm, est assez grande pour que l'incertitude sur l'estimation d'un coefficient de couplage C,, soit 25 faible, alors la valeur de ce coefficient est validée. Sinon, ce coefficient C,, est calculé après avoir augmenté l'amplitude du champ de référence BG1 pour permettre de réduire l'incertitude de mesure sur Bm, et donc l'incertitude sur C,,. 3056761 9 Ainsi, pour des raisons de précision de mesure des coefficients de couplage et de dynamique de mesure des capteurs en boucle ouverte, l'amplitude des champs magnétiques à générer BG1 peut être augmentée à mesure que les coefficients de couplage concernent des magnétomètres de plus en plus éloignés. En effet, l'amplitude de 5 BG1 doit être faible vis-à-vis de la dynamique de mesure des magnétomètres opérés en boucle ouverte, mais doit être par ailleurs suffisamment grande pour que l'incertitude sur les champs magnétiques mesurés Eim, soit compatible avec l'objectif de précision sur le calcul des C. Ces deux contraintes sont parfois antinomiques, et l'invention propose de faire varier l'amplitude de l'amplitude de BGI, en adoptant une valeur initiale faible, par 10 exemple de 10 nT, et en venant l'augmenter pour permettre le calcul de coefficients de couplage suffisamment certains pour des magnétomètres lointains. Ainsi, dans ce mode de réalisation, la détermination des coefficients de couplage comprend en outre la mesure d'une incertitude INC ? portant sur chaque coefficient de couplage calculé pour ledit un des m champs de référence et une validation 15 d'un coefficient de couplage calculé lorsque la mesure de l'incertitude portant sur le coefficient de couplage calculé est inférieure à un seuil. Lorsque tous les coefficients de couplage sont validés, pour chacun des m champs de référence, (« N » en sortie du bloc « INC ? » sur la figure 3) on procède au choix d'un nouveau magnétomètre de référence (bloc incrémental « i++ ») et à la réitération des différentes opérations décrites ci-dessus. 20 En revanche, lorsque la mesure de l'incertitude portant sur un coefficient de couplage calculé est supérieure audit seuil (« O » en sortie du bloc « INC ? » sur la figure 3), on procède à l'augmentation « A++ » de l'amplitude dudit un des m champs de référence, par exemple d'un facteur 10, et à la réitération du calcul du coefficient de couplage en effectuant le ratio entre une nouvelle mesure dudit un des m champs de 25 référence suivant l'un des n axes et l'amplitude augmentée dudit un des m champs de référence. La mesure de l'incertitude portant sur un coefficient de couplage calculé AC comprend le calcul d'une incertitude absolue 4C,J et/ou d'une incertitude relative 'I . Si l'une et/ou l'autre de ces incertitudes est supérieur à un seuil (seuil absolu 4C,J. pour 3056761 10 l'incertitude absolue ; seuil relatif ARCIJ pour l'incertitude relative), on procède à l'augmentation A++ de l'amplitude du champ de référence. L'incertitude absolue ACIJ est par exemple calculée selon m ACE - AB+ R BGy avec AB = bBM + RB où Bpi correspond à l'amplitude BGI 5 du champ de référence, Bm, correspond à l'amplitude du champ mesuré, ARBG correspond à l'incertitude relative sur l'amplitude du champ de référence (par exemple égale à 10-3), ARBM correspond à l'incertitude relative sur l'amplitude du champ mesuré (par exemple égale à 10-1 et bBM correspond au bruit efficace affectant les magnétomètres pendant la durée de mesure choisie (bBM est par exemple de 100 fT). Le seuil absolu ACIJ. pour 10 l'incertitude absolue est par exemple de 104, et le seuil relatif ARCIJ pour l'incertitude relative est par exemple de 5.10-3. Dans une variante de réalisation possible, l'amplitude du champ de référence n'est pas augmentée si les coefficients de couplage qui n'ont pas pu être calculés avec suffisamment de précision présentent une contribution négligeable, i.e. ils sont si 15 faibles qu'ils illustrent une contribution négligeable du champ généré au niveau du champ mesuré, typiquement parce que les magnétomètres sont distants et donc de faible couplage. Dans cette variante, qui est représentée sur la figure 3, le procédé comprend l'estimation d'une contribution au couplage CONT ? portant sur les coefficients 20 de couplage calculés dont la mesure d'incertitude est supérieure audit seuil. L'amplitude dudit un des m champs de référence n'est alors augmentée « A++ » que si ladite estimation de la contribution au couplage est supérieure à un seuil (sortie « O » du bloc « CONT ? »). Dans le cas contraire (sortie « N » du bloc « CONT ? »), on n'augmente pas l'amplitude du champ de référence, on affecte une valeur prédéterminée (par exemple zéro) aux 25 coefficients de couplage incertains, on procède au choix d'un nouveau magnétomètre de référence (bloc incrémental « i++ ») et à la réitération des différentes opérations décrites ci-dessus. 3056761 11 Afin de fixer des ordres de grandeurs, on peut relever que le champ magnétique décroit avec le cube de la distance, si bien que seuls les couplages avec les premiers voisins contribuent réellement. La Demanderesse a pu calculer, pour un réseau 2D de magnétomètres agencés selon une matrice de 10*10 magnétomètres dont le pas est 5 de 3 cm, que les 8 premiers voisins (cardinaux et diagonaux) d'un magnétomètre central fournissent 91% de l'information (poids cumulé des des 8 premiers voisins) qui permet de calculer le champ magnétique vu au niveau du magnétomètre central. En incluant la deuxième couronne et ses 16 magnétomètres supplémentaires, on atteint 95 %. Dans cette variante, on vient par exemple estimer la contribution possible 10 des coefficients de couplage incertains, et la comparer à un seuil. Cette estimation est par exemple réalisée en comptant le nombre Ncci de coefficients de couplage incertains et en leur affectant la plus faible valeur des coefficients de couplage certains et validés Cimin (les coefficients de couplage incertains sont nécessairement plus petits que cette valeur). La contribution possible des coefficients de couplage incertains peut alors s'exprimer comme 15 Ncci*Cimin. Le seuil peut correspondre à l'incertitude sur le coefficient de couplage le plus important Cijmax dans la matrice de couplage, soit ARC,. C. La valeur de ce coefficient de couplage le plus important Cijmax peut être prédéterminée : elle est par exemple égale à 1, ce qui correspond à la diagonale de la matrice de couplage (c'est-à-dire un auto-cou plage). 20 L'étape de détermination de la matrice de coefficients de couplage, dont on a décrit ci-dessus les différentes opérations, peut être réitérée après modification des amplitudes et directions des m champs magnétiques de référence générés par chacun des magnétomètres du réseau lorsqu'il est opéré comme magnétomètre du réseau. En effet, au niveau de chaque magnétomètre, la direction et l'amplitude des champs d'excitation 25 sont modifiées par la présence des magnétomètres voisins. L'impact sur la mesure est du deuxième ordre. Pour être sûr que cet impact soit négligeable, on vient donc actualiser l'amplitude et la direction des champs de référence des magnétomètres en fonction des coefficients de couplage déterminés au cours d'une itération de l'étape de détermination de la matrice de coefficients de couplage et les utiliser au cours d'une itération suivante de 30 cette étape. 3056761 12 Dans un mode de réalisation possible, si l'écart entre les deux itérations conduit à une évolution de la correction apportée aux champs de référence négligeable au regard de la précision visée, on considère que les coefficients de couplage calculés sont connus avec suffisamment de précision et qu'il n'est pas nécessaire de les mesurer de 5 nouveau en actualisant une nouvelle fois les valeurs des champs de référence. Dans ce mode de réalisation, le procédé peut comprendre une mesure d'écart entre les coefficients de couplage calculés entre deux itérations successives de l'étape de détermination de la matrice de coefficients de couplage, et un arrêt de la réitération de cette étape lorsque la mesure d'écart est inférieure à un seuil. 10 Cette mesure d'écart peut être, pour chaque coefficient de couplage calculé, une mesure absolue et/ou une mesure relative. La mesure absolue s'exprime par exemple selon C C et la mesure relative s'exprime par exemple selonC, yref Cu où Cu et Cijref correspondent aux valeurs d'un coefficient de couplage déterminée au cours de deux itérations successives de l'étape de détermination de la matrice de couplage. S'il 15 existe un couple i, j pour lequel l'une et/ou l'autre de ces mesures est supérieur à un seuil (seuil absolu 4C,J,,, pour la mesure absolue ; seuil relatif ARC'y pour l'incertitude relative), on réitère l'étape de détermination de la matrice de couplage après actualisation des valeurs des champs de référence. On notera que pour les premières valeurs de Cijref et l'initialisation de la 20 matrice de couplage, on peut choisir une matrice de zéro ou une matrice construite grâce à une connaissance a priori de la localisation relative des magnétomètres. L'invention s'étend également à un appareil de mesure de champ magnétique apte à mettre en oeuvre le procédé précédemment décrit, et notamment à un appareil comprenant un réseau de N magnétomètres vectoriels à pompage optique et un 25 calculateur, les magnétomètres du réseau étant apte à mettre en oeuvre les opérations de génération séquentielle GENi, par l'un des magnétomètres du réseau dit magnétomètre de référence, de m champs magnétiques de référence d'amplitudes connues et de directions distinctes, et de mesure MESj, par chacun des N magnétomètres du réseau, des m champs de référence suivant n axes, et le calculateur étant configuré pour calculer CALCij des 3056761 13 coefficients de couplage entre le magnétomètre de référence et chacun des N magnétomètres du réseau en effectuant les ratio entre les N mesures d'un des m champs de référence suivant l'un des n axes et l'amplitude dudit un des m champs de référence. 5 REVENDICATIONS1. Procédé de calibration d'un réseau de N magnétomètres vectoriels à pompage optique, comprenant une étape de détermination d'une matrice de coefficients de couplage entre magnétomètres du réseau, ladite étape de détermination comportant les opérations suivantes : génération séquentielle (GENi), par l'un des magnétomètres du réseau dit magnétomètre de référence, de m champs magnétiques de référence d'amplitudes connues et de directions distinctes ; mesure (MESj), par chacun des N magnétomètres du réseau, des m champs de référence suivant n axes ; détermination (CALCij) de coefficients de couplage entre le magnétomètre de référence et chacun des N magnétomètres du réseau, ladite détermination comprenant le calcul de coefficients de couplage en effectuant les ratios entre les N mesures d'un des m champs de référence suivant l'un des n axes et l'amplitude dudit un des m champs de référence. 2. Procédé selon la revendication 1, dans lequel la détermination des coefficients de couplage comprend en outre la mesure d'une incertitude (INC ?) portant sur chaque coefficient de couplage calculé et une validation d'un coefficient de couplage calculé lorsque la mesure de l'incertitude portant sur le coefficient de couplage calculé est inférieure à un seuil. 3. Procédé selon la revendication 2, comprenant, lorsque la mesure de l'incertitude portant sur un coefficient de couplage calculé est supérieure audit seuil, l'augmentation (A++) de l'amplitude dudit un des m champs de référence et la réitération du calcul du coefficient de couplage en effectuant le ratio entre la mesure dudit un des m champs de référence suivant l'un des n axes et l'amplitude augmentée dudit un des m champs de référence. 3056761 15 4. Procédé selon la revendication 3, comprenant l'estimation d'une contribution au couplage (CONT ?) des coefficients de couplage calculés dont la mesure d'incertitude est supérieure audit seuil, et dans lequel l'amplitude dudit un des m champs de référence n'est augmentée (A++) que si ladite estimation de la 5 contribution au couplage est supérieure à un seuil. 5. Procédé selon l'une des revendications 1 à 4, dans lequel les opérations de génération séquentielle, de mesure et de détermination sont réitérées en adoptant à chaque itération (i++) un nouveau magnétomètre du réseau en tant que 10 magnétomètre de référence. 6. Procédé selon la revendication 5, comprenant en outre la réitération de l'étape de détermination de la matrice de coefficients de couplage, après modification des amplitudes et directions des m champs magnétiques de référence générés par 15 chacun des magnétomètres du réseau lorsqu'il est opéré comme magnétomètre de référence. 7. Procédé selon la revendication 6, comprenant en outre une mesure d'écart entre les coefficients de couplage calculés lors de deux itérations successives de 20 l'étape de détermination de la matrice de coefficients de couplage, et un arrêt de la réitération de l'étape de détermination de la matrice de coefficients de couplage lorsque la mesure d'écart est inférieure à un seuil. 8. Appareil de mesure de champ magnétique, comprenant un réseau de N 25 magnétomètres vectoriels à pompage optique et un calculateur, les magnétomètres du réseau étant apte à mettre en oeuvre les opérations de génération séquentielle (GENi), par l'un des magnétomètres du réseau dit magnétomètre de référence, de m champs magnétiques de référence d'amplitudes connues et de directions distinctes, et de mesure (MESj), par chacun des N magnétomètres du réseau, des m champs de référence suivant n axes, et le calculateur étant configuré pour calculer (CALCij) des coefficients de couplage entre le magnétomètre de référence et chacun des N 3056761 16 magnétomètres du réseau en effectuant les ratio entre les N mesures d'un des m champs de référence suivant l'un des n axes et l'amplitude dudit un des m champs de référence. 5