Brevet FR3061787A1 — publié le 20180713

DOMAINE DE L'INVENTION La présente invention concerne un dispositif permettant l'analyse de la régularité et de la symétrie d'une séquence de cycles de marche ou de course d'une personne. ETAT DE LA TECHNIQUE Récemment, le progrès de l'électronique et de l'informatique a entraîné le développement de nouveaux capteurs qui rendent la mesure de la marche en consultation clinique de routine possible (par exemple à partir de capteurs inertiels). Ces capteurs donnent accès à des signaux qu'il faut synthétiser pour en tirer des conclusions sur la marche d'un sujet. Les paramètres de marche comme la vitesse ou la cadence constituent la forme la plus synthétisée des signaux de marche car ils les réduisent à un nombre. Ainsi, les constructeurs d'outils de mesure du pas fournissent des compte- rendus exhaustifs sous forme de liste de paramètres traduisant les différents aspects de la marche d'un sujet comme la symétrie, la régularité ou la rapidité de la mise en route. Il peut y avoir plus d'une centaine de paramètres. Ils sont représentés de manière adaptée sous forme de tableau, de diagramme en en barre ou de courbe. Néanmoins, la thématique de la visualisation des données intermédiaires entre les signaux temporels bruts ininterprétables et les paramètres est peu explorée. Certains paramètres nécessitent la définition de seuils pour être calculés. Ainsi, on définit par exemple qu'un sujet a atteint un régime établi lorsque l'amplitude de ses cycles a dépassé 67% de l'amplitude des cycles n°5 à 10 (on considère qu'il n'y a plus lieu de parler d'initiation à plus de 5 cycles). Certains cycles peuvent avoir des valeurs limites (60% ou 70% de l'amplitude des cycles à distance). On commet ainsi des erreurs liées au seuil choisi qui pourraient être évitées en regardant directement la donnée brute. 3061787 2 Les paramètres moyennent les caractéristiques de la marche sur l'ensemble de l'exercice, et occultent de ce fait le déroulé de l'exercice. Par exemple, pour évaluer la régularité de la marche d'un sujet, on utilise la déviation standard de la durée des cycles 5 de marche. Si les cycles sont globalement réguliers et qu'un cycle se trouve plus long que les autres, la déviation standard sera augmentée. On ne saura pas que cette augmentation sera due à un seul cycle isolé, ni à quel moment de la marche celui-ci a eu lieu. Or, le nombre de cycles erratiques, et leur instant de survenu sont des informations analytiques importantes. 10 Les paramètres sont des bons indicateurs de la performance de la marche d'un sujet. Ils n'apportent néanmoins pas d'information analytique sur l'origine de la modification de la marche. Par exemple, la vitesse de marche est une bonne évaluation globale, cependant, une vitesse peut être diminuée par des pas plus courts ou une cadence 15 diminuée. A leur tour, la taille des pas et la cadence peuvent être affectés par différents facteurs. La mesure de la marche présente d'importantes fluctuations interindividuelles dues principalement aux variabilités anatomiques et fonctionnelles de chaque individu qui influent sur le positionnement des capteurs, et aux variabilités de style locomoteur de chaque individu. Les paramètres sont affectés par cette variabilité et il n'existe pas de méthode de référence pour s'en affranchir. EXPOSE DE L'INVENTION Il existe un besoin dans l'état de l'art de disposer d'un procédé qui permet d'étudier la régularité et la symétrie d'une séquence de cycles de marche ou de course d'une personne. 3061787 3 Un but de l'invention est de permettre de prendre en compte la forme du signal pour rendre compte de la régularité de la marche ou de course, tout en s'affranchissant de la variabilité de la durée ou de l'amplitude des cycles de marche ou de la course. 5 Dans l'invention, il est construit une matrice ou tableau de visualisation d'un exercice de marche ou de course. Cette matrice de visualisation M(i,j) est dite « locogramme », et est une matrice de Gram au sens ou xij=xji avec xij étant la valeur située à la ième ligne et à la jème colonne dans la matrice M(i,j). 10 Grâce à la matrice de visualisation de l'invention, le médecin peut avoir une vue globale de l'exercice de marche à partir de laquelle il peut faire la sémiologie afin de se faire une idée de la qualité de la marche d'un patient et repérer aisément dans un exercice donné des pas atypiques. Les seuls paramètres ne suffisent pas pour avoir du patient un aperçu satisfaisant. La matrice de visualisation est à mi-chemin entre les 15 paramètres et les signaux temporels bruts, et donne accès au médecin aux données de la marche du patient. La matrice de visualisation est fondée sur la division de la marche ou de course en cycle. Il s'agit de comparer les cycles entre eux à l'aide de « distances » ou métriques 20 mathématiques. En fonction du choix de la distance, tel ou tel aspect du signal va être comparé. Plus précisément, l'invention présente un dispositif permettant l'analyse de la régularité et de la symétrie d'une séquence de N cycles de marche d'une personne, comportant: 25 - des capteurs de mesure de signaux bruts d'une grandeur physique d'un segment anatomique - une unité de traitement reliée aux capteurs de mesure pour : - traiter et séparer les signaux temporels bruts en signaux temporels distincts Ci, 30 chaque signal temporel Ci étant une série de points de la grandeur physique mesurée 3061787 4 et présentant une forme, une amplitude et une durée donnée, la série Ci étant associée à un cycle de marche ou de course donné i de la personne; - calculer au moins un coefficient de similarité xii entre le signal Ci associé au cycle i, avec un autre signal Cj associé à un cycle j de la même personne, 5 le coefficient de similarité xij étant : un coefficient de similarité de la forme (fij) des deux signaux Ci et Cj, un coefficient de similarité de l'amplitude (aij) des deux signaux Ci et Cj, ou un coefficient de similarité de la durée (dij) des deux signaux Ci et Cj, -ranger dans une matrice carrée M(i,j) (F(i,j), A(i,j) D(i,j)), dans chaque cellule (i,j) à la ligne i et à la colonne j, la valeur du coefficient de similarité xij (fij, aij, dij) ; 10 avec i et j entiers naturels variant de 1 à N, les N cycles de marche étudiés étant ordonnés chronologiquement ; - des moyens d'affichage reliés à l'unité de traitement et affichant la matrice carrée M(i,j) (F(i,j), A(i,j) D(i,j)) avec chaque valeur du coefficient de similarité xij (fij, 15 aij, dij) représentée dans la matrice M(i,j) (F(i,j), A(i,j) D(i,j)) par une représentation visuelle graduée pour permettre de visualiser à l'oeil nu la similarité entre les cycles de marche ou de course i et j. L'approche de la locomotion par la matrice de visualisation suggère de décrire un cycle 20 (marche, course, piétinement) selon trois aspects sur lesquels vont porter des calculs de distance de similarité : la durée, l'amplitude et la forme du cycle. L'amplitude et la forme du signal sont des aspects indépendants. Avantageusement, pour décrire les pas (marche, course, piétinement) d'un sujet, il peut 25 être utile de construire trois matrices de visualisation F(i,j), A(i,j) D(i,j) à partir de ces trois distances. La présente invention concerne aussi un procédé utilisant le dispositif défini ci-avant, 30 comportant comme étapes: 3061787 5 -une étape de détection (i) des signaux temporels bruts d'une grandeur physique d'un segment anatomique, mesurée lors d'une séquence de N cycles de marche ou de course i de la personne, 5 - une étape de traitement et de séparation (ii) des signaux temporels bruts en signaux temporels distincts Ci, chaque signal temporel Ci étant une série de points de la grandeur physique mesurée et présentant une forme, une amplitude et une durée donnée, la série Ci étant associée à un cycle de marche ou de course donné i de la personne; 10 - une étape de calcul (iii) d'au moins un coefficient de similarité xij (fij, aij, dij) entre un signal Ci associé au cycle de marche ou de course i, avec un autre signal Cj associé à un cycle de marche ou de course j, représentatif de la similarité entre les deux signaux Ci et Cj ; 15 le coefficient étant : un coefficient de similarité de la forme fij des deux signaux Ci et Cj, un coefficient de similarité de l'amplitude aj des deux signaux Ci et Cj, ou un coefficient de similarité de la durée dij des deux signaux Ci et Cj; - une étape de rangement (iv) dans une matrice carrée M(i,j), de la valeur du 20 coefficient de similarité xij à la ligne i et à la colonne j ; avec i et j entiers naturels variant de 1 à N, les N cycles de marche ou de course étant ordonnés chronologiquement; - une étape d'affichage (y) affichant la matrice carrée M(i,j), chaque valeur du 25 coefficient de similarité xij étant représentée dans la matrice M(i,j) par une représentation visuelle graduée pour permettre de visualiser à l'oeil nu la similarité entre les cycles de marche ou de course ou de piétinement étudiés i et j. DESCRIPTION DES FIGURES 30 D'autres objectifs, caractéristiques et avantages sortiront de la description détaillée qui 3061787 6 suit en référence aux dessins donnés à titre illustratif et non limitatif parmi lesquels : - la figure I a représente le dispositif selon l'invention ; - la figure 1 b représente le procédé selon l'invention ; 5 - la figure 2 représente la définition des termes amplitude, durée et forme des cycles de marche ou de course ; - la figure 3 représente l'outil de visualisation selon l'invention du cycle de marche d'un pied d'une personne sur un exercice de marche donné comportant un aller et retour sur une distance définie; 10 - les figures 4a et 4b représentent l'outil de visualisation F(i,j) qui se rapporte au coefficient de similarité de forme selon l'invention pour un sujet sain, respectivement pour un sujet parkinsonien, comparant pour chaque sujet, la forme de ses cycles de marche sur un exercice de marche donné comportant un aller et retour sur une distance définie; 15 - les figures 5a et 5b représentent l'outil de visualisation de la marche ou matrice carrée A(i,j) qui se rapporte au coefficient de similarité de l'amplitude selon l'invention pour un sujet sain, respectivement pour un sujet parkinsonien, comparant pour chaque sujet, l'amplitude de ses cycles de marche sur un exercice de marche donné comportant un aller et retour sur une distance 20 définie; - les figures 6a et 6b représentent l'outil de visualisation de la marche ou matrice carrée D(i,j) qui se rapporte au coefficient de similarité de durée selon l'invention pour un sujet sain, respectivement pour un sujet parkinsonien, comparant pour chaque sujet, la durée de ses cycles de marche sur un 25 exercice de marche donné comportant un aller et retour sur une distance définie; - la figure 7 représente, en fonction du temps, les signaux temporels bruts de l'accélération issus d'un accéléromètre triaxial selon les trois axes x, y et z, pour le pied droit d'une personne en mouvement sur un exercice de marche 30 donné comportant un aller et retour sur une distance définie; 3061787 7 - la figure 8 représente, en fonction du temps, les signaux temporels bruts de l'accélération sans gravité correspondant aux signaux bruts de la figure 1 auxquels on a soustrait les signaux temporels bruts issus de l'accéléromètre triaxial selon les trois axes x, y et z pour le pied droit d'une personne immobile ; 5 - la figure 9 représente, en fonction du temps, la norme de l'accélération sans gravité du pied droit de la personne sur un exercice de marche donné comportant un aller et retour sur une distance définie; - les figures 10a, 10b représentent, en fonction du temps, la norme de l'accélération sans gravité du pied droit de la personne, la norme de 10 l'accélération sans gravité du pied gauche de la personne et le repérage des temps -r_droit et -r_gauche correspondant aux cycles de marche C1, 02....Cnd+ng sur un exercice de marche donné comportant un aller et retour sur une distance définie; - la figure 10c représente le ré-échantillonnage des cycles de marche C1, 15 02....Cnd+ng. Les figures présentées dans ce document ont uniquement une valeur illustrative. DESCRIPTION DETAILLEE DE L'INVENTION 20 La présente invention concerne un dispositif 1 permettant l'analyse de la régularité et de la symétrie d'une séquence de N cycles de marche ou de course d'une personne, sans calcul de paramètres basés sur la considération globale des N cycles de marche 25 ou de course. L'expression 'cycle de marche ou de course' correspond aux phases définies entre deux appuis du talon d'un même pied sur le sol lors de la marche ou d'une course de la personne et qui ont lieues à certains instants entre les deux appuis du talon. Elle est 30 donc représentative du comportement du pied de la personne pendant un pas. 3061787 8 Comme il est connu de l'état de l'art, les différentes phases peuvent être définies en pourcentage, 0 % correspondant au premier appui talon et 100 % au deuxième. Le cycle peut être séparé en deux phases : la phase d'appui (0 à 60%) où le pied est 5 en contact avec le sol, et la phase oscillante ou pendulaire (60% à 100%) où le talon est en l'air. Cette expression définit tous les styles de marche et les styles de course ; par exemple le piétinement doit être considéré comme faisant partie de la catégorie définie par 10 l'expression 'cycle de marche'. Le dispositif 1 illustré sur la figure 1 a comporte des capteurs de mesure 2 de signaux bruts i d'une grandeur physique d'un segment anatomique, mesurée lors d'une séquence de N cycles de marche ou de course d'au moins un pied de la personne. 15 Les capteurs de mesure 2 peuvent être tout capteur, par exemple tel qu' : un accéléromètre, un gyroscope, une électromyographie, des semelles de pression ou des traceurs infrarouges/dispositifs d'acquisition vidéo ou IR. 20 Ils peuvent être placés par exemple sur le pied, la cheville, la ceinture, poignet, tête. La grandeur physique mesurée peut, par exemple, être la norme de l'accélération, la norme de l'accélération sans la composante de gravité terrestre, la vitesse, la vitesse angulaire, le déplacement, la position etc... 25 Par exemple, les grandeurs issues de capteurs inertiels peuvent être l'accélération linéaire ou la vitesse angulaire selon un des 3 axes de l'espace. Avec d'autres capteurs, d'autres grandeurs peuvent être la position du pied dans l'espace avec la stéréophotogrammétrie ou la pression au sol avec des semelles de force à capteur de pression. 30 3061787 9 Comme représenté sur la figure 1 a, le dispositif 1 comporte une unité de traitement et de calcul 3 agencée pour : - traiter et séparer les données brutes en signaux temporels distincts Ci, chaque signal 5 temporel Ci étant une série de points de la grandeur physique mesurée fonction du temps et présentant une forme, une amplitude et une durée donnée, la série Ci étant associée à un seul cycle donné i de la personne parmi les N cycles de marche comme représenté sur la figure 2 ; 10 - calculer au moins un coefficient de similarité entre un signal Ci associé à un cycle i, avec un autre signal Cj associé à un autre cycle j, représentatif de la similarité entre les deux signaux Ci et Cj ; -ranger dans une matrice carrée M(i,j) indexée en ligne i et en colonne j, et indiquant 15 pour chaque cellule (i, j), la valeur du coefficient de similarité entre les cycles i et j. Lors de ces calculs et du rangement des valeurs déterminées dans la matrice M(i,j), les entiers naturels avec i et j varient de 1 à N, et les N cycles de marche étudiés sont ordonnés pour le pied donné, chronologiquement. 20 Le traitement (« compend resampling ») et la séparation des données brutes est illustré sur les figures 10A à 10D, et les signaux temporels distincts Ci, après traitement - dessus sont nommés sur la figure 10D (ce qui correspond à la terminologie C'i dans le texte). 25 Le coefficient de similarité peut être : un coefficient de similarité de la forme fij des deux signaux Ci et Cj (comme représenté sur les figures 4a et 4b), un coefficient de similarité de l'amplitude aij des deux signaux Ci et Cj (comme représenté sur les figures 5a et 5b), ou un coefficient de similarité de la durée dij des deux signaux Ci et Cj 30 (comme représenté sur les figures 6a et 6b). 3061787 10 Un synonyme pour 'coefficient' est 'distance algébrique' au sens mathématique du terme par exemple, et pour 'similarité' est ressemblance. 5 L'expression `de similarité' implique la symétrie des coefficients et donc que fij=fji, aij=aji et dij=dji. Le coefficient de similarité peut être considéré comme étant un coefficient de corrélation dans un certain nombre de cas de figures. 10 D'autres coefficients de similarité peuvent être utilisés en fonction de l'aspect de la forme que l'on veut comparer entre deux signaux : le Dynamic Time Wraping (DTW), le coefficient de Spearman, les distances euclidiennes (L1, L2 et Loo). Ainsi, avec un coefficient de similarité qui varie dans un intervalle [a ; b], 15 - si deux signaux (ou les deux cycles de marche ou de course) sont identiques, le coefficient peut être égal à b, et - si deux signaux (ou les deux cycles de marche ou de course) n'ont aucune ressemblance, le coefficient peut être égal à a, la similarité entre les signaux augmentant linéairement entre a et b. 20 Avantageusement, l'unité de traitement et de calcul 3 est agencée pour : - calculer plusieurs coefficients de similarité fij, aij, dij différents, chaque coefficient de similarité fij, aij, dij étant associé à la même échelle de couleurs, chaque couleur ayant la même signification de similarité quelle que soit le coefficient utilisé fij, aij, dij, comme 25 illustré sur les figures 4a à 6b ; - présenter les valeurs des coefficients de similarité fij, aij, dijj dans des matrices carrées différentes F(i,j), A(i,j) D(i,j), une matrice carrée par type de coefficient de similarité . 30 Le dispositif 1 présente aussi des moyens d'affichage 4 reliés à l'unité de traitement et affichant la matrice M(i,j) pour i et j variant de 1 à N, chaque valeur du coefficient de 3061787 11 similarité étant représentée dans la matrice M(i,j) par une représentation visuelle graduée pour permettre de visualiser à l'oeil nu la similarité entre les cycles de marche ou de course i et j. 5 Avantageusement, les coefficients de similarité (fij, dij, aij) sont choisis : - pour qu'une fois calculées, toutes les valeurs des coefficients de similarité (fij, dij, aij) soient dans un même intervalle [a ; b] et - pour que plus les valeurs des coefficients de similarité (fij, dij, aij) soient élevées et plus les signaux Ci, Cj soient similaires. 10 Avantageusement, cette représentation visuelle est graduée pour représenter la valeur du coefficient de similarité. En d'autres termes, il y a une correspondance entre l'échelle des valeurs de coefficient de similarité et celle des couleurs. 15 Cette correspondance peut être proportionnelle ou non. De façon préférée, chaque valeur est représentée par une couleur située dans une échelle de couleurs graduées pour permettre de visualiser à l'oeil nu la similarité entre 20 les deux cycles de marche ou de course i et j (ce qui revient à étudier la similarité entre les deux signaux Ci et Cj). Avantageusement, l'échelle utilisée par les moyens d'affichage 4 est choisie sans seuils ou continue comme illustré sur les figures 3 à 6b. 25 Ici l'échelle est [0 ; 1] pour le coefficient de similarité de forme, le coefficient de similarité d'amplitude et le coefficient de similarité de durée permettant une comparaison rapide et aisée des différents coefficients de similarité entre eux. 30 Chaque case de la matrice peut porter une couleur correspondant au degré de similarité/ressemblance de deux pas entre eux. 3061787 12 Une couleur chaude (exemple rouge) est proche de 1, et traduit une grande ressemblance (similarité). Une couleur froide (exemple bleu) est proche de 0, et traduit une faible ressemblance 5 (similarité). A noter que cette échelle est inversée entre la figure 3 d'un côté et les figures 4a à 6b de l'autre : sur la figure 3 des couleurs foncées indiquent une ressemblance, ce qui est le contraire sur les figures 4a à 6b ou des couleurs claires indiquent une ressemblance. 10 Les moyens d'affichage 4 affichent une représentation visuelle de la matrice de similarité telle que : -chaque case représente un coefficient de similarité entre le signal Ci et Cj ; -la valeur du coefficent de similarité (fij, aij, dij) est représentée dans une échelle de 15 couleur graduée. Dans un premier mode de réalisation illustré sur la figure 3, il est représenté le coefficient de similarité entre le signal Ci relatif à la marche ou à la course d'un seul pied (pied droit, ou respectivement pied gauche) de la personne, et le signal Cj relatif 20 à la marche ou à la course du même pied, et nombre de cycles N=Nd (respectivement Ng) nombre de cycles du pied droit (respectivement nombres de cycles du pied gauche). Dans un second mode de réalisation illustré sur les figures 4a à 6b, les moyens 25 d'affichage 4 affichent: -les coefficients de similarité (fij, aij, dij) calculés pour un seul pied 3 ; - également le coefficient de similarité (fij, aij, dij) entre le signal Ci relatif à la marche ou à la course du pied droit de la personne, et le signal Cj relatif à la marche ou à la course du pied gauche de la personne, et N=Nd +Ng. 30 3061787 13 Dans ce cas, tous les cycles étudiés du pied droit, puis tous les cycles étudiés du pied gauche (ou inversement) sont ordonnés chronologiquement, en abscisse et en ordonnée de la matrice, comme illustré sur les figures 4a à 6b. Ils sont séparés ici dans la représentation par une zone O, ce qui donne 4 sous-carrés. 5 Ainsi, la matrice de similarité : - « pied droit/pied droit » est représentée par le sous-carré en bas à gauche ; - « pied gauche/pied gauche » est représentée par le sous-carré en haut à droite ; - « pied gauche/pied droit » est représentée par le sous-carré en bas à droite. 10 Lorsque le coefficient de similarité est un coefficient de forme fij, l'unité de traitement et de calcul 3 est agencée pour : -normaliser en durée chaque signal temporel Ci avec chaque signal temporel Cj pour que les deux signaux Ci et Cj aient la même durée ; 15 - normaliser en amplitude chaque signal temporel Ci ; - calculer un coefficient de similarité de forme fij entre chaque signal normalisé C'i, avec un autre signal normalisé C'j, et ce pour tous les cycles de marche. Dit autrement, pour étudier la ressemblance de la forme des cycles entre eux, il faut 20 s'affranchir de leur différence d'amplitude et de leur différence de durée. Par exemple, la normalisation en durée est réalisée par ré-échantillonage linéaire, ou par DTW, comme illustré sur la figure 10c. 25 Par exemple, la normalisation en amplitude est réalisée en divisant le Signal Ci par l'écart type ou par valeur efficace (« Root Mean Square »). Dans une première réalisation, le coefficient de similarité de forme fij peut être le coefficient de Pearson, comme illustré sur les figures 4a et 4b, les valeurs étant 30 ramenées à 0 en cas de valeurs du coefficient de Pearson inférieures à zéro. 3061787 14 D'autres distances mathématiques sont envisageables telles que la distance de corrélation de Spearman ou la technique du « dynamique time wraping ». Si l'on s'intéresse uniquement à la forme, il est important de s'affranchir de l'amplitude et de la durée. Pour cela, il est re-normalisé le signal en ré-échantillonant le signal à 5 100 échantillons (normalisation en durée) et divisant le signal par l'écart-type du cycle (normalisation en amplitude). La distance de corrélation de Pearson conserve la ligne temporelle et compare strictement la forme des pas. La distance de corrélation de Spearman et le « dynamic time wraping » déforment la ligne temporelle et indique s'il existe une déformation possible pour que la forme de deux cycles se ressemble ou si 10 les formes sont vraiment différentes. Dans une deuxième réalisation, le coefficient de similarité peut être un coefficient de similarité de durée dij. 15 Dans ce cas, comme illustré sur les figures 5a et 5b, l'unité de traitement et de calcul 3 est agencée par exemple pour : -calculer une durée Di, Dj du signal Ci et de la durée du signal Cj, -calculer le coefficient de similarité de durée dij qui est le rapport de la plus petite des 20 deux durées Di, Dj divisée par la plus grande durée Di, Dj, afin de ramener la valeur des coefficients entre 0 et 1. La durée du signal Ci est égal au nombre d'échantillons, chaque échantillon étant pris de façon régulière à une fréquence donnée. 25 Dans une troisième réalisation, le coefficient peut être un coefficient de similarité d'amplitude. Dans ce cas, comme illustré sur les figures 6a et 6b, l'unité de traitement et de calcul 30 3 est agencée par exemple pour : -calculer une amplitude Ai, Aj du signal Ci et du signal Cj, 3061787 15 -calculer le coefficient de similarité d'amplitude aij qui est le rapport de la plus petite amplitude Ai, Aj divisée par la plus grande amplitude Ai, Aj afin de ramener la valeur des coefficients entre 0 et 1. 5 L'amplitude Ai et Aj est par exemple l'écart type ou la valeur efficace (« Root Mean Square »). Une fois ces calculs effectués pour trouver les valeurs des coefficients de similarité, l'unité de traitement et de calcul 3 est agencée pour calculer : 10 -des paramètres permettant d'évaluer la régularité de la marche ; -des paramètres permettant d'évaluer la symétrie de la marche ; -le nombre de cycles nécessaires à établir une séquence ou une sous séquence de marche. 15 En particulier : - le paramètre évaluant la régularité de la marche : la moyenne de toutes les cellules (i, j) de la matrice M(i,j) (A(i, j) ; D(i,j) ; F(i,j) ) en excluant les cellules i=j, - le paramètre évaluant la régularité de la marche : la déviation standard de toutes les cellules (i, j) de la matrice M(i,j) en excluant les cellules i=j, 20 - le paramètre évaluant la symétrie de la marche : la moyenne des coefficients des cellules avec i allant de 1 à Nd (ou de 1 à Ng), et j allant de Nd+1 à Nd+Ng (j allant de Ng+1 à Nd+Ng). La présente invention concerne aussi un procédé illustré sur la figure 2 et utilisant le 25 dispositif 1 ci-avant décrit. Le procédé, illustré sur la figure 2, comporte : -une étape de détection (i) des signaux temporels bruts ; - une étape de traitement et de séparation (ii) des signaux temporels bruts en signaux 30 temporels distincts Ci; - une étape de calcul (iii) de coefficient de similarité ; 3061787 16 -une étape de rangement (iv) dans une matrice M(i,j), la valeur du coefficient de similarité à la ligne i et à la colonne j ; avec i et j entiers naturels variant de 1 à N, les N cycles de marche ou de course étant ordonnés chronologiquement; 5 - une étape d'affichage (y) affichant la matrice M(i,j) ; - une étape de visualisation (vi) par un opérateur. L'étape de visualisation par l'opérateur permet de déterminer : - le nombre de cycles nécessaires à la réalisation d'une séquence de marche ou de 10 course ou d'une sous-séquence de marche, la sous-séquence de marche ou de course étant: un demi-tour, l'initiation de la marche ou de la course, une accélération, une décélération, un régime établi, l'arrêt, un virage ; - le nombre de cellules ayant une couleur associée à une valeur faible de coefficient ; - si un cycle présente une ou plusieurs cellules ayant une couleur associée à une valeur 15 faible ou élevée de coefficient ; - si les deux cycles de marche consécutifs présentent une ou plusieurs cellules ayant une couleur associée à une valeur faible ou élevée de coefficient; - repérer les changements de régime de pas au cours de l'exercice par des hétérogénéités de couleur ; 20 -comparer la couleur des cellules pied-droit-pied droit avec les cellules pied droit-pied gauche. Il peut être prévu que l'unité de traitement et de calcul 3 compte le nombre de cycles ci-dessus définis, ou ne fait apparaitre que certaines valeurs de coefficients de 25 similarité entre les cycles au vu de valeurs de seuil prédéfinies. DESCRIPTION DU PROTOCOLE DE REALISATION D'UN EXERCICE DE MARCHE ET DESCRIPTION DE LA COHORTE DE PARKINSONIENS D'APRES AVC 2. 30 3061787 17 Etape 1 : Acquisition des données Un sujet portant un accéléromètre triaxial sur la face dorsale de chaque pied réalise un exercice de marche de 10 mètres aller-retour départ arrêté. Chacun des capteurs 5 enregistre un signal temporel en 3 dimensions (selon les trois axes x, y et z). On appelle respectivement accdroit et accgauche les signaux pour le pied droit et le pied gauche, comme représenté sur les courbes de la figure 7. 10 Etape 2 : Repérage des débuts des cycles de marche On repère (manuellement ou automatiquement) dans les signaux, les temps correspondant aux débuts des cycles de marche. Ces instants, selon la définition de Mariani et al. (2013), correspondent aux « heel-strikes », c'est-à-dire aux instants où le 15 talon touche le sol. On définit deux ensembles T - droit et T gauche correspondant respectivement à l'ensemble des heel-strike pour le pied droit et le pied gauche (figures 10A et 10B). On appelle Nd + 1 et N9 + 1 respectivement le nombre de heel-strike repérés pour le pied droit et pour le pied gauche. Le cycle de marche i du pied droit est défini comme la partie du signal commençant à l'instant T - droit,i et finissant au début du 20 cycle de marche suivant T droit,i+11 comme le montrent les figures 10A et 10B. Ainsi, Nd et N, sont les nombres de cycle de marche du pied droit et du pied gauche respectivement. 25 Etape 3 : Retrait de la gravité 3061787 18 Au début du protocole, ou dans une étape préliminaire, on a demandé au sujet de rester debout et immobile. On a enregistré la valeur (constante) des accélérations 7 accimmobile droit,x\ durant cette phase et on l'a stocké dans un vecteur accimmobiledroity \accimmobile Comme représenté sur les courbes de la figure 8, il est enlevé la gravité selon le 5 procédé suivant : accdroit,x 7 accimmobile droitA ( = acc droit y - accimmobile droit y accdroit,z \accimmobile droit,z i On réitère le même procédé pour le pied gauche. Etape 4 : Calcul de la norme 10 Comme représenté sur la courbe de la figure 9, il est calculé la norme du vecteur accélération à chaque instant t : sdroit [t] = ^l accfreedroit,x[t]2+ accfreedroity RF + accfreedroit,z[t]2 de même pour le pied gauche. 15 Cette fusion des trois axes de l'accélération permet de se rendre indépendant de la position du capteur qui est une source d'imprécision de mesure importante dans la mesure par accéléromètre. Etape 5 : Création des cycles de marche 20 Comme représenté sur la courbe des figures 10C à 10D, on note l'ensemble des cycles de marche de l'exercice de la manière suivante : 25 Pour Ci = S droit[T droit,i)T droit,i+1] Si 1 REVENDICATIONS1. Dispositif (1) permettant l'analyse de la régularité et de la symétrie d'une séquence de N cycles de marche ou de course d'une personne, comportant : - des capteurs de mesure (2) de signaux temporels bruts d'une grandeur physique d'un segment anatomique, - une unité de traitement et de calcul (3), reliée aux capteurs de mesure (2), et agencée pour : - traiter et séparer les signaux temporels bruts en signaux temporels distincts Ci, chaque signal temporel Ci étant une série de points de la grandeur physique mesurée et présentant une forme, une amplitude et une durée donnée, la série Ci étant associée à un cycle de marche ou de course donné i de la personne; - calculer au moins un coefficient de similarité entre le signal Ci associé au cycle de marche ou de course i, avec un autre signal Cj associé à un cycle de marche ou de 15 course j de la même personne, le coefficient de similarité étant : un coefficient de similarité de la forme fij des deux signaux Ci et Cj avec fij=fji, ou un coefficient de similarité de l'amplitude aij des deux signaux Ci et Cj avec aij=aji, ou un coefficient de similarité de la durée dij des deux signaux Ci et Cj avec dij=dji, 20 -ranger dans une matrice carrée M(i,j), la valeur du coefficient de similarité à la ligne i et à la colonne j; avec i et j entiers naturels variant de 1 à N, les N cycles de marche ou de course étant ordonnés chronologiquement, - des moyens d'affichage (4) reliés à l'unité de traitement et de calcul (3), et 25 affichant la matrice carrée M(i,j) avec les valeurs du coefficient de similarité représentées dans les cellules (i,j) de la matrice carrée M(i,j), par une représentation visuelle graduée, pour permettre de visualiser la similarité entre l'ensemble des cycles de marche ou de course i et j. 30 3061787 31 2. Dispositif (1) selon la revendication précédente 1, caractérisé en ce que les moyens d'affichage (4) représentent chaque valeur du coefficient de similarité par une couleur située dans une échelle continue graduée en correspondance avec l'échelle des coefficients de similarité. 5 3. Dispositif (1) selon l'une des revendications précédentes, caractérisé en ce que l'unité de traitement et de calcul (3) détermine le coefficient de similarité entre le signal Ci relatif à la marche ou à la course du pied droit (respectivement du pied gauche) de la personne, et le signal Cj relatif à la marche ou à la course du pied droit 10 (respectivement du pied gauche) de la même personne, et en ce que N=Nd (respectivement Ng) qui est le nombre de cycles du pied droit (respectivement le nombre de cycles du pied gauche). 4. Dispositif (1) selon la revendication précédente 3, caractérisé en ce que l'unité de traitement et de calcul (3) détermine également le coefficient de similarité entre le signal Ci relatif à la marche ou à la course du pied droit de la personne, et le signal Cj relatif à la marche ou à la course du pied gauche de la même personne, et en ce que le nombre de cycles N est égal à la somme du nombre de cycles du pied droit et du nombre de cycles du pied gauche : N=Nd +Ng. 5. Dispositif (1) selon l'une des revendications précédentes 2 à 4, caractérisé en ce que l'unité de traitement et de calcul (3) est agencée pour : - calculer les trois coefficients de similarité différents fij, aij, dij, chaque coefficient de similarité fij, aij, dij étant associé à la même échelle de couleurs que celle des autres coefficients de similarité, chaque couleur ayant la même signification de similarité quel que soit le coefficient de similarité utilisé ; - présenter les valeurs des coefficients de similarité fij, aij, dij dans des matrices carrées différentes (F(i,j), A(i,j),D(i,j)); les coefficients de similarité étant choisis : 3061787 32 - pour qu'une fois calculées, toutes les valeurs des coefficients de similarité fij, aij, dij soient dans un même intervalle [a ; b] et - pour que plus les valeurs des coefficients de similarité fij, aij, dij soient élevées et plus les signaux Ci, Cj soient similaires. 5 6. Dispositif (1) selon l'une des revendications précédentes, caractérisé en ce que le coefficient de similarité est un coefficient de forme fij, et en ce que l'unité de traitement et de calcul (3) est agencée pour : -normaliser en durée chaque signal temporel Ci avec chaque signal temporel Cj pour 10 que les deux signaux Ci et Cj aient la même durée ; - normaliser en amplitude chaque signal temporel Ci ; - calculer un coefficient de similarité de forme fij entre chaque signal normalisé C'i, avec un autre signal normalisé C'j, et ce pour tous les cycles de marche ou de course. 15 7. Dispositif (1) selon la revendication précédente, caractérisé en ce que le coefficient de similarité de forme fij est le coefficient de Pearson, les valeurs du coefficient de similarité étant ramenées à 0 en cas de valeur du coefficient de Pearson inférieure à zéro. 20 8. Dispositif (1) selon l'une des revendications précédentes, caractérisé en ce que le coefficient de similarité est un coefficient de similarité de durée dij, et en ce que l'unité de traitement et de calcul (3) est agencée pour : -calculer une durée Di, Dj du signal Ci et de la durée du signal Cj, -calculer le coefficient de similarité de durée dij qui est le rapport de la plus petite des 25 deux durées Di, Dj divisée par la plus grande durée Di, Dj. 9. Dispositif (1) selon l'une des revendications précédentes, caractérisé en ce que le coefficient de similarité est un coefficient de similarité d'amplitude aij, et en ce que l'unité de traitement et de calcul (3) est agencée pour : 30 -calculer une amplitude Ai, Aj du signal Ci et du signal Cj, 3061787 33 -calculer le coefficient de similarité d'amplitude aij qui est le rapport de la plus petite amplitude Ai, Aj divisée par la plus grande amplitude Ai, Aj. 10. Dispositif (1) selon la revendication précédente, caractérisé en ce que l'amplitude 5 Ai et Aj est l'écart type ou la valeur efficace. 11. Dispositif (1) selon l'une des revendications précédentes, caractérisé en ce que le cycle considéré est le cycle de la marche et en ce que l'unité de traitement et de calcul (3) est agencée pour calculer : 10 -des paramètres permettant d'évaluer la régularité de la marche ou de la course ; -des paramètres permettant d'évaluer la symétrie de la marche ou de la course ; -le nombre de cycles nécessaires à une séquence ou une sous séquence de marche ou de course. 15 12. Dispositif (1) selon la revendication précédente, caractérisé en ce que l'unité de traitement et de calcul (3) est agencée pour calculer : -la moyenne de toutes les cellules (i, j) de la matrice carrée M(i,j) en excluant les cellules i=j, -la déviation standard de toutes les cellules (i, j) de la matrice carrée M(i,j) en excluant 20 les cellules i=j, - la moyenne des cellules i allant de 1 à Nd (ou de 1 à Ng), et j allant de Nd+1 à Nd+Ng (j allant de Ng+1 à Nd+Ng). 13. Dispositif (1) selon l'une des revendications précédentes, caractérisé en ce que les 25 capteurs de mesure (2) sont : accéléromètre, gyroscope, électromyographie, semelles de pression ou dispositif d'acquisition cinématique infrarouge. 14. Dispositif (1) selon l'une des revendications précédentes, caractérisé en ce que la grandeur physique mesurée est la norme de l'accélération sans la composante de 30 gravité terrestre. 3061787 34 15. Procédé permettant l'analyse de la régularité et de la symétrie de la marche ou de course d'une personne et utilisant le dispositif (1) selon l'une revendications précédentes 1 à 14, comportant : 5 -une étape de détection (i) des signaux temporels bruts d'une grandeur physique d'un segment anatomique, mesurée lors d'une séquence de N cycles de marche ou de course de la personne, - une étape de traitement et de séparation (ii) des signaux temporels bruts en signaux 10 temporels distincts Ci, chaque signal temporel Ci étant une série de points de la grandeur physique mesurée et présentant une forme, une amplitude et une durée donnée, la série Ci étant associée à un cycle de marche ou de course donné i de la personne; 15 - une étape de calcul (iii) d'au moins un coefficient de similarité (aij, fij, dij) entre un signal Ci associé au cycle de marche ou de course i, avec un autre signal Cj associé à un cycle de marche ou de course j, représentatif de la similarité entre les deux signaux Ci et Cj ; le coefficient de similarité étant un coefficient de similarité de la forme fij des deux 20 signaux Ci et Cj avec fij=fji, ou un coefficient de similarité de l'amplitude aij des deux signaux Ci et Cj avec aij=aji, ou un coefficient de similarité de la durée dij des deux signaux Ci et Cj avec dij=dji; - une étape de rangement (iv) dans une matrice carrée M(i,j), de la valeur du 25 coefficient de similarité (aij, fij, dij) à la ligne i et à la colonne j ; avec i et j entiers naturels variant de 1 à N, les N cycles de marche ou de course étant ordonnés chronologiquement; - une étape d'affichage (y) affichant la matrice carrée M(i,j), chaque valeur du 30 coefficient de similarité (aij, fij, dij) étant représentée dans la matrice carrée M(i,j) par 3061787 une représentation visuelle graduée pour permettre de visualiser à l'oeil nu la similarité entre les cycles de marche ou de course i et j. 16. Procédé selon la revendication précédente 15, dans lequel les signaux temporels 5 bruts sont traités et séparés préalablement au calcul du coefficient de similarité automatiquement ou manuellement à l'étape de calcul (iii).