Brevet FR3082612A1 — publié le 20191220

Procédé de calibration d'un gyromètre équipant un objet DOMAINE TECHNIQUE GENERAL La présente invention concerne le domaine de la navigation sans GNSS. Plus précisément, elle concerne un procédé de calibration d'un gyromètre solidaire de magnétomètres. 10 ETAT DE L'ART Il est aujourd'hui commun de suivre la position d'un véhicule par GNSS (Global Navigation Satellite System, par exemple le GPS) ou en utilisant un réseau de communication sans fil (triangulation à l'aide de 15 bornes émettrices, réseau wifi ou autres). Ces méthodes s'avèrent très limitées car elles ne permettent pas d'assurer disponibilité et précision de l'information, l'une et l'autre étant affectées par d'éventuels masquages entre les sources et le récepteur. Elles s'avèrent également dépendantes de technologies extérieures comme 20 les satellites pour le GNSS qui peuvent être indisponibles voire volontairement brouillées. Alternativement, on connait aussi des méthodes « autonomes » pour suivre dans n'importe quel environnement le déplacement relatif d'un véhicule grâce par exemple à une centrale inertielle ou magnéto-inertielle. 25 Par déplacement relatif, on entend la trajectoire du véhicule dans l'espace par rapport à un point et à un repère donnés à l'initialisation. En plus de la trajectoire, ces méthodes permettent également d'obtenir l'orientation du véhicule par rapport au même repère initial. Une centrale inertielle de classe navigation est constituée au 30 minimum de trois accéléromètres et de trois gyromètres disposés en triaxe. Typiquement, les gyromètres « maintiennent » un repère, dans lequel une 2 double intégration temporelle des mesures des accéléromètres permet d'estimer le mouvement. Il est notablement connu que pour pouvoir utiliser les méthodes de navigation inertielle classique, telles que mises en oeuvre dans les applications lourdes comme la navigation des avions de chasse ou de ligne, des sous-marins, des navires, etc., il est nécessaire d'utiliser des capteurs de très haute précision. En effet la double intégration temporelle d'une mesure d'accélération fait qu'une erreur constante d'accélération crée une erreur de position qui augmente de façon proportionnelle au carré du temps. Alternativement, il est connu que si l'information de vecteur vitesse dans le repère du porteur peut être fournie par une source externe (odométrie sur une voiture, loch sur un navire, tube de Pitot sur un avion), il est possible d'obtenir la trajectoire du porteur dès lors que l'on est en capacité de connaître son orientation et notamment son cap. Dans un certain nombre de cas, le cap initial est connu (« alignement » initial d'un dispositif inertiel) et peut être obtenu par des capteurs autres qu'inertiels (par exemple magnétiques) ou déduit directement des capteurs inertiels lorsque ceux-ci permettent de mesurer précisément la rotation terrestre. La difficulté si l'on souhaite travailler avec des capteurs bas coût est d'« augmenter » leur performance en calibrant en continu leurs paramètres d'erreur, de sorte à maintenir au mieux l'information d'orientation d'un porteur, en particulier l'information de cap sur un véhicule terrestre. Il a par conséquent été proposé dans ce contexte de calibrer un 25 gyromètre grâce à des mesures magnétiques. La demande US2012/0116716 propose de regarder si l'évolution d'un champ magnétique entre deux points de mesure correspond au changement de l'attitude prédite par l'intégration des gyromètres. Elle fait cependant l'hypothèse que les mesures 30 magnétométriques/accelérométriques sont fiables, ce qui est délicat. En effet, on constate en pratique que : 3 le champ magnétique ambiant n'est pas toujours homogène et stationnaire ; les éléments métalliques (par exemple une carrosserie) ou aimants à proximité ont un impact sur la mesure magnétique 5 (effets de type fer doux et fer dur), voir également la demande US2004/0123474. Pour calibrer correctement un magnétomètre, il est nécessaire de le mettre dans un maximum d'orientation spatiale de manière à « exciter le modèle d'erreur » pour avoir de l'observabilité. Ce 10 problème est facilement soluble pour un équipement qu'on peut porter à la main et qu'on peut ainsi tourner dans tous les sens, voir par exemple la demande FR1757082, avant de s'en servir. Si au contraire le magnétomètre équipe un véhicule, on est contraint : 15 o par le fait que lors d'un trajet typique d'un véhicule, roulis et tangage ne dévient de zéro que faiblement. o par le fait que la caisse s'aimante, ou bien qu'on a chargé le véhicule, ce qui fait que la calibration peut n'être plus pertinente; ce qui va nécessiter des recalibrations 20 régulières. La demande US2011/0178707 concerne plus précisément les dispositifs de type smartphone comprenant une boussole magnétique et un gyromètre, et évoque le problème des perturbations magnétiques. Cette demande vise plutôt la calibration de la boussole magnétique à partir des données gyrométriques, mais propose également l'inverse (i.e. la calibration d'un gyromètre à partir de données magnétiques), en prenant quant à elle comme point de départ un quaternion q magnéto-accélérométrique et en obtenant la vitesse angulaire de recalage à partir de ce quaternion o=2q-1 q. Cependant, il est toujours fait l'hypothèse que les mesures magnétométriques/accelérométriques sont fiables, et dans le cas contraire il est juste proposé soit de ne pas mettre à jour le biais du gyromètre tant que 4 cette situation dure, soit d'utiliser un quaternion par défaut stocké dans une mémoire. Il serait souhaitable de disposer d'une nouvelle méthode de calibration d'un gyromètre d'un objet en vue de l'estimation du mouvement 5 de cet objet qui permette une excellente qualité de résultat et ne soit pas contraignante. PRESENTATION DE L'INVENTION 10 La présente invention se rapporte ainsi selon un premier aspect à un procédé de calibration d'un gyromètre équipant un objet évoluant dans un champ magnétique ambiant, le procédé étant caractérisé en ce qu'il comprend des étapes de : (a) Acquisition 15 - par le gyromètre, d'une vitesse angulaire mesurée de l'objet, et - par des moyens de mesure magnétique solidaires dudit objet, d'au moins deux composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique, au niveau des moyens de mesure magnétique ; 20 (b) Détermination par des moyens de traitement de données de valeurs d'au moins un paramètre de calibration du gyromètre minimisant une première expression définie par une vitesse angulaire estimée de l'objet et au moins une première équation magnétique sur les composantes du champ magnétique et/ou d'une dérivée i-ième du 25 champ magnétique, - la vitesse angulaire estimée de l'objet étant fonction de la vitesse angulaire mesurée et de paramètres de calibration du gyromètre, et - l'au moins une première équation magnétique supposant un 30 caractère uniforme et stationnaire du champ magnétique au niveau des moyens de mesure magnétique. 5 Selon d'autres caractéristiques avantageuses et non limitatives : - la vitesse angulaire estimée co.(g;ration) de l'objet est liée à la vitesse angulaire mesurée Cil(mesure) par un modèle co(estimat(on) = D - gyro syro (w9Yro(mesure) bgyro), où D et b gyro sont les paramètres de calibration du 5 gyromètre ; - le procédé comprend en outre une étape (c) d'estimation d'un paramètre représentatif d'une erreur sur les paramètres de calibration ; - le procédé comprend en outre, si ledit paramètre représentatif d'une erreur est supérieur à un seuil prédéterminé, une étape (d) d'acquisition par 10 des moyens d'acquisition supplémentaires d'une vitesse linéaire mesurée de l'objet, et de re-détermination par les moyens de traitement de données de valeurs de l'au moins un paramètre de calibration du gyromètre minimisant une deuxième expression définie par la vitesse angulaire estimée de l'objet et au moins une deuxième équation magnétique sur la 15 vitesse linéaire mesurée et les composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique, l'au moins une deuxième équation magnétique supposant seulement un caractère stationnaire du champ magnétique au niveau des moyens de mesure magnétique ; - le procédé comprend une étape (e) de nouvelle estimation dudit 20 paramètre représentatif d'une erreur sur les paramètres de calibration ; - ledit paramètre représentatif d'une erreur sur les paramètres de calibration est fonction de la valeur de ladite première expression, ou le cas échéant de ladite deuxième expression, calculée pour les valeurs déterminées des paramètres de calibration ; 25 - l'étape (c), et le cas échéant l'étape (e), comprend la détermination d'un intervalle d'erreur accumulée en fonction de la valeur dudit paramètre représentatif d'une erreur sur les paramètres de calibration et d'une durée depuis la dernière calibration du gyromètre, ledit seuil prédéterminé dépendant dudit intervalle d'erreur accumulée ; 30 - l'étape (b), et le cas échéant l'étape (d), comprend la mise en oeuvre d'un filtre récursif ou d'une optimisation ; 6 - les composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique acquises à l'étape (a) sont des composantes mesurées, et celles utilisées par lesdites première et deuxième équations magnétiques sont des composantes estimées en fonction des composantes 5 mesurées et de paramètres de calibration des moyens de mesure magnétique, l'étape (b) et, le cas échéant l'étape (d), comprenant également la détermination de valeurs d'au moins un paramètre de calibration des moyens de mesure magnétique ; - les moyens de mesure magnétique sont des magnétomètres triaxes, 10 les composantes estimées fi eestimati°n) du champ magnétique étant liées aux composantes mesurées M(mesure) par un modèle emesur e) = A - m (estimation) + /,,,i 1/2migneto, où A et bmaneto sont les paramètres de calibration des moyens de mesure magnétique ; - la ou les premières équations magnétiques sont de la forme Vn"M = 15 f.,(co) x V,M, et le cas échéant la ou les deuxièmes équations magnétiques sont de la forme V1114 fri.(co) x VnAl + 9n(Vn+1.114).V, où V1M est une dérivée n-ième du champ magnétique, cd la vitesse angulaire, V la vitesse linéaire, et !net gt, des fonctions prédéterminées ; - la première équation magnétique est M = -w x M, et le cas échéant 20 la deuxième équation magnétique est M = -CO xM+VM-V; - la première expression est fonction de lynm (estimatton) fn(wa(eysrtoimation)) x vnm (estimation)1 et l'éventuelle deuxième expression est fonction de 1,9.11111(es/in/anon) _ ( (estimation)) in CO 9Yro vniestimution) gn(vn±im (estimation)). 25 - ledit paramètre représentatif d'une erreur sur les paramètres de calibration est soit la moyenne de Ivia(esiimation) _ f Wgyro n (estimation)) V 2 ou x vnm(estimation) 1Vnfieesitmatton) fn(wg(eysrtoimation)) (estimation)) 12 v sur un intervalle de temps donné, soit tout ou 9n(V11+1M 7 partie d'un gradient spatial desdites composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique ; - le procédé comprend une étape (f) d'estimation par les moyens de traitement de données du mouvement dudit objet (1) en fonction de la 5 vitesse angulaire mesurée de l'objet, des composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique, et/ou de l'éventuelle vitesse linéaire mesurée de l'objet, et des valeurs des paramètres de calibration, dans lequel une orientation du véhicule est estimée à l'étape (f) uniquement en fonction de la vitesse angulaire 10 mesurée du véhicule et des valeurs des paramètres de calibration. Selon un deuxième aspect, est proposé un objet évoluant dans un champ magnétique ambiant, comprenant un gyromètre configuré pour acquérir une vitesse angulaire mesurée de l'objet, des moyens de mesure 15 magnétique configurés pour acquérir au moins deux composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique, l'objet étant caractérisé en ce qu'il comprend en outre des moyens de traitement de données configurés pour : déterminer des valeurs d'au moins un paramètre de calibration du 20 gyromètre minimisant une première expression définie par une vitesse angulaire estimée de l'objet et au moins une première équation magnétique sur les composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique, o la vitesse angulaire estimée de l'objet étant fonction de la 25 vitesse angulaire mesurée et de paramètres de calibration du gyromètre, et o l'au moins une première équation magnétique supposant un caractère uniforme et stationnaire du champ magnétique au niveau des moyens de mesure magnétique. 30 Selon d'autres caractéristiques avantageuses et non limitatives : 8 - l'objet comprend en outre des moyens d'acquisition supplémentaires configurés pour acquérir une vitesse linéaire mesurée de l'objet, les moyens de traitement de données étant en outre configurés pour : Estimer un paramètre représentatif d'une erreur sur les paramètres 5 de calibration ; et Si ledit paramètre représentatif d'une erreur est supérieur à un seuil prédéterminé, re-déterminer des valeurs de l'au moins un paramètre de calibration du gyromètre une deuxième expression définie par la vitesse angulaire estimée de l'objet et au moins une deuxième 10 équation magnétique sur la vitesse linéaire mesurée et les composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique, l'au moins une deuxième équation magnétique supposant seulement un caractère stationnaire du champ magnétique au niveau des moyens de mesure magnétique ; 15 - L'objet est un véhicule à roues, les moyens d'acquisition supplémentaire étant au moins deux odomètres équipant des roues du véhicule. Selon un troisième et un quatrième aspect, sont proposés un produit 20 programme d'ordinateur comprenant des instructions de code pour l'exécution d'un procédé selon le premier aspect de l'invention de calibration d'un gyromètre ; et 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é selon le premier aspect 25 de calibration d'un gyromètre. PRESENTATION DES FIGURES D'autres caractéristiques et avantages de la présente invention 30 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 : 9 la figure 1 représente un exemple d'architecture de véhicule pour la mise en oeuvre du procédé selon l'invention ; la figure 2 est un diagramme représentant les étapes d'un mode de réalisation préféré du procédé selon l'invention. 5 DESCRIPTION DETAILLEE Architecture 10 En référence à la figure 1, le présent procédé permet la calibration d'un gyromètre 11 d'un objet 1 évoluant dans un champ magnétique ambiant (typiquement le champ magnétique terrestre, le cas échéant altéré par les objets métalliques à proximité), noté M. Comme déjà expliqué, le champ magnétique est un champ vectoriel en l'espèce tridimensionnel, 15 c'est-à-dire associant un vecteur de dimension trois à chaque point de tridimensionnel dans lequel l'objet 1 est mobile. Cet objet 1 peut être n'importe quel objet mobile dont la connaissance de la position est souhaitée, par exemple un véhicule, en particulier un véhicule à roues, un drone, etc., mais également une 20 personne ou une partie du corps de cette personne (ses mains, sa tête, etc.). L'objet 1 est équipé du gyromètre 11 et de moyens de mesure magnétique 20. Les moyens 20 et le gyromètre 11 sont solidaires de l'objet 1, par exemple solidaires de la carrosserie dans le cas d'un véhicule, et de 25 façon générale fixes dans le repère de l'objet 1. Les moyens de mesure magnétique 20 sont typiquement des magnétomètres et/ou des gradiomètres magnétiques. Dans ce dernier cas, au lieu de mesurer la valeur des composantes du champ magnétique M, les gradiomètres mesurent directement la valeur des composantes du gradient 30 du champ magnétique M, i.e. la valeur des dérivées spatiales (éventuellement en plus du champ magnétique). De tels gradiomètres magnétiques 20 sont connus de l'homme du métier. Selon certains modes 10 de réalisation de l'invention, on pourra utiliser des gradiomètres mesurant directement la valeur des dérivées secondes (gradient d'ordre deux), et de façon générale des dérivées i-ièmes (gradient d'ordre O. Dans la suite de la description on prendra l'exemple de 5 magnétomètres mais l'homme du métier saura transposer à des gradiomètres. Le ou les magnétomètres 20 sont préférentiellement « axiaux », c'est à dire capables de mesurer une composante dudit champ magnétique, i.e. la projection dudit vecteur champ magnétique M selon leur axe (ou, dans le 10 cas d'un gradiomètre, une composante d'une dérivée spatiale d'ordre i dudit champ magnétique, i.e. la variation d'une composante de Vi-1M au niveau du gradiomètre selon son axe). Les magnétomètres 20 sont au moins au nombre de deux, et préférentiellement au moins au nombre de trois, de sorte à pouvoir acquérir 15 respectivement au moins deux ou trois composantes du champ magnétique. On comprendra en effet que deux magnétomètres peuvent suffire, s'ils sont disposés orthogonalement l'un par rapport à l'autre dans un plan orthogonal à l'axe de rotation vertical du véhicule, i.e. un plan horizontal. Néanmoins, utiliser trois magnétomètres 20 et acquérir trois 20 composantes permet de s'affranchir de la nécessité de cette disposition particulière. Avantageusement les magnétomètres 20 sont même d'un nombre supérieur à trois, comme six voire même neuf, avantageusement organisés par groupes de trois en « triaxes », i.e. un triplet de magnétomètres 20 deux 25 à deux orthogonaux associés à une même position spatiale et mesurant le champ magnétique selon les trois axes (dans le cas de six magnétomètres 20 on a deux triaxes et dans le cas de neuf magnétomètres 20 on a trois triaxes). Deux triaxes permettent de mesurer le champ magnétique ainsi qu'une partie de son gradient (voir plus loin), et à partir de huit 30 magnétomètres 20 qui ont été convenablement positionnés et orientés, on peut déterminer toutes les composantes du gradient. Ceci est par exemple 11 le cas pour trois triaxes positionnés aux coins d'un triangle rectangulaire isocèle. De façon préférée, le repère orthonormé associé à l'objet est choisi par convention (et par facilité pour la suite de la présente description) tel 5 que les triaxes sont avantageusement orientés conformément audit repère orthonormé, de sorte à faciliter encore les calculs. Mais l'homme du métier saura dans tous les cas transposer à n'importe quelle disposition spatiale de magnétomètres/gradiomètres. 10 Par gyromètre 11, on entend des moyens de mesure inertielle aptes à mesurer la vitesse angulaire de l'objet 1 selon un système de trois axes orthogonaux, qui définissent le repère de l'objet 1, i.e. mesurer les trois composantes d'une vecteur vitesse angulaire noté w. On comprend ainsi que le gyromètre 11 peut en fait désigner un ensemble de trois gyromètres 15 associés à un des trois axes, en particulier en triaxe (i.e. chacun apte à mesure une des trois composantes du vecteur vitesse angulaire noté w). Dans le mode de réalisation préféré où l'objet 1 est un véhicule, ce repère comporte préférentiellement un axe vertical et deux axes horizontaux, en particulier un axe longitudinal et un axe transversal. La 20 rotation autour de l'axe vertical du véhicule est décrite par l'angle sur lequel le conducteur agit en tournant le volant. Sur un sol généralement plat, les changements de direction du véhicule sont dans le plan horizontal, i.e. également selon ledit axe vertical du lacet. En réalité, des valeurs non-nulles pour le roulis (rotation selon l'axe longitudinal du véhicule) et le 25 tangage (rotation selon l'axe transversal du véhicule) peuvent être le résultat par exemple d'une route en pente mais sont typiquement faibles. Le procédé objet de la présente invention permet de calibrer les trois composantes du vecteur vitesse angulaire w. 30 L'objet 1 est en outre avantageusement équipé de moyens d'acquisition « supplémentaires » 10 d'une vitesse linéaire mesurée de l'objet 1 (notée V), c'est-à-dire du déplacement. Ces moyens 10 peuvent 12 directement ou indirectement permettre d'obtenir la vitesse linéaire, et ainsi être de nombreux types, par exemple des moyens de mesure inertielle. Ainsi le gyromètre 11 peut être complété par un ou plusieurs accéléromètres, voire l'objet 1 peut comprendre une centrale à inertie « 6 5 axes ». Les accéléromètres sont sensibles aux forces extérieures autres que gravitationnelles appliquées sur le capteur, et permettent de mesurer une accélération spécifique notée y. Alternativement, les moyens 10 peuvent consister, si l'objet 1 est un véhicule à roues, en au moins deux odomètres chacun pour une roue du 10 véhicule, par exemple les deux roues arrière, comme cela est représenté dans l'exemple de la figure 1. On note qu'un ensemble d'odomètres est un moyen simple et fiable d'obtenir une simple vitesse linéaire. Par odomètre en entend un équipement capable de mesurer la vitesse d'une roue en comptant les tours (« compte-tours »). Généralement, 15 les odomètres présentent une partie fixée à la roue (par exemple un aimant), et détectent chaque passage de cette partie fixée (appelé « top ») de sorte à compter le nombre de tours par unité de temps, qui est la fréquence de rotation. On connait d'autres techniques, par exemple la détection optique d'une marque sur la roue, ou le magnétomètre du brevet 20 FR2939514 qui détecte la rotation d'un objet métallique tel qu'une roue. Ici la « vitesse » d'une roue est un scalaire, i.e. la norme de la vitesse de la roue dans le référentiel terrestre (dans l'hypothèse d'absence de dérapage). Si le rayon r de la roue est connu, la mesure de la fréquence f de rotation permet d'estimer cette norme de la vitesse : vz2ffr f 25 Encore alternativement, les moyens supplémentaires 10 peuvent être un récepteur GNSS. On comprendra que l'homme du métier n'est limité à aucune technique particulière pour la mesure de vitesse. L'objet 1 comprend en outre comme expliqué des moyens de 30 traitement 21 (typiquement un processeur) pour la mise en oeuvre, en particulier en temps réel, des traitements du présent procédé, par exemple un ordinateur de bord d'un véhicule, et éventuellement une mémoire 22 et 13 une interface 23 pour restituer des informations du mouvement de l'objet 1 (une valeur de vitesse instantanée, un cap, une position sur une carte, etc.), et/ou envoyer des commandes à l'objet 1. L'objet 1 peut à ce titre être notamment un véhicule autonome, et les moyens de traitement 21 5 configurés pour mettre en oeuvre la navigation autonome du véhicule. Ainsi, lesdites commandes sont envoyées aux organes de contrôle du véhicule (moteur, actionneurs du volant, etc.) de sorte à simuler la conduite par le conducteur. A noter que les moyens de traitement 21 peuvent être externes à 10 l'objet 1, et par exemple connectés à ces derniers par un réseau sans fil. Alternativement, le gyromètre 11 et les moyens de mesure magnétique 20 peuvent être connectés aux moyens de traitement de données 21 de manière filaire, par exemple via Ethernet. 15 Procédé Le procédé objet de la présente invention est un procédé de calibration d'au moins le gyromètre 11. Par calibration, on entend la détermination d'un ou plusieurs paramètres de calibration, dont on verra 20 une liste plus loin. En particulier, certains paramètres de calibrations peuvent être considérés fiables et prédéterminés, et d'autres à déterminer. En ce qui concerne ceux à déterminer, on peut prévoir qu'ils présentent des valeurs « actuelles » (en d'autres termes qu'une calibration a déjà eu lieu), et que ces valeurs vont le cas échéant être modifiées (en cas de nouvelle 25 calibration). Dans un mode de réalisation particulièrement préféré, le procédé peut en outre être un procédé de calibration des moyens de mesure magnétique 20, i.e. le gyromètre 11 et les magnétomètres 20 peuvent être calibrés simultanément. Il s'agit d'un mode extrêmement avantageux, 30 puisque comme l'on verra il n'y a même plus besoin de supposer que l'une de la gyrométrie (les données gyrométriques) et de la magnétométrie (les données magnétiques) est fiable pour l'utiliser comme référence pour 14 calibrer l'autre : les deux se calibrent automatiquement l'une l'autre. Alternativement, il sera bien entendu possible de considérer les magnétomètres 20 comme bien calibrés et de calibrer le gyromètre 11 en conséquence, ce qui permet par exemple de calibrer plus de paramètres du 5 gyromètre 11. Comme l'on verra, dans un mode de réalisation avantageux le présent procédé est même un procédé d'estimation du mouvement de l'objet 1, i.e. il comprend suite à la calibration l'utilisation des mesures pour en déduire de façon fiable une ou plusieurs composantes du mouvement. 10 Dans une première étape (a), le procédé comprend l'acquisition par le gyromètre 11 d'une vitesse angulaire mesurée de l'objet 1, notée (mesure) CO.clyro , et par les moyens de mesure magnétique 20, d'au moins deux composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique. Ces composantes sont plus précisément des composantes 15 dites mesurées, formant un vecteur noté M(mesure) dans le cas du champ magnétique (dans le cas d'une dérivée i-ème, on pourra noter V1M De façon préférée, on acquiert trois composantes de la vitesse angulaire et au moins six voire huit composantes magnétiques (avantageusement deux voire trois triaxes magnétométriques positionnés aux coins d'un triangle 20 rectangulaire isocèle), de sorte à pouvoir en déduire le champ et au moins une partie de son gradient (i.e. les dérivées à l'ordre suivant selon les trois axes). Ces grandeurs sont avantageusement mesurées avec un échantillonnage dt (i.e. toutes les « dt » secondes) avec dt très petit devant 25 le temps caractéristique des mouvements de l'objet 1, typiquement 40 ms. Dans une étape (b), les moyens de traitement de données 21 déterminent les valeurs d'au moins un paramètre de calibration du gyromètre 11 minimisant une première expression définie par une vitesse angulaire estimée de l'objet 1 et au moins une première équation 30 magnétique sur les composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique. (mesure)). 15 L'idée est d'estimer séparément des données gyrométriques et des données magnétiques qui sont théoriquement liées par une équation magnétique. Ainsi, grâce à l'équation, on peut exprimer une quantité qui idéalement devrait être nulle (i.e. les données gyrométriques et 5 magnétiques vérifient exactement l'équation magnétique), sinon c'est que la calibration est à parfaire. La vitesse angulaire estimée de l'objet 1, notée cog(eysrtoimation), est fonction de la vitesse angulaire mesurée et de paramètres de calibration du . gyromètre 11, i.e. en général cog(eysrtimation) = 9 (a)g(n;esure) ) avec g une 10 fonction (application) qui est déterminée par les paramètres de calibration. Au cas le plus simple, g est affine et to(estimation) est liée à la vitesse "yr° (mesure) (estimation) angulaire mesurée cogyn, par une formule du type cogyro = D - ((mesure) yyro sure) _ b) , où D et ba représentent les paramètres de calibration yr du gyromètre 11. 15 Dans le cas d'une attitude tridimensionnelle, bgy' est un vecteur de biais et D est une matrice 3x3 = (matrice orthogonale de passage au bon repère) x (matrice triangulaire supérieure contenant les facteurs d'échelle et les calages). De façon simplifiée, on peut par exemple considérer que D est prédéterminée (elle varie en pratique très lentement) et que le seul 20 paramètre de calibration à déterminer pour le gyromètre 11 est bgyro, qui a effectivement tendance à varier dans le temps (on parle de dérive du gyromètre 11). Comme expliqué, les composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique acquises à l'étape (a) sont 25 préférentiellement des composantes mesurées (i.e. potentiellement biaisées), et celles utilisées par ladite première équation magnétiques sont des composantes estimées (le vecteur correspondant est noté Westillulti°n) dans le cas du champ magnétique - dans le cas d'une dérivée i-ème, on pourra noter Vi/t9(estimation)), en fonction des composantes mesurées et de 30 paramètres de calibration des moyens de mesure magnétique 20. 16 Dans le cas où les moyens de mesure magnétique 20 sont des magnétomètres (i.e. on travaille sur les composantes du champ magnétique et pas une dérivée), les composantes estimées M(estimation) du champ magnétique sont avantageusement liées aux composantes mesurées 5 M(mesure) par un modèle de la forme M(mesure) = A. M(estimation) bmagneto, où la matrice 3x3 A et le vecteur bmagneto sont les paramètres de calibration des moyens de mesure magnétique 20. Dans la mesure où c'est plutôt la composante estimée qu'on souhaite exprimer en fonction de la composante mesurée, on peut écrire M(estimation) = A-1 (m(mesure) agneto)- 10 bmagneto est généralement représentatif des effets « type fer dur», et la matrice A représentative des effets « type fer doux ». Ces effets correspondent aux effets fer dur/doux à proprement parler mais peuvent aussi inclure des phénomènes avec le même impact mais une cause différente (par exemple des effets dus à l'électronique/la physique du 15 magnétomètre 20). Comme pour le gyromètre, plus généralement, on peut regarder un modèle d'erreur général : M(estimation) = h(M(mesure)) avec h une fonction (application) qui n'est pas nécessairement affine. Dans le cas affine (m (mesure) = A . m(estimation) + agneto), A est une 20 matrice 3x3, que l'on peut écrire comme A=R-iii avec R orthogonale et triangulaire supérieure. A noter qu'on suppose que les moyens de mesure magnétique 20 en eux-mêmes sont a priori parfaitement calibrés, ce n'est que leur environnement direct (carrosserie, etc.) qui doit être corrigé en calibration 25 afin de pouvoir estimer le champ magnétique terrestre que l'on mesurerait en absence de l'objet 1. La matrice R est ainsi importante et à déterminer en sorte que le repère magnétique objet reste stable. Les méthodes classiques de calibration par rapport à des sphères ne déterminent R qu'à une rotation autour d'un axe près (voir le document E. 30 Dorveaux, D. Vissière, A. P. Mar-tin, N. Petit, "Iterative calibration method for inertial and magnetic sensors", in Proc. of the 48th IEEE Conf. on Decision 17 and Control 2009). Les axes capteurs ne restent pas fixes après l'application de ces méthodes. On peut en pratique déterminer huit coefficients de A en calibration (un coefficient, le facteur d'échelle global qui relie la mesure à l'unité 5 physique en Tesla (ou Gauss), ne peut pas être déterminé, mais en même temps, il n'est pas nécessaire de le connaître pour calibrer le gyromètre 11). Ainsi, dans le mode de réalisation particulièrement préféré de double calibration simultanée, il n'y a que trois paramètres de calibration à déterminer : b9y.j.,, A et bmagneto - Ce nombre restreint de paramètres rend 10 possible une telle double calibration dès lors que l'objet 1 présente une trajectoire « variée » (directions / vitesses qui présentent une grande diversité de valeurs en fonction du temps). Equation magnétique 15 L'évolution du champ magnétique dans le repère de l'objet 1 (dit body) au cours du temps est décrite par l'équation M = M + VM - V + aM at. Le premier terme -w x M décrit le changement du champ dans le 20 repère de l'objet 1 qui est dû à la rotation de ce repère par rapport au repère fixe. Le deuxième terme VM - V décrit le changement du champ mesuré issu du mouvement dans une région avec un champ inhomogène. Finalement, le troisième terme -5M prend en compte les instationnarités du champ at magnétique (par exemple les courants périodiques tels que le 50 Hz en 25 Europe, des mouvements d'un aimant/un objet en acier à proximité dudit objet, etc.). S'il n'y a pas de perturbation magnétique « mobile », i.e. le champ magnétique présente un caractère stationnaire, nia t 0, et on peut réduire l'équation à M = -w x M + VM - V. 18 S'il n'y a pas de perturbation magnétique du tout, i.e. le champ magnétique présente un caractère uniforme en plus de stationnaire, VM 0 et on peut encore réduire l'équation à M = x M. Dans le cas d'une utilisation de capteurs magnétiques pour 5 déterminer un cap, cette double hypothèse d'uniformité et stationnarité est le plus souvent vérifiée ou supposée (comme lorsqu'on utilise une boussole). En effet, pour le champ magnétique terrestre, l'uniformité est localement constatée en très bonne approximation dans nos latitudes, hors les cas de présence par exemple d'une structure métallique ou en béton 10 armé (typiquement un pont) dans les alentours. Les instationnarités dues à des déplacements d'objets à proximité sont encore plus rares car il faudrait que les sources de perturbations se déplacent, et les éventuelles instationnarités dues à des courant périodiques sont en moyenne nulle et peuvent aussi être traitées avec des techniques d'estimation, voir par 15 exemple document C.-I. Chesneau, M. Hillion, and C. Prieur, Motion estimation of a Rigid Body with an EKF using Magneto-Inertial Measurements, 7th Conf. on Indoor Positioning and Indoor Navigation (IPIN'16), Madrid, Spain, 2016. On note qu'on connait la possibilité de généraliser les équations 20 magnétiques (réduites ou non), de sorte à lier non pas les variations spatiales et temporelles du champ magnétique évaluées grâce aux mesures de magnétomètres, mais les variations spatiales et temporelles des gradients magnétiques (i.e. les variations spatiales du champ magnétique) évaluées en particulier grâce aux mesures de gradiomètres. Plus 25 précisément, l'équation M = -a) x M +'7M - V se dérive en V,.,"M = fn GO- VnM + gn (Vll-ki M). V, avec n + 1 > i > 1, où VIIM est une dérivée nième du champ magnétique, et f' et gn des fonctions prédéterminées, voir la demande FR1756675. Dans la suite de la présente description, on utilisera par commodité 30 les équations à l'ordre 0 (c'est-à-dire basées sur M x M +VM - V), mais l'homme du métier saura transposer le procédé aux ordres supérieurs. 19 En référence à la figure 2, le présent procédé de calibration propose astucieusement une calibration à plusieurs niveaux d'hypothèses. Plus précisément, dans l'étape (b) on tente initialement de calibrer le gyromètre en faisant la double hypothèse d'uniformité et de stationnarité du 5 champ, puis dans une éventuelle étape (c) d'estimation d'un paramètre représentatif d'une erreur sur les paramètres de calibration on vérifie l'absence de perturbations, i.e. la pertinence des hypothèses. Dans un mode de réalisation préféré, si cette vérification échoue (ledit paramètre représentatif d'une erreur est supérieur à un seuil 10 prédéterminé), on retente la calibration dans une étape (d) en ne faisant plus que l'hypothèse de stationnarité (et plus celle d'uniformité). On peut alors dans une étape (e) revérifier la pertinence de cette hypothèse. Alternativement, si la vérification de l'étape (c) échoue (ledit paramètre représentatif d'une erreur est supérieur à un seuil prédéterminé), on peut 15 arrêter le procédé et garder la calibration actuelle. On comprend que la présence de l'étape (d) permet de réaliser quand même la calibration dans les éventuels cas où le champ magnétique n'est pas uniforme, à condition d'avoir les moyens d'acquisition supplémentaire 10 évoqués précédemment. Néanmoins, l'invention n'est 20 pas limitée à ce mode de réalisation et on peut se contenter de n'avoir que le gyromètre 11 et les moyens de mesure magnétique 20, et de ne mettre en oeuvre que les étapes (a), (b) et éventuellement (c). L'éventuelle étape (d) peut ainsi comprendre l'acquisition par les moyens d'acquisition supplémentaires 10 d'une vitesse linéaire mesurée de 25 l'objet 1, notée V (en pratique cette acquisition peut être mise en oeuvre concomitamment avec les autres acquisitions de l'étape (a), en particulier à chaque pas de temps), puis la re-détermination par les moyens de traitement de données 21 de valeurs de l'au moins un paramètre de calibration du gyromètre 11 minimisant une deuxième expression définie 30 par la vitesse angulaire estimée de l'objet 1 et au moins une deuxième équation magnétique sur la vitesse linéaire et les composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique, l'au moins 20 une deuxième équation magnétique supposant seulement un caractère stationnaire du champ magnétique au niveau des moyens de mesure magnétique 20. On comprendra ici « acquisition » au sens large, i.e. directe ou 5 indirecte. En pratique les moyens 10 vont le plus souvent mesurer une première grandeur qui n'est pas la vitesse linéaire (des accélérations, la fréquence de rotation d'une roue, etc.), et la vitesse linéaire dite mesurée est en fait estimée à partir de ces grandeurs, éventuellement combinées. L'homme du métier saura obtenir de nombreuses façons une vitesse 10 linéaire de l'objet 1, par exemple par gradiométrie magnétique. L'au moins une deuxième équation supposant seulement un caractère stationnaire du champ magnétique au niveau des moyens de mesure magnétique 20 utilisée à l'étape (d) est ainsi préférentiellement de la forme V,,M = fn(co) x VnM + gn(V,'fil).V , avec co la vitesse angulaire 15 théorique, et comme expliqué M = -w x M + VB - V dans le mode de réalisation préféré dans lequel on mesure directement une composante du champ magnétique. On comprend donc que dans une hypothèse de stationnarité du champ magnétique au niveau des moyens de mesure magnétique 20, la 20 quantité vnm(esimiation) fn ((.09(eysrtoimation)) x vnm(estimation) grena(estimation)) V .n (et en particulier sa version à l'ordre 0 m(esttinutton) w (etimation) x M(estimation) 1Westilnuti°n).V) est égale à gyro zéro aux erreurs de calibration près. 25 Similairement l'au moins une première équation magnétique supposant un caractère uniforme et stationnaire du champ magnétique au niveau des moyens de mesure magnétique 20 utilisée à l'étape (b) est une version simplifiée de l'au moins une deuxième équation, et est préférentiellement de la forme V,.,1/ = fn(w) x VnM, avec comme expliqué 30 M = -w x M dans le mode de réalisation préféré dans lequel on mesure directement des composantes du champ magnétique. On comprend donc 21 que dans une double hypothèse d'uniformité et de stationnarité du champ magnétique au niveau des moyens de mesure magnétique 20, la quantité vnwesiumitto'-, _ inkto(estimation)) x vnm(estimation) (et en particulier sa gyro (4)(estimation) x m(estunation) version à l'ordre 0 heestilnatum) ) est égale à gYro 5 zéro aux erreurs de calibration près. On voit que la première calibration est aisée à mettre en oeuvre car elle ne nécessite pas la connaissance de la vitesse linéaire V. Avantageusement, la première expression (minimisée dans l'étape 10 (b)) est fonction de Westunatwn) s (D (fl e s t i m a t i o n ) x IVI(esttmati°n) et la deuxième i yro expression (minimisée à l'étape (d)) est fonction de 9Yro x _vm(estimation).v1. Sous des conditions génériques, la minimisation de ces expressions donnera les paramètres sans ambiguïté pour une trajectoire générique. On comprendra cependant 15 qu'alternativement l'homme du métier pourra utiliser comme première/deuxième expression toute autre fonctionnelle sensible à des écarts par rapport à la relation définie par la première/deuxième équation (norme L2, etc.). Pour mettre en oeuvre cette minimisation, les moyens de traitement 20 de données 21 peuvent travailler au cours du temps sur un intervalle d'une longueur donnée. A ce titre, de façon connue un filtre récursif (méthodes RLS, recursive least squares, etc.) ou une optimisation (méthode aux moindres carrés, etc.) peut être utilisé. Par exemple, dans l'hypothèse où l'on détermine simultanément les 25 paramètres de calibration des magnétomètres 20, i.e. les paramètres A et bnuigneto, on peut mettre en oeuvre un principe dit de calibration différentielle 2 dans lequel on minimise f Im(estunatton) ci(estimation)gYro x Ag(estimation)1 dt ou Im(estrinatzon)+ (estimation) flm(estrMatton) w(estimation) x m(estimation) vm(estimation). 2 gYro dt, avec 22 w(estemanon) = D co(mesure) b) et 'West ation) = A-1 (m(mesure) gyro 9Yro bmagnéto) aux moindre carrés. Caractérisation d'erreur 5 Comme expliqué, le procédé comprend avantageusement une étape (c) d'estimation d'un paramètre représentatif d'une erreur sur les paramètres de calibration, et encore plus avantageusement (si les étapes (c) et (d) ont eu lieu) une étape (e) de nouvelle estimation dudit paramètre 10 représentatif d'une erreur sur les paramètres de calibration. L'idée est d'estimer la qualité de l'information fournie par les moyens de mesure magnétique 20 pour valider ou non les niveaux d'hypothèses sur le champ magnétique ambiant. Plus précisément, et comme l'on voit sur la 15 figure 2 : si la première estimation (étape (c)) du paramètre représentatif d'une erreur est inférieure audit seuil prédéterminé, on valide la double hypothèse d'uniformité et stationnarité du champ magnétique ambiant et le résultat de la calibration de l'étape (b) est accepté. 20 L'étape (c) comprend alors la calibration effective du gyromètre 11 et/ou des moyens de mesure magnétique 20 avec les valeurs déterminées à l'étape (b) des paramètres de calibration, et il n'y a pas lieu de mettre en oeuvre les étapes (d) et (e) si la première estimation (étape (c)) du paramètre représentatif d'une 25 erreur est supérieure audit seuil prédéterminée, on rejette la double hypothèse d'uniformité et stationnarité du champ magnétique ambiant et le résultat de la calibration de l'étape (b) n'est pas accepté, on peut soit s'arrêter, soit mettre en oeuvre la recalibration de l'étape (d) : 30 o si la deuxième estimation (étape (e)) du paramètre représentatif d'une erreur est cette fois inférieure audit seuil 23 prédéterminée, on valide l'hypothèse de seule stationnarité du champ magnétique ambiant et le résultat de la calibration de l'étape (d) est accepté. L'étape (e) comprend alors la calibration effective du gyromètre 11 et/ou des moyens de 5 mesure magnétique 20 avec les valeurs déterminées à l'étape (d) des paramètres de calibration ; o si la deuxième estimation (étape (e)) du paramètre représentatif d'une erreur est toujours supérieure audit seuil prédéterminée, on rejette même l'hypothèse de stationnarité 10 du champ magnétique ambiant (on est en présence de perturbations) et le résultat de la calibration de l'étape (d) n'est pas accepté. On garde alors les paramètres de l'ancienne calibration comme calibration. On peut alors soit tenter d'utiliser l'équation magnétique complète (avec le terme -a:), 15 soit écarter temporairement les mesures magnétiques. A noter que des valeurs déterminées des paramètres de calibrations lors d'une occurrence de l'étape (b) et/ou (d), mais non utilisées pour la calibration effective, peuvent être stockées sur les moyens de stockage de données 12, et utilisées lors d'une future occurrence de l'étape (c) et/ou (e). 20 Par exemple, on peut prévoir que tant que le paramètre représentatif d'une erreur est au-dessus du seuil, on stocke les paramètres de calibration déterminés, et lorsqu'on passe en dessous du seuil, la calibration effective tient également compte des valeurs stockées. Par ailleurs, le résultat des étapes (c) et (e) peut dépendre de l'âge des paramètres de l'ancienne 25 calibration (par exemple via la valeur du seuil). A noter qu'un intervalle d'erreur accumulée après une durée t depuis la dernière bonne calibration peut être estimé en fonction de cette durée t et de l'erreur estimée dans cette calibration, par exemple en ajoutant à cette erreur un incrément qui est proportionnel à la durée t. L'intervalle d'erreur 30 prédit peut servir pour caractériser le résultat d'une future calibration, et notamment pour définir la valeur dudit seuil prédéterminé : par exemple plus 24 l'intervalle d'erreur est grand, plus il est urgent de recalibrer le gyromètre, et ainsi plus l'on sera « tolérant » avec la valeur du paramètre représentatif d'une erreur (i.e. le seuil prédéterminé sera plus haut). Par ailleurs, dans la mesure où l'intervalle d'erreur augmente proportionnellement plus vite 5 lorsque l'erreur estimée est élevée, cela permet de conserver longtemps les « bonnes » calibrations avec une erreur estimée très faible, tout en forçant plus rapidement à recalibrer les calibrations « moins bonnes ». De façon générale, ledit paramètre représentatif d'une erreur sur les 10 paramètres de calibration est fonction de la valeur de la première/deuxième expression (selon si c'est la première ou deuxième occurrence d'estimation du paramètre), calculée pour les valeurs déterminées des paramètres de calibration. Un premier mode de réalisation de l'étape (c)/(e) est dit intrinsèque 15 puisqu'il n'utilise que des grandeurs disponibles dans l'étape (b)/(d). De façon préférée, on utilise les résidus d'estimation de l'étape (b)/(d), i.e. ledit paramètre représentatif d'une erreur est en particulier la norme (p.ex. L2 ou CD) de la valeur de la première/deuxième expression sur un intervalle de 2 + (0(estimation) temps donné, typiquement f ki(estimatzon) x m(estimation)1 dt gyro 20 ou flwy(eysrtoimation) x Aeestimation) vm(estimation). 2 iwestr;nation) clt. Dans un tel cas, l'étape (c)/(e) peut être mise en oeuvre de façon concomitante avec l'étape (b)/(d). Dans le cas d'un filtre récursif, on peut utiliser la norme (p.ex. L2 ou LOE') de l'innovation du filtre sur une période de temps donné. Alternativement, on peut quantifier la déviation par rapport à une 25 distribution normale, i.e. regarder si la statistique de Westlinatu)n) (estimation). m(estimation) (ou m(estiMatton) (0(estimation).A4(estimation) gyro gyro vm(estimation). V) est similaire à une statistique de bruit gaussien ou non, par exemple en calculant les moments d'ordre supérieurs à 2. Dans un deuxième mode de réalisation de l'étape (c)/(e), dit 30 extrinsèque, on utilise des grandeurs d'un autre type, de sorte à calculer sur 25 la base d'une ou l'autre de ces grandeurs une vitesse angulaire dite de référence, de sorte à la comparer avec la vitesse angulaire estimée. Et l'homme du métier saura calculer une telle vitesse angulaire de référence à partir de divers type de grandeurs, et en particulier on pourra utiliser par 5 exemple : - un angle volant (si l'objet 1 est un véhicule), - des données GNSS, p.ex. GPS, (si/quand disponibles), - des données odométriques. Dans le dernier cas, si les moyens d'acquisition supplémentaire 10 10 d'une vitesse linéaire consistent en au moins deux odomètres (toujours si l'objet 1 est un véhicule), ladite vitesse angulaire théorique peut en effet être obtenue à partir des vitesses mesurées de roues. En particulier, dans le cas où les roues arrières sont chacune munie d'un odomètre, en supposant la géométrie différentielle d'Ackermann, une deuxième vitesse angulaire (est imation) 15 estimée coodo de l'objet 1 est liée aux vitesses mesurées vL, 12R respectivement pour la roue arrière gauche et la roue arrière droite par une (esimation) vo-ni, formule du type coodo = avec ci la distance entre les roues, et a ' ledit paramètre représentatif d'une erreur est alors en particulier la norme ,.,(estimation) - co(estimation)12 (p.ex. L2 ou Lee) de sur l'intervalle de temps uutgyro odo ( 20 donné, où 0)testimation) gyro est la vitesse angulaire pour la rotation seule autour de l'axe vertical de l'objet 1. Selon un troisième mode de réalisation, en ce qui concerne la première occurrence (étape (c)) de l'estimation du paramètre d'erreur, on peut tout simplement estimer tout ou partie d'un gradient spatial de ladite 25 composante du champ magnétique et/ou d'une dérivée i-ième du champ magnétique, i.e. VM ou Vi+1M selon l'ordre de la deuxième équation utilisée. En effet, cette estimation doit être proche de zéro pour que le caractère uniforme soit retenu. On comprend que l'estimation d'un sous-ensemble des composantes du gradient peut en pratique suffire pour cette 26 caractérisation, car chaque composante non-nulle est déjà une manifestation d'une inhomogénéité. Alternativement ou en complément, on peut recourir à de l'apprentissage pour améliorer l'estimation de ce paramètre d'erreur et/ou 5 développer une approche d'identification de conditions favorables (i.e. stationnaires et uniformes) de manière plus robuste et avec une disponibilité augmentée. En particulier, on peut mettre en oeuvre des mécanismes d'apprentissage tels que des réseaux de neurones, machines à vecteurs de 10 support, méthodes des plus proches voisins, forêts d'arbres décisionnels, etc. Ainsi, à chaque occurrence des étapes (b) à (e), on peut venir enrichir une base d'apprentissage dans laquelle chaque ensemble de données de mesures est « taggué » avec la valeur correspondante du paramètre représentatif d'une erreur, de sorte à progressivement (au fur et à mesure 15 des occurrences successives des étapes (b) à (e)) et automatiquement apprendre à distinguer les calibrations acceptables de celles non acceptables. Ainsi, le calibrage s'améliore de lui-même en permanence. Estimation du mouvement 20 Comme expliqué, le procédé comprend avantageusement une étape (f) d'estimation par les moyens de traitement de données 21 du mouvement dudit objet 1 en fonction de la vitesse angulaire mesurée de l'objet 1 et/ou des au moins deux composantes du champ magnétique et/ou d'une dérivée 25 i-ième du champ magnétique, et/ou de l'éventuelle vitesse linéaire mesurée de l'objet 1 (par les moyens 10), et des valeurs des paramètres de calibration, i.e. après recalibration. Cette étape (f), qui n'est pas représentée sur la figure 2 qui se concentre sur la calibration, peut être mise en oeuvre en continu parallèlement. 30 Par estimation du mouvement, on entend en particulier au moins l'estimation d'une orientation de l'objet 1 (dans le plan horizontal, i.e. un cap) et avantageusement l'estimation d'une norme de vitesse. L'orientation 27 s'obtient typiquement par intégration de la vitesse angulaire à partir d'une orientation initiale connue. De façon préférée, l'estimation par les moyens de traitement de données 21 du cap dudit objet 1 est uniquement faite en fonction de la 5 vitesse angulaire mesurée de l'objet 1 et des valeurs des paramètres de calibration, i.e. en utilisant les seules données gyrométriques. En résumé, on utilise les données magnétiques pour calibrer le gyromètre 11, et ensuite on n'utilise pour l'estimation du cap que les données gyrométriques (puisque celles-ci sont toujours disponibles). Le cas échéant, la vitesse 10 linéaire mesurée (obtenue via les moyens 10) n'est utilisée dans l'étape (f) que pour déterminer une vitesse globale de l'objet 1. A noter que dans un cas de véhicule autonome, l'étape (f) peut comprendre la génération d'une commande dudit véhicule 1 en fonction du mouvement estimé, de sorte à amener le véhicule 1 par exemple jusqu'à 15 une destination souhaitée, ou à stopper le véhicule 1 en le gardant dans une trajectoire dénuée d'obstacles. Equipements et objet 20 Selon un deuxième aspect, l'invention concerne en particulier un ensemble des équipements 11, 20, 21 et éventuellement 10 pour la mise en oeuvre de l'un ou l'autre des modes de réalisation du procédé. Cet ensemble peut être installé en kit dans un objet 1 « classique » de sorte à le transformer. Alternativement, en particulier dans le cas où 25 l'objet 1 est un véhicule, ce peut être un véhicule autonome déjà muni de moyens de traitement de données 21 pour la navigation du véhicule, ainsi que de capteurs tels que le gyromètre 11 et/ou des odomètres servant de moyens d'acquisition supplémentaire 10. En particulier est proposé un objet 1, en particulier de type véhicule à 30 roues, comprenant : un gyromètre 11 configuré pour acquérir une vitesse angulaire mesurée de l'objet 1 ; 28 des moyens de mesure magnétique 20 configurés pour acquérir au moins deux composantes du champ magnétique et/ou d'une dérivée i-iéme du champ magnétique ; - le cas échéant, des moyens d'acquisition supplémentaire 10 5 configurés pour acquérir une vitesse linéaire mesurée de l'objet 1 (avantageusement des odomètres équipant au moins deux des roues du véhicule et configurés pour acquérir des vitesses mesurées desdites deux roues) ; - des moyens de traitement de données 21 configurés pour 10 o déterminer des valeurs d'au moins un paramètre de calibration du gyromètre 11 minimisant une première expression définie par une vitesse angulaire estimée de l'objet 1 et au moins une première équation magnétique sur les composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique, 15 - la vitesse angulaire estimée de l'objet 1 étant fonction de la vitesse angulaire mesurée et de paramètres de calibration du gyromètre 11, et - l'au moins une première équation magnétique supposant un caractère uniforme et stationnaire du 20 champ magnétique au niveau des moyens de mesure magnétique 20, les moyens de traitement de données 21 pouvant en outre être configurés pour : o Estimer un paramètre représentatif d'une erreur sur les 25 paramètres de calibration ; et o Si ledit paramètre représentatif d'une erreur est supérieur à un seuil prédéterminé, re-déterminer des valeurs de l'au moins un paramètre de calibration du gyromètre 11 une deuxième expression définie par la vitesse angulaire estimée de l'objet 1 30 et au moins une deuxième équation magnétique sur la vitesse linéaire mesurée et les composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique, l'au moins 29 une deuxième équation magnétique supposant seulement un caractère stationnaire du champ magnétique au niveau des moyens de mesure magnétique 20. 5 Comme expliqué précédemment, l'objet 1 peut en outre comprendre une mémoire 22 et une interface 23, ainsi que d'autre capteurs comme un capteur angulaire du volant ou un récepteur GNSS. Par ailleurs, les moyens de traitement de données 21 peuvent être en outre configurés pour estimer un mouvement dudit objet 1 en fonction de 10 la vitesse angulaire mesurée de l'objet 1, des composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique, et/ou de l'éventuelle vitesse linéaire mesurée de l'objet 1 (le cas échéant selon le résultat d'une comparaison dudit paramètre représentatif d'une erreur avec un seuil), et des valeurs des paramètres de calibration. 15 Produit programme d'ordinateur Selon un troisième et un quatrième aspects, l'invention concerne un produit programme d'ordinateur comprenant des instructions de code pour 20 l'exécution (sur les moyens de traitement 21) d'un procédé de calibration d'un gyromètre 11 selon le premier aspect de l'invention, ainsi que des moyens de stockage lisibles par un équipement informatique (par exemple des moyens de stockage de données 22) sur lequel on trouve ce produit programme d'ordinateur. REVENDICATIONS1. Procédé de calibration d'un gyromètre (11) équipant un objet (1) évoluant dans un champ magnétique ambiant, le procédé étant caractérisé en ce qu'il comprend des étapes de : (a) Acquisition - par le gyromètre (11), d'une vitesse angulaire mesurée de l'objet (1), et - par des moyens de mesure magnétique (20) solidaires dudit objet (1), d'au moins deux composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique, au niveau des moyens de mesure magnétique (20) ; (b) Détermination par des moyens de traitement de données (21) de valeurs d'au moins un paramètre de calibration du gyromètre (11) minimisant une première expression définie par une vitesse angulaire estimée de l'objet (1) et au moins une première équation magnétique sur les composantes du champ magnétique et/ou d'une dérivée iiéme du champ magnétique, - la vitesse angulaire estimée de l'objet (1) étant fonction de la vitesse angulaire mesurée et de paramètres de calibration du gyromètre (11), et - l'au moins une première équation magnétique supposant un caractère uniforme et stationnaire du champ magnétique au niveau des moyens de mesure magnétique (20). 2. Procédé selon la revendication 1, dans laquelle la vitesse angulaire estimée co,(geysitiration) de l'objet (1) est liée à la vitesse angulaire mesurée co(me"") par un modèle co("timati") = D - (cogyro + 9Yro 9Yro gyro bgyro), où D et bgyro sont les paramètres de calibration du gyromètre (11). 31 3. Procédé selon l'une des revendications 1 et 2, comprenant en outre une étape (c) d'estimation d'un paramètre représentatif d'une erreur sur les paramètres de calibration. 5 4. Procédé selon la revendication 3, comprenant en outre, si ledit paramètre représentatif d'une erreur est supérieur à un seuil prédéterminé, une étape (d) d'acquisition par des moyens d'acquisition supplémentaires (10) d'une vitesse linéaire mesurée de l'objet (1), et de re- 10 détermination par les moyens de traitement de données (21) de valeurs de l'au moins un paramètre de calibration du gyromètre (11) minimisant une deuxième expression définie par la vitesse angulaire estimée de l'objet (1) et au moins une deuxième équation magnétique sur la vitesse linéaire mesurée et les composantes du champ magnétique et/ou d'une dérivée i- 15 ième du champ magnétique, l'au moins une deuxième équation magnétique supposant seulement un caractère stationnaire du champ magnétique au niveau des moyens de mesure magnétique (20). 5. Procédé selon la revendication 4, comprenant une 20 étape (e) de nouvelle estimation dudit paramètre représentatif d'une erreur sur les paramètres de calibration. 6. Procédé selon l'une des revendications 3 à 5, dans laquelle ledit paramètre représentatif d'une erreur sur les paramètres de 25 calibration est fonction de la valeur de ladite première expression, ou le cas échéant de ladite deuxième expression, calculée pour les valeurs déterminées des paramètres de calibration. 7. Procédé selon l'une des revendications 3 à 6, dans 30 laquelle l'étape (c), et le cas échéant l'étape (e), comprend la détermination d'un intervalle d'erreur accumulée en fonction de la valeur dudit paramètre représentatif d'une erreur sur les paramètres de calibration et d'une durée 32 depuis la dernière calibration du gyromètre (11), ledit seuil prédéterminé dépendant dudit intervalle d'erreur accumulée. 8. Procédé selon l'une des revendications 1 à 7, dans 5 lequel l'étape (b), et le cas échéant l'étape (d), comprend la mise en oeuvre d'un filtre récursif ou d'une optimisation. 9. Procédé selon l'une des revendications 1 à 8, dans lequel les composantes du champ magnétique et/ou d'une dérivée i-ième 10 du champ magnétique acquises à l'étape (a) sont des composantes mesurées, et celles utilisées par lesdites première et deuxième équations magnétiques sont des composantes estimées en fonction des composantes mesurées et de paramètres de calibration des moyens de mesure magnétique (20), l'étape (b) et, le cas échéant l'étape (d), comprenant 15 également la détermination de valeurs d'au moins un paramètre de calibration des moyens de mesure magnétique (20). 10. Procédé selon la revendication 9, dans laquelle les moyens de mesure magnétique (20) sont des magnétomètres triaxes, les 20 composantes estimées Ileestimattml) du champ magnétique étant liées aux composantes mesurées M(mesure) par un modèle M(mesure) = A - M(estimation) bmagneto, OÙ A et brruz,gneto sont les paramètres de calibration des moyens de mesure magnétique (20). 25 11. Procédé selon l'une des revendications 1 à 10, dans lequel la ou les premières équations magnétiques sont de la forme Vn/v/ = Ma)) x V'M, et le cas échéant la ou les deuxièmes équations magnétiques sont de la forme V,,M = x V,M + gn(V'iM).V, où V,,M est une dérivée n-ième du champ magnétique, w la vitesse angulaire, V la vitesse 30 linéaire, et fuet gn des fonctions prédéterminées. 33 12. Procédé selon la revendication 11, dans lequel la première équation magnétique est Ji/ = -w x M, et le cas échéant la deuxième équation magnétique est M = -w x M + VM - V. 5 13. Procédé selon l'une des revendication 11 et 12, dans lequel la première expression est fonction de Ivn (esiunatton) fri(wil(eysrtoimation)) vnwestimation)1 et l'éventuelle 1 f (w(estimation)) deuxième expression est fonction dev nm(esinnatzon) gyro y njw(estimation) gn (vn+i m (estimation)). 171 10 14. Procédé selon les revendications 6 et 13 en combinaison, dans laquelle ledit paramètre représentatif d'une erreur sur les paramètres de calibration est soit la moyenne de livnii,(esiumitton) - 2 (0(estimation)) x vllm(estimation)1 ou fn(no9(eysrtoimation)) 15 fn( gyro 1v nfi (esizmation) n (estimation) a ',(nv n+ m (estima ion)).v2 sur un intervalle de temps m rt donné, soit tout ou partie d'un gradient spatial desdites composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique. 20 15. Procédé selon l'une des revendications 1 à 14, comprenant une étape (f) d'estimation par les moyens de traitement de données (21) du mouvement dudit objet (1) en fonction de la vitesse angulaire mesurée de l'objet (1), des composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique, et/ou de l'éventuelle 25 vitesse linéaire mesurée de l'objet (1), et des valeurs des paramètres de calibration, dans lequel une orientation du véhicule (1) est estimée à l'étape (f) uniquement en fonction de la vitesse angulaire mesurée du véhicule (1) et des valeurs des paramètres de calibration. 34 16. Objet (1) évoluant dans un champ magnétique ambiant, comprenant un gyromètre (11) configuré pour acquérir une vitesse angulaire mesurée de l'objet (1), des moyens de mesure magnétique (20) configurés pour acquérir au moins deux composantes du champ 5 magnétique et/ou d'une dérivée i-ième du champ magnétique, l'objet (1) étant caractérisé en ce qu'il comprend en outre des moyens de traitement de données (21) configurés pour : - déterminer des valeurs d'au moins un paramètre de calibration du gyromètre (11) minimisant une première expression définie par une 10 vitesse angulaire estimée de l'objet (1) et au moins une première équation magnétique sur les composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique, o la vitesse angulaire estimée de l'objet (1) étant fonction de la vitesse angulaire mesurée et de paramètres de calibration du 15 gyromètre (11), et o l'au moins une première équation magnétique supposant un caractère uniforme et stationnaire du champ magnétique au niveau des moyens de mesure magnétique (20). 20 17. Objet selon la revendication 16, comprenant en outre des moyens d'acquisition supplémentaires (10) configurés pour acquérir une vitesse linéaire mesurée de l'objet (1), les moyens de traitement de données (21) étant en outre configurés pour : - Estimer un paramètre représentatif d'une erreur sur les paramètres 25 de calibration et - Si ledit paramètre représentatif d'une erreur est supérieur à un seuil prédéterminé, re-déterminer des valeurs de l'au moins un paramètre de calibration du gyromètre (11) une deuxième expression définie par la vitesse angulaire estimée de l'objet (1) et au moins une 30 deuxième équation magnétique sur la vitesse linéaire mesurée et les composantes du champ magnétique et/ou d'une dérivée i-ième du champ magnétique, l'au moins une deuxième équation magnétique supposant seulement un caractère stationnaire du champ magnétique au niveau des moyens de mesure magnétique (20). 18. Objet selon la revendication 17, étant un véhicule à 5 roues, les moyens d'acquisition supplémentaire (10) étant au moins deux odomètres équipant des roues du véhicule. 19. Produit programme d'ordinateur comprenant des instructions de code pour l'exécution d'un procédé de calibration d'un 10 gyromètre (11) selon l'une des revendications 1 à 15, lorsque ledit programme est exécuté sur un ordinateur. 20. Moyen de stockage lisible par un équipement informatique sur lequel un produit programme d'ordinateur comprend des 15 instructions de code pour l'exécution d'un procédé de calibration d'un gyromètre (11) selon l'une des revendications 1 à 15.