FR 3079611 A1 20191004 FR 1852637 A 20180327 DD17948 PROCÉDÉ D'INITIALISATION D'UN RÉSEAU DE CAPTEURS [1] L'invention concerne un procédé d'initialisation d'un réseau de capteurs. L'invention concerne également un support d'enregistrement d'informations pour la mise 5 en oeuvre de ce procédé ainsi qu'un appareil implémentant ce procédé. [2] Un réseau de capteurs comporte typiquement au moins quatre capteurs, les capteurs étant reliés mécaniquement ensemble et étant configurés pour mesurer une même grandeur physique selon au moins deux axes non-colinéaires de mesure. Généralement, chacun de ces capteurs est également associé à au moins un coefficient 10 de correction qui permet de corriger la mesure de ce capteur pour compenser des erreurs de mesure. [3] Lors de l'initialisation d'un tel réseau de capteurs, une première mesure de la grandeur physique réalisée par chacun des capteurs est acquise. Cette première mesure précède toutes les autres mesures acquises depuis la mise sous tension du réseau de 15 capteurs. Idéalement, en présence d'une grandeur physique uniformément répartie sur l'ensemble des capteurs, les premières mesures devraient être identiques les unes aux autres. Dès lors, dans ce cas idéal, les écarts de mesures entre les différents capteurs sont essentiellement dus à des erreurs de mesure. Les coefficients de correction peuvent alors être réglés à partir de ces premières mesures pour que, après ce réglage, les écarts 20 de mesure entre les différents capteurs soient presque nuls. [4] Par contre, si lors de l'initialisation du réseau de capteurs, à cause de la présence d'un perturbateur, la grandeur physique n'est pas uniforme, les écarts de mesures entre les différents capteurs ne sont pas nécessairement dus uniquement à des erreurs de mesure. Dans ce cas, les premières mesures ne doivent pas être utilisées pour régler les 25 coefficients de correction. En effet, contrairement aux erreurs de mesure qui sont causées par des caractéristiques internes des capteurs, le perturbateur peut être déplacé ou modifié. Ainsi, suite à un déplacement du perturbateur, comme son retrait, les coefficients de correction réglés ne sont alors plus corrects et restent incorrects tant qu'une nouvelle initialisation du réseau de capteurs n'est pas déclenchée. Dans cet intervalle de temps, les 30 mesures du réseau de capteurs ne sont plus fiables et l'utilisateur de ce réseau n'en est pas informé. Ainsi, il est souhaitable d'augmenter la fiabilité du procédé d'initialisation d'un réseau de capteur notamment vis-à-vis de la présence éventuelle d'un perturbateur. [5] L'invention vise à résoudre ce problème en proposant un procédé d'initialisation plus fiable. Elle a donc pour objet un tel procédé d'initialisation d'un réseau de capteurs 35 comportant au moins quatre capteurs, les capteurs étant reliés mécaniquement ensemble et étant configurés pour mesurer une même grandeur physique selon au moins deux axes 2 DD17948 non-colinéaires de mesure, ce procédé comportant l'acquisition de premières mesures de la grandeur physique réalisées par les capteurs, ces premières mesures précédant des deuxièmes mesures, dans lequel le procédé comporte : - une exécution d'un test statistique sur un ensemble constitué par les premières mesures réalisées par les capteurs, ce test statistique étant apte à détecter au moins une mesure aberrante dans les premières mesures causées par une présence d'un objet, appelé perturbateur, ce perturbateur modifiant de façon non-uniforme la grandeur physique mesurée, - si l'exécution du test statistique détecte au moins une mesure aberrante, alors le procédé comprend le signalement de la présence du perturbateur, et - si l'exécution du test statistique ne détecte pas au moins une mesure aberrante, alors le procédé comprend l'acquisition des deuxièmes mesures de la grandeur physique réalisées par les capteurs. [006] Les modes de réalisation de ce procédé de détection peuvent comporter une ou plusieurs des caractéristiques suivantes : ^ le procédé comprend l'identification et le signalement d'au moins un capteur ayant réalisé une mesure aberrante ; ^ l'exécution du test statistique comporte : - un calcul d'une valeur de référence pour les premières mesures des capteurs, - pour chaque capteur, le calcul d'une distance représentative de l'écart entre la première mesure de ce capteur et la valeur de référence calculée, et - la vérification que les distances calculées satisfont une condition prédéterminée et la détection de la présence d'au moins une mesure aberrante si les distances calculées satisfont cette condition prédéterminée et, en alternance, lorsque les distances calculées ne satisfont pas cette condition prédéterminée, l'absence de détection d'au moins une mesure aberrante ; ^ la condition prédéterminée est la condition suivante : Sdi aiS-rdi, où : - Sdi est la somme des N1 distances calculées les plus grandes, où N1 est un nombre supérieur ou égal à un et inférieur ou égal à Ndmax/2 ou Ndmax/3, où Ndmax est égal au nombre de distances calculées, - STd1 est la somme de toutes les distances calculées, et 3 DD17948 - ai est un coefficient prédéterminé supérieur à zéro et inférieur à un. ^ la distance calculée est la distance de Mahalanobis ; ^ le calcul de la valeur de référence comporte le calcul d'une valeur moyenne ou médiane des premières mesures ; 5 ^ si l'exécution du test statistique ne détecte pas au moins une mesure aberrante, le procédé comporte le réglage de coefficients de correction en fonction des premières mesures pour compenser des erreurs de mesures des capteurs, chaque capteur étant associé à au moins un de ces coefficients de correction, puis la correction des deuxièmes mesures à l'aide de ces coefficients de corrections 10 réglés, puis l'utilisation des deuxièmes mesures corrigées lors de toutes étapes ultérieures exploitant ces deuxièmes mesures. [007] L'invention a également pour objet un support d'enregistrement d'informations, lisible par un microprocesseur, qui comporte des instructions pour l'exécution du procédé objet de l'invention, lorsque ces instructions sont exécutées par le microprocesseur. 15 [008] L'invention a également pour objet une unité de traitement pour la mise en oeuvre du procédé objet de l'invention, cette unité de traitement étant configurée pour acquérir des premières mesures de la grandeur physique réalisées par les capteurs, ces premières mesures précédant des deuxièmes mesures, dans lequel l'unité de traitement est également configurée pour : 20 - exécuter un test statistique sur un ensemble constitué par les premières mesures réalisées par les capteurs, ce test statistique étant apte à détecter au moins une mesure aberrante dans les premières mesures causées par une présence d'un objet, appelé perturbateur, ce perturbateur modifiant de façon non-uniforme la grandeur physique mesurée, 25 - si l'exécution du test statistique détecte au moins une mesure aberrante, alors signaler la présence du perturbateur, et - si l'exécution du test statistique ne détecte pas au moins une mesure aberrante, alors acquérir les deuxièmes mesures de la grandeur physique réalisées par les capteurs. [009] L'invention a également pour objet un système comportant : 30 - un réseau de capteurs comportant au moins quatre capteurs, les capteurs étant reliés mécaniquement ensemble et étant configurés pour mesurer une même grandeur physique selon au moins deux axes non-colinéaires de mesure, et 4 DD17948 - l'unité (30) de traitement objet de l'invention. [0010] Les modes de réalisation de ce système peuvent présenter une ou plusieurs des caractéristiques suivantes : ^ chaque capteur est choisi dans le groupe constitué d'un magnétomètre, d'un 5 capteur de pression, d'un capteur de force, d'un accéléromètre et d'un gyromètre ; ^ chaque capteur comporte au moins un premier et un second axes de mesure, les premiers axes de mesure étant parallèles entre eux et les seconds axes de mesure étant parallèles entre eux ; ^ chaque capteur est un magnétomètre tri-axe ou un capteur de pression biaxes ; 10 ^ le système comporte au moins un aimant permanent apte à être fixé sans aucun degré de liberté sur un crayon, chaque capteur du réseau de capteurs est apte à mesurer l'amplitude du champ magnétique le long d'au moins trois axes non-colinéaires de mesure, et l'unité de traitement est apte à relever la trace du crayon dessinée sur une surface d'écriture, l'unité de traitement étant configurée à cet effet 15 pour localiser l'aimant permanent fixé sur le crayon et pour déduire de la localisation de l'aimant permanent la position de la pointe du crayon. [0011] L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple non limitatif et faite en se référant aux dessins sur lesquels : 20 - la figure 1 est une illustration schématique d'un système de relevé de la trace de la pointe d'un crayon sur un support d'écriture ; - la figure 2 est une illustration schématique d'une partie d'un réseau de magnétomètres utilisé dans le système de la figure 1 ; - la figure 3 est un graphe représentant les champs magnétiques mesurés par des 25 magnétomètres du réseau de la figure 2 ; - la figure 4 est un organigramme d'un procédé de relevé de la trace de la pointe d'un crayon sur un support d'écriture à l'aide du système de la figure 1, et - la figure 5 est un organigramme d'une variante du procédé de la figure 4. [0012] Dans ces figures, les mêmes références sont utilisées pour désigner les mêmes 30 éléments. [0013] Chapitre I : Exemples de modes de réalisation 5 DD17948 [0014] L'invention va maintenant être décrite dans le cas particulier d'un système de relevé de la trace de la pointe d'un crayon. [0015] La figure 1 représente un tel système 2 de relevé de la trace de la pointe 8 d'un crayon 4 sur un support d'écriture 6. L'architecture d'un tel système a déjà été décrite, par 5 exemple, dans le document W02013144338A1. Par conséquent, par la suite, les informations contenues dans ce document ne sont pas reprises ici. Pour plus d'informations sur l'architecture de ce système connu, le lecteur peut consulter ce document W02013144338A1. [0016] Le système 2 comporte un appareil 16 de localisation de la pointe 8. Le crayon 4 10 est équipé d'un aimant permanent 15 fixé sans aucun degré de liberté au crayon 4. [0017] L'appareil 16 permet de localiser l'aimant 15, et donc le crayon 4, dans un repère tridimensionnel XYZ. Par localisation, on entend ici la détermination de la position x, y, z de l'aimant 15 dans le repère XYZ et aussi la détermination de l'orientation du moment magnétique de l'aimant 15 par rapport aux directions X, Y et Z du repère XYZ. 15 [0018] A cet effet, l'appareil 16 comporte un réseau de N capteurs. Selon cet exemple de réalisation, chaque capteur est un magnétomètre tri-axes Mu. Ainsi, dans ce mode de réalisation, la grandeur physique mesurée par chacun des capteurs est le champ magnétique. Typiquement, N est supérieur ou égal à quatre ou cinq et, de préférence, supérieur à seize ou trente-deux. Selon l'exemple de réalisation qui va être décrit, le 20 nombre N de capteur est égal à trente-deux. [0019] Chaque magnétomètre Mu est raccordé par l'intermédiaire d'un bus 28 de transmission d'informations à une unité électronique de traitement 30. [0020] L'unité électronique de traitement 30 est capable de déterminer la position et l'orientation de l'aimant 15 et, à partir de ces informations, d'en déduire la position de la 25 pointe 8 sur le support 6. A cet effet, l'unité électronique de traitement 30 comporte un microprocesseur 32 programmable apte à exécuter des instructions enregistrées sur un support d'enregistrement d'informations. L'unité électronique de traitement 30 comporte donc aussi une mémoire 34 contenant les instructions nécessaires pour l'exécution, par le microprocesseur 32, du procédé de la figure 4 ou 5. 30 [0021] La figure 2 représente une partie des magnétomètres Mu de l'appareil 16. A titre d'illustration, les magnétomètres Mu sont alignés en lignes et en colonnes pour former une matrice. Les indices i et j identifient, respectivement, la ligne et la colonne de cette matrice à l'intersection de laquelle se trouve le magnétomètre Mu. Ici, pour fixer l'orientation des magnétomètres les uns par rapport aux autres, chaque magnétomètre Mu est fixé sans aucun degré de liberté aux autres magnétomètres. 6 DD17948 [0022] Les centres de deux magnétomètres Mu et M,,,+1 immédiatement consécutifs le long d'une ligne i sont séparés par une distance De façon similaire, le centre de deux magnétomètres Mu et M,,1', immédiatement consécutifs le long d'une même colonne j sont séparés par une distance Ici, quelle que soit la ligne i, la distance est la même. 5 Cette distance est donc notée d,. De façon similaire, quelle que soit la colonne j, la distance entre deux magnétomètres est la même. Cette distance est donc notée d,. Dans ce mode de réalisation particulier, les distances d, et d, sont toutes les deux égales à d. [0023] Typiquement, la distance d est comprise entre 1 et 4 cm notamment lorsque : 10 - la puissance de l'aimant permanent 15 est de 0,5 A.m2, - la sensibilité des magnétomètres Mu est de 4*10-7T, et - le nombre de magnétomètres Mu est de trente-deux. [0024] La figure 3 représente dans le repère XYZ par trente-deux flèches les premières mesures du réseau de magnétomètres. Chaque flèche représente la direction du champ 15 magnétique mesuré par un magnétomètre respectif Mu. La longueur de chaque flèche représente l'intensité du champ magnétique mesuré par le magnétomètre. correspondant. [0025] Les premières mesures représentées sur la figure 3 ont été réalisées seulement en présence du champ magnétique terrestre. En particulier, lors de la réalisation de ces premières mesures, aucun objet magnétique, distinct de l'appareil 16, ne se situait à 20 proximité des magnétomètres Mu. Par « à proximité des magnétomètres Mu », on désigne une distance suffisamment courte, entre l'un des magnétomètres Mu et cet objet magnétique, pour modifier de façon non négligeable les mesures d'au moins l'un de ces magnétomètres M. Par exemple, on considère ici que les mesures d'un magnétomètre Mu sont modifiées de façon non négligeable si la présence de l'objet magnétique fait varier 25 l'intensité du champ magnétique mesuré par ce magnétomètres Mu d'au moins 1 % ou d'au moins 5 % par rapport à l'intensité du champ magnétique mesuré dans les mêmes conditions mais en absence de cet objet magnétique. Dans la suite, l'objet magnétique est appelé perturbateur. [0026] Ainsi, on comprend que lors de la réalisation des mesures représentées sur la 30 figure 3, l'aimant 15 n'est pas situé à proximité des magnétomètres M. [0027] Dans ces conditions, étant donné les faibles dimensions du réseau de magnétomètres Mu par rapport aux distances sur lesquelles le champ magnétique terrestre varie, il est admis que le champ magnétique terrestre, au niveau des magnétomètres Mu, est uniforme. Par conséquent, en théorie, les champs magnétiques mesurés par les 35 magnétomètres Mu devraient tous être identiques les uns aux autres. Toutefois, comme le 7 DD17948 montre clairement la figure 3, en pratique et avant qu'une correction appropriée soit appliquée, même lorsque le réseau de magnétomètres Mu est seulement en présence d'un champ magnétique uniforme sur toute sa surface, les champs magnétiques mesurés par les différents magnétomètres Mu ne sont pas identiques entre eux. Cela provient d'erreurs 5 de mesure des magnétomètres M. Une partie de ces erreurs de mesure peut provenir d'un défaut de reproductibilité du procédé de fabrication des magnétomètres. A cause de ce défaut de reproductibilité, les magnétomètres Mu ne sont pas strictement identiques les uns aux autres. Ces erreurs proviennent également du fait que chaque magnétomètre Mu comporte des pièces ferromagnétiques. Au cours de l'utilisation de ces magnétomètres M u, 10 ces pièces ferromagnétiques ne sont pas toujours exposées au même champ magnétique, de sorte qu'il peut exister dans chacun des magnétomètres Mu des champs magnétiques rémanents qui ne sont pas identiques d'un magnétomètre Mu à un autre. Ces champs magnétiques rémanents décalent plus ou moins la mesure de chaque magnétomètre Mu. Ce décalage est connu sous le terme anglais de « offset » et il décale en particulier 15 l'amplitude du champ magnétique mesuré le long d'un axe de mesure. Par exemple, dans le cas des magnétomètres Mu de l'appareil 16, il a été observé que la distribution des offset des magnétomètres Mu est centrée entre 70 pIT et 80 pT. Il a aussi était observé que quelques rares magnétomètres Mu présentent des offset supérieurs à 150 pIT ou 500 pT. [0028] Pour corriger ces erreurs de mesure, chaque magnétomètre Mu est associé à un 20 ou plusieurs coefficients de correction. Ces coefficients de correction sont appliqués aux mesures non corrigées du magnétomètre Mu pour obtenir une mesure corrigée. Typiquement, la mesure non corrigée du magnétomètre Mu le long d'un axe de mesure est combinée, par addition et/ou multiplication, avec l'un de ces coefficients de correction pour obtenir la mesure corrigée le long de cet axe de mesure. Les valeurs de ces coefficients 25 de correction sont réglées, lors d'une phase d'initialisation, de sorte que lorsque les magnétomètres Mu sont placés dans un champ magnétique uniforme, les mesures de tous ces magnétomètres Mu sont identiques à une marge de tolérance prédéterminée près. Par la suite, lorsque l'on parle de « mesure » sans préciser s'il s'agit de mesure corrigée ou non, c'est qu'il s'agit de la mesure corrigée. 30 [0029] Dans le mode de réalisation particulier décrit ici, chaque magnétomètre Mu est associé à trois coefficients de correction distincts, c'est-à-dire à un coefficient de correction par axe de mesure. La valeur de chacun de ces coefficients de correction est déterminée pour compenser l'offset de mesure le long d'un axe de mesure respectif. Ainsi, la mesure corrigée d'un magnétomètre Mu est obtenue en ajoutant à chaque mesure non corrigée sur 35 un de ses axes de mesure, tel que l'un des axes X, Y, Z précédemment définis, le coefficient de correction associé à ce magnétomètre Mu et à cet axe particulier de mesure. Ici, les coefficients de correction associés à chaque magnétomètre Mu sont par exemple enregistrés dans la mémoire 34 et déterminés par l'unité 30 lors de la phase d'initialisation. 8 DD17948 [0030] Pour déterminer les valeurs des coefficients de correction, il faut s'assurer que les premières mesures acquises des magnétomètres lors de cette phase d'initialisation ne sont pas perturbées par un perturbateur. En effet, dans ce dernier cas, il n'est pas possible d'initialiser correctement le réseau de magnétomètres et, en particulier, de déterminer les 5 valeurs des coefficients de corrections à employer pour corriger les deuxièmes mesures suivantes. Or, comme l'illustre la figure 3, les premières mesures acquises lors de la phase d'initialisation présentent une grande dispersion même en absence de perturbateur ce qui rend très difficile la détection d'un perturbateur. Le procédé, mis en oeuvre par le système 2, et qui va maintenant être décrit en référence au procédé de la figure 4, permet de 10 remédier à cette difficulté. [0031] Le procédé débute par une phase 40 d'initialisation du réseau de magnétomètres Mu. Par exemple, la phase 40 est automatiquement déclenchée à chaque mise sous tension de l'appareil 16. Elle peut également être déclenchée à intervalle régulier ou en réponse à l'acquisition par l'unité 30 d'une commande particulière transmise, par exemple, 15 par l'utilisateur. [0032] La phase 40 débute par une étape 42 de détection automatique de la présence à proximité des magnétomètres Mu d'un perturbateur, c'est-à-dire d'un objet magnétique qui modifie de façon non uniforme le champ magnétique mesuré par les magnétomètres M. Cette modification non uniforme du champ magnétique mesuré causée par la présence du 20 perturbateur modifie les mesures réalisées par au moins un magnétomètre le long, généralement, de ses trois axes de mesures. Cela provoque donc en particulier une variation AV1 de l'amplitude mesurée par un premier magnétomètre Mu le long d'un premier axe de mesure. Cette modification provoque également une variation AV2 de l'amplitude mesurée par un second magnétomètre Mu le long de ce même premier axe de mesure. Le 25 second magnétomètre est différent du premier magnétomètre. Chaque variation AV1 et AV2 représente la différence entre les amplitudes mesurées le long du premier axe de mesure en présence et en absence du perturbateur. Puisque le perturbateur modifie de façon non-uniforme le champ magnétique mesuré, il existe dans le réseau de magnétomètres Mu au moins deux magnétomètres pour lesquels l'écart entre les variations 30 AV1 et AV2 est substantielle. Par exemple, l'écart est considéré comme étant substantiel si la condition suivante est vérifiée :IAV1- AV21> p x Max(IVil , 1V21), où : - le symbole « I ...1 » désigne la valeur absolue, - le symbole « Max(..., ...) » désigne la fonction qui calcule le maximum des valeurs entre parenthèses, 35 - le symbole « x » désigne l'opération de multiplication, - V1 et V2 désignent les amplitudes mesurées le long du premier axe de mesure, 9 DD17948 respectivement, par les premier et second magnétomètres, - p est une constante supérieure à 0,01 et, de préférence, supérieure à 0,05 ou 0,1 ou 0,3 ou 0,5. Ici, p est égal à 0,1 par exemple. 5 [0033] Lors de la phase 40 d'initialisation du réseau de magnétomètres, les premières mesures sont les mesures corrigées avec les valeurs des coefficients de correction déterminées lors d'une précédente exécution de la phase 40 ou avec des valeurs par défaut enregistrées dans la mémoire 34, par exemple lors de la fabrication de l'appareil 16. Les valeurs par défaut sont indépendantes des mesures réalisées par les 10 magnétomètres Mu. Par exemple, les valeurs par défaut des coefficients de correction sont utilisés s'il s'agit de la première exécution de la phase 40 ou suite à une commande de restauration des valeurs par défaut des coefficients de correction. [0034] L'étape 42 débute par une opération 44 lors de laquelle les magnétomètres Mu mesurent chacun le champ magnétique actuellement présent. Lors de cette étape, l'unité 15 électronique de traitement 30 acquiert ces premières mesures, puis, par exemple, les corrige en utilisant les précédentes valeurs des coefficients de correction enregistrées dans la mémoire 34. Ces premières mesures corrigées sont alors enregistrées dans une matrice A de dimension (n, p), où : - n est le nombre de lignes et il est égal ici au nombre N de magnétomètres Mu du réseau, 20 et - p est le nombre d'axes de mesure de chaque magnétomètre Mu c'est-à-dire ici trois. [0035] Chaque ligne k de la matrice A comporte les trois premières mesures acquises selon les trois axes de mesure d'un magnétomètre Mu particulier désigné, par la suite, par la référence Mk. Ainsi, dans la suite de cette description, un magnétomètre particulier du 25 réseau de l'appareil 16 peut aussi bien être désigné par la référence M;; que par la référence Mk. Ici, l'indice k est le numéro d'ordre de la ligne de la matrice A et il est donc compris entre 1 et N. Plus précisément, chaque ligne k comporte dans l'ordre : - la mesure MMk,x de la projection orthogonale du champ magnétique mesuré par le magnétomètre Mk sur l'axe de mesure parallèle à la direction X, 30 - la mesure MMk,y de la projection orthogonale du champ magnétique mesuré par le magnétomètre Mk sur l'axe de mesure parallèle à la direction Y, et - la mesure MMk,Z de la projection orthogonale du champ magnétique mesuré par le magnétomètre Mk sur l'axe de mesure parallèle à la direction Z. 10 DD17948 [0036] Ici, les mesures acquises lors de l'opération 44 sont appelées « premières mesures » car elles correspondent aux premières mesures acquises après la mise sous tension de l'appareil 16 ou suite à une commande particulière, telle qu'une commande de réinitialisation des valeurs des coefficients de corrections. Ces premières mesures sont 5 destinées à être utilisées ensuite pour régler les nouvelles valeurs des coefficients de corrections. Ainsi, ces premières mesures peuvent présenter une dispersion importante comme représenté sur la figure 3. [0037] Les opérations suivantes de l'étape 42 sont exécutées par l'unité 30 et correspondent à l'exécution d'un test statistique sur les premières mesures. 10 [0038] Dans ce mode de réalisation, lors d'une opération 46, pour chaque axe de mesure, une valeur de référence pour l'amplitude du champ magnétique mesuré le long de cet axe est calculée. Plus précisément, lors de l'opération 46, ici, un champ magnétique moyen, noté Xm par la suite, est calculé. Les composantes du champ magnétique moyen Xm le long des directions X, Y et Z sont notées, respectivement, Xm,x ; 15 Xm,y et Xm,z. Chaque composante Xm,x, Xm,y et Xm,z est une valeur moyenne des mesures, respectivement, MMk,x, MMk,yet MMk,z pour k variant de 1 à N. Ainsi, la composante Xm,x est calculée à l'aide de la relation suivante : [0039] Les composantes Xm,y et Xm,z sont calculées de façon similaire. Les composantes Xm,x, Xm,y et Xm,z sont les valeurs de références calculées lors de l'opération 46 pour les 20 axes de mesures parallèles aux directions, respectivement, X, Y et Z. [0040] Ensuite, lors d'une opération 48, pour chaque magnétomètre Mk, une distance dk représentative de l'écart entre les mesures MMk,x, MMk,y et MMk,z de ce magnétomètre Mk et les valeurs de référence, respectivement Xm,x, Xm,y et Xm,z, est calculée. Ici, la distance utilisée est la distance de Mahalanobis. 25 [0041] Pour cela, une matrice centrée X, de dimension (n, p) est déterminée en modifiant chaque composante de la matrice A à l'aide du champ magnétique moyen Xm. Plus précisément, pour construire la première colonne de la matrice X,, la valeur de référence Xm,x est soustraite à chacune des composantes MMk,x contenues dans la première colonne de la matrice A. Il est procédé de façon similaire pour construire les deuxième et troisième 30 colonnes de la matrice X, mais en utilisant, respectivement, les valeurs de référence Xm,y et Xm,z. 11 DD17948 [0042] Ensuite, une matrice S de covariance est construite à l'aide de la relation suivante : S = (XcT Xc) / (N-1). Dans ce texte, l'exposant « T » désigne l'opération matricielle de transposition. La matrice S ainsi obtenue est une matrice de dimension (p, p) qui est inversible. La distance dk de Mahalanobis de chaque magnétomètre Mk est alors calculée, 5 par exemple, à l'aide de la relation suivante : dk2 = XC,k S1 X,,kT, où X,,k désigne la k-ième ligne de la matrice Xc et l'exposant « -1» indique qu'il s'agit de l'inverse de la matrice S. [0043] Lors d'une opération 50, il est vérifié si les distances dk calculées satisfont une condition prédéterminée. Si les distances dk calculées satisfont cette condition prédéterminée, alors la présence de mesures aberrantes est détectée et donc la présence 10 d'un perturbateur est détectée. Dans ce cas, la présence d'un perturbateur à proximité du réseau de magnétomètres est signalée. Selon une réalisation particulière, il s'en suit une étape, dans laquelle le ou les capteurs dont les mesures sont perturbées sont identifiés et signalés. [0044] À l'inverse, si les distances dk calculées ne satisfont pas cette condition 15 prédéterminée, alors la présence de mesures aberrantes et donc d'un perturbateur n'est pas détectée. Autrement dit, si les distances dk ne satisfont pas cette condition prédéterminée, l'absence de perturbateur est détectée. Dans ce dernier cas, le signalement de la présence d'un perturbateur est inhibé. [0045] La condition prédéterminée est par exemple enregistrée dans la mémoire 34 lors 20 de la fabrication de l'appareil 16. Cette condition prédéterminée est conçue pour discriminer efficacement une non-uniformité intrinsèque d'une non-uniformité extrinsèque. Par « non-uniformité intrinsèque », on désigne ici une non-uniformité des premières mesures des magnétomètres Mg qui est causée par des éléments intrinsèques à l'appareil 16, c'est-à-dire des éléments qui sont fixés sans aucun degré de liberté au réseau de 25 magnétomètres M. Par exemple, typiquement, il s'agit de pièces en matériau ferromagnétique comme les pièces ferromagnétiques qui sont situées dans chacun des magnétomètres Mij. [0046] Par « non-uniformité extrinsèque », on désigne ici une non-uniformité des premières mesures des magnétomètres Mg qui est causée par le perturbateur, c'est-à-dire 30 par un élément extérieur au réseau de magnétomètres Mg et généralement extérieur à l'appareil 16. Le perturbateur est ici généralement un objet magnétique mobile qui peut librement être déplacé par rapport à l'appareil 16. Par exemple, le perturbateur peut être l'aimant 15 du crayon 4 si celui-ci est à proximité du réseau de magnétomètres Mq lors de l'exécution de la phase 40 d'initialisation du réseau de magnétomètres. 35 [0047] Dans ce mode de réalisation, la condition prédéterminée utilisée est la suivante : Sdl al STd1, où 12 DD17948 - Sd1 est la somme des N1 distances dk calculées les plus grandes, - STd1 est la somme de toutes les distances dk calculées, et - ai est un coefficient prédéterminé supérieur à 0 et inférieur à 1. [0048] Par la suite, cette condition est appelée « condition 1) ». 5 [0049] N1 est un nombre prédéterminé enregistré dans la mémoire 34. Le nombre N1 est supérieur ou égal à un et inférieur à Ndmax, OÙ Ndmax est égal au nombre de distances dk calculées c'est-à-dire ici à N. De préférence, N1 est inférieur ou égal à Ndmax/2 ou Ndmax/3. De préférence, Ni est également supérieur ou égal à deux ou trois. Ici, N1 est choisi égal à dix ou onze. 10 [0050] Le coefficient ai permet de régler la sensibilité de la détection de la présence d'un perturbateur. En pratique, lorsque N1 est choisi inférieur ou égal à Ndmax/2 ou Ndmax/3, ai est choisi supérieur ou égal à 0,5 ou à 0,7. Ici, ai est égal à 0,7. [0051] Il a été découvert que la condition 1) est, parmi l'ensemble des conditions possibles présentées à la fin de cette description, celle qui permet de détecter, de la façon 15 la plus fiable, la présence ou l'absence d'un perturbateur. [0052] À ce sujet, les essais suivants ont été réalisés en utilisant le crayon 4 en tant que perturbateur : - Essai 1) : le crayon 4 est posé horizontalement à un endroit quelconque sur la face du réseau de magnétomètres puis l'étape 42 de détection de la présence d'un perturbateur 20 est exécutée. - Essai 2) : le crayon 4 est disposé le long d'un bord du réseau de magnétomètres parallèle à la direction X ou Y puis l'étape 42 de détection de la présence d'un perturbateur est exécutée. - Essai 3) : le crayon 4 est maintenu vertical, la pointe 8 en appui à un endroit quelconque 25 sur le support 6, puis l'étape 42 de détection de la présence d'un perturbateur est exécutée. [0053] Chacun des essais 1) à 3) a été réitéré plusieurs fois en modifiant à chaque fois la position du crayon 4 tout en respectant les conditions de l'essai. Dans le cas de l'essai 1, le taux de détection correcte, également connu sous l'acronyme TPR (« True Positive 30 Rate »), et le taux de fausse détection, connu sous l'acronyme FPR (« False Positive Rate ») sont tous les deux égaux à 100 % . Dans le cas de l'essai 2, le taux TPR est supérieur à 99,6 % et le taux FPR est inférieur à 0,3 %. Dans le cas de l'essai 3, le taux TPR est supérieur à 50 % ou 70 % selon l'appareil 16 utilisé et le taux FPR est inférieur à 13 DD17948 0,2 %. [0054] Si la condition 1) est satisfaite, les magnétomètres Mk identifiés comme étant perturbés par la présence du perturbateur sont, par exemple, les N1 magnétomètres dont les distances dk sont les plus grandes. 5 [0055] Si à l'issue de l'étape 42 de détection de la présence d'un perturbateur, la présence d'aucun perturbateur n'a été signalée, le réglage des coefficients de correction à partir des premières mesures acquises lors de l'opération 44 de mesure est autorisé. Ainsi, le procédé se poursuit par une étape 52 de réglage automatique des coefficients de correction des magnétomètres M. 10 [0056] Par exemple, lors de l'étape 52, pour chaque magnétomètre Mk, les nouvelles valeurs des trois coefficients de correction associés respectivement aux directions X, Y et Z sont corrigés en y retranchant les mesures, respectivement, MM k,x, MMk,y et MMk,Z. Ainsi, par la suite, chaque nouvelle mesure corrigée est obtenue en utilisant ces nouvelles valeurs des coefficients de correction. Ces nouvelles mesures sont appelées par la suite 15 « deuxièmes mesures ». Cela a pour conséquence de compenser non seulement les erreurs de mesure mais aussi le champ magnétique terrestre mesuré lors de l'opération 44 de mesure. Ainsi, après l'étape 52 de réglage des coefficients de correction, en absence de perturbateur, les deuxièmes mesures corrigées des magnétomètres Mu sont toutes nulles à une marge de tolérance prédéterminée près. 20 [0057] À l'issue de l'étape 52 de réglage des coefficients de correction, la phase 40 d'initialisation est achevée et le procédé se poursuit par une phase 58 d'utilisation de deuxièmes mesures, ici, pour relever la trace de la pointe 8 du crayon 4. [0058] Si à l'issue de l'étape 42 de détection de la présence d'un perturbateur, la présence d'un perturbateur a été signalée, en réponse à ce signalement, le réglage des 25 coefficients de correction à partir des premières mesures acquises lors de l'opération 44 de mesure est interdit. Pour cela, l'exécution de l'étape 52 de réglage des coefficients de correction est inhibée. Par exemple, dans un premier mode de réalisation, l'utilisation du réseau de magnétomètres est bloquée et l'exécution de la phase 58 d'utilisation des deuxièmes mesures est interdite. Dans un deuxième mode de réalisation, en réponse au 30 signalement de la présence d'un perturbateur : - le procédé retourne à l'étape 42 de détection de la présence d'un perturbateur, ou - le procédé se poursuit par la phase 58 d'utilisation des deuxièmes mesures mais en utilisant les coefficients de correction réglés lors d'une précédente exécution de l'étape 52 de réglage des coefficients de correction ou en utilisant les valeurs par défaut de ces 35 coefficients de correction. 14 DD17948 [0059] La phase 58 d'utilisation des deuxièmes mesures est par exemple réalisée comme décrit dans la demande W02013144338. Par conséquent, cette phase 58 n'est pas décrite ici plus en détail et le lecteur est renvoyé au texte de cette demande pour plus d'informations. Lors de cette phase 58, des deuxièmes mesures des magnétomètres Mu 5 sont acquises à intervalle régulier par l'unité 30 puis corrigées à l'aide des coefficients de correction. Ensuite, seules les deuxièmes mesures ainsi corrigées sont utilisées pour déterminer la position et l'orientation de l'aimant 15 et donc la position de la pointe 8. [0060] La figure 5 représente un autre procédé de relevé de la trace du crayon 4 susceptible d'être implémenté dans l'appareil 16 à la place du procédé de la figure 4. Le 10 procédé de la figure 5 est identique au procédé de la figure 4 sauf que la phase 40 d'initialisation est remplacée par une phase 60. La phase 60 est identique à la phase 40 d'initialisation du réseau de magnétomètres sauf que l'étape 42 de détection de la présence d'un perturbateur est remplacée par une étape 62 de détection d'un perturbateur. L'étape 62 est identique à l'étape 42 de détection de la présence d'un 15 perturbateur sauf que : - les opérations 46 de calcul d'une valeur de référence et 48 de calcul de la distance dk sont remplacées par une opération 66, et - l'opération 50 de vérification de la condition prédéterminée est remplacée par une opération 70. 20 [0061] Lors de l'opération 66, pour chaque magnétomètre Mk, l'unité 30 calcule l'intensité Bk du champ magnétique mesuré par ce magnétomètre à partir des mesures mm MM ----k,y5 Mmk,Z acquises lors de l'opération 44 de mesure. Classiquement, l'intensité Bk est calculée à l'aide de la relation suivante : Bk = sqrt(MMk,x2 mmk,y2 + Mmk,Z2), où sqrt(..) est la fonction qui retourne la racine carrée du nombre entre parenthèses. 25 [0062] L'opération 70 est identique à l'opération 50 de vérification de la condition prédéterminée sauf que la condition 1) s'exprime alors, par exemple, sous la forme suivante : aiSTBi, où : - Sin est la somme des N1 intensités Bk les plus grandes parmi l'ensemble des intensités Bk calculées, 30 - STB1 est la somme de toutes les intensités Bk calculées. [0063] Chapitre II : Variantes : [0064] Chapitre 11.1. Variantes de la distance dl : [0065] D'autres méthodes sont possibles pour calculer une distance représentative de l'écart entre les mesures d'un capteur et les mesures des autres capteurs. Par exemple, 15 DD17948 lors de l'opération 48 de calcul de la distance dk, le champ magnétique moyen Xm peut être remplacé par un champ magnétique médian. Dans ce cas, les valeurs de référence Xm,x, Xm,y et Xm,z sont égales à la valeur médiane, au lieu de la valeur moyenne, des mesures réalisées selon les axes de mesure X, Y et Z. 5 [0066] La distance peut également être calculée en utilisant d'autres valeurs de référence qu'une valeur moyenne ou médiane. Par exemple, une valeur de référence peut être prise égale au maximum ou au minimum des mesures des N capteurs selon un axe de mesure. [0067] Il existe aussi d'autres méthodes de calcul de la distance dk de Mahalanobis. Par exemple, l'unité 30 réalise une analyse en composantes principales connue sous 10 l'acronyme PCA (« Principal Component Analysis ») pour déterminer les trois directions u, et w de variance maximale qui forme une base orthonormée. Ensuite, la distance dk est calculée à l'aide de la relation suivante dk2-=ky 2 L-1, où - yk est la k-ième ligne d'une matrice Y de dimension (n, p) qui contient chacune des projections orthogonales du vecteur Xc,ksur, respectivement, les axes u, v et w, et 15 - L est le vecteur dont les composantes sont formées par les valeurs propres de la matrice S de covariance. [0068] En variante, les premières mesures utilisées sont directement les mesures brutes des magnétomètres sans leur appliquer une quelconque correction. [0069] Chapitre 11.2. Variantes de la condition prédéterminée : 20 [0070] D'autres conditions prédéterminées que la condition 1) précédemment décrite peuvent être utilisées. Par exemple, la condition 1) peut être remplacée par l'une des conditions suivantes : - Condition 2) : si N2 distances dk sont supérieures à un seuil S2 prédéterminé, alors la présence d'un perturbateur est détectée, où N2 est supérieur ou égal à un et inférieur ou 25 égal à Ndmax. De préférence, le seuil S2 est calculé en fonction de 'id et ed, où : - µd est la valeur moyenne des distances dk calculées, et - ed est l'écart-type de la distribution des distances dk calculées autour de la valeur Par exemple, le seuil S2 est calculé à l'aide de la relation suivante : S2 = µd + [3d Ed, où [3d 30 est une constante prédéterminée. - Condition 3) : si N3 distances dk sont supérieures à un seuil S3 prédéterminé, alors la 16 DD17948 présence d'un perturbateur est détectée, où N3 est supérieur ou égal à un et inférieur ou égal à Ndmax. De préférence, le seuil S3 est calculé en fonction de µm et Em, où : - pm est la valeur médiane des distances dk calculées, et - Em est la déviation maximale en valeur absolue également connus sous l'acronyme 5 MAD (Maximum Absolute Deviation). Par exemple, le seuil S3 est calculé à l'aide de la relation suivante : S3 = Prn + [3m Em, où pm est une constante prédéterminée. [0071] Les différentes conditions décrites ici peuvent être utilisées quel que soit le mode de calcul de la distance. Elles peuvent aussi facilement être adaptées au cas où la 10 distance dk est remplacée par une autre grandeur comme l'amplitude Bk. [0072] D'autres conditions prédéterminées peuvent être générées par exemple en combinant à l'aide d'opérateurs logiques plusieurs des conditions 1) à 3) précédemment décrites. Par exemple, en variante, la condition prédéterminée est : « condition 1) » OU « condition 2) » où « condition 1) » et « condition 2) » sont, respectivement, les conditions 15 1) et 2)-précédemment décrites. [0073] Chapitre 11.3. Autres variantes : [0074] Dans une version simplifiée, lors de l'opération 44, les premières mesures ne subissent pas une première correction en utilisant les précédentes valeurs des coefficients de correction ou les valeurs par défaut de ces coefficients de correction. 20 [0075] Dans un autre mode de réalisation, lorsqu'un perturbateur est détecté à l'issue de l'étape 42 de détection de la présence d'un perturbateur, le procédé se poursuit par une étape automatique de signalement de la présence de ce perturbateur. Par exemple, lors de cette étape de signalement, l'unité 30 commande une interface homme-machine pour signaler à un utilisateur de l'appareil 16 la présence d'un perturbateur à proximité du 25 réseau de capteurs. Par exemple, l'unité 30 commande l'affichage d'un message écrit ou le jeu d'un message sonore pour signaler la présence d'un perturbateur. Ce message peut également demander à l'utilisateur d'éloigner le perturbateur de l'appareil 16 pendant la phase d'initialisation. [0076] Pour la mise en oeuvre de la phase d'initialisation, il n'est pas nécessaire que les 30 distances entre les magnétomètres soient connues. Il n'est pas non plus nécessaire que les magnétomètres soient fixés les uns aux autres sans aucun degré de liberté. Par exemple, en variante, les liaisons mécaniques entre les magnétomètres autorisent un déplacement en translation des magnétomètres les uns par rapport aux autres tout en interdisant tout déplacement en rotation des magnétomètres les uns par rapport aux 17 DD17948 autres. [0077] En variante, si les orientations les unes par rapport aux autres des axes de mesure des magnétomètres Mu sont connues, alors il n'est pas nécessaire que les axes de mesure de ces différents magnétomètres Mu soient parallèles entre eux. En effet, par un 5 simple changement de repère, l'unité 30 est capable de convertir les mesures le long des axes de mesure d'un magnétomètre donné en mesures réalisées, par exemple, le long des directions X, Y et Z du repère XYZ. [0078] Les différentes variantes de la structure du système décrites dans le document W02013144338 sont applicables à la structure du système 2 décrit ici. De même les 10 différentes variantes du procédé de relevé de la trace d'un crayon décrites le document W02013144338 sont applicables ici. [0079] Le réseau de magnétomètres Mu, et donc la phase 40 ou 60 d'initialisation, peut être utilisé dans d'autres appareils qu'un appareil de relevé de la trace d'un crayon. Par exemple, le procédé d'initialisation décrit ici peut aussi être implémenté dans les appareils 15 décrits dans les documents suivants : W02014135421, W02014079740 et W02014180796. [0080] De nombreux autres modes de réalisation de l'étape 52 de réglage automatique des coefficients de correction sont possibles. Par exemple, dans un autre mode de réalisation, les coefficients de correction des magnétomètres Mu sont réglés de manière à 20 compenser seulement les erreurs de mesure sans compenser le champ magnétique terrestre. Pour cela, à titre d'illustration, pour chaque magnétomètre Mu, les trois coefficients de correction associés, respectivement aux directions X, Y et Z, sont pris égaux à, respectivement, MMk,,- Xm,x, MMk,y - Xm,y et MMk,z - Xm,z. [0081] Ce qui a été décrit ici dans le cas particulier où les capteurs sont des 25 magnétomètres s'applique à tout réseau de capteurs à au moins deux axes de mesure. En particulier, cela s'applique à des réseaux de capteurs mesurant une autre grandeur physique que le champ magnétique. Par exemple, la grandeur physique mesurée par les capteurs de cet autre réseau peut être une force, une pression, une vitesse ou une accélération. Ainsi, ce qui a été décrit s'applique à des réseaux de capteurs de force ou de 30 pression, à des réseaux d'accéléromètres ou à des réseaux de gyromètres ou tout autre réseau de capteurs. [0082] Ce qui a été décrit ici s'applique aussi à tout réseau de capteurs dans lequel chaque capteur comporte plus de trois axes de mesure non colinéaires ou, au contraire, seulement deux axes de mesure non colinéaires. En particulier, ce qui a été décrit ici 35 s'applique à des réseaux de capteurs de pression dans lequel chaque capteur est apte à mesurer la pression selon une direction normale et une force de cisaillement selon une 18 DD17948 direction orthogonale à cette direction normale. Par exemple, un tel capteur est décrit dans le document US20100162832A1. [00831 Chapitre III : Avantages des modes de réalisation décrits: [0084] Le fait de détecter des mesures aberrantes dans les premières mesures permet 5 de signaler la présence d'un perturbateur. Ainsi, l'absence de fiabilité des premières mesures est signalée de sorte que des mesures correctives ou préventives peuvent être ensuite mises en oeuvre pour ne pas exploiter des premières mesures non-fiables. Par exemple, comme décrit ci-dessus, lorsque la présence d'un perturbateur est signalée, le réglage des coefficients de corrections à partir des premières mesures est inhibé. Le 10 fonctionnement du réseau de capteur peut aussi être interdit. Ainsi, le procédé d'initialisation décrit ici est plus fiable. [0085] Le fait d'utiliser une distance entre ce que mesure un capteur et une valeur de référence pour cette mesure permet de limiter le taux de fausse détection et le taux de fausse absence de détection par rapport à un procédé identique mais utilisant directement 15 l'amplitude de la grandeur physique mesurée et non pas une distance par rapport à une valeur de référence. [0086] Le fait d'utiliser la condition 1) permet de limiter le taux de fausse détection et le taux de fausse absence de détection par rapport à d'autres conditions possibles, comme par exemple les conditions 2) et 3) précédemment décrites. 20 [0087] Le fait d'utiliser la distance de Mahalanobis permet aussi d'augmenter la fiabilité de la détection de la présence d'un perturbateur par rapport à d'autres procédés identiques mais dans lesquels une autre distance, comme par exemple la distance euclidienne, est utilisée. [0088] Le fait d'utiliser des valeurs préenregistrées pour régler les coefficients de 25 correction lorsqu'un perturbateur est détecté permet d'obtenir un réglage des capteurs qui est meilleur que celui qui serait obtenu si l'étape REVENDICATIONS1. Procédé d'initialisation d'un réseau de capteurs comportant au moins quatre capteurs, les capteurs étant reliés mécaniquement ensemble et étant configurés pour mesurer une même grandeur physique selon au moins deux axes non-colinéaires de mesure, ce procédé comportant l'acquisition (44) de premières mesures de la grandeur physique réalisées par les capteurs, ces premières mesures précédant des deuxièmes mesures, caractérisé en ce que le procédé comporte : - une exécution (46, 48, 50) d'un test statistique sur un ensemble constitué par les premières mesures réalisées par les capteurs, ce test statistique étant apte à détecter au moins une mesure aberrante dans les premières mesures causées par une présence d'un objet, appelé perturbateur, ce perturbateur modifiant de façon non-uniforme la grandeur physique mesurée, - si l'exécution du test statistique détecte au moins une mesure aberrante, alors le procédé comprend le signalement (50) de la présence du perturbateur, et - si l'exécution du test statistique ne détecte pas au moins une mesure aberrante, alors le procédé comprend l'acquisition (58) des deuxièmes mesures de la grandeur physique réalisées par les capteurs. 2. Procédé selon la revendication 1, dans lequel le procédé comprend l'identification et le signalement (50) d'au moins un capteur ayant réalisé une mesure aberrante. 3. Procédé selon la revendication 1 ou 2, dans lequel l'exécution du test statistique comporte : - un calcul (46) d'une valeur de référence pour les premières mesures des capteurs, - pour chaque capteur, le calcul (48) d'une distance représentative de l'écart entre la première mesure de ce capteur et la valeur de référence calculée, et - la vérification (50) que les distances calculées satisfont une condition prédéterminée et 30 la détection de la présence d'au moins une mesure aberrante si les distances calculées satisfont cette condition prédéterminée et, en alternance, lorsque les distances calculées 20 DD17948 ne satisfont pas cette condition prédéterminée, l'absence de détection d'au moins une mesure aberrante. 4. Procédé selon la revendication 3, dans lequel la condition prédéterminée est la 5 condition suivante : Sdi aiS-rdi, où : - Sdi est la somme des N1 distances calculées les plus grandes, où N1 est un nombre supérieur ou égal à un et inférieur ou égal à Ndmax/2 ou Ndmax/3, où Ndmax est égal au nombre de distances calculées, - STd1 est la somme de toutes les distances calculées, et 10 - ai est un coefficient prédéterminé supérieur à zéro et inférieur à un. 5. Procédé selon l'une quelconque des revendications 3 à 4, dans lequel la distance calculée est la distance de Mahalanobis. 15 6. Procédé selon l'une quelconque des revendications 3 à 5, dans lequel le calcul de la valeur de référence comporte le calcul d'une valeur moyenne ou médiane des premières mesures. 7. Procédé selon l'une quelconque des revendications précédentes, dans lequel : 20 - si l'exécution du test statistique ne détecte pas au moins une mesure aberrante, le procédé comporte le réglage (52) de coefficients de correction en fonction des premières mesures pour compenser des erreurs de mesures des capteurs, chaque capteur étant associé à au moins un de ces coefficients de correction, puis - la correction (58) des deuxièmes mesures à l'aide de ces coefficients de corrections 25 réglés, puis - l'utilisation (58) des deuxièmes mesures corrigées lors de toutes étapes ultérieures exploitant ces deuxièmes mesures. 8. Support (34) d'enregistrement d'informations, lisible par un microprocesseur (32), 21 DD17948 caractérisé en ce qu'il comporte des instructions pour l'exécution d'un procédé conforme à l'une quelconque des revendications précédentes, lorsque ces instructions sont exécutées par le microprocesseur (32). 5 9. Unité (30) de traitement pour la mise en oeuvre d'un procédé d'initialisation conforme à l'une quelconque des revendications 1 à 7, cette unité de traitement étant configurée pour acquérir des premières mesures de la grandeur physique réalisées par les capteurs, ces premières mesures précédant des deuxièmes mesures, caractérisée en ce que l'unité de traitement est également configurée pour : 10 - exécuter un test statistique sur un ensemble constitué par les premières mesures réalisées par les capteurs, ce test statistique étant apte à détecter au moins une mesure aberrante dans les premières mesures causées par une présence d'un objet, appelé perturbateur, ce perturbateur modifiant de façon non-uniforme la grandeur physique mesurée, 15 - si l'exécution du test statistique détecte au moins une mesure aberrante, alors signaler la présence du perturbateur, et - si l'exécution du test statistique ne détecte pas au moins une mesure aberrante, alors acquérir les deuxièmes mesures de la grandeur physique réalisées par les capteurs. 20 10. Système comportant : - un réseau de capteurs comportant au moins quatre capteurs, les capteurs étant reliés mécaniquement ensemble et étant configurés pour mesurer une même grandeur physique selon au moins deux axes non-colinéaires de mesure, et - une unité (30) de traitement, 25 caractérisé en ce que l'unité (30) de traitement est conforme à la revendication 9. 11. Système selon la revendication 10, dans lequel chaque capteur est choisi dans le groupe constitué d'un magnétomètre, d'un capteur de pression, d'un capteur de force, d'un accéléromètre et d'un gyromètre. 22 DD17948 12. Système selon la revendication 10, dans lequel chaque capteur comporte au moins un premier et un second axes de mesure, les premiers axes de mesure étant parallèles entre eux et les seconds axes de mesure étant parallèles entre eux. 5 13. Système selon la revendication 10, dans lequel chaque capteur est un magnétomètre tri-axe ou un capteur de pression biaxes. 14. Système selon la revendication 13, dans lequel : - le système comporte au moins un aimant permanent (15) apte à être fixé sans aucun 10 degré de liberté sur un crayon, - chaque capteur du réseau de capteurs est apte à mesurer l'amplitude du champ magnétique le long d'au moins trois axes non-colinéaires de mesure, - l'unité (30) de traitement est apte à relever la trace du crayon dessinée sur une surface d'écriture, l'unité (30) de traitement étant configurée à cet effet pour localiser l'aimant 15 permanent fixé sur le crayon et pour déduire de la localisation de l'aimant permanent la position de la pointe (8) du crayon.