FR 3069329 A1 20190125 FR 1756974 A 20170721 PROCEDE ET DISPOSITIF DE MESURE DU CHAMP MAGNETIQUE PAR DES MAGNETOMETRES DOMAINE TECHNIQUE GÉNÉRAL ET ART ANTÉRIEUR La présente invention concerne de façon générale la mesure du champ magnétique par des magnétomètres. Elle trouve en particulier avantageusement application dans le domaine de la navigation sans GPS, notamment pour l'estimation du mouvement d'un objet par des techniques magnéto-inertielles. Généralités Les techniques dites de navigation magnéto-inertielles permettent d'estimer précisément la vitesse (et de là, par intégration, la position) d'un solide dans un lieu où le champ magnétique est perturbé. Notamment, il a été proposé dans le brevet FR2914739 d'utiliser, en complément des équations déjà utilisées en techniques inertielles classiques, l'équation É = -f2 A B + V'B.V pour estimer le vecteur vitesse (linéaire) V, à partir du vecteur vitesse angulaire .0 (généralement mesuré par des gyromètres de la centrale inertielle), de la mesure d'un champ magnétique B, de ses dérivées par rapport au temps É et de ses dérivées spatiales VB (gradient), mesurées ou estimées. Grâce à cette information sur la vitesse du solide, l'erreur sur l'estimation de position sera dans le pire des cas proportionnelle au temps, contrairement aux techniques inertielles classiques où elle est quadratique 25 en temps. Ces techniques magnéto-inertielles apportent entière satisfaction, et requièrent des capteurs inertiels de moindre performance que les méthodes inertielles classiques (basées sur l'intégration au travers d'un filtre de Kalman ou d'un autre observateur d'informations en provenance 30 de gyromètres et d'accéléromètres embarqués sur le solide dont on souhaite estimer la position et la vitesse : typiquement, les gyromètres 2 « maintiennent » un repère, dans lequel une double intégration temporelle des mesures des accéléromètres permet d'estimer le mouvement) et donc moins encombrants. Ainsi, par exemple, au moins une centrale inertielle est couplée à un 5 ou plusieurs magnétomètres répartis spatialement (typiquement trois magnétomètres tri-axes disposés sur les sommets d'un trièdre direct) et solidaires de la centrale inertielle (typiquement placée à l'origine du trièdre). Les données issues de ces capteurs sont transmises à une carte électronique d'acquisition et de calcul qui délivre à l'utilisateur 10 l'information de position et de vitesse par résolution de l'équation susmentionnée. L'utilisation de plusieurs magnétomètres permet de mesurer directement les gradients de champ magnétique par exemple par différence finie (développement de Taylor), au lieu de les estimer. On pourra par exemple à cet égard se référer à la demande FR 1653493 de la 15 demanderesse. Dans ce cas, les informations des magnétomètres sont donc utilisées pour calculer des informations de vitesse linéaire ou de vitesse angulaire. L'utilisation la plus classique des magnétomètres dans les systèmes embarqués concerne le calcul du cap. Dans ce cas, on fait l'hypothèse que 20 le champ magnétique mesuré par le capteur est le champ magnétique terrestre qui pointe environ vers le nord. L'information d'attitude du système calculée par ailleurs (avec des accéléromètres et des gyromètres par exemple) est utilisée pour déterminer les composantes horizontales My et Mx du champ magnétique. Le cap psi est alors déterminé comme suit : psi = atan PM ) Mx Ce calcul du cap peut être réalisé en même temps que le calcul de l'attitude dans le cadre d'algorithmes AHRS (« Attitude and Heading Reference System » selon la terminologie anglosaxonne). Le calcul du cap peut également être fait en même temps que les techniques de navigation magnéto-inertielles. 3 Calibra tion Pour que les calculs réalisés à partir des mesures de champ magnétique (les calculs de cap ou les calculs de navigation magnéto-inertielle) soient justes, il convient d'abord de retirer de la mesure des 5 capteurs magnétiques un maximum d'erreurs intrinsèques. Les erreurs classiquement considérées sont les suivantes : le biais, les facteurs d'échelles, les calages d'axes (erreurs géométriques d'orientation propre de chaque axe capteur). La suppression des erreurs des capteurs se fait habituellement au 10 moyen d'une calibration du capteur, cette calibration se déroulant en deux phases. Dans une première phase, on établit- en fonction des erreurs que l'on souhaite calibrer -un modèle dit "de calibration" reliant la mesure brute du capteurMB à la mesure du capteur calibré Mc (ce que mesurerait 15 le capteur s'il était parfait) et les erreurs capteurs. Par exemple, dans le cas d'une calibration du biais B et du facteur d'échelle matriciel A d'un magnétomètre tri-axe, le modèle utilisé peut être le suivant : Mc = A. MB +B Dans une deuxième phase, on met en oeuvre une expérimentation 20 dédiée afin de déterminer les différents paramètres du modèle d'erreur. Dans le cas de l'exemple suscité, il convient d'identifier les différentes composantes de A et de B. De nombreuses méthodes ont déjà été proposées. Une méthode classique consiste à faire tourner le magnétomètre 25 dans un environnement dépourvu de perturbation du champ magnétique terrestre, de manière à assurer qu'il ne mesure que le champ terrestre BT qui est de norme constante. On utilise ensuite des méthodes d'optimisation pour identifier A et B de manière à ce que la norme de Mc soit égale à la norme du champ magnétique terrestre. La solution est donc 30 formellement donnée par : argmin 1111. MB + BI - IBTII A,BEM3,3 (IR) xM3 (IR) 4 On pourra par exemple se référer à l'article suivant : "Complete Triaxis Magnetometer Calibra tion in the Magnetic Domain" - V. Renaudin, M. H. Afzal, and G. Lachapelle - Journal of Sensors - Volume 2010 (2010), Article ID 967245, 10 pages. 5 D'autres techniques encore peuvent utiliser des machines de calibration permettant de générer un champ magnétique maitrisé. Un exemple en ce sens a notamment été proposé dans l'article : "An Automated Indoor Scalar Calibra tion Method for Three-Axis Vector - Magnetometers" - K.Mohamadabadi, M. Hillion - IEEE Sensors 10 Journal, Institute of Electrical and Electronics Engineers (IEEE), 2014, 14, 3076-3083 Le système qui y est décrit est composé de bobines de Helmholtz triaxiales pilotées en courant de sorte que le champ généré soit de norme constante. L'ensemble est enfermé dans un blindage magnétique isolant le 15 système de l'extérieur, permettant le calcul simultané des biais des capteurs. Non linéarités Comme tous les capteurs, les magnétomètres connaissent des 20 problèmes de non linéarités. Le traitement par calibration des erreurs de non linéarités est possible, mais est complexe et non nécessairement totalement fiable. On comprend en effet que plus le nombre de paramètres à identifier est important, plus l'optimisation devient complexe à faire converger. 25 En outre, les non linéarités ont un impact sur la mesure plus faible que d'autres paramètres. Typiquement, par exemple, pour des magnétomètres AMR, les biais pourront avoir un impact de l'ordre de 10 à 100mG sur la mesure, tandis que les facteurs d'échelle impacteront la mesure de l'ordre de 1 à 10mG et les non linéarités uniquement de l'ordre 30 de 0.1 à 1mG. 5 PRÉSENTATION GÉNÉRALE DE L'INVENTION Un but général de l'invention est de proposer une technique de mesures par magnétomètres permettant de s'affranchir d'une calibration des non linéarités, tout en offrant une grande précision et une grande 5 fiabilité de mesure. A cet effet, l'invention propose un procédé de mesure de champ magnétique, dans lequel un champ magnétique est mesuré par un ensemble de magnétomètres qui sont liés à un même support mobile et qui présentent différentes orientations de directions propres par rapport à 10 ce support, caractérisé en ce que des moyens de traitement mettent en oeuvre, lors d'une mesure dans une position et orientation donnée du support et de l'ensemble de magnétomètres, - une détermination, pour chaque magnétomètre d'un champ magnétique "candidat" (en considérant que les capteurs magnétiques 15 n'ont pas de non-linéarités), - une détermination, pour chaque magnétomètre, de l'écart d'orientation entre les directions propres du magnétomètre et le champ magnétique « candidat », - une détermination du magnétomètre pour lequel cet écart est 20 minimal, le champ magnétique mesuré par ce magnétomètre étant retenu comme champ magnétique mesuré par l'ensemble de magnétomètres. Le champ magnétique « candidat » déterminé pour un capteur est par exemple le champ magnétique du capteur ou une combinaison des champs magnétiques de différents capteurs. 25 Les orientations des directions propres dépendent de chaque magnétomètre. La répartition des magnétomètres avec un maximum d'orientations différentes dans l'espace permet de trouver toujours un magnétomètre qui a une orientation la plus optimale et d'utiliser sa mesure comme mesure 30 du champ magnétique. On s'affranchit ainsi des non linéarités des magnétomètres sans avoir à recourir à une calibration spécifique pour traiter les non linéarités. 6 L'invention propose également un dispositif de mesure de champ magnétique, ainsi qu'un système de navigation magnéto-inertielle. Elle concerne également un produit programme d'ordinateur comprenant des instructions de code pour l'exécution d'un tel procédé de 5 mesure de champ magnétique, ainsi qu'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 tel procédé. PRÉSENTATION DES FIGURES 10 D'autres caractéristiques et avantages de l'invention ressortiront encore de la description qui suit, laquelle est purement illustrative et non limitative, et doit être lue en regard des figures annexées sur lesquelles : -la figure 1 est un schéma d'équipements pour la mise en oeuvre du procédé selon l'invention ; 15 - la figure 2 représente plus en détail un exemple de boitier pour la mise en oeuvre du procédé selon l'invention ; - la figure 3a illustre une configuration d'un dispositif conforme à un mode de réalisation de l'invention dans laquelle deux magnétomètres sont répartis dans un même plan ; 20 - la figure 3b illustre une configuration plane avec trois magnétomètres ; - la figure 4 illustre une configuration avec quatre magnétomètres répartis dans l'espace. 25 DESCRIPTION D'UN OU PLUSIEURS MODES DE MISE EN OEUVRE ET DE RÉALISATION Généralités En référence à la figure 1, le dispositif de mesure proposé est par exemple utilisé pour l'estimation du mouvement d'un objet 1 évoluant 30 dans un champ magnétique ambiant (typiquement le champ magnétique 7 terrestre, le cas échéant légèrement altéré par les objets métalliques à proximité), noté Comme on le sait, le champ magnétique est un champ vectoriel tridimensionnel, c'est-à-dire associant un vecteur de dimension trois à chaque point de tridimensionnel dans lequel l'objet est mobile. 5 Cet objet 1 peut être n'importe quel objet mobile dont la connaissance de la position est souhaitée, par exemple un véhicule à roues, un drone, etc., mais également un piéton. L'objet 1 est comporte dans un boitier 2 (support) une pluralité de capteurs de mesure magnétique 20, i.e. des magnétomètres axiaux 20. 10 Par magnétomètre axial, on entend un élément capable de mesurer une composante dudit champ magnétique, i.e. la projection dudit vecteur champ magnétique fi au niveau dudit magnétomètre 20 selon son axe. Plus précisément, les magnétomètres 20 sont solidaires du boitier 2. 15 Ils présentent un mouvement sensiblement identique au boitier 2 et à l'objet 1 dans le référentiel terrestre. De façon préférée, le référentiel de l'objet 1 est muni d'un repère cartésien orthonormé dans lequel les magnétomètres 20 présentent une position prédéterminée dans ce repère. 20 Sur la figure 2, le boitier 2 est fixé sur l'objet 1 (par exemple un membre du piéton) par des moyens d'attache 23. Ces moyens d'attache 23 consistent par exemple en un bracelet par exemple à bande auto-agrippante qui enserre le membre et permet la liaison solidaire. Bien entendu, l'invention n'est pas limitée à l'estimation du 25 mouvement d'un piéton, mais elle est particulièrement avantageuse dans une telle utilisation car elle permet un encombrement très réduit, ce qui est nécessaire pour que le boitier soit portable par un humain de façon ergonomique. 30 Le boitier 2 peut comprendre des moyens de traitement 21 (typiquement un processeur) pour la mise en oeuvre directement en temps réel des traitements du présent procédé, ou bien les mesures peuvent être émises via des moyens de communication 25 vers un dispositif externe tel 8 qu'un terminal mobile (smartphone) 3, voire un serveur distant 4, ou encore les mesures peuvent être enregistrées dans des moyens de stockage de données 22 locaux (une mémoire par exemple de type flash) mémoire locale pour un traitement a posteriori par exemple sur le 5 serveur 4. Les moyens de communication 25 peuvent mettre en oeuvre une communication sans fil à courte portée par exemple Bluetooth ou Wifi (en particulier dans un mode de réalisation avec un terminal mobile 3) voire être des moyens de connexion à un réseau mobile (typiquement 10 UMTS/LTE) pour une communication à longue distance. Il est à noter que les moyens de communication 25 peuvent être par exemple une connectique filaire (typiquement USB) pour transférer les données des moyens de stockage de données 22 locaux à ceux d'un terminal mobile 3 ou d'un serveur 4. 15 Si c'est un terminal mobile 3 (respectivement un serveur 4) qui héberge « l'intelligence », il comprend des moyens de traitement 31 (respectivement 41) tels qu'un processeur pour la mise en oeuvre des traitements du présent procédé qui vont être décrits. Lorsque les moyens de traitement utilisés sont ceux 21 du boitier 2, celui-ci peut encore 20 inclure des moyens de communication 25 pour transmettre la position estimée. Par exemple la position du porteur peut être envoyée au terminal mobile 3 pour afficher la position dans une interface d'un logiciel de navigation. Les moyens de traitement de données 21, 31, 41 respectivement du 25 boitier 2, d'un smartphone 3 et d'un serveur distant 4 peuvent indifféremment et selon les applications réaliser tout ou partie des étapes du procédé. Ils comportent à cet effet chacun des moyens de stockage dans lesquels sont mémorisés tout ou partie de séquences d'instructions de 30 code pour l'exécution du procédé. 9 Disposition et orientation des magnétomètres dans l'espace Chaque capteur magnétomètre, quelle que soit sa structure, 5 présente des directions appelées directions propres dans la suite du texte, pour lequel l'impact des linéarités est minimal. Les magnétomètres 20 sont disposés dans le boitier 2 de façon à couvrir une multitude d'orientations différentes pour ces directions propres. 10 Dans le cas illustré sur les figures 3a et 3b, on a représenté le cas de magnétomètres bi-axiaux 20 présentant chacun deux axes propres X et Y (en l'occurrence perpendiculaires). Les plans formés par les axes propres X et Y de chacun de ces magnétomètres 20 sont parallèles (coplanaires dans les exemples de ces 15 figures). Sur la figure 3a, le capteur comporte deux magnétomètres bi-axiaux 20 co-planaires. Ces magnétomètres 20 sont répartis de façon que leurs axes propres soient séparés angulairement de 45°. On comprend en effet qu'une telle répartition optimise la probabilité 20 qu'un champ magnétique dans le plan XY soit angulairement proche de la direction propre de l'un des capteurs. Dans le cas illustré sur la figure 3b, le système comporte trois magnétomètres bi-axiaux répartis selon une disposition coplanaire des différents axes propres X et Y des magnétomètres 20. Ces magnétomètres 25 20 sont orientés de façon que les axes propres de deux magnétomètres successifs soient séparés angulairement de 30°. On comprend là aussi qu'une telle répartition optimise la probabilité qu'un champ magnétique dans le plan XY soit angulairement proche de la direction propre de l'un des capteurs. 30 Bien entendu ces deux exemples se généralisent au cas de m magnétomètres bi-axiaux coplanaires ou répartis dans des plans parallèles du dispositif (avec m, nombre entier). 10 Dans ce cas la répartition optimale est une répartition de magnétomètres séparés par des angles de 90°/m, mais on comprendra bien que l'invention ne se limite pas à cette répartition optimale. La figure 4 illustre quant à elle le cas de quatre magnétomètres tri- 5 axiaux présentant chacun trois axes propres X, Y et Z perpendiculaires. Elle illustre une répartition des orientations de ces quatre magnétomètres dans l'espace qui permet une optimisation de la proximité angulaire entre un champ magnétique B à mesurer et au moins une direction propre de l'un des magnétomètres 20. 10 De façon plus générale, dans le cas de magnétomètres tri-axes, la répartition de leur orientation est choisie avec un maximum d'orientations dans l'espace afin de maximiser la probabilité de trouver un magnétomètre qui a une orientation la plus optimale (optimisation de la répartition de leur directions propres). 15 Choix du magnétomètre pour une mesure donnée Pour chaque magnétomètre j, on désigne par (d) les n ie[[i..n]] directions propres de mesure de ce capteur, c'est-à-dire les directions 20 liées à ce capteur pour lesquelles l'impact des non-linéarités est minimal. Le nombre de ces directions propres dépend du type capteur. Les moyens de traitement calculent pour chaque magnétomètre un champ magnétique candidat en considérant que les magnétomètres ne 25 sont pas perturbés par les non-linéarités. En effet, comme l'effet des non-linéarités est faible, on peut supposer en première approximation que le champ magnétique fourni par chaque capteur n'est pas trop mauvais, en tout cas, qu'il suffit pour la suite du procédé. Le champ magnétique candidat d'un capteur peut être le champ magnétique mesuré par ce 30 capteur ou une combinaison astucieuse (comme la moyenne par exemple) des champs magnétiques mesurés par tous les capteurs. Dans la suite, on 11 considèrera le cas, plus simple, où le champ magnétique candidat d'un capteur est le champ magnétique mesuré par ledit capteur. Les moyens de traitement (selon le cas 21, 31 ou 41) calculent pour chaque capteur l'écart d'orientation entre le champ magnétique mesuré 5 par le capteur et ses directions propres. Si BBJ est le champ magnétique « candidat » pour le capteur j, l'écart el entre les directions propres et est donné par = min (i',) iEft1..fl IL B où (c/',BBI) est l'angle entre les deux vecteur ",et BBJ. 10 Les moyens de traitement 21, 31 ou 41 déterminent ensuite le capteur] pour lequel l'écart el est minimal. C'est de ce capteur qu'on obtient la meilleure mesure, au sens où l'impact des non-linéarités est le plus faible sur la mesure. Le champ magnétique mesuré par le magnétomètre j ainsi 15 déterminé est donc retenu par les moyens de traitement 21, 31, 41 comme champ magnétique mesuré par l'ensemble de magnétomètres. La mesure de champ du magnétomètre j est stockée par lesdits moyens de traitement21, 31, 41 et/ou utilisée par lesdits moyens de traitement dans les calculs d'informations de navigation magnéto- 20 inertielles (vitesse linéaire, vitesse angulaire, position, cap, etc...). Dans ces calculs, la mesure du champ j est alors utilisée comme mesure du champ magnétique à l'instant considéré, correspondant à une position / orientation de l'ensemble de magnétomètres. Les traitements sont alors avantageusement du type de ceux décrits 25 dans le brevet FR2914739 ou dans la demande FR 1653493. Comme on l'aura compris, le fait d'utiliser le champ magnétique « candidat » du magnétomètre pour lequel l'écart est minimal permet de disposer d'une mesure non perturbée par les non linéarités des magnétomètres et ce sans avoir eu à recourir à une calibration spécifique 30 pour traiter celles-ci. Bien entendu, le procédé et le dispositif de calibration proposés peuvent être utilisés en combinaisons avec des calibrations destinées à 12 corriger non pas les erreurs de non-linéarité, mais d'autres paramètres d'erreur du capteur selon les méthodes connues. Directions propres 5 Les directions propres sont déterminées préalablement pour chaque capteur. Cette détermination peut se faire par exemple à partir de la modélisation des termes d'erreurs du capteur considéré, ou par une analyse de la physique du capteur et/ou de sa mise en oeuvre. 10 On donne ci-après un exemple modélisation et de détermination de directions propres, d'autres modélisations et détermination étant bien entendu possibles. Cas des magnétomètres AMR - modélisation des non linéarités 15 Les magnétomètres 20 sont par exemple chacun composé de trois puces identiques dont le fonctionnement repose sur la technologie dite AMR (« Anisotropic Magneto-Resistance » selon la terminologie anglosaxonne). 20 Lorsque ces puces sont plongées dans un champ magnétique, elles renvoient une tension qui s'approche d'une fonction affine de la composante du champ magnétique dans la direction de leur axe de sensibilité. Classiquement, les mesures magnétométriques d'un capteur d'axe 25 de sensibilité X et de plan principal XY (avec Y direction perpendiculaire à X) sont modélisées comme suit (voir par exemple Handbook of Magnetic Measurements, S. Tumanski) : 1+ mi.mo Où M'est le champ magnétique mesuré par le capteur et oùMi.' et 30 /le sont les champs magnétiques selon les axes X et Y. Ho est le champ rémanent, supposé grand devant les champs magnétiques en jeu. 13 On peut former un magnétomètre triaxe AMR. Un tel magnétomètre s'analyse comme formé de 3 capteurs mono-axes, dont le modèle de non linéarité est 5 Cas des magnétomètres AMR -Détermination des directions propres En fonction du champ magnétique mesuré, l'impact des non linéarité n'est pas le même sur la mesure du capteur. Ainsi, on lit directement sur le modèle que : 10 Lorsque le champ est aligné avec les axes capteurs, les non linéarités ne gênent pas la mesure du capteur. En effet : Si le champ Mic est a peu près aligné avec l'axe capteur X, on a le Mc4 0, ce qui donne me - O me 0 15 Si le champ le est a peu près aligné avec l'axe capteur Y, on a Mic Mc4 0, ce qui donne me O me O Si le champ Mc4 est a peu près aligné avec l'axe capteur Z, on a Mic 0, ce qui donne me 0 20 MO 111,6 Mc4 Lorsque le champ n'est pas aligné avec les axes capteurs, les non linéarités gênent la mesure du capteur. Dans le cas d'un champ magnétique orienté à 45° par rapport aux axes X et Y, on a alors 14 V2 Mê, , V2 Mi. Mc4 ,-,' 0 Ce qui donne un champ mesuré 1Mc2 mx = V: Mc V2 V2 - 2 Ho 1+-2 Mc/Ho 2 Mc 5 V2 Mc V2 1 Mc2 MC + i Ho V2 2 1- -2MC/HO 0 Les directions optimales (axes REVENDICATIONS1. Procédé de mesure de champ magnétique, dans lequel un champ magnétique est mesuré par un ensemble de magnétomètres qui sont liés à un même support mobile et qui présentent des différentes orientations de directions propres par rapport à ce support, caractérisé en ce que des moyens de traitement mettent en oeuvre, lors d'une mesure dans une position et orientation donnée du support et de l'ensemble de magnétomètres, - une détermination, pour chaque magnétomètre d'un champ magnétique « candidat », - une détermination, pour chaque magnétomètre, de l'écart d'orientation entre les directions propres du magnétomètre et le champ magnétique « candidat », - une détermination du magnétomètre pour lequel cet écart est minimal, le champ magnétique mesuré par ce capteur étant retenu comme champ magnétique mesuré par l'ensemble de magnétomètres. 2. Procédé selon la revendication 1, dans lequel le champ magnétique « candidat » pour un magnétomètre est le champ magnétique mesuré par ledit magnétomètre. 3. Procédé selon l'une des revendications précédentes dans lequel, pour chaque magnétomètre, les moyens de traitement calculent = min (i',) iEft1..fl IL B OÙ - j est un indice désignant le magnétomètre, mesure de ce - par (d') les n directions propres de ie[[1..n]] magnétomètre j, 16 - .6', est le champ magnétique « candidat » déterminé pour ledit magnétomètre, - (d',BBJ) est l'angle entre les deux vecteur d', et .6',.3 4. Procédé selon l'une des revendications précédentes, dans lequel 5 les directions propres d'un magnétomètre sont déterminées par modélisation, préalablement aux mesures. 5. Procédé selon l'une des revendications précédentes, dans lequel les moyens de traitement calculent des informations de navigation magnétiques en utilisant le champ magnétique mesuré par le 10 magnétomètre pour lequel l'écart aux directions propres est minimal, comme mesure du champ magnétique de l'ensemble de magnétomètres. 6. Dispositif de mesure de champ magnétique comportant un support mobile et un ensemble de magnétomètres liés audit support, 15 lesdits magnétomètres présentant différentes orientations de direction propres dans celui-ci, ledit dispositif comportant en outre des moyens de traitement adaptés pour mettre en oeuvre, lors d'une mesure dans une position et orientation donnée de l'ensemble de magnétomètres, 20 - une détermination, pour chaque magnétomètre d'un champ magnétique « candidat », - une détermination, pour chaque magnétomètre, de l'écart d'orientation entre les directions propres du magnétomètre et le champ magnétique « candidat », 25 - une détermination du magnétomètre pour lequel cet écart est minimal, le champ magnétique mesuré par ce capteur étant retenu comme champ magnétique mesuré par l'ensemble de magnétomètres. 17 7. Dispositif selon la revendication 6, comportant une pluralité de magnétomètres plans à deux axes propres, ces magnétomètres étant parallèles, leurs orientations étant réparties pour que les directions de leurs axes propres soient séparées angulairement 5 d'un angle de 90°/m, ou m est le nombre de ces magnétomètres. 8. Dispositif selon la revendication 6, dans lequel lesdits magnétomètres sont parallèles. 9. Dispositif selon la revendication 6, dans lequel les magnétomètres sont de type tri-axes, les orientations de ces magnétomètres étant 10 répartis afin d'optimiser la répartition de leurs directions propres dans l'espace. 10. Système de navigation comportant au moins un dispositif de mesure de champ magnétique selon l'une des revendications 6 à 9. 15 11. Produit programme d'ordinateur comprenant des instructions de code pour l'exécution d'un procédé de mesure de champ magnétique selon l'une des revendications 1 à 5, lorsque ledit programme est exécuté sur un ordinateur. 12. Moyen de stockage lisible par un équipement informatique sur 20 lequel un produit programme d'ordinateur comprend des instructions de code pour l'exécution d'un procédé de mesure de champ magnétique selon l'une des revendications 1 à 5. 25