FR 3080186 A1 20191018 FR 1853319 A 20180416 PROCÉDÉ DE CALIBRATION D'UN RÉSEAU DE MAGNÉTOMÈTRES DESCRIPTION DOMAINE TECHNIQUE Le domaine de l'invention est celui du suivi du mouvement d'un objet magnétique mobile au moyen d'un ensemble de magnétomètres agencés en réseau. L'invention concerne plus particulièrement la calibration du réseau de magnétomètres pour connaître les positions et orientations des magnétomètres dans un repère associé au réseau. ÉTAT DE LA TECHNIQUE ANTÉRIEURE Le suivi du mouvement d'un ou plusieurs objets magnétiques peut être réalisé au moyen d'un ensemble de magnétomètres solidaires d'un même support mécanique de sorte que leurs positions et orientations relatives sont fixes dans le temps. Les magnétomètres sont ainsi agencés en réseau, ce réseau étant parfois abusivement qualifié de matriciel alors que les positions et orientations relatives des magnétomètres ne sont pas nécessairement ordonnées de manière régulière. Afin d'estimer des paramètres de l'objet magnétique tels que sa position et son orientation dans un repère associé au réseau de magnétomètres, typiquement un repère solidaire du support portant les magnétomètres, on peut mettre en oeuvre un procédé de localisation qui exploite : - un modèle permettant de simuler le champ magnétique rayonné par l'objet en tout point de l'espace ; - des mesures issues des magnétomètres du réseau ; et - un algorithme capable d'estimer la position et l'orientation de l'objet magnétique en exploitant l'écart entre les mesures attendues et les mesures observées. La mise en oeuvre d'un tel procédé suppose de connaitre au préalable les positions et les orientations des magnétomètres dans le repère du réseau. Or cette connaissance peut être erronée du fait des tolérances de fabrication et/ou d'assemblage 2 du réseau, ainsi que du fait des tolérances de fabrication des magnétomètres (la position et l'orientation réelles des axes de mesures d'un magnétomètre ne sont généralement pas fournies, ou quand elles le sont, c'est au mieux avec une tolérance de l'ordre de 500um et 1° respectivement). Il apparaît donc primordial de procéder à une calibration 5 des magnétomètres dans le repère du réseau (donc post-fabrication et post-assemblage) avant de pouvoir exploiter les mesures pour disposer d'une localisation suffisamment précise de l'objet magnétique. A titre d'exemple de technique de calibration connue d'un réseau de magnétomètres, on peut citer le document FR 2 915 568 B1 qui divulgue un procédé qui 10 permet de calibrer les sensibilités et les défauts d'orthogonalité des axes de mesure d'un ou plusieurs magnétomètres en effectuant des rotations dans un champ magnétique uniforme. La démarche peut être étendue à l'ensemble des magnétomètres d'un réseau à condition que le champ magnétique soit homogène et uniforme sur l'ensemble du volume occupé par le réseau au cours des rotations. Cette technique permet de remonter 15 aux orientations des magnétomètres du réseau. En revanche, elle ne permet pas d'accéder à la position des magnétomètres au sein du réseau. Un autre exemple de technique de calibration connue d'un réseau de magnétomètres est fourni par le document US 7,133,793 B2 qui propose de calibrer une dalle de magnétomètres (digitaliseur X-Y) destinée à localiser la position XY d'une source 20 magnétique de type stylet. Le système de calibration est composé d'une grille de générateurs de champ (bobines) positionnée au-dessus de la dalle. L'alignement de la grille avec la dalle est mesuré par imagerie. Les générateurs, dont la position est connue au sein de la grille, sont activés séquentiellement. Un algorithme compare les positions connues des générateurs avec celles estimées par la dalle, de sorte à générer une carte de 25 calibration. Cette carte est ensuite utilisée pour corriger les décalages de positionnement X-Y et les non linéarités de la dalle. Dans ce document US 7,133,793 B2, la grille est utilisée comme référence et doit donc être parfaitement réalisée et/ou caractérisée. Sa maille et sa surface doivent en outre être adaptée à la surface de la dalle et à la résolution spatiale souhaitée. La calibration nécessite par ailleurs un système d'imagerie de positionnement 3 ainsi qu'une électronique de pilotage des bobines. Cette technique de calibration apparaît donc complexe et onéreuse, pour ne traiter en outre qu'une localisation 2D. EXPOSÉ DE L'INVENTION L'invention a pour objectif d'améliorer la précision d'une localisation 5 (position et/ou orientation) d'un objet mobile par un réseau de magnétomètres en proposant une technique simple, rapide et bas coût pour calibrer collectivement les magnétomètres du réseau avec un minimum d'équipements externes. A cet effet, l'invention propose un procédé de calibration d'un réseau de magnétomètres, comportant : 10 - une acquisition de mesures réalisées par les magnétomètres du réseau au cours d'un déplacement en survol du réseau d'un outil portant au moins deux sources de champ magnétique ; - une reconstruction, à partir des mesures acquises et pour une localisation possible des magnétomètres du réseau, d'une trajectoire des sources de champ 15 magnétique au cours du déplacement ; - un calcul d'une estimation d'un attribut de l'outil à partir de la trajectoire reconstruite; et - une comparaison entre l'attribut de l'outil et l'estimation calculée dudit attribut. 20 Certains aspects préférés mais non limitatifs de ce procédé sont les suivants : - les étapes de reconstruction, de calcul et de comparaison sont réitérées pour une ou plusieurs autres localisations possibles des magnétomètres jusqu'à ce que le résultat de l'étape de comparaison soit inférieur à un seuil ; 25 - l'étape de reconstruction comprend, pour chaque instant d'échantillonnage des mesures acquises, l'évaluation d'une ou plusieurs caractérisations possibles des sources de champ magnétique, l'évaluation d'une caractérisation possible comprenant le calcul d'un écart entre les mesures acquises et une estimation du 4 champ magnétique généré à ladite localisation possible des magnétomètres par les sources de champ magnétique présentant ladite caractérisation possible ; - une caractérisation possible des sources de champ magnétique comprend des caractéristiques géométriques et magnétiques constituées des positions et des 5 moments magnétiques des sources de champ magnétique ; - l'attribut de l'outil est une distance entre sources de champ magnétique portées par l'outil ; - l'attribut de l'outil est une amplitude du moment magnétique de chacune des sources de champ magnétique ; 10 - l'attribut de l'outil est une orientation relative des moments magnétiques de chacune des sources de champ magnétique ; - l'amplitude du moment magnétique de chacune des sources de champ magnétique est identique ; - l'outil comporte deux sources de champ magnétique dont les vecteurs directeurs 15 de moment magnétique sont orientés selon des directions opposées ; - le déplacement de l'outil comprend une pluralité de phases statiques ; - le déplacement de l'outil est réalisé à main levée par un opérateur ; - le réseau de magnétomètres est constitué de magnétomètres appartenant à une pluralité de structures de magnétomètres associées ensemble. 20 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 dessins annexés sur lesquels : 25 - les figures 1 et 2 sont des schémas représentant deux variantes d'un outil porte- aimants pouvant être utilisé dans le cadre de l'invention ; - la figure 3 est un schéma illustrant la définition d'un repère associé au réseau de magnétomètres ; 5 - la figure 4 représente un repère associé au réseau et l'outil porte-aimants dont les coordonnées sont exprimées dans le repère associé au réseau ; - la figure 5 est un ordinogramme illustrant l'enchaînement des étapes du procédé selon l'invention. 5 EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS L'invention concerne un procédé de calibration collective de magnétomètres portés par un support mécanique et formant un réseau de magnétomètres. L'objectif de la calibration est de déterminer une localisation, c'est-à-dire une position et/ou une orientation, de chaque magnétomètre formant le réseau, 10 dans un repère associé au réseau, par exemple un repère ayant pour origine l'un des magnétomètres. Le procédé de calibration selon l'invention exploite un outil portant au moins deux sources de champ magnétique, telles que des aimants ou des électroaimants dont un attribut (c'est-à-dire une caractéristique ou encore une propriété) est connu, par 15 exemple un attribut géométrique ou magnétique. On prend dans ce qui suit l'exemple d'aimants et l'outil est dit porte-aimants. Comme représenté sur les figures 1 et 2, l'outil porte-aimants 10 peut prendre la forme d'un barreau en un matériau amagnétique, par exemple de l'aluminium, sur lequel sont fixés deux aimants 1, 2, avec typiquement un aimant agencé à chacune des extrémités du barreau. Bien entendu, l'outil porte-aimants 20 10 pourrait prendre toute autre forme, tel qu'un disque. Un exemple d'attribut géométrique est celui de la distance relative D entre le barycentre des aimants 1, 2. On suppose que chacun de ces aimants 1, 2 est assimilable à un dipôle magnétique, et est donc ponctuel. Pour connaitre la distance relative D, il faut donc dans la pratique connaitre les dimensions des aimants et mesurer 25 la distance entre les aimants, de manière à pouvoir en déduire la distance entre les deux barycentres. Un exemple d'attribut magnétique est celui de l'amplitude des moments magnétiques de chacun des aimants. Un autre exemple d'attribut magnétique est celui de l'orientation relative des moments magnétiques des aimants 1, 2. 6 On prévoit de préférence que les amplitudes ml, m2 des moments magnétiques M1 = m1u1, M2 = m2u2 des aimants 1, 2 sont identiques afin de limiter les effets de masquage magnétique de l'un par rapport à l'autre. Par ailleurs, les vecteurs directeurs u1, u2 des moments magnétiques M1, M2 des aimants 1, 2 sont de préférence 5 orientés selon des directions opposées de manière à accentuer les gradients magnétiques au niveau des magnétomètres du réseau et ainsi enrichir l'information. Ces vecteurs ou orthogonaux (figure 2) à l'axe directeurs sont par exemple parallèles (figure 1) longitudinal du barreau. Le ou les attributs de l'outil porte-aimants 10, par exemple les 10 amplitudes des moments magnétiques, l'orientation relative des moments magnétiques et la distance relative D, sont connus avec un niveau de précision cohérent avec celui de la calibration attendue. Pour une calibration fine, l'amplitude des moments magnétiques est connue typiquement à mieux que 0.5% près, leur orientation relative à mieux que 1° et la distance D à mieux que 100u.m. Des précisions respectivement supérieures à 10%, 15 10° et 1mm sont jugées non suffisantes pour aboutir à une calibration qualitativement significative. L'amplitude ml, m2 du moment magnétique des aimants 1, 2 qui dépend de la nature du matériau et du volume des aimants, est adaptée à la résolution des magnétomètres du réseau et à la portée visée, par exemple 0.3A.m2 pour un 20 déplacement de l'outil 10 à quelques centimètres au-dessus d'un réseau comportant des magnétomètres sensibles à quelques centaines de nano Tesla (0.3A.m2 correspond à environ 300mm3 de matériau magnétique Néodyme-Fer-Bore). La distance D entre les aimants est suffisante pour que ces deux sources magnétiques soient séparables du point de vue des magnétomètres. Cette distance peut 25 avantageusement être réduite en positionnant les aimants tête-bêche afin de générer des gradients magnétiques locaux forts. En toute rigueur, cette distance D ne doit pas être inférieure à la plus petite distance séparant deux magnétomètres du réseau. Cette distance D est telle qu'une projection des aimants 1, 2 sur le réseau, par exemple sur un plan horizontal si le réseau est à l'horizontal, est inscrite à l'intérieur de l'aire occupée par 30 l'ensemble des magnétomètres formant le réseau. Elle peut typiquement être fixée à 7 10cm pour un survol à quelques centimètres de magnétomètres espacés de quelques centimètres et occupant la surface d'une feuille A4. Une distance de survol trop faible et/ou un moment magnétique trop fort au regard de la dynamique de mesure des magnétomètres conduisent à un champ 5 magnétique trop élevé et donc à la saturation des magnétomètres. A l'inverse un survol à une distance trop élevée et/ou un moment magnétique trop faible conduisent à un champ magnétique trop faible au regard du bruit de mesure des magnétomètres. Ainsi, on comprend que la distance de survol est choisie en fonction des performances de mesure des magnétomètres et du moment magnétique des aimants. En pratique, la 10 valeur du champ au niveau des magnétomètres les plus proches des aimants doit être typiquement 300 à 3000 fois supérieure à leur bruit de mesure. En référence à la figure 5, le procédé selon l'invention se décompose essentiellement en deux phases : une première phase ENR d'acquisition de mesures réalisées par les magnétomètres du réseau au cours d'un déplacement de l'outil 10 en 15 survol du réseau de magnétomètres, et une deuxième phase RES au cours de laquelle les mesures acquises sont traitées. Ce procédé se résume comme suit. L'outil porte-aimants 10 est déplacé au-dessus du réseau de magnétomètres, ceux-ci enregistrant l'évolution de la distribution du champ magnétique dont notamment celle générée par l'outil porte-aimants 10. Parallèlement, un modèle 20 permet de décrire le comportement magnétique des aimants et de prédire la distribution spatiale du champ magnétique générée par le porte-aimants 10. En effet, en connaissant les caractéristiques géométriques et magnétiques des aimants, un modèle de prédiction peut être élaboré, comme cela sera décrit plus loin. En comparant les mesures effectives à des mesures prédites, la position et l'orientation des magnétomètres dans le repère du 25 réseau peuvent être estimées. On prend ici l'exemple d'un outil portant deux sources de champ magnétique 1, 2. L'invention n'est cependant pas limitée à cet exemple et s'étend à un outil portant plus que deux sources de champ magnétique, par exemple trois aimants agencés pour former un dièdre, voire une matrice d'aimants, ce qui peut s'avérer 30 avantageux car les signaux enregistrés par les magnétomètres sont alors plus riches, avec 8 un rapport signal à bruit logiquement plus favorable, ce qui permet de réduire le nombre de mesures et donc la durée de l'enregistrement. On dispose alors de plusieurs distances et de plusieurs angles (distances deux à deux et angles deux à deux). Par exemple, on peut disposer les aimants aux sommets de triangles, de telle manière que l'angle entre les 5 moments magnétiques est de 27r/3. Phase d'acquisition ENR L'outil 10 est déplacé au-dessus des magnétomètres qui chacun enregistre pendant ce temps l'évolution temporelle du champ magnétique perçu. Ce déplacement peut être effectué à main levée par un opérateur, de manière aléatoire. 10 Toute la surface du réseau peut ainsi être balayée avec une hauteur de vol et une orientation de l'outil potentiellement variables. Les signaux sont d'autant plus riches en information que la trajectoire est variée. Les hauteurs de vol respectives des deux aimants sont de préférence maintenues proches l'une de l'autre de manière à ce qu'un aimant ne prenne pas le pas sur l'autre, les aimants générant un rapport signal à bruit 15 similaire au niveau des magnétomètres. Autrement dit, l'outil 10 est déplacé de manière horizontale et, de préférence, en recouvrant successivement toutes les parties formant le réseau de magnétomètres dans un plan horizontal parallèle au réseau de magnétomètres. La vitesse du déplacement doit être en adéquation avec la fréquence 20 d'échantillonnage des magnétomètres, typiquement de l'ordre de 1cm/s pour une fréquence d'échantillonnage de 100Hz. Au total, le déplacement peut ainsi durer plusieurs dizaines de secondes. Le déplacement peut aussi comprendre une pluralité de phases statiques, l'outil restant fixe pendant plusieurs secondes pour différentes positions au- 25 dessus de la matrice. L'enregistrement est alors moins riche en termes d'information mais les phases statiques permettent de mieux filtrer le bruit de mesure. Pour ce faire, un support apte à porter l'outil peut être posé sur le réseau en différentes positions successives, prédéfinies ou non. 30 9 Phase de traitement RES Cette phase RES comprend la résolution d'un problème d'optimisation pour déterminer une localisation des magnétomètres du réseau. Plus particulièrement, ladite localisation vient minimiser, au sens d'un critère d'optimisation, l'écart entre un 5 attribut de l'outil 10 et une estimation de cet attribut calculée, sans connaissance du déplacement, à partir de ladite localisation et des mesures acquises lors de la phase d'acquisition ENR. Par localisation venant minimiser l'écart au sens d'un critère d'optimisation, on entend typiquement une localisation associée à un écart inférieur à un 10 seuil prédéterminé ou une localisation associée à un écart minimal après un nombre maximum d'itérations autorisées dépassé. Dans un mode de réalisation possible, la localisation déterminée par la résolution du problème d'optimisation comprend simultanément les positions et orientations des magnétomètres du réseau. 15 Dans un autre mode de réalisation possible, le procédé comprend une étape préalable de détermination des orientations de chacun des magnétomètres du réseau, par exemple mise en oeuvre selon la technique divulguée dans le document FR 2 9515 568 B1 discuté précédemment. La localisation déterminée par la résolution du problème d'optimisation comprend alors les positions des magnétomètres du réseau. Ce 20 mode de réalisation permet de réduire le nombre de degrés de liberté du problème d'optimisation et ainsi d'éviter de conduire à une calibration erronée, notamment lorsque le rapport signal à bruit des mesures n'est pas favorable et/ou lorsque les mesures ne sont pas suffisamment fiables (par exemple du fait d'un biais variable ou de non-linéarités). 25 La résolution du problème d'optimisation comprend l'évaluation d'une pluralité de localisations possibles des magnétomètres du réseau. Chaque évaluation d'une localisation possible des magnétomètres du réseau comporte : - à partir de la localisation possible et des mesures acquises, le calcul, sans connaissance du déplacement, d'une estimation de l'attribut de l'outil 10, et 30 - le calcul d'un écart entre l'attribut et l'estimation calculée dudit attribut. 10 Comme discuté précédemment, l'attribut connu de l'outil porte-aimants peut comprendre une distance D entre aimants 1, 2 de l'outil porte-aimants 10. En complément ou indépendamment, l'attribut connu peut comprendre une amplitude du moment magnétique de chacun des aimants de l'outil porte-aimants, ou encore leur 5 orientation relative. On notera que l'estimation de l'attribut peut permettre d'identifier les aimants et donc, par exemple, de vérifier qu'il s'agit du bon outil 10. Selon une approximation, chacun des aimants 1, 2 de l'outil porte-aimants 10 peut être assimilé à un dipôle magnétique ponctuel. Cette approximation devient valable si la distance entre un magnétomètre et l'aimant est au moins supérieure 10 à deux fois la plus grande dimension de l'aimant (par exemple sa longueur pour un aimant cylindrique allongé, ou son diamètre pour un aimant en forme de pièce de monnaie). Les lois magnétostatiques permettent alors d'exprimer le champ magnétique en tout point de l'espace selon la formule : m10-1 B (P , u, m) en le = d5 (3 P - d2u) Avec x 15 - P = y la position de l'aimant par rapport au magnétomètre ; z - c/ = Pli la distance entre l'aimant et le magnétomètre, telle que d2 = x2 + y2 + Z2 ; - u est le vecteur directeur du moment magnétique de l'aimant ; et - m est l'amplitude du moment magnétique de l'aimant. 20 On notera que des modèles magnétiques plus détaillés, comme par exemple le modèle ampérien décrit dans le document WO 2017/005915 Al, peuvent aussi être utilisés pour modéliser la distribution spatiale du champ magnétique généré par les deux aimants 1, 2, en particulier dans le cas où l'aimant passe très près des magnétomètres, l'approximation dipolaire n'étant alors plus valable. 25 Les coordonnées x, y, z sont ici exprimées dans le repère associé au réseau. Comme représenté sur la figure 3, ce repère peut être obtenu en choisissant, par exemple arbitrairement, un premier capteur Cl comme origine du repère, un deuxième 11 capteur C2 pour définir un premier axe de référence, et un troisième capteur C3 pour définir un deuxième axe de référence, de sorte que ces deux axes de référence soient orthogonaux. Ce troisième capteur C3 permet de définir le sens du second axe mais n'est pas nécessairement aligné sur celui-ci. Un troisième axe peut être défini à l'aide d'un 5 quatrième capteur, afin d'obtenir un repère orthogonal en trois dimensions pour un réseau de magnétomètres non plan. La figure 4 représente l'outil porte-aimants 10 dont les coordonnées sont exprimées dans le repère associé au réseau, ayant pour origine le capteur Cl. Pour un outil comportant deux aimants, le champ vu par chacun des magnétomètres est la 10 somme des champs générés par chacun des deux aimants. Le champ vu par le capteur Cl se calcule comme suit : m1.10-1 ,45 il 1 rn2.10-1 di 2u2), C1 = P1 ue11"1) (3 P2 - "12'12), d5 die 12 avec dii = 11/3111 et d12 = 11P211. Le champ vu par le capteur Ci s'exprime comme suit : m1.10-1 é 13c., = 5 (Ci - P1) - c/Pui) dii m2. 10-1 (3 (C, - P2) - eu2) dit 15 avec dit = II CiP111 et die = II CiP211. La distance D entre les deux aimants se calcule comme D = 11P1 - P211. Le calcul, à partir des mesures acquises lors du déplacement de l'outil porte-aimants 10, de l'estimation de l'attribut de l'outil dépend de la connaissance des positions et orientations des divers magnétomètres du réseau. Ainsi une erreur portant 20 sur cette connaissance des positions et orientations des magnétomètres se traduit par une erreur sur l'estimation de l'attribut et donc un écart entre la valeur réelle de l'attribut et son estimation. Partant d'une localisation possible des magnétomètres, il est ainsi possible de déterminer une erreur d'estimation de l'attribut, par exemple une erreur 25 d'estimation de la distance entre les aimants ou encore une erreur d'estimation de l'amplitude du moment magnétique de chaque aimant. Ainsi, l'invention permet de 12 s'affranchir d'un référentiel extérieur au réseau pour estimer des erreurs de positionnement des magnétomètres du réseau et autorise une auto-calibration facile à mettre en oeuvre à l'aide d'un outil porte-aimants 10. Pour chaque instant d'échantillonnage, on peut ainsi calculer un vecteur 5 d'erreur entre les attributs connus et l'estimation de ces attributs. On note a le vecteur des attributs connus et â le vecteur des attributs estimés pour la localisation possible des magnétomètres. Une erreur quadratique peut être calculée sur l'ensemble des n instants d'échantillonnage selon l'équation n°1 suivante : E =E7,2_1(â(loc,[Bo..'1],) - a)2, avec loc le vecteur contenant les paramètres de localisation que l'on cherche à trouver 10 (positions et, le cas échéant, orientations des magnétomètres ; qui restent fixes au cours du temps), et [Bo.,_1], le vecteur des mesures magnétiques réalisées par les magnétomètres 0 à m-1 à l'instant d'échantillonnage i. Ainsi, pour trouver le vecteur loc, on cherche à résoudre le problème d'optimisation consistant à minimiser cette erreur quadratique selon l'équation n°2 15 suivante. loc = argminzoc E7,',,(e/f/oc, [B0 - a)2. Différents algorithmes d'optimisation peuvent être utilisés pour minimiser cette fonction. Par exemple, il est possible d'utiliser une descente de gradient sous réserve de dérivabilité, c'est à dire de pourvoir calculer une variation de â en 20 fonction d'une variation de loc. Dans ce cas, la dérivée est calculée numériquement. Il est aussi possible d'utiliser un algorithme génétique qui permet de s'affranchir des problèmes liés au calcul de variation mentionné précédemment. Un algorithme hybride, utilisant des mutations aléatoires lorsqu'un problème de dérivabilité survient et une descente de gradient classique en l'absence de problèmes de dérivabilité, peut être 25 préféré. La résolution de l'équation n°2 peut également être réalisée au moyen d'un filtre de Kalman, d'un filtre particulaire ou d'un réseau de neurones. Le calcul du vecteur â des attributs estimés pour une localisation possible des magnétomètres peut être effectué en plusieurs étapes. Par exemple, si â porte sur la distance D entre les aimants, une première étape consiste à estimer la 30 trajectoire suivie par les aimants lors du déplacement (les positions et moments 13 magnétiques des deux aimants à chaque instant i) selon une technique similaire d'optimisation : on cherche les positions et moments magnétiques des aimants qui minimisent l'erreur, par exemple quadratique, entre les valeurs de champs magnétiques mesurées à l'instant i et celles estimées en fixant les positions et orientations des 5 magnétomètres. Une fois les positions et moments magnétiques des aimants estimés sur l'ensemble des instants i, on calcule la distance D grâce à la formule donnée précédemment. Si le vecteur â porte sur les normes des deux moments magnétiques, l'étape mentionnée ci-dessus est effectuée, puis les normes sont déduites des valeurs des moments (vecteurs 3D) estimés. Si le vecteur â porte sur l'orientation relative des deux 10 aimants, l'étape ci-dessus est effectuée, puis l'orientation relative est estimée par exemple à partir du produit scalaire entre les 2 vecteurs directeurs u1 et u2. Ainsi, l'estimation de l'attribut de l'outil porte-aimants (i.e., le calcul du vecteur â pour une localisation possible des magnétomètres) est précédée d'une étape de reconstruction, à partir des mesures acquises et pour une localisation possible des 15 magnétomètres du réseau, d'une trajectoire des aimants de l'outil porte-aimants au cours du déplacement. Cette étape de reconstruction comprend, pour chaque instant d'échantillonnage des mesures acquises, l'évaluation d'une ou plusieurs caractérisations possibles des aimants de l'outil porte-aimants. Par caractérisation des aimants, on entend 20 notamment les caractéristiques géométriques et magnétiques des aimants constituées par leurs positions et leurs moments magnétiques. L'évaluation d'une caractérisation possible comprend le calcul d'un écart entre les mesures acquises et une estimation du champ magnétique généré à ladite localisation possible des magnétomètres par les aimants présentant ladite caractérisation 25 possible. Puis, pour le calcul de l'estimation de l'attribut, on retient pour chacun des instants d'échantillonnage, la caractérisation possible associée à l'écart le plus faible. Au global, il y a donc une optimisation récursive, à double étage : la première vise à estimer la trajectoire des aimants lors du déplacement de l'outil (les positions et moments magnétiques des aimants) pour une localisation des 30 magnétomètres donnée (pour le calcul de â et de l'erreur e selon l'équation n°1), la 14 seconde vise à converger vers les positions et/ou orientations des magnétomètres qui minimisent l'erreur selon l'équation n°2. La première optimisation peut être réalisée selon l'un quelconque des algorithmes, filtres ou réseaux mentionnés précédemment en lien avec la seconde optimisation, sans nécessairement que les moyens mis en oeuvre 5 pour réaliser ces deux optimisations soient identiques. A titre d'exemple, un filtre de Kalman peut être mis en oeuvre pour la première optimisation tandis qu'un algorithme de descente de gradient peut être mis en oeuvre pour la seconde optimisation. Ainsi, la phase de traitement RES comprend : - une reconstruction, à partir des mesures acquises et pour une localisation 10 possible des magnétomètres du réseau, d'une trajectoire des aimants de l'outil porte-aimants au cours du déplacement ; - un calcul d'une estimation d'un attribut de l'outil porte-aimants à partir de la trajectoire reconstruite; et - une comparaison entre l'attribut de l'outil porte-aimants et l'estimation 15 calculée dudit attribut. Les étapes de reconstruction, de calcul et de comparaison sont réitérées pour une ou plusieurs autres localisations possibles des magnétomètres jusqu'à ce que le résultat de l'étape de comparaison soit inférieur à un seuil. La localisation possible associée à cet égard minimal constitue alors la localisation optimale. 20 Dans tous les cas de figure, il est préférable d'initialiser l'algorithme avec le maximum de connaissances a priori, comme les positions et orientations relatives des magnétomètres (généralement, on sait par conception comment ils sont disposés dans le réseau), ainsi que leur sensibilité (données constructeur). L'algorithme est ainsi initialisé avec des données réalistes qu'il va affiner, ce qui facilite sa convergence. 25 La calibration réalisée par l'invention est simple et rapide. Elle ne nécessite ni d'équipements ni d'instruments externes lourds et couteux. Une calibration complète des magnétomètres (position, orientation et sensibilité) en une seule opération peut potentiellement être réalisée. L'invention peut par ailleurs être employée quelle que soit la distribution des magnétomètres, à condition que cette distribution reste fixe 15 pendant et après la calibration (une connaissance a priori de cette distribution étant par ailleurs souhaitable pour faciliter la convergence de l'algorithme de localisation). L'invention peut donc aussi être utilisée pour les systèmes à plusieurs structures de magnétomètres en venant exprimer les coordonnées de l'ensemble des 5 magnétomètres dans un repère associé à l'une des structures prise comme référence. De tels systèmes comprennent par exemple des structures accolées en 2D ou 3D, de manière contiguë ou non, pour étendre le volume de mesure, densifier ce volume de capture ou changer la géométrie du volume de capture, comme par exemple des structures planes disposés en V. Dans un tel cas de figure, le réseau de magnétomètres au sens de 10 l'invention est constitué de l'ensemble des magnétomètres des différentes structures associées ensemble. L'invention n'est pas limitée au procédé tel que précédemment décrit, mais s'étend également à un dispositif de magnétométrie comprenant un réseau de magnétomètres et un calculateur configuré pour mettre en oeuvre ce procédé ainsi qu'à 15 un programme d'ordinateur comprenant des instructions qui, lorsque le programme est exécuté par un ordinateur, REVENDICATIONS1. Procédé de calibration d'un réseau de magnétomètres (C1, C2, C3, Ci), comportant : - une acquisition (ENR) de mesures réalisées par les magnétomètres du réseau au cours d'un déplacement en survol du réseau d'un outil portant au moins deux sources de champ magnétiques ; - une reconstruction, à partir des mesures acquises et pour une localisation possible des magnétomètres du réseau, d'une trajectoire des sources de champ magnétique au cours du déplacement ; - un calcul d'une estimation d'un attribut (D, m1, m2) de l'outil à partir de la trajectoire reconstruite; et - une comparaison entre l'attribut de l'outil et l'estimation calculée dudit attribut. 15 2. Procédé de calibration selon la revendication 1, dans lequel les étapes de reconstruction, de calcul et de comparaison sont réitérées pour une ou plusieurs autres localisations possibles des magnétomètres jusqu'à ce que le résultat de l'étape de comparaison soit inférieur à un seuil. 20 3. Procédé de calibration selon l'une des revendications 1 et 2, dans lequel l'étape de reconstruction comprend, pour chaque instant d'échantillonnage des mesures acquises, l'évaluation d'une ou plusieurs caractérisations possibles des sources de champ magnétique de l'outil, l'évaluation d'une caractérisation possible 25 comprenant le calcul d'un écart entre les mesures acquises et une estimation du champ magnétique généré à ladite localisation possible des magnétomètres par les sources de champ magnétique présentant ladite caractérisation possible. 17 4. Procédé selon la revendication 3, dans lequel une caractérisation possible des sources de champ magnétique comprend des positions et des moments magnétiques des sources de champ magnétique. 5 5. Procédé de calibration selon l'une des revendications 1 à 4, dans lequel l'attribut de l'outil est une distance (D) entre sources de champ magnétique (1, 2) de l'outil (10). 6. Procédé de calibration selon l'une des revendications 1 à 4, dans lequel 10 l'attribut de l'outil est une amplitude (m1, m2) du moment magnétique (M1, M2) de chacune des sources de champ magnétique (1, 2) de l'outil (10). 7. Procédé de calibration selon l'une des revendications 1 à 4, dans lequel l'attribut de l'outil est une orientation relative des moments magnétiques (M1, M2) de chacune des sources de champ magnétique (1, 2) de l'outil (10). 8. Procédé de calibration selon l'une des revendications 1 à 7, dans lequel l'amplitude (m1, m2) du moment magnétique de chacune des sources de champ magnétique (1, 2) de l'outil porte-aimants (10) est identique. 9. Procédé de calibration selon l'une des revendications 1 à 8, dans lequel l'outil (10) comporte deux sources de champ magnétique (1, 2) dont les vecteurs directeurs de moment magnétique (u1, u2) sont orientés selon des directions opposées. 10. Procédé selon l'une des revendications 1 à 9, dans lequel le déplacement de l'outil comprend une pluralité de phases statiques. 11. Procédé selon l'une des revendications 1 à 10, dans lequel le déplacement de l'outil est réalisé à main levée par un opérateur. 18 12. Procédé selon l'une des revendications 1 à 11, dans lequel le réseau de magnétomètres est constitué de magnétomètres appartenant à une pluralité de structures de magnétomètres associées ensemble. 5 13. Dispositif de magnétométrie comprenant un réseau de magnétomètres et un calculateur configuré pour mettre en oeuvre le procédé de calibration selon l'une des revendications 1 à 12. 14. Programme d'ordinateur comprenant des instructions qui, lorsque le 10 programme est exécuté par un ordinateur, conduisent celui-ci à mettre en oeuvre le procédé selon l'une des revendications 1 à 12.