Brevet FR3069634A1 — publié le 20190201

PROCEDE ET DISPOSITIF DE CARACTERISATION D'UN CAP DETERMINE A PARTIR DE LA MESURE DU CHAMP MAGNETIQUE DOMAINE TECHNIQUE GÉNÉRAL ET ART ANTÉRIEUR La présente invention concerne de façon générale les techniques 5 magnéto-inertielles. Plus précisément, elle concerne la détermination de cap par des magnétomètres. Elle trouve en particulier avantageusement application dans le cas de mesures en milieu urbain ou « indoor », c'est-à-dire à l'intérieur de 10 bâtiments. On utilise classiquement les magnétomètres pour le calcul du cap dans des systèmes embarqués. Dans ce cas, on fait l'hypothèse que le champ magnétique mesuré 15 par le capteur est le champ magnétique terrestre qui pointe, pour sa composante horizontale, vers le nord magnétique. La différence entre la direction du nord magnétique et du nord géographique (appelée déclinaison magnétique) est connue et tabulée. Sans perte de généralités, on considèrera donc dans la suite que le nord magnétique et le nord 20 géographique sont confondus, et que les magnétomètres indiquent donc ce qu'on appellera le nord. Typiquement, les méthodes de détermination utilisées pour le calcul 25 du cap à partir des mesures de magnétomètres s'appuient sur : - une modélisation qui relie la mesure fournie par le magnétomètre et une information relative au cap du système, - une caractérisation de la pertinence de cette relation pour les mesures réalisées. 30 Dans une approche classique, le calcul mis en oeuvre peut être un filtrage de type Kalman avec le champ magnétique comme mesure, 2 permettant de recaler le cap qui est contenu dans l'état, avec un bruit de mesure fonction de la caractérisation de la pertinence. La modélisation consiste à écrire que le champ magnétique mesuré M 5 contient une information de cap, par exemple par l'équation M - R op, 0, (p)TmTERRE où R est la matrice de rotation permettant de passer du référentiel de l'objet vers le référentiel inertiel terrestre, lp,O,(f) les angles d'Euler et MTERRE le champ magnétique terrestre. 10 La pertinence de cette équation est caractérisée par une variance de mesure, c'est-à-dire que l'on suppose que l'erreur de cette égalité est une variable aléatoire gaussienne d'espérance nulle. Cette variance sert à calculer automatiquement le gain de Kalman qui pondère le recalage pour tenir compte des différents bruits (bruits 15 dynamiques liés à l'environnement externe et bruits de mesure des magnétomètres). Dans une autre approche encore, on met en oeuvre un calcul par filtrage linéaire utilisant le même type de modélisation et de caractérisation. Dans cette approche, c'est le réglage relatif de l'amplitude du gain de 20 l'équation qui en caractérise la pertinence. En général, ce gain est réglé à la main. L'homme de l'art saura alors le pondérer en fonction de la pertinence de la modélisation. Ainsi, dans ces deux approches, les paramètres caractérisant la pertinence du modèle (variance du bruit gaussien dans le cas du filtre de 25 Kalman, gain dans le cas du filtrage linéaire) sont généralement des paramètres constants, indépendants du champ magnétique mesuré. Il a récemment été proposé, par exemple dans les publications W. T. Faulkner, R. Alwood, W. T. David, and J. Bohlin, "Gps-denied pedestrian tracking in indoor environments using an imu and 30 magnetic compass," in Proceedings of the 2010 International Technical Meeting of The Institute of Navigation, (San Diego, CA), pp. 198 - 204, January 2010. 3 - M. H. Afzal, V. Renaudin, and G. Lachapelle, "Magnetic field based heading estimation for pedestrian navigation environments," in 2011 International Conference on Indoor Positioning and Indoor Navigation, (Guimaraes, Portugal), Sep 2011. des caractérisations de la pertinence de la modélisation prenant en compte le champ magnétique mesuré. Les caractérisations ainsi réalisées ne donnent toutefois pas pleinement satisfaction. C'est le cas en particulier dans des environnements où le champ magnétique est fortement perturbé, notamment en « indoor » ou en milieu urbain, où des perturbations magnétiques issues des courants électriques ou des éléments métalliques dans les bâtiments se superposent au champ magnétique terrestre. PRÉSENTATION GÉNÉRALE DE L'INVENTION Un but général de l'invention est de proposer une solution permettant une meilleure caractérisation de la pertinence de la modélisation utilisée, notamment dans les environnements fortement perturbés. Notamment, l'invention propose un procédé de détermination de cap, dans lequel un champ magnétique est mesuré par au moins un magnétomètre et un cap magnétique est déterminé à partir de ladite mesure et dans lequel des moyens de calcul mettent en oeuvre un traitement calculant, pour chaque mesure de cap magnétique, une valeur de caractérisation du cap ainsi déterminé, caractérisé en ce que ladite valeur de caractérisation est fonction du gradient du champ magnétique mesuré, ladite valeur étant utilisée dans la détermination du cap. Les moyens de calcul déterminent pour un instant d'échantillonnage donné : - une estimation d'une prédiction de cap qui est fonction du cap déterminé à l'instant d'échantillonnage précédent, 4 - un recalage du cap prédit ainsi estimé en fonction d'un cap magnétique déterminé à partir de mesures de magnétomètres, et dans lequel le recalage est déterminé en fonction de ladite valeur de caractérisation. 5 Ledit recalage est avantageusement fonction du gradient du champ magnétique mesuré. De cette façon, on prend en compte le fait que, plus l'environnement est perturbé, plus le gradient est important et donc plus la mesure de cap est faussée. Dans un mode de mise en oeuvre préférentiel, les moyens de calcul 10 estiment l'amplitude du recalage en fonction : - d'un modèle d'évolution des perturbations magnétiques entre deux instants d'échantillonnage et - d'une estimation a priori de l'amplitude des perturbations. 15 De cette façon, le traitement prend en compte l'évolution de la perturbation de cap magnétique entre deux instants d'échantillonnage et tient compte de l'éventuelle corrélation spatiale de la perturbation due à l'environnement dans le cas d'instants d'échantillonnage successifs. 20 Le calcul du qui en résulte est plus fiable que dans l'art antérieur. Notamment, pour un instant d'échantillonnage donné k +1, les moyens de calcul estiment une perturbation de cap liée en calculant tp(d)[k + 1] = a[k]ap(d)[k] + û[k]) + (1 - a[k]){:c] + vP(a) [k] 2 ' 25 Où a[k] = 1 2a[k]2 et vP (d) est une variable aléatoire gaussienne de variance I a[k]2(1 - a[k]2) 5 - k étant l'instant d'échantillonnage précédent, - a[k] étant un paramètre qui représente l'amplitude a priori de la perturbation magnétique, - ii et citiétant deux paramètres calculés par les moyens de 5 calculs comme une estimation de l'espérance et de la variance de l'évolution de la perturbation. De façon avantageuse, les moyens de calcul estiment le paramètre a[k] comme une fonction linéaire de la norme du gradient du champ magnétique. 10 Le paramètre ii[k] est quant à lui par exemple calculé par les moyens de calcul comme la différence entre les caps magnétiques déterminés directement à partir des sorties magnétomètres pour les instants k+1 et k, à laquelle est soustraite l'évolution de cap prédite (co[k]dt). Le paramètre au peut être estimé comme fonction de la vitesse de 15 déplacement ou du déplacement entre deux instants d'échantillonnage successifs. Le traitement met par exemple en oeuvre un filtrage Kalman dont l'état a au moins comme paramètre le cap réel et la perturbation de cap magnétique (1p, ti)(d)) 20 La prédiction de cap peut être déterminée en fonction de mesures d'un ou plusieurs capteurs d'une centrale inertielle. L'invention concerne en outre un dispositif détermination de cap par des capteurs magnétiques, comportant des magnétomètres et des moyens de calcul pour le calcul du cap à partir du champ magnétique 25 mesuré par lesdits magnétomètres, caractérisé en les moyens de calcul mettent en oeuvre, pour différents instants d'échantillonnage successifs, le traitement précité. Elle propose également un système de navigation magnéto-inertielles comportant au moins un tel dispositif de mesure de cap. 30 On utilise avantageusement un tel système en environnement urbain ou à l'intérieur de bâtiments. L'invention concerne également : 6 - un produit programme d'ordinateur comprenant des instructions de code pour l'exécution d'un procédé du type précité, lorsque ledit programme est exécuté sur un ordinateur ; 5 - 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 3 illustre les étapes principales d'un procédé conforme à un mode de mise en oeuvre de l'invention. 20 DESCRIPTION D'UN OU PLUSIEURS MODES DE MISE EN OEUVRE ET DE RÉALISATION Généralités - Dispositif de mesure 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 dans 25 un champ magnétique ambiant (typiquement le champ magnétique 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. 7 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 5 capteurs de mesure magnétique 20, i.e. des magnétomètres axiaux 20. 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 à au niveau dudit magnétomètre 20 selon son axe. 10 Plus précisément, les magnétomètres 20 sont solidaires du boitier 2. 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 15 position prédéterminée dans ce repère. 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 autoagrippante qui enserre le membre et permet la liaison solidaire. 20 Bien entendu, l'invention n'est pas limitée à l'estimation du 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. 25 Le boitier 2 peut comprendre des moyens de calcul 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 qu'un 30 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 serveur 4. 8 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 UMTS/LTE) 5 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. Si c'est un terminal mobile 3 (respectivement un serveur 4) qui 10 héberge « l'intelligence », il comprend des moyens de calcul 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 calcul utilisés sont ceux 21 du boitier 2, celui-ci peut encore inclure des moyens de communication 25 pour transmettre la position estimée. Par 15 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 calcul de données 21, 31, 41 respectivement du 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é. 20 Ils comportent à cet effet chacun des moyens de stockage dans lesquels sont mémorisés tout ou partie de séquences d'instructions de code pour l'exécution du procédé. Prédiction et recalage 25 Les moyens de calcul mettent en oeuvre (figure 3) un filtrage 100 calculant d'une part une estimation de la valeur de cap par prédiction (étape 101) et d'autre part mettant en oeuvre un recalage en fonction d'une estimation d'erreur (étape 102). Notamment, connaissant la vitesse angulaire co donnée par exemple 30 par un gyromètre, l'étape 101 calcule le cap pk+1 à un instant k+1 comme égal à Ok+1 = 11/k + W. At 9 où Ipk est le cap à l'instant k précédent et où A t est la durée séparant ces deux instants d'échantillonnage. Le recalage 102 prend en compte les mesures réalisées par les magnétomètres 20. 5 Dans ce qui suit, la mesure du cap magnétique (issue des mesures des magnétomètres) pour un instant k donné est notée z,p[k], avec z[k] = Typiquement, le cap magnétique est donné par z,p[k] = atan (111Y[kl) Mx[k] 10 où My et Mx sont deux composantes horizontales du champ magnétique dans un référentiel terrestre, ces deux composantes étant calculées en fonction de l'attitude déterminée pour l'objet 1 par une centrale inertielle de celui-ci. Les moyens de calcul calculent le cap recalé Ihrze+cia_ le correspondant à 15 l'instant k+1 comme égal à la somme du cap estimé pk+1 à cet instant et d'un recalage qui est avantageusement fonction - d'un gain Kk calculé ou ajusté par rapport à l'instant k précédent - d'une estimation de l'erreur Err(1pk+1,z,p) entre le cap prédit pk+1 et la mesure 20 Ainsi, les moyens de calcul utilisent le cap magnétique issu de la mesure par les magnétomètres pour recaler l'état, et en particulier le cap, en calculant ,recale_ +1 V'k+1 Kk- E rr(V) k+ z,p) Typiquement, l'erreur Err(1pk+1,z,p) peut être une simple différence 25 entre le cap prédit, estimé à l'instant k+ 1 (avant recalage), et le cap magnétique z, issu des mesures des magnétomètres. D'autres fonctions d'erreurs sont néanmoins possibles, notamment dans le cas d'un filtrage non linéaire. En particulier, le recalage Kk. Err(IPk+i, Zip) et plus particulièrement le 30 gain Kk, peut avantageusement dépendre du gradient du champ magnétique mesuré par les magnétomètres 20. 10 De cette façon, la correction de recalage est d'autant plus importante que le champ magnétique varie fortement et est donc susceptible d'induire des erreurs importantes sur la mesure de cap. La valeur de cap recalée ainsi obtenue est stockée par lesdits moyens 5 de calcul 21, 31, 41 et/ou utilisée par lesdits moyens pour la suite du traitement et pour les calculs d'informations de navigation magnéto-inertielles (vitesse linéaire, vitesse angulaire, position, cap, etc...). Egalement, elle peut être transmise par les moyens de calcul à des moyens d'interface, par exemple sur le téléphone pour être affichée sur 10 l'écran de celui-ci. Modèle d'évolution des perturbations liées à l'environnement On détaille dans ce qui suit un exemple de calcul possible du recalage dans lequel la mesure de cap magnétique z,l, est corrigée d'une estimation 15 de la perturbation du champ magnétique liée à l'environnement. Le cap magnétique z,l, issu de la mesure par les magnétomètres peut en effet être considéré comme se décomposant comme suit : = lp + 1p(d) + vzw 20 OÙ: -1p correspond au cap réel que l'on cherche à déterminer, - v correspond à une erreur gaussienne de mesure et - 1p(d) (« d » pour « disturbance » en anglais) correspond à une perturbation de cap magnétique lié à l'environnement (typiquement, 25 perturbations liées aux infrastructures métalliques et aux câbles électriques dans des environnements urbains ou dans des bâtiments). La perturbation de cap magnétique liée à l'environnement est fortement corrélée spatialement. Le champ magnétique lié aux perturbations de l'environnement est en effet un champ de vecteur continu 30 et les champs magnétiques en deux points A et B donnés sont d'autant plus proches que ces deux points A et B sont voisins dans l'espace. 11 Dans un mode de mise en oeuvre, la perturbation de cap 1p(d) lié à l'environnement est estimée les moyens de calcul au moyen d'une formulation prenant en compte les corrélations d'évolution que l'on peut attendre pour les perturbations de champ magnétique entre deux instants 5 d'échantillonnage successifs (corrélation temporelle, spatiale ou plus complexe). Cette estimation est construite pour permettre un modèle markovien (qui peut être utilisé dans un filtre récursif) pour l'évolution de Ipd et - permet de prendre en compte un modèle de l'évolution de la 10 perturbation de cap magnétique (corrélation temporelle, spatiale ou plus complexe) - permet de construire un modèle de filtrage pour lequel le cap est observable (sinon, on perd tout espoir de construire un estimateur de cap) 15 Dans l'estimation proposée, l'amplitude du recalage est estimée en fonction d'un modèle d'évolution des perturbations magnétiques entre deux instants d'échantillonnage et d'une estimation a priori de l'amplitude des perturbations. Les inventeurs ont trouvé (et vérifié mathématiquement) que dans le 20 cas d'une corrélation spatiale (cas de perturbations en milieu urbain ou « indoor »), une estimation adaptée de l'évolution de la perturbation 1p(d) entre deux instants d'échantillonnage séparés par un pas de temps avantageusement comme suit : tp(d)[k + 1] = a[k]ap(cO[k] + û[k]) + (1 - a[k])û[k] + v -2 tp(d)[Ic] 25 Où a[k] = 1 2a[k]2 et vP (d) est une variable aléatoire gaussienne de variance I a[k]2(1 - a[k]2) au[k] 2 12 - a[k] étant un paramètre qui représente l'amplitude a priori de la perturbation magnétique, - ii et ciu étant l'espérance et la variance d'une variable aléatoire qui est un modèle (sous la forme d'une variable aléatoire gaussienne) 5 de la variation de la perturbation de cap magnétique entre 2 instants d'échantillonnage. est une estimation prenant en compte le bruit des magnétomètres. Cette estimation est retranchée de z,l, pour le calcul d'erreur, celle-ci pouvant être directement calculée entre pk+1 et (z,p - tp (d) [k + 1]) , ou confiée 10 à un filtre de type Kalman qui a à la fois 1p(d) et lp dans son état. Détermination du paramètre a[k] Le paramètre a[k] représente l'amplitude a priori de la perturbation du cap magnétique. Il caractérise le cap et est par exemple calculé comme 15 une fonction linéaire de la norme du gradient du champ magnétique. Par exemple, a[k] = cto + aiN [k] avec N[k] = r3[1c111 20 cto et alétant deux paramètres préalablement fixés avant mise en oeuvre du traitement de filtrage, fi [k] tant le champ magnétique à l'instant k. De cette façon, comme indiqué plus haut, le terme de recalage Err(1pk+1,z,p) est fonction du gradient du champ magnétique mesuré par les 25 magnétomètres 20, ce qui permet une correction efficace. D'autres méthodes sont possibles pour déterminer le paramètre a[k] . Notamment, a[k] peut également être déterminé par comparaison avec un modèle du champ magnétique terrestre, par exemple en mettant en oeuvre les techniques telles que proposées dans la publication : 30 "Assessment of Indoor Magnetic Field Anomalies using Multiple Magnetometers Assessment of Indoor" - M.H. Afzal, V. Renaudin, G. Lachapelle - Proceedings of the 23rd International Technical Meeting of The 13 Satellite Division of the Institute of Navigation (ION GNSS 2010) September 21 - 24, 2010. Les solutions dans lesquelles le paramètre a[k] est fonction du gradient du champ magnétique permettent toutefois de meilleures 5 corrections, et ne nécessitent pas de modèles du champ magnétique terrestre. Détermination du paramètre ii Le terme ved) peut être considéré en première approximation comme 10 négligeable et être considéré comme nul dans le calcul. Le paramètre i2[k] représente la valeur la plus probable pour l'évolution de la perturbation. Ce paramètre est_par exemple calculé en utilisant les gyromètres qu'on suppose bons à court terme. L'évolution la plus probable de la 15 perturbation est alors donnée par la différence entre l'évolution de la mesure de cap magnétique et de l'évolution du cap gyrométrique : fi[k] = E [tp(d)[k + 1] - tp(d)[k]] = E [(d) [k +11+ tp[k + 1] - (te (d)[k] + tp[k]) - (tp[k +1] - tp[kpl cap magnétique mesuréà l'instantk+1 cap magnétique mesuré àl'instantk évolution de cap = Ztp[k + 1] - z»[k] co[k]dt où : E[.] est l'espérance mathématique, w[k] est la vitesse de rotation gyrométrique et dt est le pas d'échantillonnage. 20 Ainsi, le paramètre ii[k] est calculé par les moyens de calcul comme la différence entre les caps magnétiques déterminés directement à partir des sorties magnétomètres pour les instants k+1 et k, à laquelle est soustraite la rotation prédite co[k]dt. 25 Détermination du paramètre 0-' Le paramètre au représente l'image de la corrélation entre la perturbation au pas k et celle au pas k+1. On peut choisir avantageusement de l'indexer sur la vitesse de 30 déplacement (ou sur le déplacement entre deux instants d'échantillonnage 14 successifs, ce qui est similaire) de manière à rendre compte de la corrélation spatiale des perturbations magnétiques. On pourra par exemple prendre = c. 11v11 où c est un coefficient de réglage. 5 Autre mode de mise en oeuvre Dans un autre mode de mise en oeuvre, on combine différentes valeurs de cap magnétique, afin d'obtenir une valeur de variable dépendant du cap la moins possible entachée d'erreur liées aux perturbations 10 magnétiques. La variable que l'on cherche à déterminer est par exemple le cap moyen d'un morceau de trajectoire, d'autres variables dépendant du cap pouvant bien entendu être envisagées. Pour cela, les moyens de calcul enregistrent, pour chaque pas 15 d'échantillonnage k, la valeur du cap magnétique mesuré z,/,[k] et la valeur de la caractérisation du cap magnétique appelée a[k], calculée comme une fonction du gradient. Par exemple, a[k] = cto + G/111'7BM Le cap17;ftk0k1le plus probable pour le morceau de trajectoire défini 20 entre deux instants d'échantillonnage k = /cc) et k = /cl peut ainsi être posé sous la forme d'une optimisation - argmin - {/j112 17) k où II. II est une norme (L1, L2 ou autre). L'algorithme d'optimisation mis en oeuvre par les moyens de calcul 25 peut être de tout type adapté : descente de gradient, algorithme de Levenberg-Marquardt, ou algorithme de Metropolis-Hastings. Par ailleurs, dans le cas de la norme L2, le problème est alors posé sous la forme d'un moindre carré et les moyens de calcul déterminent directement vk ztp[k] I L'ko a[k] [[ko,k 111 = klo (1 / a[k]) 30 15 On comprendra bien que cette valeur est beaucoup plus avantageuse que la valeur moyenne non pondérée par la valeur de caractérisation puisque les instants k où la perturbation magnétique est la plus forte ont beaucoup moins d'impact REVENDICATIONS1. Procédé de détermination de cap, dans lequel un champ magnétique est mesuré par au moins un magnétomètre et un cap magnétique est déterminé à partir de ladite mesure et dans lequel des moyens de calcul mettent en oeuvre un traitement calculant, pour chaque mesure de cap magnétique, une valeur de caractérisation du cap ainsi déterminé, caractérisé en ce que ladite valeur de caractérisation est fonction du gradient du champ magnétique mesuré, ladite valeur étant utilisée dans la détermination du cap. 2. Procédé selon la revendication 1, dans lequel les moyens de calcul déterminent pour un instant d'échantillonnage donné - une estimation d'une prédiction de cap qui est fonction du cap déterminé à l'instant d'échantillonnage précédent, - un recalage du cap prédit ainsi estimé en fonction d'un cap magnétique déterminé à partir de mesures de magnétomètres, et dans lequel le recalage est déterminé en fonction de ladite valeur de caractérisation. 3. Procédé selon la revendication 1, dans lequel les moyens de calcul déterminent une valeur d'une variable dépendant du cap en fonction d'une pluralité de caps magnétiques mesurés et de valeurs de caractérisation. 4. Procédé selon la revendication 1, dans lequel, lors dudit recalage, les moyens de calcul estiment l'amplitude du recalage en fonction : - d'un modèle d'évolution des perturbations magnétiques entre deux instants d'échantillonnage et - d'une estimation a priori de l'amplitude des perturbations. 17 5. Procédé selon la revendication 2, dans lequel, pour un instant d'échantillonnage donné k+1, les moyens de calcul estiment une perturbation de cap liée en calculant a[k] (tp (d) [k] + û[k]) + (1- a[k]) û 21c1 + v ed)[k] 5 Où au[k] 2 a[k] = 1 2a[k]2 O ù Vtp(ct) est une variable aléatoire gaussienne de variance a[k]2(1 - a[k]2). - k étant l'instant d'échantillonnage précédent, 10 - a[k] étant un paramètre qui représente l'amplitude a priori de la perturbation magnétique û et atiétant deux paramètres calculés par les moyens de calculs comme une estimation de l'espérance et de la variance de l'évolution de la perturbation. 15 6. Procédé selon la revendication 5, dans lequel les moyens de calcul estiment le paramètre a[k] comme une fonction linéaire de la norme du gradient du champ magnétique. 7. Procédé selon la revendication 5, dans lequel le paramètre ii[k] est 20 calculé par les moyens de calcul comme la différence entre les caps magnétiques déterminés directement à partir des sorties magnétomètres pour les instants k+1 et k, à laquelle est soustraite l'évolution de cap prédite (co[k]dt). 25 8. Procédé selon la revendication 5, dans lequel le paramètre au est estimé comme fonction de la vitesse de déplacement ou du déplacement entre deux instants d'échantillonnage successifs. 9. Procédé selon la revendication 5, dans lequel le traitement met en 30 oeuvre un filtrage Kalman dont l'état a au moins comme paramètre le cap réel et la perturbation de cap magnétique (tp, Ip(d)). 18 10. Procédé selon l'une des revendications précédentes, dans lequel la prédiction de cap est déterminée en fonction de mesures d'un ou plusieurs capteurs d'une centrale inertielle. 5 11. Dispositif de détermination de cap par des capteurs magnétiques, comportant des magnétomètres et des moyens de calcul pour le calcul du cap à partir du champ magnétique mesuré par lesdits magnétomètres, caractérisé en les moyens de calcul mettent en oeuvre, 10 pour différents instants d'échantillonnage successifs, le traitement du procédé selon l'une des revendications précédentes. 12. Système de navigation magnéto-inertielles comportant au moins un dispositif de mesure de champ magnétique selon la revendication 11. 15 13. Utilisation d'un système selon la revendication 12 dans des applications en environnement urbain ou à l'intérieur de bâtiments. 14. Produit programme d'ordinateur comprenant des instructions de code 20 pour l'exécution d'un procédé de mesure de champ magnétique selon l'une des revendications 1 à 10, lorsque ledit programme est exécuté sur un ordinateur. 15. Moyen de stockage lisible par un équipement informatique sur lequel 25 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 à 10. 30