FR 3054040 A1 20180119 FR 1656785 A 20160713 1 PROCEDE D'ESTIMATION D'UNE DEVIATION ANGULAIRE MOYENNE ENTRE L'AXE MAGNETIQUE ET UN AXE DE REFERENCE D'UN OBJET MAGNETIQUE DOMAINE TECHNIQUE [1] L'invention porte sur un procédé et un dispositif d'estimation d'une déviation angulaire moyenne entre l'axe magnétique du moment magnétique d'un objet magnétique et un axe de référence dudit objet magnétique. ÉTAT DE LA TECHNIQUE ANTÉRIEURE [2] Il est connu d'utiliser des objets magnétiques, notamment dans le cadre d'un système de relevé de la trace d'un crayon magnétique sur un support d'écriture. L'objet magnétique s'entend ici comme un objet auquel est associé un moment magnétique non nul, par exemple un aimant permanent fixé à un crayon non magnétique. [3] A titre d'exemple, le document W02014/053526 décrit un système de relevé de la trace d'un crayon auquel est fixé un aimant permanent annulaire. L'objet magnétique, ici l'aimant permanent, comporte un matériau magnétique, par exemple ferromagnétique ou ferrimagnétique, uniformément réparti autour d'un axe mécanique, dit axe de référence, qui correspond à son axe de révolution. L'aimant est conçu pour que son moment magnétique soit sensiblement colinéaire à l'axe de référence. [4] Le système de relevé de la trace du crayon muni de l'aimant permanent comporte un réseau de magnétomètres aptes à mesurer le champ magnétique généré par l'aimant permanent. Les magnétomètres sont fixés à un support d'écriture. [5] Cependant, le procédé de relevé de trace suppose que l'axe magnétique de l'aimant permanent, défini comme l'axe passant par le moment magnétique, soit effectivement colinéaire à l'axe de référence, ou présente une déviation angulaire acceptable entre l'axe magnétique et l'axe de référence. En effet, une déviation angulaire de quelques dixièmes de degrés peut entraîner une erreur sur le relevé de la trace qui peut alors présenter un manque de précision préjudiciable. Il peut alors être nécessaire de procéder à une 3054040 2 estimation préalable de la déviation angulaire entre l'axe magnétique et l'axe de référence de l'objet magnétique. EXPOSÉ DE L'INVENTION 5 [006] L'invention a pour objectif de proposer un procédé d'estimation d'une déviation angulaire moyenne entre un axe de référence d'un objet magnétique et un axe magnétique colinéaire à un moment magnétique dudit objet magnétique. Pour cela, le procédé d'estimation comporte les étapes suivantes : a) positionnement dudit objet magnétique vis-à-vis d'un plan de mesure défini par un 10 réseau de magnétomètres aptes à mesurer un champ magnétique en présence de l'objet magnétique ; b) rotation dudit objet magnétique autour dudit axe de référence ; c) mesure, au cours de la rotation, du champ magnétique, pour différents instants d'une durée de mesure, par lesdits magnétomètres ; 15 d) estimation d'un moment magnétique, dit instantané, de l'objet magnétique, auxdits différents instants de mesure, à partir du champ magnétique mesuré ; e) estimation de la déviation angulaire moyenne à partir desdits moments magnétiques instantanés. [7] Certains aspects préférés, mais non limitatifs de ce procédé d'estimation sont les 20 suivants. [8] L'étape d'estimation de la déviation angulaire moyenne peut comporter une sous-étape d'estimation d'un vecteur invariant en rotation autour de l'axe de référence à partir des moments magnétiques instantanés, l'estimation de la déviation angulaire moyenne étant effectuée, en outre, à partir dudit vecteur invariant. 25 [009] L'estimation de la déviation angulaire moyenne peut comporter le calcul d'une amplitude d'écart angulaire entre les moments magnétiques instantanés par rapport au vecteur invariant. [0010] Lors de l'étape de mesure du champ magnétique, l'objet magnétique peut réaliser un nombre entier de tours autour de l'axe de référence supérieur ou égal à 1. 3054040 3 [0011] Le procédé d'estimation peut comporter une étape d'homogénéisation d'une répartition angulaire des moments magnétiques instantanés autour de l'axe de référence au cours de la durée de mesure, cette étape d'homogénéisation comportant un calcul par interpolation d'une série temporelle dite homogénéisée de moments magnétiques 5 instantanés, à partir d'une série temporelle dite initiale de moments magnétiques instantanés obtenue préalablement, de sorte que les moments magnétiques instantanés successifs de la série temporelle homogénéisée présentent un écart angulaire sensiblement constant. [0012] L'étape d'homogénéisation peut comporter une phase itérative de calcul de 10 moments magnétiques interpolés dans laquelle un moment magnétique interpolé à une itération considérée est obtenu à partir d'un moment magnétique interpolé à une itération précédente, d'une valeur seuil prédéterminée et d'un vecteur unitaire défini à partir de deux moments magnétiques instantanés successifs de la série temporelle initiale. [0013] L'étape de positionnement peut comporter : 15 - une sous-étape d'estimation d'un premier paramètre de dispersion représentatif d'une dispersion de l'intensité des moments magnétiques instantanés de l'objet magnétique, pour différentes positions dites verticales de l'objet magnétique suivant un axe sensiblement orthogonal au plan de mesure ; - une sous-étape de positionnement de l'objet magnétique à une position verticale 20 pour laquelle la valeur du premier paramètre de dispersion est inférieure ou égale à une valeur seuil prédéterminée. [0014] Le procédé d'estimation peut comporter, préalablement à l'étape de positionnement, une étape de mesure, par lesdits magnétomètres, d'un champ magnétique ambiant en l'absence de l'objet magnétique, et peut comporter une sous- 25 étape de soustraction du champ magnétique ambiant aux champs magnétiques préalablement mesurés, de manière à obtenir le champ magnétique généré par l'objet magnétique. [0015] Le procédé d'estimation peut comporter une étape d'application d'un filtre passe-bas aux valeurs de composantes des moments magnétiques instantanés préalablement 30 estimés. 3054040 4 [0016] Le filtre passe-bas peut être une moyenne glissante sur un nombre donné d'échantillons. [0017] Le nombre d'échantillons peut être prédéterminé de sorte qu'un biais défini comme un écart entre une déviation angulaire dite réelle et la déviation angulaire moyenne 5 estimée soit minimal, en valeur absolue, pour une valeur dite seuil de déviation angulaire moyenne. [0018] Le procédé d'estimation peut comporter une étape de tri de l'objet magnétique selon l'écart de la valeur estimée de la déviation angulaire moyenne vis-à-vis d'une valeur dite seuil de déviation angulaire moyenne. 10 [0019] Le procédé d'estimation peut comporter une étape d'estimation d'un paramètre dit indicateur de qualité, à partir du calcul d'un deuxième paramètre de dispersion représentatif d'une dispersion des valeurs d'un paramètre dit rayon instantané calculé comme étant la norme d'un vecteur instantané défini entre chaque moment magnétique instantané et le vecteur invariant estimés préalablement. 15 [0020] Le procédé d'estimation peut comporter en outre une étape de comparaison du deuxième paramètre de dispersion à une valeur seuil prédéterminée. [0021] La prédétermination de la valeur seuil peut comporter les sous-étapes suivantes : - obtention, Q fois avec Q un entier naturel non nul, d'une distribution temporelle dudit rayon instantané, pour un même objet magnétique ; 20 - estimation, pour chaque distribution temporelle, d'une valeur moyenne dite rayon moyen définie comme la moyenne temporelle du rayon instantané, de manière à obtenir une distribution de Q rayons moyens estimés ; - estimation, à partir de la distribution des Q rayons moyens estimés, d'une valeur moyenne et un troisième paramètre de dispersion, dite valeur seuil, représentatif 25 d'une dispersion des valeurs des rayons moyens estimés. [0022] L'invention porte également sur un système d'estimation d'une déviation angulaire moyenne entre un axe magnétique et un axe de référence d'un objet magnétique, comportant : - un réseau de magnétomètres aptes à mesurer, à chaque instant de mesure d'une durée, un champ magnétique perturbé par l'objet magnétique selon au moins un 3054040 5 axe de mesure, les magnétomètres étant répartis spatialement de manière à définir un plan de mesure ; - un organe de maintien et de rotation, apte à maintenir l'objet magnétique dans une position vis-à-vis du plan de mesure, et à mettre en rotation l'objet magnétique 5 suivant son axe de référence ; - un processeur configuré pour mettre en oeuvre le procédé d'estimation selon l'une quelconque des caractéristiques précédentes. [0023] L'invention porte également sur un support d'enregistrement d'informations, comportant des instructions pour la mise en oeuvre du procédé d'estimation selon l'une 10 quelconque des caractéristiques précédentes, ces instructions étant aptes à être exécutées par un processeur. BRÈVE DESCRIPTION DES DESSINS [0024] D'autres aspects, buts, avantages et caractéristiques de l'invention apparaîtront 15 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 référence aux dessins annexés sur lesquels : la figure 1A est une vue schématique en perspective d'un dispositif d'estimation d'une déviation angulaire moyenne d'un objet magnétique, et la figure 1B illustre 20 schématiquement l'évolution temporelle du moment magnétique de l'objet magnétique au cours d'une durée T de mesure, formant ainsi un cercle qui s'étend autour de l'axe de référence ; la figure 2 est un organigramme d'un exemple de procédé d'estimation d'une déviation angulaire moyenne de l'objet magnétique ; 25 la figure 3A illustre un exemple d'évolution temporelle de l'écart angulaire élémentaire entre deux moments magnétiques instantanés successifs au cours de la durée T de mesure ; les figures 3B et 3C illustrent de manière schématique un exemple de répartition angulaire non homogène des moments magnétiques instantanés autour de l'axe de référence ; et les figures 3D et 3E illustrent de manière schématique un exemple de 3054040 6 répartition angulaire homogène, autour de l'axe de référence, des moments magnétiques instantanés obtenus après une étape d'homogénéisation angulaire ; la figure 4 est un organigramme d'un exemple partiel d'un procédé d'estimation d'une déviation angulaire moyenne comportant une étape d'homogénéisation angulaire ; 5 la figure 5A illustre un exemple schématique d'évolution temporelle des composantes du moment magnétique instantané et la figure 5B illustre un exemple d'évolution de la valeur d'un biais sur l'estimation de la déviation angulaire moyenne estimée, en fonction de la valeur de la déviation angulaire réelle, pour différentes valeurs K du nombre d'échantillons pris en compte lors d'une étape de filtrage ; 10 la figure 6 est un organigramme d'un exemple partiel d'un procédé d'estimation d'une déviation angulaire moyenne comportant une étape de filtrage des moments magnétiques instantanés ; la figure 7A illustre un exemple d'évolution de l'écart-type de la norme des moments magnétiques instantanés en fonction de la position verticale de l'objet magnétique vis-à- 15 vis du plan de mesure ; et la figure 7B est un organigramme d'un exemple partiel d'un procédé d'estimation d'une déviation angulaire moyenne comportant une étape de positionnement préalable de l'objet magnétique. EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS 20 [0025] Sur les figures et dans la suite de la description, les mêmes références représentent les éléments identiques ou similaires. De plus, les différents éléments ne sont pas représentés à l'échelle de manière à privilégier la clarté des figures. Par ailleurs, les différents modes de réalisation et variantes ne sont pas exclusifs les uns des autres et peuvent être combinés entre eux. Sauf indication contraire, les termes « sensiblement », 25 « environ », « de l'ordre de » signifient à 10% près, et de préférence à 5% près, voire à 1% près. 3054040 7 [0026] L'invention porte sur un dispositif et un procédé d'estimation de la déviation angulaire moyenne entre l'axe magnétique associé au moment magnétique d'un objet magnétique et un axe mécanique de référence de cet objet. [0027] L'objet magnétique comporte un matériau présentant un moment magnétique, par 5 exemple spontané. Il s'agit de préférence du matériau d'un aimant permanent. L'objet magnétique peut être un aimant permanent cylindrique, par exemple annulaire, tel qu'illustré dans le document W02014/053526 cité précédemment, auquel cas l'axe de référence correspond à l'axe de révolution de l'aimant. Il peut également s'agir d'un crayon équipé d'un tel aimant ou comportant un aimant permanent différent, par exemple intégré 10 au corps du crayon, auquel cas l'axe de référence peut correspondre à l'axe longitudinal suivant lequel s'étend le crayon passant par la pointe ou la mine d'écriture du crayon. Le terme crayon est à comprendre au sens large et peut englober les stylos, feutres, pinceaux ou tout autre organe d'écriture ou de dessin. [0028] Le matériau magnétique est de préférence ferrimagnétique ou ferromagnétique. Il 15 présente un moment magnétique spontané non nul même en l'absence d'un champ magnétique extérieur. Il peut présenter un champ magnétique coercitif supérieur à 100 A.m-1 ou 500 A.m-1 et son intensité est de préférence supérieure à 0,01 A.m2 voire à 0,1 A.m2. On considère par la suite que l'aimant permanent peut être approximé par un dipôle magnétique. L'axe magnétique de l'objet est défini comme étant l'axe colinéaire au 20 moment magnétique de l'objet. L'axe de référence correspond ici à l'axe de révolution de l'objet magnétique. L'écart angulaire entre l'axe de référence et l'axe magnétique est noté déviation angulaire. Dans la mesure où l'orientation de l'axe magnétique varie dans le temps au cours de la durée T de mesure, on nomme moment magnétique instantané le vecteur moment magnétique à un instant donné, et on nomme déviation angulaire 25 instantanée la valeur de la déviation angulaire à un instant donné. [0029] La figure 1 est une vue en perspective, schématique et partielle, d'un dispositif d'estimation 1 d'une déviation angulaire moyenne d'un objet magnétique 2, ici un aimant permanent cylindrique, par exemple annulaire. 3054040 8 [0030] Le dispositif d'estimation 1 est apte à estimer le moment magnétique instantané à différents instants de mesure au cours d'une durée T de mesure, dans un repère XYZ. Plus précisément, le dispositif 1 permet d'estimer la position de l'aimant permanent 2, et son moment magnétique, à différents instants, dans le repère XYZ. En d'autres termes, le 5 dispositif 1 permet de localiser la position et l'orientation de l'aimant permanent 2 à différents instants dans le repère XYZ. [0031] On définit ici et pour la suite de la description un repère direct tridimensionnel (X,Y,Z), où les axes X et Y forment un plan parallèle au plan de mesure du réseau de magnétomètres, et où l'axe Z est orienté de manière sensiblement orthogonale au plan de 10 mesure. Dans la suite de la description, les termes « vertical » et « verticalement » s'étendent comme étant relatifs à une orientation sensiblement parallèle à l'axe Z, et les termes « horizontal » et « horizontalement » comme étant relatifs à une orientation sensiblement parallèle au plan (X,Y). Par ailleurs, les termes « inférieur » et « supérieur » s'étendent comme étant relatifs à un positionnement croissant lorsqu'on s'éloigne du plan 15 de mesure suivant la direction +Z. [0032] La position Pd de l'aimant permanent 2 correspond aux coordonnées du centre géométrique de l'aimant 2. Le centre géométrique est le barycentre non pondéré de l'ensemble des points de l'aimant permanent 2. Le moment magnétique m de l'aimant 2 présente les composantes (mx, my, mz) dans le repère XYZ. Sa norme, ou intensité, est notée 20 II mll [0033] Le dispositif 1 comporte un réseau de magnétomètres M, répartis les uns vis-à-vis des autres de manière à former un plan de mesure Pmes. Le nombre de magnétomètres M, peut être, par exemple supérieur ou égal à 2, de préférence supérieur ou égal à 16, par exemple égal à 32, notamment lorsqu'il s'agit de magnétomètres tri-axes. Le réseau de 25 magnétomètres comporte cependant au moins 3 axes de mesure distants les uns des autres et non parallèles deux à deux. [0034] Les magnétomètres M, sont fixés à une plaque de protection 3 et peuvent être situés au niveau de la face arrière de la plaque 3, celle-ci étant réalisée dans un matériau non magnétique. Par fixés, on entend qu'ils sont assemblés à la plaque 3 sans aucun degré 30 de liberté. Ils sont ici alignés en lignes et colonnes, mais peuvent être positionnés mutuellement de manière sensiblement aléatoire. Les distances entre chaque 3054040 9 magnétomètre et ses voisins sont connues et constantes dans le temps. Par exemple, elles peuvent être comprises entre 1cm et 4cm. [0035] Les magnétomètres M, présentent chacun au moins un axe de mesure, par exemple trois axes, notés x,, y,, z,. Chaque magnétomètre mesure donc l'amplitude et la direction du 5 champ magnétique B perturbé par l'aimant permanent. Plus précisément, chaque magnétomètre M, mesure la norme de la projection orthogonale du champ magnétique B suivant les axes x,, y,, z, du magnétomètre. La sensibilité des magnétomètres M, peut être de 4.10-7T. Par champ magnétique perturbé B, on entend le champ magnétique ambiant Ba, c'est-à-dire non perturbé par l'aimant, auquel s'ajoute le champ magnétique Bd généré 10 par l'aimant. [0036] Le dispositif d'estimation 1 comporte en outre une unité de calcul 4 apte à calculer la position et l'orientation du moment magnétique de l'aimant 2 dans le repère XYZ à partir des mesures des magnétomètres M. Il permet en outre de déterminer la déviation angulaire moyenne de l'aimant permanent 2 à partir des mesures du moment magnétique. 15 [0037] Pour cela, chaque magnétomètre M, est électriquement connecté à l'unité de calcul 4 par un bus de transmission d'informations (non représenté). L'unité de calcul 4 comporte un processeur programmable 5 apte à exécuter des instructions enregistrées sur un support d'enregistrement d'informations. Il comporte en outre une mémoire 6 contenant les instructions nécessaires pour la mise en oeuvre de certaines étapes d'un procédé 20 d'estimation de la déviation angulaire moyenne par le processeur 5. La mémoire 6 est également adaptée à stocker les informations calculées à chaque instant de mesure. [0038] L'unité de calcul 4 implémente un modèle mathématique associant la position de l'aimant permanent 2 dans le repère XYZ, ainsi que l'orientation et l'intensité du moment magnétique, aux mesures des magnétomètres M. Ce modèle mathématique est construit 25 à partir des équations de l'électromagnétisme, en particulier de la magnétostatique, et est paramétré notamment par les positions et orientations des magnétomètres dans le repère XYZ. [0039] De préférence, pour pouvoir approximer l'aimant permanent 2 à un dipôle magnétique, la distance entre l'aimant permanent 2 et chaque magnétomètre M, est 3054040 10 supérieure à 2, voire 3 fois la plus grande dimension de l'aimant permanent 2. Cette dimension peut être inférieure à 20cm, voire inférieure à 10cm, voire à 5cm. [0040] Le dispositif d'estimation 1 comporte également un organe 7 de maintien et de rotation, apte à maintenir l'aimant permanent 2 vis-à-vis du plan de mesure, dans une 5 position quelconque vis-à-vis du plan de mesure. L'aimant permanent peut être situé au- dessus, c'est-à-dire en regard, du réseau de magnétomètres M,, dans une position verticale suivant l'axe Z constante pendant la durée T de mesure. Il est en outre apte à assurer la rotation du l'aimant 2 suivant son axe de référence Aref pendant la durée de mesure. Aussi, pendant la durée de mesure, la position verticale de l'aimant permanent 2 est fixe et seule 10 sa position angulaire varie du fait de la rotation suivant l'axe de référence Aref. L'organe 7 de maintien et de rotation comporte un moteur 8 associé à un bras 9. Le bras 9 peut ainsi recevoir et maintenir l'aimant permanent 2 et assure la rotation de ce dernier. Le bras 9 est réalisé dans un matériau non magnétique et le moteur 8 est suffisamment éloigné du plan de mesure Pmes et de l'aimant permanent 2 pour ne pas induire de perturbation du 15 champ magnétique mesuré. [0041] La figure 2 est un organigramme d'un exemple de procédé d'estimation de la déviation angulaire moyenne que présente un objet magnétique, le procédé étant mis en oeuvre par le dispositif d'estimation de la figure 1A. 20 [0042] Le procédé comporte une étape 100 de mesure préalable du champ magnétique ambiant Ba, c'est-à-dire ici le champ magnétique non perturbé par la présence de l'aimant 2. Pour cela, chacun des magnétomètres M, mesure le champ magnétique Bia en l'absence de l'aimant 2, c'est-à-dire la projection du champ magnétique Ba sur chaque axe d'acquisition x,, y,, z, des différents magnétomètres M. 25 [0043] Le procédé comporte ensuite une étape 110 de mesure du champ magnétique perturbé par l'aimant 2 en rotation autour de son axe de référence Aref au cours d'une durée T de mesure. [0044] Lors d'une sous-étape 111, l'aimant permanent 2 est positionné vis-à-vis du plan de mesure Pmes par l'organe 7 de maintien et de rotation, qui définit l'orientation de l'axe de 30 référence Aref de l'aimant permanent 2 vis-à-vis du plan de mesure Pmes dans le repère XYZ. 3054040 11 Dans cet exemple, l'aimant permanent 2 est positionné au-dessus du plan de mesure Pmes mais tout autre positionnement est possible. L'orientation de l'axe Aref peut être quelconque, mais peut être avantageusement sensiblement orthogonale au plan de mesure Pmes. Le moment magnétique m de l'aimant permanent 2 est non colinéaire à l'axe 5 de référence Aref et forme, vis-à-vis de cet axe, une déviation angulaire a à déterminer. [0045] Lors d'une sous-étape 112, l'aimant permanent 2 est mis en rotation par l'organe 7 autour de l'axe de référence Aref. L'axe de référence Aref est statique au cours de la durée T de mesure, autrement dit sa position et son orientation dans le repère XYZ ne varient pas lors de la durée T. La vitesse de rotation de l'aimant permanent 2, notée Ô, est de 10 préférence constante au cours de la durée T. [0046] Au cours de la rotation, pendant la durée T=[ti; tN], chaque magnétomètre M, mesure le champ magnétique Bei) perturbé par la présence de l'aimant permanent 2, à une fréquence d'échantillonnage fe=N/T. On obtient ainsi N instants tj de mesure. Chaque magnétomètre M,, à l'instant tj, mesure la projection du champ magnétique B suivant le ou 15 les axes d'acquisition x,, y,, z,. On obtient ainsi une série temporelle {Bi(ti)}N de N mesures du champ magnétique perturbé. [0047] Le procédé comporte ensuite une étape 120 d'estimation du moment magnétique instantané m(tj) de l'aimant permanent 2, pour chaque instant t de mesure, à partir des mesures du champ magnétique perturbé Bei). 20 [0048] Lors d'une sous-étape 121, l'unité de calcul 4 déduit le champ magnétique généré Bid(ti) par l'aimant permanent 2, pour chaque magnétomètre M, et à chaque instant tj, à partir des mesures du champ magnétique ambiant Bi' et du champ magnétique perturbé Bei). Pour cela, on calcule Bid(tj) = Bei) - Bi'. [0049] Lors d'une sous-étape 122, l'unité de calcul 4 estime la position Pd de l'aimant 25 permanent 2 ainsi que son moment magnétique instantané m(tj), à chaque instant tj, à partir du champ magnétique Bid(tj) préalablement calculé. Pour cela, l'unité de calcul 4 résout un modèle mathématique d'équations de l'électromagnétisme associant la position et le moment magnétique de l'aimant permanent 2 au champ magnétique qu'il génère Bid(tj) = f(Pd ; m(tj)). Ainsi, à chaque instant de mesure, l'unité de calcul 4 détermine les 30 coordonnées de la position de l'aimant permanent 2 ainsi que les composantes mx, my, mz 3054040 12 du moment magnétique instantané, dans le référentiel XYZ. On obtient ainsi une série temporelle {m(tJ)}N de N vecteurs instantanés du moment magnétique m(tj). [0050] Comme le montrent les figures 1A et 1B, la série temporelle {m(tJ)}N forme un cercle C qui s'étend autour de l'axe de référence Aref. Au cercle C peuvent être associés des rayons 5 instantanés R(tj) relatifs aux moments magnétiques instantanés m(tj), ainsi qu'un rayon moyen R dont la valeur dépend de la déviation angulaire moyenne a à déterminer. [0051] Le procédé comporte ensuite une étape 130 d'estimation de la déviation angulaire moyenne a formée par l'écart angulaire entre le moment magnétique m du dipôle magnétique et l'axe de référence Aref. 10 [0052] Lors d'une sous-étape 131, l'unité de calcul estime un vecteur ma invariant lors de la rotation de l'aimant permanent 2. Pour cela, on réalise la moyenne temporelle, sur les N instants tj de mesure, de chaque coordonnée mx(tj), my(tj), mz(ti) du moment magnétique instantané m(tj) dans le repère XYZ. Ainsi, ma = ( N ; N ; N) = N. L'opérateur N correspond ici à la moyenne temporelle arithmétique, éventuellement 15 pondérée, sur les N instants de mesure. On obtient ainsi un vecteur ma sensiblement colinéaire à l'axe de référence Aref qui s'étend entre la position Pd de l'aimant 2 et la position Po du centre du cercle C. [0053] Lors d'une sous-étape 132, l'unité de calcul estime un paramètre représentatif d'une amplitude moyenne de l'écart angulaire des moments magnétiques instantanés par 20 rapport à l'axe de référence. Ici, ce paramètre est le rayon moyen R du cercle C, tel que = Ilin(t i) - m oui > N , c'est-à-dire la moyenne temporelle de la norme des vecteurs instantanés m(tj)-mo. [0054] Lors d'une sous-étape 133, l'unité détermine la déviation angulaire moyenne â à partir du rayon moyen estimé R et, dans cet exemple, du vecteur invariant mo. Pour cela, 25 il calcule l'arc-tangente du rapport du rayon moyen R sur la norme du vecteur invariant ma : â = tan-l(R/Ilmo II). On obtient ainsi une estimation de la déviation angulaire entre l'axe de référence de l'objet magnétique et son axe magnétique. [0055] Bien entendu, d'autres calculs équivalents peuvent être effectués. Ainsi, en variante, lors de la sous-étape 133, l'estimation de la déviation angulaire moyenne peut 30 être obtenue à partir de l'arc-sinus du rapport entre le rayon moyen R estimé 3054040 13 préalablement et la moyenne temporelle de la norme du moment magnétique instantané N sur les N instants tj de mesure. [0056] Il est également possible, en variante, d'estimer, lors de la sous-étape 132, le rayon instantané R(tj) du cercle à chaque instant tj de mesure et de calculer lors de la sous-étape 5 133, la déviation angulaire instantanée cc(ti) correspondante puis de déterminer la déviation angulaire moyenne â comme la moyenne temporelle N des déviations angulaires instantanées sur les N instants tj de mesure. [0057] Il est également possible, en variante, d'estimer, lors de la sous-étape 132, le rayon moyen R par ajustement d'au moins un paramètre d'une équation d'un cercle dans le 10 repère XYZ parcourant la série temporelle des moments magnétiques instantanés, par exemple la série {m(tJ)}N. Ainsi, l'équation du cercle peut être (mx(ti)-P0,x)2- (my(ti)-Po,y)2- (mz(ti)-P0,z)2 = R, où Po,x, Po,y, Po,z sont les coordonnées du centre Po du cercle C déduites à partir de la position Pd du moment magnétique estimée lors de la sous-étape 122 et du vecteur invariant mo estimé lors de la sous-étape 131, et où mx(tj), my(tj), mz(tj) sont les 15 valeurs des composantes du moment magnétique instantané de la série temporelle {m(tJ)}N considérée. L'ajustement de la valeur du rayon moyen R peut être réalisé de manière connue, par des méthodes classiques de régression, par exemple multilinéaire ou polynomiale. [0058] A titre d'exemple, l'objet magnétique est un aimant permanent à symétrie de 20 révolution autour de son axe de référence Aref. Il est positionné par exemple à 4 cm à la verticale d'un réseau de 32 magnétomètres tri-axes d'une sensibilité de 4.10-7T espacés par exemple de 4 cm. La vitesse de rotation est de ni5 radis (1 tour en 10s) et la fréquence d'échantillonnage est de 140Hz. Lors d'une durée T de 10s, l'objet magnétique réalise un unique tour de rotation et le dispositif de détermination acquiert 1400 instants de mesure. 25 L'unité de calcul estime à 0,165 A.m2 l'intensité moyenne du moment magnétique pour un bruit RMS (Roof Mean Square, en anglais) de 5.10-4 A.m2. Le dispositif de détermination estime à 0,01° la déviation angulaire moyenne entre l'axe de référence de l'aimant et son axe magnétique. [0059] De préférence, la durée T de mesure et la vitesse de rotation Ô sont choisies de 30 sorte que, pendant la durée T, l'objet magnétique ait effectué un nombre entier de tours 3054040 14 complets, par exemple un unique tour. De plus, dans le but d'obtenir une répartition angulaire sensiblement homogène des données de mesure autour de l'axe de référence Aref au cours de la durée T, la rotation élémentaire 4ej entre deux instants de mesure tj et ti+1 est sensiblement constante. 5 [0060] Comme détaillé maintenant, dans le cas où les données de mesure ne présentent pas une répartition angulaire sensiblement homogène autour de l'axe de référence Aref pendant la durée T, une étape d'homogénéisation angulaire peut être mise en oeuvre. [0061] Les figures 3A à 3E et la figure 4 illustrent une étape 140 d'homogénéisation de la 10 répartition angulaire des moments magnétiques instantanés autour de l'axe de référence Aref au cours de la durée T de mesure. Cette étape repose sur le ré-échantillonnage des moments magnétiques instantanés au cours de la durée T de mesure par interpolation, ici au moyen d'une méthode d'interpolation linéaire, de telle sorte que les écarts angulaires entre deux moments magnétiques instantanés successifs soient sensiblement constants. 15 En d'autres termes, on ré-échantillonne une série temporelle initiale de N moments magnétiques instantanés qui présentent une répartition angulaire inhomogène autour de l'axe de référence Aref au cours de la durée T, par exemple la série temporelle {m(tJ)}N estimée par exemple lors de la sous-étape 122, pour obtenir une série temporelle dite homogénéisée {mh(ti)}N, de N' moments magnétiques instantanés qui présentent une 20 répartition angulaire sensiblement homogène. [0062] Cette étape est avantageuse lorsque, comme l'illustrent les figures 3A à 3C, la répartition angulaire des moments magnétiques instantanés autour de l'axe de référence Aref n'est pas homogène au cours de la durée T. Cela se traduit par le fait que les écarts angulaires 413j définis entre deux moments magnétiques instantanés m(tj) et m(41), à des 25 instants de mesure tj et ti.,1 successifs, présentent une forte dispersion autour de la valeur moyenne d/3. [0063] Dans ce cas de figure, le vecteur mo, calculé lors de la sous-étape 131 décrite ci-dessus, peut ne pas être invariant en rotation, c'est-à-dire peut ne pas être colinéaire à l'axe de référence Aref. Cette erreur sur le calcul du vecteur mo peut introduire un biais sur 3054040 15 la valeur estimée de la déviation angulaire moyenne a. Cette étape 140 d'homogénéisation permet ainsi de réduire l'erreur sur l'estimation du vecteur invariant mo. [0064] En référence à la figure 4, cette étape 140 d'homogénéisation peut comporter une sous-étape 141 facultative dans laquelle on estime tout d'abord les écarts angulaires APJ 5 définis par les moments magnétiques estimés m(tj) et m(ti.n) à des instants de mesure successifs de la série temporelle initiale, par exemple ici {m(tJ)}N. Pour cela, on peut calculer l'écart angulaire APJ à partir de l'arc sinus de la norme du produit vectoriel de ces deux vecteurs : zl,8 = Ilm(ti) x m(ti+1)11/(11m(ti)11.11m(ti+1)11) où le symbole « x » représente le produit vectoriel et le symbole « . » le produit de scalaires. 10 [0065] Ensuite, lors d'une sous-étape 142 également facultative, l'estimation de la dispersion des écarts angulaires APJ peut être obtenue en calculant l'écart-type bop, ou un paramètre équivalent, sur les écarts angulaires APJ estimés. Si la valeur de l'écart-type bop est inférieure à un seuil prédéterminé, les données de mesure présentent une répartition angulaire sensiblement homogène. Dans le cas contraire, on procède à l'homogénéisation 15 de la répartition angulaire des moments magnétiques instantanés autour de l'axe de référence Aref au cours de la durée T de mesure. [0066] Lors d'une sous-étape 143, on ré-échantillonne la série temporelle initiale {m(tJ)}N par interpolation pour obtenir une série temporelle dite homogénéisée {mh(tk)}N, de N' moments magnétiques interpolés, dont les écarts angulaires A130 sont sensiblement 20 constants. Il s'agit d'un sur-échantillonnage lorsque N'>N mais le nombre N' peut être inférieur ou égal à N. [0067] Pour cela, une approche illustrée sur la figure 3D consiste à discrétiser l'écart angulaire APJ entre deux moments magnétiques instantanés successifs m(tj) et m(41) de la série temporelle initiale {m(tJ)}N lorsque sa valeur est supérieure à une valeur seuil 25 prédéterminée A130, par exemple par interpolation linéaire (comme représenté sur la figure 3D) ou par un autre d'interpolation, par exemple polynomiale. [0068] Comme l'illustre la figure 3C, on considère deux moments magnétiques m(ti,i) et m(tj) mesurés aux instants successifs tj et tpn, dont l'écart angulaire APJ présente ici une valeur supérieure au seuil A(30. On détermine un vecteur unitaire ej,j+Mm(ti+1)- 30 m(tj))/II m(41)-m(ti)11. Puis, on calcule un voire plusieurs moments magnétiques interpolés 3054040 16 tels que : = mj,k + 41:30.ej,j+1, avec ik,k.0 = m(tj). L'itération sur k se poursuit tant que l'écart angulaire entre m(41) et fiii,k+1 est supérieur à 4130. A titre illustratif, on considère les moments magnétiques m(t1) et m(t2) mesurés aux instants successifs ti et t2, pour lesquels l'écart angulaire 4131 présente ici une valeur supérieure au seuil 4130. On détermine 5 le vecteur unitaire el,2=(m(t2)-m(t1))/11 m(t2)-m(t1)II . Puis, on calcule un voire plusieurs moments magnétiques interpolés tels que : fi1,k+1= m1,k + 4130.e1,2, avec = m(tJ=1) et tant que l'écart angulaire entre m(t2) et fi1,k+1 est supérieur à 4130. Cette sous-étape 143 peut être appliquée à tous les couples de moments magnétiques m(tj) et m(ti,i) mesurés à des instants successifs, pour lesquels 413j>4130. A titre d'exemple, sur la figure 3D, trois 10 vecteurs interpolés 71.21, mzz, 17123 sont calculés entre les moments magnétiques instantanés m(t2) et m(t3) de telle sorte que les écarts angulaires entre ces moments magnétiques m(t2), 7'11'21, fii22, fh-'23 et m(t3) sont sensiblement constants et égaux à 4130. [0069] On obtient ainsi une série temporelle homogénéisée imilk de N' moments magnétiques instantanés imilk correspondant aux N moments magnétiques de la série 15 temporelle initiale {m(tJ)}N estimés lors de la sous-étape 122 auxquels s'ajoutent les moments magnétiques interpolés m. [0070] En variante (non représentée), on peut construire une série temporelle homogénéisée {mh}Ny sans avoir à y incorporer, comme précédemment, la série temporelle initiale {m(tJ)}N. Ainsi, on peut calculer N' moments magnétiques interpolés tels que : fit k+1 20 = 4130.eitn, avec fiik.0 = m(t J.1) et avec le vecteur unitaire ej,j+1=(m(ti+1)-m(ti))/11m(ti+1)- m(ti)11. L'itération sur k se poursuit tant que l'écart angulaire entre fitk+1 et m(ti-n) est supérieur à 4130, et tant que le produit du nombre k d'itérations et de 4130 est inférieur à une ou plusieurs fois 27c. De plus, le vecteur unitaire ejel est défini en j et j+1 de sorte que les écarts angulaires entre m(tj) et fit k+1 d'une part, et entre fitk+1 et m(ti,i) d'autre part, 25 sont inférieurs à l'écart angulaire 413j entre m(tj) et m(ti+1). On obtient ainsi la série temporelle homogénéisée imilk de N' moments magnétiques interpolés. [0071] Au préalable de la sous-étape 143 de calcul de la série temporelle homogénéisée {le}N-, il est possible d'effectuer un sur-échantillonnage de la série temporelle initiale, par exemple {m(tJ)}N, de préférence par interpolation polynomiale, éventuellement à l'aide de 30 splines, pour ainsi obtenir une nouvelle série temporelle à homogénéiser au moyen de la sous-étape 143 décrite précédemment. 3054040 17 [0072] Ainsi, l'écart-type calculé sur les écarts angulaires 413' associés à la série temporelle homogénéisée {W}N, est alors minimal ou sensiblement nul dans la mesure où ils sont sensiblement égaux à la valeur 4130. Les moments magnétiques instantanés mh(tk) présentent alors une répartition angulaire sensiblement homogène autour de l'axe de 5 référence Aref au cours de la durée T de mesure. [0073] La série temporelle homogénéisée {mh}w correspond à un ré-échantillonnage des données de mesure à une nouvelle fréquence fe'= N'/T, de sorte que chaque moment magnétique instantané mh(tk) de la série temporelle homogénéisée {W}N, peut être vue comme une estimation du moment magnétique à différents instants tk de mesure, avec 10 k={1, 2, ..., N'}. [0074] Le procédé d'estimation de la déviation angulaire moyenne â entre l'axe magnétique et l'axe de référence est ensuite poursuivi avec l'étape 130 décrite précédemment. Lors de la sous-étape 131, le vecteur invariant ma est déduit à partir des N' moments magnétiques instantanés mh(tk) de la série temporelle homogénéisée {mil}w 15 et non plus à partir des N moments magnétiques instantanés m(tj) de la série temporelle {m}N obtenue lors de la sous-étape 122. De plus, lors de la sous-étape 132, le rayon moyen du cercle C est calculé à partir des moments magnétiques mh(tk) de la série temporelle homogénéisée et/ou à partir des moments magnétiques fit(ti) de la série temporelle filtrée {fit(ti)}N décrite plus loin, et non pas à partir des moments magnétiques m(tj) de la série 20 temporelle obtenue lors de la sous-étape 122. [0075] Comme détaillé maintenant, il peut être avantageux d'effectuer une étape 150 de filtrage passe-bas des valeurs des composantes des moments magnétiques instantanés dans le repère XYZ, dans le but de diminuer la dispersion temporelle, autrement dit le bruit 25 de mesure, qu'elles peuvent présenter, par exemple à l'aide d'une moyenne mobile arithmétique ou exponentielle calculée sur un nombre K d'échantillons. C'est particulièrement avantageux lorsque la déviation angulaire moyenne est inférieure à 0,2° et/ou lorsque le bruit associé aux composantes des moments magnétiques est du même ordre de grandeur que la valeur du rayon moyen R du cercle C. 3054040 18 [0076] La figure 5A est un diagramme qui illustre de manière schématique la dispersion temporelle des valeurs des composantes m(tj) des moments magnétiques instantanés m(tj) dans le repère XYZ, ici issus de la série temporelle {m(tJ)}N obtenue lors de la sous-étape 122. 5 [0077] Cette dispersion des valeurs des composantes m(ti) au cours du temps peut induire un biais dans l'estimation de la déviation angulaire moyenne â. Le biais est ici noté b et est défini comme la différence entre la déviation angulaire réelle cc, de l'aimant permanent 2 et la déviation angulaire moyenne â, soit b = â - [0078] Dans le but de diminuer la dispersion temporelle des valeurs des composantes m(tj), 10 l'étape 150 de filtrage passe-bas par une moyenne glissante, ou moyenne mobile, peut être effectuée. [0079] Lors d'une sous-étape 151 facultative, l'unité calcule l'écart-type cy(Il m(ti)11 ) sur les N valeurs de la norme II m(ti)11 des moments magnétiques instantanés, ou un paramètre équivalent représentatif de la dispersion des valeurs de la norme m(ti)11 des moments 15 magnétiques instantanés. Si la valeur de l'écart-type cy(Il m(ti)11 ) est supérieure à une valeur seuil prédéterminée, on procède à un filtrage des 3N valeurs de composantes m(tj). [0080] Lors d'une sous-étape 152, l'unité applique un filtre HK(tj) de type passe-bas, ici par moyenne glissante sur K échantillons, par exemple une moyenne arithmétique ou exponentielle, ou un type de filtre équivalent, sur la série temporelle de 3N valeurs {m(tJ)}N 20 des composantes des moments magnétiques instantanés dans le repère XYZ. On a donc une nouvelle série temporelle telle que {fit(ti)} = Hdtj).{m(ti)}, où fit(ti) est les valeurs filtrées des composantes du moment magnétique instantané dans le repère XYZ, à l'instant tj. On obtient ainsi une série temporelle dite filtrée {fit(ti)}N de moments magnétiques instantanés qui est ensuite prise en compte pour l'exécution de l'étape 130 d'estimation 25 de la déviation angulaire moyenne â. [0081] Bien entendu, cette étape 150 de filtrage passe-bas des composantes des moments magnétiques instantanés peut être appliquée à la série temporelle {m(tJ)}N obtenue lors de la sous-étape 122, comme à la série temporelle homogénéisée {mh(tk)}N-. [0082] De préférence, le nombre K d'échantillons est choisi de sorte que le produit de K 30 avec l'écart angulaire moyen de la série temporelle considérée, par exemple la série 3054040 19 {m(tJ)}N issue de la sous-étape 122, est inférieur à une valeur donnée, par exemple 45°, et de préférence 10°, voire 5° et de préférence 1°. Alternativement, l'écart angulaire moyen peut être obtenu par la vitesse de rotation de l'organe de rotation multipliée par la fréquence d'échantillonnage. On évite ainsi d'appliquer un filtrage trop contraignant à la 5 série temporelle considérée qui risquerait de masquer la dynamique du signal de rotation. [0083] Il est avantageux d'effectuer l'étape 140 d'homogénéisation au préalable de l'étape 131 d'estimation du vecteur invariant mo. Par ailleurs, il est avantageux d'effectuer l'étape 150 de filtrage au préalable de l'étape 132 d'estimation du rayon moyen R du cercle C. Il est enfin avantageux d'effectuer l'étape 150 de filtrage avant l'étape 140 10 d'homogénéisation pour que la série temporelle prise en compte dans l'étape 140 d'homogénéisation soit moins bruitée. [0084] Comme l'illustre schématiquement la figure 5B, la valeur du biais b = â - cc, peut varier en fonction de la valeur de la déviation angulaire réelle ar. Elle est illustrée sur la 15 figure 5B par une fonction affine, mais elle peut présenter un autre type de fonction. Ainsi, il apparaît que la valeur du biais est positive pour les faibles valeurs de déviation angulaire et négative pour les fortes valeurs. Il apparaît également que la valeur de la déviation angulaire réelle pour laquelle la valeur du biais b est sensiblement nulle dépend notamment du nombre K d'échantillons pris en compte dans la moyenne glissante, 20 notamment dans le cas d'une moyenne glissante arithmétique. [0085] Lorsque le procédé d'estimation de la déviation angulaire moyenne â comporte une sous-étape supplémentaire de classification de l'aimant permanent 2 en vue d'opérer un tri en fonction d'une valeur de référence luth de la déviation angulaire, il peut être avantageux de déterminer le nombre K d'échantillons à prendre en compte dans la 25 moyenne glissante pour que la valeur du biais soit sensiblement nulle à la valeur de référence luth. A titre purement illustratif, pour une valeur de référence ad-, égale à 0,5° et pour une série temporelle de 1400 moments magnétiques instantanés dont l'intensité moyenne est de 0,165 A.m2 bruitée par un signal RMS de 5.10-4 A.m2, une moyenne glissante sur K=10 échantillons peut convenir. On diminue, voire annule ainsi la valeur du 30 biais pour la valeur de référence luth. 3054040 20 [0086] La sous-étape de classification comporte la comparaison de la valeur estimée de la déviation angulaire moyenne â avec la valeur prédéterminée de référence luth. Ainsi, on peut opérer une classification et donc un tri des objets magnétiques selon que la déviation angulaire moyenne â est supérieure ou inférieure à la valeur de référence luth. 5 [0087] A titre illustratif, pour une valeur de référence de 0,5°, on peut considérer que l'aimant permanent 2 est utilisable lorsque sa déviation angulaire estimée est inférieure ou égale à la valeur de référence luth et qu'il est à écarter lorsqu'elle est supérieure à la valeur de référence. En effet, pour une valeur estimée â inférieure à la valeur seuil luth, le biais est positif, c'est-à-dire que la valeur réelle cc, est inférieure à la valeur estimée. L'aimant 10 permanent est alors utilisable. Par contre, pour une valeur estimée â supérieure à la valeur seuil luth, le biais est négatif, c'est-à-dire que la valeur réelle cc, est supérieure à la valeur estimée. L'aimant permanent est alors à écarter. Bien entendu, la valeur de référence et la décision de tri dépendent des applications prévues pour l'aimant permanent. 15 [0088] Par ailleurs, le procédé d'estimation de la déviation angulaire peut comporter une étape 90 préalable de positionnement vertical de l'aimant permanent 2 vis-à-vis du plan de mesure Pmes, décrite en référence aux figures 7A et 7B. [0089] Cette étape de positionnement comporte une sous-étape 91 d'estimation d'un paramètre représentatif de la dispersion des valeurs de l'intensité du moment magnétique 20 instantané, pour différentes positions verticales suivant l'axe Z au-dessus du plan de mesure Pmes. [0090] Pour cela, les étapes 110 et 120 sont effectuées pour différentes valeurs de positionnement vertical de l'objet suivant l'axe Z. Pour chaque position verticale, on obtient ainsi une série temporelle de moments magnétiques instantanés, par exemple la 25 série temporelle {m(tJ)}N obtenue lors de la sous-étape 122, et on déduit par exemple l'écart-type cy(Il m(ti)11 ) sur les N valeurs de l'intensité m(ti)I1 des moments magnétiques instantanés. En référence à la figure 7A, on obtient ainsi une série spatiale qui donne la valeur de l'écart-type cy(Il m(ti)II ) pour différentes positions verticales Z occupées par l'objet magnétique. 3054040 21 [0091] Lors d'une sous-étape 92, on sélectionne une valeur de position verticale Zp pour laquelle la valeur de l'écart-type cy(Il m(ti)11) sur l'intensité du moment magnétique instantané est inférieure à une valeur seuil are. Puis, on positionne l'objet magnétique de sorte qu'il occupe la position verticale Zp. 5 [0092] Le procédé d'estimation de la déviation angulaire moyenne â est ensuite poursuivi avec l'étape 100 et les étapes 120 et 130 décrites précédemment, éventuellement complétées par les étapes 140 et 150. [0093] Par ailleurs, le procédé peut comporter une étape supplémentaire d'estimation 10 d'un indicateur de qualité des mesures effectuées pendant la durée T de mesure. Cette étape peut avoir lieu, de préférence, après l'exécution de l'étape 120 ou de l'étape 150. [0094] Pour cela, l'unité de calcul 4 détermine l'écart-type 6R sur les N valeurs du rayon instantané qui est estimé à partir des moments magnétiques instantanés m(tj). Le rayon instantané peut être calculé par la relation R(tj) = Il m(tj)-moll . Il est calculé ici à partir de la 15 série temporelle {m(tJ)}N obtenue lors de la sous-étape 122. L'indicateur de qualité est ici la grandeur ap/VN qui permet de caractériser la précision de la série de mesures. D'autres indicateur de qualité peuvent être utilisés à partir de l'écart-type 6R, par exemple asin( (6R/VN)/( N) ) ou atan( (6p/VN )I(IImoII) ). [0095] Dans le but de détecter une éventuellement anomalie ayant eu lieu lors de la série 20 de mesure, l'écart-type 6R sur les N valeurs du rayon instantané peut être comparé à une valeur de référence cy(R) préalablement déterminée. Cette valeur cy(R) peut être déterminée en effectuant un grand nombre Q d'estimations du rayon moyen R pour un même type d'aimant permanent 2. Le nombre Q peut être supérieur à 100 voire supérieur ou égal à 1000. Le rayon moyen R est par exemple estimé en exécutant les étapes 120 et 25 130, éventuellement complétées par les étapes 140 et 150. Ainsi, on obtient Q fois une distribution temporelle {R(ti)}N dudit rayon instantané R(tj), pour un même objet magnétique. Q est un entier naturel non nul, de préférence supérieur ou égal à 50, voire à 100. Puis on estime, pour chacune des Q distributions temporelles {R(ti)}N, une valeur moyenne dite rayon moyen R défini comme la moyenne temporelle du rayon instantané 30 N. On obtient ainsi une distribution {R}Q de Q rayons moyens estimés R. A partir de 3054040 22 cette distribution, on estime la valeur moyenne 0 ainsi qu'un paramètre d'estimation de la dispersion des valeurs des rayons moyens estimés R autour de la moyenne o, telle que l'écart-type cy(R). On compare la grandeur 6R à cy(R), par exemple en calculant le rapport 6R/6(R). Lorsque ce rapport est supérieur ou égal à une valeur seuil, par exemple 5 1, 2 voire 3, on peut estimer que les mesures effectuées pendant la durée T ne sont pas représentatives et ne doivent donc pas être prises en compte. Cette erreur liée aux mesures peut être issue d'un problème mécanique, par exemple un problème de stabilité de l'axe de référence, voire un problème de détection magnétique, par exemple une magnétisation excessive d'un ou plusieurs magnétomètres. 10 REVENDICATIONS1. Procédé d'estimation d'une déviation angulaire moyenne () entre un axe de référence (Are) d'un objet magnétique (2) et un axe magnétique colinéaire à un moment magnétique (m) dudit objet magnétique (2), comportant les étapes suivantes : a) positionnement (90 ; 110) dudit objet magnétique (2) vis-à-vis d'un plan de mesure (Pmes) défini par un réseau de magnétomètres (M,) aptes à mesurer un champ magnétique (B) en présence de l'objet magnétique (2) ; b)rotation (110) dudit objet magnétique (2) autour dudit axe de référence (Aref) ; c) mesure (110), au cours de la rotation, du champ magnétique (Bei)), pour différents instants (ti) d'une durée (T) de mesure, par lesdits magnétomètres (M,) ; d)estimation (120) d'un moment magnétique, dit instantané, (m(tj)) de l'objet magnétique (2), auxdits différents instants (ti) de mesure, à partir du champ magnétique mesuré (Bei)) ; e)estimation (130) de la déviation angulaire moyenne () à partir desdits moments magnétiques instantanés (m(tj)). 2. Procédé selon la revendication 1, dans lequel l'étape d'estimation (130) de la déviation angulaire moyenne () comporte une sous-étape d'estimation (131) d'un vecteur (ma) invariant en rotation autour de l'axe de référence (Aref) à partir des moments magnétiques instantanés (m(tj)), l'estimation (130) de la déviation angulaire moyenne () étant effectuée, en outre, à partir dudit vecteur invariant (ma). 3. Procédé selon la revendication 1 ou 2, dans lequel l'estimation (130) de la déviation angulaire moyenne () comporte le calcul d'une amplitude (R) d'écart angulaire entre les moments magnétiques instantanés (m(tj)) par rapport au vecteur invariant (ma). 4. Procédé selon l'une quelconque des revendications 1 à 3, dans lequel, lors de l'étape de mesure (110) du champ magnétique, l'objet magnétique (2) réalise un nombre entier de tours autour de l'axe de référence (Aref) supérieur ou égal à 1. 5. Procédé selon l'une quelconque des revendications 1 à 4, comportant une étape (140) d'homogénéisation d'une répartition angulaire des moments magnétiques 3054040 24 instantanés autour de l'axe de référence (Are) au cours de la durée (T) de mesure, comportant un calcul par interpolation d'une série temporelle dite homogénéisée ({mh(tk)}N,) de moments magnétiques instantanés, à partir d'une série temporelle dite initiale ({m(tj)}N) de moments magnétiques instantanés obtenue préalablement, de sorte 5 que les moments magnétiques instantanés successifs de la série temporelle homogénéisée ({mh(tk)},w) présentent un écart angulaire sensiblement constant. 6. Procédé selon la revendication 5, dans lequel l'étape d'homogénéisation comporte une phase itérative de calcul de moments magnétiques interpolés dans laquelle un moment magnétique interpolé à une itération considérée est obtenu à partir d'un moment 10 magnétique interpolé à une itération précédente, d'une valeur seuil prédéterminée (4(3o) et d'un vecteur unitaire défini à partir de deux moments magnétiques instantanés successifs de la série temporelle initiale ({m(ti)}N). 7. Procédé selon l'une quelconque des revendications 1 à 6, dans lequel l'étape de positionnement comporte : 15 - une sous-étape d'estimation d'un premier paramètre de dispersion (6(0 m(ti)11 )) représentatif d'une dispersion de l'intensité (II m(ti)11 ) des moments magnétiques instantanés (m(tj)) de l'objet magnétique (2), pour différentes positions dites verticales de l'objet magnétique (2) suivant un axe sensiblement orthogonal au plan de mesure (Pmes) ; 20 - une sous-étape de positionnement de l'objet magnétique (2) à une position verticale pour laquelle la valeur du premier paramètre de dispersion (6(0 m(ti)11 )) est inférieure ou égale à une valeur seuil prédéterminée. 8. Procédé selon l'une quelconque des revendications 1 à 7, comportant préalablement à l'étape (90, 110) de positionnement, une étape (100) de mesure, par lesdits 25 magnétomètres (M,), d'un champ magnétique ambiant (V) en l'absence de l'objet magnétique (2) ; et comportant une sous-étape (121) de soustraction du champ magnétique ambiant (V) aux champs magnétiques (Bei)) préalablement mesurés, de manière à obtenir le champ magnétique généré (13,d(ti)) par l'objet magnétique (2). 3054040 25 9. Procédé selon l'une quelconque des revendications 1 à 8, comportant une étape (150) d'application d'un filtre passe-bas aux valeurs de composantes (m(tj)) des moments magnétiques instantanés (m(tj)) préalablement estimés. 10. Procédé selon la revendication 9, dans lequel le filtre passe-bas est une moyenne 5 glissante sur un nombre (K) donné d'échantillons. 11. Procédé selon la revendication 10, dans lequel le nombre (K) d'échantillons est prédéterminé de sorte qu'un biais (b) défini comme un écart entre une déviation angulaire dite réelle (ar) et la déviation angulaire moyenne estimée () soit minimal, en valeur absolue, pour une valeur dite seuil (ami) de déviation angulaire moyenne. 10 12. Procédé selon l'une quelconque des revendications 1 à 11, comportant une étape de tri de l'objet magnétique (2) selon l'écart de la valeur estimée () de la déviation angulaire moyenne vis-à-vis d'une valeur dite seuil (ami) de déviation angulaire moyenne. 13. Procédé selon l'une quelconque des revendications 1 à 12, comportant une étape d'estimation d'un paramètre dit indicateur de qualité (6R/V1\1), à partir du calcul d'un 15 deuxième paramètre de dispersion (cyR) représentatif d'une dispersion des valeurs d'un paramètre dit rayon instantané (R(tj)) calculé comme étant la norme (II m(tj)-moll ) d'un vecteur instantané (m(tj)-mo) défini entre chaque moment magnétique instantané (m(tj)) et le vecteur invariant (ma) estimés préalablement. 14. Procédé selon la revendication 13, comportant en outre une étape de comparaison 20 du deuxième paramètre de dispersion (cyR) à une valeur seuil (cy(R)) prédéterminée. 15. Procédé selon la revendication 14, dans lequel la prédétermination de la valeur seuil (cy(R)) comporte les sous-étapes suivantes : - obtention, Q fois avec Q un entier naturel non nul, d'une distribution temporelle ({R(tj)}N) dudit rayon instantané (R(tj)), pour un même objet magnétique (2) ; 25 - estimation, pour chaque distribution temporelle ({R(tj)}N), d'une valeur moyenne dite rayon moyen (n) définie comme la moyenne temporelle du rayon instantané ( N), de manière à obtenir une distribution ({R}Q) de Q rayons moyens estimés (R) ; 3054040 26 - estimation, à partir de la distribution ({Rh) des Q rayons moyens estimés (n), d'une valeur moyenne ( Q) et un troisième paramètre de dispersion (cy(R)), dite valeur seuil, représentatif d'une dispersion des valeurs des rayons moyens estimés (R). 16. Système d'estimation (1) d'une déviation angulaire moyenne () entre un axe 5 magnétique et un axe de référence d'un objet magnétique (2), comportant : - un réseau de magnétomètres (M,) apte à mesurer, à chaque instant de mesure d'une durée (T), un champ magnétique perturbé par l'objet magnétique (2) selon au moins un axe de mesure, les magnétomètres (M,) étant répartis spatialement de manière à définir un plan de mesure (Pmes) ; 10 - un organe (7) de maintien et de rotation, apte à maintenir l'objet magnétique dans une position vis-à-vis du plan de mesure, et à mettre en rotation l'objet magnétique (2) suivant son axe de référence (Aret) ; - un processeur (5) configuré pour mettre en oeuvre le procédé selon l'une quelconque des revendications 1 à 15. 15 17. Support d'enregistrement d'informations, comportant des instructions pour la mise en oeuvre du procédé selon l'une quelconque des revendications 1 à 15, ces instructions étant aptes à être exécutées par un processeur.