FR 3069650 A1 20190201 FR 1757084 A 20170726 Procédé de calibration d'un magnétomètre DOMAINE TECHNIQUE GENERAL La présente invention concerne le domaine de la calibration de capteurs. Plus précisément, elle concerne un procédé de calibration d'au moins un magnétomètre. Un magnétomètre peut être utilisé dans pour diverses applications, 10 comme par exemple la navigation sans GPS, ou la détection des métaux. ETAT DE L'ART Le champ magnétique est un champ vectoriel en l'espèce 15 tridimensionnel, c'est-à-dire associant un vecteur champ magnétique noté É et de dimension trois à chaque emplacement de l'espace. Le champ magnétique en un point se caractérise ainsi par sa norme et sa direction. Un magnétomètre permet de mesurer le champ magnétique en un point. Un magnétomètre peut être un magnétomètre mono-axial, capable de 20 mesurer une composante du champ magnétique dans une position, i.e. la projection du vecteur champ magnétique É au niveau dudit magnétomètre selon son axe. Un magnétomètre peut également être un magnétomètre triaxe, capable de mesurer le champ magnétique selon ses trois axes, qui peut être constitué de trois magnétomètres mono-axe rigidement liés entre 25 eux et orientés selon des axes différents, généralement sensiblement orthogonaux. Toutefois, la mesure du champ magnétique par un magnétomètre n'est pas parfaite, de sorte que la mesure du champ magnétique dévie du champ magnétique réel. Par exemple, des effets de fers durs peuvent être crées par 30 rémanence des certains composants à de forts champs magnétiques et entraîner des biais de mesure. Des facteurs d'échelle peuvent également distordre la mesure. Par ailleurs, des effets dus à la configuration physique 2 du magnétomètre et au montage du magnétomètre peuvent également affecter la mesure. On peut ainsi par exemple écrire la mesure d'un magnétomètre tri-axe 5 sous la forme linéaire suivante : où Bin est la mesure du champ magnétique par le magnétomètre, É est le champ magnétique réel, D est un facteur d'échelle (représenté par une matrice diagonale ou pleine), et É est un biais de mesure. Le facteur d'échelle et le biais de mesure sont donc des paramètres de mesure qui viennent polluer la mesure du champ magnétique par rapport au champ magnétique réel. Il convient donc d'estimer ces paramètres, ou au moins un, afin de pouvoir corriger la mesure du champ magnétique. La calibration d'un magnétomètre revient à déterminer au moins un de 15 ces paramètres et à en déduire des corrections à apporter aux sorties du magnétomètre afin que la mesure corrigée reflète plus fidèlement le champ magnétique réel. De plus, la calibration permet de déterminer la positon du magnétomètre dans un repère lié au support sur lequel ledit magnétomètre 20 est (l'orientation du magnétomètre et son emplacement dans l'espace). Une calibration précise du magnétomètre est notamment requise pour les applications qui requièrent une grande précision. En particulier, une centrale magnéto-inertielle peut comprendre, outre des accéléromètres, un réseau de plusieurs magnétomètres disposés à des emplacements différents 25 qui permettent l'estimation simultanée du champ et du gradient magnétique. La précision de ces mesures est essentielle à la bonne performance de reconstruction du mouvement de la centrale magnéto-inertielle. Une méthode de calibration existante, appelée « méthode en champ libre », consiste à placer les magnétomètres à calibrer dans un endroit aussi 30 éloigné que possible de toute source de perturbations magnétiques. Cette méthode est très contraignante car les endroits éloignés des sources de perturbations magnétiques sont très reculés. De plus, les conditions 3 d'isolement de toute source de perturbations magnétiques sont de plus en plus compromises par les diverses constructions humaines. Une autre méthode de calibration connue utilise des bobines de type Helmholtz, Braunbek, ou autre. Cette solution nécessite toutefois de procéder 5 au préalable à un étalonnage des bobines avec un magnétomètre étalon. Une grande part de la précision du calibrage des magnétomètres dépend de cette étape d'étalonnage des bobines qui d'une part est complexe et d'autre part nécessite une grande précision. De plus, afin d'obtenir la précision désirée, il peut alors être nécessaire d'étalonner le magnétomètre étalon, complexifiant et rallongeant d'autant plus la méthode de calibration. En outre, certaines méthodes de calibration connues utilisant des bobines de Helmholtz nécessitent de placer les magnétomètres à calibrer dans une zone sur laquelle un champ magnétique homogène est créé par les bobines de Helmholtz. Une telle méthode est par exemple décrite dans le document US2017/0067981. Une telle méthode possède l'inconvénient qu'il devient nécessaire d'utiliser des bobines avec un rayon très important pour pouvoir calibrer certains objets, ce qui pose des problèmes d'encombrement. En effet, les magnétomètres à calibrer doivent être situés dans la zone sur laquelle le champ magnétique créé par les bobines de Helmholtz est homogène, et la taille de cette zone sur laquelle le champ magnétique est homogène dépend de la taille des bobines de Helmholtz. Or, les magnétomètres sont calibrés en étant montés sur un support, comme une carte électronique par exemple, qui est généralement disposé dans un objet, comme un boîtier par exemple, car l'objet et le support peuvent interférer sur la mesure du champ magnétique par les magnétomètres. La zone sur laquelle le champ magnétique est homogène doit être suffisamment grande afin de pouvoir contenir l'objet. Il serait donc souhaitable de disposer d'une nouvelle méthode de calibration d'un magnétomètre qui est à la fois précise et simple à mettre en oeuvre, et qui ne nécessite pas d'être effectuée dans un endroit éloigné de toute source de perturbations magnétiques. PRESENTATION DE L'INVENTION La présente invention se rapporte ainsi selon un premier aspect à un procédé de calibration d'un magnétomètre, le procédé étant caractérisé en 5 ce qu'il comprend les étapes suivantes : El : placer au moins un magnétomètre à l'intérieur d'un générateur d'un champ magnétique, ledit générateur comprenant une pluralité de bobinages ; E2 : générer un champ magnétique dans le générateur en appliquant un courant électrique à chaque bobinage ; 10 E3: mesurer le champ magnétique généré dans le générateur avec le au moins un magnétomètre, et mesurer les courants électriques appliqués aux bobinages du générateur ; E4 : répéter l'étape E3 en modifiant à chaque fois une position dudit au moins un magnétomètre à l'intérieur du générateur et/ou les courants 15 électriques dans les bobinages du générateur, de sorte à obtenir des mesures du champ magnétique et des courants électriques appliqués aux bobinages pour une pluralité de positions dudit au moins un magnétomètre et/ou de courants électriques dans les bobinages ; E5: générer un modèle de transfert paramétrique à partir d'un modèle 20 de mesure paramétrique dudit au moins un magnétomètre et d'un modèle paramétrique du champ magnétique, le modèle de mesure paramétrique dudit au moins un magnétomètre comprenant des paramètres de calibration dudit au moins un magnétomètre ; E6 : calculer les paramètres de calibration du au moins un 25 magnétomètre avec un algorithme d'optimisation à partir et du modèle de transfert paramétrique généré à l'étape E5 et des mesures du champ magnétique et des courants électriques réalisées dans les étapes E3 et E4. Selon d'autres caractéristiques avantageuses et non limitatives : 30 - le procédé comprend l'étape suivante : E40: calculer un modèle de transfert mesuré entre la mesure du champ magnétique par le au moins un magnétomètre et les courants électriques 4 5 appliqués aux bobinages pour les différentes positions dudit au moins un magnétomètre et/ou les différents courants électriques dans les bobinages, à partir d'une part des mesures dudit champ magnétique et d'autre part des mesures desdits courants électrique appliqués aux bobinages, 5 le calcul des paramètres de calibration étant réalisé avec un algorithme d'optimisation à partir du modèle de transfert paramétrique et des modèles de transfert mesurés ; - le procédé comprend en outre une étape de contrôle la température dudit au moins un magnétomètre à calibrer ; 10 - le au moins un magnétomètre possède une réponse linéaire, le modèle de transfert mesuré étant un matrice de transfert mesurée qui est calculée par régression linéaire, le modèle de transfert paramétrique étant une matrice de transfert paramétrique ; - le au moins un magnétomètre possède une réponse non-linéaire, le 15 modèle de transfert mesuré étant calculé avec un modèle polynomial ou un modèle en séries de Fourier ; - le générateur est un dispositif à bobines de Helmholtz triaxes comprenant six bobinages ; - le procédé comprend une étape d'identification d'un biais dudit au moins 20 un magnétomètre, l'étape d'identification étant réalisée en plaçant ledit au moins un magnétomètre dans un dispositif isolant ledit au moins un magnétomètre de tout champ magnétique ; - le procédé est appliqué sur une pluralité de magnétomètres, les paramètres de calibration des magnétomètres comprenant un emplacement 25 des magnétomètres défini par rapport à un des magnétomètres servant de référence ; - le au moins un magnétomètre est maintenu fixe durant l'étape E2 de génération du champ magnétique à l'intérieur du générateur et durant l'étape E3 de mesure du champ magnétique et du courant électrique appliqué à 30 chacun des bobinages ; - les paramètres de calibration du au moins un magnétomètre comprennent : une position dudit au moins un magnétomètre sur un support 6 sur lequel ledit au moins un magnétomètre est monté, un facteur directionnel de mesure dudit au moins un magnétomètre, et/ou un biais dudit au moins un magnétomètre ; - le modèle paramétrique du champ magnétique est un modèle de champ 5 d'ordre 1, un modèle de champ d'ordre 2, ou un modèle à éléments finis. Selon un deuxième aspect, l'invention concerne un équipement de calibration d'au moins un magnétomètre comprenant au moins un magnétomètre, un générateur d'un champ magnétique qui comprend une 10 pluralité de bobinages, caractérisé en ce qu'il comprend également un capteur de courant configuré pour mesurer le courant dans chacun des bobinages, et des moyens de traitement qui sont configurés pour : - générer un modèle de transfert paramétrique à partir d'un modèle de mesure paramétrique dudit au moins un magnétomètre et d'un modèle 15 paramétrique du champ magnétique, le modèle de mesure paramétrique dudit au moins un magnétomètre comprenant des paramètres de calibration dudit au moins un magnétomètre ; - calculer les paramètres de calibration du au moins un magnétomètre avec un algorithme d'optimisation à partir du modèle de transfert 20 paramétrique et de mesures du champ magnétique par ledit au moins un magnétomètre et de mesures des courants dans les bobinages. Selon une autre caractéristique avantageuse et non limitative, le générateur est un dispositif à bobines de Helmholtz triaxes comprenant six 25 bobinages. Selon un troisième aspect, l'invention concerne un produit programme d'ordinateur comprenant des instructions de code pour l'exécution d'un procédé de calibration d'au moins un magnétomètre selon le premier aspect 30 de l'invention. 7 Selon un quatrième aspect, l'invention concerne un moyen de stockage lisible par un équipement informatique sur lequel un produit programme d'ordinateur comprend des instructions de code pour l'exécution d'un procédé de calibration d'au moins un magnétomètre selon le premier 5 aspect de l'invention. Un but de l'invention est de simplifier la calibration d'un magnétomètre, notamment par détermination du facteur directionnel de mesure du magnétomètre, de son biais, et/ou de son emplacement sur son support. 10 Un autre but de l'invention est de calibrer de manière simultanée un magnétomètre et un modèle paramétrique du champ magnétique. Un autre but de l'invention est de pouvoir étalonner ou optimiser les systèmes de calibration actuels. 15 PRESENTATION DES FIGURES D'autres caractéristiques et avantages de la présente invention apparaîtront à la lecture de la description qui va suivre d'un mode de réalisation préférentiel. Cette description sera donnée en référence aux dessins annexés dans lesquels : 20 - la figure 1 est un schéma représentant une variante possible d'un équipement pour la mise en oeuvre du procédé selon l'invention ; - La figure 2 représente une mise en oeuvre possible du procédé selon l'invention ; - La figure 3 représente une variante possible du procédé selon 25 l'invention. DESCRIPTION DETAILLEE Principe Le présent procédé propose d'utiliser un générateur de champ 30 magnétique, comme par exemple un dispositif à bobines de Helmholtz de Braunbek ou de Maxwell. 8 De plus, le procédé propose, au lieu de calibrer le générateur puis de calibrer les magnétomètres, de calibrer simultanément les magnétomètres et le générateur. Un aspect de l'invention repose sur le fait qu'il est possible de : 5 - déterminer des paramètres de calibration des magnétomètres, notamment un emplacement et une orientation des magnétomètres sur le support sur lequel lesdits magnétomètres sont montés, un facteur directionnel de mesure, et un biais ; et - déterminer des paramètres du modèle de champ magnétique généré 10 par le générateur par utilisation d'un algorithme d'optimisation entre des mesures du champ magnétique généré par le générateur effectuées par les magnétomètres à calibrer et un modèle de mesure paramétrique, ledit modèle de mesure paramétrique étant basé sur modèle paramétrique du champ magnétique et 15 sur un modèle paramétrique de mesure du magnétomètre qui comprend les paramètres de calibration dudit magnétomètre. Il n'y a donc pas besoin de connaître le champ magnétique généré par le générateur, comme cela est le cas dans les méthodes de calibration connues dans lesquelles le générateur est calibré afin de connaitre 20 précisément les propriétés du champ magnétique généré par ledit générateur. Dans cet aspect de l'invention, il suffit d'utiliser un modèle paramétrique du champ magnétique qui est choisi par l'utilisateur sans prendre en compte les paramètres du champ magnétique que génère véritablement le générateur. 25 Les mesures du champ magnétique du générateur sont réalisées en faisant varier la position des magnétomètres à l'intérieur du générateur. Outre les paramètres de calibration du magnétomètre, cette solution peut en outre permettre de déterminer les paramètres du modèle paramétrique du champ magnétique. 30 On entend ici par position la combinaison d'un emplacement et d'une orientation, ce qui permet de décrire complètement la configuration spatiale 9 d'un objet. Dans une écriture vectorielle d'un espace à trois dimensions, l'emplacement est défini par un vecteur à trois composantes (les coordonnées spatiales dans un repère), et l'orientation est également définie par un vecteur à trois composantes (les angles de rotation par rapport au 5 repère). Equipement En référence à la figure 1, un équipement 1 de calibration d'au moins un magnétomètre 3 comprend : 10 - un générateur 2 d'un champ magnétique comprenant des bobinages 21. Le générateur 2 génère un champ magnétique lorsque des courants électriques circulent dans les bobinages 21. Dans la variante représentée sur la figure 1 le générateur 1 est un dispositif à bobines de Helmholtz triaxes comprenant six bobinages 21. Le au moins un 15 magnétomètre 3 à calibrer, de manière générale une pluralité de magnétomètres 3 disposés sur un support (comme par exemple une carte électronique), est disposé à l'intérieur du générateur 2. Dans la variante illustrée sur la figure 1, le ou les magnétomètres 3 à calibrer sont disposé(s) entre les bobinages 21 du générateur 2 à bobines de 20 Helmholtz. - des moyens de traitement 4, par exemple un ordinateur, qui comprennent un processeur 41 et une mémoire 42 et qui sont configurés pour mettre en oeuvre un procédé selon l'invention. Les moyens de traitement 4 sont connectés au générateur 2 et audit au 25 moins un magnétomètre 3. Les moyens de traitement 4 peuvent ainsi contrôler les courants électriques dans les bobinages 21 du générateur 2, et obtenir les données mesurées par le ou les magnétomètre(s) 3. - un capteur de courant 5 configuré pour mesurer le courant électrique 30 dans chacun des bobinages 21 du générateur. Le capteur de courant 5 est connecté aux moyens de traitement 4. 10 L'équipement 1 peut également comprendre des moyens de contrôle 6 de la température du au moins un magnétomètre 3. En effet, les magnétomètres peuvent être sensibles à la température. Il est ainsi 5 avantageux de maintenir le ou les magnétomètres 3 à calibrer à la température désirée avec les moyens de contrôle 6 de la température dudit au moins un magnétomètre 3. Typiquement, les moyens de contrôle 6 de la température peuvent imposer que la température à l'intérieur du générateur 2 reste constante, par exemple avec des moyens de réchauffe et/ou des 10 moyens de refroidissement, comme par exemple une résistance chauffante ou une pompe à chaleur. Procédé Comme illustré sur la figure 2, le procédé de calibration du au moins 15 un magnétomètre 3 comprend les étapes suivantes : El : placer au moins un magnétomètre 3 à calibrer à l'intérieur du générateur 2 de champ magnétique ; E2 : générer un champ magnétique dans le générateur 2 en appliquant des courants électriques aux différents bobinages 21 dudit générateur 2 ; 20 E3: mesurer le champ magnétique généré dans le générateur 2 avec le au moins un magnétomètre 3, et mesurer les courants électriques appliqués aux bobinages 21 du générateur 2 avec le capteur de courant 6; E4: répéter l'étape E3 (mesurer ledit champ magnétique ainsi que les courants électriques dans les bobines) en modifiant à chaque fois une 25 position dudit au moins un magnétomètre 3 à l'intérieur du générateur 2, et/ou les courants électriques dans les bobinages 21 du générateur 2. Une mesure du champ magnétique et des courants appliqués aux bobinages 21 étant obtenue pour chaque positionnement du magnétomètre 3 et/ou variation des courants électriques dans les bobinages 21. Typiquement, une série de 36 30 mesures peut être effectuée avec deux magnétomètres 3 monoaxe, ou 18 avec un magnétomètre 3 triaxes, en changeant de position les 11 magnétomètres 3 de position et/ou les courants électriques dans les bobinages 21 entre chaque mesure. La variation des courants électriques dans les bobinages 21 peut par exemple être fait par une variation de l'intensité et/ou une variation de la tension. Le fait de faire varier la position 5 du au moins un magnétomètre 3 et/ou des courants dans les bobinages 21 permet d'obtenir des mesures avec dans des conditions différentes afin de pouvoir être utilisées pour un algorithme d'optimisation. E5: générer un modèle de transfert paramétrique à partir d'un modèle de mesure paramétrique dudit au moins un magnétomètre 3 et d'un modèle 10 paramétrique du champ magnétique, le modèle de mesure paramétrique dudit au moins un magnétomètre comprenant les paramètres de calibration dudit au moins un magnétomètre 3. E6: calculer les paramètres de calibration du au moins un magnétomètre 3 avec un algorithme d'optimisation à partir du modèle de 15 transfert paramétrique généré à l'étape E5, et des mesures du champ magnétique et des courants pour les différents positionnements dudit au moins un magnétomètre 3 obtenues lors des étapes E3 et E4. L'algorithme d'optimisation peut être une descente de gradient, un algorithme de Levenberg-Marquardt, ou un algorithme de Metropolis-Hastings. 20 L'algorithme d'optimisation est initialisé en générant le modèle de transfert paramétrique. Une fonction coût est ensuite évaluée à partir des mesures des courants dans les bobinages 21 et des mesures du champ magnétique par le au moins un magnétomètre 3. Une valeur des paramètres de calibration du magnétomètre, ainsi que des autres paramètres du modèle 25 paramétrique, est déterminée afin de réduire la fonction coût. Selon une variante illustrée sur la figure 3, le procédé peut comprendre une étape E40 pour calculer un modèle de transfert mesuré entre le champ magnétique généré et les courants électriques appliqués aux 30 bobinages 21 pour les différentes positions dudit au moins un magnétomètre 3 et/ou les différents courants électriques dans les bobinages, à partir d'une 12 part des mesures dudit champ magnétique et d'autre part des mesures desdits courants électrique appliqués aux bobinages 21. Le calcul d'une matrice de transfert à partir des mesures peut par exemple être réalisé par régression linéaire. D'autres méthodes peuvent 5 également être utilisés pour calculer le modèle de transfert et ainsi décrire le lien entre les mesures du champ magnétique et les courants dans les bobinages 21, comme par exemple avec un modèle polynomial, ou un modèle en séries de Fourier, qui sont des méthodes adaptées si la réponse des magnétomètres 3 n'est pas linéaire. 10 Dans cette variante, le calcul des paramètres de calibration de l'étape E6 est réalisé avec un algorithme d'optimisation à partir du modèle de transfert paramétrique et des modèles de transfert mesurés. Les étapes d'évaluation d'une fonction coût et de détermination des paramètres de calibration réduisant cette fonction coût sont répétées avec les modèles de 15 transfert mesurés des différentes positions du magnétomètre dans le générateur, jusqu'à trouver un minium local ou global. Ce procédé permet donc également de déterminer les paramètres d'un modèle de champ magnétique, lesdits paramètres étant compris dans le modèle paramétrique du champ magnétique. 20 Il va maintenant être décrit comment il est possible de calculer les paramètres de calibration du ou des magnétomètre(s) 3, et les paramètres de calibration du champ magnétique généré par le générateur 2 avec un tel procédé. 25 Le champ magnétique généré par le générateur 2 est une fonction affine par rapport aux courants électriques circulant dans les bobinages 21 dudit générateur. Ainsi, si l'on note B (X) le champ magnétique en fonction de la position X à l'intérieur du générateur 2, il est possible d'écrire : 30 B (X) = T(X)1> 13 où Î est un vecteur comprenant la liste de tous les courants électriques indépendants dans chaque bobinages 21, et où T (X) désigne une matrice de transfert (un modèle de transfert dans le cas d'un magnétomètre linéaire). Les courants électriques sont mesurés par le capteur de courant 6, qui 5 transmet les données mesurées aux moyens de traitement 4. Le champ magnétique en chaque point s'obtient donc par combinaison linéaire des courants électriques dans chaque bobinage 21 par l'intermédiaire de la matrice de transfert T (X) . Un magnétomètre 3 fournit typiquement une mesure du champ 10 magnétique dans une direction donnée. Lorsqu'un magnétomètre 3 est positionné immobile à proximité des bobinages 21, le champ magnétique mesuré par ledit magnétomètre 3 peut donc être exprimé comme fonction des courants Î, et de paramètres de calibration comprenant par exemple l'emplacement X dudit magnétomètre 3 sur le support sur lequel ledit 15 magnétomètre 3 est placé, du facteur directionnel de mesure du magnétomètre (qui comprend le paramètre d'échelle dudit magnétomètre 3), et du biais dudit magnétomètre 3. Dans des bobines de Helmholtz (et dans d'autres configuration de générateur à bobinages comme ceux du type bobines de Maxwell ou 20 Braunbek) on peut montrer théoriquement que deux ensembles de mesures physiquement équivalentes ne peuvent se distinguer que par le choix : - d'une origine du repère des bobinages 21; - d'une définition des axes du repère des bobinages 21; - d'un facteur d'échelle sur le champ magnétique émis par les bobinages 25 21; - d'un facteur d'échelle sur les longueurs. Ces choix peuvent s'écrire sous la forme de contraintes sur les paramètres identifiés. Par exemple, le facteur d'échelle sur le champ magnétique émis par les bobinages 21 peut être défini en imposant le facteur 30 d'échelle d'un magnétomètre, ou en imposant une valeur de champ à l'origine du générateur 2. 14 Pour une pluralité de magnétomètres 3 monoaxiaux, à réponse linéaire (le procédé de l'invention s'applique cependant à d'autres modèles de magnétomètres, comme par exemple des magnétomètres à modèles non- 5 linéaire), qui sont fixés de manière rigide sur leur support, comme par exemple une carte électronique (appelée carte capteur) qui est montée sur un objet, lesdits magnétomètres 3 étant considérés comme ponctuels et supposés sans biais, il est possible de définir le modèle de mesure paramétrique (ou modèle de calibration) suivant pour chaque magnétomètre lo 3: Yi(X) = c/iCi(X) (1) Où Yi(X) représente la mesure scalaire du magnétomètre i lorsqu'il est placé au point X, et les Ci(X) représentent les trois composantes du champ magnétique au point X. Et où di représente un facteur directionnel de mesure 15 associé au magnétomètre i. Le facteur directionnel de mesure di est un paramètre de calibration du magnétomètre qui comprend plusieurs informations sur le magnétomètre. En effet, le facteur directionnel de mesure indique d'une part la manière dont le magnétomètre mesure les différentes composantes du champ magnétique, 20 et d'autre part le facteur d'échelle du magnétomètre. Par exemple, si : = (1 + E ; ey ; Et), avec ex, ey, ez « 1 alors le magnétomètre associé est un capteur mesurant globalement la direction X du champ magnétique et souffre d'un problème de facteur d'échelle axial par rapport à X (ex) et de problème d'alignement par rapport aux axes Y (ey) et Z (Et). L'homme de l'art 25 saura étendre ce modèle au cas d'un biais ou de non linéarités. De plus, la relation suivante existe par ailleurs : C' (X) = S "1,13k (X) (2) où S représente la transposée de la matrice d'attitude Q de la carte capteurs 30 sur laquelle sont montés les magnétomètres 3. 15 Ainsi, dans le cas où le champ magnétique est généré par un générateur 2 à bobines de Helmholtz à six bobinages 21, comme cela est le cas dans la variante illustrée sur la figure 1, il est obtenu : 5 Bk (X) = n.i.c (X, p)In (3) où n.ic (X, p.) représente le modèle paramétrique du champ magnétique, ici en l'occurrence sous la forme d'une matrice de transfert, entre le champ magnétique généré par le générateur 2 et les courants dans les six bobinages 21. Et où Pn correspond au courant électrique dans la bobine m. Le modèle paramétrique du champ magnétique n.ic (X, p.) dépend d'un ensemble de paramètres qui sont rassemblés dans le vecteur y. Le vecteur kt comprend les paramètres de calibration du champ magnétique obtenus lors de l'étape E6. 15 De plus, la relation suivante peut être établie : X[a,p]i = O[p]' Q [p] Ji P[a]1 (4) Où O[p] est la i-ème coordonnée dans le repère du générateur 2 de l'origine du repère capteurs lors de son p-ième positionnement dans le générateur 2. Le repère capteurs est un repère permettant de déterminer l'emplacement 20 des magnétomètres 3 de manière relative en utilisant un des magnétomètres 3 comme référence, formant ainsi l'origine du repère capteurs. Où Q [p] est la matrice d'attitude des magnétomètres 3, la matrice d'attitude donnant l'orientation des magnétomètres 3. Où P[a]] est la j-ème coordonnée dans le repère capteurs du magnétomètre a. Et où X[a,M est la position du 25 magnétomètre a lors de son p-ième positionnement dans le générateur 2. Ainsi, on obtient le modèle de mesure paramétrique suivant pour les magnétomètres 3: Y' (X[i,p])= cliS[p]Jkn.ic (O[p]n Q[p]ro'llir , p.)Pn (5) j k m 16 Où Y (X [i, /A) représente la mesure scalaire du magnétomètre i lors de son p-ième positionnement dans le générateur 2. La matrice de transfert mesurée lors du pième positionnement des 5 magnétomètres 3 dans le générateur 2 est définie de la manière suivante : H[p] = (1171n(X[i, p])) (6) 15i5nb magnétomètres,15m56 Où i correspond ici au nombre de magnétomètres 3 sur la carte capteurs à calibrer, où p correspond au nombre de positionnement de la carte capteurs comprenant les magnétomètres 3 dans le générateur 2, et où m correspond 10 au nombre de bobinages 21 du générateur 2. On obtient donc l'équation suivante : 1-171n(X[i, p]) =k S[p]Jkn2c (0 [p]n + Q [pro' 13 [ir , kt) (7) j La matrice 1171n(X[i, p]) est la matrice de transfert qui est déterminée lors de 15 l'étape E4, et plus précisément lors de la p-ième mise en oeuvre de l'étape E4 par le procédé. Nous allons détailler ici comment les modèles de transferts mesurés peuvent être calculés. 20 En disposant d'une série de N échantillons de mesures simultanées du champ magnétique par n magnétomètres : (yi(ti)) 15i5n,15j5N Ces mesures dépendant linéairement à chaque instant de p variables qui représentent les mesures des courants électriques dans les bobinages du générateur : (xk(ti)) 151c5p,15j5N 25 On obtient donc pour tout i,j: y(t1) = fi (xi(ti), x2(ti), , xp(ti)) (8) Pour les cas linéaires : 17 Les fonctions, ou transferts f", sont des fonctions linéaires de (ti), x2(ti), , xp (ti)). C'est à dire qu'il est possible d'écrire MO combinaisons linéaires des (ti), x2(ti), , xp(ti)): 5 y, (ti) = Tuxi(ti) T2x2(ti) Tuxp(ti) (9) La matrice T = (Ti i REVENDICATIONS1. Procédé de calibration d'au moins un magnétomètre (3), le procédé étant caractérisé en ce qu'il comprend les étapes suivantes : El : placer au moins un magnétomètre (3) à l'intérieur d'un générateur (2) d'un champ magnétique, ledit générateur (2) comprenant une pluralité de bobinages (21) ; E2: générer un champ magnétique dans le générateur (2) en appliquant des courants électriques aux différents bobinages (21) ; E3: mesurer le champ magnétique généré dans le générateur (2) avec ledit au moins un magnétomètre (3), et mesurer les courants électriques appliqués aux bobinages du générateur (2) ; E4: répéter l'étape E3 en modifiant à chaque fois une position dudit au 15 moins un magnétomètre (3) à l'intérieur du générateur (2) et/ou les courants électriques dans les bobinages (21) du générateur (2) de sorte à obtenir des mesures du champ magnétique et des courants électriques appliqués aux bobinages pour une pluralité de positions dudit au moins un magnétomètre (3) et/ou de courants électriques dans les bobinages (21) ; 20 E5: générer un modèle de transfert paramétrique à partir d'un modèle de mesure paramétrique dudit au moins un magnétomètre (3) et d'un modèle paramétrique du champ magnétique, le modèle de mesure paramétrique dudit au moins un magnétomètre (3) comprenant des paramètres de calibration dudit au moins un magnétomètre (3) ; 25 E6: calculer les paramètres de calibration du au moins un magnétomètre (3) avec un algorithme d'optimisation à partir du modèle de transfert paramétrique généré à l'étape E5 et des mesures du champ magnétique et des courants électriques réalisées dans les étapes E3 et E4. 30 2. Procédé selon la revendication 1, dans lequel le procédé comprend l'étape suivante : 25 E40: calculer un modèle de transfert mesuré entre la mesure du champ magnétique par le au moins un magnétomètre et les courants électriques appliqués aux bobinages (21) pour les différentes positions dudit au moins un magnétomètre (3) et/ou les différents courants électriques dans les 5 bobinages (21), à partir d'une part des mesures dudit champ magnétique et d'autre part des mesures desdits courants électrique appliqués aux bobinages (21) ; le calcul des paramètres de calibration étant réalisé avec un algorithme d'optimisation à partir du modèle de transfert paramétrique et des modèles de transfert mesurés. 3. Procédé selon l'une quelconque des revendications précédentes, dans lequel le procédé comprend en outre une étape de contrôle la température dudit au moins un magnétomètre (3) à calibrer. 4. Procédé selon la revendication 2 ou la revendication 3 prise dans sa dépendance selon la revendication 2, dans lequel le au moins un magnétomètre (3) possède une réponse linéaire, le modèle de transfert mesuré étant une matrice de transfert mesurée qui est calculée par régression linéaire, le modèle de transfert paramétrique étant une matrice de transfert paramétrique. 5. Procédé selon la revendication 2 ou la revendication 3 prise dans sa dépendance selon la revendication 2, dans lequel le au moins un magnétomètre (3) possède une répondre non-linéaire, le modèle de transfert mesuré étant calculé avec un modèle polynomial ou un modèle en séries de Fourier. 6. Procédé selon l'une quelconque des revendications précédentes, dans lequel le générateur (2) est un dispositif à bobines de Helmholtz triaxes comprenant six bobinages (21). 26 7. Procédé selon l'une quelconque des revendications précédentes, dans lequel le procédé comprend une étape d'identification d'un biais dudit au moins un magnétomètre (3), l'étape d'identification étant réalisée en plaçant ledit au moins un magnétomètre (3) dans un dispositif 5 isolant ledit au moins un magnétomètre (3) de tout champ magnétique. 8. Procédé selon l'une quelconque des revendications précédentes, dans lequel le procédé est appliqué sur une pluralité de magnétomètres (3), les paramètres de calibration des magnétomètres (3) 10 comprenant un emplacement des magnétomètres (3) défini par rapport à un des magnétomètres (3) servant de référence. 9. Procédé selon l'une quelconque des revendications précédentes, dans lequel ledit au moins un magnétomètre (3) est maintenu 15 fixe durant l'étape E2 de génération du champ magnétique à l'intérieur du générateur (2) et durant l'étape E3 de mesure du champ magnétique et du courant électrique appliqué à chacun des bobinages (21). 10. Procédé selon l'une quelconque des revendications 20 précédentes, dans lequel les paramètres de calibration du au moins un magnétomètre comprennent : une position dudit au moins un magnétomètre (3) sur un support sur lequel ledit au moins un magnétomètre est monté, un facteur directionnel de mesure dudit au moins un magnétomètre (3), et/ou un biais dudit au moins un magnétomètre (3). 25 11. Procédé selon l'une quelconque des revendications précédentes, dans lequel le modèle paramétrique du champ magnétique est un modèle de champ d'ordre 1, un modèle de champ d'ordre 2, ou un modèle à éléments finis. 30 12. Equipement de calibration d'au moins un magnétomètre (3) comprenant au moins un magnétomètre (3), un générateur (2) d'un champ 27 magnétique qui comprend une pluralité de bobinages (21), caractérisé en ce qu'il comprend également un capteur de courant (5) configuré pour mesurer le courant dans chacun des bobinages (21), et des moyens de traitement (4) qui sont configurés pour : 5 - générer un modèle de transfert paramétrique à partir d'un modèle de mesure paramétrique dudit au moins un magnétomètre (3) et d'un modèle paramétrique du champ magnétique, le modèle de mesure paramétrique dudit au moins un magnétomètre (3) comprenant des paramètres de calibration dudit au moins un magnétomètre (3) ; 10 - calculer les paramètres de calibration du au moins un magnétomètre (3) avec un algorithme d'optimisation à partir du modèle de transfert paramétrique et de mesures du champ magnétique par ledit au moins un magnétomètre (3) et de mesures des courants dans les bobinages (21). 15 13. Produit programme d'ordinateur comprenant des instructions de code pour l'exécution d'un procédé de calibration d'au moins un magnétomètre (3) selon l'une des revendications 1 à 11, lorsque ledit programme est exécuté sur un ordinateur. 20 14. Moyen de stockage lisible par un équipement informatique sur lequel un produit programme d'ordinateur comprend des instructions de code pour l'exécution d'un procédé de calibration d'au moins un magnétomètre (3) selon l'une des revendications 1 à 11.