FR 3042266 A1 20170414 FR 1559591 A 20151008 1 DOMAINE TECHNIQUE GENERAL La présente invention concerne le domaine de la navigation sans GPS. Plus précisément, elle concerne un procédé d'estimation du mouvement d'un piéton en marche par des techniques magnéto-inertielles. ETAT DE L'ART Il est aujourd'hui commun pour un piéton de suivre sa position par GPS ou en utilisant un réseau de communication (triangulation à l'aide de bornes émettrices, réseau wifi ou autres). Il est possible de leur associer d'autres capteurs pour améliorer le positionnement par exemple des capteurs barométriques, de champ magnétique, d'image, de radar, etc. Ces méthodes s'avèrent très limitées car elles ne fonctionnent pas en intérieur, dans les tunnels, ou trop loin des émetteurs, et s'avèrent dépendantes de technologies extérieures comme les satellites pour le GPS qui peuvent être indisponibles voire volontairement brouillées. Alternativement, on connait aussi des méthodes « autonomes » pour suivre dans n'importe quel environnement le déplacement relatif d'un véhicule lourd tel qu'un avion de chasse ou de ligne, un sous-marin, un navire, etc., grâce à une centrale inertielle ou magnéto-inertielle. Par déplacement relatif, on entend la trajectoire du véhicule dans l'espace par rapport à un point et à un repère donnés à l'initialisation. En plus de la trajectoire, ces méthodes permettent également d'obtenir l'orientation du véhicule par rapport au même repère initial. Une centrale inertielle est constituée au minimum de trois accéléromètres et de trois gyromètres disposés en triaxe. Typiquement, les gyromètres « maintiennent » un repère, dans lequel une double intégration temporelle des mesures des accéléromètres permet d'estimer le mouvement. 3042266 2 Il est notablement connu que pour pouvoir utiliser les méthodes de navigation inertielle classique, telles que mises en oeuvre dans les applications lourdes comme la navigation des avions de chasse ou de ligne, des sous-marins, des navires, etc., il est nécessaire d'utiliser des capteurs 5 de très haute précision. En effet la double intégration temporelle d'une mesure d'accélération fait qu'une erreur constante d'accélération crée une erreur de position qui augmente de façon proportionnelle au carré du temps. Et ces capteurs de haute précision sont trop lourds, trop encombrants et trop chers pour être embarqués par un piéton. 10 Pour pouvoir estimer une trajectoire avec des capteurs inertiels léger et à bas coûts tels que ceux employés dans les téléphones portables par exemple, des méthodes différentes doivent être mises en oeuvre ne nécessitant pas l'intégration des capteurs inertiels sur des périodes temporelles longues. 15 Une première méthode très répandue consiste à compter les pas effectués. La détection de chaque pas est réalisée en repérant un motif caractéristique d'un pas dans les mesures inertielles. Une estimation de la direction de la marche est obtenue séparément de la vitesse ou de la 20 distance, par exemple en s'orientant au nord magnétique terrestre à l'aide de capteurs sensibles au champ magnétique. Dans de nombreux endroits, de fortes perturbations magnétiques rendent la détermination du cap magnétique terrestre imprécise. Ces perturbations sont particulièrement fréquentes à l'intérieur de bâtiments du fait de la présence de matériaux 25 magnétiques par exemple dans le mobilier, les murs, installations électriques et objets divers, etc. Pour cette raison des méthodes complémentaires ont été proposées pour s'orienter en dépit de ces perturbations. Par ailleurs, il est courant d'associer un filtre d'attitude par exemple de type Kalman étendu pour 30 combiner les mesures de champ magnétique et inertielles. Cela permet d'améliorer la précision de l'orientation et notamment du cap. 3042266 3 Après chaque pas, la position estimée du porteur est mise à jour en effectuant un déplacement de la longueur estimée d'un pas dans la direction de la marche estimée à partir du cap calculé par la centrale inertielle. Cette méthode a été utilisée en attachant la centrale inertielle à 5 divers emplacements sur un individu, par exemple au pied, à la ceinture, dans une poche, au poignet, à la main, sur des lunettes, en frontale, etc. La performance obtenue est limitée par l'estimation imprécise de la longueur du pas et la différence entre cap et direction de la marche. Cette estimation peut être améliorée à l'aide d'une harmonisation des repères et 10 en reliant par un modèle la longueur et la fréquence des pas. Néanmoins une incertitude significative demeure car, d'une part, un piéton n'a jamais deux pas de longueur exactement identiques, et d'autre part, un piéton effectue différents types de marches avec des pas de longueur variables non distinguables, aucun modèle ne donne aujourd'hui satisfaction. 15 Une seconde méthode consiste à intégrer les mesures d'accélération et de vitesse angulaire sur des périodes de temps très courtes pour déterminer la trajectoire de la centrale inertielle et ainsi du porteur de la centrale. Cette méthode ne nécessite pas d'estimation de la longueur de 20 pas mais est limitée par l'accumulation des erreurs d'intégration dont on a déjà parlé. Pour des capteurs de type MEMS, cela se traduit par une erreur du même ordre de grandeur que la longueur d'un pas après quelques secondes. Une solution proposée dans la littérature (« Pedestrian tracking with 25 shoe-mounted inertial sensors. Computer Graphics and Applications », Foxlin, E., 2005) consiste à réinitialiser la vitesse calculée lorsque l'on sait qu'elle est nulle, méthode souvent appelée ZUPT (Zero velocity UPdaTe). Le pied en particulier a une vitesse nulle au contact du sol. Donc en plaçant la centrale inertielle sur la chaussure et en détectant la phase où le pied est 30 au contact du sol, il est possible de réinitialiser la vitesse à zéro. L'intégration de l'accélération est alors nécessaire uniquement pendant la 3042266 4 phase où le pied est en l'air, phase qui dure environ une seconde. Dès que le pied est au sol la vitesse est connue, égale à zéro. Afin d'améliorer cette méthode, il est possible d'utiliser un filtre d'estimation par exemple de Kalman, un filtre non linéaire ou tout autre filtre 5 pour combiner les informations des différents capteurs. Ce filtre peut par exemple comporter un état avec 6 degrés de liberté pour la vitesse et l'attitude. D'autres états peuvent être ajoutés par exemple la position, le biais des capteurs, etc. Le filtre donne également une mesure de l'incertitude des états estimés avec une matrice de covariance. Cela permet 10 de combiner aisément les mesures de capteurs additionnels pour lesquels une estimation de l'incertitude est aussi disponible. Et il peut s'agir de mesures de position absolue ou relative. Lors de la phase de ZUPT, la mise à jour des états du filtre est progressive et corrige éventuellement tous les états et pas uniquement la 15 vitesse. La méthode ZUPT permet ainsi d'améliorer la qualité de l'estimation du mouvement, mais s'avère poser un certain nombre de problèmes supplémentaires dus à la position peu pratique du capteur sur le pied. Cela le rend très sensible aux chocs (à régler en fonction type de semelle), peu 20 esthétique et inconfortable car il faut l'intégrer à la chaussure. De surcroît, le cap est difficile à déterminer près du sol avec la présence de champs magnétiques parasites et la chaussure elle-même est souvent constituée de matériaux magnétiques. Enfin se pose le problème de la course, pour laquelle la phase durant laquelle le pied est immobile au sol est très réduite, 25 voire nulle, ce qui empêche la correction des états du filtre. Il serait souhaitable de disposer d'une nouvelle méthode d'estimation du mouvement d'un piéton qui présente une qualité supérieure à celle des méthodes actuelles et ne soit pas contraignante. 30 PRESENTATION DE L'INVENTION 3042266 5 La présente invention se rapporte ainsi selon un premier aspect à un procédé d'estimation du mouvement d'un piéton en marche, le procédé étant caractérisé en ce qu'il comprend des étapes de : (a) Acquisition par des moyens de mesure inertielle solidaires d'un 5 membre inférieur dudit piéton et disposés de sorte à présenter sensiblement un mouvement de rotation par rapport à une extrémité distale dudit membre inférieur, d'une accélération et d'une vitesse angulaire dudit membre inférieur ; (b) Estimation par des moyens de traitement de données d'une vitesse 10 dudit membre inférieur en fonction de ladite accélération et de ladite vitesse angulaire mesurées. (c) Détermination par les moyens de traitement de données d'un intervalle temporel de ladite marche du piéton pendant lequel ladite extrémité distale dudit membre inférieur est en contact du sol en 15 fonction de l'accélération mesurée, de la vitesse angulaire mesurée, et d'un bras de levier entre les moyens de mesure inertielle et ladite extrémité distale ; (d) Dans ledit intervalle temporel déterminé : o Calcul par les moyens de traitement de données d'une vitesse 20 attendue dudit membre inférieur en fonction de ladite vitesse angulaire mesurée et dudit bras de levier ; o Correction de la vitesse et/ou de l'orientation estimées dudit membre inférieur en fonction de la vitesse attendue ; (e) Estimation par les moyens de traitement de données du mouvement 25 du piéton en fonction de la vitesse estimée dudit membre inférieur. Selon d'autres caractéristiques avantageuses et non limitatives : - la vitesse dudit membre inférieur est estimée par intégration de l'accélération mesurée exprimée dans le référentiel terrestre en fonction de 30 la vitesse angulaire mesurée, et le mouvement dudit membre inférieur est estimé par intégration de la vitesse estimée ; 3042266 6 - ladite vitesse attendue est donnée par la formule Vatt = (73 n r, où FI; est la vitesse angulaire mesurée et r- le bras de levier ; - lesdits moyens de mesure inertielle sont disposés sur ledit membre inférieur entre une cheville et un genou ; 5 - l'étape (c) comprend le calcul par les moyens de traitement de données d'une accélération attendue en fonction de ladite vitesse angulaire mesurée et du bras de levier, ledit intervalle temporel de ladite marche du piéton pendant lequel ladite extrémité distale dudit membre inférieur est en contact du sol étant déterminé en fonction de l'accélération mesurée et de ladite 10 accélération attendue ; - ledit intervalle temporel de ladite marche du piéton pendant lequel ladite extrémité distale dudit membre inférieur est en contact du sol est déterminé à l'étape (c) comme celui pendant lequel un écart entre l'accélération mesurée et l'accélération attendue est inférieur à un seuil prédéterminé ; 15 - ladite accélération attendue est donnée par la formule yatt = -É + -d: A É> + FI; A FI; A É>, où FI; est la vitesse angulaire mesurée, r- le bras de levier et g l'accélération de la pesanteur ; - ledit intervalle temporel de ladite marche du piéton pendant lequel ladite extrémité distale dudit membre inférieur est en contact du sol est déterminé 20 à l'étape (c) comme celui pendant lequel l'accélération et/ou la vitesse angulaire mesurées correspondent à un motif prédéterminé représentatif du contact de l'extrémité du membre inférieur avec le sol ; - la correction d'une vitesse estimée dudit membre inférieur à l'étape (c) comprend la mise en oeuvre d'un filtre estimateur d'état linéaire ou non 25 linéaire ; - le procédé comprend une étape (a0) préalable de détermination dudit bras de levier ; - le bras de levier est déterminé en minimisant l'écart entre l'accélération mesurée et une accélération attendue fonction de ladite vitesse angulaire 30 mesurée et du bras de levier, lors d'un intervalle temporel prédéterminé de 3042266 7 ladite marche du piéton pendant lequel ladite extrémité distale dudit membre inférieur est en contact du sol ; - ladite détermination du bras de levier comprend l'intégration du bras de levier au filtre estimateur d'état linéaire ou non linéaire ; 5 - le procédé comprend une étape (f) subséquente d'analyse par les moyens de traitement de données du mouvement estimé de sorte à identifier un trouble de la marche dudit piéton. Selon un deuxième aspect, l'invention concerne un équipement 10 d'estimation du mouvement d'un piéton en marche, caractérisé en ce qu'il comprend des moyens de traitement de données configurés pour mettre en oeuvre : - Un module de réception d'une accélération et d'une vitesse angulaire d'un membre inférieur dudit piéton acquises par des moyens de 15 mesure inertielle solidaires dudit membre inférieur et disposés de sorte à présenter sensiblement un mouvement de rotation par rapport à une extrémité distale dudit membre inférieur ; - Un module d'estimation d'une vitesse dudit membre inférieur en fonction de ladite accélération et de ladite vitesse angulaire 20 mesurées ; - Un module de détermination d'un intervalle temporel de ladite marche du piéton pendant lequel ladite extrémité distale dudit membre inférieur est en contact du sol en fonction de l'accélération mesurée, de la vitesse angulaire mesurée, et d'un bras de levier 25 entre les moyens de mesure inertielle et ladite extrémité distale ; - Un module de calcul dans ledit intervalle temporel déterminé d'une vitesse attendue dudit membre inférieur en fonction de ladite vitesse angulaire mesurée et dudit bras de levier ; - Un module de correction dans ledit intervalle temporel de la vitesse 30 estimée dudit membre inférieur en fonction de la vitesse attendue ; - Un module d'estimation du mouvement du piéton en fonction de la vitesse estimée dudit membre inférieur. 3042266 8 Selon d'autres caractéristiques avantageuses et non limitatives : - L'équipement est un boitier comprenant les moyens de mesure inertielle ; 5 - L'équipement comprend en outre des moyens d'attache du boitier au membre inférieur, un magnétomètre, et des moyens de communication. - L'équipement est un terminal mobile ou un serveur, adapté pour communiquer avec un boitier comprenant les moyens de mesure inertielle. 10 Selon un troisième aspect, l'invention concerne un système comprenant l'équipement selon le deuxième aspect de l'invention et au moins un boitier en connexion. Selon un quatrième et un cinquième aspect, l'invention concerne un 15 produit programme d'ordinateur comprenant des instructions de code pour l'exécution d'un procédé d'estimation du mouvement d'un piéton en marche selon le premier aspect de l'invention ; et un moyen de stockage lisible par un équipement informatique sur lequel un produit programme d'ordinateur comprend des instructions de code pour l'exécution d'un procédé 20 d'estimation du mouvement d'un piéton en marche selon le premier aspect de l'invention. PRESENTATION DES FIGURES 25 D'autres caractéristiques et avantages de la présente invention apparaîtront à la lecture de la description qui va suivre d'un mode de réalisation préférentiel. Cette description sera donnée en référence aux dessins annexés dans lesquels : - la figure 1 est un schéma d'équipements pour la mise en oeuvre du 30 procédé selon l'invention ; - La figure 2 représente plus en détail un exemple de boitier pour la mise en oeuvre du procédé selon l'invention ; 3042266 9 - les figures 3a-3b représentent schématiquement des étapes successives lors d'une phase de mise en contact du pied avec le sol ; - les figures 4a-4b représentent des exemples de résultats 5 d'estimation du mouvement d'un piéton obtenus grâce à la mise en oeuvre du procédé selon l'invention. DESCRIPTION DETAILLEE 10 Architecture En référence à la figure 1, le présent procédé permet l'estimation du mouvement d'un piéton 1 en marche. Le piéton présente au moins un membre inférieur 10 (i.e. une jambe) équipée de moyens de mesure 15 inertielle 20. On comprendra que chacun des deux membres inférieurs 10 du piéton 1 peut être équipé de moyens de mesure inertielle 20. Plus précisément, les moyens de mesure inertielle 20 sont solidaires de ce membre inférieur 10, i.e. qu'ils présentent un mouvement sensiblement identique dans le référentiel terrestre, on verra comment plus 20 loin. Les moyens de mesure inertielle sont plus précisément disposés dans la région dite crurale de la jambe inférieure 10, c'est-à-dire la moitié basse qui s'étend entre une cheville 12 et un genou 13 du piéton (inclus), et de façon générale tout emplacement présentant par rapport à une extrémité distale 11 du membre inférieur 10 (i.e. le talon du pied) sensiblement 25 uniquement un mouvement de rotation, c'est à dire dû à un bras de levier. Ainsi, lorsque le talon du pied 11 est posé, les moyens 20 ne peuvent quasiment faire qu'une rotation dans le référentiel terrestre, et non une translation. On verra plus loin l'intérêt d'un tel positionnement. Pour résumer, les moyens 20 sont typiquement disposés au niveau d'un tibia du 30 piéton 1. La distance entre le point de rotation de l'extrémité distale 11 et les moyens 20 est appelée dans la suite de la présente description « bras de 3042266 10 levier ». Et on note r- le vecteur représentant ce bras de levier orienté de l'extrémité 11 aux moyens 20. Les moyens de mesure inertielle 20 sont préférentiellement ceux d'un 5 boitier 2 tel que représenté sur la figure 2 présentant des moyens 23 d'attache au membre inférieur 10. Ces moyens d'attache 23 consistent par exemple en un bracelet par exemple à bande autoagrippante qui enserre le membre et permet la liaison solidaire. Comme l'on verra plus tard, il est en effet souhaitable que les moyens de mesure inertielle 20 soient disposés le 10 plus près possible du genou 13, et ne puissent pas se déplacer le long du membre 10. Par moyens de mesure inertielle 20, on entend une centrale à inertie comprenant au moins trois accéléromètres et trois gyromètres disposés en triaxe. Les gyromètres mesurent la vitesse angulaire instantanée de la 15 centrale inertielle par rapport au référentiel terrestre, notée Les accéléromètres sont sensibles aux forces extérieures autres que gravitationnelles appliquées sur le capteur, et permettent de mesurer une accélération notée ÿ. Comme on le verra, les moyens 20 sont avantageusement accompagnés d'au moins un magnétomètre 24 de sorte 20 à former des moyens 20, 24 de mesure magnéto-inertielle. Le magnétomètre 24 mesure un champ noté B. Un tel magnétomètre 24 est utile pour indiquer un cap du piéton (i.e. une direction dans un plan horizontal), en particulier au démarrage puisque comme expliqué le mouvement est relatif. Le magnétomètre 24 n'est ensuite plus indispensable 25 mais peut servir à recaler le cap qui dévie suite à l'accumulation d'erreurs liées aux mesures de vitesse angulaire. Le boitier 2 peut comprendre des moyens de traitement 21 (typiquement un processeur) pour la mise en oeuvre directement en temps réel des traitements du présent procédé, ou bien les mesures peuvent être 30 émises via des moyens de communication 25 vers un dispositif externe tel qu'un terminal mobile (smartphone) 3, voire un serveur distant 4, ou encore les mesures peuvent être enregistrées dans des moyens de stockage de 3042266 11 données 22 locaux (une mémoire par exemple de type flash) mémoire locale pour un traitement a posteriori par exemple sur le serveur 4. Les moyens de communication 25 peuvent mettre en oeuvre une communication sans fil à courte portée par exemple Bluetooth ou Wifi (en 5 particulier dans un mode de réalisation avec un terminal mobile 3) voire être des moyens de connexion à un réseau mobile (typiquement UMTS/LTE) pour une communication à longue distance. Il est à noter que les moyens de communication 25 peuvent être par exemple une connectique filaire (typiquement USB) pour transférer les données des moyens de stockage de 10 données 22 locaux à ceux d'un terminal mobile 3 ou d'un serveur 4. Si c'est un terminal mobile 3 (respectivement un serveur 4) qui héberge « l'intelligence », il comprend des moyens de traitement 31 (respectivement 41) tels qu'un processeur pour la mise en oeuvre des traitements du présent procédé qui vont être décrits. Lorsque les moyens de 15 traitement utilisés sont ceux 21 du boitier 2, celui-ci peut encore inclure des moyens de communication 25 pour transmettre la position estimée. Par exemple la position du porteur peut être envoyée au terminal mobile 3 pour afficher la position dans une interface d'un logiciel de navigation. Dans la suite de la présente description, on verra que les moyens de 20 traitement de données 21, 31, 41 respectivement du boitier 2, d'un smartphone 3 et d'un serveur distant 4 peuvent indifféremment et selon les applications réaliser tout ou partie des étapes du procédé. Principe et notations 25 Dans une première étape (a), le procédé comprend l'acquisition par les moyens de mesure inertielle 20 de l'accélération ir> et de la vitesse angulaire FI; dudit membre inférieur 10. On note que dans la suite de la description, lorsqu'on mentionne accélération/vitesse/position du membre 30 inférieur 10, on entend au niveau des moyens de mesure inertielle 20, car des points du membre inférieur 10 situé par exemple au niveau de la cuisse auront un mouvement différent. 3042266 12 Ces grandeurs sont avantageusement mesurées avec un échantillonnage dt (i.e. toutes les « dt » secondes) avec dt très petit devant le temps caractéristique des mouvements du piéton 1, typiquement 40 ms. L'orientation des moyens 20 par rapport au référentiel inertiel 5 terrestre peut être donnée par exemple par une matrice de rotation (notée R), un quaternion d'attitude (noté q), l'attitude est synonyme d'orientation dans l'espace ou des angles d'Euler (roulis cp, tangage 0, lacet i.p). Ces trois représentations sont équivalentes, on les utilise donc indifféremment dans ce document. La vitesse et la position des moyens 20 (et donc du membre 10 inférieur 10) sont notées respectivement v et à, et s'estiment respectivement par une simple et une double intégration de l'accélération dans le référentiel terrestre, qui comme l'on verra plus loin se calcule à partir de l'accélération ir> mesurée (donnée dans le référentiel mobile des moyens de mesure inertielle 20) et de l'orientation des moyens 20 par 15 rapport au référentiel terrestre (mise à jour à partir des mesures de vitesse angulaire). L'initialisation de l'attitude peut se faire par exemple à partir des mesures d'accélération (et le cas échéant des mesures d'un éventuel magnétomètre 24) en considérant que le membre 10 et donc les moyens 20 20 sont immobiles au démarrage et que le champ magnétique est égal au champ magnétique terrestre. Dans ce cas l'accélération mesurée est égale à l'opposée du champ gravitationnel ir> = -g. On a alors le roulis et le tangage avec les formules suivantes: cp = - tan-1 Yy Yz 0 = sin-1 Yx Iyx2 yy2 yz2 Le cap magnétique peut alors se calculer à partir de la mesure de 25 champ magnétique avec la formule: Bz - sin cp - By - cos (1) 11J = tan-1 Bx - cos 0 + By - sin 0 - sin cp + Bz - sin 0 - cos cp 3042266 13 La formule donnant la matrice de passage du référentiel terrestre au référentiel des moyens 20 à partir des angles d'Euler est: Ri,b cos 0 - cos 11J - cos 0 - sin 11J [ sin 0 = cos 11J - sin 0 - sin cp + cos cp - sin 11J cos cp - cos 11J - sin 0 - sin cp - sin 11J -cos 0 - sin cp sin cp - sin 11J - cos cp - cos 11J - sin 0 cos 11J - sin cp + cos cp - sin 0 - sin 11J cos 0 - cos cp La vitesse et la position des moyens 20 par rapport au référentiel 5 terrestre sont initialisées à zéro. La position initiale ne peut pas être déterminée directement à partir des mesures d'accéléromètres, gyromètres et magnétomètres, elle peut être fournie par un autre capteur (par exemple GPS) ou indiquée par l'utilisateur. On connait donc au cours du temps avec les seuls capteurs inertiels (et le cas échéant magnétiques) la position 10 relative du dispositif 2 (et donc du piéton 1), définie par rapport à la position initiale. On désignera le référentiel inertiel terrestre par l'indice i et le référentiel des moyens 20 aussi appelé body par l'indice b. Ainsi la matrice de changement de base du référentiel inertiel terrestre vers le référentiel 15 centrale est notée Ri,b. On note Rn l'estimation de cette matrice après n pas d'échantillon. L'attitude est liée à la vitesse angulaire FI; selon l'équation différentielle sur la matrice de passage Ri,b. Les coordonnées de w sont exprimées dans la base de la centrale inertielle. 0 -Wbz Wby Ri->b = Wbz 0 -(1)b,c' Ri,b by cobx 0 20 En considérant que la période de l'échantillonnage notée dt est suffisamment petite, on peut par exemple utiliser un développement à l'ordre 1 comme approximation: Ri,b dt) = Ri,b Ri,b (t) - dt Cette approximation peut être utilisée pour mettre à jour l'estimation de la matrice Rn pour chaque mesure: 3042266 14 0 -w w y lkn+i -Rn+ wz 0 -wx - Rn - dt [ -coy cox 0 Pour mettre à jour la position et la vitesse, on utilise la mesure de l'accélération. En particulier, le procédé comprend classiquement une étape (b) d'estimation de la vitesse en fonction de ladite accélération mesurée (en pratique par intégration) et une étape (e) d'estimation du mouvement en 5 fonction de la vitesse estimée (en pratique également par intégration). Ainsi, l'accélération des moyens 20 par rapport au référentiel terrestre est donnée par: à = ir> + g Et l'approximation donnant la vitesse en fonction de l'accélération est: v(t + dt) = v(t) + a(t) - dt 10 Et l'estimation de la vitesse est ainsi mise à jour pour chaque mesure selon: Qn+1 = Qn -ân dt La mesure de l'accéléromètre est connue dans le référentiel des moyens 20 body tandis que le champ gravitationnel est connu dans le référentiel terrestre, l'expression de l'accélération dans le référentiel 15 terrestre est donc: [aix Ybx 0 I aiy1= Rb,, ' [Ybyl aiz Ybz -g Donc si on exprime la vitesse dans la base du référentiel terrestre la formule de mise à jour au nième pas d'échantillon est: [ Qix,n+1 Q t Qiy,n+1=[Q iy,n + (lk n-1 '[Yby,n+[0 dt Qiz,n+1 Qiz,n Ybz,n -g - -x,n Ybx,n 0 1) Et si on exprime la vitesse dans la base du référentiel de la centrale, la formule de mise à jour au nième pas d'échantillon est: [Qbx,n+1 Vbx,n Yhx.nl f 0 Vby,n+1 = en+1.en-1 . Vby,n +en+1. Qbz,n+1 Qbz,n Ybz,n -g [Rn-1 - Yby,n + o - dt 3042266 15 Ce qui devient après substitutions et à l'ordre 1 en dt: , / 0 Qbx,n+1 Vbx,n Wz COy Vbx,n bxn [ 0I \ Vby,n+1 = Vby,n + [ (1)z 0 (1)x ' Vby,n + Yby,n+ Rn - [ 0 - dt Qbz,n+1 Vbz,n \ Wy Wx 0 Vbz,n Y bz,n -9 / Enfin l'estimation de la position est mise à jour selon: [a ix,n+1 [a ix,n Q. a iy,n+1 = aiy,n + Qiy,n ' dt Qt.z,n a iz,n+1 a iz,n De façon préférée, avec un filtre de Kalman étendu (on reviendra la dessus plus loin), la matrice de covariance est mise à jour en linéarisant les 5 formules précédentes autour du point courant. Condition pied au sol Le présent procédé utilise astucieusement la condition « pied au 10 sol » sur laquelle est basée la méthode ZUPT, mais sans la nécessité d'avoir le capteur sur le pied. Comme l'on voit sur la figure 3a, le pied noir au contact du sol a une vitesse nulle dans la phase représentée pour les étapes 2 à 4. Les moyens de mesure inertielle 20 étaient statiques lorsque le pied est au sol dans la 15 méthode ZUPT (et donc on pouvait réinitialiser la vitesse) La figure 3b représente cette fois un placement possible des moyens de mesure inertielle 20 au-dessus de la cheville 12 pendant la phase de contact talon/sol. Dans le cas du présent procédé les moyens de mesure inertielle 20 réalisent une rotation, et donc présentent une « vitesse 20 attendue » des moyens 20 (i.e. dudit membre inférieur 10) qui va permettre un recalage. Ladite vitesse attendue est fonction de la vitesse angulaire mesurée et dudit bras de levier, en particulier se calcule à partir du modèle jambe en rotation autour de l'extrémité 11 fixe sur le sol. Elle est alors égale au produit vectoriel entre le vecteur rotation instantanée FI; et le vecteur r- : fy> = 25 3042266 16 Une difficulté par rapport à la méthode ZUPT est d'arriver à déterminer un intervalle temporel de ladite marche du piéton 1 pendant lequel ladite extrémité distale 11 dudit membre inférieur 10 est en contact avec sol, et plus exactement le moment où le modèle donnant la vitesse est 5 correct. Pour cela, une étape (c) permet de déterminer cet intervalle temporel de contact du sol en fonction de l'accélération mesurée, de la vitesse angulaire mesurée, et du bras levier supposé ou estimé. Dans un premier mode de réalisation, cette étape (c) comprend le 10 calcul par les moyens de traitement de données 21, 31, 41 d'une « accélération attendue » en fonction de ladite vitesse angulaire mesurée et du bras de levier. En effet, avec le même modèle qui suppose les moyens 20 en rotation à la distance r (i.e. la longueur du bras de levier) du point fixe 11, 15 l'accélération théorique est égale à: aatt=- dt Ar+wnwAr En ajoutant le terme dû au champ de gravité la mesure d'accélération attendue est égale à: Y att = -g> + ci+7>) + (73 A (73 A É> Cette relation donne un critère pour déterminer si le modèle de bras de levier est correct, et l'étape (c) comprend la détermination de l'intervalle 20 temporel de contact au sol en fonction de l'accélération mesurée et de ladite accélération théorique. Plus précisément, on considère le modèle valide lorsque l'écart entre la mesure attendue de l'accélération et la mesure réalisée par les moyens de mesure inertielle 20 est nul: Ecart =1 - Y att =+g-- dt Ar-wAwAr 25 Le champ de gravité est connu dans le référentiel terrestre tandis que les mesures du gyromètre et de l'accéléromètre ainsi que r- sont connus dans la base de la centrale inertielle. Pour calculer l'écart vectoriel, il faut 3042266 17 donc effectuer un changement de base. Pour ne pas dépendre de la matrice de changement de base, on peut comparer les normes des vecteurs Ecart = cl(73 - Igl y-d t - (73 A (73 A É> dc7i Le terme -dt r- est en général petit par rapport aux autres termes, il 5 est donc possible pour éviter d'avoir à calculer la dérivée du la vitesse angulaire de la négliger. La formule de l'écart se simplifie alors: Ecart = i - -Fo>AFo> A Igl Cet écart n'est en pratique pas exactement nul et on sélectionne la période durant laquelle il est inférieur à un seuil prédéterminé, par exemple 10% de la En d'autres termes, ledit intervalle temporel de ladite marche du 10 piéton 1 pendant lequel ladite extrémité distale 11 dudit membre inférieur 10 est en contact du sol est déterminé à l'étape (c) comme celui pendant lequel un écart entre l'accélération mesurée et l'accélération théorique est inférieur à ce seuil prédéterminé. Alternativement, dans un second mode de réalisation on peut 15 rechercher grossièrement la phase de contact avec le sol à partir d'un motif caractéristique sur le signal des capteurs. La rotation de la jambe lors du balancement ou le choc lorsque le pied touche le sol par exemple peuvent permettre de déterminer approximativement ledit intervalle de contact au sol. Par exemple, lorsque l'accélération d'un choc est détectée, on peut 20 décider que pendant une durée caractéristique (par exemple un quart de seconde) suivant cette détection le pied est au sol. Alternativement, on peut dans une phase préalable d'étalonnage mesurer les valeurs classiques d'accélération et de vitesse angulaire correspondant à un contact du pied au sol, ces valeurs classiques 25 (ponctuelles ou sur un intervalle) constituant le motif de référence, et identifier les intervalles de temps minimisant l'écart entre les valeurs mesurées et ce motif de référence. 3042266 18 Dans tous les cas, on peut alors chercher dans ledit intervalle le moment où l'écart défini à partir de l'accélération est minimum. De cette façon, le meilleur moment pour effectuer le recalage est identifié. Par exemple, on peut utiliser une méthode alternative basée sur un motif 5 prédéterminé pour identifier grossièrement l'intervalle de contact au sol, puis minimiser l'écart entre l'accélération mesurée et l'accélération attendue au sein de cet intervalle (i.e. utiliser la méthode principale décrite ci-avant). Dans ledit intervalle temporel déterminé (et préférentiellement au point identifié dans cet intervalle), une étape (d) permet le recalage de la 10 vitesse estimée. Plus précisément, est calculée par les moyens de traitement de données 21, 31, 41 ladite vitesse attendue du membre inférieur 10 en fonction de ladite vitesse angulaire mesurée et dudit bras de levier. En particulier, comme expliqué elle est égale au produit vectoriel entre le 15 vecteur rotation instantanée FI; et le vecteur É> : Barr = Fo> A É> La vitesse estimée dudit membre inférieur 10 est alors corrigée en fonction de la vitesse attendue. La différence entre cette vitesse attendue et la vitesse estimée Q 20 obtenue par intégrations successives est utilisée pour le recalage. Il est tout à fait possible de recaler en remplaçant simplement la valeur de la vitesse estimée par la vitesse attendue, mais de façon préférée et comme expliqué, un filtre estimateur d'état linéaire (filtre Luenberger, Filtre de Kalman, etc.) ou non-linéaire (filtre de Kalman étendu, observateur invariant, etc.) est 25 utilisé. Dans la présente description, on prendra l'exemple d'un filtre de Kalman étendu, mais l'homme du métier saura transposer à d'autres filtres. Le gain de Kalman Kn+1 est calculé à partir de la matrice de covariance selon la formule d'un filtre de Kalman étendu. On peut modéliser les erreurs dues aux capteurs et aux approximations comme un bruit de 30 distribution gaussienne. La variance est estimée en mesurant le bruit des capteurs au repos. Le recalage de l'estimation de l'état 5s( -n+1,n contenant la 3042266 19 matrice de changement de base, la vitesse et la position se fait en ajoutant le terme de correction: EWx Lx [bx,n+11\ 5(n+1,n+1 = 5(n+1,n Kn+1 - Y A ry Vby,n+1 rz Qbz,n+1 On note qu'en pratique l'erreur de vitesse permet également de corriger l'orientation. 5 Lorsque le pied quitte le contact du sol, l'attitude et la vitesse sont à nouveau uniquement mises à jour à partir des valeurs mesurées de l'accélération et de la vitesse angulaire fournies par les moyens 20. Les erreurs de calcul de la vitesse dues à la double intégration des accéléromètres sont ainsi réinitialisées à chaque pas. 10 Calcul du bras de levier On note que la longueur du bras de levier peut être une constante saisie par l'utilisateur (le cas échéant après mesure). 15 Mais alternativement, on ne connaît pas forcément la position exacte du boitier 2 sur la jambe donc r- peut ne pas être connu a priori. Dans ce cas, le procédé comprend avantageusement une étape préalable (a0) de détermination dudit bras de levier. Par ailleurs on note que le boitier 2 peut légèrement bouger le long du membre inférieur 10, donc l'étape (a0) peut 20 être mise à nouveau en oeuvre (à intervalles réguliers ou sur instructions du piéton 1) lors de la marche pour redéterminer r. Il est en effet possible d'estimer la valeur de ce vecteur soit séparément soit en intégrant directement r- à l'état du filtre de Kalman. Si on suppose que les moyens de mesure inertielle 20 sont montés de façon à ce 25 que l'axe z soit aligné avec l'axe de la jambe, le bras de levier r- peut être considéré comme aligné avec l'axe z et les coordonnées x et y sont nulles: 0 [ryx = 01 rz rZ En particulier, le produit vectoriel donnant la vitesse devient: 3042266 20 CO [rx [.ox A 0 coy - rz cy A ry = COy [0 = [-COx - rz coz rz COz rz 0 En ajoutant un champ rz à l'état du filtre de Kalman étendu, un recalage de rz est effectué en fonction de l'erreur de vitesse mesurée à chaque phase de pied au sol. Alternativement on peut aussi d'estimer r- directement lorsque le pied 5 est au contact du sol. On a expliqué que cette phase peut se déterminer par exemple à partir du choc subi lorsque le pied touche le sol, ou alors il s'agit d'un intervalle prédéterminé où l'utilisateur pose sciemment le pied au contact du sol. Durant cette phase, on a expliqué que la mesure de l'accéléromètre 10 est donnée par: [ax [a - az - d(73 a = dtA i-> + (73 A(73A i-> 7 dwy - dt dwx dt 0 - -Wx2 rz (-1-) Y- 21 - rz I On peut ainsi chercher la valeur de rz qui minimise en moyenne la quantité: 7 dcoy - dt dwx odt 21 2 - rz x rz (-1-) Y- Enfin si on ne connait pas l'orientation des moyens de mesure inertielle 20, on peut ajouter les trois coordonnées du bras de levier au filtre. 15 Le recalage concerne alors les trois coordonnées. La position des moyens 20 sur le membre inférieur 10 peut également être estimée grâce à l'utilisation complémentaire d'autres capteurs lorsqu'ils sont disponibles, par exemple un GPS donnant une vitesse et l'orientation de la centrale par rapport à la terre, l'inconnue étant 20 alors le bras de levier, ou encore un système de vision donnant une vitesse 3042266 21 et l'orientation. En fonction de la précision de l'information par rapport à la vitesse exacte au niveau du capteur, un temps de filtrage plus ou moins long sera nécessaire. 5 Résultats La figure 4a représente un exemple d'une trajectoire obtenue après intégration des mesures d'un boiter 2 porté au-dessus de la cheville 12. La marche est constituée d'un aller-retour dans une rue avec descente puis 10 montée de deux escaliers. Les accéléromètres et les gyromètres utilisés dans les moyens 20 sont des capteurs MEMS précisément calibrés. Un filtre de Kalman étendu a été utilisé avec un recalage utilisant le modèle du bras de levier lorsque le pied est au sol comme précédemment décrit. On constate l'absence de dérive sur une marche de plusieurs 15 centaine de mètres (le piéton 1 revient exactement au point de départ) et la qualité et la précision de la navigation. La figure 4b représente plus précisément la descente et montée des escaliers au cours de la marche illustrée sur la figure 4a. On montre ici la hauteur du boitier 2 en fonction de la distance parcourue. Chaque foulée sur 20 le plat se reconnaît distinctivement, la cheville 12 s'élevant d'environ 18 cm en effectuant une trajectoire courbe. Lors de la descente des escaliers entre les distances marquées 118 m et 125 m, on reconnaît chaque foulée, le pied 11 franchissant deux marches et sur une hauteur d'environ 30 cm. De même entre les distances marquées de 128 m à 132 m, on reconnaît pour 25 chaque foulée le franchissement de deux marches. Ainsi, alors même qu'on est dans un environnement très problématique pour une centrale magnéto-inertielle (montée et descente d'escaliers, au lieu de rester sur du plat), on constate que la qualité offerte par le présent procédé reste irréprochable. 30 Ce dernier s'avère optimal pour suivre le parcours d'un piéton 1 dans une zone peu ou mal couverte par un signal GPS, à l'intérieur de bâtiments, en sous-sol, près de hautes parois, en forêt, etc. 3042266 22 Comme expliqué avant, une telle vue montrant pour chaque pas l'élévation des moyens 20 permet de détecter facilement des troubles de la marche. 5 Equipements et système Selon un deuxième aspect, l'invention concerne en particulier les équipements 2, 3, 4 pour la mise en oeuvre de l'un ou l'autre des modes de réalisation du procédé. 10 Comme expliqué précédemment, selon un premier mode de réalisation l'équipement est un boitier 2 autonome comprenant les moyens de mesure inertielle 20 et les moyens de traitement de données 21 configurés pour la mise en oeuvre des étapes du procédé. Le boitier 2 comprend en outre des moyens d'attache 23 du boitier 2 15 au membre inférieur 10, et le cas échéant un magnétomètre 24, des moyens de stockage de données 22 (pour le stockage des accélérations/vitesse angulaires mesurées ou des mouvements estimés) et/ou des moyens de communication 25 pour l'exportation des résultats. Selon un deuxième mode de réalisation, l'équipement est un terminal 20 mobile 3 ou un serveur 4, adapté pour communiquer avec un boitier 2 comprenant les moyens de mesure inertielle 20. En d'autres termes, le terminal 3 ou le serveur 4 comprend les moyens de traitement 31 ou 41 configurés pour la mise en oeuvre des étapes du procédé. Chaque boitier 2 peut tout de même comprendre des moyens de traitement de données 21 25 pour le contrôle des moyens 20 et la transmission (via des moyens de communication 25) des données mesurées aux moyens de traitement de données 31, 41. Il est à noter que les moyens 21, 31, 41 peuvent le cas échéant se partager des étapes du procédé. Par exemple, en cas d'application 30 médicale, les moyens de traitement 21 du boitier 2 peuvent réaliser les étapes jusqu'à (e), et a posteriori les moyens 41 du serveur 4 mettent en 3042266 23 oeuvre l'étape (f) d'analyse de sorte à identifier un éventuel trouble de la marche dudit piéton 1. L'invention concerne dans ce cas également le système comprenant l'équipement 3, 4 selon ce mode de réalisation et le ou les boitiers 2 5 « satellites » en connexion Dans tous les cas, les moyens de traitement de données 21, 31, 41 de l'équipement 2, 3, 4 « principal » sont configurés pour mettre en oeuvre : - Un module de réception d'une accélération et d'une vitesse angulaire d'un membre inférieur 10 dudit piéton 1 acquises par des moyens de 10 mesure inertielle 20 solidaires dudit membre inférieur 10 et disposés de sorte à présenter sensiblement un mouvement de rotation par rapport à une extrémité distale 11 dudit membre inférieur 10 ; - Un module d'estimation d'une vitesse dudit membre inférieur 10 en fonction de ladite accélération et de ladite vitesse angulaire 15 mesurées ; - Un module de détermination d'un intervalle temporel de ladite marche du piéton 1 pendant lequel ladite extrémité distale 11 dudit membre inférieur 10 est en contact du sol en fonction de l'accélération mesurée, de la vitesse angulaire mesurée, et d'un bras 20 de levier entre les moyens de mesure inertielle 20 et ladite extrémité distale 11 ; - Un module de calcul dans ledit intervalle temporel déterminé d'une vitesse attendue dudit membre inférieur 10 en fonction de ladite vitesse angulaire mesurée et dudit bras de levier ; 25 - Un module de correction dans ledit intervalle temporel de la vitesse estimée dudit membre inférieur 10 en fonction de la vitesse attendue ; - Un module d'estimation du mouvement du piéton 1 en fonction de la vitesse estimée dudit membre inférieur 10. 30 Produit programme d'ordinateur 3042266 24 Selon un troisième et un quatrième aspects, l'invention concerne un produit programme d'ordinateur comprenant des instructions de code pour l'exécution (sur les moyens de traitement 21, 31, 41) d'un procédé d'estimation du mouvement d'un piéton 1 en marche selon le premier 5 aspect de l'invention, ainsi que des moyens de stockage lisibles par un équipement informatique (par exemple des moyens de stockage de données 22) sur lequel on trouve ce produit programme d'ordinateur. REVENDICATIONS1. Procédé d'estimation du mouvement d'un piéton (1) en marche, le procédé étant caractérisé en ce qu'il comprend des étapes de : (a) Acquisition par des moyens de mesure inertielle (20) solidaires d'un membre inférieur (10) dudit piéton (1) et disposés de sorte à présenter sensiblement un mouvement de rotation par rapport à une extrémité distale (11) dudit membre inférieur (10), d'une accélération et d'une vitesse angulaire dudit membre inférieur (10) ; (b) Estimation par des moyens de traitement de données (21, 31, 41) d'une vitesse dudit membre inférieur (10) en fonction de ladite accélération et de ladite vitesse angulaire mesurées. (c) Détermination par les moyens de traitement de données (21, 31, 41) d'un intervalle temporel de ladite marche du piéton (1) pendant lequel ladite extrémité distale (11) dudit membre inférieur (10) est en contact du sol en fonction de l'accélération mesurée, de la vitesse angulaire mesurée, et d'un bras de levier entre les moyens de mesure inertielle (20) et ladite extrémité distale (11) ; (d) Dans ledit intervalle temporel déterminé : o calcul par les moyens de traitement de données (21, 31, 41) d'une vitesse attendue dudit membre inférieur (10) en fonction de ladite vitesse angulaire mesurée et dudit bras de levier ; o Correction de la vitesse estimée dudit membre inférieur (10) en fonction de la vitesse attendue ; (e) Estimation par les moyens de traitement de données (21, 31, 41) du mouvement du piéton (1) en fonction de la vitesse estimée dudit membre inférieur (10). 2. Procédé selon la revendication 1, dans lequel la vitesse dudit membre inférieur (10) est estimée par intégration de l'accélération 3042266 26 mesurée exprimée dans le référentiel terrestre en fonction de la vitesse angulaire mesurée, et le mouvement dudit membre inférieur (10) est estimé par intégration de la vitesse estimée. 5 3. Procédé selon l'une des revendications 1 et 2, dans lequel ladite vitesse attendue est donnée par la formule Vatt = (73 n r, où FI; est la vitesse angulaire mesurée et r- le bras de levier. 4. Procédé selon l'une des revendications 1 à 3, dans 10 lequel lesdits moyens de mesure inertielle (20) sont disposés sur ledit membre inférieur (10) entre une cheville (12) et un genou (13). 5. Procédé selon l'une des revendications 1 à 4, dans lequel l'étape (c) comprend le calcul par les moyens de traitement de 15 données (21, 31, 41) d'une accélération attendue en fonction de ladite vitesse angulaire mesurée et du bras de levier, ledit intervalle temporel de ladite marche du piéton (1) pendant lequel ladite extrémité distale (11) dudit membre inférieur (10) est en contact du sol étant déterminé en fonction de l'accélération mesurée et de ladite accélération attendue. 20 6. Procédé selon la revendication 5, dans lequel ledit intervalle temporel de ladite marche du piéton (1) pendant lequel ladite extrémité distale (11) dudit membre inférieur (10) est en contact du sol est déterminé à l'étape (c) comme celui pendant lequel un écart entre 25 l'accélération mesurée et l'accélération attendue est inférieur à un seuil prédéterminé. 7. Procédé selon l'une des revendications 5 et 6, dans lequel ladite accélération attendue est donnée par la formule 30 v att = -É + -d: A É> + FI; A FI; A É>, où FI; est la vitesse angulaire mesurée, r- le bras de levier et g l'accélération de la pesanteur. 3042266 27 8. Procédé selon l'une des revendications 1 à 4, dans lequel ledit intervalle temporel de ladite marche du piéton (1) pendant lequel ladite extrémité distale (11) dudit membre inférieur (10) est en contact du sol est déterminé à l'étape (c) comme celui pendant ou après lequel 5 l'accélération et/ou la vitesse angulaire mesurées correspondent à un motif prédéterminé représentatif du contact de l'extrémité (11) du membre inférieur (10).. 9. Procédé selon l'une des revendications 1 à 8, dans 10 lequel la correction d'une vitesse estimée dudit membre inférieur (10) à l'étape (c) comprend la mise en oeuvre d'un filtre estimateur d'état linéaire ou non linéaire. 10. Procédé selon l'une des revendications 1 à 9, 15 comprenant une étape (a0) préalable de détermination dudit bras de levier. 11. Procédé selon la revendication 10, dans lequel le bras de levier est déterminé en minimisant l'écart entre l'accélération mesurée et une accélération attendue fonction de ladite vitesse angulaire mesurée et 20 du bras de levier, lors d'un intervalle temporel prédéterminé de ladite marche du piéton (1) pendant lequel ladite extrémité distale (11) dudit membre inférieur (10) est en contact du sol. 12. Procédé selon les revendications 9 et 10 en 25 combinaison, dans lequel ladite détermination du bras de levier comprend l'intégration du bras de levier au filtre. 13. Equipement (2, 3, 4) d'estimation du mouvement d'un piéton (1) en marche, caractérisé en ce qu'il comprend des moyens de 30 traitement de données (21, 31, 41) configurés pour mettre en oeuvre : - Un module de réception d'une accélération et d'une vitesse angulaire d'un membre inférieur (10) dudit piéton (1) acquises par des moyens 3042266 28 de mesure inertielle (20) solidaires dudit membre inférieur (10) et disposés de sorte à présenter sensiblement un mouvement de rotation par rapport à une extrémité distale (11) dudit membre inférieur (10) ; 5 - Un module d'estimation d'une vitesse dudit membre inférieur (10) en fonction de ladite accélération et de ladite vitesse angulaire mesurées ; - Un module de détermination d'un intervalle temporel de ladite marche du piéton (1) pendant lequel ladite extrémité distale (11) 10 dudit membre inférieur (10) est en contact du sol en fonction de l'accélération mesurée, de la vitesse angulaire mesurée, et d'un bras de levier entre les moyens de mesure inertielle (20) et ladite extrémité distale (11) ; - Un module de calcul dans ledit intervalle temporel déterminé d'une 15 vitesse attendue dudit membre inférieur (10) en fonction de ladite vitesse angulaire mesurée et dudit bras de levier ; - Un module de correction dans ledit intervalle temporel de la vitesse estimée dudit membre inférieur (10) en fonction de la vitesse attendue ; 20 - Un module d'estimation du mouvement du piéton (1) en fonction de la vitesse estimée dudit membre inférieur (10). 14. Equipement selon la revendication 13, étant un boitier (2) comprenant les moyens de mesure inertielle (20). 25 15. Equipement selon la revendication 14, comprenant en outre des moyens d'attache (23) du boitier (2) au membre inférieur (10), un magnétomètre (24), et des moyens de communication (25) 30 16. Equipement selon la revendication 15, étant un terminal mobile (3) ou un serveur (4), adapté pour communiquer avec un boitier (2) comprenant les moyens de mesure inertielle (20). 3042266 29 17. Système comprenant l'équipement (3, 4) selon la revendication 16 et au moins un boitier (2) en connexion. 5 18. Produit programme d'ordinateur comprenant des instructions de code pour l'exécution d'un procédé d'estimation du mouvement d'un piéton (1) en marche selon l'une des revendications 1 à 12, lorsque ledit programme est exécuté sur un ordinateur. 10 19. Moyen de stockage lisible par un équipement informatique sur lequel un produit programme d'ordinateur comprend des instructions de code pour l'exécution d'un procédé d'estimation du mouvement d'un piéton (1) en marche selon l'une des revendications 1 à 12. 15