FR 3029642 A1 20160610 FR 1462106 A 20141209 1 Procédé de localisation d'au moins un objet magnétique mobile, et système associé La présente invention porte sur un procédé de localisation d'au moins un objet magnétique mobile, et un système associé. Il est connu des documents portant sur des procédés et dispositifs de localisation d'un objet magnétique, tel le document FR 2 988 862. Le document FR 2 988 862 divulgue un procédé de localisation d'un objet magnétique a partir de mesures d'un réseau de magnétomètres tri-axes, ce réseau comportant N magnétomètres tri-axes reliés mécaniquement les uns aux autres sans aucun degré de liberté pour conserver une distance connue entre chacun de ces magnétomètres, dans lequel N est un nombre entier supérieur à cinq. Le réseau de magnétomètres est compris dans un dispositif comprenant une surface d'appui pour un support d'écriture tel une ou plusieurs feuilles, reliées ou non. La localisation d'un objet magnétique peut être considéré comme sa position (coordonnées (x, y, z) dans un repère fixe par rapport au réseau de magnétomètres) et la valeur de son moment magnétique ((Mx, My, Mz) dans le même repère) représentatif de l'orientation de l'objet magnétique. L'introduction d'un bais sur les mesures utilisées dans ce type de procédé rend la localisation imprécise, voire instable Dans le cadre d'une utilisation mobile d'un tel dispositif de localisation, i.e. lorsque l'utilisateur bouge le dispositif de localisation supportant son support d'écriture, son orientation est modifiée, ce qui entraîne directement une rotation du champ magnétique terrestre, dans le repère mesuré par le réseau de magnétomètres, et par conséquent un biais de mesure sur l'ensemble des magnétomètres. 3029642 2 Le mouvement du dispositif de localisation dans un champ magnétique homogène rend donc quasiment impossible l'utilisation d'un tel dispositif de localisation. 5 L'utilisation d'un tel dispositif dans un système de numérisation d'une trace manuscrite pour lequel la précision de la localisation recherchée est inférieure à 1 mm, implique qu'une légère variation du champ magnétique terrestre dans le repère des magnétomètres, entraîne directement un écart de mesure supérieur à la précision recherchée. 10 Il est connu de l'état de l'art, de supprimer, au démarrage, le biais de mesure sur l'ensemble des capteurs, ou en l'absence de tout objet magnétique. 15 Un tel procédé ne permet pas de mouvement pendant la localisation d'un objet magnétique et ne résout donc pas la problématique mentionnée pour une utilisation mobile. Notamment, dans le cadre d'une application de numérisation d'écriture Il arrive fréquemment que l'utilisateur bouge le dispositif de localisation afin d'obtenir un angle d'écriture possible 20 tout en se déplaçant et modifiant sa posture. Il est également connu de placer un magnétomètre à une distance suffisante de l'objet magnétique mobile, de sorte que ce magnétomètre ne soit pas affecté par le champ magnétique généré par l'objet magnétique 25 mobile. Le chant mesuré par ce magnétomètre distant sert à supprimer la composante magnétique induite par le champ terrestre. Un tel procédé ne permet pas d'utiliser des dispositifs de petite taille, et ne se prête donc pas un dispositif utilisable en mobilité. En outre, ce 30 procédé est très sensible aux erreurs de mesure de ce magnétomètre distant, notamment en cas de dysfonctionnement ou de perturbation magnétique autre que terrestre de ce dernier. Un but de l'invention est de pallier ces problèmes. 3029642 3 Il est proposé, selon un aspect de l'invention, un procédé de localisation d'au moins un objet magnétique mobile par rapport à un réseau d'au moins N magnétomètres triaxiaux reliés entre-eux mécaniquement sans aucun degré de liberté pour conserver une position relative connue de ces 5 magnétomètres, N étant un nombre entier au moins égal à 2 comprenant les étapes consistant à: - charger les mesures des magnétomètres, à un instant courant; - calculer une moyenne pondérée desdites mesures des magnétomètres représentant un champ magnétique terrestre moyen 10 et un champ magnétique moyen dû au(x) objet(s) magnétique(s) mobile(s); - soustraire ladite moyenne pondérée desdites mesures de chacune desdites mesures pour obtenir des mesures modifiées; - charger lesdites mesures modifiées, ainsi qu'une localisation du ou 15 des objets magnétiques mobiles à l'instant courant en entrée d'un filtrage de localisation du ou des objets magnétiques mobile; - effectuer le filtrage de localisation comprenant les étapes consistant à: - calculer ladite moyenne pondérée d'estimations des données délivrées par lesdits magnétomètres; 20 - soustraire ladite moyenne pondérée desdites estimations de chacune desdites données estimées; et - délivrer en sortie une localisation du ou des objets magnétiques mobile à un instant suivant. 25 Un tel procédé permet de pouvoir utiliser un tel dispositif de localisation magnétique tout en en se déplaçant et/ou modifiant l'orientation du dispositif. De plus ce procédé permet de supprimer l'influence d'un champ 30 homogène sur le dispositif, (le champ terrestre en étant un ) ainsi que de supprimer les sources de perturbation qui présentent un champ magnétique homogène sur l'ensemble des magnétomètres. La plupart des sources s'apparentent à un dipôle magnétique. Si la distance entre le dispositif et la source de perturbation est suffisamment importante, le champ magnétique 3029642 4 perçu est un champ magnétique homogène qui est donc corrigé par le procédé. Des sources connues de perturbation magnétiques peuvent être 5 des moteurs à courant continu pour lesquels les champs magnétiques générés sont importants, par exemple supérieurs à 1 Gauss. Ainsi cette méthode permet également l'utilisation du dispositif dans un train ou un tram. On entend par localisation, la détermination de la position et de 10 l'orientation de l'objet magnétique mobile, par exemple en déterminant la position et le moment magnétique de l'objet magnétique avec une représentation cartésienne (composante selon trois axes orthogonaux). En variante, la détermination de la position et de l'orientation de l'objet magnétique mobile peut être effectuée en exprimant le moment magnétique 15 dans un repère sphérique c'est-à-dire selon une norme et deux angles. Cette alternative peut permettre de fixer le moment magnétique de l'objet magnétique mobile et donc estimer uniquement sa position et deux angles en initialisant le moment magnétique à une valeur connue, l'avantage et de diminuer le nombre d'inconnues à estimer pour le filtre de localisation. 20 Dans un mode de mise en oeuvre, ladite moyenne pondérée utilise des coefficients de pondération identiques pour chacun des N magnétomètres. Ce procédé offre l'avantage de nécessiter une puissance de calcul moins importante qu'en utilisant des pondérations différentes. En effet 25 on s'affranchit de 3n multiplications (avec N le nombre de capteurs). Selon un mode de mise en oeuvre, ladite moyenne pondérée utilise des coefficients de pondération plus faibles pour P magnétomètres, parmi les N magnétomètres, les plus proches d'un objet magnétique mobile 30 respectif, P étant un nombre entier inférieur à N. Ainsi, on limite l'influence des P magnétomètres les plus perturbés, voire saturés, par les objets magnétiques mobiles sont ignorés. 3029642 5 Par exemple, s'il y a plusieurs objets magnétiques mobiles, on détermine, pour chaque magnétomètre du réseau, les distances entre le magnétomètre et chaque objet magnétique mobile, puis on classe les distances déterminées pour l'ensemble des magnétomètres et des objets 5 magnétiques mobiles dans un ordre croissant, et on utilise des coefficients de corrélation plus faibles pour les P magnétomètres dont les distances classées sont les plus faibles. Par exemple, ladite moyenne pondérée utilise des coefficients de 10 pondération nuls pour lesdits P magnétomètres. Ainsi, les P magnétomètres les plus perturbés, voire saturés, par les objets magnétiques mobiles sont ignorés. 15 Dans un mode de réalisation, ladite moyenne pondérée utilise des coefficients de pondération nuls sauf pour le magnétomètre, parmi les N magnétomètres, le plus éloigné du ou des objets magnétiques mobiles. Ainsi, on prend pour moyenne du champ magnétique terrestre 20 mesuré, le champ magnétique mesuré par magnétomètre le plus éloigné du ou des objets magnétiques mobiles, donc le moins perturbé par le ou les objets magnétiques mobiles. Il est également proposé, selon un autre aspect de l'invention un 25 système de localisation d'au moins un objet magnétique mobile comprenant: - un réseau d'au moins magnétomètres triaxiaux reliés entre-eux mécaniquement sans aucun degré de liberté pour conserver une position relative connue de ces magnétomètres, N étant un nombre entier au moins égal à 2; 30 - une unité électronique de traitement adapté pour déterminer la localisation du ou des objets magnétiques mobiles à partir de mesures fournies par les magnétomètres du réseau et d'un filtre de localisation; ladite unité électronique de traitement comprenant: - des moyens de chargement des mesures des magnétomètres, à un 35 instant courant; 3029642 6 - des moyens de calcul d'une moyenne pondérée desdites mesures des magnétomètres représentant un champ magnétique terrestre moyen et un champ magnétique moyen dû au(x) objet(s) magnétique(s) mobile(s); et 5 - des moyens de calcul d'une soustraction de ladite moyenne pondérée desdites mesures de chacune desdites mesures pour obtenir des mesures modifiées; le filtre de localisation comprenant : - des moyens de chargement desdites mesures modifiées, ainsi que 10 d'une localisation du ou des objets magnétiques mobiles à l'instant courant en entrée d'un filtrage de localisation de l'objet magnétique mobile; - des moyens de calcul d'une moyenne pondérée d'estimations des données délivrées par lesdits magnétomètres; 15 - des moyens de calcul d'une soustraction de ladite moyenne pondérée desdites estimations de chacune desdites données estimées; et - des moyens de délivrance en sortie d'une localisation du ou des objets magnétiques mobile à un instant suivant. 20 L'invention sera mieux comprise à l'étude de quelques modes de réalisation décrits à titre d'exemples nullement limitatifs et illustrés par les dessins annexés sur lesquels : - les figures 1 et 2 illustrent schématiquement un dispositif de localisation selon l'état de la technique connu; et 25 - la figure 3 illustre schématiquement un procédé selon un aspect de l'invention. Sur l'ensemble des figures, les éléments ayant des références identiques sont similaires. 30 Dans la présente description, les caractéristiques et fonctions bien connues de l'homme du métier ne sont pas décrites en détails. 3029642 7 La figure 1 représente un dispositif de localisation 1 d'un objet magnétique mobile OMM selon l'état de l'art, par exemple décrit dans le document FR 2988862 A1. 5 Le dispositif de localisation 1, en l'espèce, commande un appareil électronique 2, par exemple un écran 2 raccordé à une unité de commande 3 capable de commander l'affichage d'une image sur cet écran 2. Sur l'exemple représenté, l'unité 3 commande la position et 10 l'orientation d'un curseur 4 sur l'écran 2. Par exemple, le curseur 4 a une forme tridimensionnelle, telle un parallélépipède rectangle. Un objet magnétique mobile OMMk comprend un aimant permanent 5 présentant un moment magnétique non nul même en absence 15 de champ magnétique extérieur et un ustensile 6 non magnétique. Par non magnétique, on entend un ustensile réalisé dans un matériau ne présentant aucune propriété magnétique mesurable. Par exemple, l'ustensile 6 est un crayon, un stylo ou un effaceur en bois ou en plastique. L'ustensile 6 peut être localisé à partir de la localisation de l'aimant permanent 5. L'ustensile 6 20 permet aussi de faciliter la manipulation de l'aimant permanent 5, notamment lorsque celui-ci est de petite taille. Typiquement, la combinaison de l'aimant permanent 5 et de l'ustensile 6 est librement déplaçable directement par la main d'un être humain dans un repère orthogonal XYZ fixé sans aucun degré de liberté au dispositif 1. En l'espèce, les directions X et Y sont horizontales 25 (plan du support d'écriture) et la direction Z est verticale. A cet effet, cette combinaison pèse moins d'un kilo et, de préférence, moins de 200g. Les dimensions de cette combinaison sont suffisamment réduites pour pouvoir être saisie et déplacée par une seule main d'un utilisateur. Le ou les objets magnétiques mobiles OMMk sont librement déplaçables dans le repère 30 orthogonal XYZ. Par exemple, le champ magnétique coercitif d'un aimant est supérieur à 100 A.m-1 ou 500 A.m-1. Par exemple, il est réalisé en matériau ferro- ou ferrimagnétique. L'aimant permanent 5 a une forme oblongue. Sur 35 la figure 1, la direction du moment magnétique de l'aimant permanent 5 est 3029642 8 représentée par une flèche parallèle à la direction longitudinale de cet objet. En variante, l'aimant permanent 5 peut être en forme d'anneau entourant l'ustensile 6. La plus grande longueur de l'aimant permanent 5 est notée L par la suite. 5 La puissance de l'aimant permanent 5 est typiquement supérieure à 0,01 A.m2 ou 0,1 A.m2. Dans ce mode de réalisation, l'aimant permanent 5 est fixé sans aucun degré de liberté sur l'ustensile 6. 10 Le dispositif de localisation 1 permet de localiser l'aimant permanent 5 dans le repère XYZ. Par localisation, on entend ici la détermination de la position x, y, z de l'aimant permanent 5 dans le repère XYZ et aussi la détermination de l'orientation de l'aimant permanent 5 par rapport aux axes X, Y et Z du repère XYZ. Par exemple, l'orientation de 15 l'aimant permanent 5 peut être représentée par les coordonnées cartésiennes Mx, My et M7, du moment magnétique de l'aimant permanent 5 dans le repère XYZ, ou par les coordonnées sphériques M, 0, cp du moment magnétique de l'aimant permanent 5 dans le repère XYZ, respectivement, par rapport aux axes X, Y et Z du repère. 20 Le dispositif de localisation 1 comporte un réseau de N magnétomètres tri-axes Mg. Sur la figure 1, les traits ondulés verticaux indiquent qu'une partie du dispositif de localisation 1 n'a pas été représentée. 25 Typiquement, N peut être supérieur à cinq et, de préférence, supérieur à seize ou trente-deux, voire égal à soixante-quatre. Dans ce mode de réalisation, les magnétomètres Mg sont alignés en lignes et en colonnes pour former une matrice ou réseau. En l'espèce, 30 cette matrice comporte huit lignes et huit colonnes. Les indices i et j identifient, respectivement, la ligne et la colonne de cette matrice à l'intersection de laquelle se trouve le magnétomètre Mg. 3029642 9 Sur la figure 1, seuls les magnétomètres Mil, Mi2, Mi3, Mi4 et Mi8 d'une ligne i sont visibles. La position des magnétomètres Mg les uns par rapport aux autres est décrite plus en détail sur la figure 2. 5 Chaque magnétomètre Mg est fixé sans aucun degré de liberté aux autres magnétomètres. A cet effet, les magnétomètres Mg sont fixés sans aucun degré de liberté sur une face arrière 7 d'une plaque rigide 8. Cette plaque rigide 8 présente une face avant 9 tournée vers l'aimant permanent 5. La plaque 8 est réalisée dans un matériau non magnétique 10 rigide. Par exemple, la plaque 8 peut être en verre. Chaque magnétomètre Mg mesure la direction et l'intensité du champ magnétique généré par l'aimant permanent 5. Pour cela, chaque magnétomètre Mg mesure la norme de la projection orthogonale du champ 15 magnétique généré par l'aimant permanent 5 au niveau de ce magnétomètre Mg sur trois axes de mesure de ce magnétomètre. En l'espèce, ces trois axes de mesure sont orthogonaux entre eux. Par exemple, les axes de mesure de chacun des magnétomètres Mg sont, respectivement, parallèles aux axes X, Y et Z du repère. 20 Chaque magnétomètre Mg est raccordé par l'intermédiaire d'un bus 10 de transmission d'informations à une unité de traitement 11. L'unité de traitement 11 est capable de déterminer la position et 25 l'orientation de l'aimant permanent 5 dans le repère XYZ à partir des mesures des magnétomètres Mg. A cet effet, l'unité 11 comporte un calculateur électronique 12 programmable apte à exécuter des instructions enregistrées sur un support d'enregistrement d'informations. L'unité 11 comprend donc également une mémoire 13 contenant les instructions 30 nécessaires pour l'exécution par le calculateur 12 du procédé de la figure 3. En particulier, l'unité 11 implémente un modèle physico-mathématique associant des mesures de plusieurs magnétomètres Mg à des paramètres représentatifs de la position et de l'orientation de l'aimant 35 permanent 5, i.e. également de l'ustensile 6, dans le repère XYZ. 3029642 10 Ce modèle est implémenté sous la forme d'un filtre d'estimation non linéaire, par exemple un filtre de Kalman étendu. 5 Ce modèle est typiquement construit a partir des équations physiques de magnétostatique. Pour construire ce modèle, l'aimant permanent 5 est approximé par un dipôle magnétique. Cette approximation n'introduit que très peu d'erreurs si la distance entre l'aimant permanent 5 et le magnétomètre Mg est supérieure à 2L et, de préférence, supérieure à 8L; L 10 étant la plus grande dimension de l'aimant permanent 5. Typiquement, L est inférieur à 20 cm et, de préférence, inférieur à 10 ou 5cm. L'unité 11 est également capable de restituer la position et l'orientation mesurée de l'aimant permanent 5 sur une interface 14. 15 L'unité de commande 3 est raccordée à l'unité 11 par l'intermédiaire de cette interface 14. La figure 2 représente une partie des magnétomètres Mg du 20 dispositif de localisation 1. Ces magnétomètres Mg sont alignés en lignes i parallèles à la direction X. Ces magnétomètres sont également alignés en colonnes j parallèles à la direction Y pour former une matrice ou réseau. Les lignes i et les colonnes j sont disposées dans l'ordre des indices croissants. 25 Le centre du magnétomètre Mg se trouve à l'intersection de la ligne i et de la colonne j. Le centre du magnétomètre correspond au point auquel le champ magnétique est mesuré par ce magnétomètre. En l'espèce, les indices i et j appartiennent à l'intervalle [1 ; 8]. Les centres de deux magnétomètres Mg et M1,1+1 immédiatement 30 consécutifs le long d'une ligne i sont séparés par une distance connue dij,i±i. De façon similaire, le centre de deux magnétomètres Mq et Mi+1,1 immédiatement consécutifs le long d'une même colonne j sont séparés par une distance connue di,i,i±i. 3029642 11 Dans le cas particulier présentement décrit, quelle que soit la ligne i, la distance di,1,1+1 est la même. Cette distance est donc notée di. De façon similaire, quelle que soit la colonne j, la distance di,i,i+i entre deux magnétomètres est la même. Cette distance est donc notée di. Sur l'exemple 5 décrit, les distances di et di sont toutes les deux égales à d. Typiquement, la distance d est comprise entre 1 et 4 cm lorsque 1 - la puissance de l'aimant permanent est de 0,5 A.m2, - la sensibilité des magnétomètres est de 4.10-7T, et 10 - le nombre de magnétomètres Mg est de soixante-quatre. La figure 3 illustre un mode de réalisation d'un procédé selon un aspect de l'invention. Les étapes du procédé sont mises en oeuvre par des moyens de l'unité électronique de traitement 11. A chaque étape de mesure, on calcule une moyenne pondérée desdites mesures Bir" des magnétomètres Mi j respectifs représentant la somme d'un un champ magnétique terrestre moyen et d'un champ magnétique moyen dû au(x) objet(s) magnétique(s) mobile(s) OMMk. La moyenne pondérée peut être, de manière non limitative, la moyenne arithmétique pour laquelle les coefficients de pondération sont tous égaux. 25 Ainsi, le champ magnétique terrestre est supprimé en même temps que le champ magnétique moyen généré par le ou les objets magnétiques mobiles OMMk. La suppression du champ magnétique terrestre ou perturbations 30 magnétique homogène sur le système constitue la correction recherchée tandis que la suppression de la moyenne du champ magnétique généré par la masse ferromagnétique constitue une erreur introduite. On prend ensuite en compte l'erreur introduite pour ne pas 35 dégrader la précision de l'algorithme de localisation. 15 20 3029642 12 [(e")1 On calcule donc la moyenne pondérée = (Br) des (Mnes) [Bmes x ii mesures Bi jmes= lees ii des magnétomètres Mi,1 représentant un champ Mnesii t magnétique terrestre moyen et un champ magnétique moyen dû au(x) 5 objet(s) magnétique(s) mobile(s) OMMk. [Bmes x ii On soustrait ensuite, de chacune des mesures Bires= 13j/nes ii , la Brest] t [(eles) moyenne pondérée = (Br) des mesures Bi gimes , pour obtenir des (Mnes) [B;nes ii - (le es) mesures modifiées Bijmes- = Brii - (B;nes) B;nes ii - (Br es) 10 Par exemple, dans le cas de la moyenne arithmétique comme moyenne pondérée, on a : Brel (lees) (Bres) &nes ij [ = 1x-, x,p Bmes 7 L1=1 j=1 y ij B;nesij avec axb=N Un filtrage de localisation (FL) utilise ensuite en entrée les Brest] (mies) mesures modifiées 13; jmes- = Bytes (Byies^ / et une localisation Bfies (Bfies) 15 3029642 13 (Post, Mt)= Pos t+1 du ou des objets magnétiques mobiles (OMMk) à M (y z l'instant courant t. Le filtrage de localisation utilise un calcul d'une moyenne est (B:st) pe 5 pondérée estimée = (West) d'estimations Bi jest= mstii des données (Bfst) I Brt ii t délivrées par lesdits magnétomètres (M1). En l'espèce la moyenne pondérée utilisée doit être identique à celle utilisée pour la moyenne pondérée mesurée (i.e. mêmes coefficients de pondérations respectivement associées aux magnétomètres dans la pondération des mesures et dans la pondération 10 des estimations, un coefficient de pondération est associé à un magnétomètre). Les algorithmes d'estimation reposent sur l'utilisation d'un modèle de mesures qui lie position, orientation, moment magnétique de l'aimant avec 15 le champ magnétique généré en tout point de l'espace. Le calcul de la moyenne pondérée estimée, dans ce mode de réalisation, utilise la même moyenne pondérée, en l'espèce la moyenne arithmétique: 20 [B:st ij = 7.11 Lia ali 3 1 B3est ij West ij [(B:st)1 (133e,st) avec axb=N (Bfst) (B:st) On soustrait ensuite ladite moyenne pondérée = (mst) de (West) 25 chacune desdites données estimées Birt. 3029642 14 Le filtre de localisation FL délivre ensuite en sortie une localisation 1 Pos (y (Post z o, Mt+i)= x du ou des objets magnétiques mobile OMMk à M (y z - t+i un instant suivant t+1. 5 En variante, ladite moyenne pondérée utilise des coefficients de pondération plus faibles pour P magnétomètres, parmi les N magnétomètres Mij, les plus proches d'un objet magnétique mobile OMMk respectif, P étant un nombre entier inférieur à N. Bien entendu, la moyenne pondérée des mesures est toujours identique à la moyenne pondérée des estimations. 10 Ainsi, on limite l'influence des P magnétomètres les plus perturbés, voire saturés, par les objets magnétiques mobiles sont ignorés. La moyenne pondérée peut même utiliser des coefficients de 15 pondération nuls pour lesdits P magnétomètres, de manière à éliminer totalement l'influence des P magnétomètres les plus perturbés, voire saturés, par les objets magnétiques mobiles sont ignorés. En variante, la moyenne pondérée peut utiliser des coefficients de 20 pondération nuls sauf pour le magnétomètre, parmi les N magnétomètres, le plus éloigné du ou des objets magnétiques mobiles, pour prendre comme moyenne du champ magnétique terrestre, la mesure du magnétomètre le moins perturbé par le ou les objets magnétiques mobiles. 25 Dans ce mode de réalisation, la moyenne calculée correspond exactement à la mesure du magnétomètre. Les étapes du procédé décrit ci-dessus peuvent être effectuées par un ou plusieurs processeurs programmables exécutant un programme 30 informatique pour exécuter les fonctions de l'invention en opérant sur des données d'entrée et générant des données de sortie. 3029642 15 Un programme d'ordinateur peut être écrit dans n'importe quelle forme de langage de programmation, y compris les langages compilés ou interprétés, et le programme d'ordinateur peut être déployé dans n'importe quelle forme, y compris en tant que programme autonome ou comme un 5 sous-programme, élément ou autre unité appropriée pour une utilisation dans un environnement informatique. Un programme d'ordinateur peut être déployée pour être exécuté sur un ordinateur ou sur plusieurs ordinateurs à un seul site ou répartis sur plusieurs sites et reliées entre elles par un réseau de communication. 10 Le mode de réalisation préféré de la présente invention a été décrit. Diverses modifications peuvent être apportées sans s'écarter de l'esprit et de la portée de l'invention. Par conséquent, d'autres mises en oeuvre sont dans la portée des revendications suivantes. 15 REVENDICATIONS1. Procédé de localisation d'au moins un objet magnétique mobile (OMMk) par rapport à un réseau d'au moins N magnétomètres triaxiaux (Mi,j) reliés entre-eux mécaniquement sans aucun degré de liberté pour conserver une position relative connue de ces magnétomètres N étant un nombre entier au moins égal à 2, comprenant les étapes consistant à: - charger les mesures (13i,imes) des magnétomètres (Mi,i), à un instant courant (t); - calculer une moyenne pondérée ( ) desdites mesures (13i,imes) des magnétomètres (Mi,j) représentant un champ magnétique terrestre moyen et un champ magnétique moyen dû au(x) objet(s) magnétique(s) mobile(s) (OMMk); - soustraire ladite moyenne pondérée ( ) de chacune desdites mesures (Bir") pour obtenir des mesures modifiées (13i,imes- ); - charger lesdites mesures modifiées (Biomes- ), ainsi qu'une localisation (Post, Mt) du ou des objets magnétiques mobiles (OMMk) à l'instant courant (t) en entrée d'un filtrage de localisation (FL) du ou des objets magnétiques mobile (OMMk); - effectuer le filtrage de localisation (FL) comprenant les étapes consistant à: - calculer une moyenne pondérée ( ) d'estimations (Bir) des données délivrées par lesdits magnétomètres (Mi); - soustraire ladite moyenne pondérée ( ) de chacune desdites données estimées (13i,jest); et - délivrer en sortie une localisation (Pasto, Mt+1) du ou des objets magnétiques mobile (OMMk) à un instant suivant (t+1). 30 2. Procédé selon la revendication 1, dans lequel ladite moyenne pondérée utilise des coefficients de pondération identiques pour chacun des N magnétomètres (Mi j). 3029642 17 3. Procédé selon la revendication 1, dans lequel ladite moyenne pondérée utilise des coefficients de pondération plus faibles pour P magnétomètres, parmi les N magnétomètres (Mi,j), les plus proches du ou des objets magnétiques mobiles (OMMk), P étant un nombre entier 5 inférieur à N. 4. Procédé selon la revendication 3, dans lequel ladite moyenne pondérée utilise des coefficients de pondération nuls pour lesdits P magnétomètres. 10 5. Procédé selon la revendication 1, dans lequel ladite moyenne pondérée utilise des coefficients de pondération nuls sauf pour le magnétomètre, parmi les N magnétomètres (Mi j), le plus éloigné du ou des objets magnétiques mobiles (OMMk). 15 6. Système de localisation d'au moins un objet magnétique mobile (OMMk) comprenant: - un réseau d'au moins N magnétomètres triaxiaux (ivi,d) reliés entre-eux mécaniquement sans aucun degré de liberté pour conserver une position 20 relative connue de ces magnétomètres (ivi,d), N étant un nombre entier au moins égal à 2; - une unité électronique de traitement (11) adapté pour déterminer la localisation (Post, Mt) du ou des objets magnétiques mobiles (OMMk) à partir de mesures (Bries) fournies par les magnétomètres (ivi,d) du réseau et 25 d'un filtre de localisation (FL); ladite unité électronique de traitement (11) comprenant: - des moyens de chargement des mesures (Bir") des magnétomètres (ivi,d), à un instant courant (t); - des moyens de calcul d'une moyenne pondérée ( ) desdites mesures (Bir") des magnétomètres (ivi,d) représentant un champ magnétique terrestre moyen et un champ magnétique moyen dû au(x) objet(s) magnétique(s) mobile(s) (OMMk); et - des moyens de calcul d'une soustraction de ladite moyenne pondérée ( ) de chacune desdites mesures (Bir") pour obtenir des mesures modifiées (13i,im"- ); le filtre de localisation (FL) comprenant : 3029642 18 - des moyens de chargement desdites mesures modifiées (13i,r"- ), ainsi que d'une localisation (Post, Mt) du ou des objets magnétiques mobiles (OMMk) à l'instant courant (t) en entrée d'un filtrage de localisation (FL) de l'objet magnétique mobile (OMM); 5 - des moyens de calcul d'une moyenne pondérée ( ) d'estimations (Brest) des données délivrées par lesdits magnétomètres (vi,d); - des moyens de calcul d'une soustraction de ladite moyenne pondérée ( ) de chacune desdites données estimées (13i,i"t); et - des moyens de délivrance en sortie d'une localisation (Pasto, Mt+1) du 10 ou des objets magnétiques mobile (OMMk) à un instant suivant (t+1). 15