FR 2943424 A1 20100924 FR 0951841 A 20090323 Procédé d'identification d'un axe de rotation unique d'un mouvement de rotation Domaine de l'invention La présente invention a pour objet un procédé d'identification d'un axe de rotation unique destiné à reconstruire un mouvement de rotation imposé par un utilisateur à un système d'application mobile. Plus particulièrement, la présente invention concerne un procédé de reconstitution d'un geste imposé à un téléphone mobile par un utilisateur. L'invention a également pour objet un téléphone mobile comportant des moyens pour mettre en oeuvre le procédé d'identification de l'axe de rotation relatif au mouvement de rotation imposé par l'utilisateur. Un but de l'invention est de déclencher une application à partir du mouvement de rotation reconstitué. Etat de la technique Aujourd'hui, les accéléromètres sont des composants électroniques de plus en plus courants dans des objets du quotidien tels que notamment les manettes de console de jeux vidéo où les différents mouvements appliqués à la manette se traduisent en action au niveau du jeu. Il est possible de trouver actuellement des téléphones mobiles également équipés d'accéléromètres. L'accéléromètre embarqué dans les téléphones mobiles permet en particulier de basculer automatiquement l'écran pour afficher un contenu dans le format le plus adapté. Aujourd'hui, des études pour inventer de nouveaux usages du téléphone mobile via l'accéléromètre se multiplient. Toutefois, ces études sont encore à l'état expérimental. En effet, par rapport aux manettes de console de jeux qui se trouvent dans un environnement stable, l'accéléromètre embarqué dans un téléphone mobile mesure un signal comportant, outre le mouvement imposé par l'utilisateur, des accélérations parasites. Ces accélérations parasites comportent principalement notamment l'accélération due à la gravitation terrestre et bien entendu les mouvements du corps de l'utilisateur et les vibrations de l'environnement où se situe l'utilisateur. Aujourd'hui, même s'il existe des procédés pour corriger l'accélération parasite due à la gravitation terrestre, le rapport signal/sur bruit reste relativement faible. Il est difficile d'identifier le mouvement relatif à l'utilisateur imposé à un téléphone mobile par rapport à d'autres mouvements parasites. De ce fait, l'accéléromètre embarqué dans les téléphones mobiles n'est utilisé que pour des mesures de déplacements par rapport à l'environnement, de pentes et d'inclinaisons, de détection de chocs ou des mesures de positionnement. Les mouvements associés à ces mesures sont relativement simples à identifier. Du fait que le nombre de mouvements faciles à identifier est limité, les utilisations des applications du téléphone mobile via l'accéléromètre sont également restreintes. Ainsi actuellement, il s'est fait sentir le besoin d'appliquer des mouvements complexes au téléphone mobile, tels que la rotation, afin d'utiliser les fonctions du téléphone sans avoir à manipuler des touches et à chercher dans des menus. Toutefois, les mouvements de rotation effectués par le bras de l'utilisateur ne sont pas parfaits du point de vue géométrique. En effet, ils sont liés à la chaîne squelettique. De ce fait, la reconstruction de ces mouvements de rotation est compliquée. Car d'une part le nombre de degrés de liberté de la chaîne squelettique participant au mouvement est généralement important et d'autre part le nombre de muscles impliqués dans le mouvement de chaque articulation est beaucoup plus grand que ce qui est nécessaire mécaniquement. Exposé de l'invention L'invention a justement pour but de répondre à ce besoin en remédiant aux inconvénients des techniques exposées précédemment. Pour cela, l'invention propose un procédé d'identification d'un axe de rotation fixe dans l'espace et unique afin de reconstruire le mouvement de rotation imposé par un utilisateur à un dispositif d'application mobile. Pour ce faire, l'invention met en oeuvre un algorithme d'estimation de l'axe de rotation le plus probable à partir des données fournies par l'accéléromètre embarqué. L'algorithme de l'invention permet d'améliorer les performances en termes d'estimation de l'axe de rotation, de coût calculatoire et de précision et robustesse des résultats. A cet effet, la présente invention vise un procédé d'estimation du mouvement d'un système mobile dans un milieu susceptible de générer des perturbations, le mouvement étant défini par un vecteur à trois variables et ledit système étant équipé d'au moins un capteur sensible à l'accélération ayant trois axes de sensibilité. Dans l'invention, au lieu de supprimer l'accélération parasite due à la gravitation terrestre, l'algorithme de l'invention l'utilise pour estimer l'axe de rotation. En effet, les accélérations dues au mouvement imposé par l'utilisateur sont largement inférieures à l'accélération due à la gravitation terrestre. De ce fait, la différence entre l'accélération due à la gravitation terrestre et l'accélération due au mouvement imposé est quasi constante tout au cours du mouvement et sa direction reste sensiblement inchangée. L'algorithme de l'invention effectue une projection de cette différence afin d'obtenir un nuage de points dans un plan. A partir de ce nuage de points, l'algorithme de l'invention calcule une estimation de la position de l'axe de rotation la plus probable. Puis une application associée à cet axe de rotation estimé est déclenchée. L'application peut consister entre autres en : - l'envoi d'un SMS, - l'activation du mode � mains libres �, - l'activation du renvoi d'appel, - la réponse à un appel entrant, - la réalisation d'un appel sortant, - l'ouverture d'une application telle que notamment un SMS, un MMS, une messagerie vocale ou Internet, etc.... - la navigation dans les menus du téléphone, - l'utilisation comme manette de jeu, - la mise en vibreur. Ce mode de déclenchement des applications d'un système mobile tel qu'un téléphone mobile est appelé à faire partie des fonctions standard dans les mobiles de prochaine génération. La présente invention présente ainsi l'avantage de proposer un algorithme apte à identifier des mouvements complexes du bras d'un utilisateur avec un taux de reconnaissance suffisant. L'invention permet ainsi d'augmenter les possibilités d'utilisation des accéléromètres embarqués dans un système mobile. La présente invention, peut par ailleurs être combinée avec des algorithmes existants destinés à estimer une position en translation du téléphone mobile par rapport à une position de départ. L'invention a donc pour objet un procédé de détermination d'un axe de rotation unique d'un mouvement de rotation imposé par un utilisateur à un système d'application mobile dans lequel : - on mesure des accélérations d'un accéléromètre à trois axes du 5 système pendant une durée donnée, - on détermine les orientations des accélérations mesurées, - on projette les accélérations mesurées à un module prédéfini dans le prolongement de leurs orientations obtenant un nuage de points, - on détermine un cercle médian en fonction du nuage de points, le 10 cercle médian formant un cône avec les orientations des accélérations mesurées, - on détermine une direction de l'axe de rotation unique correspondant à l'axe de révolution du cône, la direction de l'axe de rotation est définie par les coordonnées dudit axe dans un repère cartésien solidaire du système. 15 Avantageusement, l'invention est également caractérisée en ce que la détermination des orientations des accélérations mesurées comporte les étapes suivantes : - on exprime les accélérations mesurées selon des vecteurs d'accélération par rapport au repère cartésien solidaire au système, 20 - on détermine à chaque instant de mesure, une variation temporelle d'un vecteur d'accélération due au champ de gravitation local du système par rapport aux vecteurs d'accélération exprimés. Avantageusement, l'invention est également caractérisée en ce que la détermination du cercle médian comporte l'étape suivante : 25 - on détermine un rayon du cercle en fonction de l'application de l'algorithme des moindres carrés aux nuages de points, le centre du cercle passant par l'axe de révolution, - on détermine l'axe de révolution à partir du rayon du cercle médian et du centre du cercle. 30 Avantageusement, l'invention est également caractérisée en ce que - on détermine un plan comportant un maximum de points du nuage de points, - on détermine une droite perpendiculaire au plan correspondant à l'axe de rotation unique. 35 Avantageusement, l'invention est également caractérisée en ce que - on applique un filtre passe bas aux accélérations mesurées. Avantageusement, l'invention est également caractérisée en ce que le filtre passe bas est une moyenne glissante. Avantageusement, l'invention est également caractérisée en ce que - on affecte un attribut de confiance à l'axe de rotation unique déterminé en fonction de son évolution au cours du temps et/ou de la durée de la mesure. Avantageusement, l'invention est également caractérisée en ce que - on transmet à l'utilisateur une requête de confirmation du mouvement de rotation imposé, lorsque l'attribut de confiance est inférieur à un seuil de confiance prédéfini. Avantageusement, l'invention est également caractérisée en ce que - on associe dans une mémoire du système une direction donnée d'un axe de rotation à une application du système, - on déclenche l'application associée à cette direction de l'axe rotation déterminé, lorsque l'attribut de confiance est supérieur au seuil de confiance prédéfini. L'invention a également pour objet un téléphone mobile comportant des moyens aptes à mettre en oeuvre le procédé de l'invention de détermination d'un axe de rotation unique d'un mouvement de rotation imposée par un utilisateur audit téléphone mobile. . Brève description des dessins L'invention sera mieux comprise à la lecture de la description qui suit et à l'examen des figures qui l'accompagnent. Celles-ci sont présentées à titre indicatif et nullement limitatif de l'invention. La figure 1 montre une représentation schématique d'un système d'application mobile muni des moyens perfectionnés de l'invention. La figure 2 montre une représentation schématique en 3D du mouvement imposé au système mobile de la figure 1. La figure 3 montre une représentation schématique du résultat de projection des vecteurs accélérations à un module donné formant un nuage de points, selon l'invention. La figure 4 montre une représentation schématique d'un plan comportant le nuage de points. La figure 5 montre une illustration d'étapes correspondant à une identification d'un axe de rotation du mouvement imposé par un utilisateur au système mobile. La figure 6 montre un téléphone mobile muni des moyens perfectionnés de l'invention. Description détaillée de modes de réalisation de l'invention La figure 1 montre un système 10 d'application mobile. Ce système 10 permet de traduire en instructions un mouvement imposé audit système. Le système 10 peut être un téléphone mobile comme illustré à la figure 5, une manette de jeu, un terminal Internet mobile plus connu en anglais sous le nom de Mobile Internet Device (MID), une tablette PC, un livre électronique ou tout autre dispositif mobile embarquant un accéléromètre. Le système 10 comporte un capteur 11 de mouvement qui est un accéléromètre comportant trois axes de sensibilité. Les accélérations fournies par le capteur 11 sont de faible niveau. Elles sont largement inférieures à g, g étant une unité d'accélération correspondant approximativement à l'accélération de la pesanteur à la surface de la terre. La précision du capteur 11 requise est importante. Elle peut varier de 0,01 0/0 à 2 % de l'échelle de mesure du capteur. Le capteur 11 de mouvement peut être une jauge, un capteur capacitif, un capteur à induction, un capteur optique ou un potentiomètre. La liste n'est pas exhaustive. Le capteur 11 peut être couplé à un gyroscope pour obtenir des mesures complémentaires permettant de gagner en précision. Le capteur 11 de mouvement enregistre une accélération a d'un mouvement 12 du système 10 à des instants t; prédéfinis. Le mouvement 12 correspond à un mouvement de rotation, cette rotation étant représentée en pointillé à la figure 2. L'accélération a; est une grandeur vectorielle qui indique la modification affectant la vitesse du mouvement 12 en fonction du temps. Le système 10 est ici considéré comme ponctuel et est repéré par une position 1'4 à un instant t; comme indiqué sur la figure 2, Son orientation est repérée par le vecteur 9, qui correspond à l'orientation d'un repère R (x, y, z) dans un repère fixe Ro (Xo, Zo), Le repère fixe Ro correspond à un repère initial dans une phase de repos du système 10 d'applications au point 0, Dans un mode de réalisation préféré, le repère fixe Ro est un repère cartésien solidaire au système 10. Le repère fixe Ro est caractérisé par un triplet de vecteurs unitaires (i, j et k), qui représentent respectivement les directions des axes Ox, Oy et Oz. Le capteur 1 1 enregistre également à tout instant t; l'accélération de la pesanteur notée Go. Ici, à titre d'exemple non limitatif, l'axe Zo du repère fixe Ro est aligné avec la verticale et donc avec Co. Ainsi, lorsque e système 10 est dans une position initiale au repos, on peut définir un vecteur Co de dimension trois, correspondant à l'accélération de la pesanteur. Le mouvement 12 imposé au système 10 est défini par e vecteur a à trois variables correspondant au mouvement propre dudit système et aux mouvements parasites. Les trois variables correspondent à a position du système mobile 10 dans le repère R, correspondant ainsi aux coordonnées Y, Z dans le repère R. Lorsque le système 10 entre en mouvement, le vecteur a; est ajouté au vecteur initiai Co. Ainsi, le champ total mesuré par le dispositif de mesure 15 correspond à la somme des vecteurs Go + a. Le système 10 comporte un calculateur 13 représenté en pointillé à la figure 1. Le calculateur 13 est connecté au capteur 11. Le système 10 comporte également un dispositif de traitement 14 connecté au calculateur 13, Les calculateurs 13 et 14 peuvent éventuellement être mis en oeuvre 20 comme des processus exécutés sur la même entité physique. Dans la description, on prête des actions à des appareils ou à des programmes, cela signifie que ces actions sont exécutées par un microprocesseur de cet appareil ou de l'appareil comportant le programme, ledit microprocesseur étant alors commandé par des codes instructions 25 enregistrés dans une mémoire de l'appareil. Ces codes instructions permettent de mettre en oeuvre les moyens de l'appareil et donc de réaliser l'action entreprise. Le calculateur 13 comporte notamment un microprocesseur et une mémoire de programme (non représentés). La mémoire de programme est 30 divisée en plusieurs zones, chaque zone correspondant à une fonction ou à un mode de fonctionnement du programme du calculateur 13. Les zones de la mémoire de programme du calculateur 13 sont illustrées à la figure 1 par des rectangles 15 à 20. Une zone 15 comporte des codes instructions pour recevoir du 35 capteur 11 des coordonnées du vecteur accélération a. Du fait que les accélérations imposées par l'utilisateur au capteur 11 sont de faible niveau par rapport au vecteur Go, ces composantes peuvent être négligées lors des rotations typiques effectuées comme gestes avec le système 10. Par conséquent, il est considéré avec l'algorithme de l'invention que, au cours du geste, l'accélération à traiter AG;, correspondant à Go - a, conserve quasiment la direction verticale pointée vers le centre de la Terre et la valeur mesurées par rapport à un repère solidaire à la terre du vecteur accélération Go. Dans ce sens, la configuration adressée par cette invention diffère significativement des conditions imposées, par exemple, à une manette de jeu dont l'accélération imposée peut être très significative par rapport au vecteur Go. Les composantes du vecteur d'accélération AG sont perturbées par plusieurs types de bruit tel que notamment ceux susceptibles d'être générés par le milieu où se situe le système 10 (voiture, train ...) et les mouvements involontaires de l'utilisateur. La zone 15 comporte de ce fait des codes instructions pour filtrer les composantes reçues du capteur 11. Un filtre passe-bas est appliqué aux données reçues afin de supprimer les composantes à haute fréquence. Dans un mode de réalisation simple à mettre en oeuvre, bien qu'avantageux et efficace, le filtre passe-bas appliqué aux données reçues est une moyenne glissante. Bien entendu, d'autres types de filtres passe-bas peuvent être utilisés. La moyenne glissante, ou moyenne mobile, est un type de moyenne statistique permettant de supprimer les fluctuations transitoires des données reçues, de façon à en souligner les tendances à plus long terme. Cette moyenne est dite mobile parce qu'elle est recalculée de façon continue, en utilisant à chaque calcul un sous-ensemble d'éléments dans lequel un nouvel élément remplace le plus ancien ou s'ajoute au sous-ensemble. En sortie de la zone 15 de filtrage, le vecteur G;,m est un macro-vecteur dont chacune de ses composantes est remplacée par une macro- composante (Gx;,,, Gy;rr:, Gz;m). Une macro-composante est une moyenne d'un sous-ensemble de N composantes consécutives. Dans un exemple, N égale cinq. Une zone 16 comporte des codes instructions pour calculer une puissance d'accélération à partir des données reçues de la zone 15. La puissance d'accélération correspond à une différence entre les modules de deux macro-vecteurs G,-, consécutifs. Une zone 17 comporte des codes instructions pour comparer la puissance d'accélération fournie par la zone 16 à un premier seuil S1 de cinétique prédéfini. Dès que la puissance d'accélération est supérieure au seuil S1, alors des codes instructions d'une zone 18 et d'une zone 19 sont exécutés. Les codes instructions de la zone 18 traitent les données fournies par la zone de filtrage 15 afin de déterminer l'axe de rotation du mouvement imposé au système 10. Les codes instructions de la zone 19 traitent les données fournies par la zone de filtrage 15 afin de déterminer un attribut de confiance à attribuer à l'axe de rotation déterminé. Cet attribut de confiance peut être exprimé en pourcentage. Les figures 3 et 4 montrent des exemples de déterminations de l'axe de rotation par les codes instructions de la zone 18. Une zone 20 comporte des codes instructions pour comparer la puissance d'accélération fournie par la zone 16 à un deuxième seuil S2 de cinétique prédéfini. Dès que la puissance d'accélération est inférieure au seuil S2, alors l'exécution des codes instructions de la zone 18 et de la zone 19 est interrompue. Ainsi, le début du geste et la fin du geste de l'utilisateur sont détectés par rapport à l'augmentation ou la diminution de la dynamique. Le calculateur 13 fournit en entrée au dispositif 14 de traitement les données de sortie des zones 18 et 19. Le dispositif 14 de traitement comporte notamment un microprocesseur 21, une mémoire 22 de programme et une mémoire 23 de données connectés par un bus 24. La mémoire 22 de programme est divisée en plusieurs zones, chaque zone correspondant à une fonction ou à un mode de fonctionnement du programme du dispositif 14 de traitement. Une zone 25 comporte des codes instructions pour recevoir du calculateur 13 des données correspondant à un axe de rotation et au pourcentage attribué à cet axe. Une zone 26 comporte des codes instructions pour comparer le pourcentage attribué à un seuil S3 de confiance. Lorsque le pourcentage attribué est supérieur ou éventuellement égal au seuil S3 alors des codes instructions d'une zone 27 sont exécutés. Les codes instructions de la zone 27 déclenchent une application via la mémoire 23 de données. Sinon, des codes instructions d'une zone 28 sont exécutés. Les codes instructions de la zone 28 transmettent une requête de confirmation du mouvement imposé à l'utilisateur. Cette requête peut être un SMS ou un message vocal. La mémoire 23 de données est, par exemple, structurée en une table. Par exemple, chaque ligne de la table correspond à une position de l'axe de rotation, chaque colonne de la table correspond à un renseignement sur cet axe de rotation. Ainsi, la mémoire 23 de données comporte une ligne 29a correspondant à une position donnée d'un axe de rotation et une colonne 29b correspondant à une adresse AD1 d'une application de services du système 10 à déclencher. L'adresse AD est une adresse mémoire et peut être un simple identifiant alphanumérique Les figures 3 et 4 montrent deux modes de détermination de l'axe de rotation unique du mouvement 12 imposé par l'utilisateur au système 10, via les codes instructions de la zone 18. La zone 18 du calculateur 13 projette les vecteurs d'accélérations Gim délivrés par la zone 15 de filtrage. Les macro-vecteurs d'accélérations Gim sont projetés dans le sens de leurs orientations. Les vecteurs d'accélérations Gim projetés forment un nuage de points Ti. Puisque le module du vecteur accélération Gim filtré est quasiment constant alors le nuage de points représentés dans le repère local se trouvera dans des circonférences 32;. Soit P le vecteur de module unitaire qui représente l'axe de rotation du système au cours du geste. Sa position peut être représentée par ses trois composantes : P = PX.i + Py.j + PZ.k. Chaque type de geste de l'utilisateur est associé à une certaine position de P. Les coordonnées (Px, P,, PZ) sont par conséquent la signature du mouvement. Les codes instructions de la zone 18 sont aptes à identifier ces coordonnées. Et les codes instructions de la zone 19 sont aptes à définir un degré de confiance mesurable à attribuer aux coordonnées identifiés. Puisque le vecteur accélération Gim est quasiment toujours vertical et le geste considéré est une rotation autour d'un axe unique P fixe, l'angle formé par G;, et l'axe de rotation P est considéré comme constant au cours de la rotation. Par conséquent, la ligne droite parallèle au vecteur G;rr: passant par le centre du repère fixe décrira une surface conique 34 de révolution. Le problème résolu par l'invention est de déterminer la direction de l'axe 35 du cône correspondant à l'axe P de rotation unique. Un premier mode de détermination de l'axe 35 du cône 34 consiste à estimer la position dans l'espace d'un cercle 32; médian le plus probable par rapport au nuage de points Ti. Il est important de souligner que, inévitablement, il y a toujours des dérives par rapport à la rotation idéale. Ceci dépend, en grande mesure, des déplacements secondaires de la main au cours de la rotation et des mouvements de l'utilisateur par rapport au milieu dans lequel il se situe et des perturbations de ce milieu. Afin de tenir compte de ces écarts par rapport au mouvement idéal, l'algorithme applique des méthodes de minimisation d'erreurs permettant d'établir la meilleure estimation de P en fonction des mesures réalisées. Pour déterminer le cercle 32; le plus probable, le calculateur 13 calcule pour chaque point Ti le rayon r; le plus probable. Dans un mode de réalisation préféré, le rayon r; attribué à chaque point Ti du nuage de points projetés est déterminé selon l'algorithme des moindres carrés. D'autres modes de calculs existants peuvent être utilisés pour déterminer le rayon r du cercle 32 médian le plus probable. L'axe 35 passant par le centre du cercle 32 médian le plus probable correspond à l'axe P de rotation unique du mouvement 12. La figure 4 montre un autre mode de détermination préféré de l'axe 35 du cône 34. Puisque le module du vecteur accélération Gin, est quasiment constant alors le nuage de points forme des circonférences 32; situées dans un plan 36 perpendiculaire à l'axe P de rotation. Pour déterminer l'axe P de rotation unique, il faudra d'abord identifier ce plan 36 et en déduire le vecteur normal. L'équation du plan 36 est définie comme suit : x = a + by + cz. Les codes instructions de la zone 18 du calculateur 13 calculent les coefficients a, b et c du plan 36. Dans un mode de réalisation préféré, ces coefficients sont déterminés selon l'algorithme des moindres carrés appliqué aux fonctions linéaires à deux variables : f(y, z) = a + by + cz. Les mesures du vecteur accélération sont appelées : Gim =Gx;m.i + Gyim.j +Gzim.k où i allant de 1 à n et n représentant le nombre total de mesures réalisées au cours de la rotation. Les équations ci-après montrent un exemple de calcul des coefficients du plan 36. L'erreur de l'ensemble des mesures de i allant de 1 à n élevé au carré est : E =E [Gxim.i ù (a + b. Gyim.j +c. Gzim.k)]2 Les conditions nécessaires pour minimiser l'erreur sont : aE/aa = 0, aE/ab = 0 et aE/ac = 0, ce qui se traduit par les relations matricielles suivantes : m11a+m12b+m13c=q1 m21a + m22b + m23c = q2 m31a + m32b +m33C = q3 où les coefficients mil sont donnés par les expressions suivantes : m11=n m12=EGyim m13=EGzim m21 = m12 m22 = > Gy2im m23 = > Gyim Gzim m31 = m13 m32 = m23 m33 = E Gz2im q1 = I Gxim q2 = > Gxim Gyim q3 = > Gxim Gzim Le système d'équations linéaires précédent peut ainsi se représenter comme un rapport entre matrices : m11 m12 m13 m21 m22 m23 m31 m32 m33 M= ql q2 q3 a b X= Q= c M=XQ Sous réserve que M ≠ 0, ce système permet de récupérer les coefficients a, b et c. Le plan déterminé par le calculateur 13 correspond au plan le plus probable comportant le maximum de points possible du nuage de points. La direction de l'axe de rotation est donnée par normalisation en utilisant les expressions suivantes : 1 ùb ùc py V1+b2 +c2 V1+b2 +c2 +b2+c2 px pz = La figure 5 montre une illustration d'étapes correspondant à une identification, par le système 10, d'un axe de rotation du mouvement 12 imposé par un utilisateur. La figure 5 montre une première étape préliminaire 100 dans laquelle le système 10 est en mode veille. Lorsque le système 10 est dans cet état de veille, le calculateur 13 est activé en permanence. Le système 10 ne sortira de cet état de veille que lorsque la puissance de l'accélération fournie par le calculateur 13 sera supérieure au premier seuil S1. Dans une variante, le calculateur 13 est activé manuellement lors d'une étape 101. Dans ce cas, les mesures du capteur 11 ne sont traitées que si l'utilisateur démarre le fonctionnement du calculateur 13. Ce démarrage peut être effectué en cliquant sur un bouton ou en appuyant sur une touche du système 10. Même si l'implémentation de l'étape 100 est préférée à celle de l'étape 101, du fait qu'elle adhère à une notion d'interface utilisateur naturelle, la consommation d'énergie liée à l'implémentation de l'étape 100 peut obliger à adopter l'implémentation de l'étape 101. A une étape 102, le système 10 est activé. Le calculateur 13 applique, à une étape 103, un filtre passe bas aux mesures d'accélérations fournies par le capteur 13 de mouvement. Ce filtrage permet de supprimer, des composantes du vecteur accélération, des perturbations dues par plusieurs types de bruit. Ce filtrage permet de filtrer les composantes à haute fréquence. A une étape 104, le calculateur 13 calcule une puissance d'accélération pour chaque vecteur d'accélération G;m, filtré. La puissance d'accélération à un instant t; correspond à la différence entre un vecteur d'accélération à l'instant t;_1 et un vecteur d'accélération à l'instant t;. A une étape 105, dès qu'une puissance d'accélération est supérieure au seuil S1 alors, le vecteur d'accélération G4, correspondant est enregistré dans une mémoire de données du calculateur à une étape 106. Si la puissance d'accélération est comprise entre le seuil S1 et le seuil S2, à une étape 107, alors le vecteur d'accélération G;rn correspondant est enregistré dans la mémoire de données du calculateur à l'étape 106. Le calculateur 13 affecte un attribut de début de rotation au premier vecteur d'accélération G;rr: d'une série de données enregistrées dans la mémoire de données du calculateur 13. Un attribut est un paramètre interne du calculateur 13. Il peut être un nombre binaire ou alphanumérique. Le calculateur 13 affecte un attribut de fin de rotation au dernier vecteur d'accélération d'une série de données enregistrées dans la mémoire de données du calculateur 13. Le calculateur 13, selon le mode de réalisation préféré de l'invention, surveille ainsi en permanence les signaux en provenance de l'accéléromètre 11 pour déclencher l'évaluation de la rotation quand un niveau cinétique suffisamment élevé est identifié. Les vecteurs d'accélération G;,,, filtrés et enregistrés dans la mémoire de données à l'étape 106 représentent les données permettant de reconstituer le mouvement imposé. A une étape 108, le calculateur 13 projette les vecteurs d'accélération Ge enregistrés à une longueur de module prédéfinie afin d'obtenir un nuage de points Ti. A une étape 109, le calculateur 13 détermine soit le rayon r le plus probable du cercle formé par le nuage de points où les cordonnées du plan 36 comportant le maximum de points du nuage de points. Cette détermination est par exemple basée sur l'algorithme des moindres carrés. L'étape 109 permet de fournir en sortie à une étape 110, un axe de rotation estimé. Puisque l'utilisateur peut imposer au système 10 un mouvement qui n'est pas une rotation autour d'un axe fixe, l'estimation de l'algorithme doit être validée. A une étape 111, le calculateur 13 détermine un attribut de confiance à affecter à l'axe de rotation estimé. Si le mouvement de rotation imposé est une rotation faible, l'axe de rotation devrait rester fixe dans l'espace pour toute la durée de la rotation. Pour déterminer l'attribut de confiance, le calculateur détermine jusqu'à quel point l'axe de rotation du système 10 ne change pas au cours du temps. Dans un exemple préféré, le critère adopté consiste à calculer une moyenne de l'écart des extrémités des vecteurs d'accélération par rapport aux projections sur le plan estimé ou le cercle estimé. Le résultat de cette opération correspond à un niveau d'erreur à attribuer à l'axe de rotation estimé. Ce niveau d'erreur varie de 0 à 1. Plus ce niveau d'erreur est proche de zéro, plus la confiance à accorder à l'axe de rotation estimée est importante. A une étape 112, le calculateur 13 transmet au dispositif 14 de traitement l'axe de rotation estimé et le niveau de confiance attribué. A une étape 113, le dispositif 14 de traitement compare le niveau de confiance attribué à un seuil S3 de confiance prédéfini. Si le niveau de confiance est supérieur au seuil S3, alors l'estimation de l'axe s'avère suffisant. Cet axe de rotation estimé est validé à une étape 114 comme l'axe P de rotation unique du mouvement 12 de rotation. A une étape 115, le dispositif 14 de traitement exécute l'application associée dans la mémoire de données à l'axe de rotation déterminé. Si, par contre, le niveau de confiance attribué reste inférieur au seuil S3, alors l'estimation de l'axe de rotation est invalidée et l'utilisateur est prié de répéter le geste, à une étape 116. Dans ce cas, une requête de confirmation du mouvement 12 est transmise par le dispositif 14 de traitement à l'interface de l'utilisateur. Le seuil S3 de confiance peut être par exemple égal à 90%. En général le niveau de confiance attribué est inférieur au seuil S3, lorsque la durée du mouvement 12 est faible, par exemple inférieure à 1 seconde. Dans ce cas, l'accéléromètre n'a pas récupéré assez de données pour permettre au dispositif de traitement de prendre une décision. De même, le niveau de confiance attribué est inférieur au seuil S3, lorsque le milieu environnant est particulièrement vibrant. En général pour une bonne estimation de l'axe de rotation P, la durée du mouvement 12 doit être sensiblement supérieure à 3 secondes. La figure 6 montre un téléphone mobile150 d'un utilisateur, client d'un opérateur, en liaison avec un réseau 151 d'un opérateur de téléphonie mobile. Le téléphone mobile 150 est muni du système 10 de la figure 1. Le terminal 150 comporte, outre une carte SIM (Subscriber Identity Module) contenant l'IMSI (International Mobile Subscriber Identity) et le MSISDN (Mobile Subscriber Integrated Services Digital Number), un microprocesseur 152, une première mémoire 153 de programme, une mémoire 154 de données, des circuits 155 d'interface radio avec le réseau 150 de l'opérateur de téléphonie, une deuxième mémoire 156 de programme et un capteur de mouvement 157 à trois axes. Les éléments 152 à 157 sont connectés par un bus 158. La mémoire 153 de programme correspond à la mémoire de programme du dispositif de traitement et la mémoire 156 correspond à la mémoire de programme du calculateur 13. La représentation des mémoires de données et des mémoires de programmes dans le téléphone mobile n'est qu'une illustration d'implantation de composants et d'enregistrements de données. Dans la pratique, ces mémoires sont unifiées ou distribuées selon des contraintes de taille (capacité) et/ou de rapidité des traitements souhaités. En outre, la mémoire programme du dispositif de traitement peut être répartie sur d'autres serveurs de l'opérateur de téléphonie mobile, éventuellement selon le type d'applications à déclencher. Dans un mode de réalisation préféré, les mémoires de données et les mémoires de programmes sont enregistrées dans la carte SIM ou dans la 15 mémoire interne du téléphone mobile. REVENDICATIONS1 ù Procédé de détermination d'un axe de rotation (P) unique d'un mouvement (12) de rotation imposé par un utilisateur à un système d'application (10) mobile dans lequel : - on mesure des accélérations (Go ù a;) d'un accéléromètre (11) à trois axes du système pendant une durée donnée, - on détermine les orientations (8) des accélérations mesurées, - on projette (108) les accélérations mesurées à un module prédéfini dans le prolongement de leurs orientations obtenant un nuage de points (Ti), - on détermine un cercle médian (32;) en fonction du nuage de points, le cercle médian formant un cône (34) avec les orientations des accélérations mesurées, - on détermine une direction de l'axe de rotation unique correspondant à l'axe de révolution du cône, la direction de l'axe de rotation est définie par les coordonnées dudit axe dans un repère cartésien (Ro) solidaire du système. 2 ù Procédé selon la revendication 1, caractérisé en ce que la détermination des orientations des accélérations mesurées comporte les étapes suivantes : - on exprime les accélérations mesurées selon des vecteurs d'accélération par rapport au repère cartésien solidaire au système, - on détermine à chaque instant de mesure, une variation temporelle (AG;) d'un vecteur d'accélération due au champ de gravitation local du système par rapport aux vecteurs d'accélération exprimés. 3 ù Procédé selon l'une des revendications 1 à 2, caractérisé en ce que la détermination du cercle médian comporte l'étape suivante : - on détermine un rayon (r;) du cercle en fonction de l'application de l'algorithme des moindres carrés aux nuages de points, le centre du cercle passant par l'axe de révolution, - on détermine l'axe de révolution à partir du rayon du cercle médian et du centre du cercle. 4 ù Procédé selon l'une des revendications 1 à 2, caractérisé en ce que : - on détermine un plan (36) comportant un maximum de points dunuage de points, - on détermine une droite perpendiculaire au plan correspondant à l'axe de rotation unique. ù Procédé selon l'une des revendications 1 à 4, caractérisé en ce 5 que : - on applique (15, 103) un filtre passe bas aux accélérations mesurées. 6 ù Procédé selon la revendication 5, caractérisé en ce que le filtre passe bas est une moyenne glissante. 7 ù Procédé selon l'une des revendications 1 à 6, caractérisé en ce que : - on affecte (19, 111) un attribut de confiance à l'axe de rotation unique déterminé en fonction de son évolution au cours du temps et/ou de la durée de la mesure. 8 ù Procédé selon la revendication 7, caractérisé en ce que : - on transmet (28, 116) à l'utilisateur une requête de confirmation du mouvement de rotation imposé, lorsque l'attribut de confiance est inférieur à un seuil (S3) de confiance prédéfini. 9 ù Procédé selon l'une des revendications 7 à 8, caractérisé en ce que : - on associe dans une mémoire (23) du système une direction donnée d'un axe de rotation (29a) à une application (29b) du système, - on déclenche l'application associée à cette direction de l'axe rotation déterminé, lorsque l'attribut de confiance est supérieur au seuil de confiance prédéfini. 10 ù Téléphone mobile (150) caractérisé en ce qu'il comporte des moyens aptes à mettre en oeuvre un procédé de détermination d'un axe de rotation (P) unique d'un mouvement de rotation imposée par un utilisateur audit téléphone mobile selon l'une quelconque des revendications précédentes.