FR 2959112 A1 20111028 FR 1053092 A 20100423 SYSTEME D'ANALYSE DE FOULEES D'UN UTILISATEUR L'invention porte sur un système d'analyse de foulées d'un utilisateur. De manière générale, l'invention est liée aux domaines de la capture de mouvement, et de la géolocalisation d'individus en déplacement. La capture de mouvement concerne des applications grand public telles des applications de loisir (consoles de jeu interactives, suivi du geste sportif, réalité virtuelle ou réalité augmentée), des applications d'aide à la navigation 1 o de piétons (les plus utilisés à l'heure actuelle étant les systèmes de navigation par satellites tel le GPS), des applications d'aide à la mobilité de personnes vulnérables ou temporairement fragilisées par leur environnement (personnes handicapées ou plongées dans l'obscurité), et des applications de fitness (podomètre, calcul de la dépense énergétique ou de la distance 15 parcourue). La capture de mouvement concerne également des applications médicales (suivi des personnes âgées et/ou dépendantes, analyse de la marche pour la rééducation posturale ou aide au diagnostic), des applications de sécurité ou de secours (localisation de pompiers à l'intérieur d'un bâtiment en feu, suivi opérationnel de militaires, ou surveillance de 20 prisonniers), mais aussi des applications à visée commerciale (statistiques sur les trajectoires suivies par les consommateurs dans les centres commerciaux ou les supermarchés, définition d'archétypes d'usage, ou proposition de services commerciaux topo-dépendants). Il est notamment connu de reconstituer le déplacement d'un objet 25 ou d'une personne munie d'un objet comprenant un émetteur de signaux reconnus par un système d'aide à la navigation par satellites, tel le système GPS, en extérieur, et par un système de radiolocalisation à base de transmissions Ultra Large Bande (ULB) ou WiFi en intérieur. Dans de nombreuses zones géographiques, externes ou internes, 30 la navigation à base de système d'aide à la navigation par satellites ou de radiolocalisation s'avère néanmoins très délicate, en raison de l'obstruction des signaux radio nécessaires à la mesure de métriques topo-dépendantes, par exemple en cas d'indisponibilité des signaux émis par un ou plusieurs satellites du système d'aide à la navigation, nécessaires à la mesure de 35 pseudo-distances. En pratique, ces situations de blocage ou d'obstruction dégradent fortement la précision de l'information de localisation, et rendent même parfois le service d'aide à la navigation indisponible, comme illustré dans le document "An evaluation of UWB localization under non line-of-sight (NLOS) propagation" (A. Maali, A. Ouldali, H. Mimoun, et G. Baudoin), publié dans Wireless Pervasive Computing, ISWPC, 3rd International Symposium on, pages 379û382, en mai 2008. Généralement, les réalisations citées précédemment présentent un coût, une consommation, un encombrement, et/ou une infrastructure inadaptés, voire rédhibitoires pour des applications grand public. 1 o D'autres systèmes, ont pour but d'élaborer un trajet à partir d'un point de départ donné, en utilisant des capteurs d'attitude (accéléromètre(s), magnétomètre(s), gyromètre) délivrant le plus souvent des mesures inertielles. Ces mesures permettent d'assurer la navigation par des moyens plus ou moins précis et compliqués, tels que décrits, par exemple dans les 15 documents "Inertial head-tracker sensor fusion by a complementary separate-bias kalman filter" (E. Foxlin) publié en mars 1996 dans Virtual Reality Annual International Symposium, Proceedings of the IEEE 1996, pages 185û194; "Detection of spatio-temporal gait parameters by using wearable motion sensors" (K. Kogure L. Seon-Woo, et K. Mase), publié en 20 2005 dans Engineering in Medicine and Biology Society. IEEE-EMBS 2005. 27th Annual International Conference of the, pages 6836û6839, 2005; "Pedestrian tracking with shoe-mounted inertial sensors" (E. Foxlin), publié en novembre 2005 dans Computer Graphics and Applications, IEEE, 25 :38û46, Nov-Dec 2005; ou "Integration of foot-mounted inertial sensors into a 25 Bayesian location estimation framework" (P. Robertson B. Krach), publié en mars 2008 dans Positioning, Navigation and Communication, 2008. WPNC 2008. 5th Workshop on, pages 55û61. De tels systèmes sont de coût et de complexité élevés, et manquent souvent de précision. Des consoles de jeu vidéo, telle la Wii, utilisent des capteurs 30 optiques et/ou à ultrasons pour déterminer la trajectoire d'un élément de commande de jeu. Ces systèmes sont couteux et limités. Il existe également, comme présentés dans les documents "Assessment of walking features from foot inertial sensing" (S. Scapellato F. Cavallo A.M. Sabatini, et C. Martelloni), publié en mars 2005 dans 35 Biomedical Engineering, IEEE Transactions on, pages 486û 494, ou "Multisensor approach to walking distance estimation with foot inertial sensing" (D. Alvarez A.M Lopez J. Rodriguez-Uria J.C. Alvarez, et R.C. Gonzalez), publié en août 2007 dans Engineering in Medicine and Biology Society, 2007, EMBS 2007. 29th Annual International Conference of the IEEE, pages 5719-5722, des systèmes dans lesquels on suppose deux accéléromètres et un gyromètre disposés dans le plan sagittal d'un utilisateur, qui permettent de déterminer l'attitude de ce capteur dans le plan. L'accélération est intégrée dans l'axe du pied. Ces systèmes ont un coût réduit, mais présentent d'autres inconvénients. D'une part, ils supposent que les capteurs sont parfaitement placés dans le plan sagittal, ce qui est quasiment impossible à réaliser concrètement par l'utilisateur, et provoque une erreur d'estimation liée au mauvais positionnement du ou des capteurs. D'autre part, ils supposent que la marche a lieu dans le plan sagittal, ce qui peut ne pas être le cas, lors d'une marche sur le côté, par exemple. Le document, "An innovative shoe-mounted pedestrian navigation system" (K. Fyfe, Gérard Lachapelle, R. Stirling, et J. Collin), publié en avril 2003 dans Proc. European Navigation Conf. (GNSS), CD-ROM, Austrian lnst. of Navigation, divulgue un système muni de trois accéléromètres et d'un quatrième accéléromètre afin de calculer l'angle subit par le capteur. Cet accéléromètre supplémentaire est placé sur le même élément mobile, mais à distance de ce dernier. Les deux capteurs voient donc la même rotation et le même déplacement. Une telle structure permet d'estimer la vitesse de rotation de l'élément mobile, en rotation autour d'une direction sensiblement constante, en se passant d'un gyromètre. Ce système a un coût réduit du au remplacement d'un gyromètre mono-axe par un couple d'accéléromètres par rapport aux systèmes avec des gyromètres proposés par Alvarez et al et Scapellato et al. Cependant, un tel système présente les mêmes inconvénients que l'exemple cité ci-dessus, car seule la rotation autour de l'axe orthogonal au plan sagittal est prise en compte. En outre, l'estimation de la rotation par un couple d'accéléromètres, tout comme celle réalisée à l'aide d'un gyromètre, présente des dérives importantes. Un cap est également déterminé à l'aide d'un magnétomètre, par exemple lorsque le pied est posé, ce qui donne la direction du pied et non pas la direction du déplacement recherchée. La demande de brevet français FR 0857181 porte sur un dispositif de détermination d'une caractéristique trajectoire formée de positions successives d'un accéléromètre triaxial lié de manière solidaire à un élément mobile, entre un premier instant d'immobilité ou de quasi-immobilité et un deuxième instant d'immobilité de quasi immobilité de l'accéléromètre triaxial, consécutif audit premier instant, ledit dispositif comprenant, en outre, un capteur additionnel triaxial de mesure d'un vecteur d'un champ vectoriel sensiblement constant entre lesdits premier et deuxième instants d'immobilité, dans un repère global fixe lié au repère terrestre. Le capteur 1 o additionnel, par exemple un magnétomètre triaxial, est lié de manière solidaire audit élément mobile et fixe dans le repère de l'accéléromètre, et le dispositif comprend des moyens de commande. Lesdits moyens de commande comprennent : - des premiers moyens de détermination desdits premier 15 et deuxième instants d'immobilité de l'accéléromètre triaxial ; - des deuxièmes moyens de détermination d'un axe de rotation sensiblement invariant de l'accéléromètre triaxial entre lesdits premier et deuxième instants d'immobilité, et d'un plan orthogonal audit axe de rotation sensiblement invariant, dans un repère mobile lié à 20 l'accéléromètre ou au capteur additionnel à partir des vecteurs délivrés par l'accéléromètre triaxial, ou à partir des vecteurs délivrés par le capteur additionnel ; - des premiers moyens de calcul, auxdits instants successifs, de premières projections orthogonales sur ledit plan des vecteurs 25 délivrés par l'accéléromètre triaxial et des vecteurs délivré par le capteur additionnel dans ledit repère mobile ; - des troisièmes moyens de détermination, auxdits instants successifs, de la rotation de passage dudit repère mobile audit repère global fixe, à partir desdites premières projections orthogonales des 30 vecteurs délivrés par le capteur additionnel ; - des deuxièmes moyens de calcul, auxdits instants successifs, de deuxièmes projections orthogonales dans ledit repère global fixe, des premières projections orthogonales dans ledit plan des vecteurs délivrés par l'accéléromètre triaxial, fournies par lesdits premiers moyens de 35 calcul ; - des troisièmes moyens de calcul pour soustraire à chaque deuxième projection orthogonale dans ledit repère global fixe le vecteur moyen sur lesdits instants successifs, de manière à obtenir les accélérations centrées dans ledit plan, dépourvues de l'influence de la gravité terrestre et de dérives dudit dispositif, dans ledit repère global fixe ; et - des quatrièmes moyens de calcul d'une caractéristique de la trajectoire à partir des accélérations centrées. Un tel dispositif permet de déterminer une trajectoire planaire d'un accéléromètre triaxial lié à un élément mobile, notamment un pied, pour une 1 o foulée, pour ainsi reconstruire une trajectoire d'un ensemble de foulées successives. L'élément mobile peut être une autre partie du corps humain, comme une main, une partie d'un corps animal, ou une partie d'un corps artificiel tel un robot ou une souris informatique. Un tel dispositif permet d'obtenir des informations globales liées 15 au déplacement de l'utilisateur dans un plan sagittal. Il ne permet pas d'obtenir des informations, en trois dimensions, relatives à la foulée d'un utilisateur, notamment d'un coureur. Un but de l'invention est de pallier ces problèmes. Selon un aspect de l'invention, il est proposé un système 20 d'analyse de foulées d'un utilisateur, comprenant : - un premier ensemble capteur muni d'un boîtier, d'un magnétomètre triaxial lié fixement audit boîtier, et d'un référentiel mobile ; - des premiers moyens de fixation adaptés pour fixer le premier ensemble capteur à un segment de jambe de l'utilisateur ; et 25 - des moyens de traitement adaptés pour calculer l'angle de lacet et/ou de tangage et/ou de roulis dudit ensemble capteur, et pour calculer la ou les variations temporelles angulaires correspondantes, dans un référentiel global fixe lié au référentiel terrestre, à partir des mesures dudit magnétomètre triaxial délivrées dans ledit référentiel mobile. 30 Un tel système permet d'obtenir des informations, en trois dimensions, relatives à la foulée d'un utilisateur. En outre, les moyens de traitement peuvent être internes ou externes au boîtier, et lorsqu'ils sont externes, ils peuvent être embarqués à bord d'un téléphone portable ou d'un ordinateur portable. On entend par jambe, un membre inférieur complet. De manière préférentielle, le segment de jambe auquel peut être fixé le premier ensemble capteur EC1, est le cou-de-pied ou le métatarse du pied, ou le tibia. Le premier ensemble capteur peut être calibré, si nécessaire, pour que le référentiel mobile soit tel que deux axes forment un plan vertical dans un référentiel global fixe lié au référentiel terrestre. Dans la suite de la description, le terme référentiel peut être associé à un repère associé au référentiel. 1 o Selon un mode de réalisation, lesdits moyens de traitement comprennent : - des premiers moyens de calcul de la variation temporelle de l'angle de tangage du premier ensemble capteur à partir des mesures dudit magnétomètre triaxial délivrées dans ledit référentiel mobile ; et 15 - des premiers moyens de détermination d'un instant de début d'une foulée et d'un instant de fin de ladite foulée et de phases de ladite foulée, à partir de ladite variation temporelle de l'angle de tangage du premier ensemble capteur fournie par lesdits premiers moyens de calcul. Il est ainsi possible, par exemple, de déterminer avec précision le 20 temps de contact au sol (entre la pose du pied au sol et l'instant où le pied va quitter le sol) ou l'accélération du pied à l'instant instant où le pied va quitter le sol, qui sont d'importantes caractéristiques de la foulée d'un coureur. Dans un mode de réalisation, ledit premier ensemble capteur est, 25 en outre, muni d'un accéléromètre triaxial lié fixement audit boîtier, et lesdits moyens de traitement comprennent des deuxièmes moyens de détermination de valeurs de paramètres représentatifs d'une phase de ladite foulée ou de la foulée complète, en trois dimensions, dans le référentiel global fixe. 30 Ainsi la précision du système est améliorée, en obtenant des informations liées à des phases de la foulée. Selon un mode de réalisation, une foulée comprenant une phase d'appui du pied auquel est fixé le premier ensemble capteur, suivi d'une phase de cycle arrière dans laquelle ledit pied remonte du sol en direction 35 du dos de l'utilisateur, et d'une phase de cycle avant dans laquelle ledit pied se dirige devant l'utilisateur, lesdits deuxièmes moyens de détermination comprenant des premiers moyens d'estimation d'un référentiel lié au premier ensemble capteur, pour la foulée et/ou par phase de la foulée, à partir d'une décomposition en valeurs singulières, de la matrice à trois colonnes formée par des vecteurs délivrés par ledit magnétomètre triaxial à des instants successifs, dans le référentiel mobile lié au premier ensemble capteur. L'invention permet de déterminer un référentiel variable lié au premier ensemble capteur, pour l'ensemble de la foulée et par phase de la foulée, les phases d'une foulée sont définies plus loin dans la description. 1 o Les deux vecteurs associés aux deux plus grandes valeurs singulières forment le plan principal de la foulée. Dans un mode de réalisation, lesdits moyens de traitement comprennent des deuxièmes moyens de calcul de projections orthogonales des vecteurs de mesures délivrés par l'accéléromètre triaxial et par le 15 magnétomètre triaxial sur les axes du référentiel lié au premier ensemble capteur sur la foulée. Ainsi, il est possible d'obtenir l'évolution des axes du magnétomètre vu d'un référentiel réputé fixe durant la foulée, et d'en déduire l'évolution des angles décrivant l'attitude spatiale de l'ensemble capteur, 20 donc du segment de jambe auquel il est fixé, par exemple le pied. Selon un mode de réalisation, le référentiel global fixe lié au référentiel terrestre est ladite projection orthogonale des vecteurs délivrés par le magnétomètre triaxial dans le référentiel lié au premier ensemble capteur pour la foulée. 25 En effet, cela est possible car sur une foulée, le champ magnétique reste sensiblement constant dans un référentiel global fixe lié au référentiel terrestre. Dans un mode de réalisation, lesdits deuxièmes moyens de détermination comprennent des troisièmes moyens de détermination des 30 angles de tangage, de roulis et de lacet du premier ensemble capteur, à partir de ladite projection orthogonale des vecteurs délivrés par le magnétomètre triaxial dans le référentiel lié au premier ensemble capteur pour la foulée , et de détermination des vitesses angulaires correspondantes de tangage, de roulis et de lacet. Par défaut on détermine ainsi les six degrés de liberté, puis ensuite on teste les quels peuvent être négligés. Selon un mode de réalisation, lesdits deuxièmes moyens de détermination comprennent des moyens de comparaison desdites vitesses 5 angulaires de tangage, de roulis et de lacet avec un seuil. On peut ainsi tester quels degrés de liberté peuvent être négligés. Dans un mode de réalisation, lesdits deuxièmes moyens de détermination comprennent des moyens d'utilisation du référentiel lié au premier ensemble capteur pour la phase en cours, lorsqu'au moins deux 1 o desdites valeurs absolues des trois vitesses angulaires sont supérieures audit seuil, pour redresser l'orientation par compensation des déviations de la trajectoire, dues aux rotations, par rapport au référentiel global en utilisant ledit référentiel lié audit premier ensemble capteur pour la phase en cours à la place du référentiel lié au premier ensemble capteur pour la foulée. 15 Le redressement d'orientation s'entend par déterminer la trajectoire dans un référentiel adéquat, différent du référentiel calculé pour l'ensemble de la foulée, à cause des vitesses angulaires non négligeables qui font dévier la trajectoire du référentiel de la foulée complète. Ainsi, la trajectoire reconstruite est composée de trajectoires 20 unitaires, pour chaque phase de la foulée, dont les référentiels respectifs sont adaptés au nombre de degrés de liberté à prendre en compte pour chacune des phases. Selon un mode de réalisation, lesdits premiers moyens de détermination sont adaptés pour identifier des instants caractéristiques de la 25 foulée de l'utilisateur par détection d'extrema locaux de la variation temporelle de l'angle de tangage du premier ensemble capteur. Il est ainsi aisé de déterminer ces instants caractéristiques de la foulée. Dans un mode de réalisation, les deuxièmes moyens de 30 détermination comprennent des troisièmes moyens de calcul pour soustraire, au vecteur correspondant à la projection orthogonale dans le référentiel global du vecteur de mesures de l'accéléromètre triaxial délivré par les deuxièmes moyens de calcul, le vecteur moyen sur la foulée de la projection orthogonale dans le référentiel global du vecteur de mesures de 35 l'accéléromètre triaxial délivré par les deuxièmes moyens de calcul, de manière à obtenir le vecteur centré des accélérations dans le référentiel global. L'obtention des accélérations centrées dans le référentiel global, permet de s'abstraire de l'influence de la gravité. Selon un mode de réalisation, les deuxièmes moyens de détermination comprennent des quatrièmes moyens de calcul d'une double intégration temporelle dudit vecteur centré des accélérations dans le référentiel global délivré par lesdits troisièmes moyens de calcul, de manière à obtenir la trajectoire du premier ensemble capteur dans ledit référentiel 1 o global, sur la foulée. Ainsi, on obtient une caractéristique en trois dimensions de la foulée de l'utilisateur, pour chaque phase de la foulée. En outre, les deuxièmes moyens de détermination comprennent des moyens de correction des dérives, en supposant que le pied de 15 l'utilisateur a une orientation identique en début et fin de foulée, par redressement, de manière à avoir des coordonnées verticales et transversales identiques en début et fin de foulée. Ainsi, la trajectoire est redressée, sans avoir recours à des calculs d'angles entre le référentiel global et le référentiel terrestre. 20 Le système peut également comprendre un deuxième ensemble capteur semblable au premier ensemble capteur, et des deuxièmes moyens de fixation adaptés pour fixer le deuxième ensemble capteur à un segment de l'autre jambe de l'utilisateur. Le système peut alors, en outre, comprendre des moyens de 25 synchronisation pour comparer les mesures effectuées par le ou lesdits ensembles capteurs d'une jambe par rapport à l'autre. Ainsi, il est d'une part possible de caractériser la fin de cycle avant d'un pied en mesurant la contribution de l'accélération à l'instant où l'autre pied va quitter le sol. 30 D'autre part, la mesure synchrone des deux trajectoires des ensembles capteurs liés aux différentes jambes de l'utilisateur permet de déterminer la présence d'asymétries temporelles et spatiales de la foulée du coureur. En outre, le système peut comprendre un troisième ensemble 35 capteur similaire au premier ensemble capteur, et des troisièmes moyens de fixation adaptés pour fixer le troisième ensemble capteur à mi-hauteur de l'utilisateur. Ainsi, une mesure des trajectoires des deux segments de jambe, notamment des deux pieds, effectuée de manière synchrone avec la mesure de l'oscillation du centre de masse du coureur dans le plan sagittal permet de déterminer la phase de course à modifier de manière à diminuer l'amplitude de ladite oscillation. Selon un autre aspect de l'invention, il est également proposé un procédé d'analyse de foulées d'un utilisateur, dans lequel : on calcule la variation temporaire de l'angle de tangage d'un premier ensemble capteur, comprenant un magnétomètre triaxial, fixé à un segment de jambe de l'utilisateur, à partir des mesures dudit magnétomètre triaxial ; - on détermine un instant de début d'une foulée et un instant de fin de ladite foulée et des phases de ladite foulée, à partir de calcul de la variation temporelle de l'angle de tangage du premier ensemble capteur ; et - on détermine des valeurs de paramètres représentatifs d'une phase de ladite foulée, en trois dimensions, dans un référentiel global fixe lié au 20 référentiel terrestre. L'invention sera mieux comprise à l'étude de quelques modes de réalisation décrits à titre d'exemples nullement limitatifs et illustrés par les dessins annexés sur lesquels : 25 - les figures 1 a, 1 b, et 1c illustrent schématiquement la décomposition d'une foulée, selon un aspect de l'invention ; la figure 2 illustre schématiquement le déplacement d'un ensemble capteur fixé à un pied d'un utilisateur, selon un aspect de l'invention ; 30 les figures 3a et 3b illustrent schématiquement un système selon un aspect de l'invention ; la figure 4 illustre plus en détails un module de traitement selon un aspect de l'invention ; et les figures 5a, 5b, 5c, 5d et 5e illustrent des éléments 35 du module de traitement de la figure 4. Dans l'ensemble de figures, les éléments ayants les mêmes références sont similaires. Dans la présente description, on décrit un système complexe, notamment illustré par la figure 4, mais l'invention couvre également une réalisation moins complexe comprenant uniquement les éléments permettant de calculer l'angle de lacet et/ou de tangage et/ou de roulis, ainsi que la ou les variations temporelles angulaires correspondantes, dans un référentiel global fixe lié au référentiel terrestre, à partir des mesures d'un magnétomètre triaxial délivrées dans un référentiel lié au 1 o magnétomètre. En effet, toute la partie du système illustré sur la figure 4, dont les flèches entre éléments sont représentées en pointillés, est optionnelle, pour obtenir le système complexe, le reste de la figure représente un mode de réalisation moins complexe permettant de déterminer l'angle de lacet (YawGF) et/ou de tangage (PitchGF) et/ou de roulis (RoIIGF) 15 dudit premier ensemble capteur (EC1), et pour calculer la ou les variations d YawGF d PitchGF d Roll" temporelles angulaires correspondantes ( , , ), dans dt dt dt un référentiel global fixe lié au référentiel terrestre, à partir des mesures dudit magnétomètre triaxial délivrées dans ledit référentiel mobile. Tel qu'illustré sur la figure 1 a, une foulée d'un coureur est 20 décomposée en deux pas successifs (pose du pied gauche et du pied droit). Une foulée est une phase de course qui se déroule entre deux contacts au sol successifs d'un même pied. Sur la figure 1 b, sont représentées les phases d'une foulée. La foulée se décompose en deux phases distinctes. La phase d'appui pendant 25 laquelle un des deux pieds est en contact avec le sol, et La phase de suspension pendant laquelle les deux pieds sont en l'air. Pendant une phase de vol d'un pied (suspension et appui de l'autre pied) se divise elle-même en deux phases : une phase de cycle arrière durant laquelle le pied remonte du sol en direction du dos du coureur, et une phase de cycle avant durant 30 laquelle le pied se dirige dans le sens de la course, ou devant l'utilisateur. La figure 1c représente les événements ou instants caractéristiques d'une foulée considérés. Au cours de la foulée, la trajectoire du pied est caractérisée par quatre instants caractéristiques notables : - maintien du pied à plat ou FF pour "Foot Flat" en langue anglaise correspondant au maintien du bassin au dessus du pied ; - instant où le pied va quitter le sol ou TO pour "Toe Off" en langue anglaise ; - instant auquel le pied repart vers l'avant ou MSW pour "Mid Swing" durant 5 la phase de suspension ; et - pose du pied au sol ou HS pour "Heel Strike" en langue anglaise correspondant au début de la phase d'amortissement qui s'achève lorsque le bassin passe au dessus du pied. La figure 2 illustre le plan principal Pf et le référentiel Rf de la 1 o foulée. Au cours de la foulée, la trajectoire du pied se déroule au voisinage d'un plan Pf étant le plan principal de la foulée. Ce plan présente un angle en général non nul vis-à-vis du plan sagittal qui est défini par la direction de la foulée au sol et la verticale. Le référentiel de la foulée est construit à partir de deux axes du plan principal Pf de la foulée et d'un axe normal. 15 La figure 3a illustre un système d'analyse de foulées d'un utilisateur, comprenant un premier ensemble capteur EC1 muni d'un boîtier BT, d'un magnétomètre triaxial 3M lié fixement au boîtier BT. L'ensemble capteur EC1 est muni d'un référentiel mobile LF, tel que deux axes (y, z) forment un plan vertical dans un référentiel global fixe GF lié au référentiel 20 terrestre. Le système comprend également des premiers moyens de fixation MF1, par exemple des éléments élastiques munis de scratchs, adaptés pour fixer le premier ensemble capteur EC1 à un segment d'une jambe de l'utilisateur, par exemple un pied de l'utilisateur. En outre, le premier ensemble capteur EC1 comprend un accéléromètre triaxial 3A optionnel lié 25 fixement au boîtier BT. En effet, dans le système décrit, outre l'extraction des angles et vitesses angulaires, on détermine la trajectoire en trois dimensions du premier ensemble capteur EC1, et donc du segment de jambe auquel il est lié fixement. Le système comprend, en outre, un module de traitement MT. 30 En l'espèce, le module de traitement MT est interne au boîtier BT, et peut transmettre des résultats à un terminal extérieur, tel un ordinateur ou un téléphone portable, pour les afficher. En variante, un élément d'affichage des résultats pourrait être directement intégré au boîtier BT. En variante, tel qu'illustré sur la figure 3b, le module de traitement MT peut être déporté dans un terminal extérieur, tel un ordinateur ou un téléphone portable. La figure 4 illustre plus en détails un mode de réalisation du 5 module de traitement MT. Un premier module de détermination DET1 permet de déterminer un instant de début d'une foulée, un instant de fin de ladite foulée, et des phases de ladite foulée, à partir d'une variation temporelle de l'angle de tangage Pitch du premier ensemble capteur EC1. Un deuxième module de 10 détermination DET2 permet de déterminer des valeurs de paramètres représentatifs d'une phase de ladite foulée, ou de la foulée complète, en trois dimensions, dans le référentiel global fixe GF. Le module de traitement MT comprend un premier module de calcul CALC1 de la variation temporelle de l'angle de tangage Pitch du 15 premier ensemble capteur EC1 dans le référentiel mobile LF à partir des mesures du magnétomètre triaxial 3M. Le deuxième module de détermination DET2 comprend un premier module d'estimation EST1 d'un référentiel lié au premier ensemble capteur EC1, pour la foulée Rf et par phase Rap, Rar, Rav de la foulée, à partir 20 d'une décomposition en valeurs singulières ou SVD pour "singular value decomposition" en langue anglaise, dans ledit référentiel mobile LF lié au premier ensemble capteur EC1, de la matrice à trois colonnes formée par des vecteurs délivrés par le magnétomètre triaxial 3M à des instants successifs dans le référentiel mobile LF. 25 Le module de traitement MT comprend un deuxième module de calcul CALC2 de projections orthogonales PO des vecteurs de mesures délivrés par l'accéléromètre triaxial 3A et par le magnétomètre triaxial 3M sur les axes du référentiel lié au premier ensemble capteur sur la foulée Rf. Le référentiel global fixe GF, pour "global frame" en langue 30 anglaise, lié au référentiel terrestre, est la projection orthogonale PO des vecteurs délivrés par le magnétomètre triaxial 3M dans le référentiel lié au premier ensemble capteur EC1 pour la foulée Rf. Le deuxième module de détermination DET2 comprend un troisième module de détermination DET3 des angles de tangage PitchRf , de 35 roulis RollRf et de lacet YawRf du premier ensemble capteur EC1, dans le référentiel Rf, à partir de la projection orthogonale PO des vecteurs délivrés par le magnétomètre triaxial 3M dans le référentiel LF lié au premier ensemble capteur EC1 pour la foulée , et de détermination des vitesses d PitchR d RollR angulaires correspondantes de tangage , de roulis et de dt dt d YawRf lacet dt Le deuxième module de détermination DET2 comprennent des moyens de comparaison COMP des vitesses angulaires de tangage d PitchR d RollR d YawR dt , de roulis dt f et de lacet dt f avec un seuil S. Le deuxième module de détermination DET2 comprend un module d'utilisation UTIL du référentiel Rap, Rar, Rav lié au premier ensemble capteur EC1 pour la phase en cours, lorsqu'au moins deux desdites valeurs absolues des trois vitesses angulaires sont supérieures audit seuil S, pour redresser l'orientation, c'est-à-dire pour replacer la trajectoire durant la phase en cours de la foulée, en utilisant le référentiel Rap, Rar, ou Rav lié à l'ensemble capteur EC1 pour la phase en cours à la place du référentiel Rf lié au premier ensemble capteur EC1 pour l'ensemble de la foulée. Le premier module de détermination DET1 est adapté pour identifier des instants caractéristiques HS, FF, TO, MSW de la foulée de l'utilisateur par détection d'extrema locaux de la variation temporelle d Pitch dt de l'angle de tangage Pitch dans le référentiel lié au premier ensemble capteur EC1. En outre, la séparation en phases par le premier ensemble capteur EC1 permet de distinguer les moments de la foulée pendant lesquels la trajectoire en trois dimensions peut être théoriquement estimée, des moments où le nombre de degrés de liberté trop important l'en empêche (ces nombres sont obtenus par extraction des vitesses angulaires). La phase d'appui va de la pose du pied au sol HS à l'instant où le pied va quitter le sol TO : le nombre de degrés de liberté majoritaires est deux, car le roulis Roll et le tangage Pitch, les accélérations et le lacet Yaw y sont négligeables. La phase de cycle arrière va de l'instant ou le pied va quitter le sol TO à l'instant séparant le cycle arrière du cycle avant MSW : le nombre de degrés de liberté est généralement de six, car les trois accélérations et les trois angles ne sont pas négligeables. La phase de cycle avant va de l'instant séparant le cycle arrière du cycle avant MSW à la pose du pied au sol HS : le roulis Roll et le tangage Yaw y sont négligeables, donc il est quatre degrés de liberté majoritaires. L'avantage de la connaissance de ce nombre de degrés de liberté est de pouvoir valider ou non la trajectoire reconstruite obtenue. Dans les phases d'appui et de cycle avant, le fait d'avoir moins de cinq degrés de liberté permet de reconstruire la trajectoire en trois dimensions. Dans la 1 o phase de cycle arrière, le nombre de degrés de liberté étant généralement de six, une estimation plus approfondie des paramètres doit être mise en oeuvre afin d'avoir la meilleure estimation possible de la trajectoire en trois dimensions au cours de cette phase. Le découpage de la foulée permet de traiter séparément les 15 signaux, on peut calculer un référentiel de la foulée Rf différent pour chacune des phases. Il y aura ainsi trois décompositions en valeurs singulières SVD effectuées au lieu d'une, et les projections sont faites phase par phase. En effet, dans la phase d'appui au sol, par exemple, le référentiel de la foulée Rf peut être très différent du référentiel global GF estimé sur une foulée entière. 20 Le fait de changer de référentiel permet de donner plus de précisions dans l'orientation du mouvement. Le deuxième module de détermination DET2 comprend, un troisième module de calcul CALC3 pour soustraire, au vecteur correspondant à la projection orthogonale PO dans le référentiel global GF du vecteur de 25 mesures de l'accéléromètre triaxial 3A délivré par le deuxième module de calcul CALC2, le vecteur moyen sur la foulée de la projection orthogonale PO dans le référentiel global GF du vecteur de mesures de l'accéléromètre triaxial 3A délivré par le deuxième module de calcul CALC2, de manière à obtenir le vecteur centré des accélérations dans le référentiel global GF. 30 Ce centrage permet d'obtenir des accélérations centrées, i.e. dépourvues de l'influence de la gravité terrestre dans ledit repère global fixe. Le deuxième module de détermination DET2 comprend des quatrièmes moyens de calcul CALC4 d'une double intégration temporelle dudit vecteur centré des accélérations dans le référentiel global GF délivré 35 par le troisième module de calcul CALC3, de manière à obtenir la trajectoire du premier ensemble capteur EC1 dans ledit référentiel global GF, sur la foulée. En outre, un module de correction CORR des dérives en supposant que le pied a une orientation identique en début et fin de foulée. Ce module de correction effectue une déviation de la trajectoire de manière à ce que les coordonnées en hauteur ou altitude dans le référentiel global fixe GF soient les mêmes au début et à la fin d'une foulée, et les coordonnées perpendiculaires ou dans le plan du sol également. Le système peut, en outre, comprendre un deuxième ensemble 1 o capteur EC2 semblable au premier ensemble capteur EC1, et un deuxième module de fixation MF2 adaptés pour fixer le deuxième ensemble capteur EC2 à l'autre pied de l'utilisateur. Un module de synchronisation SYNC permet d'observer les évènements provenant des ensembles EC1 et EC2 avec une base de temps commune. On peut ainsi, à tout instant, connaître 15 l'écart temporel entre un événement mesuré par le premier ensemble capteur EC1 (l'instant où le pied gauche va quitter le sol, par exemple) et un autre événement mesuré par le deuxième ensemble capteur EC2 (la pose du pied droit au sol, par exemple). En outre, le système peut comprendre un troisième ensemble 20 capteur EC3 similaire au premier ensemble capteur EC1, et un troisième module de fixation MF3 adapté pour fixer le troisième ensemble capteur EC3 à mi-hauteur de l'utilisateur, par exemple au niveau de l'abdomen. La figure 5a représente un pied muni d'un premier ensemble capteur EC1. La vitesse angulaire, ou variation temporelle de l'angle de 25 tangage Pitch du premier ensemble capteur EC1, à l'instant t, est donnée i par (Pitch(tn+1) ù Pitch(tn )) , avec Pitch(t) = arctan mA/ , où mAP et mvT m17 sont les mesures du magnétomètre 3M enregistrées dans le plan principal Pf de la foulée (directions y et z). La figure 5b représente la détection, par le premier module de 30 détection DET1, des instants caractéristiques HS, FF, TO, MSW de la foulée. En début de pas FF, l'accélération est supposée nulle. Les différents instants caractéristiques cités précédemment HS, FF, TO, et MSW de la foulée sont détectés aux extrema locaux. Ce type de détection est proche de c la méthode décrite dans la demande de brevet FR 2926971. La décomposition en valeurs singulières SVD permet de déterminer les axes principaux du mouvement. Les étapes de calcul sont : - une étape d'autocorrélation des trois composantes du magnétomètre triaxial 3M pour obtenir la matrice de covariance, à trois lignes et trois colonnes, des signaux du magnétomètre 3M ; - une étape de décomposition en éléments propres de la matrice de covariance, chaque valeur propre mesurant la part de variance selon le vecteur propre associé ; et - une étape d'analyse en composantes principales, dans laquelle les 1 o données sont projetées sur chacun des axes définis par les vecteurs propres. On détermine le plan [171,172] de la trajectoire à partir des vecteurs délivrés par le magnétomètre triaxial 3M. Le plan [171,172] correspond au plan contenant au maximum les mesures, considérées comme des points ou 15 vecteurs dans un espace à trois dimensions en considérant les mesure transmises pour chaque axe du magnétomètre triaxial, et l'axe porté par le vecteur 173 est l'axe sur lequel la projection des mesures varie le moins. Le vecteur 173 est le vecteur normal au plan de la trajectoire mais aussi l'axe de rotation du premier ensemble capteur EC1. 20 Pour déterminer ces trois vecteurs 171,172,173, on forme une matrice ALF à n+1 lignes, n+1 étant le nombre d'échantillons pris entre to et tn compris, et à 3 colonnes correspondant aux trois mesures des trois axes x, y, z du premier ensemble capteur EC1. On a : Ax LF (t0) A(t0) Az LF \t0 ALF (t0 ,...,tn )_ A:F (tn) AyLF (tn) kir (tn 25 Les mesures permettant de former la matrice A' peuvent être les mesures du magnétomètre triaxial 3M ou de l'accéléromètre triaxial 3A. En, variante, on peut ne prendre qu'un sous-ensemble des échantillons aux instants successifs to, t1, ...tn. Les trois vecteurs û1,à2,û3,sont déterminés à partir de la matrice 30 ALF(=ALF(t0,...,tn)), par décomposition en valeurs singulières SVD de la matrice ALF : ALF =USV. V est une matrice à trois lignes et trois colonnes dont chaque colonne correspond à un des trois vecteurs 171,172,173. La matrice S contient les trois valeurs singulières si, s2, s3. Ces trois valeurs singulières si, s2, s3, doivent êtres classées par ordre décroissant. Le vecteur ii est la colonne de la matrice V correspondant à la plus petite valeur singulière. Les deux autres colonnes donnent les coordonnées des vecteurs formant le plan Pf principal de la foulée [û1,û2]. Le troisième module de calcul CALC3 soustrait la gravité après avoir projeté l'accélération mesurée dans le référentiel fixe GF, pour extraire l'accélération propre. A ces fins, on centre les données. En effet, on suppose avoir une vitesse nulle à chaque pied à plat FF, ie au début et à la fin de chaque 1 o foulée. Or v(t)=f a(u) du + v(0) donc si T est la durée de la foulée, 0 v(T) = f a(u) du + v(0) = 0 . Puisque v(0) = v(T) = 0, et en discrétisant la somme o N intégrale, on obtient que : a(tt) = 0 â = 0. i=o La moyenne globale de l'accélération est donc nulle sur une 15 foulée, et il suffit donc de soustraire sa moyenne à chacune des composantes projetées sur les 3 axes du référentiel fixe, la contribution de la gravité y étant constante. Le quatrième module de calcul CALC4 effectue une double intégration temporelle du vecteur centré des accélérations dans le référentiel 20 global GF délivré par le troisième module de calcul CALC3, de manière à obtenir la trajectoire du premier ensemble capteur EC1 dans le référentiel global GF, sur la foulée. La méthode d'intégration utilisée est une méthode de Simpson adaptée ou méthodes des rectangles. La première intégration permet de 25 passer de l'accélération à la vitesse échantillon par échantillon. La deuxième intégration donne la distance à partir de la vitesse, dans les 3 dimensions. Or cette vitesse est supposée nulle à l'instant t=0 (FF), ce qui n'est jamais le cas en pratique. Il faut rajouter une vitesse tangentielle d PitchR R. , R étant la distance entre le capteur et le métatarse du pied (où dt 30 la vitesse est réellement nulle au posé au sol) et d PitchR d PitchR dt = PitchRf (tn+1) ù PitchRf (tn )), la vitesse de tangage dt f étant celle fournie par le troisième module de détermination DET3. Des données expérimentales ont permis de faire une estimation de R ainsi que de l'orientation de cette vitesse dans le référentiel du premier ensemble capteur, afin de donner des contributions suffisamment précises à chacun des trois axes d'intégration. Le module de correction CORR des dérives, après intégration, trois trajectoires sont obtenues de façon indépendante, selon les trois axes du référentiel fixe. Soit x l'axe horizontal dans la direction de la foulée, Z l'axe 1 o vertical et y l'axe horizontal perpendiculaire au plan de foulée. On fait l'hypothèse que le pied a la même orientation au début et à la fin de la foulée, et donc on redresse d'une part la trajectoire z = f(x), d'autre part y = f(x), en établissant y(T) = y(0) et z(T) = z(0). Par exemple, comme illustré sur la figure 5c, sont représentées 15 les courbes z = f(x), et y = f(x) : cette dernière, projetée dans le plan de la figure, est redressée pour que le dernier point soit commun aux deux trajectoires. En d'autres termes, il s'agit de la trajectoire reconstruite dans un plan perpendiculaire au plan donné par les deux directions principales de la trajectoire avec ici une représentation simple en deux dimensions, 20 qui est redressée, i.e. dont les coordonnées du posé de pied à la fin de la foulée sont supposées être les mêmes que celles au début de la foulée, verticalement et transversalement (le pied doit retoucher le sol et on suppose que le coureur court selon un segment de droite pour faciliter la représentation finale de la trajectoire. 25 Le troisième module de détermination DET3 des angles de tangage PitchRf , de roulis RollRf et de lacet YawRf du premier ensemble capteur EC1, à partir de ladite projection orthogonale PO des vecteurs délivrés par le magnétomètre triaxial 3M dans le référentiel Rf lié au premier ensemble capteur EC1 pour la foulée, et de détermination des vitesses d PitchR d RollR 30 angulaires correspondantes de tangage f , de roulis f et de dt dt d YawR lacet f dt Les mesures du magnétomètre triaxial 3M, projetées sur le référentiel de la foulée Rf, sont utilisées pour calculer les trois angles du mouvement : - le tangage PitchRf : rotation autour de l'axe perpendiculaire au mouvement obtenu par arctan (Bu2/Bu1) ; - le roulis RollRf : rotation autour de l'axe du sens du mouvement obtenu par arctan (Bu1/BU3) ; et - le lacet YawRf : rotation autour de l'axe vertical obtenu par arctan (Bu3/Bu2) ; 1 o avec BU1 / BU2 / Bu3 les composantes du champ magnétique sur les trois axes du référentiel de la foulée Rf : mUl / mU2 / mU3. En variante, ces angles peuvent être obtenus différemment, par exemple en obtenant la matrice de rotation mettant en jeu les trois angles, à partir des projections But, Bu2 et Bu3. 15 Connaissant ces trois angles, on peut obtenir les vitesses d PitchR d RollR d YawR angulaires f , f et f de manière instantanée. dt dt dt Ces valeurs sont généralement faibles si les angles ne varient pas de façon significative. Par exemple, au cours de la phase de cycle avant, seule la variation du tangage Pitch n'est pas négligeable. Le nombre de 20 degrés de liberté dans le mouvement y est donc inférieur ou égal à quatre. La trajectoire reconstruite peut être validée lors de cette phase. Par contre, si deux ou trois de ces valeurs sont supérieures en valeur absolue au seuil S, et que les accélérations selon les trois dimensions sont de plus importantes, on a plus de quatre degrés de liberté, et la validité 25 de la trajectoire est réduite. Dans ce cas, on peut corriger la trajectoire reconstruite en affinant l'orientation de la foulée grâce à une SVD effectuée dans cette phase, qui nous donnera le plan du mouvement, différent du plan obtenu sur l'ensemble de la foulée. 30 La figure 5d représente les valeurs des trois angles de tangage PitchRf , de roulis RollRf et de lacet YawRf exprimés en radians, en fonction du temps, exprimé en secondes, au cours d'une foulée, ainsi que la trajectoire verticale de cette foulée afin de visualiser les évolutions de ces angles en fonction des phases de la foulée. 5 Un tel système permet d'obtenir des informations, en trois dimensions, relatives à la foulée d'un utilisateur, notamment d'un coureur. 10 REVENDICATIONS1. Système d'analyse de foulées d'un utilisateur, comprenant : - un premier ensemble capteur (EC1) muni d'un boîtier, d'un magnétomètre triaxial (3M) lié fixement audit boîtier (BT), et d'un référentiel mobile (LF) ; et - des premiers moyens de fixation (MF1) adaptés pour fixer le premier 1 o ensemble capteur (EC1) à un segment de jambe de l'utilisateur ; caractérisé en ce qu'il comprend des moyens de traitement (MT) adaptés pour calculer l'angle de lacet (YawGF) et/ou de tangage (PitchGF) et/ou de roulis (RoIIGF) dudit premier ensemble capteur (EC1), et pour calculer la ou les variations temporelles angulaires correspondantes d YawGF d PitchGF d RollGF 15 ( , , ), dans un référentiel global fixe (GF) lié au dt dt dt référentiel terrestre, à partir des mesures dudit magnétomètre triaxial (3M) délivrées dans ledit référentiel mobile (LF). 2. Système selon la revendication 1, dans lequel lesdits 20 moyens de traitement (MT) comprennent : - des premiers moyens de calcul (CALC1) de la variation temporelle de l'angle de tangage (Pitch) du premier ensemble capteur (EC1) à partir des mesures dudit magnétomètre triaxial (3M) délivrées dans ledit référentiel mobile (LF); et 25 - des premiers moyens de détermination (DET1) d'un instant de début d'une foulée et d'un instant de fin de ladite foulée et de phases de ladite foulée, à partir de ladite variation temporelle de l'angle de tangage (Pitch) du premier ensemble capteur (EC1) fournie par lesdits premiers moyens de calcul. 30 3. Système selon la revendication 2, dans lequel : - ledit premier ensemble capteur (EC1) est, en outre, muni d'un accéléromètre triaxial (3A) lié fixement audit boîtier (BT) ; etlesdits moyens de traitement (MT) comprennent des deuxièmes moyens de détermination (DET2) de valeurs de paramètres représentatifs d'une phase de ladite foulée ou de la foulée complète, en trois dimensions, dans le référentiel global fixe (GF). 4. Système selon la revendication 3, dans lequel, une foulée comprenant une phase d'appui du pied auquel est fixé le premier ensemble capteur (EC1), suivi d'une phase de cycle arrière dans laquelle ledit pied remonte du sol en direction du dos de l'utilisateur, et d'une phase 1 o de cycle avant dans laquelle ledit pied se dirige devant l'utilisateur, lesdits deuxièmes moyens de détermination (DET2) comprenant des premiers moyens d'estimation (EST1) d'un référentiel lié au premier ensemble capteur (EC1), pour la foulée (Rf) et/ou par phase (Rap, Rar, Rav) de la foulée, à partir d'une décomposition en valeurs singulières (SVD), de la matrice à trois 15 colonnes formée par des vecteurs délivrés par ledit magnétomètre triaxial (3M) à des instants successifs, dans le référentiel mobile (LF) lié au premier ensemble capteur (EC1). 5. Système selon la revendication 3 ou 4, dans lequel 20 lesdits moyens de traitement (MT) comprennent des deuxièmes moyens de calcul (CALC2) de projections orthogonales (PO) des vecteurs de mesures délivrés par l'accéléromètre triaxial (3A) et par le magnétomètre triaxial (3M) sur les axes du référentiel lié au premier ensemble capteur sur la foulée (Rf). 25 6. Système selon la revendication 5, dans lequel le référentiel global fixe (GF) lié au référentiel terrestre est ladite projection orthogonale (PO) des vecteurs délivrés par le magnétomètre triaxial (3M) dans le référentiel lié au premier ensemble capteur (EC1) pour la foulée (Rf). 30 7. Système selon la revendication 5 ou 6, dans lequel lesdits deuxièmes moyens de détermination (DET2) comprennent des troisièmes moyens de détermination (DET3) des angles de tangage (PitchRf ), de roulis (RollRf ) et de lacet (YawRf) du premier ensemble capteur (EC1), à partir de ladite projection orthogonale (PO) des vecteurs délivrés 35 par le magnétomètre triaxial (3M) dans le référentiel (Rf) lié au premierensemble capteur (EC1) pour la foulée , et de détermination des vitesses d Pitch / d Roll angulaires correspondantes de tangage Rf , de roulis Rf dt dt d YawR f de lacet dt 8. Système selon la revendication 7, dans lequel lesdits deuxièmes moyens de détermination (DET2) comprennent des moyens de d PitchR f comparaison (COMP) desdites vitesses angulaires de tangage dt /d Roll /dYawR de roulis Rf et de lacet f avec un seuil (S). dt dt 10 9. Système selon la revendication 8, dans lequel, lesdits deuxièmes moyens de détermination (DET2) comprennent des moyens d'utilisation (UTIL) du référentiel (Rap, Rar, Rav) lié au premier ensemble capteur (EC1) pour la phase en cours, lorsqu'au moins deux desdites valeurs absolues des trois vitesses angulaires sont supérieures audit seuil (S), pour 15 redresser l'orientation par compensation des déviations de la trajectoire, dues aux rotations, par rapport au référentiel global (GF) en utilisant ledit référentiel (Rap, Rar, Rav) lié audit premier ensemble capteur (EC1) pour la phase en cours à la place du référentiel (Rf) lié au premier ensemble capteur (EC1) pour la foulée. 20 10. Système selon l'une des revendications 6 à 9, dans lequel lesdits premiers moyens de détermination (DET1) sont adaptés pour identifier des instants caractéristiques (HS, FF, TO, MSW) de la foulée de l'utilisateur par détection d'extrema locaux de la variation temporelle 25 /d pitch de l'angle de tangage (Pitch) du premier ensemble capteur (EC1). dt J 11. Système selon l'une des revendications 5 à 10, dans lequel les deuxièmes moyens de détermination (DET2) comprennent des troisièmes moyens de calcul (CALC3) pour soustraire, au vecteur5correspondant à la projection orthogonale (PO) dans le référentiel global (GF) du vecteur de mesures de l'accéléromètre triaxial (3A) délivré par les deuxièmes moyens de calcul (CALC2), le vecteur moyen sur la foulée de la projection orthogonale (PO) dans le référentiel global (GF) du vecteur de mesures de l'accéléromètre triaxial (3A) délivré par les deuxièmes moyens de calcul (CALC2), de manière à obtenir le vecteur centré des accélérations dans le référentiel global (GF). 12. Système selon l'une des revendications 3 à 11, dans lequel les deuxièmes moyens de détermination (DET2) comprennent des quatrièmes moyens de calcul (CALC4) d'une double intégration temporelle dudit vecteur centré des accélérations dans le référentiel global (GF) délivré par lesdits troisièmes moyens de calcul (CALC3), de manière à obtenir la trajectoire du premier ensemble capteur (EC1) dans ledit référentiel global (GF), sur la foulée. 13. Système selon l'une des revendications 3 à 12, dans lequel les deuxièmes moyens de détermination (DET2) comprennent des moyens de correction (CORR) des dérives, en supposant que le pied de l'utilisateur a une orientation identique en début et fin de foulée, par redressement, de manière à avoir des coordonnées verticales et transversales identiques en début et fin de foulée. 14. Système selon l'une des revendications précédentes, comprenant un deuxième ensemble capteur (EC2) semblable au premier ensemble capteur (EC1), et des deuxièmes moyens de fixation (MF2) adaptés pour fixer le deuxième ensemble capteur (EC2) à un segment de l'autre jambe de l'utilisateur. 15. Système selon la revendication 12, comprenant, en outre des moyens de synchronisation (SYNC) pour comparer les mesures effectuées par le ou lesdits ensembles capteurs d'une jambe par rapport à l'autre. 16. Système selon la revendication 13 ou 14, comprenant, en outre un troisième ensemble capteur (EC3) similaire au premier ensemble capteur, et des troisièmes moyens de fixation (MF3) adaptés pour fixer le troisième ensemble capteur (EC3) à mi-hauteur de l'utilisateur. 17. Procédé d'analyse de foulées d'un utilisateur, dans lequel : - on calcule la variation temporaire de l'angle de tangage d'un premier ensemble capteur (EC1), comprenant un magnétomètre triaxial (3M), fixé à 1 o un segment de jambe de l'utilisateur, à partir des mesures dudit magnétomètre triaxial (3M) ; caractérisé en ce que : - on détermine un instant de début d'une foulée et un instant de fin de ladite foulée et des phases de ladite foulée, à partir de calcul de la variation 15 temporelle /d pitch de l'angle de tangage (Pitch) du premier ensemble dt capteur (EC1) ; et - on détermine des valeurs de paramètres représentatifs d'une phase de ladite foulée ou de la foulée complète, en trois dimensions, dans un référentiel global fixe (GF) lié au référentiel terrestre. 20