Brevet FR3077877A1 — publié le 20190816

PROCÉDÉ DE CALIBRATION D'UN CAPTEUR TRI-AXE AVEC SELECTION D'UNE MÉTHODE DE CALIBRATION EN FONCTION DE LA DISTRIBUTION SPATIALE DES MESURES DESCRIPTION DOMAINE TECHNIQUE Le domaine de l'invention est celui de la calibration d'un capteur tri-axe. L'invention concerne plus particulièrement une calibration automatique réalisée de manière continue alors que le capteur est en cours d'utilisation. ÉTAT DE LA TECHNIQUE ANTÉRIEURE Un magnétomètre peut être intégré à une centrale inertielle de navigation pour réaliser une mesure du champ magnétique terrestre, ce qui permet de déterminer le cap d'un dispositif embarquant ladite centrale. En l'absence d'erreurs et de perturbation locales, le magnétomètre mesure uniquement le vecteur de champ magnétique terrestre. Toutefois, lorsque la centrale inertielle est contaminée par certaines erreurs, la mesure du vecteur de champ magnétique terrestre est différente de la valeur réelle du vecteur de champ magnétique terrestre, appelée valeur de référence. Les erreurs présentes sur les mesures du magnétomètre ne se limitent pas qu'aux erreurs intrinsèques du capteur, comme le facteur d'échelle, le mauvais alignement et le biais, mais peuvent également être causées par l'environnement extérieur, notamment les erreurs de fer dur qui proviennent d'une présence de matériaux magnétisés, tels que ceux présents dans les composants électroniques, le circuit d'alimentation, etc., et les erreurs de fer doux qui sont générées par la déformation des lignes de champ magnétique dans l'environnement extérieur au capteur, par exemple due à la présence de câbles, de vis et d'écrous métalliques, etc. Le magnétomètre doit donc être calibré pour corriger les mesures brutes des erreurs mentionnées ci-dessus. Cette problématique se pose également pour les autres capteurs de la centrale inertielle, à savoir pour l'accéléromètre et le gyromètre. 2 La plupart des algorithmes de calibration reposent sur une méthode d'ajustement des mesures brutes sur une surface paramétrique, typiquement un ellipsoïde, l'écart entre cet ellipsoïde et une sphère de rayon égal à la norme du champ magnétique terrestre dans le cas d'un magnétomètre étant représentatif des erreurs 5 affectant le magnétomètre. La transformation qui permet de passer de l'ellipsoïde à la sphère est alors appliquée aux mesures brutes pour les corriger des erreurs susmentionnées. Pour que ces algorithmes de calibration soient efficaces et à même de fournir des paramètres de calibration suffisamment précis, il est nécessaire de pouvoir 10 acquérir des mesures bien distribuées, à savoir des points couvrant l'intégralité de l'ellipsoïde. Or en pratique, en cours d'utilisation du magnétomètre, les mesures qui peuvent être collectées sont généralement mal distribuées et ne couvrent pas certaines parties de l'ellipsoïde. Cette constatation est notamment faite dans l'article de A. Olivares et 15 al. "Automatic Determination of Validity of Input Data Used in Ellipsoid Fitting MARG Calibration Algorithms," Sensors, vol. 13, no. 9, pp. 11797-11817, Sep. 2013. Les auteurs de cet article montrent ainsi qu'une erreur des algorithmes d'ajustement à un ellipsoïde peut être déterminée en analysant la distribution spatiale des mesures brutes du magnétomètre. Ces auteurs proposent alors de ne mettre en oeuvre ces algorithmes de 20 calibration que lorsque l'erreur attendue est inférieure à un seuil. Cet article indique que le meilleur indicateur de validité de la distribution des mesures brutes est celui de l'enveloppe convexe. Cet indicateur présente toutefois l'inconvénient de nécessiter d'importantes ressources de calcul. EXPOSÉ DE L'INVENTION 25 L'invention a pour objectif de perfectionner les techniques connues de calibration, et propose pour ce faire un procédé de calibration d'un capteur tri-axe comprenant des étapes de calcul d'un indicateur de distribution spatiale d'un ensemble de mesures délivrées par le capteur tri-axe et de détermination, selon une méthode d'ajustement, de paramètres d'une surface paramétrique permettant d'ajuster ladite 3 surface audit ensemble de mesures. La méthode d'ajustement est plus particulièrement sélectionnée parmi plusieurs méthodes d'ajustement prédéterminées en fonction dudit indicateur de distribution spatiale. Certains aspects préférés mais non limitatifs de ce procédé sont les 5 suivants : - les méthodes d'ajustement prédéterminées comprennent ou consistent en une méthode d'ajustement à un ellipsoïde en une étape, une méthode d'ajustement à une sphère en une étape et une méthode d'ajustement en deux étapes consistant en un premier ajustement à une sphère suivi d'un deuxième ajustement à un ellipsoïde dont le 10 centre est confondu avec le centre de la sphère du premier ajustement ; - lorsque l'indicateur de distribution spatiale est inférieur à un premier seuil, la méthode d'ajustement sélectionnée est la méthode d'ajustement à un ellipsoïde en une étape, lorsque l'indicateur de distribution spatiale est supérieur au premier seuil et inférieur à un deuxième seuil, la méthode d'ajustement sélectionnée est la méthode 15 d'ajustement en deux étapes, et lorsque l'indicateur de distribution spatiale est supérieur au deuxième seuil, la méthode d'ajustement sélectionnée est la méthode d'ajustement à une sphère en une étape ; - l'étape de calcul de l'indicateur de distribution spatiale comprend des sous-étapes de détermination d'une matrice de corrélation des mesures dudit ensemble et de 20 détermination du nombre de conditionnement de la matrice de corrélation ; - il comprend une étape préalable de sélection des mesures dudit ensemble parmi les mesures délivrées par le capteur tri-axe, ladite étape préalable comprenant une estimation d'un changement d'orientation du capteur tri-axe entre les mesures délivrées et la sélection de mesures associées à des orientations du capteur jugées, au 25 moyen d'une comparaison à un seuil, comme étant suffisamment différentes ; - le capteur tri-axe est un magnétomètre et l'estimation d'un changement d'orientation du capteur entre les mesures comprend la détermination de l'écart entre l'amplitude d'une composante d'une mesure avec les amplitudes de composantes de mesures préalablement sélectionnées ; 4 - le capteur tri-axe est un magnétomètre et l'estimation d'un changement d'orientation du magnétomètre entre les mesures exploite des mesures délivrées par un gyromètre, ou les mesures délivrées par le magnétomètre et des mesures délivrées par un accéléromètre ; 5 - il comprend la détermination d'une transformation test permettant de transformer la surface paramétrique ajustée audit ensemble de mesures en une sphère unité, l'application de la transformation test aux mesures dudit ensemble pour fournir des mesures transformées et la détermination d'un degré d'ajustement des mesures transformées à la sphère unité ; 10 - il comprend en outre l'application d'une transformation de calibration courante aux mesures dudit ensemble pour fournir des mesures calibrées, la détermination d'un degré d'ajustement des mesures calibrées à la sphère unité et la sélection de la transformation test en tant que nouvelle transformation de calibration courante si le degré d'ajustement des mesures transformées à la sphère unité est supérieur au degré 15 d'ajustement des mesures calibrées à la sphère unité. BRÈVE DESCRIPTION DES DESSINS D'autres aspects, buts, avantages et caractéristiques de l'invention apparaîtront mieux à la lecture de la description détaillée suivante de formes de réalisation préférées de celle-ci, donnée à titre d'exemple non limitatif, et faite en 20 référence aux dessins annexés sur lesquels : - les figures la et lb représentent des exemples de distribution spatiale de mesures délivrées par un magnétomètre ; - la figure 2 est un ordinogramme représentant les différentes étapes du procédé selon l'invention ; 25 - les figures 3a et 3b représentent respectivement une surface paramétrique ajustée à un ensemble de mesures d'un magnétomètre selon une méthode d'ajustement à un ellipsoïde en une étape et les mesures calibrées par application de la transformation permettant de passer de cette surface à la sphère unité ; 5 - les figures 4a et 4b représentent respectivement une surface paramétrique ajustée à un ensemble de mesures du magnétomètre selon une méthode d'ajustement en deux étapes et les mesures calibrées par application de la transformation permettant de passer de cette surface à la sphère unité, l'ensemble de mesures étant 5 identique à celui des figures 3a et 3b. EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS L'invention porte sur un procédé de calibration d'un capteur tri-axe et plus particulièrement sur une calibration automatique réalisée alors que le capteur tri-axe est en cours d'utilisation. Le capteur est typiquement intégré au sein d'une centrale 10 inertielle de navigation. On prendra dans ce qui suit l'exemple non limitatif d'un magnétomètre, l'invention s'étendant à d'autres types de capteur, notamment à un accéléromètre et/ou à un gyromètre. La mesure X des trois composantes du champ magnétique terrestre Br réalisée par un magnétomètre est affectée d'erreurs qui peuvent être modélisées sous la 15 forme matricielle suivante : X=W.Br+b, où b représente un biais principalement introduit par les interférences de fer dur et W est une matrice 3*3 pouvant se décomposer selon W=Wgain.Wmis.Wsoft avec Wgain une matrice diagonale modélisant les erreurs de facteur d'échelle, Wol,s une matrice modélisant les erreurs de défauts d'alignement (non-orthogonalité) des axes de mesure et Wsoft une matrice modélisant les erreurs introduites gx 20 par les interférences de fer doux. On a ainsi Wgain = 0 gy 0 et 0 gz ( 1 - ez Ey WIIIIS = e Z 1 -Ex J. La calibration du magnétomètre vise à identifier W et b pour -ez EX 1 ensuite calculer des mesures calibrées débarrassées des erreurs : X,=W-1(X-b). En l'absence de perturbations, le magnétomètre mesure un champ 3D constant. En modifiant l'orientation du magnétomètre dans ce champ constant, les 25 mesures devraient être portées par une sphère. Mais les biais, les facteurs d'échelle et les défauts d'alignement transforment cette sphère en une surface quadratique, typiquement un ellipsoïde. Une technique de calibration consiste alors à déterminer les 6 paramètres d'une surface paramétrique qui permettent d'ajuster ladite surface aux mesures, puis de déterminer la transformation inverse qui va transformer cette surface en la sphère unité. Cette transformation inverse est alors appliquée aux mesures, de manière à fournir des mesures calibrées débarrassées des erreurs susmentionnées. 5 Comme indiqué précédemment, pour que cette technique de calibration permette de fournir des paramètres de calibration (D=W-1 et b) suffisamment précis, il est nécessaire que les mesures soient spatialement bien distribuées en ce sens qu'elles couvrent l'intégralité de l'ellipsoïde. L'invention propose de calculer un indicateur de la bonne ou mauvaise distribution spatiale des mesures, et d'utiliser cet indicateur 10 comme critère pour sélectionner une méthode d'ajustement d'une surface paramétrique aux mesures parmi plusieurs méthodes d'ajustement prédéterminées. Comme rappelé dans l'article précité de A. Olivares et al., une équation de la forme X=W.Br+b s'applique également aux autres capteurs de la centrale inertielle, en remplaçant Br par le champ gravitationnel terrestre dans le cas de l'accéléromètre (sa 15 norme est généralement de 1 si l'unité de mesure de l'accélération est g) et en remplaçant Br par une vitesse angulaire constante dans le cas du gyromètre. On a représenté sur les figures la et lb des exemples de distribution spatiale de mesures délivrées par le magnétomètre. Sur la figure la, les mesures sont quasiment coplanaires et donc particulièrement mal distribuées. Au contraire, sur la 20 figure lb, les mesures couvrent correctement la surface d'un ellipsoïde et apparaissent donc bien distribuées. Selon l'invention, la méthode d'ajustement utilisée pour déterminer les paramètres de la surface paramétrique à ajuster aux mesures sera différente pour les mesures des figures la et lb. En référence à la figure 2, le procédé selon l'invention comprend plus 25 particulièrement une étape de calcul DIS d'un indicateur de distribution spatiale d'un ensemble de mesures délivrées par le magnétomètre. Le procédé comprend en outre une étape de détermination AJS, selon une méthode d'ajustement, de paramètres d'une surface paramétrique permettant d'ajuster la surface paramétrique à l'ensemble de mesures délivrées par le magnétomètre. La méthode d'ajustement est sélectionnée parmi 7 plusieurs méthodes d'ajustement prédéterminées en fonction de l'indicateur de distribution spatiale calculé lors de l'étape de calcul DIS. Les méthodes d'ajustement prédéterminées peuvent notamment comprendre une méthode d'ajustement à un ellipsoïde en une étape, une méthode 5 d'ajustement à une sphère en une étape et une méthode d'ajustement en deux étapes consistant en un premier ajustement à une sphère suivi d'un deuxième ajustement à un ellipsoïde dont le centre est confondu avec le centre de la sphère du premier ajustement. Dans une variante de réalisation, les méthodes d'ajustement prédéterminées peuvent consister en ces trois méthodes précitées. 10 La méthode d'ajustement à un ellipsoïde peut comprendre un ajustement à un ellipsoïde tourné (« rotated ellipsoid ») nécessitant la détermination de 9 paramètres ou un ajustement à un ellipsoïde non tourné (« non-rotated ellipsoid ») nécessitant la détermination de 6 paramètres. Un exemple d'une méthode d'ajustement en deux étapes est 15 notamment décrit dans l'article de M. Gietzelt et al., "Performance comparison of accelerometer calibration algorithms based on 3D-ellipsoid fitting methods" Comput. Methods Programs Biomed., vol. 111, no. 1, pp. 62-71, Jul. 2013. La méthode d'ajustement à une sphère nécessite la détermination de 4 paramètres. Elle est utilisée pour calculer les biais liés aux interférences de fer dur et 20 approximer les facteurs d'échelle. Elle fournit une calibration considérée comme un mode dégradé. A titre d'exemple, l'indicateur de distribution spatiale est d'autant plus faible que les mesures sont bien distribuées. Alors, dans une réalisation possible, lorsque l'indicateur de distribution spatiale est inférieur à un premier seuil, la méthode 25 d'ajustement sélectionnée est la méthode d'ajustement à une ellipsoïde en une étape, lorsque l'indicateur de distribution spatiale est supérieur au premier seuil et inférieur à un deuxième seuil, la méthode d'ajustement sélectionnée est la méthode d'ajustement en deux étapes, et lorsque l'indicateur de distribution spatiale est supérieur au deuxième seuil, la méthode d'ajustement sélectionnée est la méthode d'ajustement à une sphère en une étape. Bien entendu, d'autres seuils peuvent être prévus, comme par exemple un 8 seuil pour départager un ajustement à un ellipsoïde tourné d'un ajustement à un ellipsoïde non tourné. Dans un mode de réalisation de l'invention, l'invention adopte une approche de type analyse en composantes principales pour évaluer la dispersion des 5 mesures. Plus précisément, le calcul DIS de l'indicateur de distribution spatiale comprend la détermination de la matrice de corrélation des mesures de l'ensemble de mesures et la détermination du nombre de conditionnement de la matrice de corrélation. On considère un ensemble de N mesures des trois composantes du champ magnétique et Xm la matrice 3*N des mesures centrées. Le centre c des mesures 10 est calculé par la méthode "min-max-average". Par exemple pour l'axe x, c, = (mxmin + mx,m')/2, où mx,min et mx,m' représentent respectivement le minimum et le maximum de l'ensemble des mesures sur l'axe x. Alternativement, le centre des mesures peut être calculé en utilisant la moyenne ou la médiane des mesures, mais en pratique, l'ensemble des mesures est souvent mal distribué, et le centre peut se trouver trop éloigné de la 15 moyenne (ou médiane). La matrice de corrélation S s'exprime selon S =-Nli(XmX). Les valeurs propres Ài (iE{1,2,3}) de la matrice S sont représentatives des variances des mesures le long de chaque vecteur propre. Les valeurs propres Ài contiennent donc de l'information utile à propos de la dispersion des mesures. Plus particulièrement, le 20 nombre de conditionnement de la matrice S, à savoir le ratio de la valeur propre la plus grande sur la valeur propre la plus faible, est d'autant plus proche de 1 que les mesures sont bien distribuées et d'autant plus grand que les mesures sont mal distribuées. Ce nombre de conditionnement s'écrit r = max'. Ai et est utilisé pour sélectionner la méthode mini d'ajustement à utiliser. Par exemple, lorsque r est inférieur à un premier seuil r1, la 25 méthode d'ajustement sélectionnée est la méthode d'ajustement à un ellipsoïde en une étape, lorsque r est supérieur au premier seuil r1 et inférieur à un deuxième seuil r2, la méthode d'ajustement sélectionnée est la méthode d'ajustement en deux étapes, et lorsque r est supérieur au deuxième seuil r2, la méthode d'ajustement sélectionnée est la méthode d'ajustement à une sphère en une étape. 9 Le procédé selon l'invention peut comprendre une étape préalable SLT de sélection des mesures dudit ensemble parmi les mesures délivrées par le magnétomètre. Cette étape préalable comprend une estimation d'un changement d'orientation du magnétomètre entre les mesures et la sélection de mesures associées à 5 des orientations du magnétomètre jugées, par exemple au moyen d'une comparaison à un seuil, comme étant suffisamment différentes. En s'assurant que les mesures sélectionnées sont réalisées pour différentes orientations du magnétomètre, on favorise le fait que ces mesures puissent couvrir un ellipsoïde. L'estimation d'un changement d'orientation du magnétomètre peut 10 exploiter les mesures délivrées par le magnétomètre. Cette estimation peut comprendre l'évaluation d'une mesure en venant déterminer, pour chaque composante de la mesure évaluée, l'écart K - .721 entre l'amplitude de la composante de la mesure évaluée et l'amplitude .72 d'une composante d'une mesure préalablement sélectionnée dans ledit ensemble. Si pour toutes les mesures préalablement sélectionnées dans ledit ensemble, 15 cet écart est supérieur à un seuil prédéterminé xmag, par exemple égal à 10uT, la mesure évaluée est sélectionnée. En pratique, il suffit de comparer avec l'amplitude de la composante de la dernière mesure sélectionnée (ici désignée par l'indice n). L'estimation d'un changement d'orientation du magnétomètre peut exploiter les mesures délivrées par un gyromètre, par exemple lorsque le magnétomètre 20 est intégré dans une centrale inertielle. Une mesure magnétométrique associée à une mesure du vecteur vitesse angulaire délivrée par le gyromètre dont une composante co E fco, co co 1 est supérieure à un seuil Tgyr , par exemple de 30°/s, est alors x y, z sélectionnée. En variante, les mesures délivrées par le gyromètre sont utilisées pour 25 estimer une différence d'orientation du gyromètre entre un instant actuel et un instant précédent correspondant à la dernière sélection d'une mesure. Selon cette variante, on procède à l'estimation de l'attitude relative du gyromètre au moyen d'une intégration temporelle des mesures délivrées par le gyromètre, en considérant que les biais du gyromètre sont compensés ou simplement en les négligeant. La formule de dérivation du 30 quaternion, (t) = -21 q(t) 0 Co(t), peut être discrétisée selon la méthode d'Euler pour 10 calculer l'orientation actuelle q(t) à partir de (.00i) pour n E [ta, t], avec ti, l'instant correspondant à la dernière sélection d'une mesure. Une erreur ge est calculée, représentative d'une différence d'orientation entre le quaternion actuel q(t) et le quaternion associée à la dernière mesure sélectionnée q(t), selon ge = q(t) 0 go* . 5 L'angle de ge dans la représentation axe-angle (« axis-angle ») du quaternion ge est donné par a = cos-1- ge,o, où go:, est la partie réelle de l'erreur qe. La mesure magnétométrique actuelle est sélectionnée si a > TgyrInt, avec TgyrInt un seuil par exemple de 300 . L'estimation d'un changement d'orientation du magnétomètre peut 10 également exploiter les mesures délivrées par le magnétomètre et par un accéléromètre, par exemple lorsque le magnétomètre est intégré dans une centrale inertielle. La méthode TRIAD peut à ce titre être utilisée pour rechercher des mesures magnétométriques réalisées à des angles suffisamment différents. Comme pour l'approche précédente, on calcule un angle a en tant que différence entre l'orientation 15 actuelle q(t) et l'orientation correspondant associée à la dernière mesure sélectionnée, a = cos-1- ge,o, et la mesure magnétométrique actuelle est sélectionnée si a > T - acc, avec Tacc un seuil par exemple de 30°. Cette approche ne nécessite pas de stocker les mesures gyromètriques pour réaliser l'intégration temporelle puisque des orientations absolues sont calculées à partir des mesures délivrées par le magnétomètre et par l'accéléromètre. 20 Cette approche repose toutefois sur l'hypothèse d'une absence d'accélération linéaire. Une ou plusieurs autres contraintes peuvent être utilisées lors de cette étape préalable de sélection des mesures magnétométriques pour constituer ledit ensemble servant de référence pour l'ajustement d'une surface paramétrique. Une première contrainte consiste par exemple en une durée minimum entre deux mesures 25 successivement sélectionnées, par exemple 0,1s. Une autre contrainte consiste par exemple une durée maximum entre deux mesures successivement sélectionnées, par exemple 1s. Si cette durée est dépassée, l'étape préalable est réinitialisée. Encore une autre contrainte concerne le nombre de mesures à sélectionner, avec par exemple un nombre minimal et un nombre maximal de mesures, par exemple entre 100 et 200 30 mesures devant être sélectionnées lors de l'étape préalable. 11 Suite à l'étape de détermination AJS, selon la méthode d'ajustement sélectionnée, des paramètres de la surface paramétrique, le procédé selon l'invention peut comprendre une étape CHK de vérification de la qualité de l'ajustement. Cette étape de vérification CHK peut comprendre l'étude de propriétés 5 analytiques des matrices impliquées dans l'ajustement. Cette étude peut ainsi venir vérifier si la matrice représentant les mesures est une matrice de plein rang au sens des moindres carrés. En effet, des mesures strictement coplanaires ne forment pas une matrice de plein rang. Cette étude peut également venir étudier la matrice représentant la surface paramétrée à l'aide des paramètres déterminés à l'étape AJS. Les méthodes 10 d'ajustement à un ellipsoïde ont effectivement pour objectif d'identifier une matrice symétrique définie positive A apparaissant dans l'équation (X - b)T A(X - b) = c. Si la matrice identifiée n'est pas définie positive, la qualité de l'ajustement est jugée insuffisante. L'étape de vérification CHK de vérification de la qualité de l'ajustement 15 peut également comprendre la détermination d'une transformation test permettant de transformer la surface paramétrique ajustée audit ensemble de mesures (i.e. la surface paramétrée à l'aide des paramètres déterminés à l'étape AJS), en une sphère, l'application de la transformation test aux mesures dudit ensemble pour fournir des mesures transformées et la détermination d'un degré d'ajustement des mesures 20 transformées à la sphère. La transformation test vise par exemple à transformer la surface paramétrée en une sphère unité et la détermination d'un degré d'ajustement des données transformées à la sphère comprend le calcul d'une erreur quadratique moyenne des distances entre la sphère unité et les mesures transformées. La norme de la i-ème 25 mesure transformée X (sa distance algébrique à la sphère unité) s'exprime selon IIX/,c112 = - 0112. La somme des erreurs quadratiques pour les N mesures de l'ensemble considéré s'écrit X2 = - b)112 - 1)2, et l'erreur quadratique moyenne s'exprime selon p=.'lx2/N. Cette erreur est confrontée à un seuil (par exemple de 5%, ce qui correspond à une erreur d'environ 2 uT pour le champ 12 magnétique terrestre d'amplitude 47 pT), pour décider si la qualité de l'ajustement est suffisante ou non. Dans un mode de réalisation possible, on vient évaluer une éventuelle amélioration de la calibration. Pour ce faire, le procédé comprend l'application d'une 5 transformation de calibration courante (i.e., la transformation actuellement appliquée aux mesures brutes) aux mesures dudit ensemble pour fournir des mesures calibrées, la détermination d'un degré d'ajustement des mesures calibrées à la sphère unité et la sélection UPD de la transformation test à la sphère unité en tant que nouvelle transformation de calibration courante si le degré d'ajustement des mesures 10 transformées à la sphère unité est supérieur au degré d'ajustement des mesures calibrées à la sphère unité. La transformation courante est appliquée à l'intégralité des mesures délivrées par le magnétomètre, jusqu'à ce que lui soit substituée une nouvelle transformation associée à une surface paramétrique jugée mieux ajustée aux mesures. La transformation test peut alternativement transformer la surface 15 paramétrée en une sphère de rayon Bo, avec Bo la norme du champ géomagnétique, et la détermination d'un degré d'ajustement des mesures transformées à la sphère peut comprendre le calcul d'une différence de norme entre les mesures transformées X, et la sphère de rayon Bo: I IlXclI2 - Bo I. On évalue une caractéristique de cette différence de norme, par exemple son maximum ou l'un de ses percentiles (par exemple à 95%), que 20 l'on compare à un seuil choisi au-dessus du niveau de bruit, par exemple de 5uT, pour décider si la qualité de l'ajustement est suffisante ou non. Par ailleurs, cette caractéristique peut être exploitée pour autoriser la mise en oeuvre de la calibration automatique selon l'invention en cas d'observation d'une différence de norme trop importante, par exemple supérieure à 2p.T. 25 On a représenté sur les figures 3a et 4a une surface paramétrique ajustée à un ensemble de mesures du magnétomètre respectivement selon une méthode d'ajustement à un ellipsoïde en une étape et selon une méthode d'ajustement en deux étapes, l'indicateur r de distribution spatiale des mesures étant ici égal à 9.98. On a représenté sur la figure 3b, respectivement sur la figure 4b, les mesures calibrées par 30 application de la transformation permettant de passer de la surface de la figure 3a, 13 respectivement de la figure 4a, à la sphère unité. On constate que la méthode d'ajustement en deux étapes est ici plus appropriée que la méthode d'ajustement à un ellipsoïde en une étape. L'invention n'est pas limitée au procédé tel que précédemment décrit, 5 mais s'étend également à un produit programme d'ordinateur comprenant des instructions pour l'exécution de ce procédé lorsque ledit programme est exécuté sur un ordinateur. Et l'invention s'étend également à un dispositif, par exemple une centrale de navigation inertielle, REVENDICATIONS1. Procédé de calibration d'un capteur tri-axe, comprenant des étapes de calcul (DIS) d'un indicateur de distribution spatiale d'un ensemble de mesures délivrées par le capteur tri-axe et de détermination (AJS), selon une méthode d'ajustement, de paramètres d'une surface paramétrique permettant d'ajuster ladite surface paramétrique audit ensemble de mesures délivrées, caractérisé en ce que la méthode d'ajustement est sélectionnée parmi plusieurs méthodes d'ajustement prédéterminées en fonction dudit indicateur de distribution spatiale. 2. Procédé selon la revendication 1, dans lequel les méthodes d'ajustement prédéterminées comprennent une méthode d'ajustement à un ellipsoïde en une étape, une méthode d'ajustement à une sphère en une étape et une méthode d'ajustement en deux étapes consistant en un premier ajustement à une sphère suivi d'un deuxième ajustement à un ellipsoïde dont le centre est confondu avec le centre de la sphère du premier ajustement. 3. Procédé selon la revendication 2, dans lequel lorsque l'indicateur de distribution spatiale est inférieur à un premier seuil, la méthode d'ajustement sélectionnée est la méthode d'ajustement à un ellipsoïde en une étape, lorsque l'indicateur de distribution spatiale est supérieur au premier seuil et inférieur à un deuxième seuil, la méthode d'ajustement sélectionnée est la méthode d'ajustement en deux étapes, et lorsque l'indicateur de distribution spatiale est supérieur au deuxième seuil, la méthode d'ajustement sélectionnée est la méthode d'ajustement à une sphère en une étape. 4. Procédé selon l'une des revendications 1 à 3, dans lequel l'étape de calcul (DIS) de l'indicateur de distribution spatiale comprend des sous-étapes de détermination 15 d'une matrice de corrélation des mesures dudit ensemble et de détermination du nombre de conditionnement de la matrice de corrélation. 5. Procédé selon l'une des revendications 1 à 4, comprenant une étape préalable 5 de sélection (SLT) des mesures dudit ensemble parmi les mesures délivrées par le capteur tri-axe, ladite étape préalable comprenant une estimation d'un changement d'orientation du capteur tri-axe entre les mesures délivrées et la sélection de mesures associées à des orientations du capteur jugées, au moyen d'une comparaison à un seuil, comme étant suffisamment différentes. 10 6. Procédé selon la revendication 5, dans lequel le capteur tri-axe est un magnétomètre et l'estimation d'un changement d'orientation du capteur entre les mesures délivrées comprend une détermination d'un écart entre l'amplitude d'une composante d'une mesure délivrée avec l'amplitude d'une composante d'au moins 15 une des mesures préalablement sélectionnées. 7. Procédé selon la revendication 5, dans lequel dans lequel le capteur tri-axe est un magnétomètre et l'estimation d'un changement d'orientation du magnétomètre entre les mesures délivrées exploite des mesures délivrées par un gyromètre ou 20 exploite à la fois les mesures délivrées par le magnétomètre et des mesures délivrées par un accéléromètre. 8. Procédé selon l'une des revendications 1 à 7, comprenant une détermination d'une transformation test permettant de transformer la surface paramétrique ajustée 25 audit ensemble de mesures délivrées en une sphère unité, l'application de la transformation test aux mesures dudit ensemble pour fournir des mesures transformées et une détermination d'un degré d'ajustement des mesures transformées à la sphère unité. 9. Procédé selon la revendication 8, comprenant en outre une application d'une transformation de calibration courante aux mesures dudit ensemble pour fournir des mesures calibrées, une détermination d'un degré d'ajustement des mesures calibrées 16 à la sphère unité et la sélection (UPD) de la transformation test en tant que nouvelle transformation de calibration courante si le degré d'ajustement des mesures transformées à la sphère unité est supérieur au degré d'ajustement des mesures calibrées à la sphère unité. 5 10. Programme d'ordinateur comprenant des instructions pour l'exécution du procédé selon l'une des revendications 1 à 9 lorsque ledit programme est exécuté sur un ordinateur. 10 11. Dispositif comprenant une unité de traitement de données configurée pour mettre en oeuvre le procédé selon l'une des revendications 1 à 9. 15