FR 3047070 A1 20170728 FR 1650652 A 20160127 1 PROCEDE D'ESTIMATION DE L'ACTIVITE PHYSIQUE D'UN MEMBRE SUPERIEUR DOMAINE TECHNIQUE GENERAL La présente invention concerne le domaine de la mesure d'activité physique du corps humain. Plus précisément, elle concerne un procédé d'estimation de l'activité physique d'un membre supérieur d'une personne via des techniques magnéto- inertielles. ETAT DE L'ART De nombreux dispositifs sont proposés en vente pour mesurer l'activité physique d'une personne en utilisant des capteurs sensibles à certains types de mouvements. Un exemple existant depuis de nombreuses années est le podomètre, dont le principe est de compter le nombre de pas effectués par le porteur. Plus récemment, des dispositifs ont été développés fournissant des informations plus développées sur l'activité physique que le simple comptage des pas. On peut par exemple souhaiter distinguer marche et course ou bien connaître le type de marche et l'effort associé. Au-delà de l'analyse de la marche, certains dispositifs visent à prendre en compte l'ensemble des mouvements effectués par la personne pour par exemple estimer la quantité d'énergie utilisée lors d'une activité. Pour cela, un ou plusieurs objets sont portés par la personne. De nombreuses configurations ont été proposées concernant le port de ces objets. Par exemple, ils peuvent être portés au poignet comme une montre, placés dans une poche ou à la ceinture, accrochés au torse, à une jambe, à une chaussure, etc. Mesurer l'activité physique d'une personne a des applications variées. De nombreuses personnes souhaitent par exemple augmenter ou réguler leurs dépenses caloriques notamment pour le contrôle de la masse corporelle. Dans 3047070 2 ce cadre, il peut être utile d'estimer l'activité physique effectuée par ces personnes et de réaliser un suivi en continu dans le temps. Certains sportifs peuvent chercher à optimiser leur entraînement ou simplement à suivre l'évolution de leur performance. Et indépendamment d'un objectif particulier, 5 chacun peut s'intéresser à son niveau de condition physique. Enfin dans un cadre médical, il peut être essentiel d'estimer et de suivre l'évolution de l'activité physique d'un patient. Pour des pathologies associées à des difficultés musculaires par exemple, la mesure de l'activité physique est un indicateur de l'état du patient. Le médecin peut alors suivre l'évolution de la maladie et évaluer 10 l'efficacité d'un traitement. Cela peut servir à démontrer la validité d'un traitement ou à l'adapter en fonction de la réponse du patient. Pour la mesure des mouvements, quelques capteurs fréquemment utilisés dans les dispositifs existants sont les capteurs d'accélération, de vitesse 15 angulaire, de champ magnétique, de pression atmosphérique, de position GPS, etc. En outre des mesures physiologiques sont parfois effectuées comme par exemple les pulsations cardiaques ou la température corporelle. Les premiers podomètres utilisaient ainsi un dispositif mécanique sensible aux mouvements créés lors de chaque pas. La plupart des dispositifs modernes 20 intègrent cependant des capteurs disposant d'interfaces électroniques et souvent de type MEMS. Les composants MEMS sont issus de développements récents et leur utilisation est largement répandue dans l'industrie électronique grand public. Les avantages de ces capteurs sont le faible encombrement et une consommation électrique réduite ce qui permet de les intégrer au sein d'objets 25 portés par une personne sans causer de gêne particulière. De plus ces composants sont disponibles à des coûts compatibles avec l'électronique grand public et très inférieurs à ceux de composants tactiques. Une difficulté avec de tels dispositifs concerne la fiabilité des mesures effectuées. De nombreux indicateurs sont donnés mais leur calcul n'est souvent 30 pas documenté. Et même si la méthode de calcul est expliquée, la représentativité de la variable n'est pas connue. Même pour des variables simples comme le comptage de pas, la performance est difficile à évaluer car le 3047070 3 lien avec un effort physique n'est pas direct, le pas d'un piétinement n'a rien à voir celui d'une course en sprint. On connait des méthodes pour calculer de façon nettement plus précise le déplacement relatif par intégration des capteurs d'une centrale inertielle. Une 5 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. Cependant les centrales inertielles ou magnéto-inertielles nécessaires pour cette technique sont lourdes et 10 chères telles celles embarquées dans un avion ou un sous-marin et ne peuvent pas être portées par une personne. Avec des capteurs miniaturisés de type MEMS, la performance permet de calculer un mouvement relatif pendant une durée maximale de quelques secondes. Pour cette raison des méthodes complémentaires ou alternatives ont été 15 proposées pour se repérer en dépit de la performance limitée des capteurs. Il s'agit alors d'utiliser l'information de capteurs additionnels, on peut citer par exemple la mesure du gradient de champ magnétique pour estimer la vitesse ou l'utilisation des caméras vidéos pour s'orienter par rapport à l'environnement. L'orientation par rapport à la direction verticale peut s'obtenir par la mesure du 20 champ de gravité ce dispositif est parfois appelé inclinomètre. Le cap peut s'obtenir par la mesure du champ de gravité terrestre. On constate cependant que ces mesures restent bruitées et mal exploitées. 25 La demande FR1559591 a ainsi proposé un procédé d'estimation du mouvement d'un piéton en marche grâce à des moyens de mesure inertielle solidaires d'un membre inférieur dudit piéton. Il serait souhaitable de disposer d'une méthode similaire pour estimer les efforts physiques réalisés par cette fois par un membre supérieur, et c'est l'objet 30 de la présente invention. PRESENTATION DE L'INVENTION 3047070 4 La présente invention se rapporte ainsi selon un premier aspect à un procédé d'estimation de l'activité physique exercée par un membre supérieur d'une personne, le procédé étant caractérisé en ce qu'il comprend des étapes de : 5 (a) Acquisition par des moyens de mesure inertielle solidaires d'un avant-bras dudit membre supérieur de ladite personne d'une vitesse angulaire dudit avant-bras ; (b) Estimation par des moyens de traitement de données d'un couple exercé par les muscles d'un bras du membre supérieur sur ledit avant-bras en 10 fonction de la vitesse angulaire mesurées, et de paramètres physiques dudit avant-bras ; (c) Détermination par les moyens de traitement de données d'une puissance exercée par le membre supérieur en fonction du couple estimé et de la vitesse angulaire mesurée. 15 Selon d'autres caractéristiques avantageuses et non limitatives : - Le procédé comprend suite à l'étape (a) une étape (al ) de détermination de ladite orientation de l'avant-bras ; - l'étape (a) comprend en outre l'acquisition par les moyens de mesure 20 inertielle d'une accélération spécifique dudit membre supérieur, l'orientation de l'avant-bras étant déterminée en fonction et l'accélération et la vitesse angulaire mesurées ; - la détermination de l'orientation de l'avant-bras comprend la mise en oeuvre d'un filtre estimateur d'état linéaire ou non linéaire ; 25 - lesdits moyens de mesure inertielle sont disposés sur ledit avant-bras entre un coude et un poignet ; - les paramètres physiques de l'avant-bras comprennent une longueur de l'avant-bras et un rayon de l'avant-bras ; - le couple estimé à l'étape (b) est un couple par unité de masse dudit avant- 30 bras donné par la formule 3047070 X cloL)x- 5 où (73 est la vitesse - -2 g cos 0, n, dt dwy dt doàz dt - 0 R2 Cm = R2 L2 - - 4 3 o o RZ 2 4 3 o\ 2 angulaire mesurée, R le rayon de l'avant-bras, L la longueur de l'avant-bras, g l'accélération de la pesanteur, (De l'angle d'élévation de l'avant-bras et n le vecteur unitaire normal à l'avant-bras et à la direction de l'accélération de la 5 pesanteur ; - la puissance exercée par le membre supérieur est déterminée à l'étape (c) par le produit scalaire du couple estimé et de la vitesse angulaire mesurée ; - l'étape (c) comprend l'identification du caractère moteur ou de retenue de la puissance déterminée ; 10 - le procédé est un procédé d'estimation d'une grandeur représentative de l'activité physique exercée par le membre supérieur de la personne, le procédé comprenant une étape (d) d'estimation de ladite grandeur physique à partir de la puissance déterminée ; - ladite grandeur physique est l'énergie exercée par le membre supérieur sur 15 un intervalle de temps, estimée à l'étape (d) en sommant la valeur absolue des puissances déterminées à chaque instant dudit intervalle de temps ; - le procédé comprend en outre une étape (e) de comparaison de ladite grandeur physique estimée avec au moins un seuil prédéterminé de sorte à identifier un problème musculaire du membre supérieur de la personne. 20 Selon un deuxième aspect, l'invention concerne un équipement de l'activité physique exercée par un membre supérieur d'une personne, caractérisé en ce qu'il comprend des moyens de traitement de données configurés pour mettre en oeuvre : 25 - Un module de réception d'une vitesse angulaire d'un avant-bras dudit membre supérieur de ladite personne acquise par des moyens de mesure inertielle solidaires dudit avant-bras ; - Un module d'estimation d'un couple exercé par les muscles d'un bras du membre supérieur sur ledit avant-bras en fonction de la vitesse angulaire 3047070 6 mesurées, d'une orientation dudit avant-bras, et de paramètres physiques dudit avant-bras ; - Un module de détermination d'une puissance exercée par le membre supérieur en fonction du couple estimé et de la vitesse angulaire mesurée. 5 Selon d'autres caractéristiques avantageuses et non limitatives : - L'équipement est un boitier comprenant les moyens de mesure inertielle ; - L'équipement comprend en outre des moyens d'attache du boitier au membre supérieur, et des moyens de communication. 10 - L'équipement est un terminal mobile ou un serveur, adapté pour communiquer avec un boitier comprenant les moyens de mesure inertielle. 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 15 connexion. Selon un quatrième et un cinquième aspect, l'invention concerne un produit programme d'ordinateur comprenant des instructions de code pour l'exécution d'un procédé d'estimation de l'activité physique exercée par un 20 membre supérieur d'une personne 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é d'estimation de l'activité physique exercée par un membre supérieur d'une personne selon le premier aspect de l'invention. 25 PRESENTATION DES FIGURES 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 30 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 procédé selon l'invention ; 3047070 7 - La figure 2 représente plus en détail un exemple de boitier pour la mise en oeuvre du procédé selon l'invention ; - la figure 3 représente schématiquement un membre supérieur équipé d'un boitier pour la mise en oeuvre du procédé selon l'invention. 5 DESCRIPTION DETAILLEE Architecture 10 En référence à la figure 1, le présent procédé permet l'estimation de l'activité physique exercée par un membre supérieur 10 d'une personne 1. Par « activité physique », on entend une grandeur physique représentative des efforts exercés par les muscles du membre supérieur, la grandeur physique étant choisie parmi un ensemble de grandeurs mécaniques telles que la 15 puissance instantanée, la puissance moyenne, le travail, l'énergie, etc. Ces notions seront explicitées plus loin. De façon générale, le présent procédé peut être défini comme un procédé d'estimation d'au moins une grandeur physique représentative de l'activité physique exercée par un membre supérieur 10 d'une personne 1. 20 La personne 1 présente au moins un membre supérieur 10 (i.e. un « bras » dans le langage courant, même si comme on le verra plus bas cette appellation est impropre, un bras désignant anatomiquement « l'arrière-bras », c'est-à-dire la partie entre l'épaule et le coude) équipée de moyens de mesure inertielle 20. On comprendra que chacun des deux membres supérieurs 10 de la 25 personne 1 peut être équipé de moyens de mesure inertielle 20. Les membres supérieurs 10 sont chacun constitués de trois segments : le bras 12, l'avant-bras 11 et la main. L'articulation du coude 13 relie le bras à l'avant-bras, et le poignet 14 relie l'avant-bras à la main. Plus précisément, les moyens de mesure inertielle 20 sont solidaires de 30 l'avant-bras 11, i.e. qu'ils présentent un mouvement sensiblement identique dans le référentiel terrestre, on verra comment plus loin. Les moyens de mesure inertielle sont plus précisément disposés entre le coude 13 et le poignet 14 de la personne (inclus), et de façon générale tout emplacement présentant par rapport 3047070 8 à au coude 13 du membre supérieur 10 sensiblement uniquement un mouvement de rotation, c'est à dire dû à un bras de levier. La distance entre le coude 13 et les moyens 20 importe peu comme le verra, la personne pouvant le fixer à proximité du poignet 14 pour plus de 5 distance entre le coude 13 et les moyens 20 et donc plus de précision, mais quelle que soit sa position le présent procédé fonctionne. Comme on le verra plus loin, des paramètres physiques dudit avant-bras 11 sont nécessaires, en particulier des paramètres anthropométriques dont une longueur de l'avant-bras 11 et un rayon de l'avant-bras 11. 10 Ces grandeurs peuvent être des paramètres d'entrée (saisis par la personne 1) ou alors une valeur moyenne peut être utilisée. On note que le rayon intervient dans le calcul mais comme l'on verra il peut être négligé car le terme lié au rayon est petit par rapport aux autres termes. Par ailleurs la masse de l'avant-bras 11 peut être utilisée, mais il est 15 possible de fournir des valeurs normalisées des grandeurs représentatives de l'activité physique (i.e. par unité de masse de l'avant-bras 11), ce qui fait que sa connaissance est optionnelle. Les moyens de mesure inertielle 20 sont préférentiellement ceux d'un 20 boitier 2 tel que représenté sur la figure 2 présentant des moyens 23 d'attache au membre supé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 plus près possible du poignet 25 14, et si possible 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 centrale inertielle par rapport au référentiel terrestre, notéeci. Les accéléromètres sont sensibles 30 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 peuvent comprendre des moyens 24 d'acquisition de l'orientation (dont la vitesse angulaire est la dérivée) voire de manière plus limitée seulement de 3047070 9 l'angle d'élévation (i.e. une seule composante de l'orientation) de l'avant-bras 11. Par exemple un inclinomètre qui est souvent un simple accéléromètre ne permet pas d'estimer l'orientation ni la vitesse angulaire de l'avant-bras 11 mais seulement l'angle d'élévation et avec des erreurs plus importantes qu'avec une 5 centrale inertielle complète. Il est à noter que dans un mode de réalisation où seulement l'angle d'élévation est mesuré grâce à un simple inclinomètre, cette donnée peut être « débruitée » grâce aux mesures de vitesse angulaire des moyens 20. On comprend toutefois que l'orientation complète peut être mesurée par des moyens 24 adaptés, par exemple optiques ou magnétiques. Comme on 10 le verra les moyens 24 (qu'ils mesurent une composante ou toute l'orientation) ne sont pas indispensables et seuls des moyens de mesure inertielle de la vitesse angulaire sont essentiels. 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 15 du présent procédé, ou bien les mesures peuvent être é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 données 22 locaux (une mémoire par exemple de type flash) mémoire locale pour un traitement a posteriori par 20 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 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 25 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 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 30 « 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 traitement utilisés sont ceux 21 du boitier 2, celui-ci peut encore inclure des moyens de communication 25 pour 3047070 10 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 5 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 10 Dans une première étape (a), le procédé comprend l'acquisition par les moyens de mesure inertielle 20 d'au moins la vitesse angulaire FI; dudit membre supérieur 10, et plus précisément de l'avant-bras 11. De façon préférée, l'accélération spécifique yr> de l'avant-bras 11 est 15 également mesurée. Avantageusement, le procédé comprend suite à l'étape (a) une étape (al ) de détermination par les moyens de traitement de données 21, 31, 41 de ladite orientation de l'avant-bras 11, et ce grâce à cette accélération mesurée (et la vitesse angulaire mesurée). On note que dans la suite de la description, lorsqu'on mentionne 20 accélération/vitesse/orientation du membre supérieur 10, on entend au niveau des moyens de mesure inertielle 20, car des points du membre supérieur 10 situés par exemple au niveau du bras 12 auront un mouvement différent. On introduit un référentiel inertiel terrestre (R) ainsi qu'un référentiel (Rb) ayant pour origine le centre du boitier 2 et dont les axes sont liés à ceux du 25 boitier 2. La vitesse angulaire FI; de l'avant bras 11 est égal à la rotation de (Rb) par rapport à (R). L'accélération spécifique yr> est égal à la somme de l'accélération du boitier dans le référentiel inertiel terrestre moins la valeur du champ de gravité g. Les capteurs inertiels gyromètre et accéléromètre (moyens de mesure inertielle 20) mesurent respectivement FI; et yr> exprimés dans le repère 30 (Rb). 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 de la personne 1, typiquement 10 ms. 3047070 11 L'orientation des moyens 20 par rapport au référentiel inertiel 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 5 sont équivalentes, on les utilise donc indifféremment dans ce document. 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) en considérant que le membre 10 et donc les moyens 20 sont immobiles au démarrage et que le champ magnétique est égal au champ magnétique terrestre. 10 Dans ce cas l'accélération mesurée est égale à l'opposée du champ gravitationnel ir> = On a alors le roulis et le tangage avec les formules suivantes: cp = - tan-1 Yy Yz 0 = sin-1- Yx Iyx2 yy2 yz 2 Le cap magnétique peut alors se calculer à partir de la mesure de champ magnétique avec la formule: 11J = tan-1 Bx - cos 0 + By - sin 0 - sin cp + Bz - sin 0 - cos cp 15 La formule donnant la matrice de passage du référentiel terrestre (R) au référentiel (Rb) des moyens 20 à partir des angles d'Euler est: RR1_)Rb 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 formule ci-avant permet d'estimer l'orientation du boitier 2 à tout moment à partir des mesures d'accélération et de vitesse angulaire mais une 20 limitation importante vient de l'hypothèse d'immobilité (yr> = -a). En pratique l'avant-bras 11 bouge et l'accélération n'est pas négligeable devant la valeur du champ de gravité. Pour cette raison il est opportun d'utiliser les mesures de Bz - sin cp - By - cos (1) 3047070 12 vitesse angulaire pour filtrer les erreurs liées aux accélérations subies par le boitier 2. En d'autres termes, la détermination de l'orientation de l'avant-bras 11 comprend avantageusement la mise en oeuvre d'un filtre estimateur d'état 5 linéaire ou non linéaire Les méthodes permettant d'estimer l'orientation d'une centrale inertielle sont connues. Pour cela, on peut utiliser un filtre estimateur d'état linéaire (filtre Luenberger, filtre de Kalman, etc.) ou non-linéaire (filtre de Kalman étendu, observateur invariant, etc.). Dans la présente description, on rappelle les 10 principales étapes du filtre de Kalman étendu, mais l'homme du métier saura transposer à d'autres filtres. Dans un filtre de Kalman étendu, l'état est représenté par un vecteur par exemple un quaternion d'attitude de dimension 4. Il est possible d'ajouter d'autres états à ce vecteur pour améliorer l'estimation, par exemple le biais des capteurs. 15 Une matrice de covariance est utilisée pour estimer la covariance entre chaque état du filtre, elle est de dimension n2 donc de dimension 16 dans l'exemple du quaternion d'attitude. Le filtre de Kalman s'effectue en deux étapes, une étape de prédiction et une étape mise à jour. Lors de la mise à jour, un gain de Kalman K, est calculé à partir de la matrice de covariance. 20 L'étape de prédiction repose sur un modèle lié à la vitesse angulaireci. On utilise l'équation différentielle sur la matrice de passage RR0,R, et les coordonnées de w exprimées dans la base du boitier 2. kRi->Rb =[Wbz -Wbz Wby -Wby Wbx -Wbx X RRi->Rb 0 En considérant que la période de l'échantillonnage notée dt est suffisamment petite, on utilise l'approximation d'un développement limité à l'ordre 25 1. RRi,Rb(t + dt) = Ri,b fti,b dt En notant Rn l'estimation de la matrice RRi,Rb après n pas d'échantillons, l'estimation de la matrice Rn pour chaque mesure n est: 3047070 13 -Wn,z n,y n +1 = Wn,z -Wn,x X Rn dt n,y n,x 0 L'équation est donnée ici pour une matrice de rotation mais une équation différentielle équivalente existe entre quaternion d'attitude et vitesse angulaire. En parallèle, la matrice de covariance est modifiée en linéarisant l'équation différentielle. 5 L'étape de mise à jour repose sur la mesure de l'accéléromètre. On utilise le fait que l'accélération du boitier 2 est égale à la somme de l'accélération spécifique mesurée par les moyens 20 et du champ de gravité. à = ir> + g On sait de plus que l'accélération à est en moyenne nulle (on comprend 10 que le membre supérieur 10 revient régulièrement à sa position de départ et que les mouvements dans un sens ou l'autre se compensant statistiquement). On en déduit que ir> exprimée dans le repère inertiel terrestre est égale en moyenne à -g. A partir de la matrice de rotation Rn et des mesures des moyens 20 dans 15 le repère du boitier 2, on peut exprimer l'accélération spécifique dans le repère inertiel terrestre: Yn,x [yn,y Yn,z La différence entre le champ de gravité estimé et le champ de gravité terrestre réel s'écrit alors: r n,x [ 01 erreur = Rn 1 Yn,y - 0 Yn,z g Cette différence est utilisée pour le recalage de l'attitude lors de l'étape de 20 mise à jour. En linéarisant cette formule, on parvient à calculer le gain de Kalman et par suite de mettre à jour l'état et la matrice de covariance. La mise à jour de la matrice de covariance repose sur l'hypothèse que les erreurs dues aux capteurs et aux approximations sont modélisées comme un bruit de distribution gaussienne. La variance est estimée en mesurant le bruit des 3047070 14 capteurs au repos et à partir d'hypothèses sur les mouvements effectués par l'avant-bras 11. Il est ainsi possible aux moyens de traitement de données 21, 31, 41 d'estimer l'orientation du boitier 2 (grâce aux mesures d'accélération et de vitesse 5 angulaire des moyens 20) et par conséquent celle de l'avant-bras 11. On comprendre tout autre dispositif 24 permettant cette mesure serait également utilisable, en notant toutefois qu'un inclinomètre seul mesurant l'angle d'élévation ne donne pas toute l'orientation (par exemple la torsion de l'avant-bras 11 n'est pas donnée), pour cela il faut au moins un gyromètre des moyens 20. Dans la 10 suite de la description de cette invention, il suffit que l'orientation de l'avant-bras 11 soit disponible ou même que l'angle d'élévation entre l'avant-bras 11 et la direction horizontale soit connu (avec précision). Couple contra gravitaire massique 15 Dans une étape (b), les moyens de traitement de données 21, 31, 41 estiment un couple exercé par les muscles d'un bras 12 du membre supérieur 1 sur ledit avant-bras 11 en fonction de la vitesse angulaire mesurées, d'une orientation dudit avant-bras 11, et de paramètres physiques dudit avant-bras 11. 20 Comme expliqué ce couple peut être un couple massique. Pour cela, on remarque qu'en position assise, il est naturel d'avoir le coude 13 posé sur un support par exemple une table, un accoudoir, etc., comme représenté sur la figure 3. Une grande partie de l'activité consiste ensuite à déplacer l'avant-bras 11 dans un mouvement de rotation autour du point de 25 contact entre le coude et le support. On peut noter que les personnes les plus affaiblies notamment du fait de troubles musculaires passent la majorité de leur temps assis en fauteuil avec le coude posé sur un support. Les personnes plus mobiles passent aussi beaucoup de temps avec les coudes posés sur un support par exemple un bureau, de plus il est possible d'identifier ces périodes à partir 30 des mesures inertielles. L'avant-bras 11 est modélisé comme un cylindre de longueur L, de rayon R et de masse volumique p, formant comme expliqué les paramètres physiques de l'avant-bras 11. On suppose que le coude 13 est posé sur l'accoudoir du 3047070 15 fauteuil et on note O le point de contact. On considère les rotations de l'avant-bras 11 autour du point de contact. On utilise le référentiel (Ro) ayant pour origine O et dont les axes sont ceux d'un référentiel terrestre ainsi que le référentiel (R') ayant pour origine O et dont les axes sont liés à ceux de l'avant- 5 bras 11. La vitesse angulaire FI; de l'avant-bras 11 est égal à la rotation de (Ro) par rapport à (R'). La vitesse d'un point de l'avant-bras 11 M dans (Ro) est donné par : v-> m = A OM Et par définition, le moment cinétique de l'avant-bras par rapport à O est: Lo = OM A p MdT=pf 0 M A (Fo> A OM) cit En notant p la masse volumique et en intégrant sur le volume de l'avant- 10 bras 11. En coordonnées cylindriques (é,,>, Z), avec Fz> dans l'axe de l'avant-bras 11, cela donne: L 2.1-[ R Q>,=pf f f (rF'.>+zFz>)APA(rF'.>+zFz>))rdrdodz o o o Lo = p L 2.1-[ R 10 f. ((r2 + z2) -> - CO (FO> - (r Fr> + z Fz>))(r Fr> + z Fz>)) r d, d0 dz En utilisant é, = cos 0 X + sin 0 ÿ on obtient: L f2.-rt fR Lo = p fo jo jo ((r2 + z2) Fi) - (w r cos 0 + coy r sin 0 + coz z)(r cos 0 X + r sin 0 eÿ +z)) r d, do dz En projetant sur l'axe Fz> L f2.1-[ R 10 JO 1 LOz = P ((r2 + z2) coz - (w r cos 0 + coy r sin 0 + coz z) z) r d, d0 dz 0 L 2.-rt R LOz = p (r3 coz - cox r2 z cos 0 - coy r2 z sin 0) d, d0 dz 15 Après intégration, il reste: Loz = p2iiLR4coz m R2 coz 4 2 10 10 10 3047070 16 De même en projetant sur l'axe L R Lox =P f f ((r2 + z2) cox - (r cox cos e r coy sin 0 + w z) r cos 0) r dr do dz o o 0 L 2.-rt R Lox = p f f f ((r3 sin2 6 + r z2) coX - r3 - coy sin 0 cos 0 - wz r2 z cos 0) dr de dz Après intégration, il reste: Lox - P frix ( 4 + 6 = + 1T L R4 2 1T L3 R2) (R2 L2 m COx Par symétrie, le même calcul donne pour l'axe R2 L2 Loy = (4 + -3m wY Dans le repère d'origine O et d'axes (>' ÿ,Fz>) avec Fz> dans l'axe de 5 l'avant-bras 11 (voir figure 3), on peut donc écrire le tenseur d'inertie massique: /R2 L2 4 + -3 0 01 R2 L2 o + 7 o R2 ° 0 2 Lo = m ici Par définition, le moment des forces de gravitation sur l'avant-bras 11 par rapport à O s'écrit: MID° OM A g p JJJ En coordonnées cylindriques (.>.,)>,>c) avec Fz> dans l'axe de l'avant-bras 11, cela donne: L 2-rt R MPO f f (rF,.>+zFz>)prdrdodz)Ag o o o L L f2-r[ R IVVDc = (p TC R2 f z dz) Fz> A g+ (p f f r2 .. dr do dz) A g Jo \ o o o 10 Le second terme de la somme est nul (en utilisant un argument de symétrie) et on obtient après intégration et simplifications: L Mpo = 111 -2 Fz> A g 3047070 17 En introduisant l'angle d'élévation (De et le vecteur unitaire n normal à l'avant-bras 11 et à la direction du champ de gravité tel que défini sur la figure, on a Fz> A g = COS ee Mpo = m -2 g cos Oe On utilise alors le théorème du moment cinétique en considérant que les 5 seules forces appliquées à l'avant-bras sont son poids et les forces des muscles. En notant M -Bras° le moment des forces exercées par les muscles du bras 12 sur l'avant-bras 11, on a: dLo dt MPO MBraso dLo d(73 L MBraso = dt 1\113° = 1 dt m -2 g zAg Finalement, on note Cm le couple à exercer pour lever l'avant-bras 11 divisé par sa masse et on obtient: d(73 L CM =I dt2 - g cos Oe /R2 L2 0 01 - -dcox- L 4 + -3 R2 L2 - - g cos 0e 2 0 4 + -3 0 ee ° R2 0 2 -(R2 L2) dcox- dt dw Y CM= dt dwz dt - g COS 4 + 3 dt L2) dcoy (R2 Cm = 4 + 3 dt R2 dwz 2 dt Comme expliqué, il est possible de négliger le terme portant sur le rayon devant le terme portant sur la longueur, et d'estimer le couple plus simplement par la formule : 10 3047070 18 L2 dcox 3 dt Cm = L2 dwy 3 dt 0 - - g cos ee Cette formule peut être encore simplifiée en Cm = - l=" g cos Oe en négligeant le terme portant sur la longueur, par exemple 2 dans un mode de réalisation où l'on disposerait seulement d'un inclinomètre 24 et 5 de moyens 20 capables de mesurer la vitesse angulaire (gyromètres), de sorte à calculer le couple uniquement à partir de l'angle d'élévation. Une telle valeur serait toutefois approximative et il apparait impératif de filtrer les erreurs de l'angle d'élévation avec les vitesses angulaires mesurées pour utiliser cette formule. 10 Dans une étape (c), les moyens de traitement de données 21, 31, 41 déterminent une puissance exercée par le membre supérieur 10 en fonction du couple estimé et de la vitesse angulaire mesurée. Plus précisément la puissance (instantanée) est donnée par le produit 15 scalaire du couple estimé et de la vitesse angulaire. En d'autres termes, elle est égale au couple multiplié par la vitesse angulaire. On peut le cas échéant de nouveau normaliser par la masse de l'avant-bras 11 et la puissance massique vaut: di L Pm = Cm - Fo> = 1 -dt co - -2 g cos 0 e - FO> -di i L dOe Pm = I -dt co + -2 g cos Oe dt Et en utilisant les projections sur les axes liés au boitier: /R2 L2 0 \ / -4 + -3 0 R2 L2 o .71 + 7 o R2 d(co,2 + coy2 + co,2) + L (g cox - gx coy) dt 2 Y Pm = ° 0 2 3047070 19 Comme avant cette puissance peut être calculée plus facilement en négligeant R devant L, voire en négligeant L. On remarque que dans tous les cas la puissance associée au couple exercé par le bras 12 sur l'avant-bras 11 est la somme de deux termes. Le 5 premier terme vient de l'énergie cinétique permettant de mettre en mouvement l'avant-bras 11. Le second terme vient de l'énergie potentielle liée à la hauteur de l'avant-bras 11. On remarque que le second terme qui dérive d'une énergie potentielle peut s'exprimer uniquement à partir de l'angle d'élévation. Les dispositifs fournissant uniquement l'élévation peuvent donc tirer parti de cette 10 méthode. La puissance peut être de signe négatif lorsque l'on l'abaisse l'avant-bras 11 ou lorsque sa vitesse angulaire diminue. Même dans ces les muscles doivent fournir un effort pour retenir le membre. De façon préférée, l'étape (c) comprend l'identification du caractère moteur ou de retenue de la puissance déterminée, en 15 particulier en fonction du signe de la vitesse angulaire et de l'orientation (fermeture du membre supérieur 10 pour un mouvement moteur, et ouverture du membre supérieur 10 pour un mouvement de retenue). Plusieurs modèles sont envisageables pour estimer au mieux l'effort exercé par les muscles dans ce type de mouvement. On peut considérer 20 uniquement les mouvements ayant une puissance positive car ce sont les mouvements demandant le plus d'énergie aux muscles. Ou bien on peut considérer la valeur absolue de la puissance en appliquant éventuellement un coefficient différent en fonction du signe de chaque terme constituant la puissance. Cela permet de refléter les différents types d'efforts effectués par les 25 muscles. La puissance peut être la grandeur demandée, et celle-ci peut être directement restituée/stockée par le boitier 2 ou un équipement autre tel que le terminal 3. On note que la puissance moyenne peut être calculée sur un parcours, le cas échéant la puissance moyenne des puissances positives, ou 30 encore la moyenne de la valeur absolue des puissances. Dans une étape (d), les moyens de traitement de données 21, 31, 41 peuvent déterminer une autre grandeur représentative de l'activité physique du membre supérieur 10 à partir de la puissance. 3047070 20 Par exemple peut être estimé le travail effectué par les muscles du bras lors de chaque mouvement à partir de la puissance du couple exercé par le bras 12 sur l'avant-bras 11, en sommant pour chaque échantillon la puissance massique pour connaitre le travail massique des muscles pendant la période 5 considérée (le travail peut être négatif). Alternativement, une énergie exercée par le membre supérieur sur un intervalle de temps en sommant la valeur absolue des puissances déterminées à chaque instant dudit intervalle de temps. De façon préférée le procédé comprend en outre une étape (e) de 10 comparaison de ladite grandeur physique estimée avec au moins un seuil prédéterminé de sorte à identifier un problème musculaire du membre supérieur 10 de la personne 1. Ce peut être un seuil minimal (par exemple, si le muscle de l'avant-bras n'est pas capable de dépasser une valeur de puissance massique, le patient peut souffrir d'une dégénérescence musculaire) mais également un seuil 15 maximal (par exemple, des maladies de type chorées sont définies par la survenue de mouvements incontrôlables, brusques et irréguliers). Equipements et système 20 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é. Comme expliqué précédemment, selon un premier mode de réalisation l'équipement est un boitier 2 autonome comprenant les moyens de mesure 25 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 au membre supé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 30 mesurées ou des grandeurs physiques estimées) 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 mobile 3 ou un serveur 4, adapté pour communiquer avec un boitier 2 3047070 21 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 pour le contrôle des 5 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 médicale, les moyens de traitement 21 du boitier 2 peuvent réaliser les étapes jusqu'à (c), et a 10 posteriori les moyens 41 du serveur 4 mettent en oeuvre l'étape (d) de détermination d'une grandeur physique adéquate de sorte à identifier un éventuel trouble du mouvement de la personne 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 « satellites » 15 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 vitesse angulaire d'un avant-bras 11 dudit membre supérieur 10 de ladite personne 1 (et le cas échéant de 20 l'accélération de l'avant-bras 11) acquise(s) par des moyens de mesure inertielle 20 solidaires dudit avant-bras 10 ; - Optionnellement, un module de détermination d'une orientation de l'avant-bras 11, en particulier en fonction et l'accélération et la vitesse angulaire mesurées ; 25 - Un module d'estimation d'un couple exercé par les muscles d'un bras 12 du membre supérieur 10 sur ledit avant-bras 11 en fonction de la vitesse angulaire mesurées, de l'orientation dudit avant-bras 11, et de paramètres physiques dudit avant-bras 11 ; - Un module de détermination d'une puissance exercée par le membre 30 supérieur 10 en fonction du couple estimé et de la vitesse angulaire mesurée. Produit programme d'ordinateur 3047070 22 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) de l'activité physique 5 exercée par un membre supérieur 10 d'une personne 1 selon le premier 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 de l'activité physique exercée par un membre supérieur (10) d'une personne (1), 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 avant-bras (11) dudit membre supérieur (10) de ladite personne (1) d'une vitesse angulaire dudit avant-bras (10) ; (b) Estimation par des moyens de traitement de données (21, 31, 41) d'un couple exercé par les muscles d'un bras (12) du membre supérieur (10) sur ledit avant-bras (11) en fonction de la vitesse angulaire mesurées, d'une orientation dudit avant-bras (11), et de paramètres physiques dudit avant-bras (11) ; (c) Détermination par les moyens de traitement de données (21, 31, 41) d'une puissance exercée par le membre supérieur (10) en fonction du couple estimé et de la vitesse angulaire mesurée. 2. Procédé selon la revendication 1, comprenant suite à l'étape (a) une étape (al) de détermination de ladite orientation de l'avant-bras (11). 3. Procédé selon la revendication 2, dans lequel l'étape (a) comprend en outre l'acquisition par les moyens de mesure inertielle (20) d'une accélération dudit membre supérieur (10), l'orientation de l'avant-bras (11) étant déterminée en fonction et l'accélération et la vitesse angulaire mesurées. 4. Procédé selon la revendication 3, dans lequel la détermination de l'orientation de l'avant-bras (11) comprend la mise en oeuvre d'un filtre estimateur d'état linéaire ou non linéaire. 5. Procédé selon l'une des revendications 1 à 4, dans lequel lesdits moyens de mesure inertielle (20) sont disposés sur ledit avant-bras (10) entre un coude (13) et un poignet (14). 3047070 24 6. Procédé selon l'une des revendications 1 à 5, dans lequel les paramètres physiques de l'avant-bras (11) comprennent une longueur de l'avant-bras (11) et un rayon de l'avant-bras (11). 5 7. Procédé selon la revendication 6, dans lequel le couple estimé à l'étape (b) est un couple par unité de masse dudit avant-bras (11) donné par RZ 0 2 R2 0 L2 \ la où (73 est la formule 4 0 3 4 0 3 0 - -2 g cos 0, n, vitesse R2 2 dt dwy Cm = dt doàz dt - angulaire mesurée, R le rayon de l'avant-bras (11), L la longueur de l'avant-bras 10 (11), g l'accélération de la pesanteur, (De l'angle d'élévation de l'avant-bras (11) et le vecteur unitaire normal à l'avant-bras (11) et à la direction de l'accélération de la pesanteur. 8. Procédé selon l'une des revendications 1 à 7, dans lequel la 15 puissance exercée par le membre supérieur (10) est déterminée à l'étape (c) par le produit scalaire du couple estimé et de la vitesse angulaire mesurée. 9. Procédé selon l'une des revendications 1 à 8, dans lequel l'étape (c) comprend l'identification du caractère moteur ou de retenue de la 20 puissance déterminée. 10. Procédé selon l'une des revendications 1 à 9, étant un procédé d'estimation d'une grandeur représentative de l'activité physique exercée par le membre supérieur (10) de la personne (1), le procédé comprenant 25 une étape (d) d'estimation de ladite grandeur physique à partir de la puissance déterminée. 11. Procédé selon la revendication 10, dans lequel ladite grandeur physique est l'énergie exercée par le membre supérieur (10) sur un 3047070 25 intervalle de temps, déterminée à l'étape (d) en sommant la valeur absolue des puissances déterminées à chaque instant dudit intervalle de temps. 12. Procédé selon l'une des revendications 10 et 11, comprenant 5 en outre une étape (e) de comparaison de ladite grandeur physique estimée avec au moins un seuil prédéterminé de sorte à identifier un problème musculaire du membre supérieur (10) de la personne (1). 13. Equipement (2, 3, 4) de l'activité physique exercée par un 10 membre supérieur (10) d'une personne (1), caractérisé en ce qu'il comprend des moyens de traitement de données (21, 31, 41) configurés pour mettre en oeuvre : - Un module de réception d'une vitesse angulaire d'un avant-bras (11) dudit membre supérieur (10) de ladite personne (1) acquise par des moyens de mesure inertielle (20) solidaires dudit avant-bras (10) ; 15 - Un module d'estimation d'un couple exercé par les muscles d'un bras (12) du membre supérieur (10) sur ledit avant-bras (11) en fonction de la vitesse angulaire mesurées, d'une orientation dudit avant-bras (11), et de paramètres physiques dudit avant-bras (11) ; - Un module de détermination d'une puissance exercée par le membre 20 supérieur (10) en fonction du couple estimé et de la vitesse angulaire mesurée. 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 supérieur (10) et des moyens de communication (25) 30 16. Equipement selon la revendication 13, étant un terminal mobile (3) ou un serveur (4), adapté pour communiquer avec un boitier (2) comprenant les moyens de mesure inertielle (20). 3047070 26 17. Système comprenant l'équipement (3, 4) selon la revendication 16 et au moins un boitier (2) en connexion. 18. Produit programme d'ordinateur comprenant des instructions 5 de code pour l'exécution d'un procédé d'estimation de l'activité physique exercée par un membre supérieur (10) d'une personne (1) selon l'une des revendications 1 à 12, lorsque ledit programme est exécuté sur un ordinateur. 19. Moyen de stockage lisible par un équipement informatique 10 sur lequel un produit programme d'ordinateur comprend des instructions de code pour l'exécution d'un procédé d'estimation de l'activité physique exercée par un membre supérieur (10) d'une personne (1) selon l'une des revendications 1 à 12.