Brevet FR3015053A1 — publié le 20150619

PROCEDES DE LOCALISATION D'USTENSILES MOBILES PRESENTES DEVANT UN RÉSEAU DE MAGNÉTOMÈTRES [1] L'invention concerne un procédé et un dispositif de localisation d'ustensiles 5 mobiles présentés devant un réseau de magnétomètres. L'invention concerne également un support d'enregistrement d'informations pour la mise en oeuvre de ce procédé. [2] Les demandes de brevet FR2 988 862 et FR2 988 874 divulguent des procédés de localisation d'ustensiles mobiles présentés devant un réseau de 10 magnétomètres comportant N magnétomètres tri-axes reliés mécaniquement les uns aux autres sans aucun degré de liberté pour conserver une distance connue entre chacun des magnétomètres, chaque ustensile comportant au moins un objet magnétique, où N est un nombre entier supérieur ou égal à cinq. Ces procédés connus comportent : 15 a) la mesure par chacun des magnétomètres de l'amplitude du champ magnétique le long de chacun de ses axes de mesure, b) l'estimation à partir des mesures des magnétomètres, des valeurs de plusieurs variables en résolvant un système actuel d'équations reliant ces variables à chaque mesure d'un magnétomètre tri-axes du réseau, ces variables correspondant aux 20 positions, orientations et amplitudes des moments magnétiques des objets magnétiques des ustensiles à localiser, - l'itération des étapes a) et b) à des instants successifs dans le temps pour obtenir les valeurs estimées de chaque variable à ces différents instants successifs. [3] Ces procédés fonctionnent particulièrement bien. En particulier, ces procédés 25 sont capables de localiser de nombreux objets magnétiques mobiles différents et donc de nombreux ustensiles. De plus, ils sont aussi capables de localiser simultanément plusieurs ustensiles mobiles mécaniquement indépendants les uns des autres. A cause de ces capacités, on dit que ces procédés sont « flexibles ». [4] Pour être capable de localiser simultanément de nombreux ustensiles 30 mobiles, le système d'équations comporte de nombreuses variables à estimer. Typiquement, pour chaque objet magnétique de l'ustensile à localiser, le système d'équations comporte six variables à estimer. [5] L'estimation des valeurs des nombreuses variables de ce système nécessite l'exécution d'un grand nombre d'opérations par le calculateur car ce système 35 d'équations est non linéaire. Or, plus le nombre d'opérations à exécuter par le calculateur est important, plus la consommation électrique de ce calculateur augmente. [6] L'invention vise à réduire la consommation électrique du calculateur qui exécute un tel procédé de localisation sans réduire la flexibilité du procédé. [007] Elle a donc pour objet un tel procédé dans lequel le procédé comporte : c) l'établissement des mesures des magnétomètres lorsqu'ils sont uniquement en présence d'un champ magnétique de référence généré seulement par les objets magnétiques de tous les ustensiles immobiles présents devant le réseau de 5 magnétomètres, et d) le calcul de l'amplitude d'une perturbation à partir des écarts entre les mesures réalisées lors de l'étape a) et les mesures établies lors de l'étape c), et - si l'amplitude de la perturbation franchit un seuil prédéterminé, l'étape b) est exécuté, sinon l'exécution de l'étape b) est omise et le procédé retourne à l'étape a) 10 sans avoir exécuté l'étape b). [008] Le procédé ci-dessus permet d'éviter de déclencher des exécutions inutiles de l'étape b), c'est-à-dire des exécutions de l'étape b) alors qu'aucun ustensile n'est déplacé. Il permet donc de réduire la consommation électrique du calculateur exécutant ce procédé. Ce procédé est aussi plus rapide à exécuter. 15 [009] De plus, il permet de détecter qu'une exécution de l'étape b) est nécessaire sans avoir à résoudre après chaque itération de l'étape a) le système actuel d'équations, par exemple, pour déterminer une augmentation de l'erreur d'estimation. [cm cl En effet, l'amplitude de la perturbation se calcule directement à partir des valeurs établies pour les mesures des magnétomètres lorsque les ustensiles sont 20 immobiles et des mesures réalisées lors de la dernière l'exécution de l'étape a). Il n'est donc pas nécessaire d'exécuter l'étape b) pour s'apercevoir que son exécution était inutile. Cela réduit également le nombre d'opérations réalisées par le calculateur et donc sa consommation électrique lorsqu'il exécute ce procédé. [0011] Les modes de réalisation de ce procédé peuvent comporter une ou plusieurs 25 des caractéristiques suivantes : ^ si l'amplitude de la perturbation franchit un seuil prédéterminé, le procédé comporte l'identification de l'ustensile qui a été déplacé parmi les ustensiles contenus dans une liste d'ustensiles immobiles, cette liste d'ustensiles immobiles contenant pour chaque objet magnétique de chaque ustensile 30 immobile : sa position, son orientation et l'amplitude de son moment magnétique, cette identification de l'ustensile comportant : - le calcul du barycentre de la perturbation, ce barycentre étant le barycentre des positions des différents magnétomètres du réseau, dans lequel la position de chaque magnétomètre est pondérée par l'écart entre les mesures réalisées 35 par ce magnétomètre lors de l'étape a) et les mesures établies pour ce magnétomètre lors de l'étape c), et - la sélection, dans la liste d'ustensiles immobiles, de l'ustensile dont les objets magnétiques sont les plus proches de ce barycentre, l'ustensile ainsi sélectionné étant l'ustensile identifié ; ^ si la distance qui sépare l'ustensile le plus proche du barycentre est supérieure à un seuil prédéterminé, alors le procédé comporte : - l'ajout au système actuel d'équations d'au moins une nouvelle équation reliant de nouvelles variables codant la position, l'orientation et l'amplitude du moment magnétique de chaque objet magnétique d'un nouvel ustensile à chaque mesure d'un magnétomètre tri-axes du réseau, pour obtenir un nouveau système d'équations, puis, - l'utilisation de ce nouveau système d'équations à la place du système actuel d'équations lors des itérations suivantes des étapes a) et b), le nouveau système d'équations devenant ainsi, pour les itérations suivantes des étapes a) et b), le système actuel d'équations ; ^ si l'amplitude de la perturbation franchit un seuil prédéterminé, le procédé comporte l'identification de l'ustensile qui a été déplacé parmi les ustensiles contenus dans une liste d'ustensiles immobiles, cette liste d'ustensiles immobiles contenant pour chaque objet magnétique de chaque ustensile immobile : sa position, son orientation et l'amplitude de son moment magnétique, cette identification de l'ustensile comportant : e) l'estimation du champ magnétique mesuré par chaque magnétomètre en présence d'un seul ustensile immobile et en absence de tout autre ustensile, à partir de la position, de l'orientation et de l'amplitude du moment magnétique de chaque objet magnétique de cet ustensile immobile enregistrées dans la liste d'ustensiles immobiles, f) le calcul d'une corrélation entre les mesures estimées lors de l'étape e) et les écarts entre les mesures réalisées par les magnétomètres lors de l'étape a) et les mesures établies pour ces mêmes magnétomètres lors de l'étape c), et g) la sélection dans la liste d'ustensiles immobiles de celui qui comporte le ou les objets magnétiques qui présentent la corrélation, calculée lors de l'étape f), la plus forte, l'ustensile ainsi sélectionné étant l'ustensile identifié ; ^ si la corrélation la plus forte est inférieure à un seuil prédéterminé, alors le procédé comporte : - l'ajout au système actuel d'équations d'au moins une nouvelle équation reliant de nouvelles variables codant la position, l'orientation et l'amplitude du moment magnétique de chaque objet magnétique d'un nouvel ustensile à chaque mesure d'un magnétomètre tri-axes du réseau, pour obtenir un nouveau système d'équations, puis, - l'utilisation de ce nouveau système d'équations à la place du système actuel d'équations lors des itérations suivantes des étapes a) et b), le nouveau système d'équations devenant ainsi, pour les itérations suivantes des étapes a) et b), le système actuel d'équations ; ^ le procédé comporte la détection d'un ustensile immobile, et en réponse, l'ajout de cet ustensile dans la liste des ustensiles immobiles et le remplacement dans le système actuel d'équations des variables codant la position, l'orientation et l'amplitude du moment magnétique de chaque objet magnétique de cet ustensile par des valeurs numériques constantes égales à la position actuelle, l'orientation actuelle et l'amplitude actuelle du moment magnétique de chaque objet magnétique de cet ustensile pour obtenir un nouveau système d'équations dans lequel le nombre de variables dont la valeur doit être estimée est plus petit que dans le système actuel d'équations, puis - l'utilisation du nouveau système d'équations à la place du système actuel d'équations lors des itérations suivantes des étapes a) et b), le nouveau système d'équations devenant ainsi pour les itérations suivantes des étapes a) et b), le système actuel d'équations, puis - en réponse à l'identification de cet ustensile lors de l'étape d'identification, le remplacement dans le système actuel d'équations des valeurs numériques constantes codant la position et l'orientation de chaque objet magnétique de cet ustensile par des variables pour obtenir un nouveau système d'équations dans lequel le nombre de variables, dont la valeur doit être estimée, est plus grand que dans le système actuel d'équations et l'utilisation du nouveau système d'équations à la place du système actuel d'équations lors des itérations suivantes des étapes a) et b) ; ^ le procédé comporte la suppression, dans le système actuel d'équations, des variables codant la position, l'orientation et l'amplitude du moment magnétique de chaque objet magnétique d'un ustensile si la distance estimée entre cet ustensile et le réseau de magnétomètres dépasse un seuil prédéterminé pour obtenir un nouveau système d'équations, puis l'utilisation de ce nouveau système d'équations à la place du système actuel d'équations lors des itérations suivantes des étapes a) et b), le nouveau système d'équations devenant ainsi, pour les itérations suivantes des étapes a) et b), le système actuel d'équations. [0012] Ces modes de réalisation du procédé présentent en outre les avantages suivants : - calculer le barycentre de la perturbation permet d'identifier simplement, sans exécution de l'étape b), l'ustensile qui a été déplacé parmi l'ensemble des ustensiles contenus dans la liste d'ustensiles immobiles ; - déduire la présence d'un nouvel ustensile à partir de la distance qui sépare les ustensiles immobiles du barycentre de la perturbation permet d'adapter automatiquement le système d'équations en réponse à l'ajout d'un nouvel ustensile sans avoir à exécuter l'étape b) ; - identifier l'ustensile qui a été déplacé parmi ceux présents dans la liste d'ustensiles immobiles à l'aide d'une corrélation entre le champ magnétique généré par les objets magnétiques de cet ustensile et la perturbation permet d'identifier simplement, sans exécuter l'étape b), l'ustensile qui a été déplacé ; - déduire du fait que la corrélation entre le champ magnétique généré par les objets magnétiques des ustensiles et la perturbation actuelle sont toutes inférieures à un seuil prédéterminé qu'un nouvel ustensile a été ajouté, permet d'adapter dynamiquement et automatiquement le système d'équations en réponse à l'ajout d'un nouvel ustensile ; - remplacer les variables du système actuel d'équations codant les positions, orientations et amplitudes des objets magnétiques des ustensiles immobiles par des constantes puis, à chaque fois qu'un ustensile est de nouveau déplacé, réaliser l'opération inverse permet d'adapter automatiquement le système d'équations en fonction du déplacement des ustensiles réalisés et donc d'augmenter la précision de la localisation de l'ustensile déplacé ; - supprimer du système d'équations les variables codant la position, l'orientation et l'amplitude des objets magnétiques des ustensiles qui sont éloignés d'une distance prédéterminée des magnétomètres permet de simplifier le système d'équations au fur et à mesure que des ustensiles sont retirés. [0013] L'invention a également pour objet un support d'enregistrement d'informations comportant des instructions pour l'exécution du procédé ci-dessus, lorsque ces instructions sont exécutées par un calculateur électronique. [0014] Enfin, l'invention a également pour objet un dispositif de localisation d'ustensiles mobiles, chaque ustensile comportant au moins un objet magnétique, ce dispositif comportant : - un réseau de magnétomètres comportant N magnétomètres tri-axes reliés mécaniquement les uns aux autres sans aucun degré de liberté pour conserver une distance connue entre chacun des magnétomètres, où N est un nombre entier supérieur ou égal à cinq, - une unité de traitement programmée pour : a) acquérir la mesure par chacun des magnétomètres de l'amplitude du champ magnétique le long de chacun de ses axes de mesure, b) estimer à partir des mesures acquises des magnétomètres, des valeurs de plusieurs variables en résolvant un système actuel d'équations reliant ces variables à chaque mesure d'un magnétomètre tri-axes du réseau, ces variables correspondant aux positions, orientations et amplitudes des moments magnétiques des objets magnétiques des ustensiles à localiser, - itérer les étapes a) et b) à des instants successifs dans le temps pour obtenir les valeurs estimées de chaque variable à ces différents instants successifs, 5 dans lequel l'unité de traitement est également programmée pour : c) établir des mesures des magnétomètres lorsqu'ils sont uniquement en présence d'un champ magnétique de référence généré seulement par les objets magnétiques de tous les ustensiles immobiles présents devant le réseau de magnétomètres, et d) calculer l'amplitude d'une perturbation à partir des écarts entre les mesures 10 réalisées lors de l'étape a) et les mesures établies lors de l'étape c), et - si l'amplitude de la perturbation franchit un seuil prédéterminé, exécuter l'étape b), sinon omettre l'exécution de l'étape b) et retourner à l'étape a) sans avoir exécuté l'étape b). [0015] L'invention sera mieux comprise à la lecture de la description qui va suivre, 15 donnée uniquement à titre d'exemple non limitatif et faite en se référant aux dessins, sur lesquels : - la figure 1 est une illustration schématique d'une interface homme-machine permettant de commander un appareil électrique. - la figure 2 est une illustration partielle, en vue de dessus, d'un réseau de 20 magnétomètres mis en oeuvre dans l'interface homme-machine de la figure 1 ; - la figure 3 est une illustration schématique d'une base de données utilisée dans l'interface de la figure 1 ; - la figure 4 est un organigramme d'un procédé de commande d'un appareil électrique à l'aide de l'interface homme-machine de la figure 1 ; 25 - la figure 5 est un organigramme partiel d'un autre procédé de commande d'un appareil électrique à l'aide de l'interface homme-machine de la figure 1. [0016] Dans ces figures, les mêmes références sont utilisées pour désigner les mêmes éléments. [0017] Dans la suite de cette description, les caractéristiques et fonctions bien 30 connues de l'homme du métier ne sont pas décrites en détails. [0018] La figure 1 représente une interface homme-machine 2 permettant de commander un appareil électrique 4. Ici, l'appareil électrique comporte un écran et une unité 5 de commande capable de commander l'affichage d'une image sur cet écran. 35 [0019] Par la suite, le fonctionnement de l'interface 2 est illustré dans le cas où l'unité 5 est une console de jeux vidéo. Par exemple, l'unité 5 commande le déplacement de pions sur l'écran. Pour simplifier l'illustration, seul un pion 6 est représenté. Toutefois, l'interface 2 est utilisable dans de nombreuses autres applications comme cela est décrit à la fin de cette description. [0020] L'interface 2 comporte plusieurs ustensiles déplaçables directement à la main par un être humain, dénommé « utilisateur » par la suite. Chacun de ces ustensiles comporte au moins un objet magnétique. Ici chaque objet magnétique est un aimant permanent. Pour simplifier la figure 1, seul un ustensile 10 est représenté sur cette 5 figure. [0021] Chaque ustensile est librement déplaçable, directement par la main de l'utilisateur, dans un repère orthogonal Rm fixé sans aucun degré de liberté à un dispositif 12. Ici, les directions X et Y du repère Rm sont horizontales et la direction Z est verticale. A cet effet, chaque ustensile pèse moins d'un kilo et, de préférence, 10 moins de 200g ou 100g. Les dimensions de chaque ustensile sont suffisamment réduites pour qu'il puisse être saisi et déplacé par une seule main de l'utilisateur. [0022] Dans ce mode de réalisation, l'ustensile 10 comprend un corps 14 et un aimant permanent 16. Par exemple, le corps 14 s'étend le long d'un axe longitudinal 17. Le corps 14 présente dans sa partie inférieure un pied apte à immobiliser, dans le 15 repère Rm, le corps 14 lorsque ce pied repose sur une face horizontale. Ici l'ustensile 10 est un pion d'un jeu d'échec. A titre d'illustration, le corps 14 a la forme d'une tour facilement saisissable par l'utilisateur. [0023] Par exemple, le corps 14 est entièrement réalisé dans un matériau non- magnétique, c'est-à-dire un matériau ne présentant pas de propriété magnétique 20 mesurable par le dispositif 12. Ce matériau est par exemple du plastique. L'ustensile 10 est localisé, dans le repère Rm, à partir de la position et l'orientation de l'aimant 16. A cet effet, l'aimant 16 est fixé sans aucun degré de liberté à l'ustensile 10. La direction du moment magnétique de l'aimant 16 est ici confondue avec l'axe longitudinal 17 du corps 14. Sur la figure 1 et les suivantes, la direction du moment 25 magnétique d'un aimant est représentée par une flèche. [0024] L'aimant 16 présente un moment magnétique non nul même en absence de champ magnétique extérieur. Typiquement, dans cette description, le champ magnétique coercitif de chaque aimant est supérieur à 100 A.m-1 ou 500 A.m-1. Par exemple, il est réalisé en matériau ferro- ou ferrimagnétique. La puissance de chaque 30 aimant permanent est typiquement supérieure à 0,01 A.m2 ou 0,1 A.m2. [0025] La plus grande longueur de cet aimant est notée L par la suite. [0026] Pour simplifier la description, on suppose également que la structure des autres ustensiles utilisables dans l'interface 2 est identique à celle de l'ustensile 10 ou ne diffère de celle de l'ustensile 10 que par la forme du corps 14. 35 [0027] Le dispositif 12 permet de localiser un ou plusieurs aimants dans le repère Rm. Par localisation d'un aimant, on entend ici la détermination de la position et de l'orientation de l'aimant 16 dans le repère Rm. La position est définie sans ambiguïté par les valeurs de trois variables, par exemple les coordonnées x, y et z dans le repère Rm. Plus précisément, dans le cas d'un aimant, les variables x, y et z sont les 40 coordonnées du centre géométrique de l'aimant. Le centre de géométrie d'un objet est le barycentre de l'ensemble des points de cet objet en affectant le même poids à chacun de ces points. L'orientation du moment magnétique d'un aimant est définie dans le repère Rm par les valeurs de deux variables ey et 0,. Ici, les variables ey et 0, sont les angles du moment magnétique de l'aimant, respectivement, par rapport aux axes Y et Z du repère Rm. Le dispositif 12 détermine également un sixième variable A. La variable A est l'amplitude du moment magnétique de l'aimant. [0028] Le dispositif 12 comporte à cet effet un réseau de N magnétomètres tri-axes Mi. Sur la figure 1, les traits ondulés verticaux indiquent qu'une partie du dispositif 12 n'a pas été représentée. [0029] Typiquement, N est supérieur à cinq et, de préférence, supérieur à seize ou trente-deux. Ici, N est supérieur ou égal à soixante-quatre. [0030] Dans ce mode de réalisation, les magnétomètres Mi; sont alignés en lignes et en colonnes pour former une matrice. Ici, cette matrice comporte huit lignes et huit colonnes. Les indices i et j identifient, respectivement, la ligne et la colonne de cette 15 matrice à l'intersection de laquelle se trouve le magnétomètre ni. Sur la figure 1, seuls les magnétomètres Mil , -M2, -i4 -MM Mi8 d'une ligne i sont visibles. La position - des magnétomètres Mi; les uns par rapport aux autres est décrite plus en détail en référence à la figure 2. [0031] Chaque magnétomètre Mi; est fixé sans aucun degré de liberté aux autres 20 magnétomètres. A cet effet, les magnétomètres Mi; sont fixés sans aucun degré de liberté sur une face arrière 22 d'une plaque rigide 20. Cette plaque rigide présente une face avant 24 tournée vers l'aimant 16. La plaque 20 est réalisée dans un matériau non magnétique rigide. Par exemple, la plaque 20 est en verre. [0032] Chaque magnétomètre Mi; mesure la direction et l'intensité du champ 25 magnétique généré par les aimants présents devant la face 24. Pour cela, chaque magnétomètre Mi; mesure la norme de la projection orthogonale du champ magnétique au niveau de ce magnétomètre Mi; sur trois axes de mesure de ce magnétomètre. Ici, ces trois axes de mesure sont orthogonaux entre eux. Par exemple, les axes de mesure de chacun des magnétomètres Mi; sont, 30 respectivement, parallèles aux axes X, Y et Z du repère. La sensibilité du magnétomètre Mi; est par exemple de 410-7T. [0033] Chaque magnétomètre Mi; est raccordé par l'intermédiaire d'un bus 28 de transmission d'informations à une unité de traitement 30. [0034] L'unité de traitement 30 est capable, pour chaque aimant, de le localiser dans 35 le repère Rm et d'estimer l'amplitude de son moment magnétique à partir des mesures des magnétomètres Mi. A cet effet, l'unité 30 comporte un calculateur électronique 32 programmable apte à exécuter des instructions enregistrées sur un support d'enregistrement d'informations. L'unité 30 comporte donc aussi une mémoire 34 contenant les instructions nécessaires pour l'exécution par le calculateur 40 32 du procédé de la figure 4 ou 5. L'unité 30 implémente un modèle mathématique Mp associant chaque mesure d'un magnétomètre Mu aux positions, orientations et aux amplitudes des moments magnétiques de P objets magnétiques dans le repère Rm. Chaque modèle Mp se présente sous la forme d'un système d'équations dans lequel un premier jeu de variables représente les positions et orientations des P objets magnétiques ainsi que les amplitudes des moments magnétiques de ces objets. Un second jeu de variables représente les mesures des magnétomètres Mu. Pour obtenir, les positions, orientations et amplitudes des moments magnétiques des P objets magnétiques, les variables du premier jeu sont les inconnues et les valeurs des variables du second jeu sont connues. A l'inverse, pour obtenir une estimations des valeurs mesurées par les magnétomètres Mu, les variables du premier jeu sont connues et les valeurs des variables du second jeu sont les inconnues. Chaque modèle est typiquement construit à partir des équations physiques du magnétisme. Chaque modèle est paramétré par les distances connues entres les magnétomètres Mu. Ici, les objets magnétiques sont les aimants permanents. Par exemple, pour construire ce modèle, chaque aimant permanent est approximé par un dipôle magnétique. Cette approximation n'introduit que très peu d'erreurs si la distance entre l'aimant permanent et le magnétomètre Mu est supérieure à 2L et, de préférence, supérieure à 3L, où L est la plus grande dimension de l'aimant permanent. Typiquement, L est inférieur à 20 cm et, de préférence, inférieur à 10 ou 5cm ou 1cm. [0035] Ici, le modèle MP est non linéaire. Par conséquent, pour localiser les aimants permanents, l'unité 30 le résout en mettant en oeuvre un algorithme d'estimation de sa solution. Par exemple, l'algorithme utilisé est un filtre de Kalman d'ensemble plus connu sous le terme anglais de « Unscented Kalman Filter ». On notera aussi que pour estimer les valeurs mesurées par les magnétomètres à partir de positions, orientations et amplitudes connus d'aimants permanents et à l'aide du modèle Mp, il n'est pas nécessaire de mettre en oeuvre de tels algorithmes d'estimation de la solution. Par exemple, il suffit de remplacer les variables codant les positions, orientations et amplitudes par leurs valeurs respectives. [0036] Étant donné que chaque objet magnétique est caractérisé par trois variables pour connaître sa position, deux variables pour connaître son orientation et une variable pour connaître l'amplitude de son moment magnétique, le nombre maximum d'objets magnétiques simultanément localisables par le réseau de N magnétomètres est inférieur à N/2. Dès lors, la valeur du nombre P est inférieure ou égale à N/2 et, de préférence, inférieure à N/5 ou N/10 ou N/20 pour avoir des mesures redondantes. La redondance des mesures permet d'améliorer la précision de la localisation des objets magnétiques. [0037] De plus, dans ce mode de réalisation, on suppose qu'à chaque instant, il n'y a jamais plus d'un seul ustensile qui est déplacé. Par conséquent, à chaque instant il existe au plus un aimant permanent qui se déplace. Dans ces conditions, initialement, 40 l'unité 30 implémente le modèle M1. Toutefois, dans cette exemple, les modèles M1 jusqu'à M32 sont aussi enregistrés dans la mémoire 34 pour estimer les mesures des magnétomètres Mu. [0038] L'unité 30 est également capable de transmettre une commande à l'appareil 4 par l'intermédiaire d'une interface 36 reliée à cet appareil 4. [0039] La mémoire 34 comporté aussi une base 38 de données dans laquelle sont enregistrés les ustensiles immobiles et leurs positions respectives. Cette base 38 est décrite plus en détail en référence à la figure 3. [0040] La figure 2 représente une partie des magnétomètres M,, du dispositif 12. Ces magnétomètres M,, sont alignés en lignes i parallèles à la direction X. Ces magnétomètres sont également alignés en colonnes j parallèles à la direction Y pour former une matrice. Les lignes i et les colonnes j sont disposées dans l'ordre des indices croissants. [0041] Le centre du magnétomètre M,, se trouve à l'intersection de la ligne i et de la colonne j. Le centre du magnétomètre correspond au point où le champ magnétique 15 est mesuré par ce magnétomètre. Ici, les indices i et j appartiennent à l'intervalle [1 ; 8]. [0042] Les centres de deux magnétomètres M,, et immédiatement consécutifs le long d'une ligne i sont séparés par une distance connue d,,,,,.1. De façon similaire, le centre de deux magnétomètres M,, et M,.,', immédiatement consécutifs le long d'une 20 même colonne j sont séparés par une distance connue d,,',.i. [0043] Dans le cas particulier décrit ici, quelle que soit la ligne i, la distance est la même. Cette distance est donc notée d,. De façon similaire, quelle que soit la colonne j, la distance d,,',., entre deux magnétomètres est la même. Cette distance est donc notée d,. 25 [0044] Ici, les distances d, et d, sont toutes les deux égales à d. [0045] Typiquement, la distance d est inférieure, et de préférence deux fois inférieur, à la plus petite distance qui peut exister entre deux objets magnétiques simultanément présent devant la face 24 lors d'une utilisation normale de l'interface 2. Ici, la distance d est comprise entre 1 et 4 cm lorsque : 30 - la puissance de l'aimant permanent est de 0,5 A.m2, - la sensibilité des magnétomètres est de 4*10-7T, et - le nombre de magnétomètres M,, est de soixante-quatre. [0046] La figure 3 présente plus en détail la base de données 38. Cette base de données 38 contient un identifiant Uk pour chaque ustensile différent immobile dans le 35 repère Rm. L'identifiant Uk permet de distinguer cet ustensile de tous les autres ustensiles simultanément présents devant la face 24. A chaque identifiant U k, la base de données 38 associe la position Pk de son aimant permanent, l'orientation Ok du moment magnétique de cet aimant permanent et l'amplitude Ak du moment magnétique de cet aimant permanent. Pour simplifier l'illustration, seules trois lignes 40 de la base de données 38, respectivement pour les identifiants U1, U2 et U3, sont représentées. Par la suite, étant donné que cette base de données 38 contient la liste des ustensiles immobiles dans leur repère Rm, cette base de données 38 est également appelée liste des ustensiles immobiles. [0047] Le fonctionnement de l'interface 2 va maintenant être décrit en référence au procédé de la figure 4. [0048] Dans ce mode de réalisation, on suppose qu'initialement aucun ustensile n'est présent sur la face 24 du dispositif 12. Par conséquent, la base 38 est vide. Dans cette description, on considère également que le champ magnétique terrestre est négligeable dans les mesures des magnétomètres Mu. Dans l'hypothèse où l'amplitude du champ magnétique terrestre ne serait pas négligeable, on considère alors par la suite que la mesure du champ magnétique terrestre est systématiquement déduite de la mesure actuelle des magnétomètres de sorte que seul l'écart entre ces deux mesures est pris en compte par le calculateur 32 lors de l'exécution du procédé suivant. [0049] Lors d'une étape 100, le magnétomètre Mu mesure le champ magnétique présent. [0050] Ensuite, lors d'une étape 102, le calculateur 32 acquiert les mesures des magnétomètres Mu. [0051] Lors d'une étape 104, le calculateur 32 détecte la présence et, en alternance, l'absence d'une perturbation magnétique. [0052] Pour cela, lors d'une opération 106, le calculateur 32 établit la valeur qui serait mesurée par chaque magnétomètre Mu si le seul champ magnétique présent était le champ magnétique généré par les aimants permanents des ustensiles recensés dans la base 38. Ce champ magnétique est par la suite appelé « champ magnétique de références ». [0053] Ici, le calculateur 32 calcule ces valeurs à partir des positions, orientations et amplitudes des moments magnétiques enregistrées pour chaque aimant permanent des ustensiles contenus dans la base 38. [0054] Par exemple, pour cela, le calculateur 32 sélectionne un modèle Mp où P est égal au nombre d'aimants permanents dont les positions, orientations et amplitudes sont enregistrées dans la base 38. Ensuite, il remplace les variables codant les positions, orientations et amplitudes de chaque dipôle magnétique par les valeurs correspondantes contenues dans la base 38. De cette façon, le calculateur 32 obtient les valeurs que devraient mesurer les magnétomètres Mu si seuls les ustensiles immobiles étaient présents devant la face 24. Par la suite, les valeurs ainsi établies pour chaque axe de mesure de chaque magnétomètre Mu sont appelées « mesures de référence ». [0055] Dans le cas où la base 38 est vide, les mesures de référence sont nulles, car le champ magnétique de référence est également nul. [0056] Par « perturbation » on désigne ici l'ensemble des écarts, pour chaque axe de mesure de chaque magnétomètre Mi;, entre : - la dernière mesure de référence établie lors de l'opération 106, et - la dernière mesure acquise lors de l'étape 102. [0057] Lors d'une opération 108, le calculateur 32 calcule l'amplitude de la perturbation à partir des écarts entre les mesures de référence établies lors de l'étape 106 et les mesures acquises lors de l'étape 102. Par exemple, l'amplitude Ap de la perturbation est calculée à l'aide de la relation suivante : Ap = E(miin - mrefe)2, où : - me est la mesure du magnétomètre Mi; sur son axe n de mesure, 10 - m retijn est la mesure de référence pour le même magnétomètre et sur le même axe de mesure établie lors de l'opération 106, - l'indice n identifie l'axe de mesure dans la relation ci-dessus. [0058] De plus, la somme ci-dessus est calculée pour i et j variant de 1 à N, et pour n variant de 1 à 3. 15 [0059] Lors d'une opération 110, le calculateur 32 compare l'amplitude Ap calculée à un seuil Sp prédéterminé. Si l'amplitude Ap est inférieure ou égale au seuil Sp, alors aucune perturbation n'est détectée, et le procédé retourne à l'étape 100. [0060] Dans le cas contraire, c'est-à-dire si l'amplitude Ap est supérieure au seuil Sp, une perturbation est détectée, ce qui signifie qu'un ustensile a été ajouté ou déplacé 20 par rapport à la situation précédente où tous les ustensiles étaient immobiles. [0061] Dans ce cas, le procédé se poursuit par une étape 120 lors de laquelle le calculateur 32 calcule les coordonnées du barycentre Bp de la perturbation. Plus précisément, lors de l'étape 120, le calculateur 32 calcule les coordonnées du barycentre de l'ensemble des points Pi; en affectant à chaque point Pi; un coefficient 25 de pondération fonction de l'écart entre les mesures de référence et les mesures acquises du magnétomètre ni. Les coordonnées du point Pi; sont égales aux coordonnées du centre géométrique du magnétomètre 1\/4. Le coefficient de pondération Ci; est par exemple donné par la relation suivante : = lm - M refijl )2 + (n1 - M ret1j2)2 (n1 - M ret1j3)2 30 [0062] Lors d'une étape 122, si la base 38 est vide c'est-à-dire qu'elle ne contient actuellement aucun ustensile immobile, alors le calculateur 32 sélectionne automatiquement comme système actuel d'équations le système d'équations du modèle M1. [0063] A l'inverse, si la base 38 contient au moins un identifiant d'au moins un 35 ustensile immobile, le calculateur 32 procède à une étape 124 d'identification de l'ustensile qui a été déplacé parmi l'ensemble des ustensiles identifiés par la base 38. [0064] Pour cela, lors d'une opération 128, pour chaque ustensile identifié dans la base 38, le calculateur 32 calcule la distance Db entre la position de l'aimant permanent de cet ustensile enregistré dans la base 38 et les coordonnées du 40 barycentre Bp. [0065] Lors d'une opération 130, le calculateur sélectionne, dans la base 38, l'ustensile pour lequel la distance Db précédemment calculée est la plus petite. [0066] Lors d'une opération 132, le calculateur 32 compare alors la distance Db calculée pour l'ustensile sélectionné à l'étape 130 à un seuil prédéterminé Sb. Si cette distance Db est inférieure au seuil Sb, le calculateur 32 identifie cet ustensile comme étant celui qui a été déplacé et qui a provoqué la perturbation. Dans ce cas, le procédé se poursuit par une étape 140. [0067] Dans le cas contraire, c'est-à-dire si la distance Db est supérieure au seuil Sb, cela signifie qu'aucun des ustensiles identifiés dans la base 38 n'a été identifié 10 comme étant celui qui a été déplacé pour créer la perturbation. Dans ce dernier cas, le procédé se poursuit par une étape 142. [0068] Lors de l'étape 140, le calculateur 32 supprime dans la base 38 l'identifiant de l'ustensile identifié ainsi que toutes les données associées à cet ustensile identifié. Ensuite, il remplace dans le système actuel d'équations les valeurs numériques 15 constantes qui codaient la position, l'orientation et l'amplitude du moment magnétique de l'aimant permanent de cet ustensile par des variables. On obtient alors un nouveau système d'équations. Ce nouveau système d'équations remplace le système actuel d'équations lors des itérations suivantes de l'étape 150. [0069] Lors de l'étape 142, le calculateur 32 ajoute au système actuel d'équations de 20 nouvelles équations contenant six nouvelles variables x, y, z, ey, e et A codant la position, l'orientation et l'amplitude du moment magnétique d'un aimant permanent d'un nouvel ustensile. Les équations ajoutées relient ces variables aux mesures réalisées par chaque magnétomètre ni. Par exemple, ces équations sont obtenues à partir de la relation suivante : B(r,r') = (p0/47)[3*(r-e)*[m.(r-r')]/ I r-r' 15 - mil r-r' I 3], où 25 - po : la perméabilité magnétique du vide à 4n*10-7 H.m1; - m : vecteur du moment magnétique de l'aimant permanent en A.m2 ; - r' : vecteur position de l'aimant dans le repère Rm; - r : vecteur position du magnétomètre Mi; dans le repère Rm; - (r-r') : vecteur liant le centre de l'aimant au point de mesure ; 30 - B(r,r') : vecteur du champ magnétique au niveau du magnétomètre Mi; en Tesla; - I ... I : fonction retournant la norme d'un vecteur ; - « . » : fonction produit scalaire ; - « *» : multiplication par un nombre réel. [0070] A partir de cette relation et connaissant la position du magnétomètre Mi; dans 35 le repère Rm, on obtient trois équations reliant les variables x, y, z, ey, ez, A aux mesures, respectivement, sur les trois axes de mesure du magnétomètre ni. [0071] Ainsi, à l'issue de l'étape 142, on obtient un nouveau système d'équations qui est utilisé à la place du système actuel d'équations lors des itérations suivantes de l'étape 150. [0072] A l'issue des étapes 122 et 142, le calculateur 32 dispose d'un système d'équations, appelé système actuel d'équations, reliant les variables codant la position, l'orientation et l'amplitude du moment magnétique de l'aimant permanent qui se déplace aux mesures des magnétomètres. [0073] Lors de l'étape 150, le calculateur 32 estime la position, l'orientation et l'amplitude du moment magnétique de l'aimant permanent de l'ustensile qui se déplace à l'aide de ce système actuel d'équations. Pour cela, il résout le système actuel d'équations en utilisant les dernières mesures acquises lors de l'étape 102. Par exemple, la solution de ce système actuel d'équations est estimée en utilisant un filtre de Kalman d'ensemble ou « unscented Kalman filter » en anglais. Lors de l'étape 150, les valeurs des variables à estimer codant la position sont initialisées en les prenant égales aux coordonnées du barycentre Bp calculées lors de l'étape 120. Ceci permet d'accélérer la convergence du procédé vers l'estimation des valeurs qui résolvent ce système actuel d'équations. [0074] Ensuite, lors d'une étape 152, le calculateur 32 commande l'appareil 4 en fonction de la position, de l'orientation et de l'amplitude estimées lors de l'étape 150. Par exemple, cette position, cette orientation et cette amplitude estimées sont transmises à l'unité 5 qui commande en réponse un déplacement d'une représentation graphique de l'ustensile sur l'écran. [0075] A l'issue de l'étape 152, le procédé retourne à l'étape 100. [0076] En parallèle des étapes 104 à 152, lors d'une étape 160, le calculateur 32 détecte si un ustensile s'immobilise, c'est-à-dire que cet ustensile ne se déplace plus dans le repère Rm. Pour cela, il utilise les précédentes estimations des positions et orientations de l'aimant permanent de cet ustensile. [0077] Par exemple, lors d'une opération 162, le calculateur 32 calcule la variance Vd des X dernières positions estimées pour l'aimant permanent de cet ustensile, où X est un nombre entier qui, par exemple, dépend de la fréquence Fedi de rafraîchissement des mesures des magnétomètres. On peut prendre un nombre d'échantillons correspondant à plus d'une seconde, c'est-à-dire un nombre X supérieur à Fedi x 1 secondes. Par exemple, X est supérieur à 10, 100 ou 1000. [0078] Lors d'une opération 164, le calculateur 32 compare la variance Vd à un seuil prédéterminé Sv. [0079] Si la variance Vd est supérieure à ce seuil Sv, alors le calculateur 32 considère que l'ustensile n'est pas immobile et le procédé retourne à l'étape 100. posol Dans le cas contraire, c'est-à-dire si la variance de l'ustensile est inférieure au seuil Sv, le calculateur 32 considère alors que cet ustensile est maintenant immobile. Dans ce cas, il procède à une étape 166 lors de laquelle le calculateur 32 ajoute l'identifiant de cet ustensile dans la base 38 et enregistre, associé à cet identifiant, les dernières position, orientation et amplitude du moment magnétique estimées pour l'aimant permanent de cet ustensile. [00811 De plus, lors de l'étape 166, le calculateur 32 remplace, dans le système actuel d'équations, les variables codant la position, l'orientation et l'amplitude du moment magnétique de l'aimant permanent de cet ustensile par des valeurs numériques constantes. Ces valeurs numériques constantes sont égales aux dernières position, orientation et amplitude estimées pour cet aimant permanent lors de la dernière exécution de l'étape 150. Ainsi, après l'exécution de l'étape 166, le système actuel d'équations comporte moins de variables à estimer qu'avant. [0082] Après l'étape 166, le procédé retourne à l'étape 100. [0083] En parallèle des étapes 104 à 152, lors d'une étape 170, le calculateur 32 détecte également qu'un ustensile a été retiré. Pour cela, lors d'une opération 172, le calculateur 32 calcule la distance Di, qui sépare l'aimant permanent de l'ustensile du magnétomètre Mu le plus proche. Cette distance Di, est calculée à partir de la dernière position estimée pour l'aimant permanent de cet ustensile. [0084] Ensuite, lors d'une opération 174, le calculateur 32 compare la distance Di, à un seuil prédéterminé Su. Si la distance Di, est inférieure ou égale au seuil Su, alors le calculateur 32 conserve les variables codant les position, orientation et amplitude de l'aimant de cet ustensile dans le système actuel d'équations. [0085] Dans le cas contraire, c'est-à-dire si la distance Di, est supérieure au seuil Su, lors d'une opération 176, le calculateur 32 supprime du système actuel d'équations les équations faisant intervenir les variables codant la position, l'orientation et l'amplitude de l'aimant permanent de cet ustensile. Le calculateur supprime ainsi la contribution de l'ustensile qui a été supprimé. Il n'est ensuite plus pris en compte dans le calcul. [0086] A l'issue de l'étape 170, le procédé retourne à l'étape 100. [0087] La mise en oeuvre de ce procédé va maintenant être illustrée dans le cas particulier où les ustensiles sont les pions d'un jeu d'échecs et dans le cas où un échiquier a été dessiné sur la face 24. Initialement, on suppose qu'aucun pion n'est présent sur la face 24. [0088] Lorsqu'un joueur vient disposer le premier pion sur la face 24, le calculateur détecte une perturbation. Étant donné qu'il s'agit du premier pion déposé, la base 38 est vide. Par conséquent, le calculateur 32 localise ce premier pion, tant qu'il est déplacé, en utilisant le système d'équations du modèle M1. [0089] Ensuite, le joueur dépose le premier pion sur la face 24 et le laisse immobile à cet emplacement. Le calculateur 32 détecte alors automatiquement que ce pion est immobile. En réponse, il ajoute un identifiant de ce pion dans la base 38 et les variables codant la position, l'orientation et l'amplitude de l'aimant permanent de ce pion sont remplacées par des valeurs constantes dans le système actuel d'équations. [0090] Le joueur dépose alors un second pion sur la face 24. Le procédé décrit ci- dessus conduit alors cette fois le calculateur 32 à ajouter l'identifiant d'un second 40 ustensile dans la base 38 et à remplacer, dans le système actuel d'équations, les variables codant les position, orientation et amplitude de l'aimant permanent de ce deuxième pion par des constantes. A ce stade, le système actuel d'équations comporte donc des constantes codant les positions, orientations et amplitudes à la fois des premier et deuxième pions. [0091] Ce qui vient d'être décrit précédemment est réitéré pour chaque pion déposé, les uns après les autres sur la face 24. Une fois que l'ensemble des pions ont été déposés sur la face 24, avant que les joueurs ne commencent à jouer, la base 38 comporte alors les identifiants de trente-deux pions. [0092] Ensuite, à chaque fois qu'un joueur déplace un pion, le calculateur détecte une perturbation puis identifie le pion déplacé dans la base 38. Tant qu'un pion est déplacé, le calculateur 32 le localise et, par exemple, commande l'appareil 4 pour déplacer sur l'écran un pion correspondant. Quand l'utilisateur repose le pion sur une autre case, le calculateur 32 détecte que ce pion s'immobilise. Il ajoute donc à nouveau ce pion dans la base 38 des ustensiles immobiles et enregistre dans cette même base sa nouvelle position. [0093] Quand un joueur retire un pion de l'échiquier, le calculateur 32 détecte que le pion s'éloigne de l'échiquier et que la distance qui le sépare du magnétomètre Mu le plus proche devient supérieure au seuil Su. En réponse, il supprime les variables codant la position, l'orientation et l'amplitude de l'aimant permanent de ce pion dans le système actuel d'équations. En même temps, il fait, par exemple, disparaître sur l'écran de l'appareil 4 la représentation graphique de ce pion. [0094] Ainsi, le procédé mis en oeuvre par le calculateur 32 résout seulement un système d'équations avec au plus six variables à estimer pour localiser un pion alors que 32 pions peuvent être simultanément présents sur l'échiquier. Pour cela, le calculateur 32 exploite le fait que tous les pions sont déplacés uniquement les uns après les autres. Le fait d'utiliser un système d'équations ne comportant que peu de variables pour estimer la position d'un pion qui se déplace, permet d'augmenter la précision de la localisation du déplacement de ce pion. Cela limite également le nombre d'opérations à exécuter par le calculateur 32, ce qui décroît sa consommation électrique. [0095] La figure 5 représente un autre mode de réalisation du procédé de la figure 4. Le procédé de la figure 5 est identique au procédé de la figure 4, sauf que l'étape 124 est remplacée par une étape 180. Pour simplifier la figure, seule l'étape 180 a été représentée. [0096] Plus précisément, l'étape 180 consiste à identifier l'ustensile dont un identifiant est contenu dans la base 38 et qui a été déplacé. Pour cela, ici, le calculateur 32 utilise une corrélation entre la perturbation et le champ magnétique généré par les ustensiles identifiés par la base 38. [0097] Par exemple, lors d'une opération 182, pour chaque ustensile identifié dans la base 38, le calculateur 32 estime le champ magnétique qui serait mesuré par chaque magnétomètre Mu si cet ustensile était le seul à être présent devant la face 24. Par exemple, à cet effet, le calculateur 32 utilise le système d'équations du modèle M1 et les position, orientation et amplitude de l'aimant permanent de cet ustensile enregistrées dans la base 38. [0098] Ensuite, lors d'une opération 184, le calculateur 32 calcule une valeur représentative de la corrélation Cp qui existe entre les mesures calculées lors de l'opération 182 et la perturbation. [0099] Lors d'une opération 186, le calculateur 32 sélectionne l'ustensile pour lequel la corrélation est maximale. [cm oo] Enfin, lors d'une opération 188, le calculateur 32 compare la corrélation Cp calculée pour l'ustensile sélectionné lors de l'opération 186 à un seuil prédéterminé Sc. [00101] Si la corrélation Cp est supérieure au seuil Sc, l'ustensile qui se déplace est identifié dans la base 38 et le procédé se poursuit par l'étape 140. [00102] A l'inverse, si la corrélation Cp est inférieure au seuil Sc, l'ustensile qui a causé la perturbation n'est pas identifié comme étant l'un de ceux de la base 38 et le procédé se poursuit par l'étape 142. [00103] De nombreux autres modes de réalisation du dispositif 12 sont possibles. Par exemple, la face avant 24 n'est pas nécessairement plane. Par exemple, en variante, il s'agit d'une forme en relief. Typiquement, cette forme en relief peut comporter des logements en creux destinés à recevoir et à immobiliser les ustensiles dans ces logements. [00104] Il est possible de limiter le nombre d'emplacements où un ustensile peut être disposé par rapport à la face avant 24. Pour cela, l'interface 2 comporte plusieurs emplacements, équipé chacun de détrompeur propre à empêcher la fixation à cet emplacement d'un certain nombre d'ustensiles et, au contraire, à autoriser la fixation à cet emplacement d'autres ustensiles différents. [00105] D'autres ustensiles que ceux précédemment décrits sont réalisables et utilisables dans l'interface 2. Par exemple, le corps 14 de l'ustensile peut être omis. 30 Dans ce cas, un ustensile se réduit à un aimant permanent. [00106] L'approximation utilisée pour construire le modèle Mp peut aussi être une approximation quaternaire ou supérieure, c'est-à-dire que les équations du magnétisme sont approximées à un ordre supérieur à celui correspondant à l'approximation dipolaire. 35 [00107] Les magnétomètres du réseau de magnétomètres ne sont pas nécessairement rangés en colonnes et en lignes. Ils peuvent être agencés aussi selon d'autres motifs. Par exemple, les magnétomètres sont disposés sur chaque sommet de chaque maille triangulaire ou hexagonale d'un maillage d'un plan. [cm 08] La disposition des magnétomètres les uns par rapport aux autres peut aussi 40 être aléatoire ou non régulière. Ainsi, la distance entre deux magnétomètres immédiatement consécutifs dans le réseau n'est pas nécessairement la même pour toutes les paires de deux magnétomètres immédiatement consécutifs. Par exemple, la densité de magnétomètres dans une zone donnée du réseau peut être plus importante qu'ailleurs. Augmenter la densité dans une zone donnée peut permettre d'augmenter la précision de la mesure dans cette zone. [00109] Le réseau de magnétomètres peut aussi s'étendre dans trois directions non colinéaires de l'espace. Dans ces conditions, les magnétomètres sont répartis à l'intérieur d'un volume à trois dimensions. [00110] Le nombre N de magnétomètres peut également être supérieur ou égal à 10 soixante-quatre ou quatre-vingt dix. [00111] Tous les magnétomètres du réseau de magnétomètres ne sont pas nécessairement identiques les uns aux autres. En variante, les magnétomètres n'ont pas tous la même sensibilité. [00112] L'appareil 4 peut être remplacé par tout type d'appareil électrique qui doit être 15 commandé en réponse à une action d'un être humain. Par exemple, l'appareil commandé peut être un robot, une machine-outil ou autre. [00113] Dans tous les modes de réalisation décrits ici, l'aimant permanent peut être remplacé par un objet magnétique non alimenté en permanence et qui acquiert un moment magnétique en présent d'un champ magnétique continu extérieur, tel que le 20 champ magnétique terrestre. Par exemple, l'aimant permanent est remplacé par une pièce en matériau magnétique doux. Un matériau magnétique est considéré comme doux si son champ magnétique coercitif est inférieur à 10 ou 1 A.m-1. Une telle pièce présente un moment magnétique créé par l'interaction entre le champ magnétique terrestre et la pièce en matériau magnétique doux. 25 [00114] De nombreux autres modes de réalisation du procédé de la figure 4 ou 5 sont possibles. Par exemple, le nombre de mesures utilisées pour calculer la variance Vd peut être fonction de la durée d'une fenêtre temporelle glissante. Typiquement, seules les mesures réalisées à l'intérieur de cette fenêtre temporelle glissante sont prises en compte pour calculer la variance. 30 [00115] Il existe aussi d'autres moyens que la variance pour déterminer si un ustensile est immobile. Par exemple on peut incrémenter un compteur tant que l'ustensile ne s'est pas ou peu déplacé. Dès que ce compteur atteint un seuil prédéterminé Scmax l'ustensile est dit immobile. Plus précisément, initialement, l'unité 30 mémorise la position actuelle Pi, de l'ustensile comme étant la position de 35 référence Ref. A chaque itération suivante on vérifie que l'ustensile s'est peu déplacé par rapport à la position de référence Ref. Pour cela, par exemple, l'unité 30 vérifie que l'écart entre la position actuelle P, et la position Pref est inférieur à un seuil prédéterminé S,mmo. Dans l'affirmative, l'unité 30 incrémente le compteur. Dans le cas contraire, l'unité 30 réinitialise la position de référence Pr& en la prenant égale à la 40 position actuelle P. L'unité 30 réinitialise en même temps le compteur. On reproduit le même schéma à chaque itération et pour chaque ustensile. Cette méthode est plus facile a implémenter en embarquer et requiert moins de puissance de calcul que le calcul de la variance Vd. [00116] Lors de l'opération 164, le seuil Sv peut être spécifique à la variable estimée. Par exemple, un premier seuil So est utilisé lorsque la variable code une position d'un premier type d'aimant permanent et un second seuil Sv2 est utilisé lorsque la variable code la position d'un second type d'aimant permanent, par exemple, moins puissant. L'ustensile est dit immobile si la variance sur chaque variable de chacun de ses aimants permanents est inférieure au seuil correspondant. [00117] Une fois estimée, l'amplitude des moments magnétiques des aimants permanents peut être considérée comme une constante. Ainsi, le calculateur 32 ne remplace jamais, dans le système actuel d'équations, la valeur numérique codant l'amplitude d'un moment magnétique d'un aimant permanent par une variable. Autrement dit, seules les valeurs codant la position et l'orientation de l'aimant permanent sont remplacées par des variables si une nouvelle estimation de cette position et de cette orientation de l'aimant permanent doit être réalisée. L'amplitude des aimants permanents de chaque ustensile peut aussi être une donnée connue dès l'initialisation du dispositif 12 de sorte qu'il n'est pas nécessaire d'estimer cette donnée à partir de mesures des magnétomètres Mu. [00118] Les position, orientation et amplitude enregistrées dans la base 38 ne sont pas nécessairement prises égales à la dernière estimation de ces valeurs. En variante, ces valeurs peuvent être calculées à partir de plusieurs estimations de ces valeurs obtenues lors des X précédentes estimations. Par exemple, ces valeurs sont calculées en les prenant égales à la moyenne des positions et orientations estimées pendant les X dernières itérations de l'étape 150. [00119] D'autres méthodes sont possibles pour déterminer qu'un ustensile est immobile. Par exemple, la hauteur de l'ustensile par rapport à la face 24 est calculée à partir de la position estimée de l'aimant permanent de cet ustensile. Si cette hauteur est inférieure à un seuil prédéterminé, cela signifie que l'ustensile a été posé sur la face 24 et donc qu'il est immobile. Ce mode de réalisation est particulièrement pertinent lorsque le support est équipé de logements qui ne permettent pas un déplacement horizontal lorsque l'ustensile y est inséré. [00120] D'autres méthodes sont possibles pour identifier l'ustensile de la base 38 qui a été déplacé. Par exemple, il est possible de mettre en oeuvre une méthode appelée « méthode de la convergence partielle » par la suite. Cette méthode consiste à identifier les variables du ou des ustensiles de la base 38 qui, si elles étaient réintroduites dans le système actuel d'équations permettraient de minimiser l'erreur sur les positions et amplitudes estimées. Pour cela, le calculateur 32 procède comme suit : Opération 1) : Pour le premier ustensile de la base 38, le calculateur 32 réintroduit dans le système actuel d'équations les variables codant la position et l'amplitude des aimants permanents de cet ustensile. Lors de cette opération, seule les variables de ce premier ustensile sont réintroduites dans le système d'équations. Les variables des autres ustensiles de la base 38 ne sont pas réintroduites dans le système actuel d'équations. Ainsi, à l'issue de l'opération 1), on obtient un nouveau système d'équations contenant uniquement les variables du premier ustensile. Opération 2) : Ensuite, le calculateur 32 résout ce nouveau système d'équations pour obtenir une estimation de la position et de l'amplitude de chaque aimant permanent 10 du premier ustensile. Opération 3) : Ensuite, le calculateur 32 calcule une erreur Eerr d'estimation représentative de l'écart entre : - les valeurs estimées des mesures des magnétomètres Mu lorsque l'ustensile occupe la position et l'orientation estimées lors de l'opération 2), et 15 - les valeurs réellement mesurées par les magnétomètres Mu au même instant. [00121] Les opérations 1) à 3) sont itérées pour le premier, puis le deuxième jusqu'au dernier ustensile de la base 38. Ensuite, le calculateur 32 sélectionne l'ustensile pour lequel l'erreur d'estimation Eerr calculée lors de l'opération 3) est la plus petite. Si pour cet ustensile sélectionné, l'erreur Eerr est supérieure à un seuil prédéterminé Serr, 20 alors, l'ustensile qui a été déplacé n'est pas identifié dans la base 38. Si à l'inverse, l'erreur d'estimation Eerr est inférieure à ce seuil Serr, l'ustensile qui a provoqué la perturbation est identifié parmi ceux contenus dans la base 38. Par exemple, la valeur du seuil Serr est prise égale à la valeur de l'erreur Eerr calculée à partir des positions et amplitudes estimées lors de la dernière itération de l'étape 150. 25 [00122] En variante, si la base 38 de données est composée d'ustensiles contenant chacun plusieurs aimants fixes à l'intérieur de cet ustensile, les variables réintroduites lors de l'opération 1) correspondent à la position et l'orientation de l'ustensile et non à la position, l'orientation et le moment magnétique de chaque aimant de cet ustensile. [00123] Le calcul des coordonnées du barycentre Bp peut être omis notamment si 30 l'identification de l'ustensile déplacé, parmi ceux dont les identifiants sont contenus dans la base 38, est réalisée sans utiliser les coordonnées de ce barycentre B p. Dans ce cas, lorsqu'un nouvel ustensile est ajouté, les position et orientation de l'aimant permanent de ce nouvel ustensile sont initialisées à des valeurs prédéfinies. Par exemple, ces valeurs prédéfinies correspondent au centre de la face 24. 35 [00124] D'autres méthodes sont possibles pour détecter qu'un ustensile est retiré. Par exemple, le calculateur 32 détecte que l'ustensile s'est immobilisé dans une position qui se trouve en dehors d'une zone prédéfinie du repère Rm, alors, l'ustensile est considéré comme retiré. Toutefois, si l'ustensile reste suffisamment proche des magnétomètres Mu pour que le champ magnétique généré par son aimant permanent 40 puisse influer de façon non négligeable la mesure d'au moins l'un des magnétomètres Mu, alors les constantes codant, dans le système actuel d'équations la position, l'orientation et l'amplitude de l'aimant permanent de cet ustensile ne sont pas éliminées de ce système actuel d'équations. [00125] Un ustensile peut comporter U aimants permanents, où U est un nombre entier strictement supérieur à un. Dans ce cas, cet ustensile est non pas associé à six variables, mais à 6*U variables. Dès lors, quand cet ustensile est ajouté dans la base 38, ce sont les valeurs de ces 6*U variables qui sont enregistrées dans cette base 38. De plus, quand cet ustensile se déplace, le système actuel d'équations comporte non pas six variables à estimer, mais 6*U variables à estimer. [00126] Les procédés ci-dessus ont été décrits dans le contexte particulier où au plus un seul ustensile se déplace à un instant donné. Toutefois, pour détecter la présence et, en alternance, l'absence d'une perturbation, cette hypothèse n'est pas nécessaire. Plus précisément, les opérations décrites ci-dessus pour détecter la présence d'une perturbation fonctionnent aussi dans le cas où la perturbation est causée par le déplacement simultané de deux ustensiles ou plus. Dans le cas où la perturbation est causée par le déplacement simultané de plusieurs ustensiles, les étapes 124 et 180 ne permettent pas d'identifier dans la base 38 les ustensiles simultanément déplacés. Ainsi, dans ce cas, ces étapes 124, 140 et 180 sont omises. L'étape 142 est alors remplacée par une étape lors de laquelle le système actuel d'équations est réinitialisé. Par exemple, le système d'équations est réinitialisé en remplaçant toutes les constantes numériques codant les positions, orientations et amplitudes des aimants des ustensiles immobiles par des variables à estimer. Le système d'équations peut aussi être réinitialisé à un état par défaut prédéterminé. Par exemple, dans le cas précédemment décrit d'un jeu d'échecs, la réinitialisation consiste à remplacer le système actuel d'équations par le système d'équations du modèle M32. Dans le cas où la perturbation peut être causée par le déplacement simultané de plusieurs ustensiles, il est aussi possible d'identifier dans la base 38 les ustensiles qui sont déplacés. Par exemple, pour cela, la méthode de la convergence partielle décrite ci-dessus peut être utilisée. Dans ce cas, tous les ustensiles de la base 38 qui permettent de réduire l'erreur Eerr sont considérés comme étant les ustensiles simultanément déplacés. Le système actuel d'équations est alors remplacé pour les itérations suivantes par un nouveau système d'équations dans lequel les variables codant la position et l'amplitude du moment magnétique de chaque ustensile identifié comme étant en cours de déplacement sont réintroduites. [00127] En variante, initialement, la base 38 n'est pas vide, mais comporte déjà les identifiants et les positions, orientations et amplitudes des aimants permanents de plusieurs ustensiles immobiles. Une telle initialisation de la base 38 peut être réalisée sans utiliser les mesures des magnétomètres Mu, par exemple, si la position initiale de chaque ustensile est considérée comme connue à l'avance. C'est par exemple le cas pour des ustensiles comme des pions d'un jeu d'échecs ou de dame. [00128] L'établissement des mesures du champ magnétique de référence peut consister à acquérir les mesures des différents magnétomètres Mi; à l'instant où tous les ustensiles sont immobiles. [00129] L'amplitude de la perturbation peut être calculée de façons différentes. Par 5 exemple, l'amplitude Ap est donnée par la relation suivante : Ap = Max[(mijn Mretljn)2 (rnij2 mrefij2)2 (rnij3 mrefij3)2] calculée pour tous les indices i et j variant de 1 à N, où Max[...] est la fonction qui retourne le maximum des termes contenus entre crochets. [00130] Lors de l'identification de l'ustensile immobile qui a été déplacé, l'opération 132 ou 188 peut être omise. Dans ce cas, à chaque itération de l'étape 124 ou 180, 10 un ustensile est systématiquement identifié comme étant celui qui a été déplacé parmi ceux contenus dans la base 38. [00131] De nombreuses méthodes différentes sont utilisables pour déterminer la position et l'orientation de l'objet magnétique. Par exemple, la méthode décrite dans US6269324 est utilisable. Ces méthodes n'utilisent pas nécessairement un filtre de 15 Kalman. Par exemple, les méthodes décrites dans US2002/171427A1 ou US6263230B1 sont possibles. REVENDICATIONS1. Procédé de localisation d'ustensiles mobiles présentés devant un réseau de magnétomètres comportant N magnétomètres tri-axes (Mu) reliés mécaniquement les uns aux autres sans aucun degré de liberté pour conserver une distance connue entre chacun des magnétomètres, chaque ustensile comportant au moins un objet magnétique, où N est un nombre entier supérieur ou égal à cinq, ce procédé comportant : a) la mesure (100) par chacun des magnétomètres de l'amplitude du champ 10 magnétique le long de chacun de ses axes de mesure, b) l'estimation (150) à partir des mesures des magnétomètres, des valeurs de plusieurs variables en résolvant un système actuel d'équations reliant ces variables à chaque mesure d'un magnétomètre tri-axes du réseau, ces variables correspondant aux positions, orientations et amplitudes des moments magnétiques des objets 15 magnétiques des ustensiles à localiser, - l'itération des étapes a) et b) à des instants successifs dans le temps pour obtenir les valeurs estimées de chaque variable à ces différents instants successifs, caractérisé en ce que le procédé comporte : c) l'établissement (106) des mesures des magnétomètres lorsqu'ils sont uniquement 20 en présence d'un champ magnétique de référence généré seulement par les objets magnétiques de tous les ustensiles immobiles présents devant le réseau de magnétomètres, et d) le calcul (108) de l'amplitude d'une perturbation à partir des écarts entre les mesures réalisées lors de l'étape a) et les mesures établies lors de l'étape c), et 25 - si l'amplitude de la perturbation franchit un seuil prédéterminé, l'étape b) est exécutée, sinon l'exécution de l'étape b) est omise et le procédé retourne à l'étape a) sans avoir exécuté l'étape b). 2. Procédé selon la revendication 1, dans lequel si l'amplitude de la perturbation 30 franchit un seuil prédéterminé, le procédé comporte l'identification de l'ustensile qui a été déplacé parmi les ustensiles contenus dans une liste d'ustensiles immobiles, cette liste d'ustensiles immobiles contenant pour chaque objet magnétique de chaque ustensile immobile : sa position, son orientation et l'amplitude de son moment magnétique, cette identification de l'ustensile comportant : 35 - le calcul (120) du barycentre de la perturbation, ce barycentre étant le barycentre des positions des différents magnétomètres du réseau, dans lequel la position de chaque magnétomètre est pondérée par l'écart entre les mesures réalisées par cemagnétomètre lors de l'étape a) et les mesures établies pour ce magnétomètre lors de l'étape c), et - la sélection (130), dans la liste d'ustensiles immobiles, de l'ustensile dont les objets magnétiques sont les plus proches de ce barycentre, l'ustensile ainsi sélectionné 5 étant l'ustensile identifié. 3. Procédé selon la revendication 2, dans lequel si la distance qui sépare l'ustensile le plus proche du barycentre est supérieure à un seuil prédéterminé, alors le procédé comporte : 10 - l'ajout (142) au système actuel d'équations d'au moins une nouvelle équation reliant de nouvelles variables codant la position, l'orientation et l'amplitude du moment magnétique de chaque objet magnétique d'un nouvel ustensile à chaque mesure d'un magnétomètre tri-axes du réseau, pour obtenir un nouveau système d'équations, puis, 15 - l'utilisation de ce nouveau système d'équations à la place du système actuel d'équations lors des itérations suivantes des étapes a) et b), le nouveau système d'équations devenant ainsi, pour les itérations suivantes des étapes a) et b), le système actuel d'équations. 20 4. Procédé selon l'une quelconque des revendications précédentes, dans lequel si l'amplitude de la perturbation franchit un seuil prédéterminé, le procédé comporte l'identification (180) de l'ustensile qui a été déplacé parmi les ustensiles contenus dans une liste d'ustensiles immobiles, cette liste d'ustensiles immobiles contenant pour chaque objet magnétique de chaque ustensile immobile : sa position, son 25 orientation et l'amplitude de son moment magnétique, cette identification de l'ustensile comportant : e) l'estimation (182) du champ magnétique mesuré par chaque magnétomètre en présence d'un seul ustensile immobile et en absence de tout autre ustensile, à partir de la position, de l'orientation et de l'amplitude du moment magnétique de chaque 30 objet magnétique de cet ustensile immobile enregistrées dans la liste d'ustensiles immobiles, f) le calcul (184) d'une corrélation entre les mesures estimées lors de l'étape e) et les écarts entre les mesures réalisées par les magnétomètres lors de l'étape a) et les mesures établies pour ces mêmes magnétomètres lors de l'étape c), et 35 g) la sélection (186) dans la liste d'ustensiles immobiles de celui qui comporte le ou les objets magnétiques qui présentent la corrélation, calculée lors de l'étape f), la plus forte, l'ustensile ainsi sélectionné étant l'ustensile identifié. 5. Procédé selon la revendication 4, dans lequel si la corrélation la plus forte est inférieure à un seuil prédéterminé, alors le procédé comporte : - l'ajout (142) au système actuel d'équations d'au moins une nouvelle équation reliant de nouvelles variables codant la position, l'orientation et l'amplitude du moment magnétique de chaque objet magnétique d'un nouvel ustensile à chaque mesure d'un magnétomètre tri-axes du réseau, pour obtenir un nouveau système d'équations, puis, - l'utilisation de ce nouveau système d'équations à la place du système actuel d'équations lors des itérations suivantes des étapes a) et b), le nouveau système 10 d'équations devenant ainsi, pour les itérations suivantes des étapes a) et b), le système actuel d'équations. 6. Procédé selon la revendication 2 ou 4, dans lequel le procédé comporte la détection (160) d'un ustensile immobile, et en réponse, l'ajout de cet ustensile dans la 15 liste des ustensiles immobiles et le remplacement (166) dans le système actuel d'équations des variables codant la position, l'orientation et l'amplitude du moment magnétique de chaque objet magnétique de cet ustensile par des valeurs numériques constantes égales à la position actuelle, l'orientation actuelle et l'amplitude actuelle du moment magnétique de chaque objet magnétique de cet ustensile pour obtenir un 20 nouveau système d'équations dans lequel le nombre de variables dont la valeur doit être estimée est plus petit que dans le système actuel d'équations, puis - l'utilisation du nouveau système d'équations à la place du système actuel d'équations lors des itérations suivantes des étapes a) et b), le nouveau système d'équations devenant ainsi pour les itérations suivantes des étapes a) et b), le 25 système actuel d'équations, puis - en réponse à l'identification de cet ustensile lors de l'étape d'identification, le remplacement (140) dans le système actuel d'équations des valeurs numériques constantes codant la position et l'orientation de chaque objet magnétique de cet ustensile par des variables pour obtenir un nouveau système d'équations dans lequel 30 le nombre de variables, dont la valeur doit être estimée, est plus grand que dans le système actuel d'équations et l'utilisation du nouveau système d'équations à la place du système actuel d'équations lors des itérations suivantes des étapes a) et b). 7. Procédé selon l'une quelconque des revendications précédentes, dans lequel le 35 procédé comporte la suppression (176), dans le système actuel d'équations, des variables codant la position, l'orientation et l'amplitude du moment magnétique de chaque objet magnétique d'un ustensile si la distance estimée entre cet ustensile et le réseau de magnétomètres dépasse un seuil prédéterminé pour obtenir un nouveausystème d'équations, puis l'utilisation de ce nouveau système d'équations à la place du système actuel d'équations lors des itérations suivantes des étapes a) et b), le nouveau système d'équations devenant ainsi, pour les itérations suivantes des étapes a) et b), le système actuel d'équations. 8. Support d'enregistrement d'informations, 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 un calculateur électronique. 9. Dispositif de localisation d'ustensiles mobiles, chaque ustensile comportant au moins un objet magnétique, ce dispositif comportant : - un réseau de magnétomètres comportant N magnétomètres tri-axes (M,J) reliés mécaniquement les uns aux autres sans aucun degré de liberté pour conserver une 15 distance connue entre chacun des magnétomètres, où N est un nombre entier supérieur ou égal à cinq, - une unité (30) de traitement programmée pour : a) acquérir la mesure par chacun des magnétomètres de l'amplitude du champ magnétique le long de chacun de ses axes de mesure, 20 b) estimer à partir des mesures acquises des magnétomètres, des valeurs de plusieurs variables en résolvant un système actuel d'équations reliant ces variables à chaque mesure d'un magnétomètre tri-axes du réseau, ces variables correspondant aux positions, orientations et amplitudes des moments magnétiques des objets magnétiques des ustensiles à localiser, 25 - itérer les étapes a) et b) à des instants successifs dans le temps pour obtenir les valeurs estimées de chaque variable à ces différents instants successifs, caractérisé en ce que l'unité (30) de traitement est également programmée pour : c) établir des mesures des magnétomètres lorsqu'ils sont uniquement en présence d'un champ magnétique de référence généré seulement par les objets magnétiques 30 de tous les ustensiles immobiles présents devant le réseau de magnétomètres, et d) calculer l'amplitude d'une perturbation à partir des écarts entre les mesures réalisées lors de l'étape a) et les mesures établies lors de l'étape c), et - si l'amplitude de la perturbation franchit un seuil prédéterminé, exécuter l'étape b), sinon omettre l'exécution de l'étape b) et retourner à l'étape a) sans avoir exécuté 35 l'étape b).